library/html/CWB__Plugin__Skymap_8C_source.html

 #define XIFO 4


 #pragma GCC system_header


 #include "cwb.hh"

 #include "config.hh"

 #include "network.hh"

 #include "wavearray.hh"

 #include "TString.h"

 #include "TObjArray.h"

 #include "TObjString.h"

 #include "TRandom.h"

 #include "TComplex.h"

 #include "TMath.h"

 #include "mdc.hh"

 #include <vector>


 void

 CWB_Plugin(TFile* jfile, CWB::config* cfg, network* net, WSeries<double>* x, TString ifo, int type)  {

 //!EVENT_REPORT

 // This plugin save skymap to output root file


   cout << endl;

   cout << "-----> CWB_Plugin_Skymap.C" << endl;

   cout << "ifo " << ifo.Data() << endl;

   cout << "type " << type << endl;

   cout << endl;


   if(type==CWB_PLUGIN_OLIKELIHOOD) {


     int i,j,k,n,m,l;

     int nIFO = net->ifoListSize();

     int K = net->nLag;

     int M = net->mdc__ID.size();

     int ID;

     char search = net->tYPe;


     wavearray<double> id;


     bool batch = gROOT->IsBatch();

     gROOT->SetBatch(true);


     watplot SMS(const_cast<char*>("sms"),200,20,800,400);


     char ifostr[20]  = "";

     char strtime[1024]  = "";

     char fname[1024];


     int    minTimeDet=nIFO;

     double minTime=1e40;

     double eventTime[NIFO];

     double lagTime[NIFO];

     int ifoid[NIFO],sortid[NIFO];

     double factor = 1;


     //Fill in all skymaps

     double old_cc = net->netCC;

     double old_rho = net->netRHO;

     net->netCC = -1;

     net->netRHO = 0;


     // search output root file in the system list

     TFile* froot = NULL;

     TList *files = (TList*)gROOT->GetListOfFiles();

     if (files) {

       TIter next(files);

       TSystemFile *file;

       TString fname;

       bool check=false;

       while ((file=(TSystemFile*)next())) {

          fname = file->GetName();

          // set output root file as the current file

          if(fname.Contains("wave_")) {froot=(TFile*)file;froot->cd();}

       }

       if(!froot) {

         cout << "CWB_Plugin_skymap.C : Error - output root file not found" << endl;

         exit(1);

       }

     } else {

       cout << "CWB_Plugin_skymap.C : Error - output root file not found" << endl;

       exit(1);

     }


     TDirectory* cdskymap = (TDirectory*)froot->Get("skymap");

     if(cdskymap==NULL) cdskymap = froot->mkdir("skymap");


     netevent* EVT = new netevent(nIFO);


     int rate = int(2*net->getifo(0)->TFmap.resolution(0)+0.5);


     for(k=0; k<K; k++) {                                // loop over the lags


       id = net->getwc(k)->get(const_cast<char*>("ID"), 0, 'L', rate);


       for(j=0; j<(int)id.size(); j++) {                 // loop over cluster index


         ID = size_t(id.data[j]+0.5);                    // cluster id


         EVT->output(NULL,net,factor,ID,k);              // get reconstructed parameters

         if(EVT->rho[1] < cfg->cedRHO) continue;         // skip events under rho threshold


         int masterDet=0;

         int lagMin=2147483647;

         for(n=0; n<nIFO; n++) if(EVT->lag[n]<lagMin) {lagMin=int(EVT->lag[n]);masterDet=n;}


         net->likelihood(search, net->acor, ID, k);      // exec likelihood search


         // get event network time

         double gps_start = EVT->time[masterDet]-EVT->duration[masterDet];

         double gps_stop  = EVT->time[masterDet]+EVT->duration[masterDet];


         // create a unique label

         for(n=0; n<nIFO; n++) sprintf(strtime, "%s_%.3f",strtime,EVT->start[n]);


         //likelihood skymaps

         TMarker inj(EVT->phi[1],90.-EVT->theta[1], 29);  // injected pos     (white star)

         TMarker rec(EVT->phi[0],90.-EVT->theta[0], 29);  // recostructed pos (black star)

         TMarker det(EVT->phi[3],90.-EVT->theta[3], 20);  // detected pos     (black dot )

         inj.SetMarkerSize(2.0); inj.SetMarkerColor(kWhite);

         rec.SetMarkerSize(2.0); rec.SetMarkerColor(kBlack);

         det.SetMarkerSize(1.0); det.SetMarkerColor(kBlack);


         char evtdir[64];

         sprintf(evtdir, "run_%d_id_%d",net->nRun,ID);

         TDirectory* cdevtdir = (TDirectory*)cdskymap->Get(evtdir);

         if(cdevtdir==NULL) cdevtdir = cdskymap->mkdir(evtdir);

         cdevtdir->cd();


         SMS.canvas->cd();

         sprintf(fname, "sensitivity_plus");

         SMS.plot(net->nSensitivity, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "sensitivity_cross");

         SMS.plot(net->nAlignment, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "skystat");

         SMS.plot(net->nSkyStat, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "likelihood");

         SMS.plot(net->nLikelihood, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "null_energy");

         SMS.plot(net->nNullEnergy, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "corr_energy");

         SMS.plot(net->nCorrEnergy, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "penalty");

         SMS.plot(net->nPenalty, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "disbalance");

         SMS.plot(net->nDisbalance, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "correlation");

         SMS.plot(net->nCorrelation, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "netindex");

         SMS.plot(net->nNetIndex, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "ellipticity");

         SMS.plot(net->nEllipticity, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "polarisation");

         SMS.plot(net->nPolarisation, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();

         sprintf(fname, "probability");

         SMS.plot(net->nProbability, 0);

         rec.Draw(); det.Draw(); if(M) inj.Draw();

         SMS.canvas->Write(fname);

         SMS.clear();


       }                                                 // End loop on found events

     }                                                   // End loop on lags


     // restore NET thresholds

     net->netCC = old_cc;

     net->netRHO = old_rho;

   }


   return;

 }


network::getifo
detector * getifo(size_t n)
param: detector index
Definition: network.hh:418

cfg
CWB::config * cfg
Definition: TestCWB_Plugin.C:5

network::netCC
double netCC
Definition: network.hh:578

NIFO
#define NIFO
Definition: wat.hh:56

network::nLag
size_t nLag
Definition: network.hh:555

netevent::rho
Float_t * rho
biased null statistics
Definition: netevent.hh:93

netevent::start
Double_t * start
cluster duration = stopW-startW
Definition: netevent.hh:80

netevent::duration
Float_t * duration
max cluster time relative to segment start
Definition: netevent.hh:79

CWB::config::cedRHO
double cedRHO
Definition: config.hh:280

CWB_Plugin
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
COHERENCE.
Definition: CWB_Plugin_Skymap.C:20

wavearray< double >

n
int n
Definition: cwb_net.C:10

network::nSkyStat
skymap nSkyStat
Definition: network.hh:610

const_cast< char * >
cout<< "skymap size : "<< L<< endl;for(int l=0;l< L;l++) sm.set(l, l);sm > const_cast< char * >("skymap.dat")

TString
TString("c")
Definition: cwb_report_skymap.C:111

cwb.hh

network::nRun
size_t nRun
Definition: network.hh:554

ID
int ID
Definition: TestMDC.C:70

network::nPenalty
skymap nPenalty
Definition: network.hh:606

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

netcluster::get
wavearray< double > get(char *name, size_t index=0, char atype='R', int type=1, bool=true)
param: string with parameter name param: index in the amplitude array, which define detector param: c...
Definition: netcluster.cc:2188

ifostr
char ifostr[64]
Definition: cwb_report_skymap.C:124

SMS
#define SMS
Definition: FrDisplay.cc:108

network::nPolarisation
skymap nPolarisation
Definition: network.hh:612

size
Long_t size
Definition: cwb_condor_check.C:48

config.hh

M
#define M
Definition: UniqSLagsList.C:3

m
int m
Definition: cwb_net.C:10

network::mdc__ID
std::vector< size_t > mdc__ID
Definition: network.hh:597

j
int j
Definition: cwb_net.C:10

i
i drho i
Definition: cwb_epparameters.C:85

search
search
Definition: user_parameters.C:3

watplot::plot
void plot(wavearray< double > &, char *=NULL, int=1, double=0., double=0., bool=false, float=0., float=0., bool=false, float=0., bool=false)
Definition: watplot.cc:132

network::nCorrEnergy
skymap nCorrEnergy
Definition: network.hh:607

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

network::ifoListSize
size_t ifoListSize()
Definition: network.hh:413

netevent
Definition: netevent.hh:39

CWB_PLUGIN_OLIKELIHOOD
Definition: cwb.hh:88

network::nLikelihood
skymap nLikelihood
Definition: network.hh:604

watplot::clear
void clear()
Definition: watplot.hh:40

nIFO
#define nIFO
Definition: DrawSensitivitiesS5.C:11

watplot::canvas
TCanvas * canvas
Definition: watplot.hh:174

watplot
Definition: watplot.hh:27

factor
double factor
Definition: ReadCStrainFromJobFile.C:18

CWB::config
Definition: config.hh:73

jfile
jfile
Definition: cwb_job_obj.C:25

network
Definition: network.hh:41

int
i() int(T_cor *100))

next
TIter next(twave->GetListOfBranches())

netevent::lag
Float_t * lag
time between consecutive events
Definition: netevent.hh:65

fname
char fname[1024]
Definition: cwb_report_prod_1.C:110

WSeries< double >

k
int k
Definition: cwb_report_prod_2.C:150

network::acor
double acor
Definition: network.hh:567

network.hh

netevent::output
void output(TTree *=NULL, network *=NULL, double=0., size_t=0, int=-1)
Definition: netevent.cc:1166

network::likelihood
long likelihood(char='E', double=sqrt(2.), int=0, size_t=0, int=-1, bool=false)
Definition: network.cc:4415

mdc.hh

froot
TFile * froot
Definition: TestReadWriteDetectorObject.C:11

detector::TFmap
WSeries< double > TFmap
Definition: detector.hh:336

netevent::phi
Float_t * phi
sqrt(h+*h+ + hx*hx)
Definition: netevent.hh:68

network::nAlignment
skymap nAlignment
Definition: network.hh:602

EVT
netevent EVT(itree, nifo)

netevent::theta
Float_t * theta
[0]-reconstructed, [1]-injected phi angle, [2]-RA
Definition: netevent.hh:69

netevent::time
Double_t * time
beam pattern coefficients for hx
Definition: netevent.hh:75

network::nDisbalance
skymap nDisbalance
Definition: network.hh:609

cwb_online.file
string file
Definition: cwb_online.py:385

WSeries::resolution
double resolution(int=0)
Definition: wseries.hh:137

network::getwc
netcluster * getwc(size_t n)
param: delay index
Definition: network.hh:421

l
int l
Definition: cbc_plots.C:434

network::nNetIndex
skymap nNetIndex
Definition: network.hh:608

sprintf
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)

network::nSensitivity
skymap nSensitivity
list of wdm tranformations
Definition: network.hh:601

wavearray.hh

rate
pix rate
Definition: DrawClusterWithSkyplot.C:35

K
int K
Definition: CWB_Plugin_Sim4_BRST_Config.C:87

id
Long_t id
Definition: cwb_condor_check.C:48

network::nNullEnergy
skymap nNullEnergy
Definition: network.hh:605

network::nCorrelation
skymap nCorrelation
Definition: network.hh:603

cwb_online.files
list files
Definition: cwb_online.py:529

check
int check
Definition: CWB_configPlugin_Test.C:5

network::tYPe
char tYPe
Definition: network.hh:570

network::netRHO
double netRHO
Definition: network.hh:579

network::nProbability
skymap nProbability
Definition: network.hh:613

exit
exit(0)

network::nEllipticity
skymap nEllipticity
Definition: network.hh:611