Logo Coherent WaveBurst  
Reference Guide 		Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Macros | Functions | Variables
CreateMNGDDistribution.C File Reference

Macros

#define A1   5.81
 
#define A2   17.43
 
#define A3   34.86
 
#define B   0.3
 
#define DENOMINATOR   "(pow(((x-[3])*(x-[3])+y*y)+pow([0]+sqrt(z*z+[1]*[1]),2),5./2.)*pow(z*z+[1]*[1],3./2.))"
 
#define DRAW_H2
 
#define Md1   6.6e10
 
#define Md2   -2.9e10
 
#define Md3   3.3e9
 
#define MNGD_FILENAME   "MNGDDistribution.txt"
 
#define NEVENTS   1000000
 
#define NUMERATOR   "([0]*((x-[3])*(x-[3])+y*y)+([0]+3*sqrt(z*z+[1]*[1]))*pow([0]+sqrt(z*z+[1]*[1]),2))"
 
#define SOLAR_SISTEM_DISTANCE_FROM_GC   7.62
 
#define XMAX   40
 
#define YMAX   4
 

Functions

canvas Clear ()
 
h2gd Draw ("colfz")
 
h2gd Fill (x, y)
 
 GalacticToEquatorial (ilongitude, ilatitude, olongitude, olatitude)
 
cout<< "gc_phi : "<< gc_phi
<< " gc_theta : "<< gc_theta
<< " "<< gc_rho<< endl;for(int
n=0;n< NEVENTS;n++){if(n%10000==0)
cout<< n<< endl;double x, y, z;gd-> 		GetRandom3 (x, y, z)
 
h2gd GetXaxis () -> SetNdivisions(-804)
 
h2gd GetYaxis () -> SetNdivisions(-804)
 
h2gd GetZaxis () -> SetLabelFont(42)
 
gStyle SetCanvasColor (kWhite)
 
canvas SetFillColor (kWhite)
 
canvas SetGridx ()
 
canvas SetGridy ()
 
gStyle SetLineColor (kWhite)
 
gStyle SetMarkerSize (2)
 
gStyle SetMarkerStyle (7)
 
gd SetNpx (100)
 
gd SetNpy (100)
 
gd SetNpz (100)
 
gStyle SetNumberContours (256)
 
gd1 SetParameter (0, A1)
 
gd1 SetParameter (1, B)
 
gd1 SetParameter (2, Md1)
 
gd1 SetParameter (3, SOLAR_SISTEM_DISTANCE_FROM_GC)
 
gd2 SetParameter (0, A2)
 
gd2 SetParameter (2, Md2)
 
gd3 SetParameter (0, A3)
 
gd3 SetParameter (2, Md3)
 
gStyle SetStatBorderSize (1)
 
h2gd SetStats (kFALSE)
 
gStyle SetTextFont (12)
 
h2gd SetTitle ("h2gd")
 
gStyle SetTitleColor (kBlue,"D")
 
gStyle SetTitleFillColor (kWhite)
 
gStyle SetTitleFont (12,"D")
 
gStyle SetTitleH (0.050)
 
gStyle SetTitleW (0.95)
 
gStyle SetTitleY (0.98)
 
xyz SetXYZ (xgc, ygc, zgc)
 
xyz SetXYZ (x, y, z)
 
 sprintf (formula,"[1]*[1]*[2]*%s/%s", NUMERATOR, DENOMINATOR)
 
canvas ToggleEventStatus ()
 

Variables

TCanvas * canvas = new TCanvas("AUNA", "GD", 300,40, 800, 800)
 
char formula [256]
 
double gc_phi = olongitude
 
double gc_rho = sqrt(xyz.Mag2())
 
double gc_theta = olatitude
 
TF3 * gd = new TF3("gd","gd1+gd2+gd3",-XMAX,XMAX,-XMAX,XMAX,-YMAX,YMAX)
 
TF3 * gd1 = new TF3("gd1",formula,-XMAX,XMAX,-XMAX,XMAX,-YMAX,YMAX)
 
TF3 * gd2 = new TF3("gd2",formula,-XMAX,XMAX,-XMAX,XMAX,-YMAX,YMAX)
 
TF3 * gd3 = new TF3("gd3",formula,-XMAX,XMAX,-XMAX,XMAX,-YMAX,YMAX)
 
double gd_phi [NEVENTS] = olongitude
 
double gd_rho [NEVENTS] = sqrt(xyz.Mag2())
 
double gd_theta [NEVENTS] = olatitude
 
double grad2rad = TMath::Pi()/180.
 
TH2D * h2gd = new TH2D("h2gd","h2gd", 100, -20, 20, 100, -20, 20)
 
double ilatitude = -(xyz.Theta()-TMath::Pi()/2.)*rad2grad
 
double ilongitude = xyz.Phi()*rad2grad
 
double olatitude
 
double olongitude
 
double rad2grad = 180./TMath::Pi()
 
double xgc =SOLAR_SISTEM_DISTANCE_FROM_GC
 
TVector3 xyz
 
double ygc =0.
 
double zgc =0.
 

Macro Definition Documentation

#define A1   5.81
#define A2   17.43
#define A3   34.86
#define B   0.3
#define DENOMINATOR   "(pow(((x-[3])*(x-[3])+y*y)+pow([0]+sqrt(z*z+[1]*[1]),2),5./2.)*pow(z*z+[1]*[1],3./2.))"
#define DRAW_H2
#define Md1   6.6e10
#define Md2   -2.9e10
#define Md3   3.3e9
#define MNGD_FILENAME   "MNGDDistribution.txt"
#define NEVENTS   1000000
#define NUMERATOR   "([0]*((x-[3])*(x-[3])+y*y)+([0]+3*sqrt(z*z+[1]*[1]))*pow([0]+sqrt(z*z+[1]*[1]),2))"
#define SOLAR_SISTEM_DISTANCE_FROM_GC   7.62
#define XMAX   40
#define YMAX   4

Function Documentation

canvas Clear ( )

Definition at line 1261 of file CWB_Plugin_xWRC.C.

h2gd Draw ( "colfz"  )
h2gd Fill ( x  ,
y   
)
GalacticToEquatorial ( ilongitude  ,
ilatitude  ,
olongitude  ,
olatitude   
)
cout<< "gc_phi : " << gc_phi << " gc_theta : " << gc_theta << " " << gc_rho << endl; for (int n=0;n<NEVENTS;n++) { if(n%10000==0) cout << n << endl; double x,y,z; gd-> GetRandom3 ( x  ,
y  ,
z   
)
h2gd GetXaxis ( ) -> SetNdivisions(-804)

Definition at line 321 of file cwb_report_prod_2.C.

h2gd GetYaxis ( ) -> SetNdivisions(-804)

Definition at line 213 of file cbc_plots.C.

h2gd GetZaxis ( ) -> SetLabelFont(42)
gStyle SetCanvasColor ( kWhite  )
canvas SetFillColor ( kWhite  )
canvas SetGridx ( )
canvas SetGridy ( )
gStyle SetLineColor ( kWhite  )
gStyle SetMarkerSize ( 2  )
gStyle SetMarkerStyle ( 7  )
gd SetNpx ( 100  )
gd SetNpy ( 100  )
gd SetNpz ( 100  )
gStyle SetNumberContours ( 256  )
gd1 SetParameter ( 0  ,
A1   
)
gd3 SetParameter ( 1  ,
B   
)
gd1 SetParameter ( 2  ,
Md1   
)
gd3 SetParameter ( 3  ,
SOLAR_SISTEM_DISTANCE_FROM_GC   
)
gd2 SetParameter ( 0  ,
A2   
)
gd2 SetParameter ( 2  ,
Md2   
)
gd3 SetParameter ( 0  ,
A3   
)
gd3 SetParameter ( 2  ,
Md3   
)
gStyle SetStatBorderSize ( 1  )
h2gd SetStats ( kFALSE  )
gStyle SetTextFont ( 12  )
h2gd SetTitle ( "h2gd"  )
gStyle SetTitleColor ( kBlue  ,
"D"   
)
gStyle SetTitleFillColor ( kWhite  )
gStyle SetTitleFont ( 12  ,
"D"   
)
gStyle SetTitleH ( 0.  050)
gStyle SetTitleW ( 0.  95)
gStyle SetTitleY ( 0.  98)
xyz SetXYZ ( xgc  ,
ygc  ,
zgc   
)
xyz SetXYZ ( x  ,
y  ,
z   
)
sprintf ( formula  ,
"***%s/%s"  [1][1][2],
NUMERATOR  ,
DENOMINATOR   
)
canvas ToggleEventStatus ( )

Variable Documentation

TCanvas* canvas = new TCanvas("AUNA", "GD", 300,40, 800, 800)

Definition at line 61 of file CreateMNGDDistribution.C.

char formula[256]

Definition at line 112 of file CreateMNGDDistribution.C.

double gc_phi = olongitude

Definition at line 169 of file CreateMNGDDistribution.C.

double gc_rho = sqrt(xyz.Mag2())

Definition at line 171 of file CreateMNGDDistribution.C.

double gc_theta = olatitude

Definition at line 170 of file CreateMNGDDistribution.C.

TF3* gd = new TF3("gd","gd1+gd2+gd3",-XMAX,XMAX,-XMAX,XMAX,-YMAX,YMAX)

Definition at line 133 of file CreateMNGDDistribution.C.

TF3* gd1 = new TF3("gd1",formula,-XMAX,XMAX,-XMAX,XMAX,-YMAX,YMAX)

Definition at line 115 of file CreateMNGDDistribution.C.

TF3* gd2 = new TF3("gd2",formula,-XMAX,XMAX,-XMAX,XMAX,-YMAX,YMAX)

Definition at line 121 of file CreateMNGDDistribution.C.

TF3* gd3 = new TF3("gd3",formula,-XMAX,XMAX,-XMAX,XMAX,-YMAX,YMAX)

Definition at line 127 of file CreateMNGDDistribution.C.

gd_phi[n] = olongitude

Definition at line 141 of file CreateMNGDDistribution.C.

gd_rho[n] = sqrt(xyz.Mag2())

Definition at line 143 of file CreateMNGDDistribution.C.

gd_theta[n] = olatitude

Definition at line 142 of file CreateMNGDDistribution.C.

double grad2rad = TMath::Pi()/180.

Definition at line 43 of file CreateMNGDDistribution.C.

TH2D* h2gd = new TH2D("h2gd","h2gd", 100, -20, 20, 100, -20, 20)

Definition at line 87 of file CreateMNGDDistribution.C.

ilatitude = -(xyz.Theta()-TMath::Pi()/2.)*rad2grad

Definition at line 160 of file CreateMNGDDistribution.C.

ilongitude = xyz.Phi()*rad2grad

Definition at line 159 of file CreateMNGDDistribution.C.

double olatitude

Definition at line 162 of file CreateMNGDDistribution.C.

double olongitude

Definition at line 161 of file CreateMNGDDistribution.C.

double rad2grad = 180./TMath::Pi()

Definition at line 44 of file CreateMNGDDistribution.C.

double xgc =SOLAR_SISTEM_DISTANCE_FROM_GC

Definition at line 156 of file CreateMNGDDistribution.C.

TVector3 xyz

Definition at line 155 of file CreateMNGDDistribution.C.

double ygc =0.

Definition at line 156 of file CreateMNGDDistribution.C.

double zgc =0.

Definition at line 156 of file CreateMNGDDistribution.C.