8 #define JOB_CSTRAIN_FILE "data/cstrain_XXXX_jobYY.root" // define cstrain root file
11 #define SCRATCH_TIME 10 // is the segEdge parameter (sec)
22 TFile*
jfile =
new TFile(jname);
23 if(jfile==NULL||!jfile->IsOpen())
24 {cout <<
"Error : file " << jname <<
" not found" << endl;
exit(1);}
30 if(px==NULL) {cout<<
"Error : "<<cstrain_dir<<
" not present"<<endl;
exit(1);}
31 hot[
i] = *
px;
delete px;
34 cout << hot[
i].
size() <<
" " << hot[
i].
rate() <<
" " << hot[
i].
start() << endl;
37 int jE = hot[
i].
size()-jS;
38 cout <<
"DEB SIZE " << hot[
i].
size() <<
" " << hot[
i].
rate() <<
" " <<
int(hot[i].
start()) << endl;
41 for(
int j=jS;
j<jE;
j++) {
43 rms+=pow(hot[i].data[
j],2);
46 cout <<
"DEB RMS = " << sqrt(rms) << endl;
50 #if defined (DISPLAY_TIME) || defined (DISPLAY_PSD)
56 if(
n==0)
sprintf(gtitle,
"%s %s (black) ", gtitle, ifo[
n].Data());
57 if(
n==1)
sprintf(gtitle,
"%s %s (red) ", gtitle, ifo[
n].Data());
58 if(
n==2)
sprintf(gtitle,
"%s %s (green) ", gtitle, ifo[
n].Data());
62 plot.gtitle(gtitle,
"time (sec)",
"strain");
71 plot.gtitle(gtitle,
"frequency (Hz)",
"strain/#sqrt{Hz}");
75 for(
int n=0;n<nIFO;n++) hot[n] >>
plot;
virtual size_t size() const
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")
wavearray< double > x
temporary time series
plot gtitle(gtitle,"frequency (Hz)","strain/#sqrt{Hz}")
virtual void start(double s)
wavearray< double > hot[2]
sprintf(cstrain_dir,"cstrain/cstrain-f%d", FACTOR_INDEX)