Logo coherent WaveBurst  
Library Reference Guide
Logo
LoopChirpMass.C
Go to the documentation of this file.
1 #define IDIR "/home/shubhanshu.tiwari/EBBH/catalogue/O1/O1_C01_SIM_eBBH_10_10_snr_20_e_0_500/output/"
2 #define OFILE "ebbh_10_10_snr_20_e_0_500.txt"
3 
4 //#define IDIR "/home/shubhanshu.tiwari/EBBH/catalogue/O1/O1_C01_SIM_eBBH_10_10_snr_20_e_6_500/output/"
5 //#define OFILE "ebbh_10_10_snr_20_e_6_500.txt"
6 
7 //#define IDIR "/home/shubhanshu.tiwari/EBBH/catalogue/O1/O1_C01_SIM_eBBH_10_10_snr_20_e_8_500/output/"
8 //#define OFILE "ebbh_10_10_snr_20_e_8_500.txt"
9 
10 {
11 /*
12  TTree itree("ebbh","ebbh")
13  itree.ReadFile("ebbh_10_10_snr_20_e_0_500.txt","m1/F:m2:s1/F:s2/F:e/F:rp/F:dist/F:redshift/F:mch/F:ech/F:ei/F")
14  itree.Draw("mch")
15  itree.Draw("ei")
16  itree.Draw("ech")
17  itree.Draw("e")
18 */
19 // TString ifile = TString::Format("%s/*_1_job4000.root",IDIR);
20 // ChirpMass(ifile);
21 
22  gROOT->LoadMacro("/home/vedovato/WP/eBBH/MACROS/ChirpMass.C");
23 
24  EBBH* ebbh = new EBBH;
25 
26  TString TAG = "";
27 
28  vector<TString> fList = CWB::Toolbox::getFileListFromDir(IDIR, ".root", "wave_", TAG, true);
29 
30  cout << "Number of files : " << fList.size() << endl;
31 
32  ofstream out;
33  out.open(OFILE);
34  for(int n=0; n<fList.size(); n++) {
35  cout << fList[n] << endl;
36  double status = ChirpMass(fList[n],0,0,0,ebbh);
37 
38  cout << n << " -> Eccentricity Index = " << ebbh->ei << " status : " << status << endl;
39 
40  out << ebbh->mass[0] << " " << ebbh->mass[0] << " "
41  << ebbh->spin[2] << " " << ebbh->spin[5] << " "
42  << ebbh->e0 << " " << ebbh->rp0 << " " << ebbh->dist << " " << ebbh->redshift << " "
43  << ebbh->ch_mass[0] << " " << ebbh->ch_energy[0] << " " << ebbh->ei << endl;
44  }
45  out.close();
46 
47  exit(0);
48 }
int ChirpMass(TString ifroot, int gtype=0, int entry=0, int ifo=0, EBBH *xebbh=NULL)
Definition: ChirpMass.C:134
static vector< TString > getFileListFromDir(TString dir_name, TString endString="", TString beginString="", TString containString="", bool fast=false)
Definition: Toolbox.cc:5108
int n
Definition: cwb_net.C:28
TString("c")
ofstream out
Definition: cwb_merge.C:214
TString TAG
Definition: LoopChirpMass.C:26
#define OFILE
Definition: LoopChirpMass.C:2
EBBH * ebbh
Definition: LoopChirpMass.C:24
vector< TString > fList
Definition: LoopChirpMass.C:28
#define IDIR
Definition: LoopChirpMass.C:1
exit(0)