21 #pragma GCC system_header 30 #include "TObjArray.h" 31 #include "TObjString.h" 60 int gISTAGE=-1;
IMPORT(
int,gISTAGE)
61 int gJSTAGE=-1;
IMPORT(
int,gJSTAGE)
64 {cout <<
"CWB_Plugin_nRMS - Error : input stage must be CSTRAIN" << endl;gSystem->Exit(1);}
67 {cout <<
"CWB_Plugin_nRMS - Error : simulation must be 0" << endl;gSystem->Exit(1);}
70 {cout <<
"CWB_Plugin_nRMS - Error : lagSize must be 1" << endl;gSystem->Exit(1);}
73 {cout <<
"CWB_Plugin_nRMS - Error : slagSize must be 1" << endl;gSystem->Exit(1);}
96 TDirectory* drms = (TDirectory*)jfile->Get(
"rms");
97 if(drms==NULL) drms=jfile->mkdir(
"rms");
99 char cdrms_name[32];
sprintf(cdrms_name,
"rms-f%d",gIFACTOR);
100 TDirectory* cdrms = (TDirectory*)drms->Get(cdrms_name);
101 if(cdrms==NULL) cdrms=drms->mkdir(cdrms_name);
104 cdrms->cd();pD->
nRMS.Write(pD->
Name);
detector * getifo(size_t n)
param: detector index
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
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
DATA_CONDITIONING.
network ** net
NOISE_MDC_SIMULATION.
#define IMPORT(TYPE, VAR)
unsigned int jobfOptions
history object
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)