6 #define WAVE_FILE "merge/wave_bbh_low_lh_wp10_run2.M1.C_salvo.root"
12 int nevt =
event.GetEntries();
14 cout <<
"no events are presents in the file : " <<
WAVE_FILE << endl;
16 }
else cout <<
"nevt : " << nevt << endl;
19 TList* ifoList =
event.fChain->GetUserInfo();
20 int nIFO = ifoList->GetSize();
24 cout << dParams.
name << endl;
34 cout <<
"-----------------------------------" << endl;
35 cout <<
"run : " <<
event.run << endl;
36 cout <<
"gps : " <<
int(event.time[0]) << endl;
37 cout <<
"rho : " <<
event.rho[1] << endl;
38 cout <<
"cc : " <<
event.netcc[0] << endl;
39 cout <<
"network snr : " << sqrt(event.likelihood) << endl;
40 cout <<
"-----------------------------------" << endl<<endl;
46 if((event.rho[1]>
NETRHO) && (
event.netcc[0]>
NETCC)) {
65 char fits_name[256];
sprintf(fits_name,
"%s/probability_skymap_%d.fits.gz",
ODIR,
int(event.time[0]));
66 cout <<
"Dump fits file : " << fits_name << endl;
70 for(
int l=0;
l<skyprobcc.
size();
l++) prob+=skyprobcc.
get(
l);
71 cout <<
"prob : " << prob << endl;
detectorParams getDetectorParams()
cout<< "skymap size : "<< L<< endl;for(int l=0;l< L;l++) sm.set(l, l);sm > const_cast< char * >("skymap.dat")
double getTheta(size_t i)
size_t getSkyIndex(double th, double ph)
param: theta param: phi
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
void set(size_t i, double a)
param: sky index param: value to set
double get(size_t i)
param: sky index