3 #pragma GCC system_header
10 #include "TObjArray.h"
11 #include "TObjString.h"
22 cout <<
"-----> CWB_Plugin_TShiftMisCal.C" << endl;
23 cout <<
"ifo " << ifo.Data() << endl;
24 cout <<
"type " << type << endl;
41 TObjArray*
token = mdcstring.Tokenize(
' ');
42 TObjString* iname = (TObjString*)token->At(11);
43 TString wavename = iname->GetString();
44 TObjString*
itime = (TObjString*)token->At(10);
45 TString wavetime = itime->GetString();
46 double mdctime = wavetime.Atof();
51 int cshift = 2*
int(gRandom->Integer(3)-1);
53 int starti = (mdctime - x->
start()-1.)*x->
rate();
54 int stopi = (mdctime - x->
start()+1.)*x->
rate();
56 if (starti<cshift) starti=cshift;
57 if (stopi>(
int)x->
size()-cshift) stopi=x->
size()-cshift;
59 if (starti<0) starti=0;
60 if (stopi>(
int)x->
size()+cshift) stopi=x->
size()+cshift;
63 cout <<
" Start: " << starti <<
" Stop : " << stopi <<
" cshift : " << cshift << endl;
64 for (
int jj=starti; jj<stopi; jj++) {
virtual void resize(unsigned int)
virtual size_t size() const
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
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
COHERENCE.
string getmdcList(size_t n)
virtual void start(double s)
network ** net
NOISE_MDC_SIMULATION.