6 #pragma GCC system_header
13 #include "TObjArray.h"
14 #include "TObjString.h"
32 cout <<
"Execute CWB_Plugin_MDC_Frames.C : Inject On The Fly MDC ..." << endl;
45 static bool cfg_init =
false;
48 sprintf(cmd,
"CWB::config* cfg = (CWB::config*)%p;",cfg);
49 gROOT->ProcessLine(cmd);
55 std::vector<TString> mdcPar;
58 std::vector<std::string> mdcType;
61 std::vector<std::string> channelName;
69 cout <<
"injectionList : " << injectionList << endl;
70 cout <<
"frFiles : " << frFiles << endl;
74 fr.
open(frFiles.Data());
79 cout <<
"MDC " <<
" -> nfrFiles : " << nfrFiles << endl;
83 cout <<
"mdc_range : " << mdc_range.
start <<
" " << mdc_range.
stop << endl;
88 cout <<
"mdcShift : " << mdcShift << endl;
94 vector<TString>
ifos(nIFO);
98 printf(
"GPS: %16.6f saved, injections: %d\n",
double(
long(gCWB2G->
Tb)),
i);
106 {cout <<
"CWB_Plugin_MDC_Frames.C - input rate from frame " << x->
rate()
107 <<
" do not match the one defined in config : " << cfg->
inRate << endl;gSystem->Exit(1);}
111 if(ifo.CompareTo(net->
ifoName[0])==0) {
std::vector< char * > ifoName
detector * getifo(size_t n)
param: detector index
printf("total live time: non-zero lags = %10.1f \n", liveTot)
size_t readMDClog(char *, double=0., int=11, int=12)
param: MDC log file param: approximate gps time
virtual void rate(double r)
cout<< "skymap size : "<< L<< endl;for(int l=0;l< L;l++) sm.set(l, l);sm > const_cast< char * >("skymap.dat")
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
CWB::frame fr(FRLIST_NAME)
char frFiles[NIFO_MAX+1][256]
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
COHERENCE.
std::vector< std::string > mdcType
virtual void start(double s)
std::vector< double > mdcTime
network ** net
NOISE_MDC_SIMULATION.
#define CWB_PLUGIN_EXPORT(VAR)
#define IMPORT(TYPE, VAR)
std::vector< std::string > mdcList
void setVerbose(bool verbose=true)
void setSRIndex(int srIndex)
virtual void stop(double s)
void readFrames(char *filename, char *channel, wavearray< double > &w)
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
frfile getFrList(int istart, int istop, int segEdge)
void setRetryTime(int frRetryTime=60)
void open(TString ioFile, TString chName="", Option_t *option="", bool onDisk=false, TString label=".gwf", unsigned int mode=0)