17 TFile *
experimentalFile =
new TFile(
"../../../experimentalData/proton/BraggPeak/62MeVInWater.root",
"READ");
28 TNtuple *
ntupleExperimental =
new TNtuple(
"ntupleExperimental",
"Protons, exp. data",
"depthExp:EdepExp");
35 ntupleExperimental ->
Fill(depthExp, EdepExp/normFactor);
42 TFile *
simulationFile =
new TFile(
"../../../SimulationOutputs/proton/BraggPeak/protonBraggPeak.root",
"READ");
48 simulatedPeak ->
Scale(1/simulationNormalisationFactor);
51 TCanvas *
c1 =
new TCanvas (
"c1",
"c1",200,10,600,400);
57 ntupleExperimental ->
Draw(
"EdepExp:depthExp");
58 simulatedPeak->
Draw(
"same");
61 leg =
new TLegend(0.50,0.60,0.20,0.70);
64 leg ->
AddEntry(ntupleExperimental,
"Experiment",
"P");
ntupleExperimental Draw("EdepExp:depthExp")
TNtupleSim Fill(f1, f2, f3, f4)
leg AddEntry(ntupleExperimental,"Experiment","P")
Float_t simulationNormalisationFactor
experimentalTree SetBranchAddress("EdepExp",&EdepExp)
cout<< "Opened file "<< fin<< " ixs= "<< ixs<< endl;if(ixs==0) hhh=(TH1F *) fff-> Get("h1")
ntupleExperimental SetMarkerStyle(4)
simulatedPeak SetMarkerSize(2)
TNtuple * ntupleExperimental
simulatedPeak Scale(1/simulationNormalisationFactor)
ntupleExperimental GetEntry(0)