25 #define DQ_LIVETIME CWB_CAT2 // livetime after the CAT1+CAT2 29 #define OUTPUT_LIVETIME_FILE_NAME "livetime.txt" 36 TB.
checkFile(gSystem->Getenv(
"CWB_ROOTLOGON_FILE"));
37 TB.
checkFile(gSystem->Getenv(
"CWB_PARAMETERS_FILE"));
38 TB.
checkFile(gSystem->Getenv(
"CWB_UPARAMETERS_FILE"));
41 cout <<
"GetLiveTime : simulation is !=0 -> only zero lag is defined!!!" << endl << endl;
65 cout<<
"lag step: "<<
lagStep<<endl;
66 cout<<
"number of time lags: "<<lags<<endl;
82 cout<<
"Number of standard jobs : " << jobList.size() <<endl<<endl;
86 if((fP = fopen(fName.Data(),
"w")) == NULL) {
87 cout <<
"cannot open output file " << fName.Data() <<
". \n";
90 cout <<
"Write output file : " << fName.Data() << endl;
91 fprintf(fP,
"# lagID\t live\t\t\n");
97 vector<double> shiftList(nIFO);
98 for(
int m=0;
m<lags;
m++) {
107 #if ( DQ_LIVETIME == CWB_CAT2 ) 108 live = TB.
getLiveTime(jobList, cat2List, shiftList);
110 #if ( DQ_LIVETIME == CWB_CAT3 ) 111 live = TB.
getLiveTime(jobList, cat3List, shiftList);
113 #if ( DQ_LIVETIME == CWB_HVETO ) 114 live = TB.
getLiveTime(jobList, hvetoList, shiftList);
118 cout <<
"-> livetime : " <<
int(live) <<
" sec " 119 << live/3600. <<
" h " << live/86400. <<
" day" << endl;
detector * getifo(size_t n)
param: detector index
#define OUTPUT_LIVETIME_FILE_NAME
std::vector< netcluster > wc_List
size_t add(detector *)
param: detector structure return number of detectors in the network
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 Import(TString umacro="")
fprintf(stdout,"start=%f duration=%f rate=%f\, x.start(), x.size()/x.rate(), x.rate())
printf("total live time: non-zero lags = %10.1f \, liveTot)
vector< waveSegment > cat1List
wavearray< double > lagShift
cout<< "total cat1 livetime : "<< int(cat1_time)<< " sec "<< cat1_time/3600.<< " h "<< cat1_time/86400.<< " day"<< endl;cout<< endl;vector< waveSegment > cat2List
detectorParams detParms[4]
int setTimeShifts(size_t=1, double=1., size_t=0, size_t=0, const char *=NULL, const char *="w", size_t *=NULL)
param number of time lags param time shift step in seconds param first lag ID param maximum lag ID pa...