void setRBellanStyle() {
TStyle *theStyle = new TStyle("rbStyle", "Style for Bellan Thesis");
theStyle->SetOptStat(0);
theStyle->SetPadBorderMode(0);
theStyle->SetCanvasBorderMode(0);
theStyle->SetPadColor(0);
theStyle->SetCanvasColor(0);
theStyle->SetMarkerStyle(8);
theStyle->SetMarkerSize(0.7);
theStyle->SetPalette(1);
theStyle->SetStatH(0.3);
// theStyle->SetTextFont(132);
// theStyle->SetTitleFont(132);
theStyle->SetTitleBorderSize(1);
// theStyle->SetPalette(1);
theStyle->SetOptStat(0);
theStyle->SetFitFormat("4.4g");
theStyle->SetStatY(0.99);
theStyle->SetStatX(0.99);
theStyle->SetTitleYOffset(1.6);
theStyle->SetLabelSize(0.035, "XYZ");
theStyle->SetPadGridX(true);
theStyle->SetPadGridY(true);
theStyle->SetFrameBorderMode(0);
theStyle->SetTitleFillColor(0);
theStyle->SetLegendBorderSize();
// theStyle->SetCanvasDefH(600);
// theStyle->SetCanvasDefW(400);
//theStyle->SetOptLogy(); //aaa
// theStyle->SetOptLogx();
theStyle->cd();
}
void hggPaperStyle()
{
// Hgg Paper style for plots
TStyle *hggStyle = new TStyle("hggStyle","Hgg Paper Style");
//hggStyle->SetCanvasColor(0);
//hggStyle->SetPadColor(0);
hggStyle->SetPadTickX(0);
hggStyle->SetPadTickY(0);
hggStyle->SetFrameFillColor(0);
hggStyle->SetStatColor(0);
hggStyle->SetOptStat(0);
hggStyle->SetTitleFillColor(0);
hggStyle->SetCanvasBorderMode(0);
hggStyle->SetPadBorderMode(0);
hggStyle->SetFrameBorderMode(0);
hggStyle->SetFrameBorderSize(1);
hggStyle->SetPadColor(kWhite);
hggStyle->SetCanvasColor(kWhite);
hggStyle->SetCanvasDefH(600); //Height of canvas
hggStyle->SetCanvasDefW(800); //Width of canvas
hggStyle->SetCanvasDefX(0); //POsition on screen
hggStyle->SetCanvasDefY(0);
hggStyle->SetPadLeftMargin(0.16);//0.16);
hggStyle->SetPadRightMargin(0.1);//0.02);
hggStyle->SetPadTopMargin(0.085);//0.02);
hggStyle->SetPadBottomMargin(0.12);//0.02);
// For hgg axis titles:
hggStyle->SetTitleColor(1, "XYZ");
hggStyle->SetTitleFont(42, "XYZ");
hggStyle->SetTitleSize(0.05, "XYZ");
hggStyle->SetTitleYOffset(1.5); // => 1.15 if exponents
hggStyle->SetTitleXOffset(1.);//0.9);
// For hgg axis labels:
hggStyle->SetLabelColor(1, "XYZ");
hggStyle->SetLabelFont(42, "XYZ");
hggStyle->SetLabelOffset(0.007, "XYZ");
hggStyle->SetLabelSize(0.045, "XYZ");
// Legends
hggStyle->SetLegendBorderSize(0);
hggStyle->SetLegendFillColor(kWhite);
hggStyle->SetLegendFont(42);
hggStyle->SetFillColor(10);
// Nothing for now
hggStyle->SetTextFont(42);
hggStyle->SetTextSize(0.03);
hggStyle->SetOptTitle(0);
hggStyle->SetOptStat(0);
hggStyle->cd();
}
void drawTracks() {
TStyle *simpleStyle = new TStyle("simpleStyle","");
simpleStyle->SetCanvasColor(0);
simpleStyle->SetFrameFillColor(0);
simpleStyle->SetStatColor(0);
simpleStyle->SetOptStat(0);
simpleStyle->SetTitleFillColor(0);
simpleStyle->SetCanvasBorderMode(0);
simpleStyle->SetPadBorderMode(0);
simpleStyle->SetFrameBorderMode(0);
simpleStyle->cd();
TFile* dataFile = TFile::Open("Runs_123592_123596_123615_123732/output_1.root");
TTree* dataTree = (TTree*)dataFile->Get("myanalysis/pippo");
TFile* mcFile = TFile::Open("MC_MinBias_900GeV/output_12.root");
TTree* mcTree = (TTree*)dataFile->Get("myanalysis/pippo");
std::cout << "1" << std::endl;
Int_t nEventTracks_data;
dataTree->SetBranchAddress("nEventTracks", &nEventTracks_data);
Int_t nChargedHadrons_data;
dataTree->SetBranchAddress("nChargedHadrons", &nChargedHadrons_data);
std::cout << "2" << std::endl;
Int_t nEventTracks_mc;
mcTree->SetBranchAddress("nEventTracks", &nEventTracks_mc);
Int_t nChargedHadrons_mc;
mcTree->SetBranchAddress("nChargedHadrons", &nChargedHadrons_mc);
std::cout << "3" << std::endl;
TH1F* h1_data = new TH1F("data", "", 100, 0., 100.);
h1_data->SetMarkerStyle(20);
TH1F* h1_mc = new TH1F("mc", "", 100, 0., 100.);
h1_mc->SetFillColor(kGray);
std::cout << "4" << std::endl;
for( int iEntry=0; iEntry<dataTree->GetEntries(); ++iEntry) {
dataTree->GetEntry(iEntry);
Float_t trackRatio = (nChargedHadrons_data>0) ? (Float_t)(nEventTracks_data/nChargedHadrons_data) : 0.;
h1_data->Fill( trackRatio );
}
std::cout << "5" << std::endl;
for( int iEntry=0; iEntry<mcTree->GetEntries(); ++iEntry) {
Float_t trackRatio = (nChargedHadrons_mc>0) ? (Float_t)(nEventTracks_mc/nChargedHadrons_mc) : 0.;
h1_mc->Fill( trackRatio );
}
std::cout << "6" << std::endl;
TCanvas* c1 = new TCanvas("c1", "c1", 800, 600);
c1->cd();
h1_data->Draw("E");
h1_mc->Draw("histo same");
c1->SaveAs("tracks.eps");
}
void drawPtSpectrum() {
TStyle *simpleStyle = new TStyle("simpleStyle","");
simpleStyle->SetCanvasColor(0);
simpleStyle->SetFrameFillColor(0);
simpleStyle->SetStatColor(0);
simpleStyle->SetOptStat(0);
simpleStyle->SetTitleFillColor(0);
simpleStyle->SetCanvasBorderMode(0);
simpleStyle->SetPadBorderMode(0);
simpleStyle->SetFrameBorderMode(0);
simpleStyle->cd();
TChain* chain = new TChain("jetTree");
chain->Add("2ndLevelOutputFile_PhotonJet15.root");
chain->Add("2ndLevelOutputFile_PhotonJet30.root");
chain->Add("2ndLevelOutputFile_PhotonJet80.root");
chain->Add("2ndLevelOutputFile_PhotonJet170.root");
chain->Add("2ndLevelOutputFile_PhotonJet300.root");
chain->Add("2ndLevelOutputFile_PhotonJet470.root");
chain->Add("2ndLevelOutputFile_PhotonJet800.root");
Float_t eventWeight;
chain->SetBranchAddress("eventWeight", &eventWeight);
Float_t ptHat;
chain->SetBranchAddress("ptHat", &ptHat);
TH1F* h1_ptSpectrum = new TH1F("ptSpectrum", "", 100, 0., 1000.);
h1_ptSpectrum->SetXTitle("#hat{p}_{T} [GeV/c]");
h1_ptSpectrum->SetMarkerStyle(20);
h1_ptSpectrum->SetMarkerSize(1.4);
Int_t nEntries = chain->GetEntries();
for( Int_t iEntry = 0; iEntry<nEntries; ++iEntry) {
if( (iEntry%100000) == 0 ) std::cout << "Entry: " << iEntry << "/" << nEntries << std::endl;
chain->GetEntry(iEntry);
h1_ptSpectrum->Fill(ptHat, eventWeight);
}
TCanvas* c1 = new TCanvas("c1", "c1", 800, 600);
c1->cd();
c1->SetLogy();
h1_ptSpectrum->Draw("P");
c1->SaveAs("ptSpectrum.eps");
}
// TStyle used for HPlus Studies and AnalysisPack.
//********************************************************************//
void HPlusStyle(void) {
TStyle *HplusStyle = new TStyle("HplusStyle", "HplusStyle");
HplusStyle->SetCanvasBorderMode(0);
HplusStyle->SetCanvasColor(kWhite);
HplusStyle->SetPadBorderMode(0);
HplusStyle->SetPadColor(kWhite);
HplusStyle->SetTitleFillColor(kWhite);
HplusStyle->SetGridColor(0);
HplusStyle->SetFrameBorderMode(0);
HplusStyle->SetFrameBorderSize(1);
HplusStyle->SetFrameFillColor(0);
HplusStyle->SetFrameFillStyle(0);
HplusStyle->SetFrameLineColor(1);
HplusStyle->SetFrameLineStyle(1);
HplusStyle->SetFrameLineWidth(1);
// HplusStyle->SetFillColor(kWhite);
HplusStyle->SetOptTitle(0);
HplusStyle->SetTitleFont(42, "XYZ");
HplusStyle->SetTitleFontSize(0.05);
HplusStyle->SetTitleSize(0.06, "XYZ");
HplusStyle->SetTextFont(42);
HplusStyle->SetTitleXOffset(0.9);
HplusStyle->SetTitleYOffset(1.25);
HplusStyle->SetLabelColor(1, "XYZ");
HplusStyle->SetLabelFont(42, "XYZ");
HplusStyle->SetLabelOffset(0.007, "XYZ");
HplusStyle->SetLabelSize(0.05, "XYZ");
HplusStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame
HplusStyle->SetPadTickY(1); // To get tick marks on the opposite side of the frame
HplusStyle->SetNdivisions(508,"XYZ");
HplusStyle->SetPadTopMargin(0.05);
HplusStyle->SetPadBottomMargin(0.12);
HplusStyle->SetPadLeftMargin(0.16);
HplusStyle->SetPadRightMargin(0.04);//0.02
HplusStyle->SetCanvasDefH(1280); //Height of canvas
HplusStyle->SetCanvasDefW(1024); //Width of canvas
HplusStyle->SetCanvasDefX(0); //POsition on screen
HplusStyle->SetCanvasDefY(0);
HplusStyle->SetPaintTextFormat("5.2f");
// HplusStyle->SetPalette(1, 0); // Pretty palette
HplusStyle->SetOptStat(0);
/// Statistics Box
HplusStyle->SetStatX(0.9);
HplusStyle->SetStatY(0.9);
// HplusStyle->SetStatW(0.19) // default
// HplusStyle->SetStatH(0.1) // default
HplusStyle->SetStatW(0.15);
HplusStyle->SetStatH(0.1);
HplusStyle->cd();
} //eof: void HPlusStyle(void){
/////////
// the Style Section
/////////
void setHLTStyle() {
setTDRStyle();
// TStyle *hltStyle = new TStyle(*tdrStyle);
TStyle *hltStyle = new TStyle("hltStyle","My HLT Styles");
gStyle->Copy(*hltStyle);
hltStyle->SetCanvasColor(-1);
hltStyle->SetCanvasDefH(600);
hltStyle->SetCanvasDefW(600);
hltStyle->SetPadColor(-1);
hltStyle->SetPadGridX(false);
hltStyle->SetPadGridY(false);
hltStyle->SetGridWidth(0.25);
hltStyle->SetFrameFillColor(-1); // Transparent
hltStyle->SetHistFillColor(-1); // Transparent
hltStyle->SetHistFillStyle(0); // None
hltStyle->SetHistLineWidth(3);
hltStyle->SetPadTopMargin(0.08);
hltStyle->SetPadBottomMargin(0.12);
hltStyle->SetPadLeftMargin(0.15);
hltStyle->SetPadRightMargin(0.04);
hltStyle->SetTitleSize(0.05);
hltStyle->SetTitleFillColor(-1); // Transparent
hltStyle->SetTitleH(0.05); // Set the height of the title box
hltStyle->SetTitleW(0.); // Set the width of the title box
hltStyle->SetTitleSize(0.04, "XYZ");
hltStyle->SetTitleOffset(1.2, "X"); // Another way to set the Offset
hltStyle->SetTitleOffset(1.8, "Y"); // Another way to set the Offset
hltStyle->SetLabelSize(0.035, "XYZ");
hltStyle->SetPalette(1,0);
hltStyle->SetFillColor(0); // White
hltStyle->SetFillStyle(4000); // Transparent
hltStyle->SetStatStyle(0);
hltStyle->SetTitleStyle(0);
hltStyle->SetCanvasBorderSize(0);
hltStyle->SetFrameBorderSize(0);
hltStyle->SetLegendBorderSize(0);
hltStyle->SetStatBorderSize(0);
hltStyle->SetTitleBorderSize(0);
hltStyle->cd();
}
// Lauri's Style (LAuri Wendland)
//********************************************************************//
void LAWStyle() {
TStyle *LAWStyle = new TStyle("LAWStyle", "Neat style");
LAWStyle->SetCanvasBorderMode(0);
LAWStyle->SetCanvasColor(kWhite);
LAWStyle->SetPadBorderMode(0);
LAWStyle->SetPadColor(kWhite);
LAWStyle->SetTitleFillColor(kWhite);
LAWStyle->SetGridColor(0);
LAWStyle->SetFrameBorderMode(0);
LAWStyle->SetFrameBorderSize(1);
LAWStyle->SetFrameFillColor(0);
LAWStyle->SetFrameFillStyle(0);
LAWStyle->SetFrameLineColor(1);
LAWStyle->SetFrameLineStyle(1);
LAWStyle->SetFrameLineWidth(1);
LAWStyle->SetFillColor(kWhite);
LAWStyle->SetFillColor(kWhite);
LAWStyle->SetOptTitle(0);
LAWStyle->SetTitleFont(42, "XYZ");
LAWStyle->SetTitleFontSize(0.05);
LAWStyle->SetTitleSize(0.06, "XYZ");
LAWStyle->SetTextFont(42);
LAWStyle->SetTitleXOffset(0.9);
LAWStyle->SetTitleYOffset(1.25);
LAWStyle->SetLabelColor(1, "XYZ");
LAWStyle->SetLabelFont(42, "XYZ");
LAWStyle->SetLabelOffset(0.007, "XYZ");
LAWStyle->SetLabelSize(0.05, "XYZ");
LAWStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame
LAWStyle->SetPadTickY(1); // To get tick marks on the opposite side of the frame
LAWStyle->SetNdivisions(508,"XYZ");
LAWStyle->SetPadTopMargin(0.05);
LAWStyle->SetPadBottomMargin(0.12);
LAWStyle->SetPadLeftMargin(0.16);
LAWStyle->SetPadRightMargin(0.04);//0.02
LAWStyle->SetCanvasDefH(600); //Height of canvas
LAWStyle->SetCanvasDefW(600); //Width of canvas
LAWStyle->SetCanvasDefX(0); //POsition on screen
LAWStyle->SetCanvasDefY(0);
LAWStyle->SetOptTitle(0);
LAWStyle->SetOptStat(0);
LAWStyle->cd();
} //end of: void LAWStyle(){
void setNiceStyle()
{
gROOT->SetStyle("Plain");
// --------------------------------------------------
// define proper canvas style
// --------------------------------------------------
TStyle *MyStyle = new TStyle ("MyStyle", "My style for nicer plots");
Float_t xoff = MyStyle->GetLabelOffset("X"),
yoff = MyStyle->GetLabelOffset("Y"),
zoff = MyStyle->GetLabelOffset("Z");
MyStyle->SetCanvasBorderMode ( 0 );
MyStyle->SetPadBorderMode ( 0 );
MyStyle->SetPadColor ( 0 );
MyStyle->SetCanvasColor ( 0 );
MyStyle->SetTitleColor ( 0 );
MyStyle->SetStatColor ( 0 );
MyStyle->SetTitleBorderSize ( 0 );
MyStyle->SetTitleFillColor ( 0 );
MyStyle->SetTitleH ( 0.07 );
MyStyle->SetTitleW ( 1.00 );
MyStyle->SetTitleFont ( 132 );
MyStyle->SetLabelOffset (1.5*xoff, "X");
MyStyle->SetLabelOffset (1.5*yoff, "Y");
MyStyle->SetLabelOffset (1.5*zoff, "Z");
MyStyle->SetTitleOffset (0.9, "X");
MyStyle->SetTitleOffset (0.9, "Y");
MyStyle->SetTitleOffset (0.9, "Z");
MyStyle->SetTitleSize (0.045, "X");
MyStyle->SetTitleSize (0.045, "Y");
MyStyle->SetTitleSize (0.045, "Z");
MyStyle->SetLabelFont (132, "X");
MyStyle->SetLabelFont (132, "Y");
MyStyle->SetLabelFont (132, "Z");
MyStyle->SetPalette(1);
MyStyle->cd();
}
void setStyle() {
TStyle *setStyle = new TStyle("setStyle","Style for GCTAnalyzer");
// Stuff from plain style
setStyle->SetFrameBorderMode(0);
setStyle->SetCanvasBorderMode(0);
setStyle->SetCanvasColor(kWhite);
setStyle->SetPadBorderMode(0);
setStyle->SetPadColor(kWhite);
setStyle->SetOptStat(0);
setStyle->SetOptTitle(1);
setStyle->SetOptFit(0);
setStyle->SetOptDate(0);
setStyle->cd();
}
void setDefaultStyle()
{
TStyle *defaultStyle = new TStyle("defaultStyle","Default Style");
defaultStyle->SetOptStat(0000);
defaultStyle->SetOptFit(000);
defaultStyle->SetPalette(1);
/////// pad ////////////
defaultStyle->SetPadBorderMode(1);
defaultStyle->SetPadBorderSize(1);
defaultStyle->SetPadColor(0);
defaultStyle->SetPadTopMargin(0.05);
defaultStyle->SetPadBottomMargin(0.13);
defaultStyle->SetPadLeftMargin(0.14);
defaultStyle->SetPadRightMargin(0.02);
/////// canvas /////////
defaultStyle->SetCanvasBorderMode(0);
defaultStyle->SetCanvasColor(0);
defaultStyle->SetCanvasDefH(600);
defaultStyle->SetCanvasDefW(600);
/////// frame //////////
defaultStyle->SetFrameBorderMode(0);
defaultStyle->SetFrameBorderSize(1);
defaultStyle->SetFrameFillColor(0);
defaultStyle->SetFrameLineColor(1);
/////// label //////////
defaultStyle->SetLabelOffset(0.005,"XY");
defaultStyle->SetLabelSize(0.05,"XY");
defaultStyle->SetLabelFont(42,"XY");
/////// title //////////
defaultStyle->SetTitleOffset(1.2,"X");
defaultStyle->SetTitleSize(0.05,"X");
defaultStyle->SetTitleOffset(1.4,"Y");
defaultStyle->SetTitleSize(0.05,"Y");
defaultStyle->SetTitleFont(42, "XYZ");
/////// various ////////
defaultStyle->SetLegendBorderSize(0);
defaultStyle->cd();
}
void drawRxResolutions(const std::string& dataset) {
TStyle *simpleStyle = new TStyle("simpleStyle","");
simpleStyle->SetCanvasColor(0);
simpleStyle->SetFrameFillColor(0);
simpleStyle->SetStatColor(0);
simpleStyle->SetOptStat(0);
simpleStyle->SetTitleFillColor(0);
simpleStyle->SetCanvasBorderMode(0);
simpleStyle->SetPadBorderMode(0);
simpleStyle->SetFrameBorderMode(0);
simpleStyle->cd();
bool logScale=true;
drawResolution(dataset, "Ech");
drawResolution(dataset, "Egamma");
drawResolution(dataset, "Egammanh");
//drawResolution(dataset, "Ee");
drawResolution(dataset, "Enh", logScale);
}
void setStyle() {
TStyle *setStyle = new TStyle("setStyle","Style for GCT Analyzer");
// Stuff from plain style
setStyle->SetCanvasColor(kWhite);
setStyle->SetFrameBorderMode(0);
setStyle->SetCanvasBorderMode(0);
setStyle->SetFrameFillColor(0);
setStyle->SetPadBorderMode(0);
setStyle->SetPadColor(kWhite);
setStyle->SetOptStat(0);
setStyle->SetOptTitle(0);
setStyle->SetOptFit(0);
setStyle->SetOptDate(0);
// Labels and borders
setStyle->SetLabelSize(0.055,"x");
setStyle->SetLabelSize(0.06,"y");
setStyle->SetLabelOffset(0.00,"x");
setStyle->SetLabelOffset(0.00,"y");
setStyle->SetTitleOffset(0.05,"x");
setStyle->SetTitleOffset(0.50,"y");
//setStyle->SetLabelFont(22,"x");
//setStyle->SetLabelFont(22,"y");
//setStyle->SetErrorX(0.0000);
//setStyle->SetTickLength(0.05,"x");
//setStyle->SetTickLength(0.05,"y");
//setStyle->SetLineWidth(0.8);
//setStyle->SetPadTickX(1);
//setStyle->SetPadTickY(1);
setStyle->SetPadLeftMargin(0.15);
setStyle->SetPadTopMargin(0.077);
setStyle->SetPadBottomMargin(0.123);
setStyle->cd();
}
void drawProjections() {
TStyle *simpleStyle = new TStyle("simpleStyle","");
simpleStyle->SetCanvasColor(0);
simpleStyle->SetFrameFillColor(0);
simpleStyle->SetStatColor(0);
simpleStyle->SetOptStat(0);
simpleStyle->SetTitleFillColor(0);
simpleStyle->SetCanvasBorderMode(0);
simpleStyle->SetPadBorderMode(0);
simpleStyle->SetFrameBorderMode(0);
simpleStyle->cd();
draw("PFItCone5_barrel");
draw("PFItCone5_endcap");
draw("caloItCone5_barrel");
draw("caloItCone5_endcap");
draw("Rch90_100_endcap");
}
//------------------------------------------------------------------------------
void PlotAlignmentValidation::setNiceStyle() {
TStyle *MyStyle = new TStyle ("MyStyle", "My style for nicer plots");
Float_t xoff = MyStyle->GetLabelOffset("X"),
yoff = MyStyle->GetLabelOffset("Y"),
zoff = MyStyle->GetLabelOffset("Z");
MyStyle->SetCanvasBorderMode ( 0 );
MyStyle->SetPadBorderMode ( 0 );
MyStyle->SetPadColor ( 0 );
MyStyle->SetCanvasColor ( 0 );
MyStyle->SetTitleColor ( 0 );
MyStyle->SetStatColor ( 0 );
MyStyle->SetTitleBorderSize ( 0 );
MyStyle->SetTitleFillColor ( 0 );
MyStyle->SetTitleH ( 0.07 );
MyStyle->SetTitleW ( 1.00 );
MyStyle->SetTitleFont ( 132 );
MyStyle->SetLabelOffset (1.5*xoff, "X");
MyStyle->SetLabelOffset (1.5*yoff, "Y");
MyStyle->SetLabelOffset (1.5*zoff, "Z");
MyStyle->SetTitleOffset (1.5, "X");
MyStyle->SetTitleOffset (1.2, "Y");
MyStyle->SetTitleOffset (0.9, "Z");
MyStyle->SetTitleSize (0.045, "X");
MyStyle->SetTitleSize (0.045, "Y");
MyStyle->SetTitleSize (0.045, "Z");
MyStyle->SetLabelFont (132, "X");
MyStyle->SetLabelFont (132, "Y");
MyStyle->SetLabelFont (132, "Z");
MyStyle->SetPalette(1);
MyStyle->cd();
}
void setStyle(){
//Style//
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");
// For the canvas:
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetCanvasColor(kWhite);
tdrStyle->SetCanvasDefH(600); //Height of canvas
tdrStyle->SetCanvasDefW(600); //Width of canvas
tdrStyle->SetCanvasDefX(0); //POsition on screen
tdrStyle->SetCanvasDefY(0);
// For the Pad:
tdrStyle->SetPadBorderMode(0);
tdrStyle->SetPadColor(kWhite);
tdrStyle->SetPadGridX(false);
tdrStyle->SetPadGridY(false);
tdrStyle->SetGridColor(0);
tdrStyle->SetGridStyle(3);
tdrStyle->SetGridWidth(1);
// For the frame:
tdrStyle->SetFrameBorderMode(0);
tdrStyle->SetFrameBorderSize(1);
tdrStyle->SetFrameFillColor(0);
tdrStyle->SetFrameFillStyle(0);
tdrStyle->SetFrameLineColor(1);
tdrStyle->SetFrameLineStyle(1);
tdrStyle->SetFrameLineWidth(1);
tdrStyle->SetPadTopMargin(0.05);
tdrStyle->SetPadBottomMargin(0.13);
tdrStyle->SetPadLeftMargin(0.16);
tdrStyle->SetPadRightMargin(0.02);
// For the Global title:
tdrStyle->SetOptTitle(0);
tdrStyle->SetTitleFont(42);
tdrStyle->SetTitleColor(1);
tdrStyle->SetTitleTextColor(1);
tdrStyle->SetTitleFillColor(10);
tdrStyle->SetTitleFontSize(0.05);
// For the axis titles:
tdrStyle->SetTitleColor(1, "XYZ");
tdrStyle->SetTitleFont(42, "XYZ");
tdrStyle->SetTitleSize(0.06, "XYZ");
tdrStyle->SetTitleXOffset(0.9);
tdrStyle->SetTitleYOffset(0.9);
// For the axis labels:
tdrStyle->SetLabelColor(1, "XYZ");
tdrStyle->SetLabelFont(42, "XYZ");
tdrStyle->SetLabelOffset(0.007, "XYZ");
tdrStyle->SetLabelSize(0.05, "XYZ");
// For the axis:
tdrStyle->SetAxisColor(1, "XYZ");
//tdrStyle->SetStripDecimals(kTRUE);
tdrStyle->SetTickLength(0.03, "XYZ");
tdrStyle->SetNdivisions(510, "XYZ");
tdrStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame
tdrStyle->SetPadTickY(1);
tdrStyle->SetPaperSize(20.,20.);
tdrStyle->SetHatchesLineWidth(5);
tdrStyle->SetHatchesSpacing(0.05);
tdrStyle->SetOptStat(0);
tdrStyle->cd();
}
void LatinoStyle2()
{
TStyle* GloStyle;
GloStyle = gStyle;
TStyle* LatinosStyle = new TStyle("LatinosStyle", "LatinosStyle");
gStyle = LatinosStyle;
//----------------------------------------------------------------------------
// Canvas
//----------------------------------------------------------------------------
LatinosStyle->SetCanvasBorderMode( 0);
LatinosStyle->SetCanvasBorderSize( 10);
LatinosStyle->SetCanvasColor ( 0);
LatinosStyle->SetCanvasDefH (794);
LatinosStyle->SetCanvasDefW (550);
LatinosStyle->SetCanvasDefX ( 10);
LatinosStyle->SetCanvasDefY ( 10);
//----------------------------------------------------------------------------
// Pad
//----------------------------------------------------------------------------
LatinosStyle->SetPadBorderMode ( 0);
LatinosStyle->SetPadBorderSize ( 10);
LatinosStyle->SetPadColor ( 0);
LatinosStyle->SetPadBottomMargin(0.20);
LatinosStyle->SetPadTopMargin (0.08);
LatinosStyle->SetPadLeftMargin (0.18);
LatinosStyle->SetPadRightMargin (0.05);
//----------------------------------------------------------------------------
// Frame
//----------------------------------------------------------------------------
LatinosStyle->SetFrameFillStyle ( 0);
LatinosStyle->SetFrameFillColor ( 0);
LatinosStyle->SetFrameLineColor ( 1);
LatinosStyle->SetFrameLineStyle ( 0);
LatinosStyle->SetFrameLineWidth ( 2);
LatinosStyle->SetFrameBorderMode( 0);
LatinosStyle->SetFrameBorderSize(10);
//----------------------------------------------------------------------------
// Hist
//----------------------------------------------------------------------------
LatinosStyle->SetHistFillColor(0);
LatinosStyle->SetHistFillStyle(1);
LatinosStyle->SetHistLineColor(1);
LatinosStyle->SetHistLineStyle(0);
LatinosStyle->SetHistLineWidth(1);
//----------------------------------------------------------------------------
// Axis
//----------------------------------------------------------------------------
LatinosStyle->SetLabelFont ( 42, "xyz");
LatinosStyle->SetLabelOffset(0.015, "xyz");
LatinosStyle->SetLabelSize (0.050, "xyz");
LatinosStyle->SetNdivisions ( 505, "xyz");
LatinosStyle->SetTitleFont ( 42, "xyz");
LatinosStyle->SetTitleSize (0.050, "xyz");
// LatinosStyle->SetNdivisions ( -503, "y");
LatinosStyle->SetTitleOffset( 1.4, "x");
LatinosStyle->SetTitleOffset( 1.2, "y");
LatinosStyle->SetPadTickX ( 1); // Tick marks on the opposite side of the frame
LatinosStyle->SetPadTickY ( 1); // Tick marks on the opposite side of the frame
//----------------------------------------------------------------------------
// Title
//----------------------------------------------------------------------------
LatinosStyle->SetTitleBorderSize( 0);
LatinosStyle->SetTitleFillColor ( 10);
LatinosStyle->SetTitleAlign ( 12);
LatinosStyle->SetTitleFontSize (0.045);
LatinosStyle->SetTitleX (0.590);
LatinosStyle->SetTitleY (0.830);
LatinosStyle->SetTitleFont(42, "");
//----------------------------------------------------------------------------
// Stat
//----------------------------------------------------------------------------
// LatinosStyle->SetOptStat (1110);
LatinosStyle->SetOptStat ( 0);
LatinosStyle->SetStatBorderSize( 0);
LatinosStyle->SetStatColor ( 10);
LatinosStyle->SetStatFont ( 42);
LatinosStyle->SetStatX (0.94);
LatinosStyle->SetStatY (0.91);
const Int_t NRGBs = 5;
const Int_t NCont = 255;
Double_t stops[NRGBs] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
Double_t red[NRGBs] = { 0.00, 0.00, 0.87, 1.00, 0.51 };
Double_t green[NRGBs] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
Double_t blue[NRGBs] = { 0.51, 1.00, 0.12, 0.00, 0.00 };
TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
LatinosStyle->SetNumberContours(NCont);
return;
}
void doStackGen() {
TStyle *simpleStyle = new TStyle("simpleStyle","");
simpleStyle->SetCanvasColor(0);
simpleStyle->SetFrameFillColor(0);
simpleStyle->SetStatColor(0);
simpleStyle->SetOptStat(0);
simpleStyle->SetTitleFillColor(0);
simpleStyle->SetCanvasBorderMode(0);
simpleStyle->SetPadBorderMode(0);
simpleStyle->SetFrameBorderMode(0);
simpleStyle->cd();
TFile* inFile = TFile::Open("RxGenOutputFile.root");
TH1D* h1_RchGen_vs_eta = (TH1D*)inFile->Get("RchGen_vs_eta");
TH1D* h1_RgammaGen_vs_eta = (TH1D*)inFile->Get("RgammaGen_vs_eta");
TH1D* h1_RnhGen_vs_eta = (TH1D*)inFile->Get("RnhGen_vs_eta");
// TH1D* h1_RmuGen_vs_eta = (TH1D*)inFile->Get("RmuGen_vs_eta_PFItCone5");
int nBinsX_RchGen = h1_RchGen_vs_eta->GetNbinsX();
int nBinsX_RgammaGen = h1_RgammaGen_vs_eta->GetNbinsX();
int nBinsX_RnhGen = h1_RnhGen_vs_eta->GetNbinsX();
// int nBinsX_RmuGen = h1_RmuGen_vs_eta->GetNbinsX();
double minX_RchGen = h1_RchGen_vs_eta->GetXaxis()->GetXmin();
double maxX_RchGen = h1_RchGen_vs_eta->GetXaxis()->GetXmax();
double minX_RgammaGen = h1_RgammaGen_vs_eta->GetXaxis()->GetXmin();
double maxX_RgammaGen = h1_RgammaGen_vs_eta->GetXaxis()->GetXmax();
double minX_RnhGen = h1_RnhGen_vs_eta->GetXaxis()->GetXmin();
double maxX_RnhGen = h1_RnhGen_vs_eta->GetXaxis()->GetXmax();
// double minX_RmuGen = h1_RmuGen_vs_eta->GetXaxis()->GetXmin();
// double maxX_RmuGen = h1_RmuGen_vs_eta->GetXaxis()->GetXmax();
TH1D* h1_RchGen_vs_eta_stack = new TH1D("RchGen_vs_eta_stack", "", nBinsX_RchGen, minX_RchGen, maxX_RchGen);
h1_RchGen_vs_eta_stack->SetFillColor(kRed);
h1_RchGen_vs_eta_stack->SetXTitle("#eta^{GEN}");
TH1D* h1_RgammaGen_vs_eta_stack = new TH1D("RgammaGen_vs_eta_stack", "", nBinsX_RgammaGen, minX_RgammaGen, maxX_RgammaGen);
h1_RgammaGen_vs_eta_stack->SetFillColor(kGreen);
h1_RgammaGen_vs_eta_stack->SetXTitle("#eta^{GEN}");
TH1D* h1_RnhGen_vs_eta_stack = new TH1D("RnhGen_vs_eta_stack", "", nBinsX_RnhGen, minX_RnhGen, maxX_RnhGen);
h1_RnhGen_vs_eta_stack->SetFillColor(kBlue);
h1_RnhGen_vs_eta_stack->SetXTitle("#eta^{GEN}");
// TH1D* h1_RmuGen_vs_eta_stack = new TH1D("RmuGen_vs_eta", "", nBinsX_RmuGen, minX_RmuGen, maxX_RmuGen);
// h1_RmuGen_vs_eta_stack->SetFillColor(kYellow);
// h1_RmuGen_vs_eta_stack->SetXTitle("#eta^{GEN}");
for(int i=1; i<101; ++i) {
h1_RchGen_vs_eta_stack->SetBinContent( i, h1_RchGen_vs_eta->GetBinContent(i) );
h1_RgammaGen_vs_eta_stack->SetBinContent( i, h1_RgammaGen_vs_eta->GetBinContent(i) );
h1_RnhGen_vs_eta_stack->SetBinContent( i, h1_RnhGen_vs_eta->GetBinContent(i) );
// h1_RmuGen_vs_eta_stack->SetBinContent( i, h1_RmuGen_vs_eta->GetBinContent(i) );
}
THStack* stack = new THStack("stack", "");
stack->Add(h1_RchGen_vs_eta_stack);
stack->Add(h1_RgammaGen_vs_eta_stack);
stack->Add(h1_RnhGen_vs_eta_stack);
// stack->Add(h1_RmuGen_vs_eta_stack);
stack->Draw();
stack->GetXaxis()->SetTitle("#eta^{GEN}");
stack->GetYaxis()->SetTitle("% of Jet Energy");
TCanvas* c1 = new TCanvas("c1", "c1", 800, 600);
c1->cd();
stack->Draw();
c1->SaveAs("Plots/genstack_vs_eta.eps");
TH1D* h1_Rnh_neutron_vs_eta = (TH1D*)inFile->Get("Rnh_neutron_vs_eta");
TH1D* h1_Rnh_K0L_vs_eta = (TH1D*)inFile->Get("Rnh_K0L_vs_eta");
TH1D* h1_Rnh_K0S_vs_eta = (TH1D*)inFile->Get("Rnh_K0S_vs_eta");
TH1D* h1_Rnh_lambda_vs_eta = (TH1D*)inFile->Get("Rnh_lambda_vs_eta");
TH1D* h1_Rnh_csi_vs_eta = (TH1D*)inFile->Get("Rnh_csi_vs_eta");
int nBinsX_Rnh_neutron = h1_Rnh_neutron_vs_eta->GetNbinsX();
int nBinsX_Rnh_K0L = h1_Rnh_K0L_vs_eta->GetNbinsX();
int nBinsX_Rnh_K0S = h1_Rnh_K0S_vs_eta->GetNbinsX();
int nBinsX_Rnh_lambda = h1_Rnh_lambda_vs_eta->GetNbinsX();
int nBinsX_Rnh_csi = h1_Rnh_csi_vs_eta->GetNbinsX();
double minX_Rnh_neutron = h1_Rnh_neutron_vs_eta->GetXaxis()->GetXmin();
double maxX_Rnh_neutron = h1_Rnh_neutron_vs_eta->GetXaxis()->GetXmax();
double minX_Rnh_K0L = h1_Rnh_K0L_vs_eta->GetXaxis()->GetXmin();
double maxX_Rnh_K0L = h1_Rnh_K0L_vs_eta->GetXaxis()->GetXmax();
double minX_Rnh_K0S = h1_Rnh_K0S_vs_eta->GetXaxis()->GetXmin();
double maxX_Rnh_K0S = h1_Rnh_K0S_vs_eta->GetXaxis()->GetXmax();
double minX_Rnh_lambda = h1_Rnh_lambda_vs_eta->GetXaxis()->GetXmin();
double maxX_Rnh_lambda = h1_Rnh_lambda_vs_eta->GetXaxis()->GetXmax();
double minX_Rnh_csi = h1_Rnh_csi_vs_eta->GetXaxis()->GetXmin();
double maxX_Rnh_csi = h1_Rnh_csi_vs_eta->GetXaxis()->GetXmax();
TH1D* h1_Rnh_neutron_vs_eta_stack = new TH1D("Rnh_neutron_vs_eta_stack", "", nBinsX_Rnh_neutron, minX_Rnh_neutron, maxX_Rnh_neutron);
h1_Rnh_neutron_vs_eta_stack->SetFillColor(kBlue+4);
h1_Rnh_neutron_vs_eta_stack->SetXTitle("#eta^{GEN}");
TH1D* h1_Rnh_K0L_vs_eta_stack = new TH1D("Rnh_K0L_vs_eta_stack", "", nBinsX_Rnh_K0L, minX_Rnh_K0L, maxX_Rnh_K0L);
h1_Rnh_K0L_vs_eta_stack->SetFillColor(kBlue);
h1_Rnh_K0L_vs_eta_stack->SetXTitle("#eta^{GEN}");
TH1D* h1_Rnh_K0S_vs_eta_stack = new TH1D("Rnh_K0S_vs_eta_stack", "", nBinsX_Rnh_K0S, minX_Rnh_K0S, maxX_Rnh_K0S);
h1_Rnh_K0S_vs_eta_stack->SetFillColor(kCyan);
h1_Rnh_K0S_vs_eta_stack->SetXTitle("#eta^{GEN}");
TH1D* h1_Rnh_lambda_vs_eta_stack = new TH1D("Rnh_lambda_vs_eta_stack", "", nBinsX_Rnh_lambda, minX_Rnh_lambda, maxX_Rnh_lambda);
h1_Rnh_lambda_vs_eta_stack->SetFillColor(kCyan+2);
h1_Rnh_lambda_vs_eta_stack->SetXTitle("#eta^{GEN}");
TH1D* h1_Rnh_csi_vs_eta_stack = new TH1D("Rnh_csi_vs_eta_stack", "", nBinsX_Rnh_csi, minX_Rnh_csi, maxX_Rnh_csi);
h1_Rnh_csi_vs_eta_stack->SetFillColor(kCyan+4);
h1_Rnh_csi_vs_eta_stack->SetXTitle("#eta^{GEN}");
for(int i=1; i<101; ++i) {
h1_Rnh_neutron_vs_eta_stack->SetBinContent( i, h1_Rnh_neutron_vs_eta->GetBinContent(i) );
h1_Rnh_K0L_vs_eta_stack->SetBinContent( i, h1_Rnh_K0L_vs_eta->GetBinContent(i) );
h1_Rnh_K0S_vs_eta_stack->SetBinContent( i, h1_Rnh_K0S_vs_eta->GetBinContent(i) );
h1_Rnh_lambda_vs_eta_stack->SetBinContent( i, h1_Rnh_lambda_vs_eta->GetBinContent(i) );
h1_Rnh_csi_vs_eta_stack->SetBinContent( i, h1_Rnh_csi_vs_eta->GetBinContent(i) );
}
THStack* stackNH = new THStack("stackNH", "");
stackNH->Add(h1_Rnh_neutron_vs_eta_stack);
stackNH->Add(h1_Rnh_K0L_vs_eta_stack);
stackNH->Add(h1_Rnh_K0S_vs_eta_stack);
stackNH->Add(h1_Rnh_lambda_vs_eta_stack);
stackNH->Add(h1_Rnh_csi_vs_eta_stack);
stackNH->Draw();
stackNH->GetXaxis()->SetTitle("#eta^{GEN}");
stackNH->GetYaxis()->SetTitle("% of Jet Neutral Hadron Energy");
TCanvas* c3 = new TCanvas("c3", "c3", 800, 600);
c3->cd();
stackNH->Draw();
c3->SaveAs("Plots/genstackNH_vs_eta.eps");
/*
int nBins = 21;
Double_t Lower[nBins];
for(int iLower=0; iLower<11; ++iLower)
Lower[iLower] = iLower*10.;
Lower[11] = 120.;
Lower[12] = 140.;
Lower[13] = 160.;
Lower[14] = 180.;
Lower[15] = 200.;
Lower[16] = 250.;
Lower[17] = 300.;
Lower[18] = 400.;
Lower[19] = 500.;
Lower[20] = 600.;
TH1D* h1_RchGen_vs_pt = (TH1D*)inFile->Get("RchGen_vs_pt_PFItCone5");
TH1D* h1_RgammaGen_vs_pt = (TH1D*)inFile->Get("RgammaGen_vs_pt_PFItCone5");
TH1D* h1_RnhGen_vs_pt = (TH1D*)inFile->Get("RnhGen_vs_pt_PFItCone5");
// TH1D* h1_RmuGen_vs_pt = (TH1D*)inFile->Get("RmuGen_vs_pt_PFItCone5");
TH1D* h1_RchGen_vs_pt_stack = new TH1D("RchGen_vs_pt", "", nBins-2, Lower);
h1_RchGen_vs_pt_stack->SetFillColor(kRed);
h1_RchGen_vs_pt_stack->SetXTitle("p_{T}^{GEN}");
TH1D* h1_RgammaGen_vs_pt_stack = new TH1D("RgammaGen_vs_pt", "", nBins-2, Lower);
h1_RgammaGen_vs_pt_stack->SetFillColor(kGreen);
h1_RgammaGen_vs_pt_stack->SetXTitle("p_{T}^{GEN}");
TH1D* h1_RnhGen_vs_pt_stack = new TH1D("RnhGen_vs_pt", "", nBins-2, Lower);
h1_RnhGen_vs_pt_stack->SetFillColor(kBlue);
h1_RnhGen_vs_pt_stack->SetXTitle("p_{T}^{GEN}");
// TH1D* h1_RmuGen_vs_pt_stack = new TH1D("RmuGen_vs_pt", "", nBins-2, Lower);
// h1_RmuGen_vs_pt_stack->SetFillColor(kYellow);
// h1_RmuGen_vs_pt_stack->SetXTitle("p_{T}^{GEN}");
for(int i=1; i<21; ++i) {
h1_RchGen_vs_pt_stack->SetBinContent( i, h1_RchGen_vs_pt->GetBinContent(i) );
h1_RgammaGen_vs_pt_stack->SetBinContent( i, h1_RgammaGen_vs_pt->GetBinContent(i) );
h1_RnhGen_vs_pt_stack->SetBinContent( i, h1_RnhGen_vs_pt->GetBinContent(i) );
// h1_RmuGen_vs_pt_stack->SetBinContent( i, h1_RmuGen_vs_pt->GetBinContent(i) );
}
THStack* stack2 = new THStack("stack2", "");
stack2->Add(h1_RchGen_vs_pt_stack);
stack2->Add(h1_RgammaGen_vs_pt_stack);
stack2->Add(h1_RnhGen_vs_pt_stack);
// stack2->Add(h1_RmuGen_vs_pt_stack);
stack2->Draw();
stack2->GetXaxis()->SetTitle("p_{T}^{GEN} [GeV/c]");
stack2->GetYaxis()->SetTitle("% of Jet Energy");
TCanvas* c2 = new TCanvas("c2", "c2", 800, 600);
c2->cd();
stack2->Draw();
c2->SaveAs("Plots/genstack_vs_pt.eps");
TH1D* h1_RchGen_vs_pt_barrel = (TH1D*)inFile->Get("RchGen_vs_pt_barrel_PFItCone5");
TH1D* h1_RgammaGen_vs_pt_barrel = (TH1D*)inFile->Get("RgammaGen_vs_pt_barrel_PFItCone5");
TH1D* h1_RnhGen_vs_pt_barrel = (TH1D*)inFile->Get("RnhGen_vs_pt_barrel_PFItCone5");
// TH1D* h1_RmuGen_vs_pt_barrel = (TH1D*)inFile->Get("RmuGen_vs_pt_barrel_PFItCone5");
TH1D* h1_RchGen_vs_pt_barrel_stack = new TH1D("RchGen_vs_pt_barrel", "", nBins-2, Lower);
h1_RchGen_vs_pt_barrel_stack->SetFillColor(kRed);
h1_RchGen_vs_pt_barrel_stack->SetXTitle("p_{T}^{GEN}");
TH1D* h1_RgammaGen_vs_pt_barrel_stack = new TH1D("RgammaGen_vs_pt_barrel", "", nBins-2, Lower);
h1_RgammaGen_vs_pt_barrel_stack->SetFillColor(kGreen);
h1_RgammaGen_vs_pt_barrel_stack->SetXTitle("p_{T}^{GEN}");
TH1D* h1_RnhGen_vs_pt_barrel_stack = new TH1D("RnhGen_vs_pt_barrel", "", nBins-2, Lower);
h1_RnhGen_vs_pt_barrel_stack->SetFillColor(kBlue);
h1_RnhGen_vs_pt_barrel_stack->SetXTitle("p_{T}^{GEN}");
// TH1D* h1_RmuGen_vs_pt_barrel_stack = new TH1D("RmuGen_vs_pt_barrel", "", nBins-2, Lower);
// h1_RmuGen_vs_pt_barrel_stack->SetFillColor(kYellow);
// h1_RmuGen_vs_pt_barrel_stack->SetXTitle("p_{T}^{GEN}");
for(int i=1; i<21; ++i) {
h1_RchGen_vs_pt_barrel_stack->SetBinContent( i, h1_RchGen_vs_pt_barrel->GetBinContent(i) );
h1_RgammaGen_vs_pt_barrel_stack->SetBinContent( i, h1_RgammaGen_vs_pt_barrel->GetBinContent(i) );
h1_RnhGen_vs_pt_barrel_stack->SetBinContent( i, h1_RnhGen_vs_pt_barrel->GetBinContent(i) );
// h1_RmuGen_vs_pt_barrel_stack->SetBinContent( i, h1_RmuGen_vs_pt_barrel->GetBinContent(i) );
}
THStack* stack3 = new THStack("stack3", "");
stack3->Add(h1_RchGen_vs_pt_barrel_stack);
stack3->Add(h1_RgammaGen_vs_pt_barrel_stack);
stack3->Add(h1_RnhGen_vs_pt_barrel_stack);
// stack3->Add(h1_RmuGen_vs_pt_barrel_stack);
stack3->Draw();
stack3->GetXaxis()->SetTitle("p_{T}^{GEN} [GeV/c]");
stack3->GetYaxis()->SetTitle("% of Jet Energy");
c2->cd();
stack3->Draw();
c2->SaveAs("Plots/genstack_vs_pt_barrel.eps");
TH1D* h1_RchGen_vs_pt_endcap = (TH1D*)inFile->Get("RchGen_vs_pt_endcap_PFItCone5");
TH1D* h1_RgammaGen_vs_pt_endcap = (TH1D*)inFile->Get("RgammaGen_vs_pt_endcap_PFItCone5");
TH1D* h1_RnhGen_vs_pt_endcap = (TH1D*)inFile->Get("RnhGen_vs_pt_endcap_PFItCone5");
// TH1D* h1_RmuGen_vs_pt_endcap = (TH1D*)inFile->Get("RmuGen_vs_pt_endcap_PFItCone5");
TH1D* h1_RchGen_vs_pt_endcap_stack = new TH1D("RchGen_vs_pt_endcap", "", nBins-2, Lower);
h1_RchGen_vs_pt_endcap_stack->SetFillColor(kRed);
h1_RchGen_vs_pt_endcap_stack->SetXTitle("p_{T}^{GEN}");
TH1D* h1_RgammaGen_vs_pt_endcap_stack = new TH1D("RgammaGen_vs_pt_endcap", "", nBins-2, Lower);
h1_RgammaGen_vs_pt_endcap_stack->SetFillColor(kGreen);
h1_RgammaGen_vs_pt_endcap_stack->SetXTitle("p_{T}^{GEN}");
TH1D* h1_RnhGen_vs_pt_endcap_stack = new TH1D("RnhGen_vs_pt_endcap", "", nBins-2, Lower);
h1_RnhGen_vs_pt_endcap_stack->SetFillColor(kBlue);
h1_RnhGen_vs_pt_endcap_stack->SetXTitle("p_{T}^{GEN}");
// TH1D* h1_RmuGen_vs_pt_endcap_stack = new TH1D("RmuGen_vs_pt_endcap", "", nBins-2, Lower);
// h1_RmuGen_vs_pt_endcap_stack->SetFillColor(kYellow);
// h1_RmuGen_vs_pt_endcap_stack->SetXTitle("p_{T}^{GEN}");
for(int i=1; i<21; ++i) {
h1_RchGen_vs_pt_endcap_stack->SetBinContent( i, h1_RchGen_vs_pt_endcap->GetBinContent(i) );
h1_RgammaGen_vs_pt_endcap_stack->SetBinContent( i, h1_RgammaGen_vs_pt_endcap->GetBinContent(i) );
h1_RnhGen_vs_pt_endcap_stack->SetBinContent( i, h1_RnhGen_vs_pt_endcap->GetBinContent(i) );
//h1_RmuGen_vs_pt_endcap_stack->SetBinContent( i, h1_RmuGen_vs_pt_endcap->GetBinContent(i) );
}
THStack* stack4 = new THStack("stack4", "");
stack4->Add(h1_RchGen_vs_pt_endcap_stack);
stack4->Add(h1_RgammaGen_vs_pt_endcap_stack);
stack4->Add(h1_RnhGen_vs_pt_endcap_stack);
// stack4->Add(h1_RmuGen_vs_pt_endcap_stack);
stack4->Draw();
stack4->GetXaxis()->SetTitle("p_{T}^{GEN} [GeV/c]");
stack4->GetYaxis()->SetTitle("% of Jet Energy");
c2->cd();
stack4->Draw();
c2->SaveAs("Plots/genstack_vs_pt_endcap.eps");
*/
}
void roc_plot(){
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");
// For the canvas:
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetCanvasColor(kWhite);
tdrStyle->SetCanvasDefH(600); //Height of canvas
tdrStyle->SetCanvasDefW(600); //Width of canvas
tdrStyle->SetCanvasDefX(0); //POsition on screen
tdrStyle->SetCanvasDefY(0);
// For the Pad:
tdrStyle->SetPadBorderMode(0);
tdrStyle->SetPadColor(kWhite);
tdrStyle->SetPadGridX(true);
tdrStyle->SetPadGridY(true);
tdrStyle->SetGridColor(0);
tdrStyle->SetGridStyle(3);
tdrStyle->SetGridWidth(1);
// For the frame:
tdrStyle->SetFrameBorderMode(0);
tdrStyle->SetFrameBorderSize(1);
tdrStyle->SetFrameFillColor(0);
tdrStyle->SetFrameFillStyle(0);
tdrStyle->SetFrameLineColor(1);
tdrStyle->SetFrameLineStyle(1);
tdrStyle->SetFrameLineWidth(1);
// For the histo:
tdrStyle->SetHistLineColor(1);
tdrStyle->SetHistLineStyle(0);
tdrStyle->SetHistLineWidth(1);
tdrStyle->SetEndErrorSize(2);
tdrStyle->SetErrorX(0.);
tdrStyle->SetMarkerStyle(20);
//For the fit/function:
tdrStyle->SetOptFit(1);
tdrStyle->SetFitFormat("5.4g");
tdrStyle->SetFuncColor(2);
tdrStyle->SetFuncStyle(1);
tdrStyle->SetFuncWidth(1);
//For the date:
tdrStyle->SetOptDate(0);
// For the statistics box:
tdrStyle->SetOptFile(0);
tdrStyle->SetOptStat("emr"); // To display the mean and RMS: SetOptStat("mr");
tdrStyle->SetStatColor(kWhite);
tdrStyle->SetStatFont(42);
tdrStyle->SetStatFontSize(0.025);
tdrStyle->SetStatTextColor(1);
tdrStyle->SetStatFormat("6.4g");
tdrStyle->SetStatBorderSize(1);
tdrStyle->SetStatH(0.1);
tdrStyle->SetStatW(0.15);
// Margins:
tdrStyle->SetPadTopMargin(0.05);
tdrStyle->SetPadBottomMargin(0.13);
tdrStyle->SetPadLeftMargin(0.13);
tdrStyle->SetPadRightMargin(0.05);
// For the Global title:
tdrStyle->SetOptTitle(0);
tdrStyle->SetTitleFont(42);
tdrStyle->SetTitleColor(1);
tdrStyle->SetTitleTextColor(1);
tdrStyle->SetTitleFillColor(10);
tdrStyle->SetTitleFontSize(0.05);
// For the axis titles:
tdrStyle->SetTitleColor(1, "XYZ");
tdrStyle->SetTitleFont(42, "XYZ");
tdrStyle->SetTitleSize(0.06, "XYZ");
tdrStyle->SetTitleXOffset(0.9);
tdrStyle->SetTitleYOffset(1.05);
// For the axis labels:
tdrStyle->SetLabelColor(1, "XYZ");
tdrStyle->SetLabelFont(42, "XYZ");
tdrStyle->SetLabelOffset(0.007, "XYZ");
tdrStyle->SetLabelSize(0.05, "XYZ");
// For the axis:
tdrStyle->SetAxisColor(1, "XYZ");
tdrStyle->SetStripDecimals(kTRUE);
tdrStyle->SetTickLength(0.03, "XYZ");
tdrStyle->SetNdivisions(510, "XYZ");
tdrStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame
tdrStyle->SetPadTickY(1);
// Change for log plots:
tdrStyle->SetOptLogx(0);
tdrStyle->SetOptLogy(0);
tdrStyle->SetOptLogz(0);
// Postscript options:
tdrStyle->SetPaperSize(20.,20.);
tdrStyle->SetPalette(1);
gROOT -> ForceStyle();
tdrStyle->cd();
gStyle->SetOptStat(0111);
gStyle->SetOptFit(1111);
int bkg = 787623.93;
int sig = 9313.;
const int n = 5;
//minmet
//Double_t cut[n] = {10,15,20,25,30,35,40,45,50,55} ;
//pvalue
//Double_t cut[n] = {0,0.5,1,1.5,2,2.5,3,3.5,4.0};
//metsig
Double_t cut[n] = {5,10,15,20,25};
//minmet
//Double_t sensitivity_bkg_n[n]={660142.26/bkg,551026.78/bkg,444547.21/bkg,354649.92/bkg,283093.63/bkg ,227934.73/bkg ,187564.23/bkg,158191.46/bkg,103089.49/bkg,68423.49/bkg};
//Double_t sensitivity_sig_n[n]={9274./sig ,9233./sig,9172./sig,9089./sig,8995./sig,8870./sig,8718./sig,8539./sig,7398./sig,6497./sig};
//pvalue
//Double_t sensitivity_bkg_n[n]={787623.93/bkg,290416.84/bkg, 166165.65/bkg, 87792.43/bkg, 49447.43/bkg, 30536.76/bkg, 19869.09/bkg, 13705.00/bkg, 10075.21/bkg};
//Double_t sensitivity_sig_n[n]={9313./sig,9005./sig, 8499./sig,7547./sig,6508./sig,5579./sig,4832./sig,4213./sig,3649./sig};
//metsig
Double_t sensitivity_bkg_n[n]={579088.98/bkg,94310.14/bkg, 13034.68/bkg,2503.93./bkg,756.13./bkg};
Double_t sensitivity_sig_n[n]={9045./sig,7492./sig,5899./sig,4654./sig,3726./sig};
sensitivity_bkg = new TGraph(n,cut,sensitivity_bkg_n);
sensitivity_sig = new TGraph(n,cut,sensitivity_sig_n);
//Double_t xdummy[2] = {10,55};
//Double_t xdummy[2] = {-0.1,3};
Double_t xdummy[2] = {0,25};
Double_t ydummy[2] = {0.0,1.0};
dummy = new TGraph(2,xdummy,ydummy);
dummy->SetLineColor(0);
dummy->SetLineWidth(0);
dummy->SetMarkerColor(0);
dummy->SetMarkerSize(0);
dummy->Draw("AP*");
dummy->GetYaxis()->SetTitle("Efficiency");
//dummy->GetXaxis()->SetTitle("Minimized #slash{E}_{T} Cut [GeV]");
//dummy->GetXaxis()->SetTitle("|Log10(P_{#Chi^{2}})|");
dummy->GetXaxis()->SetTitle("Met Significance");
dummy->GetYaxis()->SetLabelSize(0.04);
dummy->GetXaxis()->SetLabelSize(0.04);
sensitivity_bkg->SetLineColor(2);
sensitivity_bkg->SetMarkerColor(2);
sensitivity_sig->SetLineColor(4);
sensitivity_sig->SetMarkerColor(4);
sensitivity_bkg->Draw("LPSame");
sensitivity_sig->Draw("LPSame");
//legend = new TLegend(.70,.70,.90,.90,"Final Selection e/#mu+jets");
legend = new TLegend(0.1845638,0.2027972,0.3674497,0.3024476,NULL,"brNDC");
legend->AddEntry(sensitivity_bkg, "#gamma + Jets", "lp");
legend->AddEntry(sensitivity_sig, "Z (#nu#nu) #gamma" , "lp");
legend->SetShadowColor(0);
legend->SetFillColor(0);
legend->SetLineColor(0);
legend->Draw("Same");
latex2 = new TLatex();
latex2->SetNDC();
latex2->SetTextSize(0.04);
latex2->SetTextAlign(31);
latex2->DrawLatex(0.87, 0.96, "CMS Preliminary, #sqrt{s} = 8 TeV");
}
void setTDRStyle() {
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");
// For the canvas:
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetCanvasColor(kWhite);
tdrStyle->SetCanvasDefH(600); //Height of canvas
tdrStyle->SetCanvasDefW(600); //Width of canvas
tdrStyle->SetCanvasDefX(0); //POsition on screen
tdrStyle->SetCanvasDefY(0);
// For the Pad:
tdrStyle->SetPadBorderMode(0);
// tdrStyle->SetPadBorderSize(Width_t size = 1);
tdrStyle->SetPadColor(kWhite);
tdrStyle->SetPadGridX(false);
tdrStyle->SetPadGridY(false);
tdrStyle->SetGridColor(0);
tdrStyle->SetGridStyle(3);
tdrStyle->SetGridWidth(1);
// For the frame:
tdrStyle->SetFrameBorderMode(0);
tdrStyle->SetFrameBorderSize(1);
tdrStyle->SetFrameFillColor(0);
tdrStyle->SetFrameFillStyle(0);
tdrStyle->SetFrameLineColor(1);
tdrStyle->SetFrameLineStyle(1);
tdrStyle->SetFrameLineWidth(1);
//For the date:
tdrStyle->SetOptDate(0);
// tdrStyle->SetDateX(Float_t x = 0.01);
// tdrStyle->SetDateY(Float_t y = 0.01);
// For the statistics box:
tdrStyle->SetOptFile(11);
tdrStyle->SetOptStat(11111111);
// For the Global title:
tdrStyle->SetOptTitle(0);
tdrStyle->SetTitleFont(42);
tdrStyle->SetTitleColor(1);
tdrStyle->SetTitleTextColor(1);
tdrStyle->SetTitleFillColor(10);
tdrStyle->SetTitleFontSize(0.05);
// For the axis titles:
tdrStyle->SetTitleColor(1, "XYZ");
tdrStyle->SetTitleFont(42, "XYZ");
tdrStyle->SetTitleSize(0.06, "XYZ");
// tdrStyle->SetTitleXSize(Float_t size = 0.02); // Another way to set the size?
// tdrStyle->SetTitleYSize(Float_t size = 0.02);
tdrStyle->SetTitleXOffset(0.7);
tdrStyle->SetTitleYOffset(0.7);
// tdrStyle->SetTitleOffset(1.1, "Y"); // Another way to set the Offset
// For the axis labels:
tdrStyle->SetLabelColor(1, "XYZ");
tdrStyle->SetLabelFont(42, "XYZ");
tdrStyle->SetLabelOffset(0.007, "XYZ");
tdrStyle->SetLabelSize(0.03, "XYZ");
// For the axis:
tdrStyle->SetAxisColor(1, "XYZ");
tdrStyle->SetStripDecimals(kTRUE);
tdrStyle->SetTickLength(0.03, "XYZ");
tdrStyle->SetNdivisions(510, "XYZ");
tdrStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame
tdrStyle->SetPadTickY(1);
tdrStyle->cd();
}
void setTDRStyle() {
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetCanvasColor(kWhite);
tdrStyle->SetCanvasDefH(600);
tdrStyle->SetCanvasDefW(600);
tdrStyle->SetCanvasDefX(0);
tdrStyle->SetCanvasDefY(0);
tdrStyle->SetPadBorderMode(0);
tdrStyle->SetPadColor(kWhite);
tdrStyle->SetFrameBorderMode(0);
tdrStyle->SetFrameBorderSize(1);
tdrStyle->SetFrameFillColor(0);
tdrStyle->SetFrameFillStyle(0);
tdrStyle->SetFrameLineColor(1);
tdrStyle->SetFrameLineStyle(1);
tdrStyle->SetFrameLineWidth(1);
tdrStyle->SetHistLineColor(1);
tdrStyle->SetHistLineStyle(0);
tdrStyle->SetHistLineWidth(1);
tdrStyle->SetEndErrorSize(2);
tdrStyle->SetErrorX(0.);
tdrStyle->SetMarkerStyle(20);
tdrStyle->SetMarkerSize(2);
tdrStyle->SetOptFit(1);
tdrStyle->SetFitFormat("5.4g");
tdrStyle->SetFuncColor(2);
tdrStyle->SetFuncStyle(1);
tdrStyle->SetFuncWidth(1);
tdrStyle->SetOptDate(0);
tdrStyle->SetOptFile(0);
tdrStyle->SetOptStat(0);
tdrStyle->SetStatColor(kWhite);
tdrStyle->SetStatFont(42);
tdrStyle->SetStatFontSize(0.025);
tdrStyle->SetStatTextColor(1);
tdrStyle->SetStatFormat("6.4g");
tdrStyle->SetStatBorderSize(1);
tdrStyle->SetStatH(0.1);
tdrStyle->SetStatW(0.15);
tdrStyle->SetPadTopMargin(0.05);
tdrStyle->SetPadBottomMargin(0.13);
tdrStyle->SetPadLeftMargin(0.10);
tdrStyle->SetPadRightMargin(0.05);
tdrStyle->SetOptTitle(0);
tdrStyle->SetTitleFont(42);
tdrStyle->SetTitleColor(1);
tdrStyle->SetTitleTextColor(1);
tdrStyle->SetTitleFillColor(10);
tdrStyle->SetTitleFontSize(0.05);
// For the axis labels:
tdrStyle->SetLabelColor(1, "XYZ");
tdrStyle->SetLabelOffset(0.007, "XYZ");
tdrStyle->SetLabelSize(0.03, "XYZ");
tdrStyle->SetAxisColor(1, "XYZ");
tdrStyle->SetStripDecimals(kTRUE);
tdrStyle->SetTickLength(0.03, "XYZ");
tdrStyle->SetNdivisions(510, "XYZ");
tdrStyle->SetPadTickX(1);
tdrStyle->SetPadTickY(1);
tdrStyle->SetOptLogx(0);
tdrStyle->SetOptLogy(0);
tdrStyle->SetOptLogz(0);
tdrStyle->SetPaperSize(20.,20.);
tdrStyle->cd();
}
void Style_Toby_2()
{
// Add the saved style to the current ROOT session.
delete gROOT->GetStyle("Toby_2");
TStyle *tmpStyle = new TStyle("Toby_2", "Style for posters and talks");
tmpStyle->SetNdivisions(510, "x");
tmpStyle->SetNdivisions(510, "y");
tmpStyle->SetNdivisions(510, "z");
tmpStyle->SetAxisColor(1, "x");
tmpStyle->SetAxisColor(1, "y");
tmpStyle->SetAxisColor(1, "z");
tmpStyle->SetLabelColor(1, "x");
tmpStyle->SetLabelColor(1, "y");
tmpStyle->SetLabelColor(1, "z");
tmpStyle->SetLabelFont(42, "x");
tmpStyle->SetLabelFont(42, "y");
tmpStyle->SetLabelFont(42, "z");
tmpStyle->SetLabelOffset(0.005, "x");
tmpStyle->SetLabelOffset(0.005, "y");
tmpStyle->SetLabelOffset(0.005, "z");
tmpStyle->SetLabelSize(0.035, "x");
tmpStyle->SetLabelSize(0.035, "y");
tmpStyle->SetLabelSize(0.035, "z");
tmpStyle->SetTickLength(0.03, "x");
tmpStyle->SetTickLength(0.03, "y");
tmpStyle->SetTickLength(0.03, "z");
tmpStyle->SetTitleOffset(1, "x");
tmpStyle->SetTitleOffset(1, "y");
tmpStyle->SetTitleOffset(1, "z");
tmpStyle->SetTitleSize(0.035, "x");
tmpStyle->SetTitleSize(0.035, "y");
tmpStyle->SetTitleSize(0.035, "z");
tmpStyle->SetTitleColor(1, "x");
tmpStyle->SetTitleColor(1, "y");
tmpStyle->SetTitleColor(1, "z");
tmpStyle->SetTitleFont(42, "x");
tmpStyle->SetTitleFont(42, "y");
tmpStyle->SetTitleFont(42, "z");
tmpStyle->SetBarWidth(1);
tmpStyle->SetBarOffset(0);
tmpStyle->SetDrawBorder(0);
tmpStyle->SetOptLogx(0);
tmpStyle->SetOptLogy(0);
tmpStyle->SetOptLogz(0);
tmpStyle->SetOptDate(0);
tmpStyle->SetOptStat(111);
tmpStyle->SetOptTitle(kTRUE);
tmpStyle->SetOptFit(0);
tmpStyle->SetNumberContours(20);
tmpStyle->GetAttDate()->SetTextFont(62);
tmpStyle->GetAttDate()->SetTextSize(0.025);
tmpStyle->GetAttDate()->SetTextAngle(0);
tmpStyle->GetAttDate()->SetTextAlign(11);
tmpStyle->GetAttDate()->SetTextColor(1);
tmpStyle->SetDateX(0.01);
tmpStyle->SetDateY(0.01);
tmpStyle->SetEndErrorSize(2);
tmpStyle->SetErrorX(0.5);
tmpStyle->SetFuncColor(2);
tmpStyle->SetFuncStyle(1);
tmpStyle->SetFuncWidth(2);
tmpStyle->SetGridColor(0);
tmpStyle->SetGridStyle(3);
tmpStyle->SetGridWidth(1);
tmpStyle->SetLegendBorderSize(1);
tmpStyle->SetLegendFillColor(0);
tmpStyle->SetLegendFont(42);
tmpStyle->SetHatchesLineWidth(1);
tmpStyle->SetHatchesSpacing(1);
tmpStyle->SetFrameFillColor(0);
tmpStyle->SetFrameLineColor(1);
tmpStyle->SetFrameFillStyle(1001);
tmpStyle->SetFrameLineStyle(1);
tmpStyle->SetFrameLineWidth(1);
tmpStyle->SetFrameBorderSize(1);
tmpStyle->SetFrameBorderMode(0);
tmpStyle->SetHistFillColor(0);
tmpStyle->SetHistLineColor(602);
tmpStyle->SetHistFillStyle(1001);
tmpStyle->SetHistLineStyle(1);
tmpStyle->SetHistLineWidth(1);
tmpStyle->SetHistMinimumZero(kFALSE);
tmpStyle->SetCanvasPreferGL(kFALSE);
tmpStyle->SetCanvasColor(0);
tmpStyle->SetCanvasBorderSize(2);
tmpStyle->SetCanvasBorderMode(0);
tmpStyle->SetCanvasDefH(500);
tmpStyle->SetCanvasDefW(700);
tmpStyle->SetCanvasDefX(10);
tmpStyle->SetCanvasDefY(10);
tmpStyle->SetPadColor(0);
tmpStyle->SetPadBorderSize(2);
tmpStyle->SetPadBorderMode(0);
tmpStyle->SetPadBottomMargin(0.1);
tmpStyle->SetPadTopMargin(0.1);
tmpStyle->SetPadLeftMargin(0.1);
tmpStyle->SetPadRightMargin(0.1);
tmpStyle->SetPadGridX(kFALSE);
tmpStyle->SetPadGridY(kFALSE);
tmpStyle->SetPadTickX(0);
tmpStyle->SetPadTickY(0);
tmpStyle->SetPaperSize(20, 26);
tmpStyle->SetScreenFactor(1);
tmpStyle->SetStatColor(0);
tmpStyle->SetStatTextColor(1);
tmpStyle->SetStatBorderSize(1);
tmpStyle->SetStatFont(42);
tmpStyle->SetStatFontSize(0);
tmpStyle->SetStatStyle(1001);
tmpStyle->SetStatFormat("6.4g");
tmpStyle->SetStatX(0.9);
tmpStyle->SetStatY(0.97);
tmpStyle->SetStatW(0.15);
tmpStyle->SetStatH(0.16);
tmpStyle->SetStripDecimals(kTRUE);
tmpStyle->SetTitleAlign(23);
tmpStyle->SetTitleFillColor(0);
tmpStyle->SetTitleTextColor(1);
tmpStyle->SetTitleBorderSize(0);
tmpStyle->SetTitleFont(42);
tmpStyle->SetTitleFontSize(0.05);
tmpStyle->SetTitleStyle(0);
tmpStyle->SetTitleX(0.5);
tmpStyle->SetTitleY(0.995);
tmpStyle->SetTitleW(0);
tmpStyle->SetTitleH(0);
tmpStyle->SetLegoInnerR(0.5);
Int_t fPaletteColor[50] = {51, 52, 53, 54, 55, 56, 57, 58, 59,
60, 61, 62, 63, 64, 65, 66, 67, 68, 69,
70, 71, 72, 73, 74, 75, 76, 77, 78, 79,
80, 81, 82, 83, 84, 85, 86, 87, 88, 89,
90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100};
tmpStyle->SetPalette(50, fPaletteColor);
TString fLineStyleArrayTmp[30] = {"", " ", " 12 12", " 4 8",
" 12 16 4 16", " 20 12 4 12", " 20 12 4 12 4 12 4 12", " 20 20", " 20 12 4 12 4 12",
" 80 20", " 80 40 4 40", " ", " ", " ",
" ", " ", " ", " ", " ",
" ", " ", " ", " ", " ",
" ", " ", " ", " ", " ", " "};
for (Int_t i=0; i<30; i++)
tmpStyle->SetLineStyleString(i, fLineStyleArrayTmp[i]);
tmpStyle->SetHeaderPS("");
tmpStyle->SetTitlePS("");
tmpStyle->SetFitFormat("5.4g");
tmpStyle->SetPaintTextFormat("g");
tmpStyle->SetLineScalePS(3);
tmpStyle->SetColorModelPS(0);
tmpStyle->SetTimeOffset(788918400);
tmpStyle->SetLineColor(1);
tmpStyle->SetLineStyle(1);
tmpStyle->SetLineWidth(1);
tmpStyle->SetFillColor(19);
tmpStyle->SetFillStyle(1001);
tmpStyle->SetMarkerColor(1);
tmpStyle->SetMarkerSize(1);
tmpStyle->SetMarkerStyle(1);
tmpStyle->SetTextAlign(11);
tmpStyle->SetTextAngle(0);
tmpStyle->SetTextColor(1);
tmpStyle->SetTextFont(62);
tmpStyle->SetTextSize(0.05);
}
void blinding_study()
{
gSystem->CompileMacro("MitGPTree.h");
// First we define MIT Style for the plots.
TStyle *MitStyle = gStyle;
//gStyle = MitStyle;
// Canvas
MitStyle->SetCanvasColor (0);
MitStyle->SetCanvasBorderSize(10);
MitStyle->SetCanvasBorderMode(0);
MitStyle->SetCanvasDefH (700);
MitStyle->SetCanvasDefW (700);
MitStyle->SetCanvasDefX (100);
MitStyle->SetCanvasDefY (100);
// Pads
MitStyle->SetPadColor (0);
MitStyle->SetPadBorderSize (10);
MitStyle->SetPadBorderMode (0);
MitStyle->SetPadBottomMargin(0.13);
MitStyle->SetPadTopMargin (0.04);
MitStyle->SetPadLeftMargin (0.18);
MitStyle->SetPadRightMargin (0.04);
MitStyle->SetPadGridX (0);
MitStyle->SetPadGridY (0);
MitStyle->SetPadTickX (0);
MitStyle->SetPadTickY (0);
// Frames
MitStyle->SetFrameFillStyle ( 0);
MitStyle->SetFrameFillColor ( 0);
MitStyle->SetFrameLineColor ( 1);
MitStyle->SetFrameLineStyle ( 0);
MitStyle->SetFrameLineWidth ( 1);
MitStyle->SetFrameBorderSize(10);
MitStyle->SetFrameBorderMode( 0);
// Histograms
MitStyle->SetHistFillColor(2);
MitStyle->SetHistFillStyle(0);
MitStyle->SetHistLineColor(1);
MitStyle->SetHistLineStyle(0);
MitStyle->SetHistLineWidth(2);
MitStyle->SetNdivisions(505);
// Functions
MitStyle->SetFuncColor(1);
MitStyle->SetFuncStyle(0);
MitStyle->SetFuncWidth(2);
// Various
MitStyle->SetMarkerStyle(20);
MitStyle->SetMarkerColor(kBlack);
MitStyle->SetMarkerSize (1.2);
MitStyle->SetTitleSize (0.055,"X");
MitStyle->SetTitleOffset(1.200,"X");
MitStyle->SetLabelOffset(0.005,"X");
MitStyle->SetLabelSize (0.050,"X");
MitStyle->SetLabelFont (42 ,"X");
MitStyle->SetTickLength (-0.03,"X");
MitStyle->SetStripDecimals(kFALSE);
MitStyle->SetTitleSize (0.055,"Y");
MitStyle->SetTitleOffset(1.800,"Y");
MitStyle->SetLabelOffset(0.010,"Y");
MitStyle->SetLabelSize (0.050,"Y");
MitStyle->SetLabelFont (42 ,"Y");
MitStyle->SetTickLength (-0.03,"Y");
MitStyle->SetTextSize (0.055);
MitStyle->SetTextFont (42);
MitStyle->SetStatFont (42);
MitStyle->SetTitleFont (42);
MitStyle->SetTitleFont (42,"X");
MitStyle->SetTitleFont (42,"Y");
MitStyle->SetOptStat (0);
// Here the style section ends and the macro begins.
string sig_samples[] =
{
"s12-dmmpho-v_m1-v7a",
"s12-dmmpho-av_m1-v7a",
"s12-dmmpho-v_m10-v7a",
"s12-dmmpho-av_m10-v7a",
"s12-dmmpho-v_m100-v7a",
"s12-dmmpho-av_m100-v7a",
"s12-dmmpho-v_m200-v7a",
"s12-dmmpho-av_m200-v7a",
"s12-dmmpho-av_m300-v7a",
"s12-dmmpho-v_m500-v7a",
"s12-dmmpho-av_m500-v7a",
"s12-dmmpho-v_m1000-v7a",
"s12-dmmpho-av_m1000-v7a",
"s12-addmpho-md1_d2-v7a",
"s12-addmpho-md1_d3-v7a",
"s12-addmpho-md1_d4-v7a",
"s12-addmpho-md1_d5-v7a",
"s12-addmpho-md1_d6-v7a",
"s12-addmpho-md2_d2-v7a",
"s12-addmpho-md2_d3-v7a",
"s12-addmpho-md2_d5-v7a",
"s12-addmpho-md2_d6-v7a",
"s12-addmpho-md3_d2-v7a",
"s12-addmpho-md3_d3-v7a",
"s12-addmpho-md3_d4-v7a",
"s12-addmpho-md3_d5-v7a",
"s12-addmpho-md3_d6-v7a"
}
Int_t signal_num=0;
// This macro considers only one signal sample at a time, so the variable signal_num defines
// which signal sample to work on (it is the index of the sample in sig_samples).
double sig_weights[] =
{
4.81E-07,
4.79E-07,
4.80E-07,
4.81E-07,
4.77E-07,
4.26E-07,
4.23E-07,
3.18E-07,
2.19E-07,
2.01E-07,
9.52E-08,
3.00E-08,
8.22E-09,
6.48E-02,
1.73E-01,
6.93E-02,
9.17E-02,
9.63E-01,
5.87E-03,
4.91E-03,
4.37E-03,
4.26E-03,
1.45E-03,
7.93E-04,
5.53E-04,
4.26E-04,
3.24E-04
};
// The signal weights are given by (sigma_MC * lumi)/(N_processed). Background weights are defined similarly below.
string bg_samples[] = {"s12-zgptg130-v7c", "s12-wjets-ptw100-v7a", "s12-wgptg130-v7a", "s12-qcdht100-250-v7a", "s12-qcdht250-500-v7a", "s12-qcdht500-1000-v7a",
"s12-qcdht1000-v7a", "s12-2pibo10_25-v7a", "s12-2pibo25_250-v7a", "s12-2pibo250-v7a", "s12-2pibx10_25-v7a", "s12-2pibx25_250-v7a", "s12-2pibx250-v7a",
"s12-pj50_80-v7a", "s12-pj80_120-v7a", "s12-pj120_170-v7a", "s12-pj170_300-v7a", "s12-pj300_470-v7a", "s12-pj470_800-v7a", "s12-pj800_1400-v7a",
"s12-pj1400_1800-v7a", "s12-pj1800-v7a", "s12-zgllgptg130-v7a", "s12-zgllg-v7a"};
double bg_weights[]=
{
5.28E-03,
1.14E-01,
1.38E-02,
5.63E+03,
2.34E+02,
9.71E+00,
3.09E-01,
9.21E+00,
9.90E-01,
4.21E-04,
1.66E+01,
6.43E-01,
4.91E-05,
3.30E+01,
5.49E+00,
1.09E+00,
3.01E-01,
2.15E-02,
2.10E-03,
7.11E-05,
4.45E-07,
1.88E-08,
3.18E-03,
5.10E-01
};
string sigLine = sig_samples[signal_num];
Double_t sigWeight = sig_weights[signal_num];
TH1* sig_hist = new TH1F("sig_hist", TString("Phi Between Met and Photon for Signal Events [") + TString(sigLine) + TString("]"), 70, 0, 3.5);
TH1* bg_hist = new TH1F("bg_hist", TString("Phi Between Met and Photon for Background Events [") + TString(sigLine) + TString("]"), 70, 0, 3.5);
TH1* signifhist = new TH1F("signifhist", TString(" "), 24, 2, 3.2);
cout<<"The selected signal sample is: "<<TString(sigLine)<<" "<<endl;
TString sigFilename = TString("monoph-2013-July9_") + TString(sigLine) + TString("_noskim.root");
MitGPTree sigEvent;
sigEvent.LoadTree(TString("/scratch/cferko/hist/monoph-2013-July9/merged/")+sigFilename, 0);
sigEvent.InitTree(0);
int nDataSig=sigEvent.tree_->GetEntries();
for (int evt=0; evt<nDataSig; ++evt)
{
sigEvent.tree_->GetEntry(evt);
Double_t pho1Pt = sigEvent.pho1_.Pt();
Double_t jet1Pt = sigEvent.jet1_.Pt();
Double_t met = sigEvent.met_;
Double_t metPhi = sigEvent.metPhi_;
Double_t nphotons = sigEvent.nphotons_;
Double_t ncosmics = sigEvent.ncosmics_;
Double_t phoPassEleVeto = sigEvent.phoPassEleVeto_a1_;
Double_t jet1Pt = sigEvent.jet1_.Pt();
Double_t nlep = sigEvent.nlep_;
Double_t pho1Eta = sigEvent.pho1_.Eta();
Double_t pho1Phi = sigEvent.pho1_.Phi();
Double_t phoIsTrigger = sigEvent.phoIsTrigger_a1_;
Double_t phoLeadTimeSpan = sigEvent.phoLeadTimeSpan_a1_;
Double_t phoCoviEtaiEta = sigEvent.phoCoviEtaiEta_a1_;
Double_t phoCoviPhiiPhi = sigEvent.phoCoviPhiiPhi_a1_;
Double_t phoMipIsHalo = sigEvent.phoMipIsHalo_a1_;
Double_t phoSeedTime = sigEvent.phoSeedTime_a1_;
Double_t deltaPhi = TMath::ACos(TMath::Cos(pho1Phi - metPhi));
if (TMath::Abs(pho1Eta)<1.479 && met>140 && pho1Pt>160 && nphotons>0 && phoPassEleVeto>0 && phoIsTrigger==1 && TMath::Abs(phoLeadTimeSpan) < 8. && phoCoviEtaiEta > 0.001 && phoCoviPhiiPhi > 0.001 && phoMipIsHalo == 0 && phoSeedTime > -1.5 && nlep==0 && ncosmics==0 && jet1Pt<100) sig_hist->Fill(deltaPhi, sigWeight);
}
for (int i=0; i<24; i++)
{
string bgLine = bg_samples[i];
Double_t bgWeight = bg_weights[i];
TString bgFilename = TString("monoph-2013-July9_") + TString(bgLine) + TString("_noskim.root");
MitGPTree bgEvent;
bgEvent.LoadTree(TString("/scratch/cferko/hist/monoph-2013-July9/merged/") + bgFilename, 0);
bgEvent.InitTree(0);
int nDataBg=bgEvent.tree_->GetEntries();
for (int evt=0; evt<nDataBg; ++evt)
{
bgEvent.tree_->GetEntry(evt);
Double_t pho1Pt = bgEvent.pho1_.Pt();
Double_t jet1Pt = bgEvent.jet1_.Pt();
Double_t met = bgEvent.met_;
Double_t metPhi = bgEvent.metPhi_;
Double_t nphotons = bgEvent.nphotons_;
Double_t ncosmics = bgEvent.ncosmics_;
Double_t phoPassEleVeto = bgEvent.phoPassEleVeto_a1_;
Double_t jet1Pt = bgEvent.jet1_.Pt();
Double_t nlep = bgEvent.nlep_;
Double_t pho1Eta = bgEvent.pho1_.Eta();
Double_t pho1Phi = bgEvent.pho1_.Phi();
Double_t phoIsTrigger = bgEvent.phoIsTrigger_a1_;
Double_t phoLeadTimeSpan = bgEvent.phoLeadTimeSpan_a1_;
Double_t phoCoviEtaiEta = bgEvent.phoCoviEtaiEta_a1_;
Double_t phoCoviPhiiPhi = bgEvent.phoCoviPhiiPhi_a1_;
Double_t phoMipIsHalo = bgEvent.phoMipIsHalo_a1_;
Double_t phoSeedTime = bgEvent.phoSeedTime_a1_;
Double_t deltaPhi = TMath::ACos(TMath::Cos(pho1Phi - metPhi));
if (TMath::Abs(pho1Eta)<1.479 && met>140 && pho1Pt>160 && nphotons>0 && phoPassEleVeto>0 && phoIsTrigger==1 && TMath::Abs(phoLeadTimeSpan) < 8. && phoCoviEtaiEta > 0.001 && phoCoviPhiiPhi > 0.001 && phoMipIsHalo == 0 && phoSeedTime > -1.5 && nlep==0 && ncosmics==0 && jet1Pt<100) bg_hist->Fill(deltaPhi, bgWeight);
}
}
Double_t sigTot = sig_hist->Integral();
Double_t bgTot = bg_hist->Integral();
Double_t signifTot = sigTot/TMath::Sqrt(bgTot);
Double_t signifTarget = 0.25*signifTot;
Int_t reduced=0;
// The strange limits in the for loop are a messy hack to get the integral to work out. Since hist->Integral(a,b)
// integrates from BIN a to BIN b (bin numbers instead of values of the x-axis), I choose to work from bin 64
// (an angle cut of 3.2) to bin 40 (an angle cut of 2).
for (int i=64; i>=40; i--)
{
Double_t sigcount = sig_hist->Integral(0, i);
Double_t bgcount = bg_hist->Integral(0, i);
if (bgcount>0)
{
Double_t significance = sigcount/TMath::Sqrt(bgcount);
// cout<<"At a phi cut of "<<sig_hist->GetBinCenter(i)+0.025<<" the significance is "<<significance<<", which is a fraction "<<significance/signifTot<<" of max."<<endl;
if (significance<signifTarget && reduced==0)
{
cout<<"The signal significance is reduced to one-fourth of its uncut value at a phi cut of "<<sig_hist->GetBinCenter(i)+0.025<<endl;
reduced=1;
}
signifhist->Fill(sig_hist->GetBinCenter(i), significance);
}
}
signifhist->SetStats(kFALSE);
signifhist->GetXaxis()->SetTitle("Maximum #Delta#phi");
signifhist->GetYaxis()->SetTitle("Signal Significance");
signifhist->Draw();
c1->Print(TString("Blinding_") + TString(sigLine) + TString(".pdf"));
}
void SetStyle()
{
TStyle *MYStyle = new TStyle("Style","Style");
gStyle = MYStyle;
// Canvas
MYStyle->SetCanvasColor (0);
MYStyle->SetCanvasBorderSize(10);
MYStyle->SetCanvasBorderMode(0);
MYStyle->SetCanvasDefH (700);
MYStyle->SetCanvasDefW (700);
MYStyle->SetCanvasDefX (100);
MYStyle->SetCanvasDefY (100);
// Pads
MYStyle->SetPadColor (0);
MYStyle->SetPadBorderSize (10);
MYStyle->SetPadBorderMode (0);
MYStyle->SetPadBottomMargin(0.13);
MYStyle->SetPadTopMargin (0.08);
MYStyle->SetPadLeftMargin (0.15);
MYStyle->SetPadRightMargin (0.05);
MYStyle->SetPadGridX (0);
MYStyle->SetPadGridY (0);
MYStyle->SetPadTickX (0);
MYStyle->SetPadTickY (0);
// Frames
MYStyle->SetFrameFillStyle ( 0);
MYStyle->SetFrameFillColor ( 0);
MYStyle->SetFrameLineColor ( 1);
MYStyle->SetFrameLineStyle ( 0);
MYStyle->SetFrameLineWidth ( 1);
MYStyle->SetFrameBorderSize(10);
MYStyle->SetFrameBorderMode( 0);
// Histograms
MYStyle->SetHistFillColor(2);
// MYStyle->SetHistFillStyle(0);
MYStyle->SetHistLineColor(1);
MYStyle->SetHistLineStyle(0);
MYStyle->SetHistLineWidth(2);
MYStyle->SetNdivisions(505);
// Functions
MYStyle->SetFuncColor(1);
MYStyle->SetFuncStyle(0);
MYStyle->SetFuncWidth(2);
// Various
MYStyle->SetMarkerStyle(20);
MYStyle->SetMarkerColor(kBlack);
MYStyle->SetMarkerSize (1.2);
MYStyle->SetTitleSize (0.055,"X");
MYStyle->SetTitleOffset(1.050,"X");
MYStyle->SetLabelOffset(0.005,"X");
MYStyle->SetLabelSize (0.045,"X");
MYStyle->SetLabelFont (42 ,"X");
MYStyle->SetStripDecimals(kFALSE);
MYStyle->SetTitleSize (0.055,"Y");
MYStyle->SetTitleOffset(1.400,"Y");
MYStyle->SetLabelOffset(0.010,"Y");
MYStyle->SetLabelSize (0.045,"Y");
MYStyle->SetLabelFont (42 ,"Y");
MYStyle->SetTextSize (0.055);
MYStyle->SetTextFont (42);
MYStyle->SetStatFont (42);
MYStyle->SetTitleFont (42);
MYStyle->SetTitleFont (42,"X");
MYStyle->SetTitleFont (42,"Y");
MYStyle->SetOptStat (0);
MYStyle->SetPalette (1);
return;
}
void setTDRStyle() {
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-DR");
gROOT->SetStyle("Plain");
gROOT->ForceStyle();
tdrStyle->SetFrameBorderMode(0);
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetPadBorderMode(0);
tdrStyle->SetPadBorderMode(0);
// tdrStyle->SetFrameColor(0);
tdrStyle->SetPadColor(0);
tdrStyle->SetCanvasColor(0);
tdrStyle->SetStatColor(0);
tdrStyle->SetFillColor(0);
tdrStyle->SetPaperSize(20,26);
// tdrStyle->SetPadTopMargin(0.08);
// tdrStyle->SetPadBottomMargin(0.14);
tdrStyle->SetPadRightMargin(0.04);
tdrStyle->SetPadLeftMargin(0.05);
// tdrStyle->SetCanvasDefH(800);
// tdrStyle->SetCanvasDefW(800);
// tdrStyle->SetPadGridX(1);
// tdrStyle->SetPadGridY(1);
tdrStyle->SetPadTickX(1);
tdrStyle->SetPadTickY(1);
tdrStyle->SetTextFont(42); //132
tdrStyle->SetTextSize(0.09);
tdrStyle->SetLabelFont(42,"xyz");
tdrStyle->SetTitleFont(42,"xyz");
tdrStyle->SetLabelSize(0.045,"xyz"); //0.035
tdrStyle->SetTitleSize(0.045,"xyz");
tdrStyle->SetTitleOffset(1.5,"y");
tdrStyle->SetTitleX(0.16);
tdrStyle->SetTitleY(0.93);
tdrStyle->SetTitleColor(1);
tdrStyle->SetTitleTextColor(1);
tdrStyle->SetTitleFillColor(0);
tdrStyle->SetTitleBorderSize(1);
tdrStyle->SetTitleFontSize(0.04);
// tdrStyle->SetPadTopMargin(0.05);
// tdrStyle->SetPadBottomMargin(0.13);
// tdrStyle->SetPadLeftMargin(0.14);
// tdrStyle->SetPadRightMargin(0.02);
// use bold lines and markers
tdrStyle->SetMarkerStyle(8);
tdrStyle->SetHistLineWidth(3);
tdrStyle->SetLineWidth(1);
tdrStyle->SetOptTitle(kFALSE);
tdrStyle->SetOptStat(0);
tdrStyle->cd();
}
void SetTMVAStyle() {
TStyle *TMVAStyle = gROOT->GetStyle("TMVA");
if(TMVAStyle!=0) {
gROOT->SetStyle("TMVA");
return;
}
TMVAStyle = new TStyle(*gROOT->GetStyle("Plain")); // our style is based on Plain
TMVAStyle->SetName("TMVA");
TMVAStyle->SetTitle("TMVA style based on \"Plain\" with modifications defined in tmvaglob.C");
gROOT->GetListOfStyles()->Add(TMVAStyle);
gROOT->SetStyle("TMVA");
TMVAStyle->SetLineStyleString( 5, "[52 12]" );
TMVAStyle->SetLineStyleString( 6, "[22 12]" );
TMVAStyle->SetLineStyleString( 7, "[22 10 7 10]" );
// the pretty color palette of old
TMVAStyle->SetPalette((UsePaperStyle ? 18 : 1),0);
// use plain black on white colors
TMVAStyle->SetFrameBorderMode(0);
TMVAStyle->SetCanvasBorderMode(0);
TMVAStyle->SetPadBorderMode(0);
TMVAStyle->SetPadColor(0);
TMVAStyle->SetFillStyle(0);
TMVAStyle->SetLegendBorderSize(0);
// title properties
// TMVAStyle->SetTitleW(.4);
// TMVAStyle->SetTitleH(.10);
// MVAStyle->SetTitleX(.5);
// TMVAStyle->SetTitleY(.9);
TMVAStyle->SetTitleFillColor( c_TitleBox );
TMVAStyle->SetTitleTextColor( c_TitleText );
TMVAStyle->SetTitleBorderSize( 1 );
TMVAStyle->SetLineColor( c_TitleBorder );
if (!UsePaperStyle) {
TMVAStyle->SetFrameFillColor( c_FrameFill );
TMVAStyle->SetCanvasColor( c_Canvas );
}
// set the paper & margin sizes
TMVAStyle->SetPaperSize(20,26);
TMVAStyle->SetPadTopMargin(0.10);
TMVAStyle->SetPadRightMargin(0.05);
TMVAStyle->SetPadBottomMargin(0.11);
TMVAStyle->SetPadLeftMargin(0.12);
// use bold lines and markers
TMVAStyle->SetMarkerStyle(21);
TMVAStyle->SetMarkerSize(0.3);
TMVAStyle->SetHistLineWidth(2);
TMVAStyle->SetLineStyleString(2,"[12 12]"); // postscript dashes
// do not display any of the standard histogram decorations
TMVAStyle->SetOptTitle(1);
TMVAStyle->SetTitleH(0.052);
TMVAStyle->SetOptStat(0);
TMVAStyle->SetOptFit(0);
// put tick marks on top and RHS of plots
TMVAStyle->SetPadTickX(1);
TMVAStyle->SetPadTickY(1);
}
void set_plot_style() {
TStyle* style = new TStyle("style", "Style for P-TDR");
// For the canvas:
style->SetCanvasBorderMode(0);
style->SetCanvasColor(kWhite);
style->SetCanvasDefX(0); //Position on screen
style->SetCanvasDefY(0);
// For the Pad:
style->SetPadBorderMode(0);
style->SetPadColor(kWhite);
style->SetPadGridX(false);
style->SetPadGridY(false);
style->SetGridColor(kBlack);
style->SetGridStyle(3);
style->SetGridWidth(1);
// For the frame:
style->SetFrameBorderMode(0);
style->SetFrameBorderSize(1);
style->SetFrameFillColor(kWhite);
style->SetFrameFillStyle(0);
style->SetFrameLineColor(kBlack);
style->SetFrameLineStyle(1);
style->SetFrameLineWidth(1);
// For the histo:
// style->SetHistFillColor(1);
style->SetHistFillStyle(0); //
style->SetHistLineColor(kBlack);
style->SetHistLineStyle(0);
style->SetHistLineWidth(1);
style->SetEndErrorSize(2);
style->SetErrorX(0.);
style->SetMarkerStyle(20);
//For the fit/function:
style->SetOptFit(1);
style->SetFitFormat("5.4g");
style->SetFuncColor(kRed);
style->SetFuncStyle(1);
style->SetFuncWidth(1);
//For the date:
style->SetOptDate(0);
// For the statistics box:
style->SetOptFile(0);
style->SetOptStat(0); // To display the mean and RMS: SetOptStat("mr");
style->SetStatColor(kWhite);
style->SetStatFont(42);
style->SetStatFontSize(0.025);
style->SetStatTextColor(kBlack);
style->SetStatFormat("6.4g");
style->SetStatBorderSize(1);
style->SetStatH(0.1);
style->SetStatW(0.15);
// Margins:
style->SetPadTopMargin(1 - 0.95);
style->SetPadBottomMargin(0.1);
style->SetPadLeftMargin(0.12);
style->SetPadRightMargin(1 - 0.96);
// For the Global title:
style->SetOptTitle(0);
style->SetTitleFont(42);
style->SetTitleColor(kBlack);
style->SetTitleTextColor(kBlack);
style->SetTitleFillColor(kWhite); //10 is roughly kWhite, 10% grey?
style->SetTitleFontSize(0.05);
// For the axis titles:
style->SetTitleColor(kBlack, "XYZ");
style->SetTitleFont(42, "XYZ");
style->SetTitleSize(0.06, "XYZ");
style->SetTitleXOffset(0.9);
style->SetTitleYOffset(1.25);
// For the axis labels:
style->SetLabelColor(kBlack, "XYZ");
style->SetLabelFont(42, "XYZ");
style->SetLabelOffset(0.007, "XYZ");
style->SetLabelSize(0.05, "XYZ");
// For the axis:
style->SetAxisColor(kBlack, "XYZ");
style->SetStripDecimals(true);
style->SetTickLength(0.03, "XYZ");
style->SetNdivisions(510, "XYZ");
style->SetPadTickX(true); // To get tick marks on the opposite side of the frame
style->SetPadTickY(true);
// Change for log plots:
style->SetOptLogx(false);
style->SetOptLogy(false);
style->SetOptLogz(false);
// Set the style
style->cd();
}
void setTDRStyle() {
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");
// For the canvas:
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetCanvasColor(kWhite);
tdrStyle->SetCanvasDefH(600); //Height of canvas
tdrStyle->SetCanvasDefW(600); //Width of canvas
tdrStyle->SetCanvasDefX(0); //POsition on screen
tdrStyle->SetCanvasDefY(0);
// For the Pad:
tdrStyle->SetPadBorderMode(0);
// tdrStyle->SetPadBorderSize(Width_t size = 1);
tdrStyle->SetPadColor(kWhite);
tdrStyle->SetPadGridX(false);
tdrStyle->SetPadGridY(false);
tdrStyle->SetGridColor(0);
tdrStyle->SetGridStyle(3);
tdrStyle->SetGridWidth(1);
// For the frame:
tdrStyle->SetFrameBorderMode(0);
tdrStyle->SetFrameBorderSize(1);
tdrStyle->SetFrameFillColor(0);
tdrStyle->SetFrameFillStyle(0);
tdrStyle->SetFrameLineColor(1);
tdrStyle->SetFrameLineStyle(1);
tdrStyle->SetFrameLineWidth(1);
// For the histo:
// tdrStyle->SetHistFillColor(1);
// tdrStyle->SetHistFillStyle(0);
tdrStyle->SetHistLineColor(1);
tdrStyle->SetHistLineStyle(0);
tdrStyle->SetHistLineWidth(1);
// tdrStyle->SetLegoInnerR(Float_t rad = 0.5);
// tdrStyle->SetNumberContours(Int_t number = 20);
tdrStyle->SetEndErrorSize(2);
// tdrStyle->SetErrorMarker(20);
// tdrStyle->SetErrorX(0.);
tdrStyle->SetMarkerStyle(20);
//For the fit/function:
tdrStyle->SetOptFit(1);
tdrStyle->SetFitFormat("5.4g");
tdrStyle->SetFuncColor(2);
tdrStyle->SetFuncStyle(1);
tdrStyle->SetFuncWidth(1);
//For the date:
tdrStyle->SetOptDate(0);
// tdrStyle->SetDateX(Float_t x = 0.01);
// tdrStyle->SetDateY(Float_t y = 0.01);
// For the statistics box:
tdrStyle->SetOptFile(0);
tdrStyle->SetOptStat(0); // To display the mean and RMS: SetOptStat("mr");
tdrStyle->SetStatColor(kWhite);
tdrStyle->SetStatFont(42);
tdrStyle->SetStatFontSize(0.025);
tdrStyle->SetStatTextColor(1);
tdrStyle->SetStatFormat("6.4g");
tdrStyle->SetStatBorderSize(1);
tdrStyle->SetStatH(0.1);
tdrStyle->SetStatW(0.15);
// tdrStyle->SetStatStyle(Style_t style = 1001);
// tdrStyle->SetStatX(Float_t x = 0);
// tdrStyle->SetStatY(Float_t y = 0);
// Margins:
tdrStyle->SetPadTopMargin(0.05);
tdrStyle->SetPadBottomMargin(0.13);
tdrStyle->SetPadLeftMargin(0.16);
tdrStyle->SetPadRightMargin(0.02);
// For the Global title:
tdrStyle->SetOptTitle(0);
tdrStyle->SetTitleFont(42);
tdrStyle->SetTitleColor(1);
tdrStyle->SetTitleTextColor(1);
tdrStyle->SetTitleFillColor(10);
tdrStyle->SetTitleFontSize(0.05);
// tdrStyle->SetTitleH(0); // Set the height of the title box
// tdrStyle->SetTitleW(0); // Set the width of the title box
// tdrStyle->SetTitleX(0); // Set the position of the title box
// tdrStyle->SetTitleY(0.985); // Set the position of the title box
// tdrStyle->SetTitleStyle(Style_t style = 1001);
// tdrStyle->SetTitleBorderSize(2);
// For the axis titles:
tdrStyle->SetTitleColor(1, "XYZ");
tdrStyle->SetTitleFont(42, "XYZ");
tdrStyle->SetTitleSize(0.06, "XYZ");
// tdrStyle->SetTitleXSize(Float_t size = 0.02); // Another way to set the size?
// tdrStyle->SetTitleYSize(Float_t size = 0.02);
tdrStyle->SetTitleXOffset(0.9);
tdrStyle->SetTitleYOffset(1.4);
// tdrStyle->SetTitleOffset(1.1, "Y"); // Another way to set the Offset
// For the axis labels:
tdrStyle->SetLabelColor(1, "XYZ");
tdrStyle->SetLabelFont(42, "XYZ");
tdrStyle->SetTitleSize(0.06, "XYZ");
// tdrStyle->SetTitleXSize(Float_t size = 0.02); // Another way to set the size?
// tdrStyle->SetTitleYSize(Float_t size = 0.02);
tdrStyle->SetTitleXOffset(0.9);
tdrStyle->SetTitleYOffset(1.4);
// tdrStyle->SetTitleOffset(1.1, "Y"); // Another way to set the Offset
// For the axis labels:
tdrStyle->SetLabelColor(1, "XYZ");
tdrStyle->SetLabelFont(42, "XYZ");
tdrStyle->SetLabelOffset(0.007, "XYZ");
tdrStyle->SetLabelSize(0.05, "XYZ");
// For the axis:
tdrStyle->SetAxisColor(1, "XYZ");
tdrStyle->SetStripDecimals(kTRUE);
tdrStyle->SetTickLength(0.03, "XYZ");
tdrStyle->SetNdivisions(510, "XYZ");
tdrStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame
tdrStyle->SetPadTickY(1);
// Change for log plots:
tdrStyle->SetOptLogx(0);
tdrStyle->SetOptLogy(0);
tdrStyle->SetOptLogz(0);
// Postscript options:
tdrStyle->SetPaperSize(20.,20.);
// tdrStyle->SetLineScalePS(Float_t scale = 3);
// tdrStyle->SetLineStyleString(Int_t i, const char* text);
// tdrStyle->SetHeaderPS(const char* header);
// tdrStyle->SetTitlePS(const char* pstitle);
// tdrStyle->SetBarOffset(Float_t baroff = 0.5);
// tdrStyle->SetBarWidth(Float_t barwidth = 0.5);
// tdrStyle->SetPaintTextFormat(const char* format = "g");
// tdrStyle->SetPalette(Int_t ncolors = 0, Int_t* colors = 0);
// tdrStyle->SetTimeOffset(Double_t toffset);
// tdrStyle->SetHistMinimumZero(kTRUE);
tdrStyle->cd();
}
void Plotter::Plot(std::string filename) {
gROOT->Reset();
//gROOT->SetStyle("Plain");
TStyle *MyStyle = new TStyle("MyStyle","My Root Styles");
MyStyle->SetPaperSize(20,24);
MyStyle->SetStatColor(0);
MyStyle->SetCanvasColor(0);
MyStyle->SetPadColor(0);
MyStyle->SetPadBorderMode(0);
MyStyle->SetCanvasBorderMode(0);
MyStyle->SetFrameBorderMode(0);
MyStyle->SetOptStat(0);
MyStyle->SetStatBorderSize(2);
MyStyle->SetOptTitle(0);
MyStyle->SetPadTickX(1);
MyStyle->SetPadTickY(1);
MyStyle->SetPadBorderSize(2);
MyStyle->SetPalette(51, 0);
MyStyle->SetPadBottomMargin(0.15);
MyStyle->SetPadTopMargin(0.05);
MyStyle->SetPadLeftMargin(0.15);
MyStyle->SetPadRightMargin(0.25);
MyStyle->SetTitleColor(1);
MyStyle->SetTitleFillColor(0);
MyStyle->SetTitleFontSize(0.05);
MyStyle->SetTitleBorderSize(0);
MyStyle->SetLineWidth(1);
MyStyle->SetHistLineWidth(3);
MyStyle->SetLegendBorderSize(0);
MyStyle->SetNdivisions(502, "x");
MyStyle->SetMarkerSize(0.8);
MyStyle->SetTickLength(0.03);
MyStyle->SetTitleOffset(1.5, "x");
MyStyle->SetTitleOffset(1.5, "y");
MyStyle->SetTitleOffset(1.0, "z");
MyStyle->SetLabelSize(0.05, "x");
MyStyle->SetLabelSize(0.05, "y");
MyStyle->SetLabelSize(0.05, "z");
MyStyle->SetLabelOffset(0.03, "x");
MyStyle->SetLabelOffset(0.03, "y");
MyStyle->SetLabelOffset(0.03, "z");
MyStyle->SetTitleSize(0.05, "x");
MyStyle->SetTitleSize(0.05, "y");
MyStyle->SetTitleSize(0.05, "z");
gROOT->SetStyle("MyStyle");
//gROOT->ForceStyle();
bool DrawLog = false;
for (int i = 0; i < N_histos; ++i) {
THStack *hs;
TLegend *l;
int Nset = data.size() + bg.size() + signal.size();
if (Nset > 20)
Nset = 20.;
l = new TLegend(0.76, 0.95 - 0.8 * Nset / 20., 1.0, 0.95);
l->SetFillStyle(1001);
l->SetFillColor(kWhite);
l->SetLineColor(kWhite);
l->SetLineWidth(2);
if (bg.size() > 0) {
hs = new THStack("hs", bg.at(0).at(i)->GetName());
int j = 0;
for (std::vector<std::vector<TH1F*> >::const_iterator it = bg.begin(); it != bg.end(); ++it) {
TH1F *histogram = it->at(i);
switch (j) {
case 0:
histogram->SetFillColor(kRed);
break;
case 1:
histogram->SetFillColor(kOrange);
break;
case 2:
histogram->SetFillColor(kYellow);
break;
case 3:
histogram->SetFillColor(kGreen);
break;
case 4:
histogram->SetFillColor(kCyan);
break;
case 5:
histogram->SetFillColor(kBlue);
break;
case 6:
histogram->SetFillColor(kMagenta);
break;
case 7:
histogram->SetFillColor(kGray);
break;
case 8:
histogram->SetFillColor(kGray + 2);
break;
default:
histogram->SetFillColor(kBlack);
break;
}
hs->Add(histogram);
l->AddEntry(histogram, bg_names.at(j).c_str(), "f");
++j;
}
}
TCanvas *c = new TCanvas("c", "c", 800, 600);
c->SetLogy(DrawLog);
std::string plotname;
if (data.size() > 0) {
plotname = std::string(data.at(0).at(i)->GetName());
data.at(0).at(i)->SetMaximum(1.5 * data.at(0).at(i)->GetMaximum());
data.at(0).at(i)->GetXaxis()->SetTitleOffset(1.3);
data.at(0).at(i)->GetYaxis()->SetTitleOffset(1.3);
data.at(0).at(i)->GetYaxis()->SetTitle("Events");
data.at(0).at(i)->GetXaxis()->SetNdivisions(505);
data.at(0).at(i)->Draw("");
l->AddEntry(data.at(0).at(i), data_names.at(0).c_str(), "p");
if (bg.size() > 0)
hs->Draw("histsame");
data.at(0).at(i)->SetMarkerStyle(20);
data.at(0).at(i)->Draw("psame");
l->Draw("same");
}
if (data.size() == 0 && bg.size() > 0) {
plotname = std::string(bg.at(0).at(i)->GetName());
hs->Draw("hist");
hs->GetXaxis()->SetTitleOffset(1.3);
hs->GetXaxis()->SetNdivisions(505);
hs->GetYaxis()->SetTitleOffset(1.3);
if (bg.size() > 0)
hs->GetXaxis()->SetTitle(bg.at(0).at(i)->GetXaxis()->GetTitle());
hs->GetYaxis()->SetTitle("Events");
l->Draw("same");
}
// if (i==2) {
// TLine *line = new TLine(2, 500*1000 , 2, 0);
// line->SetLineColor(kBlack);
// line->Draw("histsame");
// }
// if (i==4) {
// TLine *line = new TLine(1, 35000 , 1, 0);
// line->SetLineColor(kBlack);
// line->Draw("histsame");
// }
// if (i==6) {
// TLine *line = new TLine(1, 900 , 1, 0);
// line->SetLineColor(kBlack);
// line->Draw("histsame");
// }
// c->Print((filename+std::string("_")+plotname+std::string(".pdf")).c_str());
if (i == 0 && N_histos > 1) {
c->Print((filename+std::string("(")).c_str());
} else if (i > 0 && i == N_histos - 1)
c->Print((filename+std::string(")")).c_str());
else
c->Print(filename.c_str());
}
}
void absoluteCS_CSnPb()
{
TStyle * style = (TStyle*)gROOT->FindObject("graphStyle");
if(!style)
{
style = new TStyle("graphStyle","graphStyle");
}
TCanvas* c = new TCanvas("c1","",1200,1200);
style->SetOptStat(0);
style->SetOptTitle(0);
style->SetPalette(1,0);
style->SetCanvasColor(10);
style->SetCanvasBorderMode(0);
style->SetFrameLineWidth(3);
style->SetFrameFillColor(10);
style->SetPadColor(10);
style->SetHistLineWidth(3);
style->SetHistLineColor(kBlue);
style->SetMarkerSize(0.9);
style->SetMarkerStyle(8);
style->SetFuncWidth(3);
style->SetFuncColor(kRed);
style->SetLabelColor(kBlack,"xyz");
style->SetTitleSize(0.06,"xyz");
style->SetTitleFillColor(10);
style->SetTitleTextColor(kBlack);
style->SetEndErrorSize(0);
gROOT->SetStyle("graphStyle");
gROOT->ForceStyle();
// read graphs
string expFileName = "/data2/analysis/total.root";
string litFileName = "/data2/analysis/literatureData.root";
TFile* expFile = new TFile(expFileName.c_str(),"READ");
TFile* litFile = new TFile(litFileName.c_str(),"READ");
string expCGraphName = "CNat";
string expSnGraphName = "SnNat";
string expPbGraphName = "PbNat";
string litCGraphName = "Natural C (n,tot)";
string litSnGraphName = "Natural Sn (n,tot)";
string litPbGraphName = "Natural Pb (n,tot)";
TGraphAsymmErrors* expCGraph = (TGraphAsymmErrors*)expFile->Get(expCGraphName.c_str());
TGraphAsymmErrors* expSnGraph = (TGraphAsymmErrors*)expFile->Get(expSnGraphName.c_str());
TGraphAsymmErrors* expPbGraph = (TGraphAsymmErrors*)expFile->Get(expPbGraphName.c_str());
TGraphAsymmErrors* litCGraph = (TGraphAsymmErrors*)litFile->Get(litCGraphName.c_str());
TGraphAsymmErrors* litSnGraph = (TGraphAsymmErrors*)litFile->Get(litSnGraphName.c_str());
TGraphAsymmErrors* litPbGraph = (TGraphAsymmErrors*)litFile->Get(litPbGraphName.c_str());
// Set graph point and line characteristics
expCGraph->SetLineColor(kRed);
expCGraph->SetLineWidth(5);
expCGraph->SetLineStyle(0);
expCGraph->SetMarkerColor(kRed);
expSnGraph->SetLineColor(kRed);
expSnGraph->SetLineWidth(5);
expSnGraph->SetLineStyle(0);
expSnGraph->SetMarkerColor(kRed);
expPbGraph->SetLineColor(kRed);
expPbGraph->SetLineWidth(5);
expPbGraph->SetLineStyle(0);
expPbGraph->SetMarkerColor(kRed);
litCGraph->SetLineColor(kBlack);
litCGraph->SetLineWidth(3);
litCGraph->SetLineStyle(2);
litCGraph->SetMarkerColor(kBlack);
litSnGraph->SetLineColor(kBlack);
litSnGraph->SetLineWidth(3);
litSnGraph->SetLineStyle(2);
litSnGraph->SetMarkerColor(kBlack);
litPbGraph->SetLineColor(kBlack);
litPbGraph->SetLineWidth(3);
litPbGraph->SetLineStyle(2);
litPbGraph->SetMarkerColor(kBlack);
// Pad dimensions and margins
gPad->SetPad(0.005, 0.995, 0.995, 0.005);
gPad->SetLeftMargin(0.15);
gPad->SetRightMargin(0.01);
gPad->SetTopMargin(0.03);
gPad->SetBottomMargin(0.15);
gPad->SetTicky(2);
// X-axis parameters
expCGraph->GetXaxis()->SetTitle("Energy (MeV)");
expCGraph->GetXaxis()->SetTitleSize(0.05);
expCGraph->GetXaxis()->SetTitleFont(2);
expCGraph->GetXaxis()->SetTitleOffset(1.4);
expCGraph->GetXaxis()->CenterTitle();
expCGraph->GetXaxis()->SetLabelOffset(0.01);
expCGraph->GetXaxis()->SetLabelSize(0.05);
expCGraph->GetXaxis()->SetLabelFont(2);
expCGraph->GetXaxis()->SetNdivisions(10);
expCGraph->GetXaxis()->SetTickLength(0.03);
// Y-axis parameters
expCGraph->GetYaxis()->SetTitle("#sigma_{tot} (barns)");
expCGraph->GetYaxis()->SetTitleSize(0.06);
expCGraph->GetYaxis()->SetTitleFont(2);
expCGraph->GetYaxis()->SetTitleOffset(0.8);
expCGraph->GetYaxis()->CenterTitle();
expCGraph->GetYaxis()->SetLabelOffset(0.01);
expCGraph->GetYaxis()->SetLabelSize(0.05);
expCGraph->GetYaxis()->SetLabelFont(2);
expCGraph->GetYaxis()->SetNdivisions(10);
expCGraph->GetYaxis()->SetTickLength(0.02);
expCGraph->Draw("");
expSnGraph->Draw("same");
expPbGraph->Draw("same");
litCGraph->Draw("same");
litSnGraph->Draw("same");
litPbGraph->Draw("same");
expCGraph->Draw("same");
expSnGraph->Draw("same");
expPbGraph->Draw("same");
gPad->SetLogx(1);
expCGraph->GetYaxis()->SetRangeUser(0,9);
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.035);
latex.SetTextAlign(13); // align at top
latex.DrawLatex(0.65,0.65,"Pb (elem.)");
latex.DrawLatex(0.35,0.52,"Sn (elem.)");
latex.DrawLatex(0.32,0.4,"C (elem.)");
// Define legend format and contents
TLegend *legend = new TLegend(0.7,0.8,0.9,0.9);
legend->AddEntry(litCGraph,"lit data (analog)","l");
legend->AddEntry(expCGraph,"new data (DSP)","l");
legend->Draw();
expFile->Close();
litFile->Close();
}
void setTDRStyle() {
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");
// For the canvas:
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetCanvasColor(kWhite);
tdrStyle->SetCanvasDefH(600); //Height of canvas
tdrStyle->SetCanvasDefW(600); //Width of canvas
tdrStyle->SetCanvasDefX(0); //POsition on screen
tdrStyle->SetCanvasDefY(0);
// For the Pad:
tdrStyle->SetPadBorderMode(0);
tdrStyle->SetPadColor(kWhite);
tdrStyle->SetPadGridX(false);
tdrStyle->SetPadGridY(false);
tdrStyle->SetGridColor(0);
tdrStyle->SetGridStyle(3);
tdrStyle->SetGridWidth(1);
// For the frame:
tdrStyle->SetFrameBorderMode(0);
tdrStyle->SetFrameBorderSize(1);
tdrStyle->SetFrameFillColor(0);
tdrStyle->SetFrameFillStyle(0);
tdrStyle->SetFrameLineColor(1);
tdrStyle->SetFrameLineStyle(1);
tdrStyle->SetFrameLineWidth(1);
// For the histo:
tdrStyle->SetHistLineColor(1);
tdrStyle->SetHistLineStyle(0);
tdrStyle->SetHistLineWidth(1);
tdrStyle->SetEndErrorSize(2);
tdrStyle->SetErrorX(0.);
tdrStyle->SetMarkerStyle(20);
//For the fit/function:
tdrStyle->SetOptFit(1);
tdrStyle->SetFitFormat("5.4g");
tdrStyle->SetFuncColor(2);
tdrStyle->SetFuncStyle(1);
tdrStyle->SetFuncWidth(1);
//For the date:
tdrStyle->SetOptDate(0);
// For the statistics box:
tdrStyle->SetOptFile(0);
tdrStyle->SetOptStat(0); // To display the mean and RMS: SetOptStat("mr");
tdrStyle->SetStatColor(kWhite);
tdrStyle->SetStatFont(42);
tdrStyle->SetStatFontSize(0.025);
tdrStyle->SetStatTextColor(1);
tdrStyle->SetStatFormat("6.4g");
tdrStyle->SetStatBorderSize(1);
tdrStyle->SetStatH(0.1);
tdrStyle->SetStatW(0.15);
// Margins:
tdrStyle->SetPadTopMargin(0.05);
tdrStyle->SetPadBottomMargin(0.13);
tdrStyle->SetPadLeftMargin(0.16);
tdrStyle->SetPadRightMargin(0.02);
// For the Global title:
tdrStyle->SetOptTitle(0);
tdrStyle->SetTitleFont(42);
tdrStyle->SetTitleColor(1);
tdrStyle->SetTitleTextColor(1);
tdrStyle->SetTitleFillColor(10);
tdrStyle->SetTitleFontSize(0.05);
// For the axis titles:
tdrStyle->SetTitleColor(1, "XYZ");
tdrStyle->SetTitleFont(42, "XYZ");
tdrStyle->SetTitleSize(0.06, "XYZ");
tdrStyle->SetTitleXOffset(0.9);
tdrStyle->SetTitleYOffset(1.25);
// For the axis labels:
tdrStyle->SetLabelColor(1, "XYZ");
tdrStyle->SetLabelFont(42, "XYZ");
tdrStyle->SetLabelOffset(0.007, "XYZ");
tdrStyle->SetLabelSize(0.05, "XYZ");
// For the axis:
tdrStyle->SetAxisColor(1, "XYZ");
tdrStyle->SetStripDecimals(kTRUE);
tdrStyle->SetTickLength(0.03, "XYZ");
tdrStyle->SetNdivisions(510, "XYZ");
tdrStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame
tdrStyle->SetPadTickY(1);
// Change for log plots:
tdrStyle->SetOptLogx(0);
tdrStyle->SetOptLogy(0);
tdrStyle->SetOptLogz(0);
// Postscript options:
tdrStyle->SetPaperSize(20.,20.);
tdrStyle->cd();
}
