3 #pragma GCC system_header
11 #include "TObjArray.h"
12 #include "TObjString.h"
26 #define SETHR 1000000 // Is the threshold which define the time slices to be cutted
29 #define TWIN 0.5 // Is the window (sec) used to integrate the energies in time
32 #define TEDG 1.5 // Is the time window (sec) to be cutted when the time integrated energy is > SETHR
40 cout <<
"-----> CWB_Plugin_Gating.C" << endl;
41 cout <<
"ifo " << ifo.Data() << endl;
42 cout <<
"type " << type << endl;
75 vector<waveSegment> glist;
82 for(
int i=N-1;
i<
size;
i++)
for(
int j=0;
j<
N;
j++) SE[
i]+=pow(hot[
n]->data[
i-
j],2);
86 int is =
i-M>X ?
i-M : X;
87 int es =
i+M<(size-X) ?
i+M : size-X;
88 for(
int k=is;
k<es;
k++) sVeto[
k]=1;
92 sVeto[0]=0;sVeto[size-1]=0;
93 for(
int i=0;
i<size-1;
i++) sVeto[
i]=sVeto[
i+1]-sVeto[
i];
103 glist.push_back(gseg);
110 for(
int j=gsize;
j<glist.size();
j++) {
111 cout <<
j <<
" -----> CWB_Plugin_Gating.C - " << net->
getifo(
n)->
Name <<
" gating time : start="
112 << glist[
j].start <<
" stop=" << glist[
j].stop <<
" len=" << glist[
j].stop-glist[
j].start << endl;
119 double gating_time = 0;
127 cout<<
"-----> CWB_Plugin_Gating.C - gating time : " << gating_time << endl << endl;
131 sprintf(info,
"-GT:%g",gating_time);
detector * getifo(size_t n)
param: detector index
virtual size_t size() const
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
COHERENCE.
virtual void rate(double r)
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
network ** net
NOISE_MDC_SIMULATION.
#define SETHR
DATA_CONDITIONING.
void PrintAnalysisInfo(CWB_STAGE stage, TString comment, TString info, bool out=true, bool log=true)
wavearray< double > hot[2]
#define IMPORT(TYPE, VAR)
wavearray< double > * getHoT()
param: no parameters
std::vector< waveSegment > segList
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)