Logo Coherent WaveBurst  
Reference Guide
Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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\n\n";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[iDumpBinary (file)
 
cout<<"Stopping the job "
<< runID<< endl;gSystem-> 
Exec ("/bin/date")
 
 for ()
 
 for (int i=0;i< gApplication->Argc();i++) sprintf(cmd_line
 
pTF[iForward (x, B, levelR)
 
pTF[iForward (x, S, levelD)
 
 fprintf (stdout,"start=%f duration=%f rate=%f\n", x.start(), x.size()/x.rate(), x.rate())
 
pTF[igetLayer (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\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 (if(=='i'||search=='I') cout<<"\n elliptical 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< TStringifos (nIFO)
 
detSegs_dq2 push_back (detSegs[0])
 

Variables

vector< waveSegmentcat1List
 
vector< waveSegmentcat2List =slagTB.readSegList(nDQF, DQF, CWB_CAT2)
 
char cmd_line [512] =""
 
TString cwb_parameters_name = TString(gSystem->Getenv("CWB_PARAMETERS_FILE"))
 
vector< waveSegmentdetSegs
 
double detSegs_ctime = slagTB.getTimeSegList(detSegs_dq2)
 
vector< waveSegmentdetSegs_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]
 
detectorpD [nIFO]
 
WSeries< double > * pTF [nIFO]
 
netclusterpwc
 
char * rootlogonBuffer = histTB.readFile("rootlogon.C")
 
Meyer< double > S (1024, 2)
 
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

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\n\n"; 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 ( )

Definition at line 66 of file cwb_net.C.

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

Definition at line 148 of file cwb_net.C.

pTF [i] Forward ( x  ,
B  ,
levelR   
)
pTF [i] Forward ( x  ,
S  ,
levelD   
)
fprintf ( stdout  )
pTF [i] getLayer ( x  ,
 
)
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 ( if(=='i' || search=='I') cout<<"\n elliptical polarisation: "<<SEARCH()<<endl  search = ='E' || search=='E')

Definition at line 12 of file cwb_net.C.

if ( rootlogonBuffer!  = NULL)

Definition at line 141 of file cwb_net.C.

if ( slagSize  ,
 
)

Definition at line 186 of file cwb_net.C.

if ( detSegs.  size() = =0)

Definition at line 248 of file cwb_net.C.

if ( simulation  )

Definition at line 257 of file cwb_net.C.

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

Definition at line 343 of file cwb_net.C.

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

Definition at line 348 of file cwb_net.C.

if ( singleDetector  )

Definition at line 375 of file cwb_net.C.

vector<TString> ifos ( nIFO  )
detSegs_dq2 push_back ( detSegs  [0])

Variable Documentation

vector<waveSegment> cat1List

Definition at line 176 of file cwb_net.C.

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

Definition at line 176 of file cwb_net.C.

s cmd_line =""

Definition at line 136 of file cwb_net.C.

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

Definition at line 147 of file cwb_net.C.

vector<waveSegment> detSegs

Definition at line 175 of file cwb_net.C.

double detSegs_ctime = slagTB.getTimeSegList(detSegs_dq2)

Definition at line 285 of file cwb_net.C.

detSegs_dq2 = slagTB.mergeSegLists(detSegs_dq2,cat2List)

Definition at line 279 of file cwb_net.C.

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

Definition at line 223 of file cwb_net.C.

frfile FRF[nIFO+1]

Definition at line 251 of file cwb_net.C.

CWB::Toolbox frTB[nIFO+1]

Definition at line 173 of file cwb_net.C.

int i

Definition at line 10 of file cwb_net.C.

int j

Definition at line 10 of file cwb_net.C.

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

Definition at line 722 of file cwb_net.C.

int job_elapsed_hour = int(job_elapsed_time/3600)

Definition at line 719 of file cwb_net.C.

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

Definition at line 720 of file cwb_net.C.

Definition at line 721 of file cwb_net.C.

int job_elapsed_time = (job_stop_time-job_start_time)

Definition at line 718 of file cwb_net.C.

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

Definition at line 723 of file cwb_net.C.

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

Definition at line 163 of file cwb_net.C.

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

Definition at line 716 of file cwb_net.C.

int jobID

Definition at line 177 of file cwb_net.C.

int m

Definition at line 10 of file cwb_net.C.

double mdcShift =0.

Definition at line 178 of file cwb_net.C.

n = pTF[i]->size()

Definition at line 10 of file cwb_net.C.

int nfrFiles[nIFO+1]

Definition at line 174 of file cwb_net.C.

detector* pD[nIFO]

Definition at line 63 of file cwb_net.C.

WSeries<double>* pTF[nIFO]

Definition at line 64 of file cwb_net.C.

netcluster* pwc

Definition at line 58 of file cwb_net.C.

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

Definition at line 140 of file cwb_net.C.

Meyer<double> S(1024, 2)
s s

Definition at line 137 of file cwb_net.C.

int segID[20]

Definition at line 177 of file cwb_net.C.

NET segList =cat2List

Definition at line 276 of file cwb_net.C.

slagID = 0

Definition at line 177 of file cwb_net.C.

CWB::Toolbox slagTB

Definition at line 172 of file cwb_net.C.

WSeries<float> v[nIFO]

Definition at line 62 of file cwb_net.C.

Definition at line 57 of file cwb_net.C.

WSeries<double> wM

Definition at line 61 of file cwb_net.C.

wavearray<double> x

Definition at line 60 of file cwb_net.C.