7 #define SKYMAP_FILE "/home/vedovato/waveburst/coherence/offline/SIMULATION/ADV_SIM_SGQ9_L1H1V1_run5/data/ced_931158378_44_ADV_SIM_SGQ9_L1H1V1_run5_30_job1/L1H1V1_931158395.102_931158395.102_931158395.102/probability.root" 9 #define SpeedOfLightInVacuo 2.99792458000000000e+08 11 #define ALPHA_FACTOR 4 18 gSystem->Load(
"libMathCore");
34 for(
int l=0;
l<
L;
l++) osm[
n]->
set(
l,0);
40 XYZVector
L1(pD[0]->Rv[0],pD[0]->Rv[1],pD[0]->Rv[2]);
41 XYZVector
H1(pD[1]->Rv[0],pD[1]->Rv[1],pD[1]->Rv[2]);
42 XYZVector
V1(pD[2]->Rv[0],pD[2]->Rv[1],pD[2]->Rv[2]);
49 XYZVector D3 = D1.Cross(
V1-L1);
50 XYZVector D2 = D1.Cross(D3);
52 XYZVector eD1 = D1.Unit();
53 XYZVector eD2 = D2.Unit();
54 XYZVector eD3 = D3.Unit();
56 Rotation3D R(eD1,eD2,eD3);
58 XYZVector H1L1 =
H1-
L1;
60 XYZVector V1L1 =
V1-
L1;
62 XYZVector V1H1 =
V1-
H1;
65 TVector3 vH1L1(H1L1.X(),H1L1.Y(),H1L1.Z());
66 TVector3 vV1L1(V1L1.X(),V1L1.Y(),V1L1.Z());
67 TVector3 vV1H1(V1H1.X(),V1H1.Y(),V1H1.Z());
68 TVector3 veD3(eD3.X(),eD3.Y(),eD3.Z());
81 double alpha =
k*dalpha;
83 oR.Rotate(alpha,veD3.Unit());
85 if(
n==0) roSD = oR*vH1L1;
86 if(
n==1) roSD = oR*vV1L1;
87 if(
n==2) roSD = oR*vV1H1;
90 for (
int i=0;
i<
N;
i++) {
94 if(
n==0) vR.Rotate(angle,vH1L1.Unit());
95 if(
n==1) vR.Rotate(angle,vV1L1.Unit());
96 if(
n==2) vR.Rotate(angle,vV1H1.Unit());
98 TVector3 rvSD = vR*roSD;
99 double theta = rvSD.Theta()*r2d;
100 double phi = rvSD.Phi()*r2d;
101 if (phi<0) phi+=360;
if (phi>360) phi-=360;
102 if (theta<0) theta+=180;
if (theta>180) theta-=180;
108 for (
int i=0;
i<
N;
i++) {
115 for (
int i=0;i<
N;i++) {
123 for(
int l=0;l<
L;l++) {
125 for (
int n=0;
n<
nIFO;
n++) like+=osm[
n]->
get(l);
detector * getifo(size_t n)
param: detector index
void set(size_t i, double a)
void Draw(int dpaletteId=1, Option_t *option="colfz")
size_t getSkyIndex(double th, double ph)
param: theta param: phi
detector V1((char *)"V1")
#define SpeedOfLightInVacuo
void set(size_t i, double a)
param: sky index param: value to set
detector H1((char *)"H1")
double get(size_t i)
param: sky index
void SetOptions(TString projection="hammer", TString coordinate="Geographic", double resolution=1, bool goff=false)