21 #pragma GCC system_header 29 #include "TObjArray.h" 30 #include "TObjString.h" 42 cout <<
"-----> CWB_Plugin_supercluster.C" << endl;
43 cout <<
"ifo " << ifo.Data() << endl;
44 cout <<
"type " << type << endl;
53 cout <<
"type==CWB_PLUGIN_ISUPERCLUSTER" << endl;
56 void* gIFILE=NULL;
IMPORT(
void*,gIFILE)
57 cout <<
"-----> CWB_Plugin_wavegraph.C -> " <<
" gIFILE : " << gIFILE << endl;
58 TFile*
ifile = (TFile*)gIFILE;
62 cout <<
"-----> CWB_Plugin_wavegraph.C -> " <<
" gIFACTOR : " << gIFACTOR << endl;
80 for(
int j=0;
j<(
int)lags;
j++) {
85 if(ifile!=NULL) wc.
read(ifile,
"coherence",
"clusters",0,cycle);
86 else wc.
read(jfile,
"coherence",
"clusters",0,cycle);
89 for(
int i=nRES-1;
i>=0;
i--)
90 wc.
read(ifile,
"coherence",
"clusters",-1,cycle,rateANA>>(
i+cfg->
l_low));
92 for(
int i=nRES-1;
i>=0;
i--)
93 wc.
read(jfile,
"coherence",
"clusters",-1,cycle,rateANA>>(
i+cfg->
l_low));
95 if(!cfg->
simulation) cout<<
"process lag " <<cycle<<
" ..."<<endl;
96 cout<<
"loaded clusters|pixels: "<<wc.
csize()<<
"|"<<wc.
size()<<endl;
100 cout<<
"super clusters|pixels: "<<wc.
esize(0)<<
"|"<<wc.
psize(0)<<endl;
115 double pfrac = ptot>0 ? double(psel)/double(ptot) : 0.;
116 cout<<
"selected pixels: "<<psel<<
", fraction: "<<pfrac<< endl;
117 if(count<10000)
break;
123 nnn += pwc->
psize(-1);
126 if(mmm) cout<<
"events in the buffer: "<<net->
events()<<
"|"<<nnn<<
"|"<<nnn/double(mmm)<<
"\n";
127 else cout<<
"events in the buffer: "<<net->
events()<<
"\n";
130 pwc->
write(jfile,
"supercluster",
"clusters",0,cycle);
131 pwc->
write(jfile,
"supercluster",
"clusters",-1,cycle);
132 cout<<cycle<<
"|"<<pwc->
csize()<<
"|"<<pwc->
size()<<
" ";cout.flush();
detector * getifo(size_t n)
param: detector index
virtual size_t defragment(double T, double F, TH2F *=NULL)
T - maximum time gap in seconds F - maximum frequency gap in Hz.
size_t write(const char *file, int app=0)
std::vector< netcluster > wc_List
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<< "\ layers : "<< nLAYERS<< "\ dF(hz) : "<< dF<< "\ 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 setcore(bool core, int id=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...
long subNetCut(int lag, float subnet=0.6, float subcut=0.33, TH2F *hist=NULL)
network ** net
NOISE_MDC_SIMULATION.
virtual size_t size() const
wavearray< double > hot[2]
#define IMPORT(TYPE, VAR)
wavearray< double > * getHoT()
param: no parameters
void setDelayIndex(double rate)
param: MDC log file
size_t cpf(const netcluster &, bool=false, int=0)
size_t loadTDampSSE(network &net, char c, size_t BATCH=10000, size_t LOUD=0)
wavearray< double > lagShift
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
SUPERCLUSTER.
netcluster * getwc(size_t n)
param: delay index
size_t read(const char *)