3 #pragma GCC system_header
11 #include "TObjArray.h"
12 #include "TObjString.h"
15 #define L1_AMP_CAL_ERR_SIGMA 0.2 // 20% amplitude uncertainty
16 #define H1_AMP_CAL_ERR_SIGMA 0.2 // 20% amplitude uncertainty
17 #define V1_AMP_CAL_ERR_SIGMA 0.1 // 10% amplitude uncertainty
41 cout <<
"-----> CWB_Plugin_AmplitudeMisCal.C" << endl;
42 cout <<
"ifo " << ifo.Data() << endl;
43 cout <<
"type " << type << endl;
47 cout <<
"APPLY AMPLITUDE MIS-CALIBRATION !!!" << endl;
54 TObjArray*
token = mdcstring.Tokenize(
' ');
55 TObjString* iname = (TObjString*)token->At(11);
56 TString wavename = iname->GetString();
57 TObjString*
itime = (TObjString*)token->At(10);
58 TString wavetime = itime->GetString();
59 double mdctime = wavetime.Atof();
66 if (starti<0) starti=0;
67 if (stopi>(
int)x->
size()) stopi=x->
size();
74 for (
int jj=starti; jj<stopi; jj++) {
76 x->
data[jj] *= cfactor;
81 TString gfile=
"CWB_Plugin_AmplitudeMisCal_Plot_MDC_"+ifo+
".root";
84 cout <<
"CWB_Plugin_AmplitudeMisCal.C : created plot file name : " << gfile << endl;
virtual size_t size() const
virtual void rate(double r)
#define L1_AMP_CAL_ERR_SIGMA
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
string getmdcList(size_t n)
virtual void start(double s)
network ** net
NOISE_MDC_SIMULATION.
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
gwavearray< double > * gx
#define V1_AMP_CAL_ERR_SIGMA
#define H1_AMP_CAL_ERR_SIGMA
void Draw(GWAT_DRAW type=GWAT_TIME, TString options="ALP", Color_t color=kBlack)