dna/splitting/plot.C

0001 void plot() {
0002     Int_t numberOfSplit = 200;
0003
0004     TFile* ref = new TFile("splitting.root");
0005     TH1F* clusterSize = (TH1F*)ref->Get("2");
0006     TH1F* energyDeposit = (TH1F*)ref->Get("1");
0007
0008     // Re-normalization by the number of split to get absolute frequency histogram
0009     // Not necessary for energy deposit, as we scored edep * weight
0010     clusterSize->Scale(1.0/numberOfSplit);
0011
0012     TCanvas* c1 = new TCanvas("c1","",1200,550);
0013     c1->Divide(2,1);
0014     c1->cd(1);
0015     gPad->SetLogy();
0016     clusterSize->Draw();
0017     c1->cd(2);
0018     energyDeposit->Draw();
0019     c1->Update();
0020 }