#!/usr/local/bin/python -i

import aa
from ROOT import *


cxx(""" // pyroot seems broken, so do this in c++
#include "Det.hh"

void set_dom_positions(Det& d)
{
  foreach_map( k,v, d.doms ) 
   {
   v.pos = Vec();
   foreach( pmt, v.pmts ) v.pos+=pmt.pos;
   v.pos/=v.pmts.size();
   v.pmts.clear();
   }
}""");


d = Det("/sps/km3net/users/tseitz/detector/km3net_115.det")
set_dom_positions( d );

f = EventFile("/sps/km3net/users/heijboer/example_files/km3_v4_nueCC_1.JTE.root");


 
j = dumper();
for evt in f :
    print evt.hits.size()
    if evt.hits.size() < 500 : continue
    d.apply(evt);
    e = evt;
    break



outf = open ("/afs/in2p3.fr/home/h/heijboer/www/detevt_km3.js","w")
outf.write( "det = "+ j(d)+ ";\n" );
outf.write( "evt = "+ j(e)+ ";\n" );



outf.close()