16 TFile
f(
"testem6_0.root");
20 TH1D*
h7 = (TH1D*)
f.Get(
"h7");
21 TH1D*
h8 = (TH1D*)
f.Get(
"h8");
22 TH1D*
h9 = (TH1D*)
f.Get(
"h9");
23 TH1D*
h10 = (TH1D*)
f.Get(
"h10");
24 TH1D*
h11 = (TH1D*)
f.Get(
"h11");
27 TH1D*
h12 = (TH1D*)
f.Get(
"h12");
28 TH1D*
h13 = (TH1D*)
f.Get(
"h13");
29 TH1D*
h14 = (TH1D*)
f.Get(
"h14");
30 TH1D*
h15 = (TH1D*)
f.Get(
"h15");
31 TH1D*
h16 = (TH1D*)
f.Get(
"h16");
34 TH1D*
h17 = (TH1D*)
f.Get(
"h17");
38 TCanvas*
c1 =
new TCanvas(
"c1",
" ");
42 h7->GetXaxis()->SetLimits(std::log10(minE),std::log10(maxE));
43 h7->SetTitle(
"Comparison of e+ Annihilation Processes(totcrsPerAtom)");
44 h7->GetXaxis()->SetTickLength(0);
45 h7->GetXaxis()->SetLabelOffset(999);
46 h7->GetYaxis()->SetTitle(
"microbarn");
48 h7->SetLineColor(kRed);
49 h9->SetLineColor(kBlue);
50 h11->SetMarkerColor(kRed);
51 h11->SetMarkerStyle(22);
52 h10->SetMarkerStyle(22);
56 h8->Draw(
"HIST SAME");
57 h9->Draw(
"HIST SAME");
58 h10->Draw(
"HIST P SAME");
59 h11->Draw(
"HIST P SAME");
61 TGaxis *
xaxis =
new TGaxis(std::log10(minE),minY,std::log10(maxE),minY,minE,maxE,510,
"G");
62 xaxis->SetTitle(
"GeV");
65 gStyle->SetOptStat(0);
67 TLegend*
leg =
new TLegend(0.78,0.59,0.98,0.76);
68 leg->AddEntry(h8,
"to 2 gammas",
"l");
69 leg->AddEntry(h7,
"to MuPair",
"l");
70 leg->AddEntry(h9,
"to Hadrons",
"l");
75 TCanvas*
c2 =
new TCanvas(
"c2",
" ");
79 h15->SetMinimum(minY1);
80 h15->SetMaximum(maxY1);
81 h15->GetXaxis()->SetLimits(std::log10(minE),std::log10(maxE));
82 h15->SetTitle(
"Comparison of EM Processes(totcrsPerVolume)");
83 h15->GetXaxis()->SetTickLength(0);
84 h15->GetXaxis()->SetLabelOffset(999);
85 h15->GetYaxis()->SetTitle(
"1/mm");
88 h13->SetLineColor(kGreen);
89 h14->SetLineColor(kRed);
90 h16->SetLineColor(kBlue);
94 h13->Draw(
"HIST SAME");
95 h14->Draw(
"HIST SAME");
96 h12->Draw(
"HIST SAME");
97 h16->Draw(
"HIST SAME");
99 TGaxis *
xaxis1 =
new TGaxis(std::log10(minE),minY1,std::log10(maxE),minY1,minE,maxE,510,
"G");
100 xaxis1->SetTitle(
"GeV");
103 gStyle->SetOptStat(0);
105 TLegend*
leg1 =
new TLegend(0.78,0.59,0.98,0.76);
106 leg1->AddEntry(h15,
"to 2 gammas",
"l");
107 leg1->AddEntry(h14,
"to MuPair",
"l");
108 leg1->AddEntry(h16,
"to Hadrons",
"l");
109 leg1->AddEntry(h12,
"Bremsstrahlung",
"l");
110 leg1->AddEntry(h13,
"Ionization",
"l");
TFile f("testem6_0.root")