Logo coherent WaveBurst  
Library Reference Guide
Logo
DrawSingleDetectorAntennaPatterns.C
Go to the documentation of this file.
1 //
2 // Draw Antenna Pattern Componets for a single detector
3 // Author : Gabriele Vedovato
4 
5 {
6  #define PROJECTION ""
7  //#define PROJECTION "hammer"
8  #define RESOLUTION 2
9  #define COORDINATES "Geographic"
10 
11  #define DISPLAY_WORLD_MAP
12  #define WORLD_MAP_DIR "$CWB_GWAT/data/"
13 
14  //#define IFO "L1"
15  //#define IFO "H1"
16  #define IFO "V1"
17 
18  #define FPLUS
19 
20  //#define OFILE_NAME "DetectorAntennaPattern.png"
21 
22  gskymap* gSM = new gskymap(int(6));
24 
25  detector D((char*)IFO);
26 
27  int L = gSM->size();
28  for(int l=0;l<L;l++) {
29  double th = gSM->getTheta(l);
30  double ph = gSM->getPhi(l);
31  wavecomplex F = D.antenna(th,ph,0);
32 #ifdef FPLUS
33  gSM->set(l,F.real());
34 #else
35  gSM->set(l,F.imag());
36 #endif
37  }
38 
39 #ifdef DISPLAY_WORLD_MAP
40  TString world_map = gSystem->ExpandPathName(WORLD_MAP_DIR);
41  gSM->SetWorldMapPath(world_map.Data());
42  gSM->SetWorldMap();
43 #endif
44 
45 #ifdef FPLUS
46  TString title = TString("F_{+}");
47 #else
48  TString title = TString("F_{x}");
49 #endif
50 
51  gSM->SetTitle(title);
52  gSM->SetGridxColor(kWhite);
53  gSM->SetGridyColor(kWhite);
54  gSM->Draw(0);
55 
56 #ifdef OFILE_NAME
57  cout << "Write : " << OFILE_NAME << endl;
58  gSM->Print(OFILE_NAME);
59  exit(0);
60 #endif
61 
62 }
gskymap * gSM
#define PROJECTION
double imag() const
Definition: wavecomplex.hh:70
TString("c")
void set(size_t i, double a)
Definition: gskymap.hh:128
TH2F * ph
double getTheta(size_t i)
Definition: skymap.hh:224
void Draw(int dpaletteId=1, Option_t *option="colfz")
Definition: gskymap.cc:460
void SetGridxColor(Color_t colorGridx=kBlack)
Definition: gskymap.hh:142
#define WORLD_MAP_DIR
#define COORDINATES
wavecomplex antenna(double, double, double=0.)
param: source theta,phi, polarization angle psi in degrees
Definition: detector.cc:490
#define OFILE_NAME
Definition: CreateSkyMask.C:8
double real() const
Definition: wavecomplex.hh:69
detector D((char *) IFO)
int l
void SetWorldMapPath(TString worldMapPath)
Definition: gskymap.hh:156
double F
void SetGridyColor(Color_t colorGridy=kBlack)
Definition: gskymap.hh:146
void SetTitle(TString title)
Definition: gskymap.hh:152
#define RESOLUTION
double getPhi(size_t i)
Definition: skymap.hh:164
void SetWorldMap(bool drawWorldMap=true)
Definition: gskymap.hh:154
size_t size()
Definition: skymap.hh:136
void Print(TString pname)
Definition: gskymap.cc:1122
#define IFO
void SetOptions(TString projection="hammer", TString coordinate="Geographic", double resolution=1, bool goff=false)
Definition: gskymap.cc:84
exit(0)