8 // Calculates some geometrical properties of PbPb collisions
9 // in the Glauber Model
11 // Wood-Saxon nuclear density function
13 TCanvas *c1 = new TCanvas("c1","Wood Saxon",400,10,600,700);
14 TF1* f_wdsx = new TF1("f_wdsx", wdsx, 0, 15., 4);
15 f_wdsx->SetParameter(0,6.624);
16 f_wdsx->SetParameter(1,0.549);
17 f_wdsx->SetParameter(2,0.000);
18 f_wdsx->SetParameter(3,7.69e-4);
21 // Wood Saxon-nuclear density (b-z)
23 TCanvas *c2 = new TCanvas("c2","Wood Saxon",400,10,600,700);
24 TF2* f_wdsx_bz = new TF2("f_wdsx_bz", wdsx_bz, 0, 15., 0., 15., 4);
25 f_wdsx_bz->SetParameter(0,6.624);
26 f_wdsx_bz->SetParameter(1,0.549);
27 f_wdsx_bz->SetParameter(2,0.000);
28 f_wdsx_bz->SetParameter(3,7.69e-4);
31 // Wood Saxon-nuclear density (z, for fixed b)
33 TCanvas *c3 = new TCanvas("c3","Wood Saxon",400,10,600,700);
34 f_wdsx_z = new TF1("f_wdsx_z", wdsx_z, 0, 15., 5);
35 f_wdsx_z->SetParameter(0,6.624);
36 f_wdsx_z->SetParameter(1,0.549);
37 f_wdsx_z->SetParameter(2,0.000);
38 f_wdsx_z->SetParameter(3,7.69e-4);
39 f_wdsx_z->SetParameter(4,0.);
44 TCanvas *c4 = new TCanvas("c4","T_A",400,10,600,700);
45 f_ta = new TF1("f_ta", ta, 0, 15., 0);
48 // Kernel of overlap function
51 TCanvas *c5 = new TCanvas("c5","T_A",400,10,600,700);
52 f_ta_rfi = new TF2("f_ta_rfi", ta_rfi, 0, 15., 0., TMath::Pi(), 1);
53 f_ta_rfi->SetParameter(0,0.);
58 TCanvas *c6 = new TCanvas("c6","T_AA",400,10,600,700);
59 TF1* f_taa = new TF1("f_taa", taa,0.,15., 0);
66 Double_t taa(Double_t* x, Double_t* dum)
68 printf("taa %f\n", x[0]);
71 Double_t y = 10./(208*208)*hijing->Profile((Float_t)b);
76 Double_t wdsx(Double_t* x, Double_t* par)
79 // Wood Saxon Parameterisation
80 // as a function of radius
88 Double_t y = n * (1.+w*(xx/r0)*(xx/r0))/(1.+TMath::Exp((xx-r0)/d));
92 Double_t wdsx_bz(Double_t* x, Double_t* par)
95 // Wood Saxon Parameterisation
96 // as a function of z and b
100 Double_t r0 = par[0];
104 Double_t xx = TMath::Sqrt(bb*bb+zz*zz);
105 Double_t y = n * (1.+w*(xx/r0)*(xx/r0))/(1.+TMath::Exp((xx-r0)/d));
109 Double_t wdsx_z(Double_t* x, Double_t* par)
112 // Wood Saxon Parameterisation
113 // as a function of z for fixed b
115 Double_t bb = par[4];
117 Double_t r0 = par[0];
121 Double_t xx = TMath::Sqrt(bb*bb+zz*zz);
122 Double_t y = n * (1.+w*(xx/r0)*(xx/r0))/(1.+TMath::Exp((xx-r0)/d));
126 Double_t ta(Double_t* x, Double_t* par)
129 // Thickness function
132 f_wdsx_z->SetParameter(4,b);
133 Double_t y = 2.*f_wdsx_z->Integral(0.,15.);
137 Double_t ta_rfi(Double_t* x, Double_t* par)
140 // Kernel for overlap function
145 Double_t r2 = TMath::Sqrt(r1*r1+b*b-2.*r1*b*TMath::Cos(phi));
146 Double_t y = r1*f_ta->Eval(r1)*f_ta->Eval(r2);
147 // Double_t y = r1*f_ta->Eval(r1);
152 Double_t taa(Double_t* x, Double_t* par)
160 f_ta_rfi->SetParameter(0,b);
164 f_ta_rfi->Integral(0.,15., 0., TMath::Pi(), 0.001);
165 printf("taa %f %f\n", x[0], y);