void QCDHT_cutflow_weighted()
{
//=========Macro generated from canvas: c1/c1
//========= (Wed Feb 24 17:48:49 2016) by ROOT version6.02/05
TCanvas *c1 = new TCanvas("c1", "c1",65,52,700,500);
c1->Range(-1.25,-2351156,16.25,2.11604e+07);
c1->SetFillColor(0);
c1->SetBorderMode(0);
c1->SetBorderSize(2);
c1->SetFrameBorderMode(0);
c1->SetFrameBorderMode(0);
TH1D *cutflow1 = new TH1D("cutflow1","cut flow",14,0.5,14.5);
cutflow1->SetBinContent(1,1.791357e+07);
cutflow1->SetBinContent(2,1.791357e+07);
cutflow1->SetBinContent(3,1.482016e+07);
cutflow1->SetBinContent(4,1.397007e+07);
cutflow1->SetBinContent(5,8159947);
cutflow1->SetBinContent(6,1368453);
cutflow1->SetBinContent(7,12035.83);
cutflow1->SetBinContent(8,10653.79);
cutflow1->SetBinContent(9,4100.401);
cutflow1->SetBinContent(10,4323.004);
cutflow1->SetBinContent(11,1820.803);
cutflow1->SetBinContent(12,374.2397);
cutflow1->SetBinContent(13,35.33991);
cutflow1->SetBinContent(14,214.6597);
cutflow1->SetEntries(1.137194e+08);
TPaveStats *ptstats = new TPaveStats(0.78,0.775,0.98,0.935,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("cutflow");
AText->SetTextSize(0.0368);
AText = ptstats->AddText("Entries = 1.137194e+08");
AText = ptstats->AddText("Mean = 0");
AText = ptstats->AddText("RMS = 0");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
ptstats->Draw();
cutflow1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(cutflow1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
cutflow1->SetLineColor(ci);
cutflow1->GetXaxis()->SetBinLabel(1,"All");
cutflow1->GetXaxis()->SetBinLabel(2,"Trigger");
cutflow1->GetXaxis()->SetBinLabel(3,"p_{T}");
cutflow1->GetXaxis()->SetBinLabel(4,"|#eta|");
cutflow1->GetXaxis()->SetBinLabel(5,"|#Delta#eta|");
cutflow1->GetXaxis()->SetBinLabel(6,"M(jet_{0},jet(1))");
cutflow1->GetXaxis()->SetBinLabel(7,"M(jets)");
cutflow1->GetXaxis()->SetBinLabel(8,"#tau_{21}");
cutflow1->GetXaxis()->SetBinLabel(9,"0b");
cutflow1->GetXaxis()->SetBinLabel(10,"1b");
cutflow1->GetXaxis()->SetBinLabel(11,"2b");
cutflow1->GetXaxis()->SetBinLabel(12,"3b");
cutflow1->GetXaxis()->SetBinLabel(13,"4b");
cutflow1->GetXaxis()->SetBinLabel(14,"3b+HPHP");
cutflow1->GetXaxis()->SetLabelFont(42);
cutflow1->GetXaxis()->SetLabelSize(0.035);
cutflow1->GetXaxis()->SetTitleSize(0.035);
cutflow1->GetXaxis()->SetTitleFont(42);
cutflow1->GetYaxis()->SetLabelFont(42);
cutflow1->GetYaxis()->SetLabelSize(0.035);
cutflow1->GetYaxis()->SetTitleSize(0.035);
cutflow1->GetYaxis()->SetTitleFont(42);
cutflow1->GetZaxis()->SetLabelFont(42);
cutflow1->GetZaxis()->SetLabelSize(0.035);
cutflow1->GetZaxis()->SetTitleSize(0.035);
cutflow1->GetZaxis()->SetTitleFont(42);
cutflow1->Draw("");
TPaveText *pt = new TPaveText(0.4246264,0.94,0.5753736,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("cut flow");
pt->Draw();
c1->Modified();
c1->cd();
c1->SetSelected(c1);
}
void fixStats()
{
TH1* ratioHist = (TH1*) gROOT->FindObject("data");
TPaveStats* stats = (TPaveStats*) ratioHist->FindObject("stats");
if(stats)stats->SetTextSize(0.040209);
}
// -----------------------------------------------------------------------------
// Create plots
void plot( std::string& path,
std::string& type,
std::string& sample,
std::string& dir,
std::string& histo,
double lumi,
bool scale ) {
std::string canvas_name = histo + "_" + dir + "_" + type + "_" + sample;
// Create canvas
TCanvas* canvas = new TCanvas(canvas_name.c_str(),"");
canvas->SetFillColor(0);
canvas->SetFrameBorderMode(0);
canvas->SetFrameFillColor(0);
canvas->SetTopMargin(0.10);
canvas->SetBottomMargin(0.12);
canvas->SetLeftMargin(0.12);
canvas->SetRightMargin(0.15);
// Retrieve histogram
std::string file_name = path + type + "_" + sample + ".root";
TFile* f = new TFile(file_name.c_str(),"READ");
TDirectory* d = (TDirectory*)f->Get(dir.c_str());
TH2D* his = (TH1*)d->Get(histo.c_str());
if ( !his ) return;
//his->Rebin2D(2,2);
if ( true ) { gPad->SetLogz(); }
if ( scale ) his->Scale( lumi / 100. );
his->SetMaximum(1.e4);
his->SetMinimum(1.e-5);
//his->SetMinimum( his->GetMinimum(1.e-12) );
// his->SetMaximum( 20000. );
// his->SetMinimum( 2.e-4 );
// his->SetMaximum( 20000. );
// his->SetMinimum( 20. );
double xmin = his->GetXaxis()->GetXmin();
double xmax = his->GetXaxis()->GetXmax();
double ymin = his->GetYaxis()->GetXmin();
double ymax = his->GetYaxis()->GetXmax();
// Reset title
std::string title = ";" + std::string(his->GetXaxis()->GetTitle()) + ";" + std::string(his->GetYaxis()->GetTitle());
his->SetTitle(title.c_str());
his->GetXaxis()->SetTitle("x_{2}");
his->GetXaxis()->SetTitleOffset(1.2);
his->GetYaxis()->SetTitle("x_{1}");
his->GetYaxis()->SetTitleOffset(1.4);
his->Draw("COLZ");
gPad->Update();
// Lumi
if (1) {
std::stringstream ss;
ss << "#int L dt = " << lumi << " pb^{-1}";
double xpos = 0.05 * (xmax-xmin)+xmin;
double ypos = 0.25 * (ymax-ymin)+ymin;
TLatex* text1 = new TLatex(xpos,ypos,ss.str().c_str());
text1->SetTextAlign(12);
text1->SetTextSize(0.035);
text1->Draw();
}
// Jet type
if (1) {
double xpos = 0.05 * (xmax-xmin)+xmin;
double ypos = 0.15 * (ymax-ymin)+ymin;
TText* text2 = new TText(xpos,ypos,type.c_str());
text2->SetTextAlign(12);
text2->SetTextSize(0.035);
text2->Draw();
}
// Sample
if (1) {
double xpos = 0.05 * (xmax-xmin)+xmin;
double ypos = 0.10 * (ymax-ymin)+ymin;
TText* text3 = new TText(xpos,ypos,sample.c_str());
text3->SetTextAlign(12);
text3->SetTextSize(0.035);
text3->Draw();
}
// Stats
gStyle->SetOptStat("i");
his->SetStats(1);
TPaveStats* stats = (TPaveStats*)his->GetListOfFunctions()->FindObject("stats");
std::string stats_pos = "br";
if ( stats ) {
stats->SetFillColor(0);
stats->SetLineColor(0);
stats->SetShadowColor(0);
if ( stats_pos == "tr" ) {
stats->SetX1NDC(0.60); stats->SetY1NDC(0.68); stats->SetX2NDC(0.83); stats->SetY2NDC(0.88);
} else if ( stats_pos == "br" ) {
stats->SetX1NDC(0.60); stats->SetY1NDC(0.18); stats->SetX2NDC(0.83); stats->SetY2NDC(0.28);
} else {
stats->SetX1NDC(0.60); stats->SetY1NDC(0.68); stats->SetX2NDC(0.83); stats->SetY2NDC(0.88);
}
}
// Scale
gStyle->SetPalette(1);
TPaletteAxis* palette = (TPaletteAxis*)his->GetListOfFunctions()->FindObject("palette");
if ( palette ) {
palette->SetY1NDC(0.2);
palette->SetY2NDC(0.8);
}
canvas->Modified();
canvas->cd();
canvas->SetSelected(canvas);
canvas->SaveAs(std::string(canvas_name+".png").c_str());
//canvas->Write();
}
void SinglePi(const TString ref_vers, const TString val_vers, bool fastsim){
//Warning!!! This rootlogon hacks the root color pallate
setColors();
TString ref_file = "pi50scan"+ref_vers+"_fullGeom_ECALHCAL_CaloTowers.root";
TString val_file = "pi50scan"+val_vers+"_fullGeom_ECALHCAL_CaloTowers.root";
TFile f1(ref_file);
TFile f2(val_file);
// service variables
//
//Profiles
const int Nprof = 12;
TProfile* f1_prof[Nprof];
TProfile* f2_prof[Nprof];
char labelp[Nprof][64];
//1D Histos
const int Nhist1 = 7;
TH1* f1_hist1[Nhist1];
TH1* f2_hist1[Nhist1];
char label1[Nhist1][64];
//Labels
//Profiles
sprintf(labelp[0], "CaloTowersTask_emean_vs_ieta_E1.gif");
sprintf(labelp[1], "CaloTowersTask_emean_vs_ieta_H1.gif");
sprintf(labelp[2], "CaloTowersTask_emean_vs_ieta_EH1.gif");
sprintf(labelp[3], "RecHitsTask_emean_vs_ieta_E.gif");
sprintf(labelp[4], "RecHitsTask_emean_vs_ieta_H.gif");
sprintf(labelp[5], "RecHitsTask_emean_vs_ieta_EH.gif");
if (!fastsim) {
sprintf(labelp[6], "SimHitsTask_emean_vs_ieta_E.gif");
sprintf(labelp[7], "SimHitsTask_emean_vs_ieta_H.gif");
sprintf(labelp[8], "SimHitsTask_emean_vs_ieta_EH.gif");
}
sprintf(labelp[9], "RecHitsTask_timing_vs_energy_profile_HB.gif");
sprintf(labelp[10], "RecHitsTask_timing_vs_energy_profile_HE.gif");
sprintf(labelp[11], "RecHitsTask_timing_vs_energy_profile_HF.gif");
//1D Histos
sprintf(label1[0], "N_calotowers_HB.gif");
sprintf(label1[1], "N_calotowers_HE.gif");
sprintf(label1[2], "N_calotowers_HF.gif");
sprintf(label1[3], "RecHits_energy_HB.gif");
sprintf(label1[4], "RecHits_energy_HE.gif");
sprintf(label1[5], "RecHits_energy_HO.gif");
sprintf(label1[6], "RecHits_energy_HF.gif");
TDirectory *td = fileDirectory(&f1, "CaloTowersTask");
//f1.cd("DQMData/CaloTowersV/CaloTowersTask");
//gDirectory->pwd();
td->pwd();
f1_prof[0] = (TProfile*)td->Get("emean_vs_ieta_E1");
f1_prof[1] = (TProfile*)td->Get("emean_vs_ieta_H1");
f1_prof[2] = (TProfile*)td->Get("emean_vs_ieta_EH1");
f1_hist1[0] = (TH1*)td->Get("CaloTowersTask_number_of_fired_towers_HB");
f1_hist1[1] = (TH1*)td->Get("CaloTowersTask_number_of_fired_towers_HE");
f1_hist1[2] = (TH1*)td->Get("CaloTowersTask_number_of_fired_towers_HF");
td = fileDirectory(&f1, "HcalRecHitTask");
//f1.cd("DQMData/HcalRecHitsV/HcalRecHitTask");
f1_prof[3] = (TProfile*)td->Get("HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_E");
f1_prof[4] = (TProfile*)td->Get("HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths");
f1_prof[5] = (TProfile*)td->Get("HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_EH");
f1_prof[9] = (TProfile*)td->Get("HcalRecHitTask_timing_vs_energy_profile_HB");
f1_prof[10] = (TProfile*)td->Get("HcalRecHitTask_timing_vs_energy_profile_Low_HE");
f1_prof[11] = (TProfile*)td->Get("HcalRecHitTask_timing_vs_energy_profile_Low_HF");
f1_hist1[3] = (TH1*)td->Get("HcalRecHitTask_energy_of_rechits_HB");
f1_hist1[4] = (TH1*)td->Get("HcalRecHitTask_energy_of_rechits_HE");
f1_hist1[5] = (TH1*)td->Get("HcalRecHitTask_energy_of_rechits_HO");
f1_hist1[6] = (TH1*)td->Get("HcalRecHitTask_energy_of_rechits_HF");
if (!fastsim) {
td = fileDirectory(&f1, "HcalSimHitTask");
//f1.cd("DQMData/HcalSimHitsV/HcalSimHitTask");
f1_prof[6] = (TProfile*)td->Get("HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths_E");
f1_prof[7] = (TProfile*)td->Get("HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths");
f1_prof[8] = (TProfile*)td->Get("HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths_EH");
}
td = fileDirectory(&f2, "CaloTowersTask");
//f2.cd("DQMData/CaloTowersV/CaloTowersTask");
//gDirectory->pwd();
td->pwd();
f2_prof[0] = (TProfile*)td->Get("emean_vs_ieta_E1");
f2_prof[1] = (TProfile*)td->Get("emean_vs_ieta_H1");
f2_prof[2] = (TProfile*)td->Get("emean_vs_ieta_EH1");
f2_hist1[0] = (TH1*)td->Get("CaloTowersTask_number_of_fired_towers_HB");
f2_hist1[1] = (TH1*)td->Get("CaloTowersTask_number_of_fired_towers_HE");
f2_hist1[2] = (TH1*)td->Get("CaloTowersTask_number_of_fired_towers_HF");
td = fileDirectory(&f2, "HcalRecHitTask");
//f2.cd("DQMData/HcalRecHitsV/HcalRecHitTask");
f2_prof[3] = (TProfile*)td->Get("HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_E");
f2_prof[4] = (TProfile*)td->Get("HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths");
f2_prof[5] = (TProfile*)td->Get("HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_EH");
f2_prof[9] = (TProfile*)td->Get("HcalRecHitTask_timing_vs_energy_profile_HB");
f2_prof[10] = (TProfile*)td->Get("HcalRecHitTask_timing_vs_energy_profile_Low_HE");
f2_prof[11] = (TProfile*)td->Get("HcalRecHitTask_timing_vs_energy_profile_Low_HF");
f2_hist1[3] = (TH1*)td->Get("HcalRecHitTask_energy_of_rechits_HB");
f2_hist1[4] = (TH1*)td->Get("HcalRecHitTask_energy_of_rechits_HE");
f2_hist1[5] = (TH1*)td->Get("HcalRecHitTask_energy_of_rechits_HO");
f2_hist1[6] = (TH1*)td->Get("HcalRecHitTask_energy_of_rechits_HF");
if (!fastsim) {
td = fileDirectory(&f2, "HcalSimHitTask");
//f2.cd("DQMData/HcalSimHitsV/HcalSimHitTask");
f2_prof[6] = (TProfile*)td->Get("HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths_E");
f2_prof[7] = (TProfile*)td->Get("HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths");
f2_prof[8] = (TProfile*)td->Get("HcalSimHitTask_En_simhits_cone_profile_vs_ieta_all_depths_EH");
}
//Profiles
f1_prof[0]->GetXaxis()->SetTitle("CaloTowers eE (GeV) vs ieta 1 Tower");
f1_prof[1]->GetXaxis()->SetTitle("CaloTowers hE (GeV) vs ieta 1 Tower");
f1_prof[2]->GetXaxis()->SetTitle("CaloTowers eE+hE (GeV) vs ieta 1 Tower");
f1_prof[3]->GetXaxis()->SetTitle("RecHits eE (GeV) vs ieta R = 0.3 Cone");
f1_prof[4]->GetXaxis()->SetTitle("RecHits hE (GeV) vs ieta R = 0.3 Cone");
f1_prof[5]->GetXaxis()->SetTitle("RecHits eE+hE (GeV) vs ieta R = 0.3 Cone");
if (!fastsim) {
f1_prof[6]->GetXaxis()->SetTitle("SimHits eE (GeV) vs ieta R = 0.3 Cone");
f1_prof[7]->GetXaxis()->SetTitle("SimHits hE (GeV) vs ieta R = 0.3 Cone");
f1_prof[8]->GetXaxis()->SetTitle("SimHits eE+hE (GeV) vs ieta R = 0.3 Cone");
}
f1_prof[9]->GetXaxis()->SetTitle("HB RecHits timing (ns) vs Energy (GeV)");
f1_prof[10]->GetXaxis()->SetTitle("HE RecHits timing (ns) vs Energy (GeV)");
f1_prof[11]->GetXaxis()->SetTitle("HF RecHits timing (ns) vs Energy (GeV)");
//1D Histos
f1_hist1[0]->GetXaxis()->SetTitle("Number of HB CaloTowers");
f1_hist1[1]->GetXaxis()->SetTitle("Number of HE CaloTowers");
f1_hist1[2]->GetXaxis()->SetTitle("Number of HF CaloTowers");
f1_hist1[3]->GetXaxis()->SetTitle("HB RecHits energy (GeV)");
f1_hist1[4]->GetXaxis()->SetTitle("HE RecHits energy (GeV)");
f1_hist1[5]->GetXaxis()->SetTitle("HO RecHits energy (GeV)");
f1_hist1[6]->GetXaxis()->SetTitle("HF RecHits energy (GeV)");
//
f1_prof[0]->SetMaximum(20.);
f1_prof[1]->SetMaximum(40.);
f1_prof[2]->SetMaximum(40.);
f1_prof[4]->SetMaximum(50.);
f1_prof[5]->SetMaximum(50.);
f1_prof[0]->SetMinimum(0.);
f1_prof[1]->SetMinimum(0.);
f1_prof[2]->SetMinimum(0.);
f1_prof[3]->SetMinimum(0.);
f1_prof[4]->SetMinimum(0.);
f1_prof[5]->SetMinimum(0.);
if (!fastsim) {
f1_prof[6]->SetMinimum(0.);
f1_prof[7]->SetMinimum(0.);
f1_prof[8]->SetMinimum(0.);
}
f1_prof[9]->GetXaxis()->SetRangeUser(0.,75.);
f1_prof[10]->GetXaxis()->SetRangeUser(0.,75.);
f1_prof[11]->GetXaxis()->SetRangeUser(0.,75.);
/*
f1_prof[9]->SetMinimum(0.);
f1_prof[10]->SetMinimum(0.);
f1_prof[11]->SetMinimum(0.);
*/
// f1_hist[2]->GetXaxis()->SetRangeUser(0.,1200.);
// f1_hist[7]->GetXaxis()->SetRangeUser(0.,160.);
// hist1->GetXaxis()->SetNdivisions(-21);
// hist1->GetYaxis()->SetNdivisions(-1003);
f1_hist1[0]->GetXaxis()->SetRangeUser(0.,100.);
f2_hist1[0]->GetXaxis()->SetRangeUser(0.,100.);
f1_hist1[1]->GetXaxis()->SetRangeUser(0.,150.);
f2_hist1[1]->GetXaxis()->SetRangeUser(0.,150.);
f1_hist1[2]->GetXaxis()->SetRangeUser(0.,100.);
f2_hist1[2]->GetXaxis()->SetRangeUser(0.,100.);
f1_hist1[3]->GetXaxis()->SetRangeUser(0.,100.);
f2_hist1[3]->GetXaxis()->SetRangeUser(0.,100.);
f1_hist1[4]->GetXaxis()->SetRangeUser(0.,100.);
f2_hist1[4]->GetXaxis()->SetRangeUser(0.,100.);
f1_hist1[5]->GetXaxis()->SetRangeUser(0.,100.);
f2_hist1[5]->GetXaxis()->SetRangeUser(0.,100.);
f1_hist1[6]->GetXaxis()->SetRangeUser(0.,100.);
f2_hist1[6]->GetXaxis()->SetRangeUser(0.,100.);
f1_hist1[3]->SetMaximum(1.e8);
f1_hist1[4]->SetMaximum(1.e8);
f1_hist1[5]->SetMaximum(1.e8);
f1_hist1[6]->SetMaximum(1.e8);
// 1D-histo
for (int i = 0; i < Nhist1; i++){
TCanvas *myc = new TCanvas("myc","",800,600);
gStyle->SetOptStat(1111);
if(i > 2) myc->SetLogy();
f1_hist1[i]->SetStats(kTRUE); // stat box
f2_hist1[i]->SetStats(kTRUE);
f1_hist1[i]->SetTitle("");
f2_hist1[i]->SetTitle("");
f1_hist1[i]->SetLineWidth(2);
f2_hist1[i]->SetLineWidth(2);
// diffferent histo colors and styles
f1_hist1[i]->SetLineColor(41);
f1_hist1[i]->SetLineStyle(1);
f2_hist1[i]->SetLineColor(43);
f2_hist1[i]->SetLineStyle(2);
//Set maximum to the larger of the two
if (f1_hist1[i]->GetMaximum() < f2_hist1[i]->GetMaximum()) f1_hist1[i]->SetMaximum(1.05 * f2_hist1[i]->GetMaximum());
TLegend *leg = new TLegend(0.2, 0.91, 0.6, 0.99, "","brNDC");
leg->SetBorderSize(2);
// leg->SetFillColor(51); // see new color definition above
leg->SetFillStyle(1001); //
leg->AddEntry(f1_hist1[i],"CMSSW_"+ref_vers,"l");
leg->AddEntry(f2_hist1[i],"CMSSW_"+val_vers,"l");
TPaveStats *ptstats = new TPaveStats(0.85,0.86,0.98,0.98,"brNDC");
ptstats->SetTextColor(41);
f1_hist1[i]->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(f1_hist1[i]->GetListOfFunctions());
TPaveStats *ptstats2 = new TPaveStats(0.85,0.74,0.98,0.86,"brNDC");
ptstats2->SetTextColor(43);
f2_hist1[i]->GetListOfFunctions()->Add(ptstats2);
ptstats2->SetParent(f2_hist1[i]->GetListOfFunctions());
f1_hist1[i]->Draw(""); // "stat"
f2_hist1[i]->Draw("histsames");
leg->Draw();
myc->SaveAs(label1[i]);
if(myc) delete myc;
}
// Profiles
for (int i = 0; i < Nprof; i++){
TCanvas *myc = new TCanvas("myc","",800,600);
bool skipHisto = false;
if (fastsim && i>=6 && i<=8) skipHisto = true;
if (!skipHisto) {
f1_prof[i]->SetStats(kFALSE);
f2_prof[i]->SetStats(kFALSE);
f1_prof[i]->SetTitle("");
f2_prof[i]->SetTitle("");
f1_prof[i]->SetLineColor(41);
f1_prof[i]->SetLineStyle(1);
f1_prof[i]->SetLineWidth(1);
f1_prof[i]->SetMarkerColor(41);
f1_prof[i]->SetMarkerStyle(21);
f1_prof[i]->SetMarkerSize(1.0);
f2_prof[i]->SetLineColor(43);
f2_prof[i]->SetLineStyle(1);
f2_prof[i]->SetLineWidth(1);
f2_prof[i]->SetMarkerColor(43);
f2_prof[i]->SetMarkerStyle(20);
f2_prof[i]->SetMarkerSize(0.8);
if(i > 8 ) {
f1_prof[i]->SetMarkerSize(0.1);
f2_prof[i]->SetMarkerSize(0.3);
}
myc->SetGrid();
if( i <= 8) {
f1_prof[i]->Draw("histpl");
f2_prof[i]->Draw("histplsame"); // esame
}
else {
f1_prof[i]->Draw("pl");
f2_prof[i]->Draw("plsame"); // esame
}
TLegend *leg = new TLegend(0.40, 0.91, 0.74, 0.99, "","brNDC");
leg->SetBorderSize(2);
leg->SetFillStyle(1001);
leg->AddEntry(f1_prof[i],"CMSSW_"+ref_vers,"pl");
leg->AddEntry(f2_prof[i],"CMSSW_"+val_vers,"pl");
leg->Draw("");
myc->SaveAs(labelp[i]);
}
if(myc) delete myc;
}
TCanvas *myc = new TCanvas("myc","",800,600);
TProfile* ratio1 = (TProfile*)f2_prof[2]->Clone();
ratio1->Divide(f1_prof[2]);
ratio1->SetMaximum(1.2);
ratio1->SetMinimum(0.8);
myc->SetGrid();
ratio1->Draw("hist pl");
TLegend *leg = new TLegend(0.20, 0.91, 0.70, 0.99, "","brNDC");
leg->SetBorderSize(2);
leg->SetFillStyle(1001);
leg->AddEntry(ratio1,"CaloTowers scale (pi50) ratio "+val_vers+"/"+ref_vers+" vs ieta","pl");
leg->Draw("");
// f1_prof[2]->Draw();
myc->SaveAs("Ratio.gif");
// close ROOT files
//
f1.Close() ;
f2.Close() ;
return ;
}
void many_plots2(int first = 5808, int last = 11688, int mask = 0x801E)
{
TH1D *h[12][2];
int i, j;
char str[64];
TCanvas *cv[2];
TPaveStats *st;
float y, dy;
gROOT->SetStyle("Plain");
gStyle->SetOptStat(1000000);
gStyle->SetOptFit();
TFile *fRoot = new TFile("many_plots.root", "RECREATE");
HPainter *p = new HPainter(mask, first, last);
p->SetFile(fRoot);
if (!p->IsOpen()) {
printf("Something wrong with data files.\n");
return;
}
TCut cVeto("gtFromVeto > 60");
TCut cIso("(gtFromPrevious > 45 && gtToNext > 80 && EventsBetween == 0) || (gtFromPrevious == gtFromVeto)");
TCut cX("PositronX[0] < 0 || (PositronX[0] > 2 && PositronX[0] < 94)");
TCut cY("PositronX[1] < 0 || (PositronX[1] > 2 && PositronX[1] < 94)");
TCut cZ("PositronX[2] > 3.5 && PositronX[2] < 95.5");
TCut cXYZ = cX && cY && cZ;
TCut cR("Distance < 100 && DistanceZ > -40 && DistanceZ < 40");
TCut cT10("gtDiff > 1");
TCut cT20("gtDiff > 2");
TCut cT200("gtDiff < 20"); // strong cut
TCut cGamma("AnnihilationEnergy < 1.5 && AnnihilationGammas < 9");
TCut cPe("PositronEnergy > 1");
TCut cXY("PositronX[0]>=0 && PositronX[1]>=0");
TCut cN4("NeutronEnergy > 4 && NeutronHits >= 5");
for (i=0; i<2; i++) {
sprintf(str, "hR%d", i);
h[0][i] = new TH1D(str, "Distance between positron and neutron;R, cm;mHz", 40, 0, 160);
sprintf(str, "hRZ%d", i);
h[1][i] = new TH1D(str, "Distance between positron and neutron, projection Z;R_{z}, cm;mHz", 100, -100, 100);
sprintf(str, "hT%d", i);
h[2][i] = new TH1D(str, "Time between positron and neutron;T, us;mHz", 50, 0, 50);
sprintf(str, "hX%d", i);
h[3][i] = new TH1D(str, "Positron vertex X;X, cm;mHz", 25, 0, 100);
sprintf(str, "hY%d", i);
h[4][i] = new TH1D(str, "Positron vertex Y;Y, cm;mHz", 25, 0, 100);
sprintf(str, "hZ%d", i);
h[5][i] = new TH1D(str, "Positron vertex Z;Z, cm;mHz", 100, 0, 100);
sprintf(str, "hNE%d", i);
h[6][i] = new TH1D(str, "Energy detected in neutron capture;E_{n}, MeV;mHz", 50, 0, 10);
sprintf(str, "hNN%d", i);
h[7][i] = new TH1D(str, "Number of hits in SiPM for neutron capture;N_{n};mHz", 20, 0, 20);
sprintf(str, "hGE%d", i);
h[8][i] = new TH1D(str, "Energy beyond positron cluster;E_{#gamma}, MeV;mHz", 15, 0, 3);
sprintf(str, "hGN%d", i);
h[9][i] = new TH1D(str, "Number of SiPM hits out of positron cluster;N_{#gamma};mHz", 10, 0, 10);
sprintf(str, "hE%d", i);
h[10][i] = new TH1D(str, "Positorn kinetic energy;E, MeV;mHz", 40, 0, 8);
sprintf(str, "hEC%d", i);
h[11][i] = new TH1D(str, "Positorn kinetic energy, strong background cuts;E, MeV;mHz", 40, 0, 8);
}
printf("Histograms are created\n");
p->Project(h[0][0], "Distance", cVeto && cIso && cT20 && cGamma && cPe && cXYZ);
p->Project(h[0][1], "Distance", !cVeto && cIso && cT20 && cGamma && cPe && cXYZ);
printf("Distance.\n");
p->Project(h[1][0], "DistanceZ", cVeto && cIso && cT20 && cGamma && cPe && cXYZ);
p->Project(h[1][1], "DistanceZ", !cVeto && cIso && cT20 && cGamma && cPe && cXYZ);
printf("DistanceZ.\n");
p->Project(h[2][0], "gtDiff", cVeto && cIso && cT10 && cGamma && cPe && cXYZ && cR);
p->Project(h[2][1], "gtDiff", !cVeto && cIso && cT10 && cGamma && cPe && cXYZ && cR);
printf("gtDiff.\n");
p->Project(h[3][0], "PositronX[0]+2", cVeto && cIso && cT20 && cGamma && cPe && cY && cZ && cR && "PositronX[0] >= 0");
p->Project(h[3][1], "PositronX[0]+2", !cVeto && cIso && cT20 && cGamma && cPe && cY && cZ && cR && "PositronX[0] >= 0");
printf("X.\n");
p->Project(h[4][0], "PositronX[1]+2", cVeto && cIso && cT20 && cGamma && cPe && cX && cZ && cR && "PositronX[1] >= 0");
p->Project(h[4][1], "PositronX[1]+2", !cVeto && cIso && cT20 && cGamma && cPe && cX && cZ && cR && "PositronX[1] >= 0");
printf("Y.\n");
p->Project(h[5][0], "PositronX[2]+0.5", cVeto && cIso && cT20 && cGamma && cPe && cX && cY && cR);
p->Project(h[5][1], "PositronX[2]+0.5", !cVeto && cIso && cT20 && cGamma && cPe && cX && cY && cR);
printf("Z.\n");
p->Project(h[6][0], "NeutronEnergy", cVeto && cIso && cT20 && cGamma && cPe && cY && cZ && cR);
p->Project(h[6][1], "NeutronEnergy", !cVeto && cIso && cT20 && cGamma && cPe && cY && cZ && cR);
printf("NE.\n");
p->Project(h[7][0], "NeutronHits", cVeto && cIso && cT20 && cGamma && cPe && cXYZ && cR);
p->Project(h[7][1], "NeutronHits", !cVeto && cIso && cT20 && cGamma && cPe && cXYZ && cR);
printf("NN.\n");
p->Project(h[8][0], "AnnihilationEnergy", cVeto && cIso && cT20 && cPe && cXYZ && cR);
p->Project(h[8][1], "AnnihilationEnergy", !cVeto && cIso && cT20 && cPe && cXYZ && cR);
printf("AE.\n");
p->Project(h[9][0], "AnnihilationGammas", cVeto && cIso && cT20 && cPe && cXYZ && cR);
p->Project(h[9][1], "AnnihilationGammas", !cVeto && cIso && cT20 && cPe && cXYZ && cR);
printf("AG.\n");
p->Project(h[10][0], "PositronEnergy", cVeto && cIso && cT20 && cGamma && cPe && cXYZ && cR);
p->Project(h[10][1], "PositronEnergy", !cVeto && cIso && cT20 && cGamma && cPe && cXYZ && cR);
p->Project(h[11][0], "PositronEnergy", cVeto && cIso && cT20 && cGamma && cPe && cXYZ && cR && cN4 && cT200);
p->Project(h[11][1], "PositronEnergy", !cVeto && cIso && cT20 && cGamma && cPe && cXYZ && cR && cN4 && cT200);
printf("Projections are done\n");
for (i=0; i<12; i++) {
for (j=0; j<2; j++) h[i][j]->SetLineWidth(4);
h[i][0]->SetLineColor(kGreen);
h[i][1]->SetLineColor(kRed);
h[i][0]->SetMinimum(0);
h[i][1]->SetMinimum(0);
}
// h[1][0]->Fit("gaus", "0");
// h[1][1]->Fit("gaus", "0");
// TF1 *fdec = new TF1("FDEC", "[0]*(exp(-x/[1]) - exp(-x/[2]))", 3);
// fdec->SetParNames("Const", "t_{CAPTURE}", "t_{THERM}");
// fdec->SetParameters(h[2][0]->Integral()/4, 15, 5);
// h[2][0]->Fit("FDEC", "0", "", 2, 50);
for (i=0; i<2; i++) {
sprintf(str, "CV%d", i);
cv[i] = new TCanvas(str, "Plots", 1800, 1200);
cv[i]->Divide(3, 2);
for (j=0; j<6; j++) {
cv[i]->cd(j+1);
h[6*i+j][0]->Draw();
h[6*i+j][1]->Draw("sames");
gPad->Update();
if (i == 0 && j == 1) {
st = (TPaveStats *) h[6*i+j][0]->FindObject("stats");
st->SetOptStat(1100);
st->SetLineColor(kGreen);
st->SetTextColor(kGreen);
y = st->GetY1NDC();
dy = st->GetY2NDC() - y;
st->SetX1NDC(0.72);
st = (TPaveStats *) h[6*i+j][1]->FindObject("stats");
st->SetOptStat(1100);
st->SetLineColor(kRed);
st->SetTextColor(kRed);
st->SetX1NDC(0.72);
st->SetY2NDC(y);
st->SetY1NDC(y - dy);
} else {
st = (TPaveStats *) h[6*i+j][0]->FindObject("stats");
st->SetLineColor(kGreen);
st->SetTextColor(kGreen);
y = st->GetY1NDC();
dy = st->GetY2NDC() - y;
st->SetX1NDC(0.72);
st->SetY1NDC(y + dy/2);
st = (TPaveStats *) h[6*i+j][1]->FindObject("stats");
st->SetLineColor(kRed);
st->SetTextColor(kRed);
st->SetX1NDC(0.72);
st->SetY2NDC(y + dy/2);
}
if (h[6*i+j][0]->GetMaximum() > h[6*i+j][1]->GetMaximum()) {
h[6*i+j][0]->Draw();
h[6*i+j][1]->Draw("sames");
} else {
h[6*i+j][1]->Draw();
h[6*i+j][0]->Draw("sames");
}
gPad->Update();
}
cv[i]->Update();
}
fRoot->cd();
for (i=0; i<12; i++) for (j=0; j<2; j++) h[i][j]->Write();
fRoot->Close();
}
void initial_z_diff()
{
//=========Macro generated from canvas: c/
//========= (Fri Jul 24 10:47:18 2015) by ROOT version6.04/00
TCanvas *c = new TCanvas("c", "",0,45,600,500);
c->SetHighLightColor(2);
c->Range(-187.5,-3568.557,187.5,32117.01);
c->SetFillColor(0);
c->SetBorderMode(0);
c->SetBorderSize(2);
c->SetFrameBorderMode(0);
c->SetFrameBorderMode(0);
TH1D *hist__1 = new TH1D("hist__1","Difference in PMT-constructed position and mc-truth start position",100,-150,150);
hist__1->SetBinContent(0,49);
hist__1->SetBinContent(1,12);
hist__1->SetBinContent(2,7);
hist__1->SetBinContent(3,9);
hist__1->SetBinContent(4,12);
hist__1->SetBinContent(5,21);
hist__1->SetBinContent(6,30);
hist__1->SetBinContent(7,36);
hist__1->SetBinContent(8,47);
hist__1->SetBinContent(9,67);
hist__1->SetBinContent(10,70);
hist__1->SetBinContent(11,90);
hist__1->SetBinContent(12,140);
hist__1->SetBinContent(13,162);
hist__1->SetBinContent(14,208);
hist__1->SetBinContent(15,208);
hist__1->SetBinContent(16,317);
hist__1->SetBinContent(17,367);
hist__1->SetBinContent(18,460);
hist__1->SetBinContent(19,581);
hist__1->SetBinContent(20,685);
hist__1->SetBinContent(21,838);
hist__1->SetBinContent(22,963);
hist__1->SetBinContent(23,1079);
hist__1->SetBinContent(24,1231);
hist__1->SetBinContent(25,1281);
hist__1->SetBinContent(26,1264);
hist__1->SetBinContent(27,1318);
hist__1->SetBinContent(28,1318);
hist__1->SetBinContent(29,1265);
hist__1->SetBinContent(30,1282);
hist__1->SetBinContent(31,1235);
hist__1->SetBinContent(32,1315);
hist__1->SetBinContent(33,1290);
hist__1->SetBinContent(34,1399);
hist__1->SetBinContent(35,1363);
hist__1->SetBinContent(36,1528);
hist__1->SetBinContent(37,1700);
hist__1->SetBinContent(38,1822);
hist__1->SetBinContent(39,2053);
hist__1->SetBinContent(40,2441);
hist__1->SetBinContent(41,2903);
hist__1->SetBinContent(42,3359);
hist__1->SetBinContent(43,4326);
hist__1->SetBinContent(44,5618);
hist__1->SetBinContent(45,7284);
hist__1->SetBinContent(46,9575);
hist__1->SetBinContent(47,12857);
hist__1->SetBinContent(48,16804);
hist__1->SetBinContent(49,21562);
hist__1->SetBinContent(50,27189);
hist__1->SetBinContent(51,26947);
hist__1->SetBinContent(52,21828);
hist__1->SetBinContent(53,17254);
hist__1->SetBinContent(54,13388);
hist__1->SetBinContent(55,10171);
hist__1->SetBinContent(56,7452);
hist__1->SetBinContent(57,5752);
hist__1->SetBinContent(58,4452);
hist__1->SetBinContent(59,3522);
hist__1->SetBinContent(60,2826);
hist__1->SetBinContent(61,2449);
hist__1->SetBinContent(62,2170);
hist__1->SetBinContent(63,1855);
hist__1->SetBinContent(64,1635);
hist__1->SetBinContent(65,1515);
hist__1->SetBinContent(66,1393);
hist__1->SetBinContent(67,1344);
hist__1->SetBinContent(68,1279);
hist__1->SetBinContent(69,1232);
hist__1->SetBinContent(70,1248);
hist__1->SetBinContent(71,1277);
hist__1->SetBinContent(72,1226);
hist__1->SetBinContent(73,1247);
hist__1->SetBinContent(74,1229);
hist__1->SetBinContent(75,1187);
hist__1->SetBinContent(76,1265);
hist__1->SetBinContent(77,1253);
hist__1->SetBinContent(78,1269);
hist__1->SetBinContent(79,1210);
hist__1->SetBinContent(80,1055);
hist__1->SetBinContent(81,914);
hist__1->SetBinContent(82,761);
hist__1->SetBinContent(83,632);
hist__1->SetBinContent(84,515);
hist__1->SetBinContent(85,450);
hist__1->SetBinContent(86,356);
hist__1->SetBinContent(87,264);
hist__1->SetBinContent(88,227);
hist__1->SetBinContent(89,176);
hist__1->SetBinContent(90,127);
hist__1->SetBinContent(91,134);
hist__1->SetBinContent(92,85);
hist__1->SetBinContent(93,60);
hist__1->SetBinContent(94,47);
hist__1->SetBinContent(95,29);
hist__1->SetBinContent(96,31);
hist__1->SetBinContent(97,18);
hist__1->SetBinContent(98,14);
hist__1->SetBinContent(99,12);
hist__1->SetBinContent(100,12);
hist__1->SetBinContent(101,63);
hist__1->SetEntries(289897);
TPaveStats *ptstats = new TPaveStats(0.78,0.775,0.98,0.935,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("hist");
AText->SetTextSize(0.0368);
AText = ptstats->AddText("Entries = 289897 ");
AText = ptstats->AddText("Mean = 0.6888");
AText = ptstats->AddText("Std Dev = 33.12");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
ptstats->Draw();
hist__1->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(hist__1);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
hist__1->SetLineColor(ci);
hist__1->GetXaxis()->SetTitle("#DeltaZ [cm]");
hist__1->GetXaxis()->SetLabelFont(42);
hist__1->GetXaxis()->SetLabelSize(0.035);
hist__1->GetXaxis()->SetTitleSize(0.035);
hist__1->GetXaxis()->SetTitleFont(42);
hist__1->GetYaxis()->SetLabelFont(42);
hist__1->GetYaxis()->SetLabelSize(0.035);
hist__1->GetYaxis()->SetTitleSize(0.035);
hist__1->GetYaxis()->SetTitleFont(42);
hist__1->GetZaxis()->SetLabelFont(42);
hist__1->GetZaxis()->SetLabelSize(0.035);
hist__1->GetZaxis()->SetTitleSize(0.035);
hist__1->GetZaxis()->SetTitleFont(42);
hist__1->Draw("");
TPaveText *pt = new TPaveText(0.15,0.9368947,0.85,0.995,"blNDC");
pt->SetName("title");
pt->SetBorderSize(0);
pt->SetFillColor(0);
pt->SetFillStyle(0);
pt->SetTextFont(42);
AText = pt->AddText("Difference in PMT-constructed position and mc-truth start position");
pt->Draw();
c->Modified();
c->cd();
c->SetSelected(c);
}
//========= Macro generated from object: h1/l1l2Pt {(((goodVtx && lepPt[0]>25 && lepPt[1]>20 && l1l2DR>0.3 && !(abs(lepEta[0]) > 1.4 && abs(lepEta[0]) < 1.6) && !(abs(lepEta[1]) > 1.4 && abs(lepEta[1]) < 1.6))*((lepID[0]*lepID[1]==-11*11)||(lepID[0]*lepID[1]==-13*13)))*(totEveW))*(l1l2M>81 && l1l2M < 101)}
//========= by ROOT version6.02/05
void initPtSpectrum()
{
TH1F *h11 = new TH1F("h11","l1l2Pt {(((goodVtx && lepPt[0]>25 && lepPt[1]>20 && l1l2DR>0.3 && !(abs(lepEta[0]) > 1.4 && abs(lepEta[0]) < 1.6) && !(abs(lepEta[1]) > 1.4 && abs(lepEta[1]) < 1.6))*((lepID[0]*lepID[1]==-11*11)||(lepID[0]*lepID[1]==-13*13)))*(totEveW))*(l1l2M>81 && l1l2M < 101)}",500,0,500);
h11->SetBinContent(1,0.00162118);
h11->SetBinContent(2,0.00445683);
h11->SetBinContent(3,0.00644035);
h11->SetBinContent(4,0.00748253);
h11->SetBinContent(5,0.00769155);
h11->SetBinContent(6,0.00751327);
h11->SetBinContent(7,0.00721233);
h11->SetBinContent(8,0.00682488);
h11->SetBinContent(9,0.00637513);
h11->SetBinContent(10,0.00593842);
h11->SetBinContent(11,0.00554317);
h11->SetBinContent(12,0.0052028);
h11->SetBinContent(13,0.00483506);
h11->SetBinContent(14,0.00455992);
h11->SetBinContent(15,0.00425223);
h11->SetBinContent(16,0.00398354);
h11->SetBinContent(17,0.00370052);
h11->SetBinContent(18,0.00348872);
h11->SetBinContent(19,0.00326733);
h11->SetBinContent(20,0.00307653);
h11->SetBinContent(21,0.00290754);
h11->SetBinContent(22,0.0027055);
h11->SetBinContent(23,0.00256185);
h11->SetBinContent(24,0.0024293);
h11->SetBinContent(25,0.00226721);
h11->SetBinContent(26,0.00216945);
h11->SetBinContent(27,0.00203008);
h11->SetBinContent(28,0.00190211);
h11->SetBinContent(29,0.00176505);
h11->SetBinContent(30,0.00170501);
h11->SetBinContent(31,0.00158584);
h11->SetBinContent(32,0.00152027);
h11->SetBinContent(33,0.00145021);
h11->SetBinContent(34,0.00136008);
h11->SetBinContent(35,0.0012743);
h11->SetBinContent(36,0.00123754);
h11->SetBinContent(37,0.0011335);
h11->SetBinContent(38,0.00108776);
h11->SetBinContent(39,0.00105849);
h11->SetBinContent(40,0.000988127);
h11->SetBinContent(41,0.000939292);
h11->SetBinContent(42,0.000887452);
h11->SetBinContent(43,0.000858611);
h11->SetBinContent(44,0.000834414);
h11->SetBinContent(45,0.000771037);
h11->SetBinContent(46,0.000734779);
h11->SetBinContent(47,0.000702866);
h11->SetBinContent(48,0.000678294);
h11->SetBinContent(49,0.00064758);
h11->SetBinContent(50,0.000604205);
h11->SetBinContent(51,0.000591994);
h11->SetBinContent(52,0.000573715);
h11->SetBinContent(53,0.000530116);
h11->SetBinContent(54,0.000523973);
h11->SetBinContent(55,0.000498652);
h11->SetBinContent(56,0.000484644);
h11->SetBinContent(57,0.00044876);
h11->SetBinContent(58,0.000427635);
h11->SetBinContent(59,0.000415274);
h11->SetBinContent(60,0.000411603);
h11->SetBinContent(61,0.000384784);
h11->SetBinContent(62,0.000367629);
h11->SetBinContent(63,0.000363284);
h11->SetBinContent(64,0.000343207);
h11->SetBinContent(65,0.000344181);
h11->SetBinContent(66,0.000322981);
h11->SetBinContent(67,0.00029736);
h11->SetBinContent(68,0.000294738);
h11->SetBinContent(69,0.00027646);
h11->SetBinContent(70,0.000265372);
h11->SetBinContent(71,0.000258181);
h11->SetBinContent(72,0.000230016);
h11->SetBinContent(73,0.000241101);
h11->SetBinContent(74,0.000233312);
h11->SetBinContent(75,0.000217358);
h11->SetBinContent(76,0.000211366);
h11->SetBinContent(77,0.00020897);
h11->SetBinContent(78,0.000200881);
h11->SetBinContent(79,0.000188897);
h11->SetBinContent(80,0.000189721);
h11->SetBinContent(81,0.000178261);
h11->SetBinContent(82,0.000170397);
h11->SetBinContent(83,0.000165828);
h11->SetBinContent(84,0.000162308);
h11->SetBinContent(85,0.000154968);
h11->SetBinContent(86,0.000145156);
h11->SetBinContent(87,0.000131225);
h11->SetBinContent(88,0.000136842);
h11->SetBinContent(89,0.000126431);
h11->SetBinContent(90,0.00012703);
h11->SetBinContent(91,0.000124709);
h11->SetBinContent(92,0.000120364);
h11->SetBinContent(93,0.000109804);
h11->SetBinContent(94,0.000113024);
h11->SetBinContent(95,0.000106883);
h11->SetBinContent(96,0.00010486);
h11->SetBinContent(97,0.000106958);
h11->SetBinContent(98,9.25021e-05);
h11->SetBinContent(99,9.61722e-05);
h11->SetBinContent(100,9.26519e-05);
h11->SetBinContent(101,8.41135e-05);
h11->SetBinContent(102,8.55366e-05);
h11->SetBinContent(103,8.95062e-05);
h11->SetBinContent(104,7.33281e-05);
h11->SetBinContent(105,7.96196e-05);
h11->SetBinContent(106,7.66985e-05);
h11->SetBinContent(107,7.43018e-05);
h11->SetBinContent(108,7.28038e-05);
h11->SetBinContent(109,7.24293e-05);
h11->SetBinContent(110,6.85346e-05);
h11->SetBinContent(111,6.20933e-05);
h11->SetBinContent(112,6.46398e-05);
h11->SetBinContent(113,5.30301e-05);
h11->SetBinContent(114,6.11196e-05);
h11->SetBinContent(115,5.72996e-05);
h11->SetBinContent(116,5.75992e-05);
h11->SetBinContent(117,5.78988e-05);
h11->SetBinContent(118,5.19066e-05);
h11->SetBinContent(119,4.84611e-05);
h11->SetBinContent(120,4.99591e-05);
h11->SetBinContent(121,5.38541e-05);
h11->SetBinContent(122,5.26556e-05);
h11->SetBinContent(123,4.56897e-05);
h11->SetBinContent(124,4.65885e-05);
h11->SetBinContent(125,4.3892e-05);
h11->SetBinContent(126,4.1645e-05);
h11->SetBinContent(127,4.25438e-05);
h11->SetBinContent(128,4.05214e-05);
h11->SetBinContent(129,3.79748e-05);
h11->SetBinContent(130,3.99222e-05);
h11->SetBinContent(131,3.84242e-05);
h11->SetBinContent(132,2.92861e-05);
h11->SetBinContent(133,3.36304e-05);
h11->SetBinContent(134,3.07092e-05);
h11->SetBinContent(135,3.66265e-05);
h11->SetBinContent(136,2.88367e-05);
h11->SetBinContent(137,2.96606e-05);
h11->SetBinContent(138,3.0859e-05);
h11->SetBinContent(139,3.13084e-05);
h11->SetBinContent(140,3.34057e-05);
h11->SetBinContent(141,2.92861e-05);
h11->SetBinContent(142,2.62901e-05);
h11->SetBinContent(143,2.90614e-05);
h11->SetBinContent(144,2.85371e-05);
h11->SetBinContent(145,2.8612e-05);
h11->SetBinContent(146,2.21706e-05);
h11->SetBinContent(147,2.65897e-05);
h11->SetBinContent(148,2.40431e-05);
h11->SetBinContent(149,2.96606e-05);
h11->SetBinContent(150,2.47921e-05);
h11->SetBinContent(151,2.32192e-05);
h11->SetBinContent(152,2.32941e-05);
h11->SetBinContent(153,1.91745e-05);
h11->SetBinContent(154,2.20208e-05);
h11->SetBinContent(155,2.08224e-05);
h11->SetBinContent(156,2.14965e-05);
h11->SetBinContent(157,1.85004e-05);
h11->SetBinContent(158,1.82757e-05);
h11->SetBinContent(159,2.04479e-05);
h11->SetBinContent(160,1.75267e-05);
h11->SetBinContent(161,1.82757e-05);
h11->SetBinContent(162,1.45307e-05);
h11->SetBinContent(163,1.8051e-05);
h11->SetBinContent(164,1.43809e-05);
h11->SetBinContent(165,1.47554e-05);
h11->SetBinContent(166,1.84255e-05);
h11->SetBinContent(167,1.61785e-05);
h11->SetBinContent(168,1.76765e-05);
h11->SetBinContent(169,1.52048e-05);
h11->SetBinContent(170,1.34072e-05);
h11->SetBinContent(171,1.49801e-05);
h11->SetBinContent(172,1.21339e-05);
h11->SetBinContent(173,1.33323e-05);
h11->SetBinContent(174,1.61036e-05);
h11->SetBinContent(175,1.27331e-05);
h11->SetBinContent(176,1.40813e-05);
h11->SetBinContent(177,1.36319e-05);
h11->SetBinContent(178,1.26582e-05);
h11->SetBinContent(179,1.03363e-05);
h11->SetBinContent(180,1.2059e-05);
h11->SetBinContent(181,1.19092e-05);
h11->SetBinContent(182,1.11602e-05);
h11->SetBinContent(183,1.00367e-05);
h11->SetBinContent(184,1.09355e-05);
h11->SetBinContent(185,8.68847e-06);
h11->SetBinContent(186,1.03363e-05);
h11->SetBinContent(187,9.88688e-06);
h11->SetBinContent(188,1.01116e-05);
h11->SetBinContent(189,8.68847e-06);
h11->SetBinContent(190,1.19092e-05);
h11->SetBinContent(191,9.28767e-06);
h11->SetBinContent(192,1.06359e-05);
h11->SetBinContent(193,9.81198e-06);
h11->SetBinContent(194,8.23907e-06);
h11->SetBinContent(195,7.04066e-06);
h11->SetBinContent(196,9.36257e-06);
h11->SetBinContent(197,8.23907e-06);
h11->SetBinContent(198,8.01437e-06);
h11->SetBinContent(199,7.86456e-06);
h11->SetBinContent(200,8.68847e-06);
h11->SetBinContent(201,9.13787e-06);
h11->SetBinContent(202,6.06695e-06);
h11->SetBinContent(203,5.61755e-06);
h11->SetBinContent(204,8.08927e-06);
h11->SetBinContent(205,8.38887e-06);
h11->SetBinContent(206,8.76337e-06);
h11->SetBinContent(207,6.66616e-06);
h11->SetBinContent(208,6.59125e-06);
h11->SetBinContent(209,4.86854e-06);
h11->SetBinContent(210,5.69245e-06);
h11->SetBinContent(211,8.23907e-06);
h11->SetBinContent(212,7.63986e-06);
h11->SetBinContent(213,5.31795e-06);
h11->SetBinContent(214,6.06695e-06);
h11->SetBinContent(215,6.36655e-06);
h11->SetBinContent(216,4.49404e-06);
h11->SetBinContent(217,4.64384e-06);
h11->SetBinContent(218,4.79364e-06);
h11->SetBinContent(219,6.59125e-06);
h11->SetBinContent(220,4.11954e-06);
h11->SetBinContent(221,3.14583e-06);
h11->SetBinContent(222,5.09324e-06);
h11->SetBinContent(223,4.19444e-06);
h11->SetBinContent(224,6.29165e-06);
h11->SetBinContent(225,4.26934e-06);
h11->SetBinContent(226,3.59523e-06);
h11->SetBinContent(227,5.09324e-06);
h11->SetBinContent(228,5.09324e-06);
h11->SetBinContent(229,3.59523e-06);
h11->SetBinContent(230,3.59523e-06);
h11->SetBinContent(231,2.99603e-06);
h11->SetBinContent(232,3.96974e-06);
h11->SetBinContent(233,2.84623e-06);
h11->SetBinContent(234,4.04464e-06);
h11->SetBinContent(235,4.94344e-06);
h11->SetBinContent(236,3.44543e-06);
h11->SetBinContent(237,2.09722e-06);
h11->SetBinContent(238,4.19444e-06);
h11->SetBinContent(239,3.52033e-06);
h11->SetBinContent(240,3.67013e-06);
h11->SetBinContent(241,4.41914e-06);
h11->SetBinContent(242,4.19444e-06);
h11->SetBinContent(243,2.62152e-06);
h11->SetBinContent(244,3.67013e-06);
h11->SetBinContent(245,3.96974e-06);
h11->SetBinContent(246,4.49404e-06);
h11->SetBinContent(247,2.84623e-06);
h11->SetBinContent(248,2.77133e-06);
h11->SetBinContent(249,1.64781e-06);
h11->SetBinContent(250,2.92113e-06);
h11->SetBinContent(251,2.99603e-06);
h11->SetBinContent(252,2.69642e-06);
h11->SetBinContent(253,2.69642e-06);
h11->SetBinContent(254,3.29563e-06);
h11->SetBinContent(255,3.52033e-06);
h11->SetBinContent(256,2.09722e-06);
h11->SetBinContent(257,2.47172e-06);
h11->SetBinContent(258,2.09722e-06);
h11->SetBinContent(259,1.34821e-06);
h11->SetBinContent(260,3.59523e-06);
h11->SetBinContent(261,2.99603e-06);
h11->SetBinContent(262,2.24702e-06);
h11->SetBinContent(263,1.87252e-06);
h11->SetBinContent(264,2.32192e-06);
h11->SetBinContent(265,2.69642e-06);
h11->SetBinContent(266,2.02232e-06);
h11->SetBinContent(267,1.49801e-06);
h11->SetBinContent(268,2.69642e-06);
h11->SetBinContent(269,3.59523e-06);
h11->SetBinContent(270,6.74106e-07);
h11->SetBinContent(271,3.81994e-06);
h11->SetBinContent(272,2.62152e-06);
h11->SetBinContent(273,1.34821e-06);
h11->SetBinContent(274,1.12351e-06);
h11->SetBinContent(275,2.99603e-07);
h11->SetBinContent(276,5.99205e-07);
h11->SetBinContent(277,2.84623e-06);
h11->SetBinContent(278,2.02232e-06);
h11->SetBinContent(279,2.02232e-06);
h11->SetBinContent(280,1.57291e-06);
h11->SetBinContent(281,1.57291e-06);
h11->SetBinContent(282,2.24702e-06);
h11->SetBinContent(283,9.73708e-07);
h11->SetBinContent(284,1.94742e-06);
h11->SetBinContent(285,2.09722e-06);
h11->SetBinContent(286,1.79762e-06);
h11->SetBinContent(287,2.09722e-06);
h11->SetBinContent(288,1.57291e-06);
h11->SetBinContent(289,2.24702e-06);
h11->SetBinContent(290,1.87252e-06);
h11->SetBinContent(291,2.39682e-06);
h11->SetBinContent(292,1.57291e-06);
h11->SetBinContent(293,2.62152e-06);
h11->SetBinContent(294,7.49007e-08);
h11->SetBinContent(295,2.99603e-07);
h11->SetBinContent(296,2.17212e-06);
h11->SetBinContent(297,2.32192e-06);
h11->SetBinContent(298,1.04861e-06);
h11->SetBinContent(299,1.94742e-06);
h11->SetBinContent(300,8.23907e-07);
h11->SetBinContent(301,9.73708e-07);
h11->SetBinContent(302,1.34821e-06);
h11->SetBinContent(303,1.42311e-06);
h11->SetBinContent(304,1.04861e-06);
h11->SetBinContent(305,1.12351e-06);
h11->SetBinContent(306,2.24702e-07);
h11->SetBinContent(307,6.74106e-07);
h11->SetBinContent(308,4.49404e-07);
h11->SetBinContent(309,1.57291e-06);
h11->SetBinContent(310,1.49801e-06);
h11->SetBinContent(311,1.27331e-06);
h11->SetBinContent(312,8.23907e-07);
h11->SetBinContent(313,9.73708e-07);
h11->SetBinContent(314,1.49801e-07);
h11->SetBinContent(315,1.87252e-06);
h11->SetBinContent(316,8.23907e-07);
h11->SetBinContent(317,9.73708e-07);
h11->SetBinContent(318,1.87252e-06);
h11->SetBinContent(319,1.19841e-06);
h11->SetBinContent(320,1.34821e-06);
h11->SetBinContent(321,1.27331e-06);
h11->SetBinContent(322,6.74106e-07);
h11->SetBinContent(323,7.49006e-07);
h11->SetBinContent(324,1.27331e-06);
h11->SetBinContent(325,9.73708e-07);
h11->SetBinContent(327,1.34821e-06);
h11->SetBinContent(328,1.12351e-06);
h11->SetBinContent(329,1.12351e-06);
h11->SetBinContent(330,7.49006e-07);
h11->SetBinContent(331,7.49006e-07);
h11->SetBinContent(332,1.64781e-06);
h11->SetBinContent(333,9.73708e-07);
h11->SetBinContent(334,7.49007e-07);
h11->SetBinContent(335,8.98808e-07);
h11->SetBinContent(336,1.04861e-06);
h11->SetBinContent(337,7.49006e-07);
h11->SetBinContent(338,6.74106e-07);
h11->SetBinContent(339,7.49006e-07);
h11->SetBinContent(340,1.49801e-07);
h11->SetBinContent(341,1.19841e-06);
h11->SetBinContent(342,1.19841e-06);
h11->SetBinContent(343,7.49006e-07);
h11->SetBinContent(344,4.49404e-07);
h11->SetBinContent(345,7.49007e-08);
h11->SetBinContent(346,8.98808e-07);
h11->SetBinContent(347,5.99205e-07);
h11->SetBinContent(348,-2.24702e-07);
h11->SetBinContent(349,8.23907e-07);
h11->SetBinContent(350,5.99205e-07);
h11->SetBinContent(351,6.74106e-07);
h11->SetBinContent(352,5.24305e-07);
h11->SetBinContent(353,3.74503e-07);
h11->SetBinContent(354,2.24702e-07);
h11->SetBinContent(355,9.73708e-07);
h11->SetBinContent(356,3.74503e-07);
h11->SetBinContent(357,7.49006e-07);
h11->SetBinContent(358,4.49404e-07);
h11->SetBinContent(359,2.24702e-07);
h11->SetBinContent(360,8.98808e-07);
h11->SetBinContent(361,2.99603e-07);
h11->SetBinContent(362,4.49404e-07);
h11->SetBinContent(363,5.24305e-07);
h11->SetBinContent(364,5.24305e-07);
h11->SetBinContent(365,8.23907e-07);
h11->SetBinContent(366,7.49007e-07);
h11->SetBinContent(367,9.73708e-07);
h11->SetBinContent(368,1.49801e-07);
h11->SetBinContent(369,5.99205e-07);
h11->SetBinContent(370,5.99205e-07);
h11->SetBinContent(371,3.74503e-07);
h11->SetBinContent(373,3.74503e-07);
h11->SetBinContent(374,5.99205e-07);
h11->SetBinContent(375,7.49006e-07);
h11->SetBinContent(376,7.49007e-08);
h11->SetBinContent(377,8.98808e-07);
h11->SetBinContent(378,3.74503e-07);
h11->SetBinContent(379,6.74106e-07);
h11->SetBinContent(380,3.74503e-07);
h11->SetBinContent(381,7.49006e-07);
h11->SetBinContent(382,-1.49801e-07);
h11->SetBinContent(383,1.49801e-07);
h11->SetBinContent(384,3.74503e-07);
h11->SetBinContent(385,8.23907e-07);
h11->SetBinContent(386,8.23907e-07);
h11->SetBinContent(387,7.49007e-08);
h11->SetBinContent(388,3.74503e-07);
h11->SetBinContent(389,4.49404e-07);
h11->SetBinContent(390,5.99205e-07);
h11->SetBinContent(391,8.23907e-07);
h11->SetBinContent(393,3.74503e-07);
h11->SetBinContent(394,5.99205e-07);
h11->SetBinContent(395,2.99603e-07);
h11->SetBinContent(396,5.24305e-07);
h11->SetBinContent(397,1.49801e-07);
h11->SetBinContent(398,7.49006e-07);
h11->SetBinContent(399,5.99205e-07);
h11->SetBinContent(401,1.49801e-07);
h11->SetBinContent(402,6.74106e-07);
h11->SetBinContent(404,4.49404e-07);
h11->SetBinContent(405,2.24702e-07);
h11->SetBinContent(406,8.23907e-07);
h11->SetBinContent(407,1.49801e-07);
h11->SetBinContent(408,2.24702e-07);
h11->SetBinContent(409,6.74106e-07);
h11->SetBinContent(411,7.49006e-07);
h11->SetBinContent(412,3.74503e-07);
h11->SetBinContent(413,4.49404e-07);
h11->SetBinContent(414,-2.24702e-07);
h11->SetBinContent(415,2.24702e-07);
h11->SetBinContent(416,4.49404e-07);
h11->SetBinContent(417,5.24305e-07);
h11->SetBinContent(418,2.99603e-07);
h11->SetBinContent(419,2.99603e-07);
h11->SetBinContent(420,2.99603e-07);
h11->SetBinContent(421,3.74503e-07);
h11->SetBinContent(422,3.74503e-07);
h11->SetBinContent(423,4.49404e-07);
h11->SetBinContent(424,1.49801e-07);
h11->SetBinContent(425,7.49007e-08);
h11->SetBinContent(426,7.49007e-08);
h11->SetBinContent(427,3.74503e-07);
h11->SetBinContent(428,6.74106e-07);
h11->SetBinContent(429,5.24305e-07);
h11->SetBinContent(430,-1.49801e-07);
h11->SetBinContent(431,5.99205e-07);
h11->SetBinContent(432,3.74503e-07);
h11->SetBinContent(433,7.49007e-08);
h11->SetBinContent(434,4.49404e-07);
h11->SetBinContent(435,2.24702e-07);
h11->SetBinContent(437,-7.49007e-08);
h11->SetBinContent(438,7.49007e-08);
h11->SetBinContent(439,3.74503e-07);
h11->SetBinContent(440,2.24702e-07);
h11->SetBinContent(441,1.49801e-07);
h11->SetBinContent(442,2.99603e-07);
h11->SetBinContent(444,4.49404e-07);
h11->SetBinContent(445,2.24702e-07);
h11->SetBinContent(446,2.99603e-07);
h11->SetBinContent(447,2.99603e-07);
h11->SetBinContent(448,3.74503e-07);
h11->SetBinContent(449,6.74106e-07);
h11->SetBinContent(450,1.49801e-07);
h11->SetBinContent(451,1.49801e-07);
h11->SetBinContent(452,7.49007e-08);
h11->SetBinContent(453,2.24702e-07);
h11->SetBinContent(454,2.24702e-07);
h11->SetBinContent(455,7.49007e-08);
h11->SetBinContent(456,2.24702e-07);
h11->SetBinContent(458,6.74106e-07);
h11->SetBinContent(459,2.24702e-07);
h11->SetBinContent(460,7.49007e-08);
h11->SetBinContent(461,-1.49801e-07);
h11->SetBinContent(462,-2.99603e-07);
h11->SetBinContent(463,1.49801e-07);
h11->SetBinContent(465,-7.49007e-08);
h11->SetBinContent(466,2.24702e-07);
h11->SetBinContent(467,2.99603e-07);
h11->SetBinContent(468,7.49007e-08);
h11->SetBinContent(469,2.24702e-07);
h11->SetBinContent(470,2.99603e-07);
h11->SetBinContent(472,2.99603e-07);
h11->SetBinContent(473,1.49801e-07);
h11->SetBinContent(475,2.99603e-07);
h11->SetBinContent(476,2.24702e-07);
h11->SetBinContent(477,3.74503e-07);
h11->SetBinContent(478,2.99603e-07);
h11->SetBinContent(479,1.49801e-07);
h11->SetBinContent(481,2.99603e-07);
h11->SetBinContent(482,1.49801e-07);
h11->SetBinContent(483,-7.49007e-08);
h11->SetBinContent(484,7.49007e-08);
h11->SetBinContent(486,-7.49007e-08);
h11->SetBinContent(487,1.49801e-07);
h11->SetBinContent(488,2.99603e-07);
h11->SetBinContent(489,7.49007e-08);
h11->SetBinContent(490,7.49007e-08);
h11->SetBinContent(491,4.49404e-07);
h11->SetBinContent(493,2.99603e-07);
h11->SetBinContent(494,7.49007e-08);
h11->SetBinContent(495,1.49801e-07);
h11->SetBinContent(496,3.74503e-07);
h11->SetBinContent(497,1.49801e-07);
h11->SetBinContent(498,7.49007e-08);
h11->SetBinContent(499,7.49007e-08);
h11->SetBinContent(500,1.49801e-07);
h11->SetBinContent(501,5.31795e-06);
h11->SetBinError(1,1.16773e-05);
h11->SetBinError(2,1.94462e-05);
h11->SetBinError(3,2.34236e-05);
h11->SetBinError(4,2.52918e-05);
h11->SetBinError(5,2.58022e-05);
h11->SetBinError(6,2.56526e-05);
h11->SetBinError(7,2.52681e-05);
h11->SetBinError(8,2.46734e-05);
h11->SetBinError(9,2.3991e-05);
h11->SetBinError(10,2.33053e-05);
h11->SetBinError(11,2.26188e-05);
h11->SetBinError(12,2.19792e-05);
h11->SetBinError(13,2.12817e-05);
h11->SetBinError(14,2.07621e-05);
h11->SetBinError(15,2.01215e-05);
h11->SetBinError(16,1.9568e-05);
h11->SetBinError(17,1.90165e-05);
h11->SetBinError(18,1.85359e-05);
h11->SetBinError(19,1.80705e-05);
h11->SetBinError(20,1.76214e-05);
h11->SetBinError(21,1.72475e-05);
h11->SetBinError(22,1.68428e-05);
h11->SetBinError(23,1.64659e-05);
h11->SetBinError(24,1.61463e-05);
h11->SetBinError(25,1.57737e-05);
h11->SetBinError(26,1.55361e-05);
h11->SetBinError(27,1.52325e-05);
h11->SetBinError(28,1.49379e-05);
h11->SetBinError(29,1.46212e-05);
h11->SetBinError(30,1.44319e-05);
h11->SetBinError(31,1.42007e-05);
h11->SetBinError(32,1.40534e-05);
h11->SetBinError(33,1.37491e-05);
h11->SetBinError(34,1.34838e-05);
h11->SetBinError(35,1.32623e-05);
h11->SetBinError(36,1.30295e-05);
h11->SetBinError(37,1.27957e-05);
h11->SetBinError(38,1.26251e-05);
h11->SetBinError(39,1.24504e-05);
h11->SetBinError(40,1.2105e-05);
h11->SetBinError(41,1.19487e-05);
h11->SetBinError(42,1.17809e-05);
h11->SetBinError(43,1.15279e-05);
h11->SetBinError(44,1.1309e-05);
h11->SetBinError(45,1.11239e-05);
h11->SetBinError(46,1.0813e-05);
h11->SetBinError(47,1.07087e-05);
h11->SetBinError(48,1.04981e-05);
h11->SetBinError(49,1.02417e-05);
h11->SetBinError(50,1.01191e-05);
h11->SetBinError(51,9.86729e-06);
h11->SetBinError(52,9.75234e-06);
h11->SetBinError(53,9.48876e-06);
h11->SetBinError(54,9.28496e-06);
h11->SetBinError(55,9.14739e-06);
h11->SetBinError(56,8.98558e-06);
h11->SetBinError(57,8.81062e-06);
h11->SetBinError(58,8.63699e-06);
h11->SetBinError(59,8.50972e-06);
h11->SetBinError(60,8.38252e-06);
h11->SetBinError(61,8.18545e-06);
h11->SetBinError(62,8.09204e-06);
h11->SetBinError(63,7.92108e-06);
h11->SetBinError(64,7.70203e-06);
h11->SetBinError(65,7.66882e-06);
h11->SetBinError(66,7.45177e-06);
h11->SetBinError(67,7.3242e-06);
h11->SetBinError(68,7.26845e-06);
h11->SetBinError(69,7.04265e-06);
h11->SetBinError(70,6.90022e-06);
h11->SetBinError(71,6.79453e-06);
h11->SetBinError(72,6.66111e-06);
h11->SetBinError(73,6.60529e-06);
h11->SetBinError(74,6.43056e-06);
h11->SetBinError(75,6.33697e-06);
h11->SetBinError(76,6.20998e-06);
h11->SetBinError(77,6.10796e-06);
h11->SetBinError(78,5.98643e-06);
h11->SetBinError(79,5.89768e-06);
h11->SetBinError(80,5.83312e-06);
h11->SetBinError(81,5.75225e-06);
h11->SetBinError(82,5.57041e-06);
h11->SetBinError(83,5.57192e-06);
h11->SetBinError(84,5.40584e-06);
h11->SetBinError(85,5.32638e-06);
h11->SetBinError(86,5.24625e-06);
h11->SetBinError(87,5.08994e-06);
h11->SetBinError(88,5.04065e-06);
h11->SetBinError(89,4.96496e-06);
h11->SetBinError(90,4.88063e-06);
h11->SetBinError(91,4.74719e-06);
h11->SetBinError(92,4.76253e-06);
h11->SetBinError(93,4.66554e-06);
h11->SetBinError(94,4.62629e-06);
h11->SetBinError(95,4.47465e-06);
h11->SetBinError(96,4.44004e-06);
h11->SetBinError(97,4.39304e-06);
h11->SetBinError(98,4.32028e-06);
h11->SetBinError(99,4.16356e-06);
h11->SetBinError(100,4.17096e-06);
h11->SetBinError(101,4.03422e-06);
h11->SetBinError(102,4.01121e-06);
h11->SetBinError(103,4.0161e-06);
h11->SetBinError(104,3.88257e-06);
h11->SetBinError(105,3.786e-06);
h11->SetBinError(106,3.73754e-06);
h11->SetBinError(107,3.78823e-06);
h11->SetBinError(108,3.77487e-06);
h11->SetBinError(109,3.68995e-06);
h11->SetBinError(110,3.64251e-06);
h11->SetBinError(111,3.50436e-06);
h11->SetBinError(112,3.54415e-06);
h11->SetBinError(113,3.47784e-06);
h11->SetBinError(114,3.41106e-06);
h11->SetBinError(115,3.42829e-06);
h11->SetBinError(116,3.33875e-06);
h11->SetBinError(117,3.28113e-06);
h11->SetBinError(118,3.24675e-06);
h11->SetBinError(119,3.15028e-06);
h11->SetBinError(120,3.11086e-06);
h11->SetBinError(121,3.11806e-06);
h11->SetBinError(122,3.05628e-06);
h11->SetBinError(123,2.9904e-06);
h11->SetBinError(124,3.02769e-06);
h11->SetBinError(125,2.89702e-06);
h11->SetBinError(126,2.93931e-06);
h11->SetBinError(127,2.84622e-06);
h11->SetBinError(128,2.82347e-06);
h11->SetBinError(129,2.79752e-06);
h11->SetBinError(130,2.84327e-06);
h11->SetBinError(131,2.78546e-06);
h11->SetBinError(132,2.67029e-06);
h11->SetBinError(133,2.62794e-06);
h11->SetBinError(134,2.55761e-06);
h11->SetBinError(135,2.61295e-06);
h11->SetBinError(136,2.53005e-06);
h11->SetBinError(137,2.43629e-06);
h11->SetBinError(138,2.53337e-06);
h11->SetBinError(139,2.46149e-06);
h11->SetBinError(140,2.47286e-06);
h11->SetBinError(141,2.44433e-06);
h11->SetBinError(142,2.32916e-06);
h11->SetBinError(143,2.29885e-06);
h11->SetBinError(144,2.29763e-06);
h11->SetBinError(145,2.28907e-06);
h11->SetBinError(146,2.23701e-06);
h11->SetBinError(147,2.27555e-06);
h11->SetBinError(148,2.23073e-06);
h11->SetBinError(149,2.27925e-06);
h11->SetBinError(150,2.20035e-06);
h11->SetBinError(151,2.10789e-06);
h11->SetBinError(152,2.17212e-06);
h11->SetBinError(153,2.08649e-06);
h11->SetBinError(154,2.04027e-06);
h11->SetBinError(155,2.05124e-06);
h11->SetBinError(156,2.02232e-06);
h11->SetBinError(157,2.0028e-06);
h11->SetBinError(158,2.0098e-06);
h11->SetBinError(159,1.95173e-06);
h11->SetBinError(160,1.94742e-06);
h11->SetBinError(161,1.95891e-06);
h11->SetBinError(162,1.84687e-06);
h11->SetBinError(163,1.85143e-06);
h11->SetBinError(164,1.91547e-06);
h11->SetBinError(165,1.76136e-06);
h11->SetBinError(166,1.87999e-06);
h11->SetBinError(167,1.69149e-06);
h11->SetBinError(168,1.79136e-06);
h11->SetBinError(169,1.80229e-06);
h11->SetBinError(170,1.71292e-06);
h11->SetBinError(171,1.62381e-06);
h11->SetBinError(172,1.64781e-06);
h11->SetBinError(173,1.61688e-06);
h11->SetBinError(174,1.67315e-06);
h11->SetBinError(175,1.63757e-06);
h11->SetBinError(176,1.58888e-06);
h11->SetBinError(177,1.60294e-06);
h11->SetBinError(178,1.53318e-06);
h11->SetBinError(179,1.61341e-06);
h11->SetBinError(180,1.60119e-06);
h11->SetBinError(181,1.54048e-06);
h11->SetBinError(182,1.50735e-06);
h11->SetBinError(183,1.51292e-06);
h11->SetBinError(184,1.51292e-06);
h11->SetBinError(185,1.46392e-06);
h11->SetBinError(186,1.50175e-06);
h11->SetBinError(187,1.41322e-06);
h11->SetBinError(188,1.41521e-06);
h11->SetBinError(189,1.43684e-06);
h11->SetBinError(190,1.46583e-06);
h11->SetBinError(191,1.42901e-06);
h11->SetBinError(192,1.37295e-06);
h11->SetBinError(193,1.42311e-06);
h11->SetBinError(194,1.35651e-06);
h11->SetBinError(195,1.2578e-06);
h11->SetBinError(196,1.38718e-06);
h11->SetBinError(197,1.24884e-06);
h11->SetBinError(198,1.35444e-06);
h11->SetBinError(199,1.29948e-06);
h11->SetBinError(200,1.2578e-06);
h11->SetBinError(201,1.25333e-06);
h11->SetBinError(202,1.21929e-06);
h11->SetBinError(203,1.22846e-06);
h11->SetBinError(204,1.23983e-06);
h11->SetBinError(205,1.15551e-06);
h11->SetBinError(206,1.25109e-06);
h11->SetBinError(207,1.23302e-06);
h11->SetBinError(208,1.14576e-06);
h11->SetBinError(209,1.14331e-06);
h11->SetBinError(210,1.18428e-06);
h11->SetBinError(211,1.12101e-06);
h11->SetBinError(212,1.11096e-06);
h11->SetBinError(213,1.08283e-06);
h11->SetBinError(214,1.14331e-06);
h11->SetBinError(215,1.08283e-06);
h11->SetBinError(216,1.14576e-06);
h11->SetBinError(217,1.016e-06);
h11->SetBinError(218,1.126e-06);
h11->SetBinError(219,1.03243e-06);
h11->SetBinError(220,9.50383e-07);
h11->SetBinError(221,1.03243e-06);
h11->SetBinError(222,1.08541e-06);
h11->SetBinError(223,9.99299e-07);
h11->SetBinError(224,9.88007e-07);
h11->SetBinError(225,9.62116e-07);
h11->SetBinError(226,9.82313e-07);
h11->SetBinError(227,9.70823e-07);
h11->SetBinError(228,9.99299e-07);
h11->SetBinError(229,9.88007e-07);
h11->SetBinError(230,9.29493e-07);
h11->SetBinError(231,9.5333e-07);
h11->SetBinError(232,9.08122e-07);
h11->SetBinError(233,8.53999e-07);
h11->SetBinError(234,9.41487e-07);
h11->SetBinError(235,9.29493e-07);
h11->SetBinError(236,9.17342e-07);
h11->SetBinError(237,8.67038e-07);
h11->SetBinError(238,9.05028e-07);
h11->SetBinError(239,8.50708e-07);
h11->SetBinError(240,8.83066e-07);
h11->SetBinError(241,8.63796e-07);
h11->SetBinError(242,8.40758e-07);
h11->SetBinError(243,8.70267e-07);
h11->SetBinError(244,8.1707e-07);
h11->SetBinError(245,8.23907e-07);
h11->SetBinError(246,8.40758e-07);
h11->SetBinError(247,8.13629e-07);
h11->SetBinError(248,8.50708e-07);
h11->SetBinError(249,8.06705e-07);
h11->SetBinError(250,7.96205e-07);
h11->SetBinError(251,7.92674e-07);
h11->SetBinError(252,7.5646e-07);
h11->SetBinError(253,7.18422e-07);
h11->SetBinError(254,7.49007e-07);
h11->SetBinError(255,7.60159e-07);
h11->SetBinError(256,7.5646e-07);
h11->SetBinError(257,7.96205e-07);
h11->SetBinError(258,7.0263e-07);
h11->SetBinError(259,7.33874e-07);
h11->SetBinError(260,6.52969e-07);
h11->SetBinError(261,7.18422e-07);
h11->SetBinError(262,7.49007e-07);
h11->SetBinError(263,6.98627e-07);
h11->SetBinError(264,6.03868e-07);
h11->SetBinError(265,6.4432e-07);
h11->SetBinError(266,7.14507e-07);
h11->SetBinError(267,7.0263e-07);
h11->SetBinError(268,7.26189e-07);
h11->SetBinError(269,7.18422e-07);
h11->SetBinError(270,6.82377e-07);
h11->SetBinError(271,6.74106e-07);
h11->SetBinError(272,6.65732e-07);
h11->SetBinError(273,6.61505e-07);
h11->SetBinError(274,6.03868e-07);
h11->SetBinError(275,6.35553e-07);
h11->SetBinError(276,6.26664e-07);
h11->SetBinError(277,6.52969e-07);
h11->SetBinError(278,5.75323e-07);
h11->SetBinError(279,6.31124e-07);
h11->SetBinError(280,6.98627e-07);
h11->SetBinError(281,6.13088e-07);
h11->SetBinError(282,6.4432e-07);
h11->SetBinError(283,5.65488e-07);
h11->SetBinError(284,5.50405e-07);
h11->SetBinError(285,6.17647e-07);
h11->SetBinError(286,5.89768e-07);
h11->SetBinError(287,6.08496e-07);
h11->SetBinError(288,6.13088e-07);
h11->SetBinError(289,5.40116e-07);
h11->SetBinError(290,5.55478e-07);
h11->SetBinError(291,5.29628e-07);
h11->SetBinError(292,6.57251e-07);
h11->SetBinError(293,5.75323e-07);
h11->SetBinError(294,5.02449e-07);
h11->SetBinError(295,5.08001e-07);
h11->SetBinError(296,5.65488e-07);
h11->SetBinError(297,5.55478e-07);
h11->SetBinError(298,5.60505e-07);
h11->SetBinError(299,5.40116e-07);
h11->SetBinError(300,5.55478e-07);
h11->SetBinError(301,5.34898e-07);
h11->SetBinError(302,4.85412e-07);
h11->SetBinError(303,6.03868e-07);
h11->SetBinError(304,5.18927e-07);
h11->SetBinError(305,5.02449e-07);
h11->SetBinError(306,5.02449e-07);
h11->SetBinError(307,4.03352e-07);
h11->SetBinError(308,4.96835e-07);
h11->SetBinError(309,5.02449e-07);
h11->SetBinError(310,5.50405e-07);
h11->SetBinError(311,4.79598e-07);
h11->SetBinError(312,4.67754e-07);
h11->SetBinError(313,5.13493e-07);
h11->SetBinError(314,4.36742e-07);
h11->SetBinError(315,4.43118e-07);
h11->SetBinError(316,4.55603e-07);
h11->SetBinError(317,5.13493e-07);
h11->SetBinError(318,5.34898e-07);
h11->SetBinError(319,4.61719e-07);
h11->SetBinError(320,4.49404e-07);
h11->SetBinError(321,3.89195e-07);
h11->SetBinError(322,3.89195e-07);
h11->SetBinError(323,3.8192e-07);
h11->SetBinError(324,4.67754e-07);
h11->SetBinError(325,4.30272e-07);
h11->SetBinError(326,4.10248e-07);
h11->SetBinError(327,5.18927e-07);
h11->SetBinError(328,4.03352e-07);
h11->SetBinError(329,4.17029e-07);
h11->SetBinError(330,3.51315e-07);
h11->SetBinError(331,3.96337e-07);
h11->SetBinError(332,4.61719e-07);
h11->SetBinError(333,4.91156e-07);
h11->SetBinError(334,3.66937e-07);
h11->SetBinError(335,4.23702e-07);
h11->SetBinError(336,4.10248e-07);
h11->SetBinError(337,3.96337e-07);
h11->SetBinError(338,3.59211e-07);
h11->SetBinError(339,3.8192e-07);
h11->SetBinError(340,3.8192e-07);
h11->SetBinError(341,4.49404e-07);
h11->SetBinError(342,3.66937e-07);
h11->SetBinError(343,4.10248e-07);
h11->SetBinError(344,3.17777e-07);
h11->SetBinError(345,3.43238e-07);
h11->SetBinError(346,4.23702e-07);
h11->SetBinError(347,3.66937e-07);
h11->SetBinError(348,4.17029e-07);
h11->SetBinError(349,3.43238e-07);
h11->SetBinError(350,3.96337e-07);
h11->SetBinError(351,3.08823e-07);
h11->SetBinError(352,4.03352e-07);
h11->SetBinError(353,3.26484e-07);
h11->SetBinError(354,3.43238e-07);
h11->SetBinError(355,4.30272e-07);
h11->SetBinError(356,2.70058e-07);
h11->SetBinError(357,3.34966e-07);
h11->SetBinError(358,2.80253e-07);
h11->SetBinError(359,3.59211e-07);
h11->SetBinError(360,3.66937e-07);
h11->SetBinError(361,3.34966e-07);
h11->SetBinError(362,3.34966e-07);
h11->SetBinError(363,3.26484e-07);
h11->SetBinError(364,3.59211e-07);
h11->SetBinError(365,3.59211e-07);
h11->SetBinError(366,3.51315e-07);
h11->SetBinError(367,3.74503e-07);
h11->SetBinError(368,3.17777e-07);
h11->SetBinError(369,3.17777e-07);
h11->SetBinError(370,3.17777e-07);
h11->SetBinError(371,2.90089e-07);
h11->SetBinError(372,2.59463e-07);
h11->SetBinError(373,3.43238e-07);
h11->SetBinError(374,3.17777e-07);
h11->SetBinError(375,2.80253e-07);
h11->SetBinError(376,2.24702e-07);
h11->SetBinError(377,3.51315e-07);
h11->SetBinError(378,2.24702e-07);
h11->SetBinError(379,2.70058e-07);
h11->SetBinError(380,2.90089e-07);
h11->SetBinError(381,3.34966e-07);
h11->SetBinError(382,2.59463e-07);
h11->SetBinError(383,1.49801e-07);
h11->SetBinError(384,3.26484e-07);
h11->SetBinError(385,3.43238e-07);
h11->SetBinError(386,3.26484e-07);
h11->SetBinError(387,2.70058e-07);
h11->SetBinError(388,2.90089e-07);
h11->SetBinError(389,2.99603e-07);
h11->SetBinError(390,2.99603e-07);
h11->SetBinError(391,3.26484e-07);
h11->SetBinError(392,3.17777e-07);
h11->SetBinError(393,2.48417e-07);
h11->SetBinError(394,2.80253e-07);
h11->SetBinError(395,3.17777e-07);
h11->SetBinError(396,3.08823e-07);
h11->SetBinError(397,3.17777e-07);
h11->SetBinError(398,2.80253e-07);
h11->SetBinError(399,2.80253e-07);
h11->SetBinError(400,2.80253e-07);
h11->SetBinError(401,2.99603e-07);
h11->SetBinError(402,3.08823e-07);
h11->SetBinError(403,2.11851e-07);
h11->SetBinError(404,2.59463e-07);
h11->SetBinError(405,1.98169e-07);
h11->SetBinError(406,3.08823e-07);
h11->SetBinError(407,2.11851e-07);
h11->SetBinError(408,2.24702e-07);
h11->SetBinError(409,2.70058e-07);
h11->SetBinError(410,2.11851e-07);
h11->SetBinError(411,2.99603e-07);
h11->SetBinError(412,2.90089e-07);
h11->SetBinError(413,2.59463e-07);
h11->SetBinError(414,2.24702e-07);
h11->SetBinError(415,1.98169e-07);
h11->SetBinError(416,3.34966e-07);
h11->SetBinError(417,2.48417e-07);
h11->SetBinError(418,2.11851e-07);
h11->SetBinError(419,1.83468e-07);
h11->SetBinError(420,2.36857e-07);
h11->SetBinError(421,2.70058e-07);
h11->SetBinError(422,1.98169e-07);
h11->SetBinError(423,2.59463e-07);
h11->SetBinError(424,2.36857e-07);
h11->SetBinError(425,1.67483e-07);
h11->SetBinError(426,2.24702e-07);
h11->SetBinError(427,3.08823e-07);
h11->SetBinError(428,2.24702e-07);
h11->SetBinError(429,2.48417e-07);
h11->SetBinError(430,2.36857e-07);
h11->SetBinError(431,2.59463e-07);
h11->SetBinError(432,2.48417e-07);
h11->SetBinError(433,2.24702e-07);
h11->SetBinError(434,1.83468e-07);
h11->SetBinError(435,1.98169e-07);
h11->SetBinError(437,1.29732e-07);
h11->SetBinError(438,2.24702e-07);
h11->SetBinError(439,2.24702e-07);
h11->SetBinError(440,2.24702e-07);
h11->SetBinError(441,2.36857e-07);
h11->SetBinError(442,2.59463e-07);
h11->SetBinError(443,1.05926e-07);
h11->SetBinError(444,2.36857e-07);
h11->SetBinError(445,1.98169e-07);
h11->SetBinError(446,1.83468e-07);
h11->SetBinError(447,2.36857e-07);
h11->SetBinError(448,2.90089e-07);
h11->SetBinError(449,2.48417e-07);
h11->SetBinError(450,2.11851e-07);
h11->SetBinError(451,1.49801e-07);
h11->SetBinError(452,1.29732e-07);
h11->SetBinError(453,1.67483e-07);
h11->SetBinError(454,1.98169e-07);
h11->SetBinError(455,1.67483e-07);
h11->SetBinError(456,1.67483e-07);
h11->SetBinError(457,1.05926e-07);
h11->SetBinError(458,2.70058e-07);
h11->SetBinError(459,1.98169e-07);
h11->SetBinError(460,1.98169e-07);
h11->SetBinError(461,1.49801e-07);
h11->SetBinError(462,1.83468e-07);
h11->SetBinError(463,1.05926e-07);
h11->SetBinError(464,1.49801e-07);
h11->SetBinError(465,1.67483e-07);
h11->SetBinError(466,1.98169e-07);
h11->SetBinError(467,1.83468e-07);
h11->SetBinError(468,1.29732e-07);
h11->SetBinError(469,1.67483e-07);
h11->SetBinError(470,2.36857e-07);
h11->SetBinError(471,2.11851e-07);
h11->SetBinError(472,2.11851e-07);
h11->SetBinError(473,1.49801e-07);
h11->SetBinError(474,2.11851e-07);
h11->SetBinError(475,1.83468e-07);
h11->SetBinError(476,1.98169e-07);
h11->SetBinError(477,1.67483e-07);
h11->SetBinError(478,1.49801e-07);
h11->SetBinError(479,1.83468e-07);
h11->SetBinError(480,1.83468e-07);
h11->SetBinError(481,1.83468e-07);
h11->SetBinError(482,1.49801e-07);
h11->SetBinError(483,1.98169e-07);
h11->SetBinError(484,1.29732e-07);
h11->SetBinError(485,1.83468e-07);
h11->SetBinError(486,7.49007e-08);
h11->SetBinError(487,1.49801e-07);
h11->SetBinError(488,1.49801e-07);
h11->SetBinError(489,1.67483e-07);
h11->SetBinError(490,1.67483e-07);
h11->SetBinError(491,2.11851e-07);
h11->SetBinError(492,1.83468e-07);
h11->SetBinError(493,1.83468e-07);
h11->SetBinError(494,1.29732e-07);
h11->SetBinError(495,1.49801e-07);
h11->SetBinError(496,1.98169e-07);
h11->SetBinError(497,1.05926e-07);
h11->SetBinError(498,1.67483e-07);
h11->SetBinError(499,1.29732e-07);
h11->SetBinError(500,2.11851e-07);
h11->SetBinError(501,1.38718e-06);
h11->SetEntries(3.17523e+06);
TPaveStats *ptstats = new TPaveStats(0.78,0.775,0.98,0.935,"brNDC");
ptstats->SetName("stats");
ptstats->SetBorderSize(1);
ptstats->SetFillColor(0);
ptstats->SetTextAlign(12);
ptstats->SetTextFont(42);
TText *AText = ptstats->AddText("h1");
AText->SetTextSize(0.0368);
AText = ptstats->AddText("Entries = 3175226");
AText = ptstats->AddText("Mean = 22.48");
AText = ptstats->AddText("RMS = 27.23");
ptstats->SetOptStat(1111);
ptstats->SetOptFit(0);
// ptstats->Draw();
h11->GetListOfFunctions()->Add(ptstats);
ptstats->SetParent(h11);
/*
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#000099");
h11->SetLineColor(ci);
h11->GetXaxis()->SetLabelFont(42);
h11->GetXaxis()->SetLabelSize(0.035);
h11->GetXaxis()->SetTitleSize(0.035);
h11->GetXaxis()->SetTitleFont(42);
h11->GetYaxis()->SetLabelFont(42);
h11->GetYaxis()->SetLabelSize(0.035);
h11->GetYaxis()->SetTitleSize(0.035);
h11->GetYaxis()->SetTitleFont(42);
h11->GetZaxis()->SetLabelFont(42);
h11->GetZaxis()->SetLabelSize(0.035);
h11->GetZaxis()->SetTitleSize(0.035);
h11->GetZaxis()->SetTitleFont(42);
h11->Draw("");
*/
h_ptSpectrum = h11;
}
