From: cnattras Date: Sat, 21 Apr 2012 14:54:13 +0000 (+0000) Subject: Changes in efficiency to use ITS standalone tracks properly X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=8663c99d78339d3dbd820ed3e78abeaddeb8c28b;p=u%2Fmrichter%2FAliRoot.git Changes in efficiency to use ITS standalone tracks properly --- diff --git a/PWGLF/totEt/macros/hadEt/CorrEfficiency.C b/PWGLF/totEt/macros/hadEt/CorrEfficiency.C index 43ba119fc1b..6c65a1c353b 100644 --- a/PWGLF/totEt/macros/hadEt/CorrEfficiency.C +++ b/PWGLF/totEt/macros/hadEt/CorrEfficiency.C @@ -3,6 +3,8 @@ //Uses the output of AliAnalysisTaskHadEt //This is not actually what gets used in the correction class AliAnalysisHadEtCorrections - that is done in the macro GetCorrections.C - but this is useful for making plots and playing around with different options +TFile *file;// = new TFile(filename); + //this function calculates the efficiency propagating errors properly TH1D* bayneseffdiv(TH1D* numerator, TH1D* denominator,Char_t* name) { @@ -47,15 +49,17 @@ TH1D* bayneseffdiv(TH1D* numerator, TH1D* denominator,Char_t* name) //This is a somewhat messy function that gets the efficiency for different particles -TH1D *GetHisto(float cut = 0.12, char *name, int mycase, bool eta, int color, int marker,bool TPC, bool ITS, bool PbPb, int cb, int cblast){ +TH1D *GetHisto(float cut = 0.12, char *name, int mycase, bool eta, int color, int marker,bool TPC, bool ITS, bool PbPb, int cb, int cblast, char *filename){ //TFile *file = new TFile("Et.ESD.new.sim.merged.root"); - TFile *file = new TFile("rootFiles/LHC11a4_bis/Et.ESD.new.sim.LHC11a4_bis.root"); + //TFile *file = new TFile("rootFiles/LHC11a4_bis/Et.ESD.new.sim.LHC11a4_bis.root"); + file->cd(); //TFile *file = new TFile("rootFiles/LHC10d7/Et.ESD.new.sim.LHC10d7.root"); TList *list = file->FindObject("out2"); char *myname = "ITS"; - if(TPC&&!ITS) myname = "TPC"; - if(TPC&&ITS) myname = "TPCITS"; - cout<<"Using tracks from "<Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"AntiProton",i))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiPlus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiMinus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KMinus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KPlus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"Proton",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"AntiProton",i))); + } } } else{ @@ -77,6 +89,14 @@ TH1D *GetHisto(float cut = 0.12, char *name, int mycase, bool eta, int color, in numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"KPlus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"Proton"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"AntiProton"))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiPlus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiMinus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KMinus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KPlus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"Proton"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"AntiProton"))); + } } break; case 1://pion @@ -85,11 +105,19 @@ TH1D *GetHisto(float cut = 0.12, char *name, int mycase, bool eta, int color, in if(i==cb) numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiPlus",i)))->Clone("RecoPion"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiPlus",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiMinus",i))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiMinus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiPlus",i))); + } } } else{ numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"PiPlus")))->Clone("RecoPion"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"PiMinus"))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiPlus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiMinus"))); + } } break; case 2://kaon @@ -98,13 +126,21 @@ TH1D *GetHisto(float cut = 0.12, char *name, int mycase, bool eta, int color, in if(i==cb) numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i)))->Clone("RecoKaon"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KMinus",i))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KPlus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KMinus",i))); + } } } else{ - cout<<"I am kaoning here !"<FindObject(Form("EtNReconstructed%s%s",myname,"KPlus")))->Clone("RecoKaon"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"KMinus"))); - cout<<"Numerator "<GetEntries()<Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KPlus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KMinus"))); + } + //cout<<"Numerator "<GetEntries()<FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i)))->Clone("RecoProton"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"AntiProton",i))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"Proton",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"AntiProton",i))); + } } } else{ - cout<<"I am protoning here !"<FindObject(Form("EtNReconstructed%s%s",myname,"Proton")))->Clone("RecoProton"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"AntiProton"))); - cout<<"Numerator "<GetEntries()<Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"Proton"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"AntiProton"))); + } + //cout<<"Numerator "<GetEntries()<FindObject(Form("EtNReconstructed%s%s%s",myname,"EPlus",cbname)))->Clone("RecoElectron"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",myname,"EMinus",cbname))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",mynamealt,"EPlus",cbname))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",mynamealt,"EMinus",cbname))); + } break; } TH2F *denominatorParent; @@ -164,10 +212,10 @@ TH1D *GetHisto(float cut = 0.12, char *name, int mycase, bool eta, int color, in } } else{ - cout<<"I am here kaoning"<FindObject("EtNSimulatedKPlus"))->Clone("RecoKaon"); denominatorParent->Add((TH2F*) out2->FindObject("EtNSimulatedKMinus")); - cout<<"Denominator "<GetEntries()<GetEntries()<FindObject("EtNSimulatedProton"))->Clone("RecoProton"); denominatorParent->Add((TH2F*) out2->FindObject("EtNSimulatedAntiProton")); - cout<<"Denominator "<GetEntries()<GetEntries()<Add((TH2F*) out2->FindObject("EtNSimulatedEMinus")); break; } - cout<<"Numerator "<GetEntries()<GetEntries()<GetEntries()<GetEntries()<Sumw2(); denominatorParent->Sumw2(); TH1D *denominator; @@ -199,14 +247,14 @@ TH1D *GetHisto(float cut = 0.12, char *name, int mycase, bool eta, int color, in if(eta){ int lowbin = numeratorParent->GetYaxis()->FindBin(-cut+.001);//make sure we don't accv0entally get the wrong bin int highbin = numeratorParent->GetYaxis()->FindBin(cut-.001); - cout<<"Projecting from "<GetYaxis()->GetBinLowEdge(lowbin)<<" to "<GetYaxis()->GetBinLowEdge(highbin+1)<GetYaxis()->GetBinLowEdge(lowbin)<<" to "<GetYaxis()->GetBinLowEdge(highbin+1)<ProjectionX(Form("garbage%s",name),lowbin,highbin); numerator = numeratorParent->ProjectionX(name,lowbin,highbin); } else{ int lowbin = denominatorParent->GetXaxis()->FindBin(cut);//make sure we don't accidentally get the wrong bin int highbin = denominatorParent->GetXaxis()->GetNbins(); - cout<<"Here Projecting from "<GetXaxis()->GetBinLowEdge(lowbin)<<" to "<GetXaxis()->GetBinLowEdge(highbin+1)<GetXaxis()->GetBinLowEdge(lowbin)<<" to "<GetXaxis()->GetBinLowEdge(highbin+1)<ProjectionY(name,lowbin,highbin); denominator = denominatorParent->ProjectionY(Form("denominator%s",name),lowbin,highbin); } @@ -227,14 +275,20 @@ TH1D *GetHisto(float cut = 0.12, char *name, int mycase, bool eta, int color, in result->SetLineColor(color); result->SetMarkerStyle(marker); //result->Draw("e"); + //file->Close(); + //delete file; + delete numerator; + delete denominator; + result->SetName(name); return result; } //this is a method that makes pretty plots -void CorrEfficiency(char *shortprodname= "LHC11a4_bis", char *prodname = "LHC11a4_bis HIJING 2.76 TeV Pb+Pb", bool TPC = true,bool ITS = true, bool eta = true, bool PbPb = true, int cb = 0, int cblast = -1){ +void CorrEfficiency(char *prodname = "LHC11a4_bis HIJING 2.76 TeV Pb+Pb",char *shortprodname= "LHC11a4_bis", char *filename = "rootFiles/LHC11a4_bis/Et.ESD.new.sim.LHC11a4_bis.root", bool TPC = false,bool ITS = true, bool eta = true, bool PbPb = true, int cb = 0, int cblast = -1){ + file = new TFile(filename); gStyle->SetOptTitle(0); gStyle->SetOptStat(0); gStyle->SetOptFit(0); @@ -263,10 +317,10 @@ void CorrEfficiency(char *shortprodname= "LHC11a4_bis", char *prodname = "LHC11a phoscut = 0.1; emcalcut = 0.15; } - TH1D *PHOStotal = GetHisto(phoscut,"PHOStotal",0,eta,colortotal,phosmarker,TPC,ITS,PbPb,cb,cblast); - TH1D *PHOSpi = GetHisto(phoscut,"PHOSpi",1,eta,colorpi,phosmarker,TPC,ITS,PbPb,cb,cblast); - TH1D *PHOSp = GetHisto(phoscut,"PHOSp",3,eta,colork,phosmarker,TPC,ITS,PbPb,cb,cblast); - TH1D *PHOSk = GetHisto(phoscut,"PHOSk",2,eta,colorp,phosmarker,TPC,ITS,PbPb,cb,cblast); + TH1D *PHOStotal = GetHisto(phoscut,"PHOStotal",0,eta,colortotal,phosmarker,TPC,ITS,PbPb,cb,cblast,filename); + TH1D *PHOSpi = GetHisto(phoscut,"PHOSpi",1,eta,colorpi,phosmarker,TPC,ITS,PbPb,cb,cblast,filename); + TH1D *PHOSp = GetHisto(phoscut,"PHOSp",3,eta,colork,phosmarker,TPC,ITS,PbPb,cb,cblast,filename); + TH1D *PHOSk = GetHisto(phoscut,"PHOSk",2,eta,colorp,phosmarker,TPC,ITS,PbPb,cb,cblast,filename); if(eta) PHOStotal->GetXaxis()->SetRange(PHOStotal->GetXaxis()->FindBin(0.05),PHOStotal->GetXaxis()->FindBin(1.0)); //if(ITS&&!TPC){PHOStotal->GetXaxis()->SetRange(PHOStotal->GetXaxis()->FindBin(0.05),PHOStotal->GetXaxis()->FindBin(1.0));} //else{PHOStotal->GetXaxis()->SetRange(PHOStotal->GetXaxis()->FindBin(0.0),PHOStotal->GetXaxis()->FindBin(3.0));} @@ -283,14 +337,14 @@ void CorrEfficiency(char *shortprodname= "LHC11a4_bis", char *prodname = "LHC11a PHOSpi->Draw("same"); PHOSp->Draw("same"); PHOSk->Draw("same"); - cout<<"Hadrons"<Draw("same"); EMCALpi->Draw("same"); EMCALp->Draw("same"); @@ -339,12 +393,12 @@ void CorrEfficiency(char *shortprodname= "LHC11a4_bis", char *prodname = "LHC11a tex5->Draw(); } if(!TPC){ - c->SaveAs("pics/CorrEfficiency.eps"); - c->SaveAs("pics/CorrEfficiency.png"); + c->SaveAs(Form("pics/%s/CorrEfficiency.eps",shortprodname)); + c->SaveAs(Form("pics/%s/CorrEfficiency.png",shortprodname)); } else{ - c->SaveAs("pics/CorrEfficiencyTPCITS.eps"); - c->SaveAs("pics/CorrEfficiencyTPCITS.png"); + c->SaveAs(Form("pics/%s/CorrEfficiencyTPCITS.eps",shortprodname)); + c->SaveAs(Form("pics/%s/CorrEfficiencyTPCITS.png",shortprodname)); } if(PbPb){//make one more plot @@ -352,41 +406,84 @@ void CorrEfficiency(char *shortprodname= "LHC11a4_bis", char *prodname = "LHC11a //three centrality bins for efficiency 0-25%, 25-50%, 50-90% //use same for unidentified hadrons //kaons & protons - centrality dependence is more significant but I don't think I can do better on the binning - int pid = 0;//h=0,pi=1,p=3,k=2 - TCanvas *c2 = new TCanvas("c2","c2",600,400); - c2->SetTopMargin(0.02); - c2->SetRightMargin(0.02); - c2->SetBorderSize(0); - c2->SetFillColor(0); - c2->SetFillColor(0); - c2->SetBorderMode(0); - c2->SetFrameFillColor(0); - c2->SetFrameBorderMode(0); - cout<SetTopMargin(0.02); + c2->SetRightMargin(0.02); + c2->SetBorderSize(0); + c2->SetFillColor(0); + c2->SetFillColor(0); + c2->SetBorderMode(0); + c2->SetFrameFillColor(0); + c2->SetFrameBorderMode(0); + //cout<GetXaxis()->SetRange(cb0->GetXaxis()->FindBin(0.05),cb0->GetXaxis()->FindBin(1.0)); - cb0->SetMaximum(1.0); - cb0->Draw(); - cb4->Draw("same"); - // cb9->Draw("same"); - cb14->Draw("same"); - TLegend *leg = new TLegend(0.124161,0.747312,0.318792,0.959677);//(left,bottom,right,top) - leg->SetTextSize(0.0537634); - leg->SetBorderSize(0); - leg->SetLineColor(0); - leg->SetLineStyle(0); - leg->SetLineWidth(0); - leg->SetFillColor(0); - leg->SetFillStyle(0); - leg->AddEntry(cb0,"0-25%"); - leg->AddEntry(cb4,"25-50%"); - //leg->AddEntry(cb9,"45-50%"); - leg->AddEntry(cb14,"50-90%"); - leg->Draw(); - c2->SaveAs(Form("pics/CorrEfficiency%i.png",pid)); + TH1D *cb0 = GetHisto(phoscut,"cb0",pid,eta,1,20,TPC,ITS,PbPb,0,4,filename); + TH1D *cb4 = GetHisto(phoscut,"cb5",pid,eta,4,20,TPC,ITS,PbPb,5,9,filename); + //TH1D *cb9 = GetHisto(phoscut,"cb9",pid,eta,TColor::kGreen+4,20,TPC,ITS,PbPb,10,14); + TH1D *cb14 = GetHisto(phoscut,"cb14",pid,eta,2,20,TPC,ITS,PbPb,10,18,filename); + cb0->GetXaxis()->SetRange(cb0->GetXaxis()->FindBin(0.05),cb0->GetXaxis()->FindBin(5.0)); + cb0->SetMaximum(1.0); + cb0->Draw(); + cb4->Draw("same"); + // cb9->Draw("same"); + cb14->Draw("same"); + TLegend *leg = new TLegend(0.124161,0.747312,0.318792,0.959677);//(left,bottom,right,top) + leg->SetTextSize(0.0537634); + leg->SetBorderSize(0); + leg->SetLineColor(0); + leg->SetLineStyle(0); + leg->SetLineWidth(0); + leg->SetFillColor(0); + leg->SetFillStyle(0); + leg->AddEntry(cb0,"0-25%"); + leg->AddEntry(cb4,"25-50%"); + //leg->AddEntry(cb9,"45-50%"); + leg->AddEntry(cb14,"50-90%"); + leg->Draw(); + TString particle; + TString unidentified = "h"; + TString proton = "p"; + TString pion = "#pi"; + TString kaon = "K"; + switch(pid){ + case 0: + particle = unidentified; + break; + case 1: + particle = pion; + break; + case 2: + particle = kaon; + break; + case 3: + particle = proton; + break; + } + TLatex *texpid = new TLatex(0.25237,0.589541,particle.Data()); + texpid->SetTextSize(0.0806452); + texpid->Draw(); + switch(pid){ + case 0: + if(!TPC){c2->SaveAs(Form("pics/%s/CorrEfficiencyITSUnidentified.eps",shortprodname,pid));} + else{c2->SaveAs(Form("pics/%s/CorrEfficiencyTPCITSUnidentified.eps",shortprodname,pid));} + break; + case 1: + if(!TPC){c2->SaveAs(Form("pics/%s/CorrEfficiencyITSPion.eps",shortprodname,pid));} + else{c2->SaveAs(Form("pics/%s/CorrEfficiencyTPCITSPion.eps",shortprodname,pid));} + break; + case 3: + if(!TPC){c2->SaveAs(Form("pics/%s/CorrEfficiencyITSProton.eps",shortprodname,pid));} + else{c2->SaveAs(Form("pics/%s/CorrEfficiencyTPCITSProton.eps",shortprodname,pid));} + break; + case 2: + if(!TPC){c2->SaveAs(Form("pics/%s/CorrEfficiencyITSKaon.eps",shortprodname,pid));} + else{c2->SaveAs(Form("pics/%s/CorrEfficiencyTPCITSKaon.eps",shortprodname,pid));} + break; + } + } } } diff --git a/PWGLF/totEt/macros/hadEt/GetCorrections.C b/PWGLF/totEt/macros/hadEt/GetCorrections.C index 09b2a16daed..93b681ea168 100644 --- a/PWGLF/totEt/macros/hadEt/GetCorrections.C +++ b/PWGLF/totEt/macros/hadEt/GetCorrections.C @@ -10,6 +10,14 @@ // #include // #include "TStopwatch.h" +//Corrections added in by hand to deal with the inadequacies of PYTHIA and HIJING +TH1D *pp276TPCBkgd(); +TH1D *pp276ITSBkgd(); +TH1D *PbPb276TPCBkgd(); +TH1D *PbPb276ITSBkgd(); + Double_t xAxis1[112] = {0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1, 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.2, 0.22, 0.24, 0.26, 0.28, 0.3, 0.32, 0.34, 0.36, 0.38, 0.4, 0.42, 0.44, 0.46, 0.48, 0.5, 0.52, 0.54, 0.56, 0.58, 0.6, 0.62, 0.64, 0.66, 0.68, 0.7, 0.72, 0.74, 0.76, 0.78, 0.8, 0.82, 0.84, 0.86, 0.88, 0.9, 0.92, 0.94, 0.96, 0.98, 1, 1.05, 1.1, 1.15, 1.2, 1.25, 1.3, 1.35, 1.4, 1.45, 1.5, 1.55, 1.6, 1.65, 1.7, 1.75, 1.8, 1.85, 1.9, 1.95, 2, 2.2, 2.4, 2.6, 2.8, 3, 3.2, 3.4, 3.6, 3.8, 4, 4.2, 4.4, 4.6, 4.8, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10, 12, 14, 16, 18, 20, 25}; + + Float_t CorrNeutral(float ptcut, char *prodname, char *shortprodname, bool ispp = true, bool forSim = true, bool TPC, bool hadronic = false, float etacut = 0.7); TH1D *GetHistoCorrNeutral(float cut, char *name, bool ispp, bool forSim, int mycase, bool eta, int color, int marker, bool hadronic = false); @@ -160,14 +168,17 @@ void GetCorrections(char *prodname = "Enter Production Name", char *shortprodnam hadCorrectionEMCAL->SetNotIDConstCorrectionITSNoIDLowBound(1./NoIDITS*.98); hadCorrectionEMCAL->SetNotIDConstCorrectionTPCNoIDHighBound(1./NoIDTPC*1.02); hadCorrectionEMCAL->SetNotIDConstCorrectionITSNoIDHighBound(1./NoIDITS*1.02); - - TH1D *efficiencyPionTPC = GetHistoEfficiency(etacut,"hEfficiencyPionTPC",1,1,20,true,true); + + //Here we're going to do a bit of a cheat. We want the efficiency for ITS standalone tracks + TPC+ITS tracks. This is returned by the function for the efficiency function if I ask for ITS only efficiency. Had I known how this worked I probably would have written the code differently... but... + //anyhow I left a switch for changing it back. + bool useITSStandalone = true; + TH1D *efficiencyPionTPC = GetHistoEfficiency(etacut,"hEfficiencyPionTPC",1,1,20,useITSStandalone,true); hadCorrectionEMCAL->SetEfficiencyPionTPC(efficiencyPionTPC); - TH1D *efficiencyKaonTPC = GetHistoEfficiency(etacut,"hEfficiencyKaonTPC",2,1,20,true,true); + TH1D *efficiencyKaonTPC = GetHistoEfficiency(etacut,"hEfficiencyKaonTPC",2,1,20,useITSStandalone,true); hadCorrectionEMCAL->SetEfficiencyKaonTPC(efficiencyKaonTPC); - TH1D *efficiencyProtonTPC = GetHistoEfficiency(etacut,"hEfficiencyProtonTPC",3,1,20,true,true); + TH1D *efficiencyProtonTPC = GetHistoEfficiency(etacut,"hEfficiencyProtonTPC",3,1,20,useITSStandalone,true); hadCorrectionEMCAL->SetEfficiencyProtonTPC(efficiencyProtonTPC); - TH1D *efficiencyHadronTPC = GetHistoEfficiency(etacut,"hEfficiencyHadronTPC",0,1,20,true,true); + TH1D *efficiencyHadronTPC = GetHistoEfficiency(etacut,"hEfficiencyHadronTPC",0,1,20,useITSStandalone,true); hadCorrectionEMCAL->SetEfficiencyHadronTPC(efficiencyHadronTPC); TH1D *efficiencyPionITS = GetHistoEfficiency(etacut,"hEfficiencyPionITS",1,1,20,false,true); hadCorrectionEMCAL->SetEfficiencyPionITS(efficiencyPionITS); @@ -179,33 +190,33 @@ void GetCorrections(char *prodname = "Enter Production Name", char *shortprodnam hadCorrectionEMCAL->SetEfficiencyHadronITS(efficiencyHadronITS); if(!ispp){ - TH1D *efficiencyPionTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB0",1,1,20,true,true,0,4); + TH1D *efficiencyPionTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB0",1,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB0->Clone(Form("Test%i",i)),i); - TH1D *efficiencyPionTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB5",1,1,20,true,true,5,9); + TH1D *efficiencyPionTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB5",1,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB5->Clone(Form("Test%i",i)),i); - TH1D *efficiencyPionTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB10",1,1,20,true,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionEMCAL->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB10->Clone(Form("Test%i",i)),i); + TH1D *efficiencyPionTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB10",1,1,20,useITSStandalone,true,10,15); + for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB10->Clone(Form("Test%i",i)),i); - TH1D *efficiencyKaonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB0",2,1,20,true,true,0,4); + TH1D *efficiencyKaonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB0",2,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB0->Clone(Form("Test%i",i)),i); - TH1D *efficiencyKaonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB5",2,1,20,true,true,5,9); + TH1D *efficiencyKaonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB5",2,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB5->Clone(Form("Test%i",i)),i); - TH1D *efficiencyKaonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB10",2,1,20,true,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionEMCAL->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB10->Clone(Form("Test%i",i)),i);//Kaon + TH1D *efficiencyKaonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB10",2,1,20,useITSStandalone,true,10,15); + for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB10->Clone(Form("Test%i",i)),i);//Kaon - TH1D *efficiencyProtonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB0",3,1,20,true,true,0,4); + TH1D *efficiencyProtonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB0",3,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB0->Clone(Form("Test%i",i)),i); - TH1D *efficiencyProtonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB5",3,1,20,true,true,5,9); + TH1D *efficiencyProtonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB5",3,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB5->Clone(Form("Test%i",i)),i); - TH1D *efficiencyProtonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB10",3,1,20,true,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionEMCAL->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB10->Clone(Form("Test%i",i)),i);//Proton + TH1D *efficiencyProtonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB10",3,1,20,useITSStandalone,true,10,15); + for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB10->Clone(Form("Test%i",i)),i);//Proton - TH1D *efficiencyHadronTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB0",0,1,20,true,true,0,4); + TH1D *efficiencyHadronTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB0",0,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB0->Clone(Form("Test%i",i)),i); - TH1D *efficiencyHadronTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB5",0,1,20,true,true,5,9); + TH1D *efficiencyHadronTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB5",0,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB5->Clone(Form("Test%i",i)),i); - TH1D *efficiencyHadronTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB10",0,1,20,true,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionEMCAL->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB10->Clone(Form("Test%i",i)),i);//Hadron + TH1D *efficiencyHadronTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB10",0,1,20,useITSStandalone,true,10,15); + for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB10->Clone(Form("Test%i",i)),i);//Hadron TH1D *efficiencyPionITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB0",1,1,20,false,true,0,4); @@ -213,39 +224,78 @@ void GetCorrections(char *prodname = "Enter Production Name", char *shortprodnam TH1D *efficiencyPionITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB5",1,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB10",1,1,20,false,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionEMCAL->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB10->Clone(Form("Test%i",i)),i);//Pion + for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB10->Clone(Form("Test%i",i)),i);//Pion TH1D *efficiencyKaonITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB0",2,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB5",2,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB10",2,1,20,false,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionEMCAL->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB10->Clone(Form("Test%i",i)),i);//Kaon + for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB10->Clone(Form("Test%i",i)),i);//Kaon TH1D *efficiencyProtonITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB0",3,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB5",3,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB10",3,1,20,false,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionEMCAL->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB10->Clone(Form("Test%i",i)),i);//Proton + for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB10->Clone(Form("Test%i",i)),i);//Proton TH1D *efficiencyHadronITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB0",0,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionEMCAL->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB5",0,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionEMCAL->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB10",0,1,20,false,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionEMCAL->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB10->Clone(Form("Test%i",i)),i);//Hadron + for(int i=10;i<=20;i++) hadCorrectionEMCAL->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB10->Clone(Form("Test%i",i)),i);//Hadron }//EMCAL + hadCorrectionEMCAL->SetEfficiencyErrorLowBound(0.99); + hadCorrectionEMCAL->SetEfficiencyErrorHighBound(1.01); //CorrEfficiencyPlots(true,prodname,shortprodname); //CorrEfficiencyPlots(false,prodname,shortprodname,infilename); //hadCorrectionEMCAL->GetEfficiencyHadronTPC()->Draw(); - TH1D *backgroundTPC = GetHistoCorrBkgd(etacut,"hBackgroundTPC",true,ispp,forSim); - TH1D *backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); + TH1D *backgroundTPC; + TH1D *backgroundITS; + if((dataset==20111 || dataset==20100) && !forSim){//2.76 TeV p+p or Pb+Pb + if(dataset==20111){ + cout<<"Fixing 2.76 TeV p+p background to be average of 900 GeV and 7 TeV scaling"<SetName("hBackgroundTPC"); + backgroundITS = pp276ITSBkgd(); + backgroundITS->SetName("hBackgroundITS"); + } + else{//PbPb + cout<<"Fixing 2.76 TeV Pb+Pb background to be average of 900 GeV and 7 TeV scaling with baryon enhancement"<SetName("hBackgroundTPC"); + //ITS background is currently a placeholder + backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); + } + } + else{ + backgroundTPC = GetHistoCorrBkgd(etacut,"hBackgroundTPC",true,ispp,forSim); + backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); + } + float bkgdpcterror = 0.0; + switch(dataset){ + case 2009: + bkgdpcterror = 0.37; + break; + case 20111: + bkgdpcterror = 0.38; + break; + case 2010: + bkgdpcterror = 0.13; + break; + case 20100: + bkgdpcterror = 0.76; + break; + } hadCorrectionEMCAL->SetBackgroundCorrectionTPC(backgroundTPC); hadCorrectionEMCAL->SetBackgroundCorrectionITS(backgroundITS); + hadCorrectionEMCAL->SetBackgroundErrorLowBound(1.0-bkgdpcterror/100.0); + hadCorrectionEMCAL->SetBackgroundErrorHighBound(1.0+bkgdpcterror/100.0); //CorrBkgdPlots(prodname,shortprodname,true,ispp,forSim); //CorrBkgdPlots(prodname,shortprodname,false,ispp,forSim); @@ -346,10 +396,10 @@ void GetCorrections(char *prodname = "Enter Production Name", char *shortprodnam hadCorrectionPHOS->SetNotIDConstCorrectionTPCNoIDHighBound(1./NoIDTPC*1.02); hadCorrectionPHOS->SetNotIDConstCorrectionITSNoIDHighBound(1./NoIDITS*1.02); - TH1D *efficiencyPionTPC = GetHistoEfficiency(etacut,"hEfficiencyPionTPC",1,1,20,true,true); - TH1D *efficiencyKaonTPC = GetHistoEfficiency(etacut,"hEfficiencyKaonTPC",2,1,20,true,true); - TH1D *efficiencyProtonTPC = GetHistoEfficiency(etacut,"hEfficiencyProtonTPC",3,1,20,true,true); - TH1D *efficiencyHadronTPC = GetHistoEfficiency(etacut,"hEfficiencyHadronTPC",0,1,20,true,true); + TH1D *efficiencyPionTPC = GetHistoEfficiency(etacut,"hEfficiencyPionTPC",1,1,20,useITSStandalone,true); + TH1D *efficiencyKaonTPC = GetHistoEfficiency(etacut,"hEfficiencyKaonTPC",2,1,20,useITSStandalone,true); + TH1D *efficiencyProtonTPC = GetHistoEfficiency(etacut,"hEfficiencyProtonTPC",3,1,20,useITSStandalone,true); + TH1D *efficiencyHadronTPC = GetHistoEfficiency(etacut,"hEfficiencyHadronTPC",0,1,20,useITSStandalone,true); TH1D *efficiencyPionITS = GetHistoEfficiency(etacut,"hEfficiencyPionITS",1,1,20,false,true); TH1D *efficiencyKaonITS = GetHistoEfficiency(etacut,"hEfficiencyKaonITS",2,1,20,false,true); TH1D *efficiencyProtonITS = GetHistoEfficiency(etacut,"hEfficiencyProtonITS",3,1,20,false,true); @@ -367,33 +417,33 @@ void GetCorrections(char *prodname = "Enter Production Name", char *shortprodnam if(!ispp){ - TH1D *efficiencyPionTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB0",1,1,20,true,true,0,4); + TH1D *efficiencyPionTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB0",1,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB0->Clone(Form("Test%i",i)),i); - TH1D *efficiencyPionTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB5",1,1,20,true,true,5,9); + TH1D *efficiencyPionTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB5",1,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB5->Clone(Form("Test%i",i)),i); - TH1D *efficiencyPionTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB10",1,1,20,true,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionPHOS->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB10->Clone(Form("Test%i",i)),i); + TH1D *efficiencyPionTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionTPCCB10",1,1,20,useITSStandalone,true,10,15); + for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyPionTPC((TH1D*)efficiencyPionTPCCB10->Clone(Form("Test%i",i)),i); - TH1D *efficiencyKaonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB0",2,1,20,true,true,0,4); + TH1D *efficiencyKaonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB0",2,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB0->Clone(Form("Test%i",i)),i); - TH1D *efficiencyKaonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB5",2,1,20,true,true,5,9); + TH1D *efficiencyKaonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB5",2,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB5->Clone(Form("Test%i",i)),i); - TH1D *efficiencyKaonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB10",2,1,20,true,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionPHOS->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB10->Clone(Form("Test%i",i)),i);//Kaon + TH1D *efficiencyKaonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonTPCCB10",2,1,20,useITSStandalone,true,10,15); + for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyKaonTPC((TH1D*)efficiencyKaonTPCCB10->Clone(Form("Test%i",i)),i);//Kaon - TH1D *efficiencyProtonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB0",3,1,20,true,true,0,4); + TH1D *efficiencyProtonTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB0",3,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB0->Clone(Form("Test%i",i)),i); - TH1D *efficiencyProtonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB5",3,1,20,true,true,5,9); + TH1D *efficiencyProtonTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB5",3,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB5->Clone(Form("Test%i",i)),i); - TH1D *efficiencyProtonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB10",3,1,20,true,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionPHOS->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB10->Clone(Form("Test%i",i)),i);//Proton + TH1D *efficiencyProtonTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonTPCCB10",3,1,20,useITSStandalone,true,10,15); + for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyProtonTPC((TH1D*)efficiencyProtonTPCCB10->Clone(Form("Test%i",i)),i);//Proton - TH1D *efficiencyHadronTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB0",0,1,20,true,true,0,4); + TH1D *efficiencyHadronTPCCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB0",0,1,20,useITSStandalone,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB0->Clone(Form("Test%i",i)),i); - TH1D *efficiencyHadronTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB5",0,1,20,true,true,5,9); + TH1D *efficiencyHadronTPCCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB5",0,1,20,useITSStandalone,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB5->Clone(Form("Test%i",i)),i); - TH1D *efficiencyHadronTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB10",0,1,20,true,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionPHOS->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB10->Clone(Form("Test%i",i)),i);//Hadron + TH1D *efficiencyHadronTPCCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronTPCCB10",0,1,20,useITSStandalone,true,10,15); + for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyHadronTPC((TH1D*)efficiencyHadronTPCCB10->Clone(Form("Test%i",i)),i);//Hadron TH1D *efficiencyPionITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB0",1,1,20,false,true,0,4); @@ -401,36 +451,59 @@ void GetCorrections(char *prodname = "Enter Production Name", char *shortprodnam TH1D *efficiencyPionITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB5",1,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyPionITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyPionITSCB10",1,1,20,false,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionPHOS->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB10->Clone(Form("Test%i",i)),i);//Pion + for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyPionITS((TH1D*)efficiencyPionITSCB10->Clone(Form("Test%i",i)),i);//Pion TH1D *efficiencyKaonITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB0",2,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB5",2,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyKaonITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyKaonITSCB10",2,1,20,false,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionPHOS->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB10->Clone(Form("Test%i",i)),i);//Kaon + for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyKaonITS((TH1D*)efficiencyKaonITSCB10->Clone(Form("Test%i",i)),i);//Kaon TH1D *efficiencyProtonITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB0",3,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB5",3,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyProtonITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyProtonITSCB10",3,1,20,false,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionPHOS->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB10->Clone(Form("Test%i",i)),i);//Proton + for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyProtonITS((TH1D*)efficiencyProtonITSCB10->Clone(Form("Test%i",i)),i);//Proton TH1D *efficiencyHadronITSCB0 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB0",0,1,20,false,true,0,4); for(int i=0;i<=4;i++) hadCorrectionPHOS->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB0->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronITSCB5 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB5",0,1,20,false,true,5,9); for(int i=5;i<=9;i++) hadCorrectionPHOS->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB5->Clone(Form("Test%i",i)),i); TH1D *efficiencyHadronITSCB10 = GetHistoEfficiency(etacut,"hEfficiencyHadronITSCB10",0,1,20,false,true,10,15); - for(int i=10;i<=19;i++) hadCorrectionPHOS->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB10->Clone(Form("Test%i",i)),i);//Hadron + for(int i=10;i<=20;i++) hadCorrectionPHOS->SetEfficiencyHadronITS((TH1D*)efficiencyHadronITSCB10->Clone(Form("Test%i",i)),i);//Hadron }//EMCAL - TH1D *backgroundTPC = GetHistoCorrBkgd(etacut,"hBackgroundTPC",true,ispp,forSim); - TH1D *backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); + + TH1D *backgroundTPC; + TH1D *backgroundITS; + if((dataset==20111 || dataset==20100) && !forSim){//2.76 TeV p+p or Pb+Pb + if(dataset==20111){ + cout<<"Fixing 2.76 TeV p+p background to be average of 900 GeV and 7 TeV scaling"<SetName("hBackgroundTPCPHOS"); + backgroundITS = pp276ITSBkgd(); + backgroundITS->SetName("hBackgroundITSPHOS"); + } + else{//PbPb + cout<<"Fixing 2.76 TeV Pb+Pb background to be average of 900 GeV and 7 TeV scaling with baryon enhancement"<SetName("hBackgroundTPCPHOS"); + //ITS background is currently a placeholder + backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); + } + } + else{ + backgroundTPC = GetHistoCorrBkgd(etacut,"hBackgroundTPC",true,ispp,forSim); + backgroundITS = GetHistoCorrBkgd(etacut,"hBackgroundITS",false,ispp,forSim); + } hadCorrectionPHOS->SetBackgroundCorrectionTPC(backgroundTPC); hadCorrectionPHOS->SetBackgroundCorrectionITS(backgroundITS); - CorrBkgdPlots(prodname,shortprodname,true,ispp,forSim); - CorrBkgdPlots(prodname,shortprodname,false,ispp,forSim); + hadCorrectionPHOS->SetBackgroundErrorLowBound(1.0-0.001); + hadCorrectionPHOS->SetBackgroundErrorHighBound(1.0+0.001); + //CorrBkgdPlots(prodname,shortprodname,true,ispp,forSim); + //CorrBkgdPlots(prodname,shortprodname,false,ispp,forSim); hadCorrectionPHOS->Report(); //Write the output @@ -877,7 +950,7 @@ Float_t CorrPtCut(float ptcut, char *prodname, char *shortprodname, bool ispp, b TH1D *GetHistoCorrNotID(float etacut,char *name, bool TPC, bool eta, bool ispp, bool forSim){ file->cd(); char *myname = mynameITS; - if(TPC) myname = mynameTPC; + if(TPC) myname = mynameTPCITS; TH2F *notid = ((TH2F*) out2->FindObject(Form("EtReconstructed%sUnidentifiedAssumingPion",myname)))->Clone("notid"); TH2F *nNotid = ((TH2F*) out2->FindObject(Form("EtNReconstructed%sUnidentified",myname)))->Clone("nNotid"); if(!eta){ @@ -1077,7 +1150,7 @@ Float_t CorrNotIDConst(float ptcut, float etacut,char *name, char *prodname, cha TH1D *GetHistoNoID(float etacut, char *name, bool eta, bool TPC, bool ispp, bool forSim){ file->cd(); char *myname = mynameITS; - if(TPC) myname = mynameTPC; + if(TPC) myname = mynameTPCITS; TH2F *notid = ((TH2F*) out2->FindObject(Form("EtReconstructed%sChargedHadronAssumingPion",myname)))->Clone("notid"); TH2F *nNotid = ((TH2F*) out2->FindObject(Form("EtNReconstructed%sChargedHadron",myname)))->Clone("nNotid"); @@ -1285,8 +1358,9 @@ TH1D *GetHistoEfficiency(float cut, char *name, int mycase, int color, int marke bool eta = true; file->cd(); char *myname = mynameITS; - if(TPC&&!ITS) myname = mynameTPC; + if(TPC && !ITS){ myname = mynameTPCITS;cout<<"I was here I should not be called"<Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"AntiProton",i))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiPlus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiMinus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KMinus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KPlus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"Proton",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"AntiProton",i))); + } } } else{ @@ -1308,6 +1390,14 @@ TH1D *GetHistoEfficiency(float cut, char *name, int mycase, int color, int marke numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"KPlus"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"Proton"))); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"AntiProton"))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiPlus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiMinus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KMinus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KPlus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"Proton"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"AntiProton"))); + } } break; case 1://pion @@ -1316,11 +1406,19 @@ TH1D *GetHistoEfficiency(float cut, char *name, int mycase, int color, int marke if(i==cb) numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiPlus",i)))->Clone("RecoPion"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiPlus",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"PiMinus",i))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiMinus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"PiPlus",i))); + } } } else{ numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"PiPlus")))->Clone("RecoPion"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"PiMinus"))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiPlus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"PiMinus"))); + } } break; case 2://kaon @@ -1329,11 +1427,21 @@ TH1D *GetHistoEfficiency(float cut, char *name, int mycase, int color, int marke if(i==cb) numeratorParent = (TH2F*)((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i)))->Clone("RecoKaon"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KPlus",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"KMinus",i))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KPlus",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"KMinus",i))); + } } } else{ + // cout<<"I am kaoning here !"<FindObject(Form("EtNReconstructed%s%s",myname,"KPlus")))->Clone("RecoKaon"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"KMinus"))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KPlus"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"KMinus"))); + } + //cout<<"Numerator "<GetEntries()<FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i)))->Clone("RecoProton"); else{numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"Proton",i)));} numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",myname,"AntiProton",i))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"Proton",i))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%sCB%i",mynamealt,"AntiProton",i))); + } } } else{ + //cout<<"I am protoning here !"<FindObject(Form("EtNReconstructed%s%s",myname,"Proton")))->Clone("RecoProton"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",myname,"AntiProton"))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"Proton"))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s",mynamealt,"AntiProton"))); + } + //cout<<"Numerator "<GetEntries()<FindObject(Form("EtNReconstructed%s%s%s",myname,"EPlus",cbname)))->Clone("RecoElectron"); numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",myname,"EMinus",cbname))); + if(!TPC&&ITS){//ITS Standalone tracks - from leftover hits, don't want to double count + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",mynamealt,"EPlus",cbname))); + numeratorParent->Add((TH2F*) out2->FindObject(Form("EtNReconstructed%s%s%s",mynamealt,"EMinus",cbname))); + } break; } TH2F *denominatorParent; @@ -1567,7 +1689,7 @@ TH1D *GetHistoCorrBkgd(float etacut,char *name, bool TPC,bool ispp,bool forSim){ char *reweightname = reweightedNo; if(!forSim) reweightname = reweightedYes; char *myname = mynameITS; - if(TPC) myname = mynameTPC; + if(TPC) myname = mynameTPCITS; TH2F *signal = ((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiPlus",myname)))->Clone("signal"); signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedPiMinus",myname))); signal->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sIdentifiedKMinus",myname))); @@ -1585,6 +1707,10 @@ TH1D *GetHistoCorrBkgd(float etacut,char *name, bool TPC,bool ispp,bool forSim){ bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiXiDaughters",myname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sOmegaDaughters",myname))); bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sAntiOmegaDaughters",myname))); + bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sConversionElectrons",myname)) ); + bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryMuons",myname)) ); + bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryPions",myname)) ); + bkgd->Add((TH2F*) out2->FindObject(Form("EtReconstructed%sSecondaryProtons",myname)) ); int lowbin = bkgd->GetYaxis()->FindBin(-etacut+.001);//make sure we don't accidentally get the wrong bin int highbin = bkgd->GetYaxis()->FindBin(etacut-.001); //cout<<"Projecting from "<GetYaxis()->GetBinLowEdge(lowbin)<<" to "<GetYaxis()->GetBinLowEdge(highbin+1)<ProjectionX(name,lowbin,highbin); TH1D *numerator = bkgd->ProjectionX("numerator",lowbin,highbin); + denominator->Add(numerator); numerator->Divide(denominator); numerator->SetYTitle("Ratio of E_{T}^{background}/E_{T}^{real}"); numerator->GetYaxis()->SetTitleOffset(1.2); @@ -1665,3 +1792,1099 @@ void CorrBkgdPlots(char *prodname, char *shortprodname, bool TPC,bool ispp,bool delete leg; } + +//TH1D *ppBkgd() +TH1D *pp276TPCBkgd(){ +//=========Macro generated from canvas: c4/c4 +//========= (Thu May 5 23:07:14 2011) by ROOT version5.28/00c + TH1D *Allpt7TeVScaling = new TH1D("Allpt7TeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); + Allpt7TeVScaling->SetBinContent(13,0.1673787); + Allpt7TeVScaling->SetBinContent(14,0.1548408); + Allpt7TeVScaling->SetBinContent(15,0.1328529); + Allpt7TeVScaling->SetBinContent(16,0.1220669); + Allpt7TeVScaling->SetBinContent(17,0.1109226); + Allpt7TeVScaling->SetBinContent(18,0.0981996); + Allpt7TeVScaling->SetBinContent(19,0.09152354); + Allpt7TeVScaling->SetBinContent(20,0.07937801); + Allpt7TeVScaling->SetBinContent(21,0.07163593); + Allpt7TeVScaling->SetBinContent(22,0.0562236); + Allpt7TeVScaling->SetBinContent(23,0.04398756); + Allpt7TeVScaling->SetBinContent(24,0.039459); + Allpt7TeVScaling->SetBinContent(25,0.03618901); + Allpt7TeVScaling->SetBinContent(26,0.0322147); + Allpt7TeVScaling->SetBinContent(27,0.03294796); + Allpt7TeVScaling->SetBinContent(28,0.03114461); + Allpt7TeVScaling->SetBinContent(29,0.03039579); + Allpt7TeVScaling->SetBinContent(30,0.03261539); + Allpt7TeVScaling->SetBinContent(31,0.03239832); + Allpt7TeVScaling->SetBinContent(32,0.03055859); + Allpt7TeVScaling->SetBinContent(33,0.03181271); + Allpt7TeVScaling->SetBinContent(34,0.02882665); + Allpt7TeVScaling->SetBinContent(35,0.03102149); + Allpt7TeVScaling->SetBinContent(36,0.0272696); + Allpt7TeVScaling->SetBinContent(37,0.02783727); + Allpt7TeVScaling->SetBinContent(38,0.02334413); + Allpt7TeVScaling->SetBinContent(39,0.02357657); + Allpt7TeVScaling->SetBinContent(40,0.02630686); + Allpt7TeVScaling->SetBinContent(41,0.02372341); + Allpt7TeVScaling->SetBinContent(42,0.02279457); + Allpt7TeVScaling->SetBinContent(43,0.01999143); + Allpt7TeVScaling->SetBinContent(44,0.02019579); + Allpt7TeVScaling->SetBinContent(45,0.02029719); + Allpt7TeVScaling->SetBinContent(46,0.0166231); + Allpt7TeVScaling->SetBinContent(47,0.01824842); + Allpt7TeVScaling->SetBinContent(48,0.0188995); + Allpt7TeVScaling->SetBinContent(49,0.01836067); + Allpt7TeVScaling->SetBinContent(50,0.02113978); + Allpt7TeVScaling->SetBinContent(51,0.02236854); + Allpt7TeVScaling->SetBinContent(52,0.0204874); + Allpt7TeVScaling->SetBinContent(53,0.02236788); + Allpt7TeVScaling->SetBinContent(54,0.01862457); + Allpt7TeVScaling->SetBinContent(55,0.0215448); + Allpt7TeVScaling->SetBinContent(56,0.01757402); + Allpt7TeVScaling->SetBinContent(57,0.01899477); + Allpt7TeVScaling->SetBinContent(58,0.01766989); + Allpt7TeVScaling->SetBinContent(59,0.01984435); + Allpt7TeVScaling->SetBinContent(60,0.01842092); + Allpt7TeVScaling->SetBinContent(61,0.01668053); + Allpt7TeVScaling->SetBinContent(62,0.01708603); + Allpt7TeVScaling->SetBinContent(63,0.01557306); + Allpt7TeVScaling->SetBinContent(64,0.01878942); + Allpt7TeVScaling->SetBinContent(65,0.01458043); + Allpt7TeVScaling->SetBinContent(66,0.0151346); + Allpt7TeVScaling->SetBinContent(67,0.01565407); + Allpt7TeVScaling->SetBinContent(68,0.01327014); + Allpt7TeVScaling->SetBinContent(69,0.01284378); + Allpt7TeVScaling->SetBinContent(70,0.01046102); + Allpt7TeVScaling->SetBinContent(71,0.009963117); + Allpt7TeVScaling->SetBinContent(72,0.01359609); + Allpt7TeVScaling->SetBinContent(73,0.01088263); + Allpt7TeVScaling->SetBinContent(74,0.01317689); + Allpt7TeVScaling->SetBinContent(75,0.01433201); + Allpt7TeVScaling->SetBinContent(76,0.01215939); + Allpt7TeVScaling->SetBinContent(77,0.01380766); + Allpt7TeVScaling->SetBinContent(78,0.01499505); + Allpt7TeVScaling->SetBinContent(79,0.015013); + Allpt7TeVScaling->SetBinContent(80,0.01140629); + Allpt7TeVScaling->SetBinContent(81,0.009461979); + Allpt7TeVScaling->SetBinContent(82,0.007775343); + Allpt7TeVScaling->SetBinContent(83,0.01092198); + Allpt7TeVScaling->SetBinContent(84,0.009924767); + Allpt7TeVScaling->SetBinContent(85,0.006399066); + Allpt7TeVScaling->SetBinContent(86,0.008370865); + Allpt7TeVScaling->SetBinContent(87,0.01115192); + Allpt7TeVScaling->SetBinContent(88,0.003781796); + Allpt7TeVScaling->SetBinContent(89,0.01596916); + Allpt7TeVScaling->SetBinContent(90,0.005959723); + Allpt7TeVScaling->SetBinContent(91,0.006885685); + Allpt7TeVScaling->SetBinContent(93,0.006011974); + Allpt7TeVScaling->SetBinContent(94,0.008862645); + Allpt7TeVScaling->SetBinContent(95,0.009963223); + Allpt7TeVScaling->SetBinContent(96,0.01457121); + Allpt7TeVScaling->SetBinContent(97,0.007500135); + Allpt7TeVScaling->SetBinContent(98,0.02636933); + Allpt7TeVScaling->SetBinContent(103,0.06353767); + Allpt7TeVScaling->SetBinError(13,0.02859277); + Allpt7TeVScaling->SetBinError(14,0.006796474); + Allpt7TeVScaling->SetBinError(15,0.004234415); + Allpt7TeVScaling->SetBinError(16,0.0036414); + Allpt7TeVScaling->SetBinError(17,0.003333835); + Allpt7TeVScaling->SetBinError(18,0.003009034); + Allpt7TeVScaling->SetBinError(19,0.002839933); + Allpt7TeVScaling->SetBinError(20,0.002588262); + Allpt7TeVScaling->SetBinError(21,0.001703954); + Allpt7TeVScaling->SetBinError(22,0.001495734); + Allpt7TeVScaling->SetBinError(23,0.001297475); + Allpt7TeVScaling->SetBinError(24,0.001213434); + Allpt7TeVScaling->SetBinError(25,0.001162398); + Allpt7TeVScaling->SetBinError(26,0.001137158); + Allpt7TeVScaling->SetBinError(27,0.001208358); + Allpt7TeVScaling->SetBinError(28,0.001183675); + Allpt7TeVScaling->SetBinError(29,0.001263691); + Allpt7TeVScaling->SetBinError(30,0.001320041); + Allpt7TeVScaling->SetBinError(31,0.001362399); + Allpt7TeVScaling->SetBinError(32,0.001338259); + Allpt7TeVScaling->SetBinError(33,0.001369813); + Allpt7TeVScaling->SetBinError(34,0.001358652); + Allpt7TeVScaling->SetBinError(35,0.00142864); + Allpt7TeVScaling->SetBinError(36,0.001378454); + Allpt7TeVScaling->SetBinError(37,0.001384729); + Allpt7TeVScaling->SetBinError(38,0.001257903); + Allpt7TeVScaling->SetBinError(39,0.001295459); + Allpt7TeVScaling->SetBinError(40,0.001546992); + Allpt7TeVScaling->SetBinError(41,0.001496938); + Allpt7TeVScaling->SetBinError(42,0.001486813); + Allpt7TeVScaling->SetBinError(43,0.00140119); + Allpt7TeVScaling->SetBinError(44,0.001478331); + Allpt7TeVScaling->SetBinError(45,0.001566594); + Allpt7TeVScaling->SetBinError(46,0.001407568); + Allpt7TeVScaling->SetBinError(47,0.001509742); + Allpt7TeVScaling->SetBinError(48,0.001583129); + Allpt7TeVScaling->SetBinError(49,0.001668789); + Allpt7TeVScaling->SetBinError(50,0.001915902); + Allpt7TeVScaling->SetBinError(51,0.001968585); + Allpt7TeVScaling->SetBinError(52,0.00195473); + Allpt7TeVScaling->SetBinError(53,0.00207896); + Allpt7TeVScaling->SetBinError(54,0.001937327); + Allpt7TeVScaling->SetBinError(55,0.002108783); + Allpt7TeVScaling->SetBinError(56,0.001816997); + Allpt7TeVScaling->SetBinError(57,0.001892795); + Allpt7TeVScaling->SetBinError(58,0.001803616); + Allpt7TeVScaling->SetBinError(59,0.001994761); + Allpt7TeVScaling->SetBinError(60,0.001964909); + Allpt7TeVScaling->SetBinError(61,0.001255034); + Allpt7TeVScaling->SetBinError(62,0.001312232); + Allpt7TeVScaling->SetBinError(63,0.001358296); + Allpt7TeVScaling->SetBinError(64,0.001569736); + Allpt7TeVScaling->SetBinError(65,0.001469072); + Allpt7TeVScaling->SetBinError(66,0.001588803); + Allpt7TeVScaling->SetBinError(67,0.001683862); + Allpt7TeVScaling->SetBinError(68,0.00164699); + Allpt7TeVScaling->SetBinError(69,0.001760645); + Allpt7TeVScaling->SetBinError(70,0.00160423); + Allpt7TeVScaling->SetBinError(71,0.001693146); + Allpt7TeVScaling->SetBinError(72,0.002042705); + Allpt7TeVScaling->SetBinError(73,0.001965856); + Allpt7TeVScaling->SetBinError(74,0.002245841); + Allpt7TeVScaling->SetBinError(75,0.002520926); + Allpt7TeVScaling->SetBinError(76,0.002402524); + Allpt7TeVScaling->SetBinError(77,0.002627656); + Allpt7TeVScaling->SetBinError(78,0.002963027); + Allpt7TeVScaling->SetBinError(79,0.003088673); + Allpt7TeVScaling->SetBinError(80,0.002958055); + Allpt7TeVScaling->SetBinError(81,0.001434729); + Allpt7TeVScaling->SetBinError(82,0.001531537); + Allpt7TeVScaling->SetBinError(83,0.002200666); + Allpt7TeVScaling->SetBinError(84,0.00242244); + Allpt7TeVScaling->SetBinError(85,0.002270338); + Allpt7TeVScaling->SetBinError(86,0.002972849); + Allpt7TeVScaling->SetBinError(87,0.003966084); + Allpt7TeVScaling->SetBinError(88,0.002679263); + Allpt7TeVScaling->SetBinError(89,0.006083804); + Allpt7TeVScaling->SetBinError(90,0.004227019); + Allpt7TeVScaling->SetBinError(91,0.004886052); + Allpt7TeVScaling->SetBinError(93,0.006030514); + Allpt7TeVScaling->SetBinError(94,0.008892003); + Allpt7TeVScaling->SetBinError(95,0.01001297); + Allpt7TeVScaling->SetBinError(96,0.008480241); + Allpt7TeVScaling->SetBinError(97,0.007527505); + Allpt7TeVScaling->SetBinError(98,0.01889477); + Allpt7TeVScaling->SetBinError(103,0.065502); + Allpt7TeVScaling->SetMinimum(0); + Allpt7TeVScaling->SetMaximum(0.04); + Allpt7TeVScaling->SetEntries(1049.145); + Allpt7TeVScaling->SetStats(0); + Allpt7TeVScaling->SetMarkerStyle(20); + Allpt7TeVScaling->GetXaxis()->SetTitle("p_{T}"); + Allpt7TeVScaling->GetXaxis()->SetRange(1,91); + Allpt7TeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); + Allpt7TeVScaling->GetYaxis()->SetTitleOffset(1.2); + TH1D *Allpt900GeVScaling = new TH1D("Allpt900GeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); + Allpt900GeVScaling->SetBinContent(13,0.1950247); + Allpt900GeVScaling->SetBinContent(14,0.147561); + Allpt900GeVScaling->SetBinContent(15,0.136378); + Allpt900GeVScaling->SetBinContent(16,0.1233744); + Allpt900GeVScaling->SetBinContent(17,0.1126807); + Allpt900GeVScaling->SetBinContent(18,0.09909555); + Allpt900GeVScaling->SetBinContent(19,0.09191239); + Allpt900GeVScaling->SetBinContent(20,0.08438053); + Allpt900GeVScaling->SetBinContent(21,0.07119319); + Allpt900GeVScaling->SetBinContent(22,0.05773434); + Allpt900GeVScaling->SetBinContent(23,0.04292845); + Allpt900GeVScaling->SetBinContent(24,0.0396127); + Allpt900GeVScaling->SetBinContent(25,0.03395077); + Allpt900GeVScaling->SetBinContent(26,0.03305251); + Allpt900GeVScaling->SetBinContent(27,0.03297093); + Allpt900GeVScaling->SetBinContent(28,0.031215); + Allpt900GeVScaling->SetBinContent(29,0.03178372); + Allpt900GeVScaling->SetBinContent(30,0.03236985); + Allpt900GeVScaling->SetBinContent(31,0.03271276); + Allpt900GeVScaling->SetBinContent(32,0.02993124); + Allpt900GeVScaling->SetBinContent(33,0.03004746); + Allpt900GeVScaling->SetBinContent(34,0.03077143); + Allpt900GeVScaling->SetBinContent(35,0.03211214); + Allpt900GeVScaling->SetBinContent(36,0.02807512); + Allpt900GeVScaling->SetBinContent(37,0.02651954); + Allpt900GeVScaling->SetBinContent(38,0.02503885); + Allpt900GeVScaling->SetBinContent(39,0.02325777); + Allpt900GeVScaling->SetBinContent(40,0.02773712); + Allpt900GeVScaling->SetBinContent(41,0.02491155); + Allpt900GeVScaling->SetBinContent(42,0.02389028); + Allpt900GeVScaling->SetBinContent(43,0.02085532); + Allpt900GeVScaling->SetBinContent(44,0.01959721); + Allpt900GeVScaling->SetBinContent(45,0.01937255); + Allpt900GeVScaling->SetBinContent(46,0.01745646); + Allpt900GeVScaling->SetBinContent(47,0.01890853); + Allpt900GeVScaling->SetBinContent(48,0.02053692); + Allpt900GeVScaling->SetBinContent(49,0.01813791); + Allpt900GeVScaling->SetBinContent(50,0.02187495); + Allpt900GeVScaling->SetBinContent(51,0.02409687); + Allpt900GeVScaling->SetBinContent(52,0.0206565); + Allpt900GeVScaling->SetBinContent(53,0.02269165); + Allpt900GeVScaling->SetBinContent(54,0.01850795); + Allpt900GeVScaling->SetBinContent(55,0.02115427); + Allpt900GeVScaling->SetBinContent(56,0.01543302); + Allpt900GeVScaling->SetBinContent(57,0.01642095); + Allpt900GeVScaling->SetBinContent(58,0.01542331); + Allpt900GeVScaling->SetBinContent(59,0.01931254); + Allpt900GeVScaling->SetBinContent(60,0.01650654); + Allpt900GeVScaling->SetBinContent(61,0.01827512); + Allpt900GeVScaling->SetBinContent(62,0.01700125); + Allpt900GeVScaling->SetBinContent(63,0.01544321); + Allpt900GeVScaling->SetBinContent(64,0.01830588); + Allpt900GeVScaling->SetBinContent(65,0.01584369); + Allpt900GeVScaling->SetBinContent(66,0.0145421); + Allpt900GeVScaling->SetBinContent(67,0.01382206); + Allpt900GeVScaling->SetBinContent(68,0.01358825); + Allpt900GeVScaling->SetBinContent(69,0.01294888); + Allpt900GeVScaling->SetBinContent(70,0.009346655); + Allpt900GeVScaling->SetBinContent(71,0.01235781); + Allpt900GeVScaling->SetBinContent(72,0.0127376); + Allpt900GeVScaling->SetBinContent(73,0.008968166); + Allpt900GeVScaling->SetBinContent(74,0.01467765); + Allpt900GeVScaling->SetBinContent(75,0.01432309); + Allpt900GeVScaling->SetBinContent(76,0.01204856); + Allpt900GeVScaling->SetBinContent(77,0.01436071); + Allpt900GeVScaling->SetBinContent(78,0.01451222); + Allpt900GeVScaling->SetBinContent(79,0.01378967); + Allpt900GeVScaling->SetBinContent(80,0.01040148); + Allpt900GeVScaling->SetBinContent(81,0.01041213); + Allpt900GeVScaling->SetBinContent(82,0.009971816); + Allpt900GeVScaling->SetBinContent(83,0.01417668); + Allpt900GeVScaling->SetBinContent(84,0.008713292); + Allpt900GeVScaling->SetBinContent(85,0.005425386); + Allpt900GeVScaling->SetBinContent(86,0.005303683); + Allpt900GeVScaling->SetBinContent(87,0.006759432); + Allpt900GeVScaling->SetBinContent(88,0.004445477); + Allpt900GeVScaling->SetBinContent(89,0.01810702); + Allpt900GeVScaling->SetBinContent(90,0.0108909); + Allpt900GeVScaling->SetBinContent(91,0.01186746); + Allpt900GeVScaling->SetBinContent(92,0.006128581); + Allpt900GeVScaling->SetBinContent(93,0.007452721); + Allpt900GeVScaling->SetBinContent(94,0.0259471); + Allpt900GeVScaling->SetBinContent(95,0.0114321); + Allpt900GeVScaling->SetBinContent(98,0.03264295); + Allpt900GeVScaling->SetBinContent(103,0.07829265); + Allpt900GeVScaling->SetBinError(13,0.03504783); + Allpt900GeVScaling->SetBinError(14,0.007336667); + Allpt900GeVScaling->SetBinError(15,0.004765958); + Allpt900GeVScaling->SetBinError(16,0.004028143); + Allpt900GeVScaling->SetBinError(17,0.003760166); + Allpt900GeVScaling->SetBinError(18,0.003357858); + Allpt900GeVScaling->SetBinError(19,0.003156502); + Allpt900GeVScaling->SetBinError(20,0.002970822); + Allpt900GeVScaling->SetBinError(21,0.001886555); + Allpt900GeVScaling->SetBinError(22,0.001676045); + Allpt900GeVScaling->SetBinError(23,0.001411794); + Allpt900GeVScaling->SetBinError(24,0.001346425); + Allpt900GeVScaling->SetBinError(25,0.001254258); + Allpt900GeVScaling->SetBinError(26,0.001299384); + Allpt900GeVScaling->SetBinError(27,0.001328066); + Allpt900GeVScaling->SetBinError(28,0.001330088); + Allpt900GeVScaling->SetBinError(29,0.001446767); + Allpt900GeVScaling->SetBinError(30,0.00145402); + Allpt900GeVScaling->SetBinError(31,0.001463994); + Allpt900GeVScaling->SetBinError(32,0.001449209); + Allpt900GeVScaling->SetBinError(33,0.001414847); + Allpt900GeVScaling->SetBinError(34,0.001598429); + Allpt900GeVScaling->SetBinError(35,0.001590854); + Allpt900GeVScaling->SetBinError(36,0.001526798); + Allpt900GeVScaling->SetBinError(37,0.001486126); + Allpt900GeVScaling->SetBinError(38,0.001457142); + Allpt900GeVScaling->SetBinError(39,0.001474238); + Allpt900GeVScaling->SetBinError(40,0.001757902); + Allpt900GeVScaling->SetBinError(41,0.001712601); + Allpt900GeVScaling->SetBinError(42,0.001671385); + Allpt900GeVScaling->SetBinError(43,0.001570556); + Allpt900GeVScaling->SetBinError(44,0.001627725); + Allpt900GeVScaling->SetBinError(45,0.001637262); + Allpt900GeVScaling->SetBinError(46,0.001641227); + Allpt900GeVScaling->SetBinError(47,0.001716345); + Allpt900GeVScaling->SetBinError(48,0.001907778); + Allpt900GeVScaling->SetBinError(49,0.001842589); + Allpt900GeVScaling->SetBinError(50,0.00215445); + Allpt900GeVScaling->SetBinError(51,0.002310602); + Allpt900GeVScaling->SetBinError(52,0.002171046); + Allpt900GeVScaling->SetBinError(53,0.002386521); + Allpt900GeVScaling->SetBinError(54,0.002149237); + Allpt900GeVScaling->SetBinError(55,0.002233898); + Allpt900GeVScaling->SetBinError(56,0.001856697); + Allpt900GeVScaling->SetBinError(57,0.001891742); + Allpt900GeVScaling->SetBinError(58,0.001840334); + Allpt900GeVScaling->SetBinError(59,0.002168741); + Allpt900GeVScaling->SetBinError(60,0.002044863); + Allpt900GeVScaling->SetBinError(61,0.001431972); + Allpt900GeVScaling->SetBinError(62,0.001452706); + Allpt900GeVScaling->SetBinError(63,0.001477499); + Allpt900GeVScaling->SetBinError(64,0.001722091); + Allpt900GeVScaling->SetBinError(65,0.001674803); + Allpt900GeVScaling->SetBinError(66,0.001737491); + Allpt900GeVScaling->SetBinError(67,0.001783912); + Allpt900GeVScaling->SetBinError(68,0.00186443); + Allpt900GeVScaling->SetBinError(69,0.001965873); + Allpt900GeVScaling->SetBinError(70,0.001687347); + Allpt900GeVScaling->SetBinError(71,0.002072538); + Allpt900GeVScaling->SetBinError(72,0.002233273); + Allpt900GeVScaling->SetBinError(73,0.001965905); + Allpt900GeVScaling->SetBinError(74,0.002669571); + Allpt900GeVScaling->SetBinError(75,0.00278613); + Allpt900GeVScaling->SetBinError(76,0.002661412); + Allpt900GeVScaling->SetBinError(77,0.002952605); + Allpt900GeVScaling->SetBinError(78,0.003268721); + Allpt900GeVScaling->SetBinError(79,0.003273694); + Allpt900GeVScaling->SetBinError(80,0.002900898); + Allpt900GeVScaling->SetBinError(81,0.001677933); + Allpt900GeVScaling->SetBinError(82,0.001933408); + Allpt900GeVScaling->SetBinError(83,0.002759328); + Allpt900GeVScaling->SetBinError(84,0.002526604); + Allpt900GeVScaling->SetBinError(85,0.002297818); + Allpt900GeVScaling->SetBinError(86,0.002659213); + Allpt900GeVScaling->SetBinError(87,0.003392141); + Allpt900GeVScaling->SetBinError(88,0.003150511); + Allpt900GeVScaling->SetBinError(89,0.007499872); + Allpt900GeVScaling->SetBinError(90,0.006322856); + Allpt900GeVScaling->SetBinError(91,0.00689283); + Allpt900GeVScaling->SetBinError(92,0.006146806); + Allpt900GeVScaling->SetBinError(93,0.007481131); + Allpt900GeVScaling->SetBinError(94,0.01853139); + Allpt900GeVScaling->SetBinError(95,0.01149767); + Allpt900GeVScaling->SetBinError(98,0.02345991); + Allpt900GeVScaling->SetBinError(103,0.08126059); + Allpt900GeVScaling->SetMinimum(0); + Allpt900GeVScaling->SetMaximum(0.04); + Allpt900GeVScaling->SetEntries(724.3671); + Allpt900GeVScaling->SetStats(0); + Allpt900GeVScaling->SetMarkerStyle(20); + Allpt900GeVScaling->GetXaxis()->SetTitle("p_{T}"); + Allpt900GeVScaling->GetXaxis()->SetRange(1,91); + Allpt900GeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); + Allpt900GeVScaling->GetYaxis()->SetTitleOffset(1.2); + Allpt7TeVScaling->Add(Allpt900GeVScaling); + Allpt7TeVScaling->Scale(0.5); + delete Allpt900GeVScaling; + return Allpt7TeVScaling; +} + +TH1D *PbPb276TPCBkgd(){ + TH1D *Allpt7TeVScaling = new TH1D("Allpt7TeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); + Allpt7TeVScaling->SetBinContent(13,0.1745755); + Allpt7TeVScaling->SetBinContent(14,0.1624186); + Allpt7TeVScaling->SetBinContent(15,0.1431048); + Allpt7TeVScaling->SetBinContent(16,0.1331614); + Allpt7TeVScaling->SetBinContent(17,0.121262); + Allpt7TeVScaling->SetBinContent(18,0.1082923); + Allpt7TeVScaling->SetBinContent(19,0.1011451); + Allpt7TeVScaling->SetBinContent(20,0.08875923); + Allpt7TeVScaling->SetBinContent(21,0.08100933); + Allpt7TeVScaling->SetBinContent(22,0.06527292); + Allpt7TeVScaling->SetBinContent(23,0.05220065); + Allpt7TeVScaling->SetBinContent(24,0.04752904); + Allpt7TeVScaling->SetBinContent(25,0.04421568); + Allpt7TeVScaling->SetBinContent(26,0.03879551); + Allpt7TeVScaling->SetBinContent(27,0.03944679); + Allpt7TeVScaling->SetBinContent(28,0.03707215); + Allpt7TeVScaling->SetBinContent(29,0.03663599); + Allpt7TeVScaling->SetBinContent(30,0.03937437); + Allpt7TeVScaling->SetBinContent(31,0.03888919); + Allpt7TeVScaling->SetBinContent(32,0.03646502); + Allpt7TeVScaling->SetBinContent(33,0.03784992); + Allpt7TeVScaling->SetBinContent(34,0.03360079); + Allpt7TeVScaling->SetBinContent(35,0.03761566); + Allpt7TeVScaling->SetBinContent(36,0.03276653); + Allpt7TeVScaling->SetBinContent(37,0.03401566); + Allpt7TeVScaling->SetBinContent(38,0.02803893); + Allpt7TeVScaling->SetBinContent(39,0.0290018); + Allpt7TeVScaling->SetBinContent(40,0.03323201); + Allpt7TeVScaling->SetBinContent(41,0.03006077); + Allpt7TeVScaling->SetBinContent(42,0.02867674); + Allpt7TeVScaling->SetBinContent(43,0.02471421); + Allpt7TeVScaling->SetBinContent(44,0.02610375); + Allpt7TeVScaling->SetBinContent(45,0.02601621); + Allpt7TeVScaling->SetBinContent(46,0.02184781); + Allpt7TeVScaling->SetBinContent(47,0.02479646); + Allpt7TeVScaling->SetBinContent(48,0.02519956); + Allpt7TeVScaling->SetBinContent(49,0.02376095); + Allpt7TeVScaling->SetBinContent(50,0.02836189); + Allpt7TeVScaling->SetBinContent(51,0.0339595); + Allpt7TeVScaling->SetBinContent(52,0.03057126); + Allpt7TeVScaling->SetBinContent(53,0.03238621); + Allpt7TeVScaling->SetBinContent(54,0.0264272); + Allpt7TeVScaling->SetBinContent(55,0.03203146); + Allpt7TeVScaling->SetBinContent(56,0.02554187); + Allpt7TeVScaling->SetBinContent(57,0.02901947); + Allpt7TeVScaling->SetBinContent(58,0.02584016); + Allpt7TeVScaling->SetBinContent(59,0.03055042); + Allpt7TeVScaling->SetBinContent(60,0.02658193); + Allpt7TeVScaling->SetBinContent(61,0.02553988); + Allpt7TeVScaling->SetBinContent(62,0.02787508); + Allpt7TeVScaling->SetBinContent(63,0.02836785); + Allpt7TeVScaling->SetBinContent(64,0.03449766); + Allpt7TeVScaling->SetBinContent(65,0.02710397); + Allpt7TeVScaling->SetBinContent(66,0.02741966); + Allpt7TeVScaling->SetBinContent(67,0.03184318); + Allpt7TeVScaling->SetBinContent(68,0.02993716); + Allpt7TeVScaling->SetBinContent(69,0.03113033); + Allpt7TeVScaling->SetBinContent(70,0.01946293); + Allpt7TeVScaling->SetBinContent(71,0.02890383); + Allpt7TeVScaling->SetBinContent(72,0.02602078); + Allpt7TeVScaling->SetBinContent(73,0.02063028); + Allpt7TeVScaling->SetBinContent(74,0.03186679); + Allpt7TeVScaling->SetBinContent(75,0.0399062); + Allpt7TeVScaling->SetBinContent(76,0.02051914); + Allpt7TeVScaling->SetBinContent(77,0.04454688); + Allpt7TeVScaling->SetBinContent(78,0.04798449); + Allpt7TeVScaling->SetBinContent(79,0.0311417); + Allpt7TeVScaling->SetBinContent(80,0.03774521); + Allpt7TeVScaling->SetBinContent(81,0.0219128); + Allpt7TeVScaling->SetBinContent(82,0.02323492); + Allpt7TeVScaling->SetBinContent(83,0.02412182); + Allpt7TeVScaling->SetBinContent(84,0.03636555); + Allpt7TeVScaling->SetBinContent(85,0.007488934); + Allpt7TeVScaling->SetBinContent(86,0.01676086); + Allpt7TeVScaling->SetBinContent(87,0.02427021); + Allpt7TeVScaling->SetBinContent(88,0.02482882); + Allpt7TeVScaling->SetBinContent(89,0.03680926); + Allpt7TeVScaling->SetBinContent(90,0.01600685); + Allpt7TeVScaling->SetBinContent(91,0.006885685); + Allpt7TeVScaling->SetBinContent(93,0.0195411); + Allpt7TeVScaling->SetBinContent(94,0.008862645); + Allpt7TeVScaling->SetBinContent(95,0.009963223); + Allpt7TeVScaling->SetBinContent(96,0.01399091); + Allpt7TeVScaling->SetBinContent(97,0.007500135); + Allpt7TeVScaling->SetBinContent(98,0.01922651); + Allpt7TeVScaling->SetBinContent(103,0.03511665); + Allpt7TeVScaling->SetBinError(13,0.03042345); + Allpt7TeVScaling->SetBinError(14,0.007201474); + Allpt7TeVScaling->SetBinError(15,0.004623859); + Allpt7TeVScaling->SetBinError(16,0.004039713); + Allpt7TeVScaling->SetBinError(17,0.003694238); + Allpt7TeVScaling->SetBinError(18,0.003377785); + Allpt7TeVScaling->SetBinError(19,0.00319956); + Allpt7TeVScaling->SetBinError(20,0.00294949); + Allpt7TeVScaling->SetBinError(21,0.001965652); + Allpt7TeVScaling->SetBinError(22,0.001771706); + Allpt7TeVScaling->SetBinError(23,0.0015726); + Allpt7TeVScaling->SetBinError(24,0.001495611); + Allpt7TeVScaling->SetBinError(25,0.00145454); + Allpt7TeVScaling->SetBinError(26,0.001398941); + Allpt7TeVScaling->SetBinError(27,0.001495129); + Allpt7TeVScaling->SetBinError(28,0.001474393); + Allpt7TeVScaling->SetBinError(29,0.001608602); + Allpt7TeVScaling->SetBinError(30,0.001684714); + Allpt7TeVScaling->SetBinError(31,0.001755251); + Allpt7TeVScaling->SetBinError(32,0.001658064); + Allpt7TeVScaling->SetBinError(33,0.001734642); + Allpt7TeVScaling->SetBinError(34,0.001624169); + Allpt7TeVScaling->SetBinError(35,0.001893747); + Allpt7TeVScaling->SetBinError(36,0.001748243); + Allpt7TeVScaling->SetBinError(37,0.001824646); + Allpt7TeVScaling->SetBinError(38,0.001587514); + Allpt7TeVScaling->SetBinError(39,0.001697203); + Allpt7TeVScaling->SetBinError(40,0.002112545); + Allpt7TeVScaling->SetBinError(41,0.002060819); + Allpt7TeVScaling->SetBinError(42,0.002026571); + Allpt7TeVScaling->SetBinError(43,0.001821759); + Allpt7TeVScaling->SetBinError(44,0.002051793); + Allpt7TeVScaling->SetBinError(45,0.002193957); + Allpt7TeVScaling->SetBinError(46,0.001991362); + Allpt7TeVScaling->SetBinError(47,0.002242771); + Allpt7TeVScaling->SetBinError(48,0.002305595); + Allpt7TeVScaling->SetBinError(49,0.002335076); + Allpt7TeVScaling->SetBinError(50,0.002798601); + Allpt7TeVScaling->SetBinError(51,0.003536788); + Allpt7TeVScaling->SetBinError(52,0.003440952); + Allpt7TeVScaling->SetBinError(53,0.003495971); + Allpt7TeVScaling->SetBinError(54,0.003067046); + Allpt7TeVScaling->SetBinError(55,0.003553133); + Allpt7TeVScaling->SetBinError(56,0.003084681); + Allpt7TeVScaling->SetBinError(57,0.003207651); + Allpt7TeVScaling->SetBinError(58,0.002958898); + Allpt7TeVScaling->SetBinError(59,0.003486212); + Allpt7TeVScaling->SetBinError(60,0.003160784); + Allpt7TeVScaling->SetBinError(61,0.00221753); + Allpt7TeVScaling->SetBinError(62,0.002573063); + Allpt7TeVScaling->SetBinError(63,0.003014171); + Allpt7TeVScaling->SetBinError(64,0.003535776); + Allpt7TeVScaling->SetBinError(65,0.003472028); + Allpt7TeVScaling->SetBinError(66,0.003649103); + Allpt7TeVScaling->SetBinError(67,0.004547596); + Allpt7TeVScaling->SetBinError(68,0.004987629); + Allpt7TeVScaling->SetBinError(69,0.005726305); + Allpt7TeVScaling->SetBinError(70,0.003978767); + Allpt7TeVScaling->SetBinError(71,0.00668463); + Allpt7TeVScaling->SetBinError(72,0.005573138); + Allpt7TeVScaling->SetBinError(73,0.005310436); + Allpt7TeVScaling->SetBinError(74,0.007698246); + Allpt7TeVScaling->SetBinError(75,0.01002912); + Allpt7TeVScaling->SetBinError(76,0.006029285); + Allpt7TeVScaling->SetBinError(77,0.0117225); + Allpt7TeVScaling->SetBinError(78,0.01333887); + Allpt7TeVScaling->SetBinError(79,0.009949805); + Allpt7TeVScaling->SetBinError(80,0.01335105); + Allpt7TeVScaling->SetBinError(81,0.004966615); + Allpt7TeVScaling->SetBinError(82,0.006669427); + Allpt7TeVScaling->SetBinError(83,0.007591192); + Allpt7TeVScaling->SetBinError(84,0.01217267); + Allpt7TeVScaling->SetBinError(85,0.002708538); + Allpt7TeVScaling->SetBinError(86,0.008143581); + Allpt7TeVScaling->SetBinError(87,0.01149673); + Allpt7TeVScaling->SetBinError(88,0.01759547); + Allpt7TeVScaling->SetBinError(89,0.01766708); + Allpt7TeVScaling->SetBinError(90,0.01339557); + Allpt7TeVScaling->SetBinError(91,0.004886052); + Allpt7TeVScaling->SetBinError(93,0.01960311); + Allpt7TeVScaling->SetBinError(94,0.008892003); + Allpt7TeVScaling->SetBinError(95,0.01001297); + Allpt7TeVScaling->SetBinError(96,0.008178781); + Allpt7TeVScaling->SetBinError(97,0.007527505); + Allpt7TeVScaling->SetBinError(98,0.01463862); + Allpt7TeVScaling->SetBinError(103,0.03621621); + Allpt7TeVScaling->SetMinimum(0); + Allpt7TeVScaling->SetMaximum(0.04); + Allpt7TeVScaling->SetEntries(1949.423); + Allpt7TeVScaling->SetStats(0); + Allpt7TeVScaling->SetMarkerStyle(20); + Allpt7TeVScaling->GetXaxis()->SetTitle("p_{T}"); + Allpt7TeVScaling->GetXaxis()->SetRange(1,91); + Allpt7TeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); + Allpt7TeVScaling->GetYaxis()->SetTitleOffset(1.2); + + TH1D *Allpt900GeVScaling = new TH1D("Allpt900GeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); + Allpt900GeVScaling->SetBinContent(13,0.2098361); + Allpt900GeVScaling->SetBinContent(14,0.1557776); + Allpt900GeVScaling->SetBinContent(15,0.1480364); + Allpt900GeVScaling->SetBinContent(16,0.1368096); + Allpt900GeVScaling->SetBinContent(17,0.1244875); + Allpt900GeVScaling->SetBinContent(18,0.1119436); + Allpt900GeVScaling->SetBinContent(19,0.104192); + Allpt900GeVScaling->SetBinContent(20,0.0938359); + Allpt900GeVScaling->SetBinContent(21,0.08228509); + Allpt900GeVScaling->SetBinContent(22,0.06783117); + Allpt900GeVScaling->SetBinContent(23,0.05079209); + Allpt900GeVScaling->SetBinContent(24,0.04795673); + Allpt900GeVScaling->SetBinContent(25,0.04181242); + Allpt900GeVScaling->SetBinContent(26,0.03994469); + Allpt900GeVScaling->SetBinContent(27,0.03981588); + Allpt900GeVScaling->SetBinContent(28,0.0377358); + Allpt900GeVScaling->SetBinContent(29,0.0399667); + Allpt900GeVScaling->SetBinContent(30,0.04095473); + Allpt900GeVScaling->SetBinContent(31,0.04058456); + Allpt900GeVScaling->SetBinContent(32,0.03641218); + Allpt900GeVScaling->SetBinContent(33,0.03599966); + Allpt900GeVScaling->SetBinContent(34,0.03853838); + Allpt900GeVScaling->SetBinContent(35,0.04105409); + Allpt900GeVScaling->SetBinContent(36,0.03461075); + Allpt900GeVScaling->SetBinContent(37,0.03478729); + Allpt900GeVScaling->SetBinContent(38,0.03150168); + Allpt900GeVScaling->SetBinContent(39,0.03121329); + Allpt900GeVScaling->SetBinContent(40,0.0370806); + Allpt900GeVScaling->SetBinContent(41,0.03322445); + Allpt900GeVScaling->SetBinContent(42,0.03200198); + Allpt900GeVScaling->SetBinContent(43,0.02753292); + Allpt900GeVScaling->SetBinContent(44,0.02776155); + Allpt900GeVScaling->SetBinContent(45,0.02671865); + Allpt900GeVScaling->SetBinContent(46,0.02612601); + Allpt900GeVScaling->SetBinContent(47,0.02779079); + Allpt900GeVScaling->SetBinContent(48,0.02762603); + Allpt900GeVScaling->SetBinContent(49,0.0288977); + Allpt900GeVScaling->SetBinContent(50,0.03181282); + Allpt900GeVScaling->SetBinContent(51,0.04194707); + Allpt900GeVScaling->SetBinContent(52,0.03287613); + Allpt900GeVScaling->SetBinContent(53,0.04363778); + Allpt900GeVScaling->SetBinContent(54,0.03003709); + Allpt900GeVScaling->SetBinContent(55,0.03286254); + Allpt900GeVScaling->SetBinContent(56,0.02585759); + Allpt900GeVScaling->SetBinContent(57,0.03390592); + Allpt900GeVScaling->SetBinContent(58,0.03212414); + Allpt900GeVScaling->SetBinContent(59,0.03559123); + Allpt900GeVScaling->SetBinContent(60,0.02464371); + Allpt900GeVScaling->SetBinContent(61,0.03680997); + Allpt900GeVScaling->SetBinContent(62,0.03569024); + Allpt900GeVScaling->SetBinContent(63,0.03662885); + Allpt900GeVScaling->SetBinContent(64,0.04690057); + Allpt900GeVScaling->SetBinContent(65,0.04036927); + Allpt900GeVScaling->SetBinContent(66,0.03237555); + Allpt900GeVScaling->SetBinContent(67,0.04285642); + Allpt900GeVScaling->SetBinContent(68,0.04778185); + Allpt900GeVScaling->SetBinContent(69,0.0411135); + Allpt900GeVScaling->SetBinContent(70,0.02465524); + Allpt900GeVScaling->SetBinContent(71,0.03686565); + Allpt900GeVScaling->SetBinContent(72,0.03694553); + Allpt900GeVScaling->SetBinContent(73,0.02206743); + Allpt900GeVScaling->SetBinContent(74,0.04854334); + Allpt900GeVScaling->SetBinContent(75,0.04190287); + Allpt900GeVScaling->SetBinContent(76,0.03255752); + Allpt900GeVScaling->SetBinContent(77,0.0715847); + Allpt900GeVScaling->SetBinContent(78,0.07329647); + Allpt900GeVScaling->SetBinContent(79,0.02727783); + Allpt900GeVScaling->SetBinContent(80,0.041527); + Allpt900GeVScaling->SetBinContent(81,0.0400001); + Allpt900GeVScaling->SetBinContent(82,0.0246785); + Allpt900GeVScaling->SetBinContent(83,0.03546198); + Allpt900GeVScaling->SetBinContent(84,0.0244789); + Allpt900GeVScaling->SetBinContent(85,0.01692187); + Allpt900GeVScaling->SetBinContent(86,0.01669025); + Allpt900GeVScaling->SetBinContent(87,0.01506698); + Allpt900GeVScaling->SetBinContent(88,0.02896953); + Allpt900GeVScaling->SetBinContent(89,0.04019321); + Allpt900GeVScaling->SetBinContent(90,0.0108909); + Allpt900GeVScaling->SetBinContent(91,0.02059646); + Allpt900GeVScaling->SetBinContent(92,0.01781294); + Allpt900GeVScaling->SetBinContent(93,0.01930956); + Allpt900GeVScaling->SetBinContent(94,0.0259471); + Allpt900GeVScaling->SetBinContent(95,0.0114321); + Allpt900GeVScaling->SetBinContent(98,0.04961994); + Allpt900GeVScaling->SetBinContent(103,0.1436233); + Allpt900GeVScaling->SetBinError(13,0.04135227); + Allpt900GeVScaling->SetBinError(14,0.007955716); + Allpt900GeVScaling->SetBinError(15,0.005353158); + Allpt900GeVScaling->SetBinError(16,0.00465994); + Allpt900GeVScaling->SetBinError(17,0.004317818); + Allpt900GeVScaling->SetBinError(18,0.003992759); + Allpt900GeVScaling->SetBinError(19,0.003794229); + Allpt900GeVScaling->SetBinError(20,0.003428406); + Allpt900GeVScaling->SetBinError(21,0.002309662); + Allpt900GeVScaling->SetBinError(22,0.002082593); + Allpt900GeVScaling->SetBinError(23,0.001768867); + Allpt900GeVScaling->SetBinError(24,0.001734781); + Allpt900GeVScaling->SetBinError(25,0.001668325); + Allpt900GeVScaling->SetBinError(26,0.001736728); + Allpt900GeVScaling->SetBinError(27,0.001793642); + Allpt900GeVScaling->SetBinError(28,0.001824914); + Allpt900GeVScaling->SetBinError(29,0.002216558); + Allpt900GeVScaling->SetBinError(30,0.00224095); + Allpt900GeVScaling->SetBinError(31,0.002244602); + Allpt900GeVScaling->SetBinError(32,0.002073202); + Allpt900GeVScaling->SetBinError(33,0.001913981); + Allpt900GeVScaling->SetBinError(34,0.002456659); + Allpt900GeVScaling->SetBinError(35,0.002596584); + Allpt900GeVScaling->SetBinError(36,0.002232065); + Allpt900GeVScaling->SetBinError(37,0.00251725); + Allpt900GeVScaling->SetBinError(38,0.002196274); + Allpt900GeVScaling->SetBinError(39,0.002513529); + Allpt900GeVScaling->SetBinError(40,0.002973807); + Allpt900GeVScaling->SetBinError(41,0.002807704); + Allpt900GeVScaling->SetBinError(42,0.002773487); + Allpt900GeVScaling->SetBinError(43,0.002527062); + Allpt900GeVScaling->SetBinError(44,0.002921618); + Allpt900GeVScaling->SetBinError(45,0.002907984); + Allpt900GeVScaling->SetBinError(46,0.003115439); + Allpt900GeVScaling->SetBinError(47,0.003224686); + Allpt900GeVScaling->SetBinError(48,0.00312276); + Allpt900GeVScaling->SetBinError(49,0.00393121); + Allpt900GeVScaling->SetBinError(50,0.003948797); + Allpt900GeVScaling->SetBinError(51,0.005565144); + Allpt900GeVScaling->SetBinError(52,0.004615559); + Allpt900GeVScaling->SetBinError(53,0.006431509); + Allpt900GeVScaling->SetBinError(54,0.004787345); + Allpt900GeVScaling->SetBinError(55,0.004506561); + Allpt900GeVScaling->SetBinError(56,0.004277351); + Allpt900GeVScaling->SetBinError(57,0.005084143); + Allpt900GeVScaling->SetBinError(58,0.005085463); + Allpt900GeVScaling->SetBinError(59,0.005312616); + Allpt900GeVScaling->SetBinError(60,0.003899113); + Allpt900GeVScaling->SetBinError(61,0.003977416); + Allpt900GeVScaling->SetBinError(62,0.004292332); + Allpt900GeVScaling->SetBinError(63,0.005072534); + Allpt900GeVScaling->SetBinError(64,0.006271891); + Allpt900GeVScaling->SetBinError(65,0.006463988); + Allpt900GeVScaling->SetBinError(66,0.005869634); + Allpt900GeVScaling->SetBinError(67,0.008405251); + Allpt900GeVScaling->SetBinError(68,0.009847454); + Allpt900GeVScaling->SetBinError(69,0.009749354); + Allpt900GeVScaling->SetBinError(70,0.007179523); + Allpt900GeVScaling->SetBinError(71,0.01004867); + Allpt900GeVScaling->SetBinError(72,0.01065873); + Allpt900GeVScaling->SetBinError(73,0.007849093); + Allpt900GeVScaling->SetBinError(74,0.01402708); + Allpt900GeVScaling->SetBinError(75,0.01330162); + Allpt900GeVScaling->SetBinError(76,0.01245402); + Allpt900GeVScaling->SetBinError(77,0.02130282); + Allpt900GeVScaling->SetBinError(78,0.02296044); + Allpt900GeVScaling->SetBinError(79,0.01175655); + Allpt900GeVScaling->SetBinError(80,0.0187519); + Allpt900GeVScaling->SetBinError(81,0.009932965); + Allpt900GeVScaling->SetBinError(82,0.007798875); + Allpt900GeVScaling->SetBinError(83,0.01126926); + Allpt900GeVScaling->SetBinError(84,0.01123177); + Allpt900GeVScaling->SetBinError(85,0.01013651); + Allpt900GeVScaling->SetBinError(86,0.01187595); + Allpt900GeVScaling->SetBinError(87,0.009213867); + Allpt900GeVScaling->SetBinError(88,0.02054123); + Allpt900GeVScaling->SetBinError(89,0.01976358); + Allpt900GeVScaling->SetBinError(90,0.006322856); + Allpt900GeVScaling->SetBinError(91,0.01400211); + Allpt900GeVScaling->SetBinError(92,0.01786717); + Allpt900GeVScaling->SetBinError(93,0.01938449); + Allpt900GeVScaling->SetBinError(94,0.01853139); + Allpt900GeVScaling->SetBinError(95,0.01149767); + Allpt900GeVScaling->SetBinError(98,0.03786105); + Allpt900GeVScaling->SetBinError(103,0.1493976); + Allpt900GeVScaling->SetMinimum(0); + Allpt900GeVScaling->SetMaximum(0.04); + Allpt900GeVScaling->SetEntries(482.6447); + Allpt900GeVScaling->SetStats(0); + Allpt900GeVScaling->SetMarkerStyle(20); + Allpt900GeVScaling->GetXaxis()->SetTitle("p_{T}"); + Allpt900GeVScaling->GetXaxis()->SetRange(1,91); + Allpt900GeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); + Allpt900GeVScaling->GetYaxis()->SetTitleOffset(1.2); + Allpt7TeVScaling->Add(Allpt900GeVScaling); + Allpt7TeVScaling->Scale(0.5); + delete Allpt900GeVScaling; + return Allpt7TeVScaling; +} + +TH1D *pp276ITSBkgd(){ + TH1D *Allpt7TeVScaling = new TH1D("Allpt7TeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); + Allpt7TeVScaling->SetBinContent(11,0.2557976); + Allpt7TeVScaling->SetBinContent(12,0.1130502); + Allpt7TeVScaling->SetBinContent(13,0.07504393); + Allpt7TeVScaling->SetBinContent(14,0.05774832); + Allpt7TeVScaling->SetBinContent(15,0.06389806); + Allpt7TeVScaling->SetBinContent(16,0.0468956); + Allpt7TeVScaling->SetBinContent(17,0.04355167); + Allpt7TeVScaling->SetBinContent(18,0.04035013); + Allpt7TeVScaling->SetBinContent(19,0.04359472); + Allpt7TeVScaling->SetBinContent(20,0.03392261); + Allpt7TeVScaling->SetBinContent(21,0.0327676); + Allpt7TeVScaling->SetBinContent(22,0.02916985); + Allpt7TeVScaling->SetBinContent(23,0.02603138); + Allpt7TeVScaling->SetBinContent(24,0.02298032); + Allpt7TeVScaling->SetBinContent(25,0.02366939); + Allpt7TeVScaling->SetBinContent(26,0.01941947); + Allpt7TeVScaling->SetBinContent(27,0.02077003); + Allpt7TeVScaling->SetBinContent(28,0.01985847); + Allpt7TeVScaling->SetBinContent(29,0.0177459); + Allpt7TeVScaling->SetBinContent(30,0.01889314); + Allpt7TeVScaling->SetBinContent(31,0.01783637); + Allpt7TeVScaling->SetBinContent(32,0.01742894); + Allpt7TeVScaling->SetBinContent(33,0.01400672); + Allpt7TeVScaling->SetBinContent(34,0.01509014); + Allpt7TeVScaling->SetBinContent(35,0.01473409); + Allpt7TeVScaling->SetBinContent(36,0.01337399); + Allpt7TeVScaling->SetBinContent(37,0.01458302); + Allpt7TeVScaling->SetBinContent(38,0.01103552); + Allpt7TeVScaling->SetBinContent(39,0.01440217); + Allpt7TeVScaling->SetBinContent(40,0.0147644); + Allpt7TeVScaling->SetBinContent(41,0.01459072); + Allpt7TeVScaling->SetBinContent(42,0.01318416); + Allpt7TeVScaling->SetBinContent(43,0.01364365); + Allpt7TeVScaling->SetBinContent(44,0.01186023); + Allpt7TeVScaling->SetBinContent(45,0.01218222); + Allpt7TeVScaling->SetBinContent(46,0.01015807); + Allpt7TeVScaling->SetBinContent(47,0.01158454); + Allpt7TeVScaling->SetBinContent(48,0.009338749); + Allpt7TeVScaling->SetBinContent(49,0.01320016); + Allpt7TeVScaling->SetBinContent(50,0.01029357); + Allpt7TeVScaling->SetBinContent(51,0.01238829); + Allpt7TeVScaling->SetBinContent(52,0.008714745); + Allpt7TeVScaling->SetBinContent(53,0.01300408); + Allpt7TeVScaling->SetBinContent(54,0.009445714); + Allpt7TeVScaling->SetBinContent(55,0.01301231); + Allpt7TeVScaling->SetBinContent(56,0.01010809); + Allpt7TeVScaling->SetBinContent(57,0.01181843); + Allpt7TeVScaling->SetBinContent(58,0.01037317); + Allpt7TeVScaling->SetBinContent(59,0.01172911); + Allpt7TeVScaling->SetBinContent(60,0.007576451); + Allpt7TeVScaling->SetBinContent(61,0.011185); + Allpt7TeVScaling->SetBinContent(62,0.01009882); + Allpt7TeVScaling->SetBinContent(63,0.01042078); + Allpt7TeVScaling->SetBinContent(64,0.01079128); + Allpt7TeVScaling->SetBinContent(65,0.008657009); + Allpt7TeVScaling->SetBinContent(66,0.009641373); + Allpt7TeVScaling->SetBinContent(67,0.01201938); + Allpt7TeVScaling->SetBinContent(68,0.01031088); + Allpt7TeVScaling->SetBinContent(69,0.009905995); + Allpt7TeVScaling->SetBinContent(70,0.009569216); + Allpt7TeVScaling->SetBinContent(71,0.009473246); + Allpt7TeVScaling->SetBinContent(72,0.00825013); + Allpt7TeVScaling->SetBinContent(73,0.008405926); + Allpt7TeVScaling->SetBinContent(74,0.006479467); + Allpt7TeVScaling->SetBinContent(75,0.01135131); + Allpt7TeVScaling->SetBinContent(76,0.007318432); + Allpt7TeVScaling->SetBinContent(77,0.01145386); + Allpt7TeVScaling->SetBinContent(78,0.01425465); + Allpt7TeVScaling->SetBinContent(79,0.01205772); + Allpt7TeVScaling->SetBinContent(80,0.009454997); + Allpt7TeVScaling->SetBinContent(81,0.006863416); + Allpt7TeVScaling->SetBinContent(82,0.006285407); + Allpt7TeVScaling->SetBinContent(83,0.004974776); + Allpt7TeVScaling->SetBinContent(84,0.00748478); + Allpt7TeVScaling->SetBinContent(85,0.002485323); + Allpt7TeVScaling->SetBinContent(86,0.009973802); + Allpt7TeVScaling->SetBinContent(87,0.00543998); + Allpt7TeVScaling->SetBinContent(89,0.01703395); + Allpt7TeVScaling->SetBinContent(90,0.004394557); + Allpt7TeVScaling->SetBinContent(96,0.01833441); + Allpt7TeVScaling->SetBinContent(98,0.02964621); + Allpt7TeVScaling->SetBinContent(101,0.0574538); + Allpt7TeVScaling->SetBinError(11,0.07407797); + Allpt7TeVScaling->SetBinError(12,0.01275005); + Allpt7TeVScaling->SetBinError(13,0.00608927); + Allpt7TeVScaling->SetBinError(14,0.00428747); + Allpt7TeVScaling->SetBinError(15,0.004071853); + Allpt7TeVScaling->SetBinError(16,0.003225169); + Allpt7TeVScaling->SetBinError(17,0.002968269); + Allpt7TeVScaling->SetBinError(18,0.002761875); + Allpt7TeVScaling->SetBinError(19,0.002800096); + Allpt7TeVScaling->SetBinError(20,0.002386005); + Allpt7TeVScaling->SetBinError(21,0.001589813); + Allpt7TeVScaling->SetBinError(22,0.001475329); + Allpt7TeVScaling->SetBinError(23,0.001364666); + Allpt7TeVScaling->SetBinError(24,0.00126007); + Allpt7TeVScaling->SetBinError(25,0.001265789); + Allpt7TeVScaling->SetBinError(26,0.001156208); + Allpt7TeVScaling->SetBinError(27,0.001182566); + Allpt7TeVScaling->SetBinError(28,0.001171098); + Allpt7TeVScaling->SetBinError(29,0.001122203); + Allpt7TeVScaling->SetBinError(30,0.001173633); + Allpt7TeVScaling->SetBinError(31,0.001154546); + Allpt7TeVScaling->SetBinError(32,0.001162445); + Allpt7TeVScaling->SetBinError(33,0.001067574); + Allpt7TeVScaling->SetBinError(34,0.001116243); + Allpt7TeVScaling->SetBinError(35,0.001144851); + Allpt7TeVScaling->SetBinError(36,0.001110251); + Allpt7TeVScaling->SetBinError(37,0.001192711); + Allpt7TeVScaling->SetBinError(38,0.001062774); + Allpt7TeVScaling->SetBinError(39,0.001255752); + Allpt7TeVScaling->SetBinError(40,0.001319021); + Allpt7TeVScaling->SetBinError(41,0.001323284); + Allpt7TeVScaling->SetBinError(42,0.001303463); + Allpt7TeVScaling->SetBinError(43,0.001388458); + Allpt7TeVScaling->SetBinError(44,0.001312344); + Allpt7TeVScaling->SetBinError(45,0.001383885); + Allpt7TeVScaling->SetBinError(46,0.001289152); + Allpt7TeVScaling->SetBinError(47,0.001430287); + Allpt7TeVScaling->SetBinError(48,0.001302716); + Allpt7TeVScaling->SetBinError(49,0.001578988); + Allpt7TeVScaling->SetBinError(50,0.0014581); + Allpt7TeVScaling->SetBinError(51,0.00165325); + Allpt7TeVScaling->SetBinError(52,0.001404793); + Allpt7TeVScaling->SetBinError(53,0.001757303); + Allpt7TeVScaling->SetBinError(54,0.001562987); + Allpt7TeVScaling->SetBinError(55,0.001911932); + Allpt7TeVScaling->SetBinError(56,0.001695077); + Allpt7TeVScaling->SetBinError(57,0.001932989); + Allpt7TeVScaling->SetBinError(58,0.001844963); + Allpt7TeVScaling->SetBinError(59,0.002000117); + Allpt7TeVScaling->SetBinError(60,0.001689391); + Allpt7TeVScaling->SetBinError(61,0.001386339); + Allpt7TeVScaling->SetBinError(62,0.001349763); + Allpt7TeVScaling->SetBinError(63,0.001454327); + Allpt7TeVScaling->SetBinError(64,0.001584791); + Allpt7TeVScaling->SetBinError(65,0.00151686); + Allpt7TeVScaling->SetBinError(66,0.001668507); + Allpt7TeVScaling->SetBinError(67,0.001963962); + Allpt7TeVScaling->SetBinError(68,0.001970941); + Allpt7TeVScaling->SetBinError(69,0.001992913); + Allpt7TeVScaling->SetBinError(70,0.002006771); + Allpt7TeVScaling->SetBinError(71,0.002187254); + Allpt7TeVScaling->SetBinError(72,0.002140473); + Allpt7TeVScaling->SetBinError(73,0.002256925); + Allpt7TeVScaling->SetBinError(74,0.002056219); + Allpt7TeVScaling->SetBinError(75,0.002854892); + Allpt7TeVScaling->SetBinError(76,0.002450159); + Allpt7TeVScaling->SetBinError(77,0.003197387); + Allpt7TeVScaling->SetBinError(78,0.003840947); + Allpt7TeVScaling->SetBinError(79,0.003660464); + Allpt7TeVScaling->SetBinError(80,0.003363956); + Allpt7TeVScaling->SetBinError(81,0.001624793); + Allpt7TeVScaling->SetBinError(82,0.00181065); + Allpt7TeVScaling->SetBinError(83,0.001888966); + Allpt7TeVScaling->SetBinError(84,0.002671173); + Allpt7TeVScaling->SetBinError(85,0.001759843); + Allpt7TeVScaling->SetBinError(86,0.004094911); + Allpt7TeVScaling->SetBinError(87,0.003871198); + Allpt7TeVScaling->SetBinError(89,0.008630108); + Allpt7TeVScaling->SetBinError(90,0.004405367); + Allpt7TeVScaling->SetBinError(96,0.01308755); + Allpt7TeVScaling->SetBinError(98,0.02127892); + Allpt7TeVScaling->SetBinError(101,0.05904744); + Allpt7TeVScaling->SetMinimum(0); + Allpt7TeVScaling->SetMaximum(0.04); + Allpt7TeVScaling->SetEntries(300.6659); + Allpt7TeVScaling->SetStats(0); + Allpt7TeVScaling->SetMarkerStyle(20); + Allpt7TeVScaling->GetXaxis()->SetTitle("p_{T}"); + Allpt7TeVScaling->GetXaxis()->SetRange(1,61); + Allpt7TeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); + Allpt7TeVScaling->GetYaxis()->SetTitleOffset(1.2); + TH1D *Allpt900GeVScaling = new TH1D("Allpt900GeVScaling","Reconstructed E_{T} from misidentified electrons",111, xAxis1); + Allpt900GeVScaling->SetBinContent(9,1); + Allpt900GeVScaling->SetBinContent(11,0.2834258); + Allpt900GeVScaling->SetBinContent(12,0.1188776); + Allpt900GeVScaling->SetBinContent(13,0.08085647); + Allpt900GeVScaling->SetBinContent(14,0.05677858); + Allpt900GeVScaling->SetBinContent(15,0.0616757); + Allpt900GeVScaling->SetBinContent(16,0.04844882); + Allpt900GeVScaling->SetBinContent(17,0.03999142); + Allpt900GeVScaling->SetBinContent(18,0.03890007); + Allpt900GeVScaling->SetBinContent(19,0.0406008); + Allpt900GeVScaling->SetBinContent(20,0.02900627); + Allpt900GeVScaling->SetBinContent(21,0.03088961); + Allpt900GeVScaling->SetBinContent(22,0.03033992); + Allpt900GeVScaling->SetBinContent(23,0.02521969); + Allpt900GeVScaling->SetBinContent(24,0.02143008); + Allpt900GeVScaling->SetBinContent(25,0.02266187); + Allpt900GeVScaling->SetBinContent(26,0.01981764); + Allpt900GeVScaling->SetBinContent(27,0.02127366); + Allpt900GeVScaling->SetBinContent(28,0.0192753); + Allpt900GeVScaling->SetBinContent(29,0.01854744); + Allpt900GeVScaling->SetBinContent(30,0.01903471); + Allpt900GeVScaling->SetBinContent(31,0.01819287); + Allpt900GeVScaling->SetBinContent(32,0.01715784); + Allpt900GeVScaling->SetBinContent(33,0.01488607); + Allpt900GeVScaling->SetBinContent(34,0.0156348); + Allpt900GeVScaling->SetBinContent(35,0.01699235); + Allpt900GeVScaling->SetBinContent(36,0.01398954); + Allpt900GeVScaling->SetBinContent(37,0.01435865); + Allpt900GeVScaling->SetBinContent(38,0.01236315); + Allpt900GeVScaling->SetBinContent(39,0.01282288); + Allpt900GeVScaling->SetBinContent(40,0.01668125); + Allpt900GeVScaling->SetBinContent(41,0.01372736); + Allpt900GeVScaling->SetBinContent(42,0.01433745); + Allpt900GeVScaling->SetBinContent(43,0.01477496); + Allpt900GeVScaling->SetBinContent(44,0.0116776); + Allpt900GeVScaling->SetBinContent(45,0.01245205); + Allpt900GeVScaling->SetBinContent(46,0.009475905); + Allpt900GeVScaling->SetBinContent(47,0.01117031); + Allpt900GeVScaling->SetBinContent(48,0.009558273); + Allpt900GeVScaling->SetBinContent(49,0.01404047); + Allpt900GeVScaling->SetBinContent(50,0.01037464); + Allpt900GeVScaling->SetBinContent(51,0.01195013); + Allpt900GeVScaling->SetBinContent(52,0.01155682); + Allpt900GeVScaling->SetBinContent(53,0.0137127); + Allpt900GeVScaling->SetBinContent(54,0.009255665); + Allpt900GeVScaling->SetBinContent(55,0.01133598); + Allpt900GeVScaling->SetBinContent(56,0.009080946); + Allpt900GeVScaling->SetBinContent(57,0.01124535); + Allpt900GeVScaling->SetBinContent(58,0.008398657); + Allpt900GeVScaling->SetBinContent(59,0.008324869); + Allpt900GeVScaling->SetBinContent(60,0.007103489); + Allpt900GeVScaling->SetBinContent(61,0.01142644); + Allpt900GeVScaling->SetBinContent(62,0.0110906); + Allpt900GeVScaling->SetBinContent(63,0.01076479); + Allpt900GeVScaling->SetBinContent(64,0.01168746); + Allpt900GeVScaling->SetBinContent(65,0.008513672); + Allpt900GeVScaling->SetBinContent(66,0.00828279); + Allpt900GeVScaling->SetBinContent(67,0.01140648); + Allpt900GeVScaling->SetBinContent(68,0.01111323); + Allpt900GeVScaling->SetBinContent(69,0.007425678); + Allpt900GeVScaling->SetBinContent(70,0.009847098); + Allpt900GeVScaling->SetBinContent(71,0.009834419); + Allpt900GeVScaling->SetBinContent(72,0.008930991); + Allpt900GeVScaling->SetBinContent(73,0.00669819); + Allpt900GeVScaling->SetBinContent(74,0.00662618); + Allpt900GeVScaling->SetBinContent(75,0.01370167); + Allpt900GeVScaling->SetBinContent(76,0.006451024); + Allpt900GeVScaling->SetBinContent(77,0.01307605); + Allpt900GeVScaling->SetBinContent(78,0.01579264); + Allpt900GeVScaling->SetBinContent(79,0.01255205); + Allpt900GeVScaling->SetBinContent(80,0.01010045); + Allpt900GeVScaling->SetBinContent(81,0.01082155); + Allpt900GeVScaling->SetBinContent(82,0.006495255); + Allpt900GeVScaling->SetBinContent(83,0.005512525); + Allpt900GeVScaling->SetBinContent(84,0.003527945); + Allpt900GeVScaling->SetBinContent(85,0.001806946); + Allpt900GeVScaling->SetBinContent(86,0.004225982); + Allpt900GeVScaling->SetBinContent(89,0.01106803); + Allpt900GeVScaling->SetBinContent(90,0.005044269); + Allpt900GeVScaling->SetBinContent(98,0.04456838); + Allpt900GeVScaling->SetBinContent(101,0.07960857); + Allpt900GeVScaling->SetBinError(9,1.414214); + Allpt900GeVScaling->SetBinError(11,0.08936516); + Allpt900GeVScaling->SetBinError(12,0.01474769); + Allpt900GeVScaling->SetBinError(13,0.006947232); + Allpt900GeVScaling->SetBinError(14,0.00470276); + Allpt900GeVScaling->SetBinError(15,0.004394222); + Allpt900GeVScaling->SetBinError(16,0.003598075); + Allpt900GeVScaling->SetBinError(17,0.00312581); + Allpt900GeVScaling->SetBinError(18,0.003002622); + Allpt900GeVScaling->SetBinError(19,0.002976063); + Allpt900GeVScaling->SetBinError(20,0.002438669); + Allpt900GeVScaling->SetBinError(21,0.00172144); + Allpt900GeVScaling->SetBinError(22,0.001661658); + Allpt900GeVScaling->SetBinError(23,0.001477901); + Allpt900GeVScaling->SetBinError(24,0.001329128); + Allpt900GeVScaling->SetBinError(25,0.00137131); + Allpt900GeVScaling->SetBinError(26,0.001299862); + Allpt900GeVScaling->SetBinError(27,0.001332285); + Allpt900GeVScaling->SetBinError(28,0.001274491); + Allpt900GeVScaling->SetBinError(29,0.001266741); + Allpt900GeVScaling->SetBinError(30,0.001305369); + Allpt900GeVScaling->SetBinError(31,0.001295113); + Allpt900GeVScaling->SetBinError(32,0.001278961); + Allpt900GeVScaling->SetBinError(33,0.00122243); + Allpt900GeVScaling->SetBinError(34,0.001256858); + Allpt900GeVScaling->SetBinError(35,0.001357681); + Allpt900GeVScaling->SetBinError(36,0.001267233); + Allpt900GeVScaling->SetBinError(37,0.00131839); + Allpt900GeVScaling->SetBinError(38,0.001256835); + Allpt900GeVScaling->SetBinError(39,0.001306764); + Allpt900GeVScaling->SetBinError(40,0.001544152); + Allpt900GeVScaling->SetBinError(41,0.001438378); + Allpt900GeVScaling->SetBinError(42,0.001507702); + Allpt900GeVScaling->SetBinError(43,0.001615741); + Allpt900GeVScaling->SetBinError(44,0.00143817); + Allpt900GeVScaling->SetBinError(45,0.001536765); + Allpt900GeVScaling->SetBinError(46,0.001377493); + Allpt900GeVScaling->SetBinError(47,0.001536932); + Allpt900GeVScaling->SetBinError(48,0.001466076); + Allpt900GeVScaling->SetBinError(49,0.001814185); + Allpt900GeVScaling->SetBinError(50,0.001620601); + Allpt900GeVScaling->SetBinError(51,0.001814658); + Allpt900GeVScaling->SetBinError(52,0.001822325); + Allpt900GeVScaling->SetBinError(53,0.001976947); + Allpt900GeVScaling->SetBinError(54,0.001729338); + Allpt900GeVScaling->SetBinError(55,0.001988634); + Allpt900GeVScaling->SetBinError(56,0.001792157); + Allpt900GeVScaling->SetBinError(57,0.00206644); + Allpt900GeVScaling->SetBinError(58,0.001842037); + Allpt900GeVScaling->SetBinError(59,0.001826512); + Allpt900GeVScaling->SetBinError(60,0.001732292); + Allpt900GeVScaling->SetBinError(61,0.001512696); + Allpt900GeVScaling->SetBinError(62,0.001543542); + Allpt900GeVScaling->SetBinError(63,0.001633358); + Allpt900GeVScaling->SetBinError(64,0.001816012); + Allpt900GeVScaling->SetBinError(65,0.001680175); + Allpt900GeVScaling->SetBinError(66,0.001692314); + Allpt900GeVScaling->SetBinError(67,0.00213108); + Allpt900GeVScaling->SetBinError(68,0.002295308); + Allpt900GeVScaling->SetBinError(69,0.001926907); + Allpt900GeVScaling->SetBinError(70,0.002272706); + Allpt900GeVScaling->SetBinError(71,0.002475109); + Allpt900GeVScaling->SetBinError(72,0.002399123); + Allpt900GeVScaling->SetBinError(73,0.002241209); + Allpt900GeVScaling->SetBinError(74,0.002365641); + Allpt900GeVScaling->SetBinError(75,0.003450507); + Allpt900GeVScaling->SetBinError(76,0.002503392); + Allpt900GeVScaling->SetBinError(77,0.003801781); + Allpt900GeVScaling->SetBinError(78,0.004417338); + Allpt900GeVScaling->SetBinError(79,0.00421245); + Allpt900GeVScaling->SetBinError(80,0.003844397); + Allpt900GeVScaling->SetBinError(81,0.002274688); + Allpt900GeVScaling->SetBinError(82,0.002043975); + Allpt900GeVScaling->SetBinError(83,0.002262563); + Allpt900GeVScaling->SetBinError(84,0.00204239); + Allpt900GeVScaling->SetBinError(85,0.00180851); + Allpt900GeVScaling->SetBinError(86,0.002994568); + Allpt900GeVScaling->SetBinError(89,0.007939965); + Allpt900GeVScaling->SetBinError(90,0.005058548); + Allpt900GeVScaling->SetBinError(98,0.03222111); + Allpt900GeVScaling->SetBinError(101,0.08264771); + Allpt900GeVScaling->SetMinimum(0); + Allpt900GeVScaling->SetMaximum(0.04); + Allpt900GeVScaling->SetEntries(3.837105); + Allpt900GeVScaling->SetStats(0); + Allpt900GeVScaling->SetMarkerStyle(20); + Allpt900GeVScaling->GetXaxis()->SetTitle("p_{T}"); + Allpt900GeVScaling->GetXaxis()->SetRange(1,61); + Allpt900GeVScaling->GetYaxis()->SetTitle("Ratio of E_{T}^{background}/E_{T}^{had, meas.}"); + Allpt900GeVScaling->GetYaxis()->SetTitleOffset(1.2); + + Allpt7TeVScaling->Add(Allpt900GeVScaling); + Allpt7TeVScaling->Scale(0.5); + delete Allpt900GeVScaling; + return Allpt7TeVScaling; +}