21 #pragma GCC system_header 28 #include "TObjArray.h" 29 #include "TObjString.h" 41 #define SGW_MDC_POLARIZATION "SCALAR" 42 #define SGW_SEARCH_POLARIZATION "SCALAR" 71 cout <<
"-----> macro/CWB_Plugin_SGW.C" << endl;
72 cout <<
"ifo " << ifo.Data() << endl;
73 cout <<
"type " << type << endl;
81 if(gOPT.search_polarization==
"SCALAR") cfg->
gamma = -1.0;
124 for(
int j=0;
j<token->GetEntries();
j++) {
126 TObjString* tok = (TObjString*)token->At(
j);
127 TString stok = tok->GetString();
129 if(stok.Contains(
"--sgw_mdc_polarization=")) {
130 gOPT.mdc_polarization=stok;
131 gOPT.mdc_polarization.Remove(0,gOPT.mdc_polarization.Last(
'=')+1);
132 gOPT.mdc_polarization.ToUpper();;
134 if(gOPT.mdc_polarization!=
"TENSOR" && gOPT.mdc_polarization!=
"SCALAR") {
135 cout <<
"CWB_Plugin_SGW Error : parameter sgw_mdc_polarization must be : SCALAR/TENSOR" << endl;
139 if(stok.Contains(
"--sgw_search_polarization=")) {
140 gOPT.search_polarization=stok;
141 gOPT.search_polarization.Remove(0,gOPT.search_polarization.Last(
'=')+1);
142 gOPT.search_polarization.ToUpper();;
144 if(gOPT.search_polarization!=
"TENSOR" && gOPT.search_polarization!=
"SCALAR") {
145 cout <<
"CWB_Plugin_SGW Error : parameter sgw_search_polarization must be : SCALAR/TENSOR" << endl;
162 cout <<
"-----------------------------------------" << endl;
163 cout <<
"SGW config options " << endl;
164 cout <<
"-----------------------------------------" << endl << endl;
165 cout <<
"SGW_MDC_POLARIZATION " << gOPT.mdc_polarization << endl;
166 cout <<
"SGW_SEARCH_POLARIZATION " << gOPT.search_polarization << endl;
detector * getifo(size_t n)
param: detector index
void setAntenna(detector *)
param: detector (use theta, phi index array)
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<< "\ layers : "<< nLAYERS<< "\ dF(hz) : "<< dF<< "\ 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
void PrintUserOptions(CWB::config *cfg)
void setPolarization(POLARIZATION polarization=TENSOR)
network ** net
NOISE_MDC_SIMULATION.
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
COHERENCE.
#define SGW_SEARCH_POLARIZATION
void setDelay(const char *="L1")
void ReadUserOptions(TString options)
void setSkyMaps(double, double=0., double=180., double=0., double=360.)
param: sky map granularity step, degrees param: theta begin, degrees param: theta end...
#define SGW_MDC_POLARIZATION