Logo coherent WaveBurst  
Library Reference Guide
Logo
gwosc_parameters_lhv_gwtc-1.C
Go to the documentation of this file.
1 {
2  strcpy(analysis,"2G");
3 
4  nIFO = 3;
5  cfg_search = 'r';
6  optim=false;
7 
8  strcpy(ifo[0],"L1");
9  strcpy(ifo[1],"H1");
10  strcpy(ifo[2],"V1");
11  strcpy(refIFO,"L1");
12 
13  //lags
14  lagSize = 1;
15  lagStep = 1.;
16  lagOff = 0;
17  lagMax = 0;
18 
19  //superlags
20  slagSize = 0; // number of super lags (simulation=1) - if slagSize=0 -> Standard Segments
21  slagMin = 0;
22  slagMax = 0;
23  slagOff = 0;
24 
25  //jobs
26  segLen = 1200;
27  segMLS = 600;
28  segTHR = 200;
29  segEdge = 10;
30 
31  //frequency
32  fLow = 16.; // low frequency of the search
33  fHigh = 1024.; // high frequency of the search
34 
35  levelR = 3;
36  l_low = 4; // low frequency resolution level
37  l_high = 10; // high frequency resolution level
38 
39  strcpy(wdmXTalk,"wdmXTalk/OverlapCatalog16-1024.bin");
40 
42 
43  bpp = 0.001;
44  subnet = 0.5;
45  subcut = 0.0;
46  netRHO = 5.5;
47  netCC = 0.5;
48  Acore = 1.7;
49  Tgap = 0.2;
50  Fgap = 128.0;
51  delta = 0.5;
52  cfg_gamma = -1.0;
53  LOUD = 300;
54 
55  pattern = 5;
56 
57  iwindow = 30.;
58 
59  //precision=GetPrecision(100,5);
60 
61 // nSky = -99; // save pixels skymap probability pixels up to cumulative prob < 0.99
62  nSky=196608; // save all pixels
63 
64 // jobfOptions = CWB_JOBF_SAVE_CSTRAIN; // save cstrain data into the final ROOT file
65 
66  //simulation
67  nfactor = 1;
69 
70  strcpy(channelNamesRaw[0],"L1:GWOSC-4KHZ_R1_STRAIN");
71  strcpy(channelNamesRaw[1],"H1:GWOSC-4KHZ_R1_STRAIN");
72  strcpy(channelNamesRaw[2],"V1:GWOSC-4KHZ_R1_STRAIN");
73 
74  strcpy(frFiles[0],"input/L1_frames.in");
75  strcpy(frFiles[1],"input/H1_frames.in");
76  strcpy(frFiles[2],"input/V1_frames.in");
77 
78  // dq file list
79  // {ifo, dqcat_file, dqcat[0/1/2], shift[sec], inverse[false/true], 4columns[true/false]}
80  nDQF=12;
81  dqfile dqf[12]={
82  {"L1" ,"input/cwb_period.txt", CWB_CAT0, 0., false, false},
83  {"H1" ,"input/cwb_period.txt", CWB_CAT0, 0., false, false},
84  {"V1" ,"input/cwb_period.txt", CWB_CAT0, 0., false, false},
85  {"L1" ,"../DQ/L1_cat0.txt", CWB_CAT0, 0., false, false},
86  {"H1" ,"../DQ/H1_cat0.txt", CWB_CAT0, 0., false, false},
87  {"V1" ,"../DQ/V1_cat0.txt", CWB_CAT0, 0., false, false},
88  {"L1" ,"../DQ/L1_cat1.txt", CWB_CAT1, 0., false, false},
89  {"H1" ,"../DQ/H1_cat1.txt", CWB_CAT1, 0., false, false},
90  {"V1" ,"../DQ/V1_cat1.txt", CWB_CAT1, 0., false, false},
91  {"L1" ,"../DQ/L1_cat2.txt", CWB_CAT2, 0., false, false},
92  {"H1" ,"../DQ/H1_cat2.txt", CWB_CAT2, 0., false, false},
93  {"V1" ,"../DQ/V1_cat2.txt", CWB_CAT2, 0., false, false}
94  };
95  for(int i=0;i<12;i++) DQF[i]=dqf[i];
96 
97 
98  plugin = TMacro("macro/CWB_Plugin_WF.C"); // Macro source
99  configPlugin = TMacro("macro/CWB_Plugin_Config.C"); // Macro config
100  //plugin.SetTitle("macro/CWB_Plugin_WF_C.so");
101 
102  TString optwf = ""; // NOTE : add space at the end of each line
103 
104  //optwf += "wf_output_disable=root "; // disable output root file (to be used when QLveto is enabled)
105  optwf += "wf_output_enable=root "; // enable output root file (to be used when QLveto is enabled)
106  optwf += "wf_output_disable=inj "; // disable save injection to the output root file
107  optwf += "wf_output_enable=rec "; // enable save reconstructed waveform to the output root file
108  //optwf += "wf_output_disable=wht "; // disable save whitened data to the output root file
109  //optwf += "wf_output_disable=dat "; // disable save rec+null data to the output root file
110  //optwf += "wf_output_disable=nul "; // disable save null data to the output root file
111  optwf += "wf_output_enable=wht "; // enable save whitened data to the output root file
112  optwf += "wf_output_enable=dat "; // enable save rec+null data to the output root file
113  optwf += "wf_output_enable=nul "; // enable save null data to the output root file
114 
115  optwf += "wf_inj_tstep=150.000 "; // is the injection step time (used only for PE simulations, must be >0)
116 
117  optwf += "sn_strain_file_name=input/L1_ASD.txt "; // add L1 strain file name to first detector
118  optwf += "sn_strain_file_name=input/H1_ASD.txt "; // add H1 strain file name to second detector
119  optwf += "sn_strain_file_name=input/V1_ASD.txt "; // add V1 strain file name to second detector
120  optwf += "sn_seed=1 "; // add seed for random noise generation to L1 detector
121  optwf += "sn_seed=2 "; // add seed for random noise generation to H1 detector
122  optwf += "sn_seed=3 "; // add seed for random noise generation to V1 detector
123 
124  strcpy(parPlugin,optwf.Data()); // set WF plugin parameters
125  strcpy(comment,"GWOSC");
126 }
char wdmXTalk[1024]
strcpy(ifo[0],"L1")
char parPlugin[1024]
char channelNamesRaw[NIFO_MAX][50]
TString("c")
char comment[1024]
char frFiles[NIFO_MAX+1][256]
Definition: test_config1.C:166
char refIFO[4]
Definition: test_config1.C:14
i drho i
char ifo[NIFO_MAX][8]
dqfile dqf[12]
dqfile DQF[12]
Definition: test_config1.C:171
TMacro plugin
TMacro configPlugin