LHAPDF veraion 5.9.1
[u/mrichter/AliRoot.git] / macros / printProfile.C
... / ...
CommitLineData
1void printProfile(char* g3 = "geant3/timing_200_woZDC.root", char* g4 ="geant4/timing_200_woZDC.root")
2{
3 // Create monitors
4 AliTransportMonitor* mon3 = AliTransportMonitor::Import(g3);
5 AliTransportMonitor* mon4 = AliTransportMonitor::Import(g4);
6 //
7 // Create timing histos and volume lists
8 mon3->Print();
9 tH3 = (TH1F*) volume_timing->Clone();
10 TObjArray* volumes3 = mon3->GetVolumes();
11 mon4->Print();
12 tH4 = (TH1F*) volume_timing->Clone();
13 TObjArray* volumes4 = mon4->GetVolumes();
14 //
15 // Normalise the geant 3 histo
16 // Fraction of time in %
17 Float_t ri = tH3->Integral();
18 tH3->Scale(100./ri);
19 //
20 Float_t cumulant = 0.; // cumulant
21 Int_t i = 1; // position in list
22 //
23 // print up to 90% of total time
24 //
25 printf(" Pos. Volume t (%) ct (%) Time_G3 Time_G4 perStepG3 perStepG4 StepsG3 StepsG4 te+e- tgam thad te+e- tgam thad\n");
26 while (cumulant < 90.) {
27 Float_t rt = tH3->GetBinContent(i);
28 cumulant += rt;
29 // extract the corresponding geant3 and geant4 volumes
30 char* volN = (tH3->GetXaxis())->GetBinLabel(i);
31 AliTransportMonitor::AliTransportMonitorVol* vol4 = (AliTransportMonitor::AliTransportMonitorVol*)
32 volumes4->FindObject(volN);
33 AliTransportMonitor::AliTransportMonitorVol* vol3 = 0;
34 if (i > 1) {
35 vol3 = (AliTransportMonitor::AliTransportMonitorVol*)
36 volumes3->FindObject(volN);
37 } else {
38 vol3 = (AliTransportMonitor::AliTransportMonitorVol*)
39 volumes3->At(1);
40 }
41 TString s(volN);
42 TString s20(s(0,20));
43 // Total time
44
45 Float_t t4[5];
46 t4[0] = vol4->GetTotalTime();
47 Float_t t3[5];
48 t3[0] = vol3->GetTotalTime();
49
50 // For e+/e-
51 t4[1] = ElectronTime(vol4);
52 t3[1] = ElectronTime(vol3);
53
54 // For photons
55 t4[2] = PhotonTime(vol4);
56 t3[2] = PhotonTime(vol3);
57
58 // Hadrons
59 t4[3] = HadronTime(vol4);
60 t3[3] = HadronTime(vol3);
61
62 printf("%5d %20s %6.2f %6.2f %8.3f %8.3f %8.3e %8.3e %8.3e %8.3e %8.3f %8.3f %8.3f %8.3f %8.3f %8.3f\n",
63 i,
64 s20.Data(),
65 rt, cumulant, t3[0], t4[0], t3[0]/(1.+vol3->GetNSteps()),
66 t4[0]/(1.+vol4->GetNSteps()), vol3->GetNSteps(), vol4->GetNSteps(),
67 t3[1]/t3[0] * 100., t3[2]/t3[0] * 100., t3[3]/t3[0] * 100.,
68 t4[1]/t4[0] * 100., t4[2]/t4[0] * 100., t4[3]/t4[0] * 100.
69
70 );
71 i++;
72 }
73}
74
75Float_t ElectronTime(AliTransportMonitor::AliTransportMonitorVol* vol)
76{
77 Int_t nt = vol->GetNtypes();
78 Float_t t = 0.;
79 for (Int_t i = 0; i < nt; i++) {
80 Int_t pdg = TMath::Abs(vol->GetPDG(i));
81 if (pdg == 11) t += vol->GetTime(i);
82 }
83
84 return t;
85}
86
87Float_t PhotonTime(AliTransportMonitor::AliTransportMonitorVol* vol)
88{
89 Int_t nt = vol->GetNtypes();
90 Float_t t = 0.;
91 for (Int_t i = 0; i < nt; i++) {
92 Int_t pdg = TMath::Abs(vol->GetPDG(i));
93 if (pdg == 22) t += vol->GetTime(i);
94 }
95 return t;
96}
97
98Float_t HadronTime(AliTransportMonitor::AliTransportMonitorVol* vol)
99{
100 Int_t nt = vol->GetNtypes();
101 Float_t t = 0.;
102 for (Int_t i = 0; i < nt; i++) {
103 Int_t pdg = TMath::Abs(vol->GetPDG(i));
104 if (pdg == 211 || pdg == 2212 || pdg == 2112 || pdg == 130 || pdg == 321) t += vol->GetTime(i);
105 }
106 return t;
107}
108
109