Logo Coherent WaveBurst  
Reference Guide 		Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
CWB_Plugin_Sim4_BRST_Config.C
Go to the documentation of this file.
1 {
2  //!NOISE_MDC_SIMULATION
3  // Config Plugin to generate simulated MDC injected 'on the fly' using simulation=4
4 
5 
6  cout << "-----> CWB_Plugin_Sim4_BRST_Config.C -> " << " run : " << net->nRun << " gIFACTOR " << gIFACTOR << endl;
7 
8  network** net;
9  CWB_PLUGIN_IMPORT(network**,net);
10 
11  int seed = 100*(*net)->nRun+gIFACTOR;
12 
13  CWB::mdc* MDC;
14  CWB_PLUGIN_IMPORT(CWB::mdc*,MDC);
15 
16  char wf_name[256];
17  waveform wf;
18  vector<mdcpar> par;
19 
20  par.resize(1);
21 
22  double F_GA[4] = {0.0001,0.001,0.0025,0.004};
23  for(int n=0;n<4;n++) {
24  par[0].name="duration"; par[0].value=F_GA[n];
25  MDC->AddWaveform(MDC_GA, par);
26  }
27 
28  par.resize(2);
29 
30  double F_Q3[4] = {70,235,849,1615};
31  for(int n=0;n<4;n++) {
32  par[0].name="frequency"; par[0].value=F_Q3[n];
33  par[1].name="Q"; par[1].value=3.;
34  MDC->AddWaveform(MDC_SG, par);
35  }
36 
37  double F_Q8d9[7] = {70,100,153,235,361,554,849};
38  for(int n=0;n<7;n++) {
39  par[0].name="frequency"; par[0].value=F_Q8d9[n];
40  par[1].name="Q"; par[1].value=8.9;
41  MDC->AddWaveform(MDC_SG, par);
42  }
43 
44  double F_Q9[4] = {1053,1304,1615,2000};
45  for(int n=0;n<4;n++) {
46  par[0].name="frequency"; par[0].value=F_Q9[n];
47  par[1].name="Q"; par[1].value=9.;
48  MDC->AddWaveform(MDC_SG, par);
49  }
50 
51  double F_Q100[4] = {70,235,849,1615};
52  for(int n=0;n<4;n++) {
53  par[0].name="frequency"; par[0].value=F_Q100[n];
54  par[1].name="Q"; par[1].value=100.;
55  MDC->AddWaveform(MDC_SG, par);
56  }
57  par.resize(6);
58 
59  double F_WNB[5] = {100,250,1000,1000,1000};
60  double B_WNB[5] = {100,100,10,1000,1000};
61  double D_WNB[5] = {0.1,0.1,0.1,0.1,0.01};
62  for(int n=0;n<5;n++) {
63  par[0].name="frequency"; par[0].value=F_WNB[n];
64  par[1].name="bandwidth"; par[1].value=B_WNB[n];
65  par[2].name="duration"; par[2].value=D_WNB[n];
66  for(int m=0;m<30;m++) {
67  par[3].name="pseed"; par[3].value=seed+100000+n*100+m;
68  par[4].name="xseed"; par[4].value=seed+100001+n*100+m;
69  par[5].name="mode"; par[5].value=1; // simmetric
70  MDC->AddWaveform(MDC_WNB, par);
71  }
72  }
73 
74  MDC->Print();
75 
76  // --------------------------------------------------------
77  // define injection parameters
78  // --------------------------------------------------------
79  MDC->SetInjHrss(2.5e-21);
80  MDC->SetInjRate(0.01);
81  MDC->SetInjJitter(10.0);
82  MDC->SetSourceListSeed(gIFACTOR);
83 
84  // --------------------------------------------------------
85  // define the MDC hrss
86  // --------------------------------------------------------
87  int K=MDC.wfList.size();
88  std::vector<double> mdcHRSS(K);
89  for(int i=0;i<K;i++) {
90  if(TString(MDC->wfList[i].name)=="GA0d1") mdcHRSS[i]=8.980e-22;
91  if(TString(MDC->wfList[i].name)=="GA1d0") mdcHRSS[i]=5.097e-22;
92  if(TString(MDC->wfList[i].name)=="GA2d5") mdcHRSS[i]=9.252e-22;
93  if(TString(MDC->wfList[i].name)=="GA4d0") mdcHRSS[i]=1.440e-21;
94  if(TString(MDC->wfList[i].name)=="SG70Q3") mdcHRSS[i]=1.041e-21;
95  if(TString(MDC->wfList[i].name)=="SG235Q3") mdcHRSS[i]=3.180e-22;
96  if(TString(MDC->wfList[i].name)=="SG849Q3") mdcHRSS[i]=9.420e-22;
97  if(TString(MDC->wfList[i].name)=="SG1615Q3") mdcHRSS[i]=1.747e-21;
98  if(TString(MDC->wfList[i].name)=="SG70Q8d9") mdcHRSS[i]=1.219e-21;
99  if(TString(MDC->wfList[i].name)=="SG100Q8d9") mdcHRSS[i]=5.879e-22;
100  if(TString(MDC->wfList[i].name)=="SG153Q8d9") mdcHRSS[i]=3.170e-22;
101  if(TString(MDC->wfList[i].name)=="SG235Q8d9") mdcHRSS[i]=3.148e-22;
102  if(TString(MDC->wfList[i].name)=="SG361Q8d9") mdcHRSS[i]=5.268e-22;
103  if(TString(MDC->wfList[i].name)=="SG554Q8d9") mdcHRSS[i]=5.772e-22;
104  if(TString(MDC->wfList[i].name)=="SG849Q8d9") mdcHRSS[i]=9.610e-22;
105  if(TString(MDC->wfList[i].name)=="SG1053Q9") mdcHRSS[i]=1.025e-21;
106  if(TString(MDC->wfList[i].name)=="SG1304Q9") mdcHRSS[i]=1.230e-21;
107  if(TString(MDC->wfList[i].name)=="SG1615Q9") mdcHRSS[i]=1.681e-21;
108  if(TString(MDC->wfList[i].name)=="SG2000Q9") mdcHRSS[i]=2.930e-21;
109  if(TString(MDC->wfList[i].name)=="SG70Q100") mdcHRSS[i]=1.650e-21;
110  if(TString(MDC->wfList[i].name)=="SG235Q100") mdcHRSS[i]=3.557e-22;
111  if(TString(MDC->wfList[i].name)=="SG849Q100") mdcHRSS[i]=8.918e-22;
112  if(TString(MDC->wfList[i].name)=="SG1615Q100") mdcHRSS[i]=1.555e-21;
113  if(TString(MDC->wfList[i].name)=="WNB100_100_0d100") mdcHRSS[i]=4.378e-22;
114  if(TString(MDC->wfList[i].name)=="WNB250_100_0d100") mdcHRSS[i]=4.405e-22;
115  if(TString(MDC->wfList[i].name)=="WNB1000_10_0d100") mdcHRSS[i]=9.035e-22;
116  if(TString(MDC->wfList[i].name)=="WNB1000_1000_0d100") mdcHRSS[i]=3.373e-21;
117  if(TString(MDC->wfList[i].name)=="WNB1000_1000_0d010") mdcHRSS[i]=1.841e-21;
118  }
119 
120 
121  // --------------------------------------------------------
122  // define event sky distribution
123  // --------------------------------------------------------
124  par.resize(3);
125  par[0].name="entries";par[0].value=100000; // pool of events
126  par[1].name="rho_min";par[1].value=.2; // min rho Kpc
127  par[2].name="rho_max";par[2].value=40.; // max rho Kpc
128  MDC->SetSkyDistribution(MDC_RANDOM,"x+50./x",par,seed);
129 
130 }
waveform
Definition: mdc.hh:189
mdcHRSS
std::vector< double > mdcHRSS(K)
wf_name
char wf_name[256]
Definition: CWB_Plugin_Sim4_BRST_Config.C:16
CWB::mdc::SetInjRate
void SetInjRate(double inj_rate=MDC_INJ_RATE)
Definition: mdc.hh:279
MDC_GA
Definition: mdc.hh:128
n
int n
Definition: cwb_net.C:10
CWB::mdc::SetSourceListSeed
void SetSourceListSeed(unsigned int srcList_seed)
Definition: mdc.hh:288
TString
TString("c")
Definition: cwb_report_skymap.C:111
network::nRun
size_t nRun
Definition: network.hh:554
CWB_PLUGIN_IMPORT
CWB_PLUGIN_IMPORT(network **, net)
seed
int seed
Definition: CWB_Plugin_Sim4_BRST_Config.C:11
m
int m
Definition: cwb_net.C:10
i
i drho i
Definition: cwb_epparameters.C:85
CWB::mdc::AddWaveform
mdcid AddWaveform(MDC_TYPE mdc_type, vector< mdcpar > par, TString uname="")
Definition: mdc.cc:445
CWB::mdc::Print
void Print(int level=0)
Definition: mdc.cc:2707
CWB::mdc::SetSkyDistribution
void SetSkyDistribution(MDC_DISTRIBUTION sky_distribution, vector< mdcpar > par, int seed=0, bool add=false)
Definition: mdc.cc:3415
CWB::mdc
Definition: mdc.hh:216
F_Q9
double F_Q9[4]
Definition: CWB_Plugin_Sim4_BRST_Config.C:44
net
network ** net
Definition: CWB_Plugin_Sim4_BRST_Config.C:1
gIFACTOR
int gIFACTOR
Definition: cbc_plots_sim4.C:115
F_WNB
double F_WNB[5]
Definition: CWB_Plugin_Sim4_BRST_Config.C:59
CWB::mdc::SetInjJitter
void SetInjJitter(double inj_jitter=MDC_INJ_JITTER)
Definition: mdc.hh:285
e
double e
Definition: CreateListEBBH.C:35
CWB::mdc::wfList
vector< waveform > wfList
Definition: mdc.hh:355
F_Q8d9
double F_Q8d9[7]
Definition: CWB_Plugin_Sim4_BRST_Config.C:37
MDC_SG
Definition: mdc.hh:119
F_Q3
double F_Q3[4]
Definition: CWB_Plugin_Sim4_BRST_Config.C:30
D_WNB
double D_WNB[5]
Definition: CWB_Plugin_Sim4_BRST_Config.C:61
F_GA
double F_GA[4]
Definition: CWB_Plugin_Sim4_BRST_Config.C:22
B_WNB
double B_WNB[5]
Definition: CWB_Plugin_Sim4_BRST_Config.C:60
F_Q100
double F_Q100[4]
Definition: CWB_Plugin_Sim4_BRST_Config.C:51
K
int K
Definition: CWB_Plugin_Sim4_BRST_Config.C:87
CWB::mdc::SetInjHrss
void SetInjHrss(double inj_hrss=MDC_INJ_HRSS)
Definition: mdc.hh:275
MDC
CWB::mdc * MDC
Definition: CWB_Plugin_Sim4_BRST_Config.C:13
MDC_WNB
Definition: mdc.hh:127
par
vector< mdcpar > par
Definition: CWB_Plugin_Sim4_BRST_Config.C:18
wf
waveform wf
Definition: CWB_Plugin_Sim4_BRST_Config.C:17