12 gStyle->SetPalette(1);
13 gROOT->SetStyle(
"Plain");
15 c1 =
new TCanvas (
"c1",
"",20,20,1500,500);
18 system (
"rm -rf dna.root");
19 system (
"hadd dna.root dna_*.root");
24 ntuple = (TNtuple*)
f.Get(
"dna");
27 gStyle->SetOptStat(000000);
30 ntuple->SetFillStyle(1001);
31 ntuple->SetFillColor(2);
32 ntuple->Draw(
"flagProcess",
"",
"B");
35 ntuple->SetFillStyle(1001);
36 ntuple->SetFillColor(3);
37 ntuple->Draw(
"flagProcess",
"flagProcess==12||flagProcess==15||flagProcess==22||flagProcess==32||flagProcess==42||flagProcess==52||flagProcess==62",
"Bsame");
40 ntuple->SetFillStyle(1001);
41 ntuple->SetFillColor(4);
42 ntuple->Draw(
"flagProcess",
"flagProcess==11||flagProcess==21||flagProcess==31||flagProcess==41||flagProcess==51||flagProcess==61||flagProcess==110||flagProcess==210||flagProcess==410||flagProcess==510||flagProcess==710||flagProcess==120||flagProcess==220||flagProcess==420||flagProcess==520||flagProcess==720",
"Bsame");
45 ntuple->SetFillStyle(1001);
46 ntuple->SetFillColor(5);
47 ntuple->Draw(
"flagProcess",
"flagProcess==13||flagProcess==23||flagProcess==33||flagProcess==43||flagProcess==53||flagProcess==63||flagProcess==73||flagProcess==130||flagProcess==230||flagProcess==430||flagProcess==530||flagProcess==730",
"Bsame");
63 ntuple->SetMarkerColor(2);
65 ntuple->Draw(
"x:y:z",
"flagParticle==1");
88 ntuple->SetBranchAddress(
"flagParticle",&flagParticle);
89 ntuple->SetBranchAddress(
"flagProcess",&flagProcess);
90 ntuple->SetBranchAddress(
"x",&x);
91 ntuple->SetBranchAddress(
"y",&y);
92 ntuple->SetBranchAddress(
"z",&z);
93 ntuple->SetBranchAddress(
"totalEnergyDeposit",&totalEnergyDeposit);
94 ntuple->SetBranchAddress(
"stepLength",&stepLength);
95 ntuple->SetBranchAddress(
"kineticEnergyDifference",&kineticEnergyDifference);
96 ntuple->SetBranchAddress(
"kineticEnergy",&kineticEnergy);
97 ntuple->SetBranchAddress(
"cosTheta",&angle);
98 ntuple->SetBranchAddress(
"eventID",&eventID);
99 ntuple->SetBranchAddress(
"trackID",&trackID);
100 ntuple->SetBranchAddress(
"parentID",&parentID);
101 ntuple->SetBranchAddress(
"stepID",&stepID);
103 TH1F*
hsolvE =
new TH1F (
"hsolvE",
"solvE",100,0,2000);
104 TH1F*
helastE =
new TH1F (
"helastE",
"elastE",100,0,2000);
105 TH1F*
hexcitE =
new TH1F (
"hexcitE",
"excitE",100,0,2000);
106 TH1F*
hioniE =
new TH1F (
"hiioniE",
"ioniE",100,0,2000);
107 TH1F*
hattE =
new TH1F (
"hattE",
"attE",100,0,2000);
108 TH1F*
hvibE =
new TH1F (
"hvibE",
"vibE",100,0,2000);
110 for (
Int_t j=0;j<ntuple->GetEntries(); j++)
114 if (flagProcess==10) hsolvE->Fill(x);
115 if (flagProcess==11) helastE->Fill(x);
116 if (flagProcess==12) hexcitE->Fill(x);
117 if (flagProcess==13) hioniE->Fill(x);
118 if (flagProcess==14) hattE->Fill(x);
119 if (flagProcess==15) hvibE->Fill(x);
123 helastE->GetXaxis()->SetTitle(
"x (nm)");
124 helastE->SetLineColor(2);
126 hexcitE->SetLineColor(3);
127 hioniE->SetLineColor(4);
128 hattE->SetLineColor(5);
129 hvibE->SetLineColor(6);
130 hsolvE->SetLineColor(7);
135 hexcitE->Draw(
"SAME");
136 hioniE->Draw(
"SAME");
139 hsolvE->Draw(
"SAME");
system("rm -rf microbeam.root")
Double_t kineticEnergyDifference
Double_t totalEnergyDeposit