9 #define JOB_STRAIN_FILE "data/strain_931158378_44_ADV_SIM_SGQ9_L1H1V1_2G_MSA_job1.root"
20 TFile*
jfile =
new TFile(jname);
21 if(jfile==NULL||!jfile->IsOpen())
22 {cout <<
"Error : file " << jname <<
" not found" << endl;
exit(1);}
27 hot[
i] = *
px;
delete px;
38 cout << hot[
i].
size() <<
" " << hot[
i].
rate() <<
" " << hot[
i].
start() << endl;
41 #if defined (DISPLAY_TIME) || defined (DISPLAY_PSD)
47 if(
n==0)
sprintf(gtitle,
"%s %s (black) ", gtitle, ifo[
n].Data());
48 if(
n==1)
sprintf(gtitle,
"%s %s (red) ", gtitle, ifo[
n].Data());
49 if(
n==2)
sprintf(gtitle,
"%s %s (green) ", gtitle, ifo[
n].Data());
53 plot.gtitle(gtitle,
"time (sec)",
"strain");
62 plot.gtitle(gtitle,
"frequency (Hz)",
"strain/#sqrt{Hz}");
66 for(
int n=0;n<nIFO;n++) hot[n] >>
plot;
virtual size_t size() const
virtual void rate(double r)
plot gtitle(gtitle,"frequency (Hz)","strain/#sqrt{Hz}")
cout<< "skymap size : "<< L<< endl;for(int l=0;l< L;l++) sm.set(l, l);sm > const_cast< char * >("skymap.dat")
void add(const wavearray< DataType_t > &, int=0, int=0, int=0)
virtual void start(double s)
wavearray< double > hot[3]
wavearray< double > x
temporary time series
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)