Logo coherent WaveBurst  
Library Reference Guide
Logo
wmdc_config_rde.C
Go to the documentation of this file.
1 {
2  // --------------------------------------------------------
3  // define job parameters
4  // --------------------------------------------------------
5  TString frDir = "frames";
6  TString frLabel = "RDE_TST1";
7 
8  // define frame segment list
9  TString segmentList = "Segments/S6a_LS-segs.txt";
10 
11  // select segments
12  int jobmin = 210; // start segment
13  int jobmax = 312; // end segment
14  int jobstep = 10; // frames per job
15 
16  // --------------------------------------------------------
17  // define network
18  // --------------------------------------------------------
19  int nIFO=3;
20  TString ifo[nIFO]={"L1","H1","V1"};
21  CWB::mdc MDC(nIFO,ifo);
22 
23  // --------------------------------------------------------
24  // define waveforms
25  // --------------------------------------------------------
26  int nWF=3;
27  double F[nWF] = {1590,2090,2590};
28  double TAU[nWF] = {0.2,0.2,0.2};
29 
30  MDC.SetInjLength(2.0);
31 
32  char wf_name[256];
34  wf.par.resize(3);
35  for (int n=0;n<nWF;n++) {
36  sprintf(wf_name,"RDE_%d_%1.1f",F[n],TAU[n]);
37  wf.name = wf_name;
38  wf.name.ReplaceAll(".","d");
39  // add 20 different iota for each waveform
40  for(int m=0;m<20;m++) {
41  double cos_iota=gRandom->Uniform(-1,1);
42  double iota=acos(cos_iota)*180./TMath::Pi();
43  wf.hp = MDC.GetRD(F[n],TAU[n],iota,0);
44  wf.hx = MDC.GetRD(F[n],TAU[n],iota,1);
45  wf.hpPath = wf.name;
46  wf.hxPath = wf.name;
47  wf.par[0].name="freq"; wf.par[0].value=F[n];
48  wf.par[1].name="tau"; wf.par[1].value=TAU[n];
49  wf.par[2].name="iota"; wf.par[2].value=iota;
50  MDC.AddWaveform(wf);
51  }
52  }
53  MDC.Print();
54 
55  // --------------------------------------------------------
56  // define injection parameters
57  // --------------------------------------------------------
58  MDC.SetInjHrss(2.5e-21);
59  MDC.SetInjRate(0.01);
60  MDC.SetInjJitter(10.0);
61 
62  // --------------------------------------------------------
63  // define sky distribution
64  // --------------------------------------------------------
65  vector<mdcpar> par; par.resize(3);
66  par[0].name="entries";par[0].value=100000; // pool of events
67  par[1].name="rho_min";par[1].value=1; // min rho // Mpc
68  par[2].name="rho_max";par[2].value=100; // max rho // Mpc
70 
71  //vector<mdcpar> par; par.resize(1);
72  //par[0].name="entries";par[0].value=100000; // number of entries
73  //MDC.SetSkyDistribution(MDC_MNGD,par);
74 
75  //vector<mdcpar> par; par.resize(1);
76  //par[0].name="distance_thr";par[0].value=100; // distance max = 100 Mpc
77  //MDC.SetSkyDistribution(MDC_GWGC,"../data/GWGCCatalog_Rev1d2.txt",par);
78 
79  //vector<mdcpar> par; par.resize(3);
80  //par[0].name="theta";par[0].value=30; // theta
81  //par[1].name="phi"; par[1].value=60; // phi
82  //par[2].name="psi"; par[2].value=90; // psi
83  //MDC.SetSkyDistribution(MDC_EARTH_FIX,par);
84 
85 }
TString segmentList
double F[nWF]
Definition: mdc.hh:219
void SetInjRate(double inj_rate=MDC_INJ_RATE)
Definition: mdc.hh:311
int n
Definition: cwb_net.C:28
TString("c")
CWB::mdc MDC(nIFO, ifo)
int jobmin
void SetInjLength(double inj_length=MDC_INJ_LENGTH)
Definition: mdc.hh:304
int nWF
wavearray< double > hp
Definition: mdc.hh:226
int m
Definition: cwb_net.C:28
mdcid AddWaveform(MDC_TYPE mdc_type, vector< mdcpar > par, TString uname="")
Definition: mdc.cc:472
void Print(int level=0)
Definition: mdc.cc:2736
void SetSkyDistribution(MDC_DISTRIBUTION sky_distribution, vector< mdcpar > par, int seed=0, bool add=false)
Definition: mdc.cc:3444
TString frLabel
int jobstep
TString name
Definition: mdc.hh:222
TString hxPath
Definition: mdc.hh:224
Definition: mdc.hh:248
vector< mdcpar > par
waveform wf
double Pi
vector< mdcpar > par
Definition: mdc.hh:225
TString frDir[NIFO_MAX]
int nIFO
char wf_name[256]
void SetInjJitter(double inj_jitter=MDC_INJ_JITTER)
Definition: mdc.hh:317
double e
wavearray< double > hx
Definition: mdc.hh:227
int jobmax
TString ifo[nIFO]
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
void SetInjHrss(double inj_hrss=MDC_INJ_HRSS)
Definition: mdc.hh:307
wavearray< double > GetRD(double frequency, double tau, double iota, bool polarization=0)
Definition: mdc.cc:3168
TString hpPath
Definition: mdc.hh:223
double TAU[nWF]