| 1 |
anderson |
1.1 |
{
|
| 2 |
|
|
//************************************************
|
| 3 |
|
|
// Variables
|
| 4 |
|
|
|
| 5 |
|
|
TFile *file1 = new TFile("PhotonJet20-200.root");
|
| 6 |
|
|
TFile *file2 = new TFile("JetET110.root");
|
| 7 |
|
|
|
| 8 |
|
|
TTree *tree1, *tree2;
|
| 9 |
|
|
file1->GetObject("TreePhotonAndJet",tree1);
|
| 10 |
|
|
file2->GetObject("TreePhotonAndJet",tree2);
|
| 11 |
|
|
|
| 12 |
|
|
TH1F* htemp1 = new TH1F("htemp1","blah",20, 0.0, 1.0);
|
| 13 |
|
|
TH1F* htemp2 = new TH1F("htemp2","blah",20, 0.0, 1.0);
|
| 14 |
|
|
|
| 15 |
|
|
TString thingToDraw1 = "recPhtnMtch.r19";
|
| 16 |
|
|
TString thingToDraw2 = thingToDraw1;
|
| 17 |
|
|
TString label1 = "#gamma";
|
| 18 |
|
|
TString label2 = "jet";
|
| 19 |
|
|
TString cutsAll = "recPhtnMtch.et>50&&recPhtnMtch.et<150&&abs(recPhtnMtch.eta)<1.5";
|
| 20 |
|
|
TString title = "R19";
|
| 21 |
|
|
TString titleX = "E(max) / E(3x3)";
|
| 22 |
|
|
|
| 23 |
|
|
TString saveFileName = title+".gif";
|
| 24 |
|
|
TString saveFileLocation = "/afs/cern.ch/user/a/anderson/www/";
|
| 25 |
|
|
|
| 26 |
|
|
//************************************************
|
| 27 |
|
|
|
| 28 |
|
|
|
| 29 |
|
|
// Draw from the TTrees
|
| 30 |
|
|
tree1->Draw(thingToDraw1+" >> htemp1", cutsAll);
|
| 31 |
|
|
tree2->Draw(thingToDraw2+" >> htemp2", cutsAll);
|
| 32 |
|
|
|
| 33 |
|
|
// Scale the histograms
|
| 34 |
|
|
Double_t scale1 = 1/htemp1->Integral();
|
| 35 |
|
|
Double_t scale2 = 1/htemp2->Integral();
|
| 36 |
|
|
htemp1->Scale(scale1);
|
| 37 |
|
|
htemp2->Scale(scale2);
|
| 38 |
|
|
|
| 39 |
|
|
// Set the histogram colors & lines
|
| 40 |
|
|
htemp1->SetLineColor(kRed);
|
| 41 |
|
|
htemp2->SetLineColor(kGreen);
|
| 42 |
|
|
htemp1->SetLineWidth(2);
|
| 43 |
|
|
htemp2->SetLineWidth(2);
|
| 44 |
|
|
|
| 45 |
|
|
// Turn off Stats Box
|
| 46 |
|
|
gStyle->SetOptStat(0);
|
| 47 |
|
|
|
| 48 |
|
|
// Create TStack but will later draw without stacking
|
| 49 |
|
|
THStack *tempStack = new THStack();
|
| 50 |
|
|
tempStack->Add(htemp1);
|
| 51 |
|
|
tempStack->Add(htemp2);
|
| 52 |
|
|
|
| 53 |
|
|
// Draw the histogram and titles
|
| 54 |
|
|
tempStack->Draw();
|
| 55 |
|
|
tempStack->SetTitle(title);
|
| 56 |
|
|
tempStack->GetXaxis()->SetTitle(titleX);
|
| 57 |
|
|
tempStack->Draw("nostack");
|
| 58 |
|
|
|
| 59 |
|
|
// Draw the legend
|
| 60 |
|
|
TLegend *infoBox = new TLegend(0.74, 0.74, 0.98, 0.98, "");
|
| 61 |
|
|
infoBox->AddEntry(htemp1,label1,"L");
|
| 62 |
|
|
infoBox->AddEntry(htemp2,label2,"L");
|
| 63 |
|
|
infoBox->Draw();
|
| 64 |
|
|
|
| 65 |
|
|
// Save the canvas
|
| 66 |
|
|
c1->SaveAs(saveFileLocation+saveFileName);
|
| 67 |
|
|
}
|
