25 jfile = gROOT->GetFile();
26 jfile->ReOpen(
"UPDATE");
27 if(jfile==NULL) {cout <<
"Error opening input root file " << endl;
exit(1);}
38 else cdprod = (TDirectory*)jfile->Get(
"prod");
45 TObjString* path_tok = (TObjString*)token->At(token->GetEntries()-1);
47 cout << jname.Data() << endl;
50 cout <<
"nIFO : " <<
nIFO << endl;
79 cout<<
"maximum time delay between detectors: "<<mTau<<endl;
80 cout<<
" maximum time delay difference: "<<dTau<<endl;
81 cout<<
" netRHO and netCC: "<<NET.
netRHO<<
", "<<NET.
netCC<<endl;
82 cout<<
" regulator delta, gamma: "<<NET.
delta <<
" " <<NET.
gamma<<endl<<endl;
84 gSystem->Exec(
"/bin/date");
85 gSystem->Exec(
"/bin/hostname");
92 for(
int iii=0; iii<
nfactor; iii++) {
94 cout<<
"factor="<<factor<<endl;
96 sprintf(sfactor,
"factor_%g",factor);
104 if(pTF[
i]==NULL) {cout <<
"Error : data not present, job terminated!!!" << endl;
return;}
115 pD[
i]->
white(60.,1,8.,20.);
122 cout<<
"After "<<pD[
i]->
Name<<
" data conditioning"<<endl;
139 for(
int j=0;
j<lags;
j++){
142 TIter nextkey(((TDirectory*)jfile->Get(wcdir))->GetListOfKeys());
143 TKey *key;
while(key=(TKey*)nextkey()) wc.
append(*(
netcluster*)key->ReadObj());
146 cout<<m<<
"|"<<pwc->
size()<<
" ";
153 TDirectory* cdced = NULL;
169 cout<<
"rejected weak pixels: "<<NET.
netcut(NET.
netRHO,
'r',0,1)<<
"\n";
170 cout<<
"rejected loud pixels: "<<NET.
netcut(NET.
netCC,
'c',0,1)<<
"\n";
171 cout<<
"events in the buffer: "<<NET.
events()<<
"\n";
175 cout<<
"dump ced into "<<end_CED<<
"\n";
176 CWB::ced ced(&NET,&netburst,end_CED);
191 cout<<
"\nSearch done\n";
192 cout<<
"reconstructed events: "<<NET.
events()<<
"\n";
198 if(
simulation) cdsim->cd();
else cdprod->cd();
200 TTree*
net_tree = netburst.setTree();
205 int rnID =
int(gRandom->Rndm(13)*1.e9);
206 UserGroup_t*
uinfo = gSystem->GetUserInfo();
208 gSystem->Exec(
TString(
"mkdir -p /dev/shm/")+uname);
209 sprintf(ofile,
"/dev/shm/%s/waveburst-%d.txt",uname.Data(),
rnID);
213 if(dump) netburst.dopen(ofile,
"w");
214 netburst.output(net_tree,&NET);
216 if(dump) netburst.dclose();
219 macro.Write(
"output.txt");
220 gSystem->Exec(
TString(
"rm ")+ofile);
223 jfile->ReOpen(
"READ");
224 gROOT->RefreshBrowsers();
225 gROOT->ForceStyle(0);
std::vector< char * > ifoName
detector * getifo(size_t n)
param: detector index
WSeries< double > * pTF[NIFO_MAX]
void white(double dT=0, int wtype=1, double offset=0., double stride=0.)
what it does: see algorithm description in wseries.hh
void Export(TString fname="")
void setAntenna(detector *)
param: detector (use theta, phi index array)
std::vector< double > * getmdcTime()
size_t setIndexMode(size_t=0)
void add(const wavearray< DataType_t > &, int=0, int=0, int=0)
virtual size_t supercluster(char atype, double S, bool core)
param: statistic: E - excess power, L - likelihood param: selection threshold T for likelihood cluste...
void Forward(size_t k)
param: number of steps
void setDelayIndex(double rate)
param: MDC log file
double getDelay(const char *c="")
virtual size_t append(netcluster &wc)
param: input netcluster return size of appended pixel list
size_t cpf(const netcluster &, bool=false, int=0)
void setDelay(const char *="L1")
void bandPass(double f1, double f2, double a=0.)
long likelihood(char='E', double=sqrt(2.), int=0, size_t=0, int=-1, bool=false)
WSeries< double > * getTFmap()
param: no parameters
virtual void lprFilter(double, int=0, double=0., double=0.)
void Forward(int n=-1)
param: wavelet - n is number of steps (-1 means full decomposition)
netcluster * getwc(size_t n)
param: delay index
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
void setDelayFilters(detector *=NULL)
size_t netcut(double, char='L', size_t=0, int=1)
param: threshold param: minimum cluster size processed by the corrcut param: cluster type return: num...
int setTimeShifts(size_t=1, double=1., size_t=0, size_t=0, const char *=NULL, const char *="w", size_t *=NULL)
param number of time lags param time shift step in seconds param first lag ID param maximum lag ID pa...
void Inverse(int n=-1)
param: n - number of steps (-1 means full reconstruction)
size_t setsim(WSeries< double > &, std::vector< double > *, double=5., double=8., bool saveWF=false)
void setSkyMaps(double, double=0., double=180., double=0., double=360.)
param: sky map granularity step, degrees param: theta begin, degrees param: theta end...