3 #pragma GCC system_header
10 #include "TObjArray.h"
11 #include "TObjString.h"
26 cout <<
"-----> plugins/CWB_Plugin_Coherence.C" << endl;
27 cout <<
"ifo " << ifo.Data() << endl;
28 cout <<
"type " << type << endl;
37 cout <<
"type==CWB_PLUGIN_ICOHERENCE" << endl;
41 cout <<
"-----> CWB_Plugin_Sim4.C -> " <<
" gIFACTOR : " << gIFACTOR << endl;
49 int upN = rateANA/1024;
if(upN<1) upN=1;
58 int layers = level>0 ? 1<<level : 0;
59 int rate = rateANA>>level;
60 cout <<
"level : " << level <<
"\t rate(hz) : " << rate
61 <<
"\t layers : " << layers <<
"\t df(hz) : " << rateANA/2./double(1<<level)
62 <<
"\t dt(ms) : " << 1000./rate << endl;
79 cout<<
"thresholds in units of noise variance: Eo="<<Eo<<
" Emax="<<Eo*2<<endl;
82 cout<<
"live time in zero lag: "<<TL<<endl;
83 if(TL <= 0.) {cout<<
"livetime is zero : exit"<<endl;
exit(1);}
86 if(cfg->
simulation) {cout<<
"ifactor|clusters|pixels ";cout.flush();}
87 else {cout<<
"lag|clusters|pixels "; cout.flush();}
88 int csize_tot=0;
int psize_tot=0;
97 pwc->
write(jfile,
"coherence",
"clusters",0,cycle);
98 pwc->
write(jfile,
"coherence",
"clusters",-1,cycle);
99 cout<<cycle<<
"|"<<pwc->
csize()<<
"|"<<pwc->
size()<<
" ";cout.flush();
100 csize_tot+=pwc->
csize(); psize_tot+=pwc->
size();
detector * getifo(size_t n)
param: detector index
size_t write(const char *file, int app=0)
long getNetworkPixels(int LAG, double Eo, double DD=1., TH1F *hist=NULL)
cout<< "skymap size : "<< L<< endl;for(int l=0;l< L;l++) sm.set(l, l);sm > const_cast< char * >("skymap.dat")
cout<< endl;cout<< "ts size = "<< ts.size()<< " ts rate = "<< ts.rate()<< endl;tf.Forward(ts, wdm);int levels=tf.getLevel();cout<< "tf size = "<< tf.size()<< endl;double dF=tf.resolution();double dT=1./(2 *dF);cout<< "rate(hz) : "<< RATE<< "\t layers : "<< nLAYERS<< "\t dF(hz) : "<< dF<< "\t dT(ms) : "<< dT *1000.<< endl;int itime=TIME_PIXEL_INDEX;int ifreq=FREQ_PIXEL_INDEX;int index=(levels+1)*itime+ifreq;double time=itime *dT;double freq=(ifreq >0)?ifreq *dF:dF/4;cout<< endl;cout<< "PIXEL TIME = "<< time<< " sec "<< endl;cout<< "PIXEL FREQ = "<< freq<< " Hz "<< endl;cout<< endl;wavearray< double > x
size_t cluster(int kt, int kf)
param: time gap in pixels return: number of reconstructed clusters
network ** net
NOISE_MDC_SIMULATION.
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
COHERENCE.
double setVeto(double=5.)
param: time window around injections
wavearray< double > hot[2]
#define IMPORT(TYPE, VAR)
double maxEnergy(wavearray< DataType_t > &ts, Wavelet &w, double=0, int=1, int=0, TH1F *=NULL)
wavearray< double > * getHoT()
param: no parameters
double getDelay(const char *c="")
std::vector< WDM< double > * > wdmList
double THRESHOLD(double bpp)
param: selected fraction of LTF pixels assuming Gaussian noise
WSeries< double > * getTFmap()
param: no parameters
netcluster * getwc(size_t n)
param: delay index