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();
}
// 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){
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();
}
/////////
// 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 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 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 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 scan_comp_engine( TString engine1_file, TString engine_file) {
engine *myengine1 = new engine(engine1_file,0);
// TString engine_file="paw/hms_dc_52949.root";
engine *myengine = new engine(engine_file,0);
Long64_t nent_engine1 = myengine1->fChain->GetEntriesFast();
Long64_t nent_engine = myengine->fChain->GetEntriesFast();
cout << " engine1 entries = " << nent_engine1 << " engine entries = " << nent_engine << endl;
Long64_t nb_engine1 = 0,nb_engine = 0;
//
TH1F *hntr_engine1 = new TH1F("hntr_engine1","; Ntracks; Counts",50,0,50);
TH1F *hntr_eng = new TH1F("hntr_eng","; Ntracks; Counts",50,0,50);
TH2F *hntr_2d = new TH2F("hntr_2d","; ENGINE1 Ntracks ; ENGINE Ntracks",40,0,40,40,0,40);
TH1F *hntr_diff = new TH1F("hntr_diff","; Ntracks; Counts",50,0,50);
TH1F *hxpfp_eng[3],*hxpfp_engine1[3],*hxpfp_diff[3],*hxpfp_sub[3];
TH1F *hypfp_eng[3],*hypfp_engine1[3],*hypfp_diff[3],*hypfp_sub[3];
TH1F *hyfp_eng[3],*hyfp_engine1[3],*hyfp_diff[3],*hyfp_sub[3];
TH1F *hxfp_eng[3],*hxfp_engine1[3],*hxfp_diff[3],*hxfp_sub[3];
//
Int_t i;
for (i=0 ; i<3 ;i++) {
hxpfp_eng[i] = new TH1F(Form("hxpfp_eng%d",i),"; Xpfp (mr); Counts",100,-100.,100.);
hxpfp_engine1[i] = new TH1F(Form("hxpfp_engine1%d",i),"; Xpfp (mr); Counts",100,-100.,100.);
hxpfp_diff[i] = new TH1F(Form("hxpfp_diff%d",i),"Difference; Xpfp (mr); Counts",100,-100.,100.);
hypfp_eng[i] = new TH1F(Form("hypfp_eng%d",i),"; Ypfp (mr); Counts",100,-50.,50.);
hypfp_engine1[i] = new TH1F(Form("hypfp_engine1%d",i),"; Ypfp (mr); Counts",100,-50.,50.);
hypfp_diff[i] = new TH1F(Form("hypfp_diff%d",i),"Difference; Ypfp (mr); Counts",100,-50.,50.);
hxfp_eng[i] = new TH1F(Form("hxfp_eng%d",i),"; Xfp (cm); Counts",100,-50,50);
hxfp_engine1[i] = new TH1F(Form("hxfp_engine1%d",i),"; Xfp (cm); Counts",100,-50,50);
hxfp_diff[i] = new TH1F(Form("hxfp_diff%d",i),"Difference; Xfp (cm); Counts",100,-50,50);
hyfp_eng[i] = new TH1F(Form("hyfp_eng%d",i),"; Yfp (cm); Counts",60,-30,30);
hyfp_engine1[i] = new TH1F(Form("hyfp_engine1%d",i),"; Yfp (cm); Counts",60,-30,30);
hyfp_diff[i] = new TH1F(Form("hyfp_diff%d",i),"Difference; Yfp (cm); Counts",60,-30,30);
hyfp_sub[i] = new TH1F(Form("hyfp_sub%d",i),"; Yfp (ENGINE1-ENGINE); Counts",60,-.5,.5);
hxfp_sub[i] = new TH1F(Form("hxfp_sub%d",i),"; Xfp (ENGINE1-ENGINE); Counts",100,-.5,.5);
hxpfp_sub[i] = new TH1F(Form("hxpfp_sub%d",i),"; Xpfp (ENGINE1-ENGINE); Counts",100,-.2,.2);
hypfp_sub[i] = new TH1F(Form("hypfp_sub%d",i),"; Ypfp (ENGINE1-ENGINE); Counts",100,-.2,.2);
}//
//
Int_t counter=0;
for (Long64_t ni=0 ; ni<nent_engine1 ;ni++) {
nb_engine1 = myengine1->fChain->GetEntry(ni);
nb_engine = myengine->fChain->GetEntry(ni);
hntr_engine1->Fill(myengine1->dc_ntr);
hntr_eng->Fill(myengine->dc_ntr);
if ((myengine1->dc_ntr) != (myengine->dc_ntr) )hntr_2d->Fill(myengine1->dc_ntr,myengine->dc_ntr);
if ( myengine1->dc_ntr==myengine->dc_ntr &&myengine1->dc_ntr<=3 ) {
for (i=0 ; i<myengine->dc_ntr ;i++) {
hxfp_engine1[i]->Fill(myengine1->dc_xfp[i]);
hyfp_engine1[i]->Fill(myengine1->dc_yfp[i]);
hxpfp_engine1[i]->Fill(myengine1->dc_xpfp[i]*1000);
hypfp_engine1[i]->Fill(myengine1->dc_ypfp[i]*1000);
hxfp_eng[i]->Fill(myengine->dc_xfp[i]);
hyfp_eng[i]->Fill(myengine->dc_yfp[i]);
hxpfp_eng[i]->Fill(myengine->dc_xpfp[i]*1000);
hypfp_eng[i]->Fill(myengine->dc_ypfp[i]*1000);
if ( TMath::Abs((myengine1->dc_yfp[i]-myengine->dc_yfp[i])) != 0.0 || (myengine1->dc_xfp[i]-myengine->dc_xfp[i]) != 0.0
|| (myengine1->dc_xpfp[i]-myengine->dc_xpfp[i]) != 0.0 || (myengine1->dc_ypfp[i]-myengine->dc_ypfp[i]) != 0.0 ) {
hyfp_sub[i]->Fill((myengine1->dc_yfp[i]-myengine->dc_yfp[i]));
hxfp_sub[i]->Fill((myengine1->dc_xfp[i]-myengine->dc_xfp[i]));
hypfp_sub[i]->Fill((myengine1->dc_ypfp[i]-myengine->dc_ypfp[i]));
hxpfp_sub[i]->Fill((myengine1->dc_xpfp[i]-myengine->dc_xpfp[i]));
}
}
}
if ( myengine->dc_ntr == 1 && myengine1->dc_ntr == 1 && -1==1) {
cout << " Counter = "<< counter++ << endl;
myengine1->PrintTrack(ni);
myengine->PrintTrack(ni);
}
}
cout << " Engine1 File= " << engine1_file << endl;
cout << " Engine File= " << engine_file << endl;
TStyle *MyStyle = new TStyle("MyStyle"," My Style");
MyStyle->SetOptStat(1000000);
MyStyle->SetTitleOffset(1.,"Y");
MyStyle->SetTitleOffset(.7,"X");
MyStyle->SetLabelSize(0.04,"XY");
MyStyle->SetTitleSize(0.06,"XY");
MyStyle->SetPadLeftMargin(0.12);
MyStyle->SetStatFontSize(0.2);
MyStyle->SetTitleFontSize(0.1);
gROOT->SetStyle("MyStyle");
//
TCanvas *cplot2[3];
for (i=0 ; i<3 ;i++) {
cplot2[i] = new TCanvas(Form("cplot2_%d",i),Form(" Focal Plane Positions ( %d track) ",i+1),800,800);
cplot2[i]->Divide(2,3);
cplot2[i]->cd(1);
hxfp_engine1[i]->Draw();
hxfp_eng[i]->Draw("same");
hxfp_eng[i]->SetLineColor(2);
cplot2[i]->cd(3);
hxfp_diff[i]->Add(hxfp_engine1[i],hxfp_eng[i],-1);
hxfp_diff[i]->Draw();
cplot2[i]->cd(5);
hxfp_sub[i]->Draw();
cplot2[i]->cd(2);
hyfp_engine1[i]->Draw();
hyfp_eng[i]->Draw("same");
hyfp_eng[i]->SetLineColor(2);
cplot2[i]->cd(4);
hyfp_diff[i]->Add(hyfp_engine1[i],hyfp_eng[i],-1);
hyfp_diff[i]->Draw();
cplot2[i]->cd(6);
hyfp_sub[i]->Draw();
}
TCanvas *cplot3[3];
for (i=0 ; i<3 ;i++) {
cplot3[i]= new TCanvas(Form("cplot3_%d",i),Form(" Focal Plane Angles (%d track) ",i+1),800,800);
cplot3[i]->Divide(2,3);
cplot3[i]->cd(1);
hxpfp_engine1[i]->Draw();
hxpfp_eng[i]->Draw("same");
hxpfp_eng[i]->SetLineColor(2);
cplot3[i]->cd(3);
hxpfp_diff[i]->Add(hxpfp_engine1[i],hxpfp_eng[i],-1);
hxpfp_diff[i]->Draw();
cplot3[i]->cd(5);
hxpfp_sub[i]->Draw();
cplot3[i]->cd(2);
hypfp_engine1[i]->Draw();
hypfp_eng[i]->Draw("same");
hypfp_eng[i]->SetLineColor(2);
cplot3[i]->cd(4);
hypfp_diff[i]->Add(hypfp_engine1[i],hypfp_eng[i],-1);
hypfp_diff[i]->Draw();
cplot3[i]->cd(6);
hypfp_sub[i]->Draw();
}
TCanvas *cplot1 = new TCanvas("cplot1","2D ntrack ",800,800);
hntr_2d->Draw("colz");
TCanvas *cplot = new TCanvas("cplot","1D ntrack ",800,800);
cplot->Divide(1,3);
cplot->cd(1);
hntr_engine1->Draw();
cplot->cd(2);
hntr_eng->Draw();
hntr_eng->SetLineColor(2);
cplot->cd(3);
hntr_diff->Add(hntr_engine1,hntr_eng,-1);
hntr_diff->Draw();
//
}
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();
}
TStyle* sPhenixStyle()
{
TStyle *sphenixStyle = new TStyle("sPHENIX","sPHENIX style");
// use plain black on white colors
Int_t icol=0; // WHITE
sphenixStyle->SetFrameBorderMode(icol);
sphenixStyle->SetFrameFillColor(icol);
sphenixStyle->SetCanvasBorderMode(icol);
sphenixStyle->SetCanvasColor(icol);
sphenixStyle->SetPadBorderMode(icol);
sphenixStyle->SetPadColor(icol);
sphenixStyle->SetStatColor(icol);
//sphenixStyle->SetFillColor(icol); // don't use: white fill color for *all* objects
// set the paper & margin sizes
sphenixStyle->SetPaperSize(20,26);
// set margin sizes
sphenixStyle->SetPadTopMargin(0.05);
sphenixStyle->SetPadRightMargin(0.05);
sphenixStyle->SetPadBottomMargin(0.16);
sphenixStyle->SetPadLeftMargin(0.16);
// set title offsets (for axis label)
sphenixStyle->SetTitleXOffset(1.4);
sphenixStyle->SetTitleYOffset(1.4);
// use large fonts
//Int_t font=72; // Helvetica italics
Int_t font=42; // Helvetica
Double_t tsize=0.05;
sphenixStyle->SetTextFont(font);
sphenixStyle->SetTextSize(tsize);
sphenixStyle->SetLabelFont(font,"x");
sphenixStyle->SetTitleFont(font,"x");
sphenixStyle->SetLabelFont(font,"y");
sphenixStyle->SetTitleFont(font,"y");
sphenixStyle->SetLabelFont(font,"z");
sphenixStyle->SetTitleFont(font,"z");
sphenixStyle->SetLabelSize(tsize,"x");
sphenixStyle->SetTitleSize(tsize,"x");
sphenixStyle->SetLabelSize(tsize,"y");
sphenixStyle->SetTitleSize(tsize,"y");
sphenixStyle->SetLabelSize(tsize,"z");
sphenixStyle->SetTitleSize(tsize,"z");
// use bold lines and markers
sphenixStyle->SetMarkerStyle(20);
sphenixStyle->SetMarkerSize(1.2);
sphenixStyle->SetHistLineWidth(2.);
sphenixStyle->SetLineStyleString(2,"[12 12]"); // postscript dashes
// get rid of X error bars
//sphenixStyle->SetErrorX(0.001);
// get rid of error bar caps
sphenixStyle->SetEndErrorSize(0.);
// do not display any of the standard histogram decorations
sphenixStyle->SetOptTitle(0);
//sphenixStyle->SetOptStat(1111);
sphenixStyle->SetOptStat(0);
//sphenixStyle->SetOptFit(1111);
sphenixStyle->SetOptFit(0);
// put tick marks on top and RHS of plots
sphenixStyle->SetPadTickX(1);
sphenixStyle->SetPadTickY(1);
// legend modificatin
sphenixStyle->SetLegendBorderSize(0);
sphenixStyle->SetLegendFillColor(0);
sphenixStyle->SetLegendFont(font);
#if ROOT_VERSION_CODE >= ROOT_VERSION(6,00,0)
std::cout << "sPhenixStyle: ROOT6 mode" << std::endl;
sphenixStyle->SetLegendTextSize(tsize);
sphenixStyle->SetPalette(kBird);
#else
std::cout << "sPhenixStyle: ROOT5 mode" << std::endl;
// color palette - manually define 'kBird' palette only available in ROOT 6
Int_t alpha = 0;
Double_t stops[9] = { 0.0000, 0.1250, 0.2500, 0.3750, 0.5000, 0.6250, 0.7500, 0.8750, 1.0000};
Double_t red[9] = { 0.2082, 0.0592, 0.0780, 0.0232, 0.1802, 0.5301, 0.8186, 0.9956, 0.9764};
Double_t green[9] = { 0.1664, 0.3599, 0.5041, 0.6419, 0.7178, 0.7492, 0.7328, 0.7862, 0.9832};
Double_t blue[9] = { 0.5293, 0.8684, 0.8385, 0.7914, 0.6425, 0.4662, 0.3499, 0.1968, 0.0539};
TColor::CreateGradientColorTable(9, stops, red, green, blue, 255, alpha);
#endif
sphenixStyle->SetNumberContours(80);
return sphenixStyle;
}
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 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();
}
// Create my custom FORMAL style
// *****************************
void SemiFormalStyle(void) {
std::cout << " *** void getSemiFormalStyle(void) *** Custom TStyle initialized!" << std::endl;
TStyle *semiFStyle = new TStyle("semiFStyle","Style for AN09");
// For the canvas:
semiFStyle->SetCanvasBorderMode(0);
semiFStyle->SetCanvasColor(kWhite);
semiFStyle->SetCanvasDefH(600); //Height of canvas
semiFStyle->SetCanvasDefW(600); //Width of canvas
semiFStyle->SetCanvasDefX(0); //POsition on screen
semiFStyle->SetCanvasDefY(0);
// For the Pad:
semiFStyle->SetPadBorderMode(0);
// semiFStyle->SetOptTitle(0); // delete histogram titles
// semiFStyle->SetPadBorderSize(Width_t size = 1);
semiFStyle->SetPadColor(kWhite);
semiFStyle->SetPadGridX(true);
semiFStyle->SetPadGridY(true);
semiFStyle->SetGridColor(0);
semiFStyle->SetGridStyle(3);
semiFStyle->SetGridWidth(1);
// For the frame:
semiFStyle->SetFrameBorderMode(0);
semiFStyle->SetFrameBorderSize(1);
semiFStyle->SetFrameFillColor(0);
semiFStyle->SetFrameFillStyle(0);
semiFStyle->SetFrameLineColor(1);
semiFStyle->SetFrameLineStyle(1);
semiFStyle->SetFrameLineWidth(1);
// For the histo:
// semiFStyle->SetHistFillColor(1);
// semiFStyle->SetHistFillStyle(0);
semiFStyle->SetHistLineColor(1);
semiFStyle->SetHistLineStyle(0);
semiFStyle->SetHistLineWidth(1);
// semiFStyle->SetLegoInnerR(Float_t rad = 0.5);
// semiFStyle->SetNumberContours(Int_t number = 20);
semiFStyle->SetEndErrorSize(2);
//semiFStyle->SetErrorMarker(20);
semiFStyle->SetErrorX(0.);
semiFStyle->SetMarkerStyle(20);
//For the fit/function:
semiFStyle->SetOptFit(1);
semiFStyle->SetFitFormat("5.4g");
semiFStyle->SetFuncColor(2);
semiFStyle->SetFuncStyle(1);
semiFStyle->SetFuncWidth(1);
//For the date:
semiFStyle->SetOptDate(1);
// For the statistics box:
semiFStyle->SetOptFile(0);
semiFStyle->SetOptStat("mnre"); // stats box
semiFStyle->SetOptFit(1111); // fit box
semiFStyle->SetStatColor(kWhite);
semiFStyle->SetStatFont(42);
semiFStyle->SetStatFontSize(0.025);
semiFStyle->SetStatTextColor(1);
semiFStyle->SetStatFormat("6.4g");
semiFStyle->SetStatBorderSize(1);
semiFStyle->SetStatH(0.1);
semiFStyle->SetStatW(0.15);
// Margins:
semiFStyle->SetPadTopMargin(0.15);
semiFStyle->SetPadBottomMargin(0.13);
semiFStyle->SetPadLeftMargin(0.13);
semiFStyle->SetPadRightMargin(0.15);
// For the Global title:
// semiFStyle->SetOptTitle(0);
semiFStyle->SetTitleFont(42);
semiFStyle->SetTitleColor(1);
semiFStyle->SetTitleTextColor(1);
semiFStyle->SetTitleFillColor(10);
semiFStyle->SetTitleFontSize(0.05);
// For the axis titles:
semiFStyle->SetTitleColor(1, "XYZ");
semiFStyle->SetTitleFont(42, "XYZ");
semiFStyle->SetTitleSize(0.06, "XYZ");
semiFStyle->SetTitleXOffset(1.5);
semiFStyle->SetTitleYOffset(1.5);
semiFStyle->SetTitleOffset(1.5, "Y"); // Another way to set the Offset
// For the axis labels:
semiFStyle->SetLabelColor(1, "XYZ");
semiFStyle->SetLabelFont(42, "XYZ");
semiFStyle->SetLabelOffset(0.007, "XYZ");
semiFStyle->SetLabelSize(0.05, "XYZ");
// For the axis:
semiFStyle->SetAxisColor(1, "XYZ");
semiFStyle->SetStripDecimals(kTRUE);
semiFStyle->SetTickLength(0.03, "XYZ");
semiFStyle->SetNdivisions(510, "XYZ");
semiFStyle->SetPadTickX(0); // To get tick marks on the opposite side of the frame
semiFStyle->SetPadTickY(0);
// Change for log plots:
semiFStyle->SetOptLogx(0);
semiFStyle->SetOptLogy(0);
semiFStyle->SetOptLogz(0);
semiFStyle->SetPalette(1,0);
semiFStyle->cd();
} //end of: void SemiFormalStyle(void){
// Create my custom FORMAL style
//********************************************************************//
void getMyStyle(void) {
std::cout << " *** void getMyStyle(void) *** Custom TStyle initialized!" << std::endl;
TStyle *myStyle = new TStyle("myStyle","Style for AN09");
// For the canvas:
myStyle->SetCanvasBorderMode(0);
myStyle->SetCanvasColor(kWhite);
myStyle->SetCanvasDefH(600); //Height of canvas
myStyle->SetCanvasDefW(600); //Width of canvas
myStyle->SetCanvasDefX(0); //POsition on screen
myStyle->SetCanvasDefY(0);
// For the Pad:
myStyle->SetPadBorderMode(0);
// myStyle->SetOptTitle(0); // delete histogram titles
// myStyle->SetPadBorderSize(Width_t size = 1);
myStyle->SetPadColor(kWhite);
myStyle->SetPadGridX(true);
myStyle->SetPadGridY(true);
myStyle->SetGridColor(0);
myStyle->SetGridStyle(3);
myStyle->SetGridWidth(1);
// For the frame:
myStyle->SetFrameBorderMode(0);
myStyle->SetFrameBorderSize(1);
myStyle->SetFrameFillColor(0);
myStyle->SetFrameFillStyle(0);
myStyle->SetFrameLineColor(1);
myStyle->SetFrameLineStyle(1);
myStyle->SetFrameLineWidth(1);
// For the histo:
// myStyle->SetHistFillColor(1);
// myStyle->SetHistFillStyle(0);
myStyle->SetHistLineColor(1);
myStyle->SetHistLineStyle(0);
myStyle->SetHistLineWidth(1);
// myStyle->SetLegoInnerR(Float_t rad = 0.5);
// myStyle->SetNumberContours(Int_t number = 20);
myStyle->SetEndErrorSize(2);
//myStyle->SetErrorMarker(20);
myStyle->SetErrorX(0.);
myStyle->SetMarkerStyle(20);
//For the fit/function:
myStyle->SetOptFit(1);
myStyle->SetFitFormat("5.4g");
myStyle->SetFuncColor(2);
myStyle->SetFuncStyle(1);
myStyle->SetFuncWidth(1);
//For the date:
myStyle->SetOptDate(1);
// myStyle->SetDateX(Float_t x = 0.01);
// myStyle->SetDateY(Float_t y = 0.01);
// For the statistics box:
myStyle->SetOptFile(0);
myStyle->SetOptStat("miourne"); // stats box
myStyle->SetOptFit(1111); // fit box
myStyle->SetStatColor(kWhite);
myStyle->SetStatFont(42);
myStyle->SetStatFontSize(0.025);
myStyle->SetStatTextColor(1);
myStyle->SetStatFormat("6.4g");
myStyle->SetStatBorderSize(1);
myStyle->SetStatH(0.1);
myStyle->SetStatW(0.15);
// myStyle->SetStatStyle(Style_t style = 1001);
// myStyle->SetStatX(Float_t x = 0);
// myStyle->SetStatY(Float_t y = 0);
// Margins:
myStyle->SetPadTopMargin(0.15);
myStyle->SetPadBottomMargin(0.13);
myStyle->SetPadLeftMargin(0.13);
myStyle->SetPadRightMargin(0.15);
// For the Global title:
// myStyle->SetOptTitle(0);
myStyle->SetTitleFont(42);
myStyle->SetTitleColor(1);
myStyle->SetTitleTextColor(1);
myStyle->SetTitleFillColor(10);
myStyle->SetTitleFontSize(0.05);
// myStyle->SetTitleH(0); // Set the height of the title box
// myStyle->SetTitleW(0); // Set the width of the title box
// myStyle->SetTitleX(0); // Set the position of the title box
// myStyle->SetTitleY(0.985); // Set the position of the title box
// myStyle->SetTitleStyle(Style_t style = 1001);
// myStyle->SetTitleBorderSize(2);
// For the axis titles:
myStyle->SetTitleColor(1, "XYZ");
myStyle->SetTitleFont(42, "XYZ");
myStyle->SetTitleSize(0.06, "XYZ");
// myStyle->SetTitleXSize(Float_t size = 0.02); // Another way to set the size?
// myStyle->SetTitleYSize(Float_t size = 0.02);
myStyle->SetTitleXOffset(1.2);
myStyle->SetTitleYOffset(1.2);
myStyle->SetTitleOffset(1.2, "Y"); // Another way to set the Offset
// For the axis labels:
myStyle->SetLabelColor(1, "XYZ");
myStyle->SetLabelFont(42, "XYZ");
myStyle->SetLabelOffset(0.007, "XYZ");
myStyle->SetLabelSize(0.05, "XYZ");
// For the axis:
myStyle->SetAxisColor(1, "XYZ");
myStyle->SetStripDecimals(kTRUE);
myStyle->SetTickLength(0.03, "XYZ");
myStyle->SetNdivisions(510, "XYZ");
myStyle->SetPadTickX(0); // To get tick marks on the opposite side of the frame
myStyle->SetPadTickY(0);
// Change for log plots:
myStyle->SetOptLogx(0);
myStyle->SetOptLogy(0);
myStyle->SetOptLogz(0);
myStyle->SetPalette(1,0);
// Postscript options:
// myStyle->SetPaperSize(15.,15.);
// myStyle->SetLineScalePS(Float_t scale = 3);
// myStyle->SetLineStyleString(Int_t i, const char* text);
// myStyle->SetHeaderPS(const char* header);
// myStyle->SetTitlePS(const char* pstitle);
// myStyle->SetBarOffset(Float_t baroff = 0.5);
// myStyle->SetBarWidth(Float_t barwidth = 0.5);
// myStyle->SetPaintTextFormat(const char* format = "g");
// myStyle->SetPalette(Int_t ncolors = 0, Int_t* colors = 0);
// myStyle->SetTimeOffset(Double_t toffset);
// myStyle->SetHistMinimumZero(kTRUE);
myStyle->cd();
} //end of: void getMyStyle(void){
/////
// Set setTDRStyle_modified (from link https://twiki.cern.ch/twiki/pub/CMS/TRK10001/setTDRStyle_modified.C)
/////
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->SetHistFillColor(0);
tdrStyle->SetHistLineColor(1);
tdrStyle->SetHistLineStyle(0);
tdrStyle->SetHistLineWidth(1);
tdrStyle->SetErrorX(0.);
tdrStyle->SetMarkerStyle(20);
//For the fit/function:
tdrStyle->SetOptFit(1);
tdrStyle->SetFitFormat("5.4g");
tdrStyle->SetFuncStyle(1);
tdrStyle->SetFuncWidth(1);
//For the date:
tdrStyle->SetOptDate(0);
// For the statistics box:
tdrStyle->SetOptFile(0);
tdrStyle->SetOptStat("e");
tdrStyle->SetStatColor(kGray);
tdrStyle->SetStatFont(42);
tdrStyle->SetStatTextColor(1);
tdrStyle->SetStatFormat("6.4g");
tdrStyle->SetStatBorderSize(0);
tdrStyle->SetStatX(1.); //Starting position on X axis
tdrStyle->SetStatY(1.); //Starting position on Y axis
tdrStyle->SetStatFontSize(0.025); //Vertical Size
tdrStyle->SetStatW(0.15); //Horizontal size
// tdrStyle->SetStatStyle(Style_t style = 1001);
// Margins:
tdrStyle->SetPadTopMargin(0.05);
tdrStyle->SetPadBottomMargin(0.125);
tdrStyle->SetPadLeftMargin(0.105);
tdrStyle->SetPadRightMargin(0.1);
// 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.05, "XYZ");
tdrStyle->SetTitleXOffset(0.9);
tdrStyle->SetTitleYOffset(0.8);
tdrStyle->SetTitleOffset(0.7, "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.045, "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);
tdrStyle->cd();
}
void comp_beam_engine_hcana( TString hcana_file, TString engine_file) {
analyze_hcana_tree *myhcana = new analyze_hcana_tree(hcana_file,0);
analyze_engine_tree *myengine = new analyze_engine_tree(engine_file,0);
Long64_t nent_hcana = myhcana->fChain->GetEntriesFast();
Long64_t nent_engine = myengine->fChain->GetEntriesFast();
Long64_t nent_loop=TMath::Min(nent_hcana,nent_engine);
cout << " hcana entries = " << nent_hcana << " engine entries = " << nent_engine << endl;
Long64_t nb_hcana = 0,nb_engine = 0;
TH1F *comp_frx_raw_adc[2],*comp_fry_raw_adc[2],*comp_frx_adc[2],*comp_fry_adc[2],*comp_frx[2],*comp_fry[2];
TH1F *pos_diff[2],*raw_adc_diff[2], *adc_diff[2];
TString comp[2] = {"eng","hcana"};
Int_t i;
for (i=0 ; i<2 ;i++) {
comp_frx_raw_adc[i] = new TH1F(Form("comp_frx_raw_adc_%s",comp[i].Data()),"frx_raw_adc; Raw ADC value; Events",1200,3200,4400);
comp_fry_raw_adc[i] = new TH1F(Form("comp_fry_raw_adc_%s",comp[i].Data()),"fry_raw_adc; Raw ADC value; Events",1200,3200,4400);
comp_frx_adc[i] = new TH1F(Form("comp_frx_adc_%s",comp[i].Data()),"frx_adc; ADC value; Events",1000,-500,500);
comp_fry_adc[i] = new TH1F(Form("comp_fry_adc_%s",comp[i].Data()),"fry_adc; ADC value; Events",1000,-500,500);
comp_frx[i] = new TH1F(Form("comp_frx_%s",comp[i].Data()),"frx; Position (cm); Events",100,-0.5,0.5);
comp_fry[i] = new TH1F(Form("comp_fry_%s",comp[i].Data()),"fry; Position (cm); Events",100,-0.5,0.5);
pos_diff[i] = new TH1F(Form("pos_diff_%d",i),"Differences in reconstructed position; Difference (cm); Events",100,-0.5,0.5);
raw_adc_diff[i] = new TH1F(Form("raw_adc_dif_%d",i),"Differences in raw ADC values; Difference ; Events",1200,3200,4400);
adc_diff[i] = new TH1F(Form("adc_diff_%d",i),"Differences in pedestal corrected ADC values; Difference; Events",1000,-500,500);
}
//
char *s = new char[1];
Long64_t engine_ent=0;
for (Long64_t ni=0 ; ni<nent_loop ;ni++) {
nb_hcana = myhcana->fChain->GetEntry(ni);
nb_engine = myengine->fChain->GetEntry(engine_ent++);
while (myhcana->g_evnum!=myengine->evnum) {
nb_engine = myengine->fChain->GetEntry(engine_ent++);
}
if (myhcana->fEvtHdr_fEvtType==1&& myengine->evtype==1 && myhcana->g_evnum==myengine->evnum) {
comp_frx_raw_adc[0]->Fill(myengine->frx_raw_adc);
comp_frx_raw_adc[1]->Fill(myhcana->RB_raster_frx_raw_adc);
comp_frx_adc[0]->Fill(myengine->frx_adc);
comp_frx_adc[1]->Fill(myhcana->RB_raster_frx_adc);
comp_frx[0]->Fill(myengine->frx);
comp_frx[1]->Fill(myhcana->RB_raster_frx);
comp_fry_raw_adc[0]->Fill(myengine->fry_raw_adc);
comp_fry_raw_adc[1]->Fill(myhcana->RB_raster_fry_raw_adc);
comp_fry_adc[0]->Fill(myengine->fry_adc);
comp_fry_adc[1]->Fill(myhcana->RB_raster_fry_adc);
comp_fry[0]->Fill(myengine->fry);
comp_fry[1]->Fill(myhcana->RB_raster_fry);
} // if
} // for loop
cout << " Hcana File= " << hcana_file << endl;
cout << " Engine File= " << engine_file << endl;
gStyle->Reset();
TStyle *MyStyle = new TStyle("MyStyle"," My Style");
// MyStyle->SetOptStat(10000010);
MyStyle->SetOptStat("i");
MyStyle->SetTitleOffset(1.5,"Y");
MyStyle->SetTitleOffset(.7,"X");
MyStyle->SetLabelSize(0.04,"XY");
MyStyle->SetTitleSize(0.06,"XY");
MyStyle->SetTitleH(0.04);
MyStyle->SetPadLeftMargin(0.12);
MyStyle->SetStatFontSize(0.2);
MyStyle->SetTitleFontSize(0.1);
gROOT->SetStyle("MyStyle");
//draw X signals
TCanvas * cx = new TCanvas("cx","Compare Raster X Signals",1000,1200);
cx->Divide(2,3);
cx->cd(1);
comp_frx_raw_adc[0]->Draw();
comp_frx_raw_adc[0]->SetLineColor(2);
comp_frx_raw_adc[1]->Draw("same");
legx = new TLegend(0.15,0.8,0.30,0.88);
legx->AddEntry(comp_frx_raw_adc[0],"ENGINE","l");
legx->AddEntry(comp_frx_raw_adc[1],"HCANA","l");
legx->Draw();
cx->cd(2);
raw_adc_diff[0]->Add(comp_frx_raw_adc[1],comp_frx_raw_adc[0],-1);
raw_adc_diff[0]->Draw();
raw_adc_diff[0]->SetLineColor(kGreen-2);
cx->cd(3);
comp_frx_adc[0]->Draw();
comp_frx_adc[0]->SetLineColor(2);
comp_frx_adc[1]->Draw("same");
cx->cd(4);
adc_diff[0]->Add(comp_frx_adc[1],comp_frx_adc[0],-1);
adc_diff[0]->Draw();
adc_diff[0]->SetLineColor(kGreen-2);
cx->cd(5);
comp_frx[0]->Draw();
comp_frx[0]->SetLineColor(2);
comp_frx[1]->Draw("same");
cx->cd(6);
pos_diff[0]->Add(comp_frx[1],comp_frx[0],-1);
pos_diff[0]->Draw();
pos_diff[0]->SetLineColor(kGreen-2);
// draw Y signals
TCanvas * cy = new TCanvas("cy","Compare Raster Y Signals",1000,1200);
cy->Divide(2,3);
cy->cd(1);
comp_fry_raw_adc[0]->Draw();
comp_fry_raw_adc[0]->SetLineColor(2);
comp_fry_raw_adc[1]->Draw("same");
legy = new TLegend(0.15,0.8,0.3,0.88);
legy->AddEntry(comp_fry_raw_adc[0],"ENGINE","l");
legy->AddEntry(comp_fry_raw_adc[1],"HCANA","l");
legy->Draw();
cy->cd(2);
raw_adc_diff[1]->Add(comp_fry_raw_adc[1],comp_fry_raw_adc[0],-1);
raw_adc_diff[1]->Draw();
raw_adc_diff[1]->SetLineColor(kGreen-2);
cy->cd(3);
comp_fry_adc[0]->Draw();
comp_fry_adc[0]->SetLineColor(2);
comp_fry_adc[1]->Draw("same");
cy->cd(4);
adc_diff[1]->Add(comp_fry_adc[1],comp_fry_adc[0],-1);
adc_diff[1]->Draw();
adc_diff[1]->SetLineColor(kGreen-2);
cy->cd(5);
comp_fry[0]->Draw();
comp_fry[0]->SetLineColor(2);
comp_fry[1]->Draw("same");
cy->cd(6);
pos_diff[1]->Add(comp_fry[1],comp_fry[0],-1);
pos_diff[1]->Draw();
pos_diff[1]->SetLineColor(kGreen-2);
}
bool // returns true if success
processStyleSection(FILE *fp,string& theline, bool& new_section)
{
vector<string> v_tokens;
string *sid = NULL;
TStyle *thestyle = NULL;
if (gl_verbose)
cout << "Processing style section" << endl;
new_section=false;
while (getLine(fp,theline,"style")) {
if (!theline.size()) continue;
if (theline[0] == '#') continue; // comments are welcome
if (theline[0] == '[') {
new_section=true;
return true;
}
string key, value;
if (!getKeyValue(theline,key,value)) continue;
if (key == "id") {
if (sid != NULL) {
cerr << "no more than one id per style section allowed " << value << endl;
break;
}
sid = new string(value);
thestyle = new TStyle(*gStyle); // Assume current attributes, let user override specifics
thestyle->SetNameTitle(sid->c_str(),sid->c_str());
glmap_id2style.insert(pair<string,TStyle*>(*sid,thestyle));
continue;
}
else if (!sid) {
// assume the style is the global style
thestyle = gStyle;
sid = new string("using global"); // so as to pass this check the next time
cerr << "No style ID defined, assuming global style" << endl;
continue;
}
if (key == "style") {
if (value == "TDR")
setTDRStyle();
else if (value == "Plain")
gROOT->SetStyle("Plain");
else
cerr << "unknown style name " << value << endl;
}
else if (key == "optstat") {
if (gl_verbose) cout << "OptStat = " << thestyle->GetOptStat() << endl;
thestyle->SetOptStat(value.c_str());
if (gl_verbose) cout << "OptStat = " << thestyle->GetOptStat() << endl;
}
else if (key == "opttitle") thestyle->SetOptTitle(str2int(value));
else if (key == "padrightmargin") thestyle->SetPadRightMargin (str2flt(value));
else if (key == "padleftmargin") thestyle->SetPadLeftMargin (str2flt(value));
else if (key == "padtopmargin") thestyle->SetPadTopMargin (str2flt(value));
else if (key == "padbottommargin") thestyle->SetPadBottomMargin(str2flt(value));
else if (key == "padgridx") thestyle->SetPadGridX (str2int(value));
else if (key == "padgridy") thestyle->SetPadGridY (str2int(value));
// Set the position/size of the title box
else if (key == "titlexndc") thestyle->SetTitleX(str2flt(value));
else if (key == "titleyndc") thestyle->SetTitleY(str2flt(value));
else if (key == "titlewndc") thestyle->SetTitleW(str2flt(value));
else if (key == "titlehndc") thestyle->SetTitleH(str2flt(value));
else if (key == "titlefont") thestyle->SetTitleFont(str2int(value));
else if (key == "titlebordersize") thestyle->SetTitleBorderSize(str2int(value));
else if (key == "markercolor") thestyle->SetMarkerColor(str2int(value));
else if (key == "markerstyle") thestyle->SetMarkerStyle(str2int(value));
else if (key == "markersize") thestyle->SetMarkerSize(str2int(value));
else if (key == "linecolor") thestyle->SetLineColor(str2int(value));
else if (key == "linestyle") thestyle->SetLineStyle(str2int(value));
else if (key == "linewidth") thestyle->SetLineWidth(str2int(value));
else if (key == "fillcolor") thestyle->SetFillColor(str2int(value));
else if (key == "fillstyle") thestyle->SetFillStyle(str2int(value));
// axes
else if (key == "xtitlesize") thestyle->SetTitleSize(str2flt(value),"X");
else if (key == "ytitlesize") thestyle->SetTitleSize(str2flt(value),"Y");
else if (key == "ztitlesize") thestyle->SetTitleSize(str2flt(value),"Z");
else if (key == "xtitleoffset") thestyle->SetTitleOffset(str2flt(value),"X");
else if (key == "ytitleoffset") thestyle->SetTitleOffset(str2flt(value),"Y");
else if (key == "ztitleoffset") thestyle->SetTitleOffset(str2flt(value),"Z");
else if (key == "xlabeloffset") thestyle->SetLabelOffset(str2flt(value),"X");
else if (key == "ylabeloffset") thestyle->SetLabelOffset(str2flt(value),"Y");
else if (key == "zlabeloffset") thestyle->SetLabelOffset(str2flt(value),"Z");
else if (key == "xtitlefont") thestyle->SetTitleFont(str2int(value),"X");
else if (key == "ytitlefont") thestyle->SetTitleFont(str2int(value),"Y");
else if (key == "ztitlefont") thestyle->SetTitleFont(str2int(value),"Z");
else if (key == "xlabelsize") thestyle->SetLabelSize(str2flt(value),"X");
else if (key == "ylabelsize") thestyle->SetLabelSize(str2flt(value),"Y");
else if (key == "zlabelsize") thestyle->SetLabelSize(str2flt(value),"Z");
else if (key == "xlabelfont") thestyle->SetLabelFont(str2int(value),"X");
else if (key == "ylabelfont") thestyle->SetLabelFont(str2int(value),"Y");
else if (key == "zlabelfont") thestyle->SetLabelFont(str2int(value),"Z");
else if (key == "xndiv") thestyle->SetNdivisions(str2int(value),"X");
else if (key == "yndiv") thestyle->SetNdivisions(str2int(value),"Y");
else if (key == "zndiv") thestyle->SetNdivisions(str2int(value),"Z");
// Set the position of the statbox
else if (key == "statx2ndc") thestyle->SetStatX(str2flt(value));
else if (key == "staty2ndc") thestyle->SetStatY(str2flt(value));
else if (key == "statwndc") thestyle->SetStatW(str2flt(value));
else if (key == "stathndc") thestyle->SetStatH(str2flt(value));
else if (key == "statfont") thestyle->SetStatFont (str2int(value));
else if (key == "statfontsize") thestyle->SetStatFontSize(str2flt(value));
else if (key == "statformat") thestyle->SetStatFormat (value.c_str());
else if (key == "statstyle") thestyle->SetStatStyle (str2int(value.c_str()));
else if (key == "fitformat") thestyle->SetFitFormat (value.c_str());
else if (key == "painttextfmt") thestyle->SetPaintTextFormat(value.c_str());
else if (key == "markersize") thestyle->SetMarkerSize(str2flt(value));
else if (key == "palette") thestyle->SetPalette(str2int(value));
else {
cerr << "Unknown key " << key << endl;
}
} // while loop
return true;
} // processStyleSection
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 setTDRStyle() {
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");
tdrStyle->SetCanvasColor (0);
tdrStyle->SetCanvasBorderSize(10);
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetCanvasDefH (700);
tdrStyle->SetCanvasDefW (700);
tdrStyle->SetCanvasDefX (100);
tdrStyle->SetCanvasDefY (100);
// color palette for 2D temperature plots
tdrStyle->SetPalette(1,0);
// Pads
tdrStyle->SetPadColor (0);
tdrStyle->SetPadBorderSize (5);
tdrStyle->SetPadBorderMode (0);
tdrStyle->SetPadBorderMode (0);
tdrStyle->SetPadBottomMargin(0.15);
tdrStyle->SetPadTopMargin (0.08);
tdrStyle->SetPadLeftMargin (0.12);
tdrStyle->SetPadRightMargin (0.07);
tdrStyle->SetPadGridX (0);
tdrStyle->SetPadGridY (0);
tdrStyle->SetPadTickX (1);
tdrStyle->SetPadTickY (1);
// Frames
tdrStyle->SetLineWidth(3);
tdrStyle->SetFrameFillStyle ( 0);
tdrStyle->SetFrameFillColor ( 0);
tdrStyle->SetFrameLineColor ( 1);
tdrStyle->SetFrameLineStyle ( 0);
tdrStyle->SetFrameLineWidth ( 2);
tdrStyle->SetFrameBorderSize(5);
tdrStyle->SetFrameBorderMode( 0);
tdrStyle->SetFillStyle( 0);
// Histograms
tdrStyle->SetHistFillColor(2);
tdrStyle->SetHistFillStyle(0);
tdrStyle->SetHistLineColor(1);
tdrStyle->SetHistLineStyle(0);
tdrStyle->SetHistLineWidth(4);
tdrStyle->SetNdivisions(510);
// Functions
tdrStyle->SetFuncColor(1);
tdrStyle->SetFuncStyle(0);
tdrStyle->SetFuncWidth(2);
// Various
tdrStyle->SetMarkerStyle(20);
tdrStyle->SetMarkerColor(kBlack);
tdrStyle->SetMarkerSize (1.4);
tdrStyle->SetTitleBorderSize(0);
tdrStyle->SetTitleFillColor (0);
tdrStyle->SetTitleX (0.2);
tdrStyle->SetTitleSize (0.055,"X");
tdrStyle->SetTitleOffset(1.200,"X");
tdrStyle->SetLabelOffset(0.005,"X");
tdrStyle->SetLabelSize (0.050,"X");
tdrStyle->SetLabelFont (42 ,"X");
tdrStyle->SetStripDecimals(kFALSE);
tdrStyle->SetTitleSize (0.055,"Y");
tdrStyle->SetTitleOffset(1.000,"Y");
tdrStyle->SetLabelOffset(0.010,"Y");
tdrStyle->SetLabelSize (0.050,"Y");
tdrStyle->SetLabelFont (42 ,"Y");
tdrStyle->SetTextSize (0.055);
tdrStyle->SetTextFont (42);
tdrStyle->SetStatFont (42);
tdrStyle->SetTitleFont (42);
tdrStyle->SetTitleFont (42,"X");
tdrStyle->SetTitleFont (42,"Y");
tdrStyle->SetOptStat (0);
tdrStyle->SetOptTitle (0);
tdrStyle->SetLineStyleString(11,"20 10");
tdrStyle->SetLegendBorderSize(0);
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 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 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 setHHStyle(TStyle& HHStyle)
{
const int fontstyle=42;
HHStyle.SetPalette(1);
// ==============
// Canvas
// ==============
HHStyle.SetCanvasBorderMode(0);
HHStyle.SetCanvasColor(kWhite);
HHStyle.SetCanvasDefH(600); //Height of canvas
HHStyle.SetCanvasDefW(600); //Width of canvas
HHStyle.SetCanvasDefX(0); //Position on screen
HHStyle.SetCanvasDefY(0);
// ==============
// Pad
// ==============
HHStyle.SetPadBorderMode(0);
// HHStyle.SetPadBorderSize(Width_t size = 1);
HHStyle.SetPadColor(kWhite);
HHStyle.SetPadGridX(false);
HHStyle.SetPadGridY(false);
HHStyle.SetGridColor(kGray);
HHStyle.SetGridStyle(3);
HHStyle.SetGridWidth(1);
// ==============
// Frame
// ==============
HHStyle.SetFrameBorderMode(0);
HHStyle.SetFrameBorderSize(1);
HHStyle.SetFrameFillColor(0);
HHStyle.SetFrameFillStyle(0);
HHStyle.SetFrameLineColor(1);
HHStyle.SetFrameLineStyle(1);
HHStyle.SetFrameLineWidth(1);
// ==============
// Histo
// ==============
HHStyle.SetErrorX(0.0);
HHStyle.SetEndErrorSize(8);
// HHStyle.SetHistFillColor(1);
// HHStyle.SetHistFillStyle(0);
// HHStyle.SetHistLineColor(1);
HHStyle.SetHistLineStyle(0);
HHStyle.SetHistLineWidth(1);
// HHStyle.SetLegoInnerR(Float_t rad = 0.5);
// HHStyle.SetNumberContours(Int_t number = 20);
// HHStyle.SetErrorMarker(20);
HHStyle.SetMarkerStyle(20);
// ==============
// Fit/function
// ==============
HHStyle.SetOptFit(0);
HHStyle.SetFitFormat("5.4g");
HHStyle.SetFuncColor(2);
HHStyle.SetFuncStyle(1);
HHStyle.SetFuncWidth(1);
// ==============
// Date
// ==============
HHStyle.SetOptDate(0);
// HHStyle.SetDateX(Float_t x = 0.01);
// HHStyle.SetDateY(Float_t y = 0.01);
// =====================
// Statistics Box
// =====================
HHStyle.SetOptFile(0);
HHStyle.SetOptStat(0); // To display the mean and RMS: SetOptStat("mr");
HHStyle.SetStatColor(kWhite);
HHStyle.SetStatFont(fontstyle);
HHStyle.SetStatFontSize(0.025);
HHStyle.SetStatTextColor(1);
HHStyle.SetStatFormat("6.4g");
HHStyle.SetStatBorderSize(1);
HHStyle.SetStatH(0.1);
HHStyle.SetStatW(0.15);
// HHStyle.SetStatStyle(Style_t style = 1001);
// HHStyle.SetStatX(Float_t x = 0);
// HHStyle.SetStatY(Float_t y = 0);
// ==============
// Margins
// ==============
HHStyle.SetPadTopMargin(0.1);
HHStyle.SetPadBottomMargin(0.15);
HHStyle.SetPadLeftMargin(0.20);
HHStyle.SetPadRightMargin(0.05);
// ==============
// Global Title
// ==============
HHStyle.SetOptTitle(0);
HHStyle.SetTitleFont(fontstyle);
HHStyle.SetTitleColor(1);
HHStyle.SetTitleTextColor(1);
HHStyle.SetTitleFillColor(10);
HHStyle.SetTitleFontSize(0.05);
// HHStyle.SetTitleH(0); // Set the height of the title box
// HHStyle.SetTitleW(0); // Set the width of the title box
// HHStyle.SetTitleX(0); // Set the position of the title box
// HHStyle.SetTitleY(0.985); // Set the position of the title box
// HHStyle.SetTitleStyle(Style_t style = 1001);
// HHStyle.SetTitleBorderSize(2);
// ==============
// Axis titles
// ==============
HHStyle.SetTitleColor(1, "XYZ");
HHStyle.SetTitleFont(fontstyle, "XYZ");
HHStyle.SetTitleSize(0.05, "XYZ");
// HHStyle.SetTitleXSize(Float_t size = 0.02); // Another way to set the size?
// HHStyle.SetTitleYSize(Float_t size = 0.02);
HHStyle.SetTitleXOffset(1.0);
HHStyle.SetTitleYOffset(1.7);
// HHStyle.SetTitleOffset(1.1, "Y"); // Another way to set the Offset
// ==============
// Axis Label
// ==============
//HHStyle.SetLabelColor(1, "XYZ");
HHStyle.SetLabelFont(fontstyle, "XYZ");
HHStyle.SetLabelOffset(0.007, "XYZ");
HHStyle.SetLabelSize(0.04, "XYZ");
// ==============
// Axis
// ==============
HHStyle.SetAxisColor(1, "XYZ");
HHStyle.SetStripDecimals(kTRUE);
HHStyle.SetTickLength(0.03, "XYZ");
HHStyle.SetNdivisions(510, "XYZ");
HHStyle.SetPadTickX(1); // To get tick marks on the opposite side of the frame
HHStyle.SetPadTickY(1);
// Change for log plots:
HHStyle.SetOptLogx(0);
HHStyle.SetOptLogy(0);
HHStyle.SetOptLogz(0);
// ==============
// Text
// ==============
HHStyle.SetTextAlign(11);
HHStyle.SetTextAngle(0);
HHStyle.SetTextColor(1);
HHStyle.SetTextFont(fontstyle);
HHStyle.SetTextSize(0.05);
// =====================
// Postscript options:
// =====================
HHStyle.SetPaperSize(20.,20.);
// HHStyle.SetLineScalePS(Float_t scale = 3);
// HHStyle.SetLineStyleString(Int_t i, const char* text);
// HHStyle.SetHeaderPS(const char* header);
// HHStyle.SetTitlePS(const char* pstitle);
// HHStyle.SetBarOffset(Float_t baroff = 0.5);
// HHStyle.SetBarWidth(Float_t barwidth = 0.5);
// HHStyle.SetPaintTextFormat(const char* format = "g");
// HHStyle.SetPalette(Int_t ncolors = 0, Int_t* colors = 0);
// HHStyle.SetTimeOffset(Double_t toffset);
// HHStyle.SetHistMinimumZero(kTRUE);
}
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 TDRStyle() {
TStyle *TDRStyle = new TStyle("TDRStyle", "TDR style");
/// helvetica, size absolute
// TDRStyle->font = 43
// TDRStyle->titleSize = 33
// TDRStyle->labelSize = 27
// TDRStyle->statSize = 14
/// 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.5); // needed for uncertainty histograms
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(43);
TDRStyle->SetStatFontSize(14);
TDRStyle->SetStatTextColor(1);
TDRStyle->SetStatFormat("6.4g");
TDRStyle->SetStatBorderSize(1);
TDRStyle->SetStatH(0.1);
TDRStyle->SetStatW(0.15);
/// Margins:
TDRStyle->SetPadTopMargin(0.05); // default
TDRStyle->SetPadBottomMargin(0.13);
TDRStyle->SetPadLeftMargin(0.16);
TDRStyle->SetPadRightMargin(0.05); // default
/// For the Global title:
TDRStyle->SetOptTitle(0);
TDRStyle->SetTitleFont(43);
TDRStyle->SetTitleColor(1);
TDRStyle->SetTitleTextColor(1);
TDRStyle->SetTitleFillColor(10);
TDRStyle->SetTitleFontSize(27);
/// For the axis titles:
TDRStyle->SetTitleColor(1, "XYZ");
TDRStyle->SetTitleFont(43, "XYZ");
TDRStyle->SetTitleSize(33, "XYZ");
TDRStyle->SetTitleXOffset(0.9);
TDRStyle->SetTitleYOffset(1.25);
/// For the axis labels:
TDRStyle->SetLabelColor(1, "XYZ");
TDRStyle->SetLabelFont(43, "XYZ");
TDRStyle->SetLabelOffset(0.007, "XYZ");
TDRStyle->SetLabelSize(27, "XYZ");
/// For the axis:
TDRStyle->SetAxisColor(1, "XYZ");
TDRStyle->SetStripDecimals(true);
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);
/// For the text boxes
TDRStyle->SetTextFont(43+20); // Bold
TDRStyle->SetTextSize(27);
/// Change for log plots:
TDRStyle->SetOptLogx(0);
TDRStyle->SetOptLogy(0);
TDRStyle->SetOptLogz(0);
/// Postscript options:
TDRStyle->SetPaperSize(20,20);
TDRStyle->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 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();
}
