Logo Coherent WaveBurst  
Reference Guide 		Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
CWB_Plugin_SimNoise.C
Go to the documentation of this file.
1 #define XIFO 4
2 
3 #pragma GCC system_header
4 
5 #include "cwb.hh"
6 #include "config.hh"
7 #include "network.hh"
8 #include "wavearray.hh"
9 #include "TString.h"
10 #include "TObjArray.h"
11 #include "TObjString.h"
12 #include "TRandom.h"
13 #include "Toolbox.hh"
14 
15 void
16 CWB_Plugin(TFile* jfile, CWB::config* cfg, network* net, WSeries<double>* x, TString ifo, int type) {
17 
18 // Plugin to generate simulated gaussian noise
19 
20  cout << endl;
21  cout << "-----> CWB_Plugin_GenNoise.C" << endl;
22  cout << "ifo " << ifo.Data() << endl;
23  cout << "type " << type << endl;
24  cout << endl;
25 
26 
27  if(type==CWB_PLUGIN_CONFIG) {
28  cfg->dataPlugin=true; // disable read data from frames
29  }
30 
31  if(type==CWB_PLUGIN_DATA) {
32 
33  CWB::Toolbox TB;
34 
35  int seed;
36  if(ifo.CompareTo("L1")==0) seed=1000;
37  if(ifo.CompareTo("H1")==0) seed=2000;
38  if(ifo.CompareTo("V1")==0) seed=3000;
39  if(ifo.CompareTo("J1")==0) seed=4000;
40  if(ifo.CompareTo("A2")==0) seed=5000;
41 
42  TString fName;
43  if(ifo.CompareTo("L1")==0) fName="plugins/strains/advLIGO_NSNS_Opt_8khz_one_side.txt";
44  if(ifo.CompareTo("H1")==0) fName="plugins/strains/advLIGO_NSNS_Opt_8khz_one_side.txt";
45  if(ifo.CompareTo("V1")==0) fName="plugins/strains/advVIRGO_sensitivity_12May09_8khz_one_side.txt";
46  if(ifo.CompareTo("J1")==0) fName="plugins/strains/LCGT_sensitivity_8khz_one_side.txt";
47  if(ifo.CompareTo("A2")==0) fName="plugins/strains/advLIGO_NSNS_Opt_8khz_one_side.txt";
48 
49  int size=x->size();
50  double start=x->start();
51  TB.getSimNoise(*x, fName, seed, net->nRun);
52  x->resize(size);
53  x->start(start);
54  }
55  return;
56 }
cfg
CWB::config * cfg
Definition: TestCWB_Plugin.C:5
WSeries::resize
virtual void resize(unsigned int)
Definition: wseries.cc:883
wavearray::size
virtual size_t size() const
Definition: wavearray.hh:127
CWB::config::dataPlugin
bool dataPlugin
Definition: config.hh:346
TString
TString("c")
Definition: cwb_report_skymap.C:111
network::nRun
size_t nRun
Definition: network.hh:554
x
cout<< endl;cout<< "ts size = "<< ts.size()<< " ts rate = "<< ts.rate()<< endl;tf.Forward(ts, wdm);int levels=tf.getLevel();cout<< "tf size = "<< tf.size()<< endl;double dF=tf.resolution();double dT=1./(2 *dF);cout<< "rate(hz) : "<< RATE<< "\t layers : "<< nLAYERS<< "\t dF(hz) : "<< dF<< "\t dT(ms) : "<< dT *1000.<< endl;int itime=TIME_PIXEL_INDEX;int ifreq=FREQ_PIXEL_INDEX;int index=(levels+1)*itime+ifreq;double time=itime *dT;double freq=(ifreq >0)?ifreq *dF:dF/4;cout<< endl;cout<< "PIXEL TIME = "<< time<< " sec "<< endl;cout<< "PIXEL FREQ = "<< freq<< " Hz "<< endl;cout<< endl;wavearray< double > x
Definition: BaseFunctionWDM.C:51
size
Long_t size
Definition: cwb_condor_check.C:48
wavearray::start
virtual void start(double s)
Definition: wavearray.hh:119
TB
CWB::Toolbox TB
Definition: ComputeSNR.C:5
ifo
char ifo[NIFO_MAX][8]
Definition: cwb1G_parameters.C:17
net
network ** net
NOISE_MDC_SIMULATION.
Definition: CWB_Plugin_HEN_BKG_Config.C:4
jfile
jfile
Definition: cwb_job_obj.C:25
CWB::Toolbox::getSimNoise
static void getSimNoise(wavearray< double > &u, TString fName, int seed, int run)
Definition: Toolbox.cc:4809
seed
int seed
Definition: CWB_Plugin_HEN_BKG_Config.C:13
CWB_Plugin
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
COHERENCE.
Definition: CWB_Plugin_SimNoise.C:16
fName
char fName[256]
Definition: CreateListEBBH.C:25