singlee_nicePlots.cpp File Reference

#include <vector.h>
#include <map>

Include dependency graph for singlee_nicePlots.cpp:

Go to the source code of this file.

Functions
void	loadPresentationStyle ()
std::vector< Float_t >	rgaufit (TH1F *h, int iter)
void	singlee_nicePlots ()
float	effRMS (TH1F *h)

---

Function Documentation

float effRMS (  TH1F *  h  )

Definition at line 474 of file singlee_nicePlots.cpp. References effRMS(). Referenced by effRMS(), and singlee_nicePlots(). { int nbin = h->GetNbinsX(); int nevent = h->GetEntries(); int integral = h->Integral(); int maxBin = h->GetMaximumBin(); int leftBin = maxBin, rightBin = maxBin; --leftBin; ++rightBin; float sum = h->GetBinContent(maxBin); float add = 0; int last = 0; // cout << nbin << " | " << maxBin << " | " << sum << " | " << integral << " | " << h->GetBinContent(leftBin) << " | " << // h->GetBinContent(rightBin) << endl; while (sum < 0.6827*integral) { if (h->GetBinContent(leftBin) > h->GetBinContent(rightBin)) { add = h->GetBinContent(leftBin); --leftBin; last = 1; } else { add = h->GetBinContent(rightBin); ++rightBin; last = 2; } sum = sum + add; } // if (integral != 0) cout << (float) sum/integral*100. << " | "; if (last = 0) return -999.; if (last == 1) ++leftBin; if (last == 2) --rightBin; float effRMS = (h->GetBinCenter(rightBin) - h->GetBinCenter(leftBin))/2.; return effRMS; }

void loadPresentationStyle (   )

Definition at line 8 of file singlee_nicePlots.cpp. Referenced by singlee_nicePlots(). { gROOT->SetStyle("Plain"); gStyle->SetOptTitle(0); gStyle->SetPadTickX(1); gStyle->SetPadTickY(1); }

std::vector<Float_t> rgaufit (  TH1F *  h, int  iter )

Definition at line 16 of file singlee_nicePlots.cpp. { h->Fit("gaus","q"); float m; float s; TF1 *f; for (int i=0; i<iter; i++) { f = (TF1*)h->GetFunction("gaus"); m = f->GetParameter(1); s = f->GetParameter(2); h->Fit("gaus","qr","",m-1.15*s,m+3.*s); //printf("%d %f\n",i,h->GetFunction("gaus")->GetParameter(2)); } //printf("%f\n",h->GetFunction("gaus")->GetParameter(2)); f = (TF1*)h->GetFunction("gaus"); std::vector<Float_t> par; for (int i=0; i<f->GetNumberFreeParameters(); i++) par.push_back(f->GetParameter(i)); for (int i=0; i<f->GetNumberFreeParameters(); i++) par.push_back(f->GetParError(i)); return par; }

void singlee_nicePlots (   )

Definition at line 37 of file singlee_nicePlots.cpp. References effRMS(), and loadPresentationStyle(). { loadPresentationStyle(); // output figures: file type suffix char suffix[] = "eps"; // output figures: directory char outDir[] = "Figs"; // temp variables char str[128]; TFile *finResolutions = new TFile("histosFrame.root"); TNtuple *ntuple = (TNtuple *) finResolutions->Get("ntu"); vector<TH1F *> histosE, histosP, histosComb; //float start = 5., stop = 50., step = 3; float start = 5., stop = 100., step = 5; float upper = (int)(stop-start)/step; char stmp[10]; for (int ih = 0; ih < upper; ih++){ sprintf(stmp,"hE%d",ih); histosE.push_back(new TH1F(stmp, stmp, 2000, -2., 2.)); sprintf(stmp,"hP%d",ih); histosP.push_back(new TH1F(stmp, stmp, 2000, -2., 2.)); sprintf(stmp,"hComb%d",ih); histosComb.push_back(new TH1F(stmp, stmp, 2000, -2., 2.)); } float iev ,ie, ne, ETrue, Euncorr, E, p, comb, etaTrue; ntuple->SetBranchAddress("iev",&iev); ntuple->SetBranchAddress("ne",&ne); ntuple->SetBranchAddress("ie",&ie); ntuple->SetBranchAddress("ETrue",&ETrue); ntuple->SetBranchAddress("Euncorr",&Euncorr); ntuple->SetBranchAddress("E",&E); ntuple->SetBranchAddress("p",&p); ntuple->SetBranchAddress("comb",&comb); ntuple->SetBranchAddress("etaTrue",&etaTrue); for (int i = 0; i < ntuple->GetEntries(); i++) { ntuple->GetEntry(i); if (fabs(etaTrue) > 1.460 || ETrue > stop) continue; int idh = (int)(ETrue-start)/step; Axis_t eRat = (E - ETrue)/ETrue; Axis_t pRat = (p - ETrue)/ETrue; Axis_t combRat = (comb-ETrue)/ETrue; histosE[idh]->Fill(eRat); histosP[idh]->Fill(pRat); histosComb[idh]->Fill(combRat); } // TCanvas *proba = new TCanvas("proba","proba"); // proba->Divide(5,5); TGraph *gEffRMSE = new TGraph(); TGraph *gEffRMSP = new TGraph(); TGraph *gEffRMSComb = new TGraph(); for (int ih = 0; ih < upper; ih++){ // proba->cd(ih); // histosComb[0]->Draw(); float effRMSE = effRMS(histosE[ih]); float effRMSP = effRMS(histosP[ih]); float effRMSComb = effRMS(histosComb[ih]); cout << ih << " | " << effRMSE << " | " << effRMSP << " | " << effRMSComb << endl; gEffRMSE->SetPoint(ih,start+ih*step+step/2.,effRMSE); gEffRMSP->SetPoint(ih,start+ih*step+step/2.,effRMSP); gEffRMSComb->SetPoint(ih,start+ih*step+step/2.,effRMSComb); } TMultiGraph *mg = new TMultiGraph(); TCanvas *cRes = new TCanvas("cEffRMS","cEffRMS"); cRes->cd(); //gEffRMSE->Draw("AP"); mg->Add(gEffRMSE); gEffRMSE->SetMarkerStyle(20); gEffRMSE->SetMarkerColor(2); gEffRMSE->SetMarkerSize(.75); //gEffRMSP->Draw("P"); mg->Add(gEffRMSP); gEffRMSP->SetMarkerStyle(20); gEffRMSP->SetMarkerColor(3); gEffRMSP->SetMarkerSize(.75); //gEffRMSComb->Draw("P"); mg->Add(gEffRMSComb); gEffRMSComb->SetMarkerStyle(20); gEffRMSComb->SetMarkerColor(4); gEffRMSComb->SetMarkerSize(.75); mg->Draw("AP"); mg->GetXaxis()->SetTitle("E [GeV]"); mg->GetYaxis()->SetTitle("eff RMS"); gPad->Modified(); gPad->Update(); sprintf(str,"%s/effRMS.%s",outDir,suffix); gPad->Print(str); // HLT electrons and classes TFile *finElectronHLTClass = new TFile("histosFrame.root"); // eta distribution for electron classes - all electrons TH1F *hElectronHLTEta = (TH1F*)finElectronHLTClass->Get("hElectronHLTEta"); TH1F *hElectronEta = (TH1F*)finElectronHLTClass->Get("hElectronEta"); TH1F *hElectronHLTEtaClass[9]; hElectronHLTEtaClass[0] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass0"); hElectronHLTEtaClass[1] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass1"); hElectronHLTEtaClass[2] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass2"); hElectronHLTEtaClass[3] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass3"); hElectronHLTEtaClass[4] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass4"); hElectronHLTEtaClass[5] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass10"); hElectronHLTEtaClass[6] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass11"); hElectronHLTEtaClass[7] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass12"); hElectronHLTEtaClass[8] = (TH1F*)finElectronHLTClass->Get("hElectronHLTEtaClass13"); TH1F *hElectronEtaClass[9]; hElectronEtaClass[0] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass0"); hElectronEtaClass[1] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass1"); hElectronEtaClass[2] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass2"); hElectronEtaClass[3] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass3"); hElectronEtaClass[4] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass4"); hElectronEtaClass[5] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass10"); hElectronEtaClass[6] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass11"); hElectronEtaClass[7] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass12"); hElectronEtaClass[8] = (TH1F*)finElectronHLTClass->Get("hElectronEtaClass13"); cout << hElectronEtaClass[0]->GetEntries() << endl; char className[5][64]; sprintf(className[0],"golden"); sprintf(className[1],"big brem"); sprintf(className[2],"narrow"); sprintf(className[3],"showering"); sprintf(className[4],"cracks"); int hColors[5] = { 3, 4, 6, 7, 2 }; TLegend *legClass = new TLegend(.125,.62,.35,.885,"","NDC"); legClass->SetFillColor(0); legClass->SetBorderSize(1); TLegend *legHLTClass = new TLegend(.625,.625,.85,.89,"","NDC"); legHLTClass->SetFillColor(0); legHLTClass->SetBorderSize(1); THStack *hsElectronHLTEtaClass = new THStack("hsElectronHLTEtaClass","hsElectronHLTEtaClass"); THStack *hsElectronEtaClass = new THStack("hsElectronEtaClass","hsElectronEtaClass"); for (int ih=0; ih<9; ih++) { if (ih==4) continue; for (int ibin=0; ibin<hElectronHLTEta->GetNbinsX(); ibin++) { if (hElectronHLTEta->GetBinContent(ibin+1)==0) continue; if (hElectronEtaClass[ih]->GetBinContent(ibin+1)!=0) hElectronHLTEtaClass[ih]->SetBinContent( ibin+1, //hElectronHLTEtaClass[ih]->GetBinContent(ibin+1)/hElectronHLTEta->GetBinContent(ibin+1); hElectronHLTEtaClass[ih]->GetBinContent(ibin+1)/hElectronEtaClass[ih]->GetBinContent(ibin+1)); hsElectronHLTEtaClass->Add(hElectronHLTEtaClass[ih]); hElectronHLTEtaClass[ih]->SetLineColor( hColors[ih%5] ); hElectronEtaClass[ih]->SetBinContent( ibin+1, //hElectronHLTEtaClass[ih]->GetBinContent(ibin+1)/hElectronHLTEta->GetBinContent(ibin+1); hElectronEtaClass[ih]->GetBinContent(ibin+1)/hElectronEta->GetBinContent(ibin+1)); hsElectronEtaClass->Add(hElectronEtaClass[ih]); hElectronEtaClass[ih]->SetLineColor( hColors[ih%5] ); } if (ih<5) { sprintf(str,"%s",className[ih]); legClass->AddEntry(hElectronHLTEtaClass[ih],str,"l"); legHLTClass->AddEntry(hElectronHLTEtaClass[ih],str,"l"); } } hsElectronHLTEtaClass->Draw("nostack"); hsElectronHLTEtaClass->SetMaximum(1.2); hsElectronHLTEtaClass->GetXaxis()->SetTitle("#eta"); hsElectronHLTEtaClass->GetYaxis()->SetTitle("fraction of HLT electrons"); hsElectronHLTEtaClass->Draw("nostack"); legHLTClass->Draw(); sprintf(str,"%s/electronHLTEtaClass.%s",outDir,suffix); gPad->Print(str); hsElectronEtaClass->Draw("nostack"); hsElectronEtaClass->SetMaximum(1.2); hsElectronEtaClass->GetXaxis()->SetRangeUser(0,2.5); hsElectronEtaClass->GetXaxis()->SetTitle("#eta"); hsElectronEtaClass->GetYaxis()->SetTitle("fraction of electrons"); hsElectronEtaClass->Draw("nostack"); legClass->Draw(); sprintf(str,"%s/electronEtaClass.%s",outDir,suffix); gPad->Print(str); // HLT electrons and classes -- Pt function TFile *finElectronHLTPtClass = new TFile("histosFrame.root"); // eta distribution for electron classes - all electrons TH1F *hElectronHLTPt = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPt"); TH1F *hElectronHLTPtClass[9]; hElectronHLTPtClass[0] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass0"); hElectronHLTPtClass[1] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass1"); hElectronHLTPtClass[2] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass2"); hElectronHLTPtClass[3] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass3"); hElectronHLTPtClass[4] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass4"); hElectronHLTPtClass[5] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass10"); hElectronHLTPtClass[6] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass11"); hElectronHLTPtClass[7] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass12"); hElectronHLTPtClass[8] = (TH1F*)finElectronHLTPtClass->Get("hElectronHLTPtClass13"); TH1F *hElectronPtClass[9]; hElectronPtClass[0] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass0"); hElectronPtClass[1] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass1"); hElectronPtClass[2] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass2"); hElectronPtClass[3] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass3"); hElectronPtClass[4] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass4"); hElectronPtClass[5] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass10"); hElectronPtClass[6] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass11"); hElectronPtClass[7] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass12"); hElectronPtClass[8] = (TH1F*)finElectronHLTPtClass->Get("hElectronPtClass13"); cout << hElectronPtClass[0]->GetEntries() << endl; char className[5][64]; sprintf(className[0],"golden"); sprintf(className[1],"big brem"); sprintf(className[2],"narrow"); sprintf(className[3],"showering"); sprintf(className[4],"cracks"); int hColors[5] = { 3, 4, 6, 7, 2 }; TLegend *legClass = new TLegend(.625,.225,.85,.49,"","NDC"); legClass->SetFillColor(0); legClass->SetBorderSize(1); THStack *hsElectronHLTPtClass = new THStack("hsElectronHLTPtClass","hsElectronHLTPtClass"); for (int ih=0; ih<4; ih++) { for (int ibin=0; ibin<hElectronHLTPt->GetNbinsX(); ibin++) { if (hElectronHLTPt->GetBinContent(ibin+1)==0) continue; if (hElectronPtClass[ih]->GetBinContent(ibin+1)!=0 || hElectronPtClass[5+ih]->GetBinContent(ibin+1)!=0) hElectronHLTPtClass[ih]->SetBinContent( ibin+1, //hElectronHLTEtaClass[ih]->GetBinContent(ibin+1)/hElectronHLTEta->GetBinContent(ibin+1); (hElectronHLTPtClass[ih]->GetBinContent(ibin+1) + hElectronHLTPtClass[5+ih]->GetBinContent(ibin+1)) / ( hElectronPtClass[ih]->GetBinContent(ibin+1) + hElectronPtClass[5+ih]->GetBinContent(ibin+1))); hsElectronHLTPtClass->Add(hElectronHLTPtClass[ih]); hElectronHLTPtClass[ih]->SetLineColor( hColors[ih%5] ); hElectronHLTPtClass[ih]->SetMarkerColor( hColors[ih%5] ); hElectronHLTPtClass[ih]->SetMarkerStyle( 20 ); hElectronHLTPtClass[ih]->SetMarkerSize( .75 ); } if (ih<5) { sprintf(str,"%s",className[ih]); legClass->AddEntry(hElectronHLTPtClass[ih],str,"l"); } } hsElectronHLTPtClass->Draw("nostack"); hsElectronHLTPtClass->SetMaximum(1.2); hsElectronHLTPtClass->GetXaxis()->SetTitle("p_{T}"); hsElectronHLTPtClass->GetYaxis()->SetTitle("fraction of HLT electrons"); hsElectronHLTPtClass->Draw("nostack"); legClass->Draw(); sprintf(str,"%s/electronHLTPtClass.%s",outDir,suffix); gPad->Print(str); // E,p,comb/Etrue vs E/p gStyle->SetPalette(1); gStyle->SetOptStat(0); TFile *finElectronEP = new TFile("histosFrame.root"); TH1F *h2ElectronEOverEtrueEPClass[9]; h2ElectronEOverEtrueEPClass[0] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass0"); h2ElectronEOverEtrueEPClass[1] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass1"); h2ElectronEOverEtrueEPClass[2] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass2"); h2ElectronEOverEtrueEPClass[3] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass3"); h2ElectronEOverEtrueEPClass[4] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass4"); h2ElectronEOverEtrueEPClass[5] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass10"); h2ElectronEOverEtrueEPClass[6] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass11"); h2ElectronEOverEtrueEPClass[7] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass12"); h2ElectronEOverEtrueEPClass[8] = (TH1F*)finElectronEP->Get("h2ElectronEOverEtrueEPClass13"); TH1F *h2ElectronPOverEtrueEPClass[9]; h2ElectronPOverEtrueEPClass[0] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass0"); h2ElectronPOverEtrueEPClass[1] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass1"); h2ElectronPOverEtrueEPClass[2] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass2"); h2ElectronPOverEtrueEPClass[3] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass3"); h2ElectronPOverEtrueEPClass[4] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass4"); h2ElectronPOverEtrueEPClass[5] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass10"); h2ElectronPOverEtrueEPClass[6] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass11"); h2ElectronPOverEtrueEPClass[7] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass12"); h2ElectronPOverEtrueEPClass[8] = (TH1F*)finElectronEP->Get("h2ElectronPOverEtrueEPClass13"); TH1F *h2ElectronCombOverEtrueEPClass[9]; h2ElectronCombOverEtrueEPClass[0] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass0"); h2ElectronCombOverEtrueEPClass[1] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass1"); h2ElectronCombOverEtrueEPClass[2] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass2"); h2ElectronCombOverEtrueEPClass[3] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass3"); h2ElectronCombOverEtrueEPClass[4] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass4"); h2ElectronCombOverEtrueEPClass[5] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass10"); h2ElectronCombOverEtrueEPClass[6] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass11"); h2ElectronCombOverEtrueEPClass[7] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass12"); h2ElectronCombOverEtrueEPClass[8] = (TH1F*)finElectronEP->Get("h2ElectronCombOverEtrueEPClass13"); // barrel h2ElectronEOverEtrueEPClass[0]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[0]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[0]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[0]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EB_golden.%s",outDir,suffix); gPad->Print(str); h2ElectronEOverEtrueEPClass[1]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[1]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[1]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[1]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EB_bigBrem.%s",outDir,suffix); gPad->Print(str); h2ElectronEOverEtrueEPClass[2]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[2]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[2]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[2]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EB_narrow.%s",outDir,suffix); gPad->Print(str); h2ElectronEOverEtrueEPClass[3]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[3]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[3]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[3]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EB_showering.%s",outDir,suffix); gPad->Print(str); // h2ElectronEOverEtrueEPClass[0]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[0]->Add(h2ElectronEOverEtrueEPClass[2]); h2ElectronEOverEtrueEPClass[0]->Add(h2ElectronEOverEtrueEPClass[3]); h2ElectronEOverEtrueEPClass[0]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[0]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[0]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EB.%s",outDir,suffix); gPad->Print(str); h2ElectronPOverEtrueEPClass[0]->Add(h2ElectronPOverEtrueEPClass[1]); h2ElectronPOverEtrueEPClass[0]->Add(h2ElectronPOverEtrueEPClass[2]); h2ElectronPOverEtrueEPClass[0]->Add(h2ElectronPOverEtrueEPClass[3]); h2ElectronPOverEtrueEPClass[0]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronPOverEtrueEPClass[0]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronPOverEtrueEPClass[0]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/POverEtrue_vs_EOverP_EB.%s",outDir,suffix); gPad->Print(str); h2ElectronCombOverEtrueEPClass[0]->Add(h2ElectronCombOverEtrueEPClass[1]); h2ElectronCombOverEtrueEPClass[0]->Add(h2ElectronCombOverEtrueEPClass[2]); h2ElectronCombOverEtrueEPClass[0]->Add(h2ElectronCombOverEtrueEPClass[3]); h2ElectronCombOverEtrueEPClass[0]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronCombOverEtrueEPClass[0]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronCombOverEtrueEPClass[0]->GetYaxis()->SetTitle("#tilde{E}_{calo}/E_{kine}"); sprintf(str,"%s/combOverEtrue_vs_EOverP_EB.%s",outDir,suffix); gPad->Print(str); // endcap h2ElectronEOverEtrueEPClass[5]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[5]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[5]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[5]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EE_golden.%s",outDir,suffix); gPad->Print(str); h2ElectronEOverEtrueEPClass[6]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[6]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[6]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[6]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EE_bigBrem.%s",outDir,suffix); gPad->Print(str); h2ElectronEOverEtrueEPClass[7]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[7]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[7]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[7]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EE_narrow.%s",outDir,suffix); gPad->Print(str); h2ElectronEOverEtrueEPClass[8]->Add(h2ElectronEOverEtrueEPClass[1]); h2ElectronEOverEtrueEPClass[8]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[8]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[8]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EE_showering.%s",outDir,suffix); gPad->Print(str); // h2ElectronEOverEtrueEPClass[5]->Add(h2ElectronEOverEtrueEPClass[6]); h2ElectronEOverEtrueEPClass[5]->Add(h2ElectronEOverEtrueEPClass[6]); h2ElectronEOverEtrueEPClass[5]->Add(h2ElectronEOverEtrueEPClass[8]); h2ElectronEOverEtrueEPClass[5]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronEOverEtrueEPClass[5]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronEOverEtrueEPClass[5]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/EOverEtrue_vs_EOverP_EE.%s",outDir,suffix); gPad->Print(str); h2ElectronPOverEtrueEPClass[5]->Add(h2ElectronPOverEtrueEPClass[6]); h2ElectronPOverEtrueEPClass[5]->Add(h2ElectronPOverEtrueEPClass[6]); h2ElectronPOverEtrueEPClass[5]->Add(h2ElectronPOverEtrueEPClass[8]); h2ElectronPOverEtrueEPClass[5]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronPOverEtrueEPClass[5]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronPOverEtrueEPClass[5]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/POverEtrue_vs_EOverP_EE.%s",outDir,suffix); gPad->Print(str); h2ElectronCombOverEtrueEPClass[5]->Add(h2ElectronCombOverEtrueEPClass[6]); h2ElectronCombOverEtrueEPClass[5]->Add(h2ElectronCombOverEtrueEPClass[6]); h2ElectronCombOverEtrueEPClass[5]->Add(h2ElectronCombOverEtrueEPClass[8]); h2ElectronCombOverEtrueEPClass[5]->Draw("colz"); gPad->SetRightMargin(0.15); h2ElectronCombOverEtrueEPClass[5]->GetXaxis()->SetTitle("E_{calo}/p_{tk}"); h2ElectronCombOverEtrueEPClass[5]->GetYaxis()->SetTitle("E_{calo}/E_{kine}"); sprintf(str,"%s/combOverEtrue_vs_EOverP_EE.%s",outDir,suffix); gPad->Print(str); /* for (int ih=0; ih<4; ih++) { for (int ibin=0; ibin<hElectronHLTPt->GetNbinsX(); ibin++) { if (hElectronHLTPt->GetBinContent(ibin+1)==0) continue; if (hElectronPtClass[ih]->GetBinContent(ibin+1)!=0 || hElectronPtClass[5+ih]->GetBinContent(ibin+1)!=0) hElectronHLTPtClass[ih]->SetBinContent( ibin+1, //hElectronHLTEtaClass[ih]->GetBinContent(ibin+1)/hElectronHLTEta->GetBinContent(ibin+1); (hElectronHLTPtClass[ih]->GetBinContent(ibin+1) + hElectronHLTPtClass[5+ih]->GetBinContent(ibin+1)) / ( hElectronPtClass[ih]->GetBinContent(ibin+1) + hElectronPtClass[5+ih]->GetBinContent(ibin+1))); hsElectronHLTPtClass->Add(hElectronHLTPtClass[ih]); hElectronHLTPtClass[ih]->SetLineColor( hColors[ih%5] ); hElectronHLTPtClass[ih]->SetMarkerColor( hColors[ih%5] ); hElectronHLTPtClass[ih]->SetMarkerStyle( 20 ); hElectronHLTPtClass[ih]->SetMarkerSize( .75 ); } if (ih<5) { sprintf(str,"%s",className[ih]); legClass->AddEntry(hElectronHLTPtClass[ih],str,"l"); } } */ }

---