Logo coherent WaveBurst  
Library Reference Guide
Logo
TestWriteFrameMDC.C
Go to the documentation of this file.
1 {
2  //
3  // Create waveforms and write to frame file
4  // Author : Gabriele Vedovato
5 
6  #define N_IFO 3
7  #define NWAVEFORM 15
8 
9  double F[NWAVEFORM];
10  double Q[NWAVEFORM];
11 
12  F[0]=70; Q[0]=8.9;
13  F[1]=100; Q[1]=8.9;
14  F[2]=153; Q[2]=8.9;
15  F[3]=235; Q[3]=8.9;
16  F[4]=361; Q[4]=8.9;
17  F[5]=554; Q[5]=8.9;
18  F[6]=849; Q[6]=8.9;
19  F[7]=945; Q[7]=9;
20  F[8]=1053; Q[8]=9;
21  F[9]=1172; Q[9]=9;
22  F[10]=1304; Q[10]=9;
23  F[11]=1451; Q[11]=9;
24  F[12]=1615; Q[12]=9;
25  F[13]=1797; Q[13]=9;
26  F[14]=2000; Q[14]=9;
27 
28 
29 
30  TString ifo[N_IFO]={"L1","H1","V1"};
31  CWB::mdc MDC(N_IFO,ifo);
32 
33  char wf_name[256];
35  for (int n=0;n<NWAVEFORM;n++) {
36  sprintf(wf_name,"SG%fQ%1.1f",F[n],Q[n]);
37  wf.name = wf_name;
38  wf.name.ReplaceAll(".","d");
39  wf.name.ReplaceAll("d0","");
40  wf.hp = MDC.GetSGQ(F[n],Q[n]);
41  wf.hx = wf.hp;
42  wf.hx = 0;
43  MDC.AddWaveform(wf);
44  }
45  MDC.Print();
46 
47 
48  vector<mdcpar> par(3);
49  par[0].name="entries"; par[0].value=100000;
50  par[1].name="rho_min"; par[1].value=5;
51  par[2].name="rho_max"; par[2].value=10;
53 
54  //par.resize(1);
55  //par[0].name="entries"; par[0].value=100000;
56  //MDC.SetSkyDistribution(MDC_MNGD,par);
57 
58  //vector<mdcpar> par(1);
59  //par[0].name="distance_thr"; par[0].value=100;
60  //MDC.SetSkyDistribution(MDC_GWGC,"../data/GWGCCatalog_Rev1d7.txt",par);
61 
62 /*
63  par.resize(3);
64  par[0].name="theta"; par[0].value=30;
65  par[1].name="phi"; par[0].value=60;
66  par[2].name="rho"; par[0].value=1;
67  MDC.SetSkyDistribution(MDC_EARTH_FIX,par);
68 */
69 
70  MDC.SetInjRate(0.1);
71  MDC.SetInjJitter(1.0);
72 
73  TString frDir = "frames";
74  TString frLabel = "TEST";
75  size_t gps = 123456789;
76  size_t length = 1000;
77  bool xlog = true;
78 
79  TString ofName = MDC.WriteFrameFile(frDir, frLabel, gps, length, xlog);
80 /*
81  for(int n=0;n<3;n++) {
82  MDC.WriteFrameFile(frDir, frLabel, gps, length, xlog);
83  gps+=length;
84  }
85 */
86  //cout << ofName.Data() << endl;
87 
88  exit(0);
89 }
#define NWAVEFORM
TString ofName
Definition: mdc.hh:219
size_t length
void SetInjRate(double inj_rate=MDC_INJ_RATE)
Definition: mdc.hh:311
int n
Definition: cwb_net.C:28
TString("c")
wavearray< double > GetSGQ(double frequency, double Q)
Definition: mdc.cc:3005
wavearray< double > hp
Definition: mdc.hh:226
double Q[NWAVEFORM]
mdcid AddWaveform(MDC_TYPE mdc_type, vector< mdcpar > par, TString uname="")
Definition: mdc.cc:472
void Print(int level=0)
Definition: mdc.cc:2736
#define N_IFO
void SetSkyDistribution(MDC_DISTRIBUTION sky_distribution, vector< mdcpar > par, int seed=0, bool add=false)
Definition: mdc.cc:3444
exit(0)
TString name
Definition: mdc.hh:222
Definition: mdc.hh:248
TString WriteFrameFile(TString frDir, TString frLabel, size_t gps, size_t length=1000, bool log=false, vector< TString > chName=vector< TString >())
Definition: mdc.cc:4283
waveform wf
F[0]
char wf_name[256]
void SetInjJitter(double inj_jitter=MDC_INJ_JITTER)
Definition: mdc.hh:317
TString frDir
wavearray< double > hx
Definition: mdc.hh:227
TString ifo[N_IFO]
size_t gps
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
TString frLabel
CWB::mdc MDC(N_IFO, ifo)
vector< mdcpar > par(3)
bool xlog