Logo coherent WaveBurst  
Library Reference Guide 		Logo
Functions | Variables
cwb_net.C File Reference
#include <vector>
Include dependency graph for cwb_net.C:

Functions

history AddHistory (job_stage, "WATVERSION", watversion('s'))
 
history AddHistory (job_stage, "WORKDIR", work_dir)
 
history AddHistory (job_stage, "DATALABEL", data_label)
 
history AddHistory (job_stage, "CMDLINE", cmd_line)
 
history AddLog (job_stage, "START JOB")
 
s gApplication Argv (i))
 
cout<<" network of ";for(i=0;i< nIFO;i++) cout<< ifo[i]<<" ";cout<<" detectors\\";Meyer< double > B (1024)
 
TB checkFile (gSystem->Getenv("CWB_ROOTLOGON_FILE"))
 
TB checkFile (gSystem->Getenv("CWB_PARAMETERS_FILE"))
 
TB checkFile (gSystem->Getenv("CWB_UPARAMETERS_FILE"))
 
pTF [i] DumpBinary (file)
 
cout<<"Stopping the job "<< runID<< endl;gSystem-> Exec ("/bin/date")
 
 for ()
 pointers to detectors More...
 
 for (int i=0;i< gApplication->Argc();i++) sprintf(cmd_line
 
pTF [i] Forward (x, B, levelR)
 
pTF [i] Forward (x, S, levelD)
 
 fprintf (stdout,"start=%f duration=%f rate=%f\, x.start(), x.size()/x.rate(), x.rate())
 
pTF [i] getLayer (x, 0)
 
cout<< "live time after cat 2 : "<< detSegs_ctime<< endl;if(detSegs_ctime< segTHR) {cout<< "job segment live time after cat2 < "<< segTHR<< " sec, job terminated !!!"<< endl;exit(1);} double Tb=detSegs[0].start;double Te=detSegs[0].stop;double dT=Te-Tb;char file[512], tdf00[512], tdf90[512], buFFer[1024];int rnID=int(gRandom->Rndm(13) *1.e9);if(simulation) { i=NET.readMDClog(injectionList, double(long(Tb)) -mdcShift);printf("GPS: %16.6f saved, injections: %d\", double(long(Tb)), i);frTB[nIFO].shiftBurstMDCLog(NET.mdcList, ifos, mdcShift);for(int i=0;i< NET.mdcTime.size();i++) NET.mdcTime[i]+=mdcShift;vector< waveSegment > mdcSegs(NET.mdcTime.size());for(int k=0;k< NET.mdcTime.size();k++) {mdcSegs[k].start=NET.mdcTime[k]-gap;mdcSegs[k].stop=NET.mdcTime[k]+gap;} vector< waveSegment > mdcSegs_dq2=slagTB.mergeSegLists(detSegs_dq2, mdcSegs);double mdcSegs_ctime=slagTB.getTimeSegList(mdcSegs_dq2);cout<< "live time in zero lag after cat2+inj : "<< mdcSegs_ctime<< endl;if(mdcSegs_ctime==0) {cout<< "job segment with zero cat2+inj live time in zero lag, job terminated !!!"<< endl;exit(1);} } if(dump_infos_and_exit) exit(0);if(mask >0.) NET.setSkyMask(mask, skyMaskFile);for(i=0;i< nIFO;i++) { frTB[i].readFrames(FRF[i], channelNamesRaw[i], x);x.start(x.start()+dataShift[i]);x.start(x.start() -segLen *(segID[i]-segID[0]));if(singleDetector) TB.resampleToPowerOfTwo(x);sprintf(file,"%s/%s_%d_%s_%d_%d.dat", nodedir, ifo[i], int(Tb), data_label, runID, rnID);if(dump_sensitivity_and_exit) { sprintf(file,"%s/sensitivity_%s_%d_%s_job%d.txt", dump_dir, ifo[i], int(Tb), data_label, runID);cout<< endl<< "Dump Sensitivity : "<< file<< endl<< endl;TB.makeSpectrum(file, x);continue;} if(dcCal[i]>0.) x *=dcCal[i];if(fResample >0) { x.FFT(1);x.resize(fResample/x.rate() *x.size());x.FFT(-1);x.rate(fResample);} pTF[i]=pD[i]-> getTFmap ()
 
 if (if(=='i'||search=='I') cout<<"\lliptical polarisation: "<< SEARCH()<< endl search=='E'||search=='E')
 
 if (rootlogonBuffer!=NULL)
 
 if (slagSize >0)
 
 if (detSegs.size()==0)
 
 if (simulation)
 
 if (i >0 &&pTF[0]->start() !=x.start())
 
 if (i >0 &&pTF[0]->rate() !=x.rate())
 
 if (singleDetector)
 
vector< TString > ifos (nIFO)
 
detSegs_dq2 push_back (detSegs[0])
 

Variables

vector< waveSegment > cat1List
 
vector< waveSegment > cat2List =slagTB.readSegList(nDQF, DQF, CWB_CAT2)
 
char cmd_line [512] =""
 
TString cwb_parameters_name = TString(gSystem->Getenv("CWB_PARAMETERS_FILE"))
 
vector< waveSegment > detSegs
 
double detSegs_ctime = slagTB.getTimeSegList(detSegs_dq2)
 
vector< waveSegment > detSegs_dq2 = slagTB.mergeSegLists(detSegs_dq2,cat2List)
 
 else
 
frfile FRF [nIFO+1]
 
CWB::Toolbox frTB [nIFO+1]
 
int i
 
int j
 
int job_data_size_sec = int(detSegs[0].stop-detSegs[0].start)
 
int job_elapsed_hour = int(job_elapsed_time/3600)
 
int job_elapsed_min = int((job_elapsed_time-3600*job_elapsed_hour)/60)
 
int job_elapsed_sec = int(job_elapsed_time-3600*job_elapsed_hour-60*job_elapsed_min)
 
int job_elapsed_time = (job_stop_time-job_start_time)
 
double job_speed_factor = double(job_data_size_sec)/double(job_elapsed_time)
 
int job_start_time =CWB::Time("now").GetGPS()
 
int job_stop_time =CWB::Time("now").GetGPS()
 
int jobID
 
int m
 
double mdcShift =0.
 
int n = pTF[i]->size()
 
int nfrFiles [nIFO+1]
 
detector * pD [nIFO]
 
WSeries< double > * pTF [nIFO]
 
netcluster * pwc
 
char * rootlogonBuffer = histTB.readFile("rootlogon.C")
 
Meyer< double > S (1024, 2)
 
s s
 
int segID [20]
 
NET segList =cat2List
 
int slagID = 0
 
CWB::Toolbox slagTB
 
WSeries< float > v [nIFO]
 
netcluster wc
 
WSeries< double > wM
 
wavearray< double > x
 

Function Documentation

◆ AddHistory() [1/4]

history AddHistory ( job_stage  ,
"WATVERSION"  ,
watversion( 's')   
)

◆ AddHistory() [2/4]

history AddHistory ( job_stage  ,
"WORKDIR"  ,
work_dir   
)

◆ AddHistory() [3/4]

history AddHistory ( job_stage  ,
"DATALABEL"  ,
data_label   
)

◆ AddHistory() [4/4]

history AddHistory ( job_stage  ,
"CMDLINE"  ,
cmd_line   
)

◆ AddLog()

history AddLog ( job_stage  ,
"START JOB"   
)

◆ Argv()

s gApplication Argv ( i  )

◆ B()

cout<<" network of "; for(i=0; i<nIFO; i++) cout<<ifo[i]<<" "; cout<<" detectors\n\n"; Meyer<double> B ( 1024  )

◆ checkFile() [1/3]

TB checkFile ( gSystem->  Getenv"CWB_ROOTLOGON_FILE")

◆ checkFile() [2/3]

TB checkFile ( gSystem->  Getenv"CWB_PARAMETERS_FILE")

◆ checkFile() [3/3]

TB checkFile ( gSystem->  Getenv"CWB_UPARAMETERS_FILE")

◆ DumpBinary()

pTF [i] DumpBinary ( file  )

◆ Exec()

cout<<"Stopping the job "<<runID<<endl; gSystem-> Exec ( "/bin/date"  )

◆ for() [1/2]

for ( )

pointers to detectors

Definition at line 84 of file cwb_net.C.

◆ for() [2/2]

for ( int  i = 0;i<gApplication->Argc();i++)

Definition at line 166 of file cwb_net.C.

◆ Forward() [1/2]

pTF [i] Forward ( x  ,
B  ,
levelR   
)

◆ Forward() [2/2]

pTF [i] Forward ( x  ,
S  ,
levelD   
)

◆ fprintf()

fprintf ( stdout  ,
"  start = %f duration=%f rate=%f\,
x.  start(),
x.  size)/x.rate(,
x.  rate() 
)

◆ getLayer()

pTF [i] getLayer ( x  ,
0   
)

◆ getTFmap()

cout<< "live time after cat 2 : " << detSegs_ctime << endl; if(detSegs_ctime<segTHR) {cout << "job segment live time after cat2 < " << segTHR << " sec, job terminated !!!" << endl;exit(1);} double Tb=detSegs[0].start; double Te=detSegs[0].stop; double dT = Te-Tb; char file[512], tdf00[512], tdf90[512], buFFer[1024]; int rnID = int(gRandom->Rndm(13)*1.e9); if(simulation) { i=NET.readMDClog(injectionList,double(long(Tb))-mdcShift); printf("GPS: %16.6f saved, injections: %d\n",double(long(Tb)),i); frTB[nIFO].shiftBurstMDCLog(NET.mdcList, ifos, mdcShift); for(int i=0;i<NET.mdcTime.size();i++) NET.mdcTime[i]+=mdcShift; vector<waveSegment> mdcSegs(NET.mdcTime.size()); for(int k=0;k<NET.mdcTime.size();k++) {mdcSegs[k].start=NET.mdcTime[k]-gap;mdcSegs[k].stop=NET.mdcTime[k]+gap;} vector<waveSegment> mdcSegs_dq2 = slagTB.mergeSegLists(detSegs_dq2,mdcSegs); double mdcSegs_ctime = slagTB.getTimeSegList(mdcSegs_dq2); cout << "live time in zero lag after cat2+inj : " << mdcSegs_ctime << endl; if(mdcSegs_ctime==0) {cout << "job segment with zero cat2+inj live time in zero lag, job terminated !!!" << endl;exit(1);} } if(dump_infos_and_exit) exit(0); if(mask>0.) NET.setSkyMask(mask,skyMaskFile); for(i=0; i<nIFO; i++) { frTB[i].readFrames(FRF[i],channelNamesRaw[i],x); x.start(x.start()+dataShift[i]); x.start(x.start()-segLen*(segID[i]-segID[0])); if(singleDetector) TB.resampleToPowerOfTwo(x); sprintf(file,"%s/%s_%d_%s_%d_%d.dat", nodedir,ifo[i],int(Tb),data_label,runID,rnID); if(dump_sensitivity_and_exit) { sprintf(file,"%s/sensitivity_%s_%d_%s_job%d.txt",dump_dir,ifo[i],int(Tb),data_label,runID); cout << endl << "Dump Sensitivity : " << file << endl << endl; TB.makeSpectrum(file, x); continue; } if(dcCal[i]>0.) x*=dcCal[i]; if(fResample>0) { x.FFT(1); x.resize(fResample/x.rate()*x.size()); x.FFT(-1); x.rate(fResample); } pTF[i] = pD[i]-> getTFmap ( )

◆ if() [1/8]

if ( if(=='i' || search=='I') cout<<"\n elliptical polarisation: "<<SEARCH()<<endl  search = ='E' || search=='E')

Definition at line 30 of file cwb_net.C.

◆ if() [2/8]

if ( rootlogonBuffer!  = NULL)

Definition at line 159 of file cwb_net.C.

◆ if() [3/8]

if ( slagSize  ,
0   
)

Definition at line 204 of file cwb_net.C.

◆ if() [4/8]

if ( detSegs.  size() = =0)

Definition at line 266 of file cwb_net.C.

◆ if() [5/8]

if ( simulation  )

Definition at line 275 of file cwb_net.C.

◆ if() [6/8]

if ( i  ,
0 &&pTF->start() !  [0] = x.start() 
)

Definition at line 361 of file cwb_net.C.

◆ if() [7/8]

if ( i  ,
0 &&pTF->rate() !  [0] = x.rate() 
)

Definition at line 366 of file cwb_net.C.

◆ if() [8/8]

if ( singleDetector  )

Definition at line 393 of file cwb_net.C.

◆ ifos()

vector<TString> ifos ( nIFO  )

◆ push_back()

detSegs_dq2 push_back ( detSegs  [0])

Variable Documentation

◆ cat1List

vector<waveSegment> cat1List

Definition at line 194 of file cwb_net.C.

◆ cat2List

cat2List =slagTB.readSegList(nDQF, DQF, CWB_CAT2)

Definition at line 194 of file cwb_net.C.

◆ cmd_line

s cmd_line =""

Definition at line 154 of file cwb_net.C.

◆ cwb_parameters_name

TString cwb_parameters_name = TString(gSystem->Getenv("CWB_PARAMETERS_FILE"))

Definition at line 165 of file cwb_net.C.

◆ detSegs

vector<waveSegment> detSegs

Definition at line 193 of file cwb_net.C.

◆ detSegs_ctime

double detSegs_ctime = slagTB.getTimeSegList(detSegs_dq2)

Definition at line 303 of file cwb_net.C.

◆ detSegs_dq2

detSegs_dq2 = slagTB.mergeSegLists(detSegs_dq2,cat2List)

Definition at line 297 of file cwb_net.C.

◆ else

else
Initial value:
{
jobID = runID
jobID
int jobID
Definition: cwb_net.C:195

Definition at line 241 of file cwb_net.C.

◆ FRF

frfile FRF[nIFO+1]

Definition at line 269 of file cwb_net.C.

◆ frTB

CWB::Toolbox frTB[nIFO+1]

Definition at line 191 of file cwb_net.C.

◆ i

int i

Definition at line 28 of file cwb_net.C.

◆ j

int j

Definition at line 28 of file cwb_net.C.

◆ job_data_size_sec

int job_data_size_sec = int(detSegs[0].stop-detSegs[0].start)

Definition at line 740 of file cwb_net.C.

◆ job_elapsed_hour

int job_elapsed_hour = int(job_elapsed_time/3600)

Definition at line 737 of file cwb_net.C.

◆ job_elapsed_min

int job_elapsed_min = int((job_elapsed_time-3600*job_elapsed_hour)/60)

Definition at line 738 of file cwb_net.C.

◆ job_elapsed_sec

int job_elapsed_sec = int(job_elapsed_time-3600*job_elapsed_hour-60*job_elapsed_min)

Definition at line 739 of file cwb_net.C.

◆ job_elapsed_time

int job_elapsed_time = (job_stop_time-job_start_time)

Definition at line 736 of file cwb_net.C.

◆ job_speed_factor

double job_speed_factor = double(job_data_size_sec)/double(job_elapsed_time)

Definition at line 741 of file cwb_net.C.

◆ job_start_time

int job_start_time =CWB::Time("now").GetGPS()

Definition at line 181 of file cwb_net.C.

◆ job_stop_time

int job_stop_time =CWB::Time("now").GetGPS()

Definition at line 734 of file cwb_net.C.

◆ jobID

int jobID

Definition at line 195 of file cwb_net.C.

◆ m

int m

Definition at line 28 of file cwb_net.C.

◆ mdcShift

double mdcShift =0.

Definition at line 196 of file cwb_net.C.

◆ n

n = pTF[i]->size()

Definition at line 28 of file cwb_net.C.

◆ nfrFiles

int nfrFiles[nIFO+1]

Definition at line 192 of file cwb_net.C.

◆ pD

detector* pD[nIFO]

Definition at line 81 of file cwb_net.C.

◆ pTF

WSeries<double>* pTF[nIFO]

Definition at line 82 of file cwb_net.C.

◆ pwc

netcluster* pwc

Definition at line 76 of file cwb_net.C.

◆ rootlogonBuffer

char* rootlogonBuffer = histTB.readFile("rootlogon.C")

Definition at line 158 of file cwb_net.C.

◆ S

Meyer<double> S(1024, 2)

◆ s

s s

Definition at line 155 of file cwb_net.C.

◆ segID

int segID[20]

Definition at line 195 of file cwb_net.C.

◆ segList

NET segList =cat2List

Definition at line 294 of file cwb_net.C.

◆ slagID

slagID = 0

Definition at line 195 of file cwb_net.C.

◆ slagTB

CWB::Toolbox slagTB

Definition at line 190 of file cwb_net.C.

◆ v

WSeries<float> v[nIFO]

Definition at line 80 of file cwb_net.C.

◆ wc

netcluster wc

Definition at line 75 of file cwb_net.C.

◆ wM

WSeries<double> wM

Definition at line 79 of file cwb_net.C.

◆ x

wavearray<double> x

Definition at line 78 of file cwb_net.C.