Skip to content

Code and combine output #12

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
danielguerrero wants to merge 11 commits into
base: master
Choose a base branch
from
Open

Code and combine output #12

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
8f977c8
New files
danielguerrero Sep 21, 2017
d626aa4
New files
danielguerrero Sep 22, 2017
0c87086
New files
danielguerrero Sep 22, 2017
cd13990
New files
danielguerrero Sep 26, 2017
25ebea1
New files
danielguerrero Oct 3, 2017
91ef8c9
New files
danielguerrero Oct 3, 2017
28ea66c
New errorbars
danielguerrero Oct 26, 2017
5136db8
New errorbars
danielguerrero Oct 26, 2017
b99148d
New errorbars
danielguerrero Oct 26, 2017
42fceb5
Updates
danielguerrero Dec 29, 2017
92e237a
Updates
danielguerrero Dec 29, 2017
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
The table of contents is too big for display.
Diff view
Diff view
  •  
  •  
  •  
606 changes: 263 additions & 343 deletions BackgroundPrediction_Kinematic_Split.c
View file
Open in desktop

Large diffs are not rendered by default.

430 changes: 158 additions & 272 deletions BackgroundPrediction_Kinematic_Split_MMR.c
View file
Open in desktop

Large diffs are not rendered by default.

887 changes: 887 additions & 0 deletions BackgroundPrediction_Kinematic_Split_MMR_blinded.c
View file
Open in desktop

Large diffs are not rendered by default.

974 changes: 974 additions & 0 deletions BackgroundPrediction_Kinematic_Split_blinded.c
View file
Open in desktop

Large diffs are not rendered by default.

382 changes: 382 additions & 0 deletions DisplayCutFlow.cc
View file
Open in desktop

Large diffs are not rendered by default.

61 changes: 25 additions & 36 deletions DisplayCutFlow_LMR.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -208,7 +208,7 @@ void pushBackHistograms(Samples &sample, TFile *file, TFile * old, bool MC=true)

void signalEff(int npoint)
{
int massPtsStd[8]={260, 270, 300, 350, 400, 450, 500, 550};
int massPtsStd[10]={260, 270, 300, 350, 400, 450, 500, 550, 600, 650};
int massPts[npoint];
for(int i=0; i<npoint; i++) massPts[i]=massPtsStd[i];
std::cout<<" === Signal Efficiencies === "<<std::endl;
Expand All @@ -229,7 +229,7 @@ void signalEff(int npoint)

void cumEff(int npoint)
{
int massPtsStd[8]={260, 270, 300, 350, 400, 450, 500, 550};
int massPtsStd[10]={260, 270, 300, 350, 400, 450, 500, 550, 600, 650};
int massPts[npoint];
for(int i=0; i<npoint; i++) massPts[i]=massPtsStd[i];

Expand Down Expand Up @@ -260,40 +260,28 @@ void DisplayCutFlow_LMR()
lumi_13TeV = "(2016)"; // default is "5.1 fb^{-1}"


//int n=10;
int massL[6]={260, 270, 300, 400, 500, 550};// 800, 1000, 1200};
// int massL[5]={260, 270, 300, 400, 600}; //, 700, 800, 900, 1000, 1100};
TFile *glugluToX260=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-260_narrow_13TeV-madgraph.root");
TFile *glugluToX270=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-270_narrow_13TeV-madgraph.root");
TFile *glugluToX300=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-300_narrow_13TeV-madgraph.root");
//TFile *glugluToX350=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-350_narrow_13TeV-madgraph.root");
TFile *glugluToX400=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-400_narrow_13TeV-madgraph.root");
// TFile *glugluToX450=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-450_narrow_13TeV-madgraph.root");
TFile *glugluToX500=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-500_narrow_13TeV-madgraph.root");
TFile *glugluToX550=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-550_narrow_13TeV-madgraph.root");
//TFile *glugluToX600=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-600_narrow_13TeV-madgraph.root");
//TFile *glugluToX650=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-650_narrow_13TeV-madgraph.root");
// TFile *glugluToX700=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-700_narrow_13TeV-madgraph.root");
//TFile *glugluToX800=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-800_narrow_13TeV-madgraph.root");
//TFile *glugluToX900=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-900_narrow_13TeV-madgraph.root");
//TFile *glugluToX1000=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-1000_narrow_13TeV-madgraph.root");
//TFile *glugluToX1200=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-1200_narrow_13TeV-madgraph.root");
int massL[10]={260, 270, 300, 350, 400, 450, 500, 550, 600, 650};
TFile *glugluToX260=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-260_narrow_13TeV-madgraph.root");
TFile *glugluToX270=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-270_narrow_13TeV-madgraph.root");
TFile *glugluToX300=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-300_narrow_13TeV-madgraph.root");
TFile *glugluToX350=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-350_narrow_13TeV-madgraph.root");
TFile *glugluToX400=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-400_narrow_13TeV-madgraph.root");
TFile *glugluToX450=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-450_narrow_13TeV-madgraph.root");
TFile *glugluToX500=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-500_narrow_13TeV-madgraph.root");
TFile *glugluToX550=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-550_narrow_13TeV-madgraph.root");
TFile *glugluToX600=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-600_narrow_13TeV-madgraph.root");
TFile *glugluToX650=new TFile("PreselectedWithRegressionDeepCSV/LMRSelection_chi2/Histograms_LMR_GluGluToBulkGravitonToHHTo4B_M-650_narrow_13TeV-madgraph.root");

pushBackHistograms(signals, glugluToX260, glugluToX260);
pushBackHistograms(signals, glugluToX270, glugluToX270);
pushBackHistograms(signals, glugluToX300, glugluToX300 );
// pushBackHistograms(signals, glugluToX350, glugluToX350 );
pushBackHistograms(signals, glugluToX300, glugluToX300);
pushBackHistograms(signals, glugluToX350, glugluToX350);
pushBackHistograms(signals, glugluToX400, glugluToX400);
// pushBackHistograms(signals, glugluToX450, glugluToX450);
pushBackHistograms(signals, glugluToX450, glugluToX450);
pushBackHistograms(signals, glugluToX500, glugluToX500);
pushBackHistograms(signals, glugluToX550, glugluToX550);
//pushBackHistograms(signals, glugluToX600, glugluToX600);
// pushBackHistograms(signals, glugluToX650, glugluToX650);
// pushBackHistograms(signals, glugluToX700, glugluToX700);
//pushBackHistograms(signals, glugluToX800, glugluToX800);
//pushBackHistograms(signals, glugluToX900, glugluToX900);
//pushBackHistograms(signals, glugluToX1000, glugluToX1000);
//pushBackHistograms(signals, glugluToX1200, glugluToX1200);
pushBackHistograms(signals, glugluToX600, glugluToX600);
pushBackHistograms(signals, glugluToX650, glugluToX650);

std::vector<TH1F *> h_Cut;
TH1F *h_Init=new TH1F("h_Init", "h_Init", 154, 260, 1800);
Expand All @@ -302,13 +290,9 @@ void DisplayCutFlow_LMR()
TH1F *h_Events8=(TH1F*)h_Init->Clone("h_Events8");
TH1F *h_Events10=(TH1F*)h_Init->Clone("h_Events10");
TH1F *h_Events16=(TH1F*)h_Init->Clone("h_Events16");
for (unsigned int i=0; i<6; ++i)
// for (unsigned int i=0; i<9; ++i)
for (unsigned int i=0; i<10; ++i)
{
h_Init->SetBinContent(h_Init->FindBin(massL[i]), signals.v_CountWithPU.at(i)->GetBinContent(1));
// std::cout<<"signals.v_CountWithPU.at(i)->GetBinContent(1) = "<<.signals.v_CountWithPU.at(i)->GetBinContent(1)<<std::endl;
// h_Events2->SetBinContent(i, signals.v_Cuts.at(i)->GetBinContent(2));
//
std::cout<<massL[i]<<" "<<signals.v_Cuts.at(i)->GetBinContent(4)<< " "<<signals.v_Cuts.at(i)->GetBinContent(8)<<" "<< signals.v_Cuts.at(i)->GetBinContent(10)<< " "<< signals.v_Cuts.at(i)->GetBinContent(12)<<std::endl;
h_Events3->SetBinContent(h_Init->FindBin(massL[i]), signals.v_Cuts.at(i)->GetBinContent(4));
h_Events8->SetBinContent(h_Init->FindBin(massL[i]), signals.v_Cuts.at(i)->GetBinContent(8));
Expand Down Expand Up @@ -372,7 +356,12 @@ void DisplayCutFlow_LMR()
// g_Ae_2->Draw("AL*");

c_Ae->SaveAs("c_Ae_LMR.pdf");
c_Ae->SaveAs("c_Ae.root");
c_Ae->SaveAs("c_Ae_LMR.png");
c_Ae->SaveAs("c_Ae_LMR.C");
c_Ae->SaveAs("c_Ae_LMR.root");





}
55 changes: 24 additions & 31 deletions DisplayCutFlow_MMR.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -207,7 +207,7 @@ void pushBackHistograms(Samples &sample, TFile *file, TFile * old, bool MC=true)

void signalEff(int npoint)
{
int massPtsStd[10]={450, 500, 550, 650, 700, 750, 800, 900, 1000, 1200};// 350, 450, 500 ,550, 650, 700, 800, 900, 1000, 1200};
int massPtsStd[11]={400, 450, 500, 550, 600, 650, 750, 800, 900, 1000, 1200};
int massPts[npoint];


Expand All @@ -230,7 +230,7 @@ void signalEff(int npoint)

void cumEff(int npoint)
{
int massPtsStd[10]={450, 500, 550, 650, 700, 750, 800, 900, 1000, 1200};// 350, 450, 550, 650, 700, 800, 900, 1000, 1200};
int massPtsStd[11]={400, 450, 500, 550, 600, 650, 750, 800, 900, 1000, 1200};

int massPts[npoint];
for(int i=0; i<npoint; i++) massPts[i]=massPtsStd[i];
Expand Down Expand Up @@ -258,41 +258,35 @@ void DisplayCutFlow_MMR()
gROOT->SetStyle("Plain");
gStyle->SetOptStat(0000);
writeExtraText = true; // if extra text
extraText = "Preliminary"; // default extra text is "Preliminary"
extraText = "Simulation"; // default extra text is "Preliminary"
lumi_13TeV = "(2016)"; // default is "5.1 fb^{-1}"


//int n=10;
int massL[]={400, 500, 550, 600, 650, 750, 800, 900,1000, 1200};
// int massL[5]={260, 270, 300, 400, 600}; //, 700, 800, 900, 1000, 1100};
TFile *glugluToX400=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-400_narrow_13TeV-madgraph.root");
//TFile *glugluToX450=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-450_narrow_13TeV-madgraph.root");
TFile *glugluToX500=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-500_narrow_13TeV-madgraph.root");
TFile *glugluToX550=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-550_narrow_13TeV-madgraph.root");
TFile *glugluToX600=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-600_narrow_13TeV-madgraph.root");
TFile *glugluToX650=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-650_narrow_13TeV-madgraph.root");
//TFile *glugluToX700=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-700_narrow_13TeV-madgraph.root");
TFile *glugluToX750=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-750_narrow_13TeV-madgraph.root");

TFile *glugluToX800=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-800_narrow_13TeV-madgraph.root");
TFile *glugluToX900=new TFile("Histograms_GluGluToBulkGravitonToHHTo4B_M-900_narrow_13TeV-madgraph.root");
TFile *glugluToX1000=new TFile("Histograms_BulkGravTohhTohbbhbb_narrow_M-1000_13TeV-madgraph.root");
TFile *glugluToX1200=new TFile("Histograms_BulkGravTohhTohbbhbb_narrow_M-1200_13TeV-madgraph.root");
// TFile *glugluToX1400=new TFile("Histograms_BulkGravTohhTohbbhbb_narrow_M-1400_13TeV-madgraph.root");

pushBackHistograms(signals, glugluToX400, glugluToX400);
// pushBackHistograms(signals, glugluToX450, glugluToX450);
int massL[]={400, 450, 500, 550, 600, 650, 750, 800, 900, 1000, 1200};

TFile *glugluToX400=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-400_narrow_13TeV-madgraph.root");
TFile *glugluToX450=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-450_narrow_13TeV-madgraph.root");
TFile *glugluToX500=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-500_narrow_13TeV-madgraph.root");
TFile *glugluToX550=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-550_narrow_13TeV-madgraph.root");
TFile *glugluToX600=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-600_narrow_13TeV-madgraph.root");
TFile *glugluToX650=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-650_narrow_13TeV-madgraph.root");
TFile *glugluToX750=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-750_narrow_13TeV-madgraph.root");
TFile *glugluToX800=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-800_narrow_13TeV-madgraph.root");
TFile *glugluToX900=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_GluGluToBulkGravitonToHHTo4B_M-900_narrow_13TeV-madgraph.root");
TFile *glugluToX1000=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_BulkGravTohhTohbbhbb_narrow_M-1000_13TeV-madgraph.root");
TFile *glugluToX1200=new TFile("PreselectedWithRegressionDeepCSV/MMRSelection_chi2/Histograms_MMR_BulkGravTohhTohbbhbb_narrow_M-1200_13TeV-madgraph.root");

pushBackHistograms(signals, glugluToX400, glugluToX400);
pushBackHistograms(signals, glugluToX450, glugluToX450);
pushBackHistograms(signals, glugluToX500, glugluToX500);
pushBackHistograms(signals, glugluToX550, glugluToX550);
pushBackHistograms(signals, glugluToX600, glugluToX600);
pushBackHistograms(signals, glugluToX650, glugluToX650);
// pushBackHistograms(signals, glugluToX700, glugluToX700);
pushBackHistograms(signals, glugluToX750, glugluToX750);
pushBackHistograms(signals, glugluToX800, glugluToX800);
pushBackHistograms(signals, glugluToX900, glugluToX900);
pushBackHistograms(signals, glugluToX1000, glugluToX1000);
pushBackHistograms(signals, glugluToX1200, glugluToX1200);
// pushBackHistograms(signals, glugluToX1400, glugluToX1400);

std::vector<TH1F *> h_Cut;
TH1F *h_Init=new TH1F("h_Init", "h_Init", 154, 260, 1800);
Expand All @@ -301,8 +295,8 @@ TH1F *h_Events3=(TH1F*)h_Init->Clone("h_Events3");
TH1F *h_Events8=(TH1F*)h_Init->Clone("h_Events8");
TH1F *h_Events10=(TH1F*)h_Init->Clone("h_Events10");
TH1F *h_Events16=(TH1F*)h_Init->Clone("h_Events16");
// for (unsigned int i=2; i<=8; ++i)
for (unsigned int i=0; i<9; ++i)

for (unsigned int i=0; i<11; ++i)
{
h_Init->SetBinContent(h_Init->FindBin(massL[i]), signals.v_CountWithPU.at(i)->GetBinContent(1));
// std::cout<<"signals.v_CountWithPU.at(i)->GetBinContent(1) = "<<.signals.v_CountWithPU.at(i)->GetBinContent(1)<<std::endl;
Expand Down Expand Up @@ -368,11 +362,10 @@ g_Ae_3->SetMinimum(0.0001);

c_Ae->Update();


// g_Ae_2->Draw("AL*");

c_Ae->SaveAs("c_Ae_MMR.pdf");
c_Ae->SaveAs("c_Ae.root");
c_Ae->SaveAs("c_Ae_MMR.png");
c_Ae->SaveAs("c_Ae_MMR.C");
c_Ae->SaveAs("c_Ae_MMR.root");


}
Loading