wmdc_config_sgq9_lf.C

Go to the documentation of this file.

{
  // -------------------------------------------------------- 
  // define job parameters
  // -------------------------------------------------------- 
  TString frDir       = "frames";
  TString frLabel     = "TestWaveMDC";

  // define frame segment list
  TString segmentList = "Segments/S6a_LS-segs.txt";

  // select segments
  int jobmin  = 210;    // start segment
  int jobmax  = 312;    // end segment
  int jobstep =  10;    // frames per job 

  // -------------------------------------------------------- 
  // define network
  // -------------------------------------------------------- 
  int nIFO=3;
  TString ifo[nIFO]={"L1","H1","V1"};
  CWB::mdc MDC(nIFO,ifo); 

  // -------------------------------------------------------- 
  // define waveforms
  // -------------------------------------------------------- 
  int nWF=15;
  double F[nWF] = {70,100,153,235,361,554,849,945,1053,1172,1304,1451,1615,1797,2000};
  double Q[nWF] = {8.9,8.9,8.9,8.9,8.9,8.9,8.9,9,9,9,9,9,9,9,9};

  char wf_name[256];
  waveform wf;
  wf.par.resize(2);
  for (int n=0;n<nWF;n++) {
    sprintf(wf_name,"SG%dQ%1.1f",F[n],Q[n]);
    wf.name = wf_name;
    wf.name.ReplaceAll(".","d");
    wf.name.ReplaceAll("d0","");
    wf.hp = MDC.GetSGQ(F[n],Q[n]);
    wf.hpPath = wf.name;
    wf.hx = wf.hp;
    wf.hx = 0;
    wf.hxPath = "";
    wf.par[0].name="freq"; wf.par[0].value=F[n];
    wf.par[1].name="Q";    wf.par[1].value=Q[n];
    MDC.AddWaveform(wf);
  }
  MDC.Print();

  // -------------------------------------------------------- 
  // define injection parameters
  // -------------------------------------------------------- 
  MDC.SetInjHrss(2.5e-21);
  MDC.SetInjRate(0.01);
  MDC.SetInjJitter(10.0);

  // --------------------------------------------------------
  // define sky distribution
  // --------------------------------------------------------
  vector<mdcpar> par; par.resize(3);
  par[0].name="entries";par[0].value=100000;  // pool of events
  par[1].name="rho_min";par[1].value=1;         // min rho // Mpc
  par[2].name="rho_max";par[2].value=100;       // max rho // Mpc
  MDC.SetSkyDistribution(MDC_RANDOM,par);

  //vector<mdcpar> par; par.resize(1);
  //par[0].name="entries";par[0].value=100000;  // number of entries
  //MDC.SetSkyDistribution(MDC_MNGD,par);

  //vector<mdcpar> par; par.resize(1);
  //par[0].name="distance_thr";par[0].value=100;  // distance max = 100 Mpc
  //MDC.SetSkyDistribution(MDC_GWGC,"../data/GWGCCatalog_Rev1d2.txt",par);

  //vector<mdcpar> par; par.resize(3);
  //par[0].name="theta";par[0].value=30;     // theta
  //par[1].name="phi";  par[1].value=60;     // phi
  //par[2].name="psi";  par[2].value=90;     // psi
  //MDC.SetSkyDistribution(MDC_EARTH_FIX,par);

}