void hggPaperStyle()
{
// Hgg Paper style for plots
TStyle *hggStyle = new TStyle("hggStyle","Hgg Paper Style");
//hggStyle->SetCanvasColor(0);
//hggStyle->SetPadColor(0);
hggStyle->SetPadTickX(0);
hggStyle->SetPadTickY(0);
hggStyle->SetFrameFillColor(0);
hggStyle->SetStatColor(0);
hggStyle->SetOptStat(0);
hggStyle->SetTitleFillColor(0);
hggStyle->SetCanvasBorderMode(0);
hggStyle->SetPadBorderMode(0);
hggStyle->SetFrameBorderMode(0);
hggStyle->SetFrameBorderSize(1);
hggStyle->SetPadColor(kWhite);
hggStyle->SetCanvasColor(kWhite);
hggStyle->SetCanvasDefH(600); //Height of canvas
hggStyle->SetCanvasDefW(800); //Width of canvas
hggStyle->SetCanvasDefX(0); //POsition on screen
hggStyle->SetCanvasDefY(0);
hggStyle->SetPadLeftMargin(0.16);//0.16);
hggStyle->SetPadRightMargin(0.1);//0.02);
hggStyle->SetPadTopMargin(0.085);//0.02);
hggStyle->SetPadBottomMargin(0.12);//0.02);
// For hgg axis titles:
hggStyle->SetTitleColor(1, "XYZ");
hggStyle->SetTitleFont(42, "XYZ");
hggStyle->SetTitleSize(0.05, "XYZ");
hggStyle->SetTitleYOffset(1.5); // => 1.15 if exponents
hggStyle->SetTitleXOffset(1.);//0.9);
// For hgg axis labels:
hggStyle->SetLabelColor(1, "XYZ");
hggStyle->SetLabelFont(42, "XYZ");
hggStyle->SetLabelOffset(0.007, "XYZ");
hggStyle->SetLabelSize(0.045, "XYZ");
// Legends
hggStyle->SetLegendBorderSize(0);
hggStyle->SetLegendFillColor(kWhite);
hggStyle->SetLegendFont(42);
hggStyle->SetFillColor(10);
// Nothing for now
hggStyle->SetTextFont(42);
hggStyle->SetTextSize(0.03);
hggStyle->SetOptTitle(0);
hggStyle->SetOptStat(0);
hggStyle->cd();
}
void setNiceStyle()
{
gROOT->SetStyle("Plain");
// --------------------------------------------------
// define proper canvas style
// --------------------------------------------------
TStyle *MyStyle = new TStyle ("MyStyle", "My style for nicer plots");
Float_t xoff = MyStyle->GetLabelOffset("X"),
yoff = MyStyle->GetLabelOffset("Y"),
zoff = MyStyle->GetLabelOffset("Z");
MyStyle->SetCanvasBorderMode ( 0 );
MyStyle->SetPadBorderMode ( 0 );
MyStyle->SetPadColor ( 0 );
MyStyle->SetCanvasColor ( 0 );
MyStyle->SetTitleColor ( 0 );
MyStyle->SetStatColor ( 0 );
MyStyle->SetTitleBorderSize ( 0 );
MyStyle->SetTitleFillColor ( 0 );
MyStyle->SetTitleH ( 0.07 );
MyStyle->SetTitleW ( 1.00 );
MyStyle->SetTitleFont ( 132 );
MyStyle->SetLabelOffset (1.5*xoff, "X");
MyStyle->SetLabelOffset (1.5*yoff, "Y");
MyStyle->SetLabelOffset (1.5*zoff, "Z");
MyStyle->SetTitleOffset (0.9, "X");
MyStyle->SetTitleOffset (0.9, "Y");
MyStyle->SetTitleOffset (0.9, "Z");
MyStyle->SetTitleSize (0.045, "X");
MyStyle->SetTitleSize (0.045, "Y");
MyStyle->SetTitleSize (0.045, "Z");
MyStyle->SetLabelFont (132, "X");
MyStyle->SetLabelFont (132, "Y");
MyStyle->SetLabelFont (132, "Z");
MyStyle->SetPalette(1);
MyStyle->cd();
}
// 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 PlotAlignmentValidation::setNiceStyle() {
TStyle *MyStyle = new TStyle ("MyStyle", "My style for nicer plots");
Float_t xoff = MyStyle->GetLabelOffset("X"),
yoff = MyStyle->GetLabelOffset("Y"),
zoff = MyStyle->GetLabelOffset("Z");
MyStyle->SetCanvasBorderMode ( 0 );
MyStyle->SetPadBorderMode ( 0 );
MyStyle->SetPadColor ( 0 );
MyStyle->SetCanvasColor ( 0 );
MyStyle->SetTitleColor ( 0 );
MyStyle->SetStatColor ( 0 );
MyStyle->SetTitleBorderSize ( 0 );
MyStyle->SetTitleFillColor ( 0 );
MyStyle->SetTitleH ( 0.07 );
MyStyle->SetTitleW ( 1.00 );
MyStyle->SetTitleFont ( 132 );
MyStyle->SetLabelOffset (1.5*xoff, "X");
MyStyle->SetLabelOffset (1.5*yoff, "Y");
MyStyle->SetLabelOffset (1.5*zoff, "Z");
MyStyle->SetTitleOffset (1.5, "X");
MyStyle->SetTitleOffset (1.2, "Y");
MyStyle->SetTitleOffset (0.9, "Z");
MyStyle->SetTitleSize (0.045, "X");
MyStyle->SetTitleSize (0.045, "Y");
MyStyle->SetTitleSize (0.045, "Z");
MyStyle->SetLabelFont (132, "X");
MyStyle->SetLabelFont (132, "Y");
MyStyle->SetLabelFont (132, "Z");
MyStyle->SetPalette(1);
MyStyle->cd();
}
// 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 *setStyle = new TStyle("setStyle","Style for GCT Analyzer");
// Stuff from plain style
setStyle->SetCanvasColor(kWhite);
setStyle->SetFrameBorderMode(0);
setStyle->SetCanvasBorderMode(0);
setStyle->SetFrameFillColor(0);
setStyle->SetPadBorderMode(0);
setStyle->SetPadColor(kWhite);
setStyle->SetOptStat(0);
setStyle->SetOptTitle(0);
setStyle->SetOptFit(0);
setStyle->SetOptDate(0);
// Labels and borders
setStyle->SetLabelSize(0.055,"x");
setStyle->SetLabelSize(0.06,"y");
setStyle->SetLabelOffset(0.00,"x");
setStyle->SetLabelOffset(0.00,"y");
setStyle->SetTitleOffset(0.05,"x");
setStyle->SetTitleOffset(0.50,"y");
//setStyle->SetLabelFont(22,"x");
//setStyle->SetLabelFont(22,"y");
//setStyle->SetErrorX(0.0000);
//setStyle->SetTickLength(0.05,"x");
//setStyle->SetTickLength(0.05,"y");
//setStyle->SetLineWidth(0.8);
//setStyle->SetPadTickX(1);
//setStyle->SetPadTickY(1);
setStyle->SetPadLeftMargin(0.15);
setStyle->SetPadTopMargin(0.077);
setStyle->SetPadBottomMargin(0.123);
setStyle->cd();
}
void blinding_study()
{
gSystem->CompileMacro("MitGPTree.h");
// First we define MIT Style for the plots.
TStyle *MitStyle = gStyle;
//gStyle = MitStyle;
// Canvas
MitStyle->SetCanvasColor (0);
MitStyle->SetCanvasBorderSize(10);
MitStyle->SetCanvasBorderMode(0);
MitStyle->SetCanvasDefH (700);
MitStyle->SetCanvasDefW (700);
MitStyle->SetCanvasDefX (100);
MitStyle->SetCanvasDefY (100);
// Pads
MitStyle->SetPadColor (0);
MitStyle->SetPadBorderSize (10);
MitStyle->SetPadBorderMode (0);
MitStyle->SetPadBottomMargin(0.13);
MitStyle->SetPadTopMargin (0.04);
MitStyle->SetPadLeftMargin (0.18);
MitStyle->SetPadRightMargin (0.04);
MitStyle->SetPadGridX (0);
MitStyle->SetPadGridY (0);
MitStyle->SetPadTickX (0);
MitStyle->SetPadTickY (0);
// Frames
MitStyle->SetFrameFillStyle ( 0);
MitStyle->SetFrameFillColor ( 0);
MitStyle->SetFrameLineColor ( 1);
MitStyle->SetFrameLineStyle ( 0);
MitStyle->SetFrameLineWidth ( 1);
MitStyle->SetFrameBorderSize(10);
MitStyle->SetFrameBorderMode( 0);
// Histograms
MitStyle->SetHistFillColor(2);
MitStyle->SetHistFillStyle(0);
MitStyle->SetHistLineColor(1);
MitStyle->SetHistLineStyle(0);
MitStyle->SetHistLineWidth(2);
MitStyle->SetNdivisions(505);
// Functions
MitStyle->SetFuncColor(1);
MitStyle->SetFuncStyle(0);
MitStyle->SetFuncWidth(2);
// Various
MitStyle->SetMarkerStyle(20);
MitStyle->SetMarkerColor(kBlack);
MitStyle->SetMarkerSize (1.2);
MitStyle->SetTitleSize (0.055,"X");
MitStyle->SetTitleOffset(1.200,"X");
MitStyle->SetLabelOffset(0.005,"X");
MitStyle->SetLabelSize (0.050,"X");
MitStyle->SetLabelFont (42 ,"X");
MitStyle->SetTickLength (-0.03,"X");
MitStyle->SetStripDecimals(kFALSE);
MitStyle->SetTitleSize (0.055,"Y");
MitStyle->SetTitleOffset(1.800,"Y");
MitStyle->SetLabelOffset(0.010,"Y");
MitStyle->SetLabelSize (0.050,"Y");
MitStyle->SetLabelFont (42 ,"Y");
MitStyle->SetTickLength (-0.03,"Y");
MitStyle->SetTextSize (0.055);
MitStyle->SetTextFont (42);
MitStyle->SetStatFont (42);
MitStyle->SetTitleFont (42);
MitStyle->SetTitleFont (42,"X");
MitStyle->SetTitleFont (42,"Y");
MitStyle->SetOptStat (0);
// Here the style section ends and the macro begins.
string sig_samples[] =
{
"s12-dmmpho-v_m1-v7a",
"s12-dmmpho-av_m1-v7a",
"s12-dmmpho-v_m10-v7a",
"s12-dmmpho-av_m10-v7a",
"s12-dmmpho-v_m100-v7a",
"s12-dmmpho-av_m100-v7a",
"s12-dmmpho-v_m200-v7a",
"s12-dmmpho-av_m200-v7a",
"s12-dmmpho-av_m300-v7a",
"s12-dmmpho-v_m500-v7a",
"s12-dmmpho-av_m500-v7a",
"s12-dmmpho-v_m1000-v7a",
"s12-dmmpho-av_m1000-v7a",
"s12-addmpho-md1_d2-v7a",
"s12-addmpho-md1_d3-v7a",
"s12-addmpho-md1_d4-v7a",
"s12-addmpho-md1_d5-v7a",
"s12-addmpho-md1_d6-v7a",
"s12-addmpho-md2_d2-v7a",
"s12-addmpho-md2_d3-v7a",
"s12-addmpho-md2_d5-v7a",
"s12-addmpho-md2_d6-v7a",
"s12-addmpho-md3_d2-v7a",
"s12-addmpho-md3_d3-v7a",
"s12-addmpho-md3_d4-v7a",
"s12-addmpho-md3_d5-v7a",
"s12-addmpho-md3_d6-v7a"
}
Int_t signal_num=0;
// This macro considers only one signal sample at a time, so the variable signal_num defines
// which signal sample to work on (it is the index of the sample in sig_samples).
double sig_weights[] =
{
4.81E-07,
4.79E-07,
4.80E-07,
4.81E-07,
4.77E-07,
4.26E-07,
4.23E-07,
3.18E-07,
2.19E-07,
2.01E-07,
9.52E-08,
3.00E-08,
8.22E-09,
6.48E-02,
1.73E-01,
6.93E-02,
9.17E-02,
9.63E-01,
5.87E-03,
4.91E-03,
4.37E-03,
4.26E-03,
1.45E-03,
7.93E-04,
5.53E-04,
4.26E-04,
3.24E-04
};
// The signal weights are given by (sigma_MC * lumi)/(N_processed). Background weights are defined similarly below.
string bg_samples[] = {"s12-zgptg130-v7c", "s12-wjets-ptw100-v7a", "s12-wgptg130-v7a", "s12-qcdht100-250-v7a", "s12-qcdht250-500-v7a", "s12-qcdht500-1000-v7a",
"s12-qcdht1000-v7a", "s12-2pibo10_25-v7a", "s12-2pibo25_250-v7a", "s12-2pibo250-v7a", "s12-2pibx10_25-v7a", "s12-2pibx25_250-v7a", "s12-2pibx250-v7a",
"s12-pj50_80-v7a", "s12-pj80_120-v7a", "s12-pj120_170-v7a", "s12-pj170_300-v7a", "s12-pj300_470-v7a", "s12-pj470_800-v7a", "s12-pj800_1400-v7a",
"s12-pj1400_1800-v7a", "s12-pj1800-v7a", "s12-zgllgptg130-v7a", "s12-zgllg-v7a"};
double bg_weights[]=
{
5.28E-03,
1.14E-01,
1.38E-02,
5.63E+03,
2.34E+02,
9.71E+00,
3.09E-01,
9.21E+00,
9.90E-01,
4.21E-04,
1.66E+01,
6.43E-01,
4.91E-05,
3.30E+01,
5.49E+00,
1.09E+00,
3.01E-01,
2.15E-02,
2.10E-03,
7.11E-05,
4.45E-07,
1.88E-08,
3.18E-03,
5.10E-01
};
string sigLine = sig_samples[signal_num];
Double_t sigWeight = sig_weights[signal_num];
TH1* sig_hist = new TH1F("sig_hist", TString("Phi Between Met and Photon for Signal Events [") + TString(sigLine) + TString("]"), 70, 0, 3.5);
TH1* bg_hist = new TH1F("bg_hist", TString("Phi Between Met and Photon for Background Events [") + TString(sigLine) + TString("]"), 70, 0, 3.5);
TH1* signifhist = new TH1F("signifhist", TString(" "), 24, 2, 3.2);
cout<<"The selected signal sample is: "<<TString(sigLine)<<" "<<endl;
TString sigFilename = TString("monoph-2013-July9_") + TString(sigLine) + TString("_noskim.root");
MitGPTree sigEvent;
sigEvent.LoadTree(TString("/scratch/cferko/hist/monoph-2013-July9/merged/")+sigFilename, 0);
sigEvent.InitTree(0);
int nDataSig=sigEvent.tree_->GetEntries();
for (int evt=0; evt<nDataSig; ++evt)
{
sigEvent.tree_->GetEntry(evt);
Double_t pho1Pt = sigEvent.pho1_.Pt();
Double_t jet1Pt = sigEvent.jet1_.Pt();
Double_t met = sigEvent.met_;
Double_t metPhi = sigEvent.metPhi_;
Double_t nphotons = sigEvent.nphotons_;
Double_t ncosmics = sigEvent.ncosmics_;
Double_t phoPassEleVeto = sigEvent.phoPassEleVeto_a1_;
Double_t jet1Pt = sigEvent.jet1_.Pt();
Double_t nlep = sigEvent.nlep_;
Double_t pho1Eta = sigEvent.pho1_.Eta();
Double_t pho1Phi = sigEvent.pho1_.Phi();
Double_t phoIsTrigger = sigEvent.phoIsTrigger_a1_;
Double_t phoLeadTimeSpan = sigEvent.phoLeadTimeSpan_a1_;
Double_t phoCoviEtaiEta = sigEvent.phoCoviEtaiEta_a1_;
Double_t phoCoviPhiiPhi = sigEvent.phoCoviPhiiPhi_a1_;
Double_t phoMipIsHalo = sigEvent.phoMipIsHalo_a1_;
Double_t phoSeedTime = sigEvent.phoSeedTime_a1_;
Double_t deltaPhi = TMath::ACos(TMath::Cos(pho1Phi - metPhi));
if (TMath::Abs(pho1Eta)<1.479 && met>140 && pho1Pt>160 && nphotons>0 && phoPassEleVeto>0 && phoIsTrigger==1 && TMath::Abs(phoLeadTimeSpan) < 8. && phoCoviEtaiEta > 0.001 && phoCoviPhiiPhi > 0.001 && phoMipIsHalo == 0 && phoSeedTime > -1.5 && nlep==0 && ncosmics==0 && jet1Pt<100) sig_hist->Fill(deltaPhi, sigWeight);
}
for (int i=0; i<24; i++)
{
string bgLine = bg_samples[i];
Double_t bgWeight = bg_weights[i];
TString bgFilename = TString("monoph-2013-July9_") + TString(bgLine) + TString("_noskim.root");
MitGPTree bgEvent;
bgEvent.LoadTree(TString("/scratch/cferko/hist/monoph-2013-July9/merged/") + bgFilename, 0);
bgEvent.InitTree(0);
int nDataBg=bgEvent.tree_->GetEntries();
for (int evt=0; evt<nDataBg; ++evt)
{
bgEvent.tree_->GetEntry(evt);
Double_t pho1Pt = bgEvent.pho1_.Pt();
Double_t jet1Pt = bgEvent.jet1_.Pt();
Double_t met = bgEvent.met_;
Double_t metPhi = bgEvent.metPhi_;
Double_t nphotons = bgEvent.nphotons_;
Double_t ncosmics = bgEvent.ncosmics_;
Double_t phoPassEleVeto = bgEvent.phoPassEleVeto_a1_;
Double_t jet1Pt = bgEvent.jet1_.Pt();
Double_t nlep = bgEvent.nlep_;
Double_t pho1Eta = bgEvent.pho1_.Eta();
Double_t pho1Phi = bgEvent.pho1_.Phi();
Double_t phoIsTrigger = bgEvent.phoIsTrigger_a1_;
Double_t phoLeadTimeSpan = bgEvent.phoLeadTimeSpan_a1_;
Double_t phoCoviEtaiEta = bgEvent.phoCoviEtaiEta_a1_;
Double_t phoCoviPhiiPhi = bgEvent.phoCoviPhiiPhi_a1_;
Double_t phoMipIsHalo = bgEvent.phoMipIsHalo_a1_;
Double_t phoSeedTime = bgEvent.phoSeedTime_a1_;
Double_t deltaPhi = TMath::ACos(TMath::Cos(pho1Phi - metPhi));
if (TMath::Abs(pho1Eta)<1.479 && met>140 && pho1Pt>160 && nphotons>0 && phoPassEleVeto>0 && phoIsTrigger==1 && TMath::Abs(phoLeadTimeSpan) < 8. && phoCoviEtaiEta > 0.001 && phoCoviPhiiPhi > 0.001 && phoMipIsHalo == 0 && phoSeedTime > -1.5 && nlep==0 && ncosmics==0 && jet1Pt<100) bg_hist->Fill(deltaPhi, bgWeight);
}
}
Double_t sigTot = sig_hist->Integral();
Double_t bgTot = bg_hist->Integral();
Double_t signifTot = sigTot/TMath::Sqrt(bgTot);
Double_t signifTarget = 0.25*signifTot;
Int_t reduced=0;
// The strange limits in the for loop are a messy hack to get the integral to work out. Since hist->Integral(a,b)
// integrates from BIN a to BIN b (bin numbers instead of values of the x-axis), I choose to work from bin 64
// (an angle cut of 3.2) to bin 40 (an angle cut of 2).
for (int i=64; i>=40; i--)
{
Double_t sigcount = sig_hist->Integral(0, i);
Double_t bgcount = bg_hist->Integral(0, i);
if (bgcount>0)
{
Double_t significance = sigcount/TMath::Sqrt(bgcount);
// cout<<"At a phi cut of "<<sig_hist->GetBinCenter(i)+0.025<<" the significance is "<<significance<<", which is a fraction "<<significance/signifTot<<" of max."<<endl;
if (significance<signifTarget && reduced==0)
{
cout<<"The signal significance is reduced to one-fourth of its uncut value at a phi cut of "<<sig_hist->GetBinCenter(i)+0.025<<endl;
reduced=1;
}
signifhist->Fill(sig_hist->GetBinCenter(i), significance);
}
}
signifhist->SetStats(kFALSE);
signifhist->GetXaxis()->SetTitle("Maximum #Delta#phi");
signifhist->GetYaxis()->SetTitle("Signal Significance");
signifhist->Draw();
c1->Print(TString("Blinding_") + TString(sigLine) + TString(".pdf"));
}
void SetStyle() {
TH1::SetDefaultSumw2();
gROOT->SetStyle("Plain");
TStyle *MITStyle = new TStyle("MIT-Style","The Perfect Style for Plots ;-)");
gStyle = MITStyle;
gStyle->SetTitleSize(0.08,"t");// gStyle->SetTitleX(0.8);
gStyle->SetTitleW(0.3);
gStyle->SetFrameFillColor(2);
gStyle->SetFrameLineColor(2);
// Canvas
MITStyle->SetCanvasColor (0);
MITStyle->SetCanvasBorderSize(10);
MITStyle->SetCanvasBorderMode(0);
MITStyle->SetCanvasDefH (500);
MITStyle->SetCanvasDefW (550);
// MITStyle->SetCanvasDefX (100);
// MITStyle->SetCanvasDefY (100);
// Pads
MITStyle->SetPadColor (0);
// MITStyle->SetPadBorderSize (10);
MITStyle->SetPadBorderMode (0);
MITStyle->SetPadBottomMargin(0.18);
MITStyle->SetPadTopMargin (0.06);
MITStyle->SetPadLeftMargin (0.22);
MITStyle->SetPadRightMargin (0.06);
MITStyle->SetPadGridX (0);
MITStyle->SetPadGridY (0);
MITStyle->SetPadTickX (1);
MITStyle->SetPadTickY (1);
// Frames
MITStyle->SetFrameFillStyle ( 0);
MITStyle->SetFrameFillColor ( 0);
MITStyle->SetFrameLineColor ( 1);
MITStyle->SetFrameLineStyle ( 0);
MITStyle->SetFrameLineWidth ( 1);
// MITStyle->SetFrameBorderSize(10);
MITStyle->SetFrameBorderMode( 0);
// Histograms
MITStyle->SetHistLineColor(1);
MITStyle->SetHistLineStyle(0);
MITStyle->SetHistLineWidth(3);
MITStyle->SetNdivisions(505,"X");
MITStyle->SetNdivisions(505,"Y");
MITStyle->SetNdivisions(505,"Z");
// Functions
MITStyle->SetFuncColor(1);
MITStyle->SetFuncStyle(0);
MITStyle->SetFuncWidth(2);
// Various
MITStyle->SetMarkerStyle(20);
MITStyle->SetMarkerColor(kBlack);
MITStyle->SetMarkerSize (1.2);
// MITStyle->SetTitleSize (0.1,"t");
MITStyle->SetTitleSize (0.070,"X");
MITStyle->SetTitleOffset(1.0,"X");
MITStyle->SetTitleFont (42,"X");
MITStyle->SetLabelOffset(0.005,"X");
MITStyle->SetLabelSize (0.070,"X");
MITStyle->SetLabelFont (42 ,"X");
MITStyle->SetTitleSize (0.070,"Y");
MITStyle->SetTitleOffset(1.0,"Y");
MITStyle->SetTitleFont (42,"Y");
MITStyle->SetLabelOffset(0.005,"Y");
MITStyle->SetLabelSize (0.070,"Y");
MITStyle->SetLabelFont (42 ,"Y");
MITStyle->SetTitleSize (0.06,"Z");
MITStyle->SetTitleOffset(1.640,"Z");
MITStyle->SetTitleFont (42,"Z");
MITStyle->SetLabelOffset(0.005,"Z");
MITStyle->SetLabelSize (0.070,"Z");
MITStyle->SetLabelFont (42 ,"Z");
MITStyle->SetStripDecimals(kFALSE);
MITStyle->SetTitleBorderSize(0);
MITStyle->SetTitleFillColor(0);
// MITStyle->SetTitleAlign(6);
MITStyle->SetTextSize (20);
MITStyle->SetTextFont (43);
MITStyle->SetStatFont (42);
MITStyle->SetOptStat (0);
MITStyle->SetLegendBorderSize(0);
MITStyle->SetEndErrorSize(0);
MITStyle->SetErrorX(0);
MITStyle->SetPalette (1,0);
}
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 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 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();
}
void WZPaperStyle()
{
TStyle* WZPaperStyle = new TStyle("WZPaperStyle", "WZPaperStyle");
gStyle = WZPaperStyle;
//----------------------------------------------------------------------------
// Canvas
//----------------------------------------------------------------------------
WZPaperStyle->SetCanvasBorderMode( 0);
WZPaperStyle->SetCanvasBorderSize( 10);
WZPaperStyle->SetCanvasColor ( 0);
WZPaperStyle->SetCanvasDefH (600);
WZPaperStyle->SetCanvasDefW (550);
WZPaperStyle->SetCanvasDefX ( 10);
WZPaperStyle->SetCanvasDefY ( 10);
//----------------------------------------------------------------------------
// Pad
//----------------------------------------------------------------------------
WZPaperStyle->SetPadBorderMode ( 0);
WZPaperStyle->SetPadBorderSize ( 10);
WZPaperStyle->SetPadColor ( 0);
WZPaperStyle->SetPadBottomMargin(0.20);
WZPaperStyle->SetPadTopMargin (0.08);
WZPaperStyle->SetPadLeftMargin (0.18);
WZPaperStyle->SetPadRightMargin (0.05);
//----------------------------------------------------------------------------
// Frame
//----------------------------------------------------------------------------
WZPaperStyle->SetFrameBorderMode( 0);
WZPaperStyle->SetFrameBorderSize(10);
WZPaperStyle->SetFrameFillStyle ( 0);
WZPaperStyle->SetFrameFillColor ( 0);
WZPaperStyle->SetFrameLineColor ( 1);
WZPaperStyle->SetFrameLineStyle ( 0);
WZPaperStyle->SetFrameLineWidth ( 2);
//----------------------------------------------------------------------------
// Hist
//----------------------------------------------------------------------------
WZPaperStyle->SetHistFillColor(0);
WZPaperStyle->SetHistFillStyle(1);
WZPaperStyle->SetHistLineColor(1);
WZPaperStyle->SetHistLineStyle(0);
WZPaperStyle->SetHistLineWidth(1);
//----------------------------------------------------------------------------
// Axis
//----------------------------------------------------------------------------
WZPaperStyle->SetLabelFont ( 42, "xyz");
WZPaperStyle->SetLabelOffset(0.015, "xyz");
WZPaperStyle->SetLabelSize (0.040, "xyz");
WZPaperStyle->SetNdivisions ( 505, "xyz");
WZPaperStyle->SetTitleFont ( 42, "xyz");
WZPaperStyle->SetTitleSize (0.040, "xyz");
WZPaperStyle->SetTitleOffset( 1.15, "x");
WZPaperStyle->SetTitleOffset( 1.15, "y");
WZPaperStyle->SetPadTickX ( 1);
WZPaperStyle->SetPadTickY ( 1);
//----------------------------------------------------------------------------
// Title
//----------------------------------------------------------------------------
WZPaperStyle->SetTitleAlign ( 33);
WZPaperStyle->SetTitleBorderSize( 0);
WZPaperStyle->SetTitleFillColor ( 10);
WZPaperStyle->SetTitleFont ( 42);
WZPaperStyle->SetTitleFontSize (0.045);
WZPaperStyle->SetTitleX (0.960);
WZPaperStyle->SetTitleY (0.990);
//----------------------------------------------------------------------------
// Stat
//----------------------------------------------------------------------------
WZPaperStyle->SetOptStat ( 0);
WZPaperStyle->SetStatBorderSize( 0);
WZPaperStyle->SetStatColor ( 10);
WZPaperStyle->SetStatFont ( 42);
WZPaperStyle->SetStatX (0.94);
WZPaperStyle->SetStatY (0.91);
}
int electronCompare()
{
TString CMP_DIR = gSystem->Getenv( "CMP_DIR" ) ;
TString CMP_URL = gSystem->Getenv( "CMP_URL" ) ;
TString CMP_TITLE = gSystem->Getenv( "CMP_TITLE" ) ;
TString CMP_RED_FILE = gSystem->Getenv( "CMP_RED_FILE" ) ;
TString CMP_BLUE_FILE = gSystem->Getenv( "CMP_BLUE_FILE" ) ;
TString CMP_RED_NAME = gSystem->Getenv( "CMP_RED_NAME" ) ;
TString CMP_BLUE_NAME = gSystem->Getenv( "CMP_BLUE_NAME" ) ;
TString CMP_RED_COMMENT = gSystem->Getenv( "CMP_RED_COMMENT" ) ;
TString CMP_BLUE_COMMENT = gSystem->Getenv( "CMP_BLUE_COMMENT" ) ;
TString CMP_CONFIG = gSystem->Getenv( "CMP_CONFIG" ) ;
TString CMP_RED_RELEASE = gSystem->Getenv( "CMP_RED_RELEASE" ) ;
TString CMP_BLUE_RELEASE = gSystem->Getenv( "CMP_BLUE_RELEASE" ) ;
//-----
// AC
// std::cout << "red_file : C : " << CMP_RED_FILE << std::endl;
// std::cout << "blue_file : C : " << CMP_BLUE_FILE << std::endl;
std::cout << "red_release : C : " << CMP_RED_RELEASE << std::endl;
std::cout << "blue_release : C : " << CMP_BLUE_RELEASE << std::endl;
//-----
// style:
TStyle *eleStyle = new TStyle("eleStyle","Style for electron validation");
eleStyle->SetCanvasBorderMode(0);
eleStyle->SetCanvasColor(kWhite);
eleStyle->SetCanvasDefH(600);
eleStyle->SetCanvasDefW(800);
eleStyle->SetCanvasDefX(0);
eleStyle->SetCanvasDefY(0);
eleStyle->SetPadBorderMode(0);
eleStyle->SetPadColor(kWhite);
eleStyle->SetPadGridX(false);
eleStyle->SetPadGridY(false);
eleStyle->SetGridColor(0);
eleStyle->SetGridStyle(3);
eleStyle->SetGridWidth(1);
eleStyle->SetOptStat(1);
eleStyle->SetPadTickX(1);
eleStyle->SetPadTickY(1);
eleStyle->SetHistLineColor(1);
eleStyle->SetHistLineStyle(0);
eleStyle->SetHistLineWidth(2);
eleStyle->SetEndErrorSize(2);
eleStyle->SetErrorX(0.);
eleStyle->SetOptStat(1);
eleStyle->SetTitleColor(1, "XYZ");
eleStyle->SetTitleFont(42, "XYZ");
eleStyle->SetTitleXOffset(1.0);
eleStyle->SetTitleYOffset(1.0);
eleStyle->SetLabelOffset(0.005, "XYZ");
eleStyle->SetTitleSize(0.05, "XYZ");
eleStyle->SetTitleFont(22,"X");
eleStyle->SetTitleFont(22,"Y");
eleStyle->SetHistLineWidth(2);
eleStyle->SetPadBottomMargin(0.13);
eleStyle->SetPadLeftMargin(0.15);
eleStyle->SetMarkerStyle(21);
eleStyle->SetMarkerSize(0.8);
//-- AC --
eleStyle->SetPadRightMargin(0.2) ;
eleStyle->cd();
gROOT->ForceStyle();
TString internal_path("DQMData/Run 1/EgammaV/Run summary/") ;
TString old_internal_path("DQMData/EgammaV/") ;
TString file_ref_dir ;
TFile * file_ref = 0 ;
if ( CMP_BLUE_FILE != "" )
{
file_ref = TFile::Open(CMP_BLUE_FILE) ;
if (file_ref!=0)
{
std::cout<<"open "<<CMP_BLUE_FILE<<std::endl ;
if (file_ref->cd(internal_path)==kTRUE)
{
std::cerr<<"cd "<<internal_path<<std::endl ;
file_ref_dir = internal_path ;
file_ref->cd() ;
}
else if (file_ref->cd(old_internal_path)==kTRUE)
{
std::cerr<<"cd "<<old_internal_path<<std::endl ;
file_ref_dir = old_internal_path ;
file_ref->cd() ;
}
else
{
std::cerr<<"Failed move to: "<<internal_path<<" and "<<old_internal_path<<std::endl ;
file_ref_dir = "" ;
}
}
else
{ std::cerr<<"Failed to open: "<<CMP_BLUE_FILE<<std::endl ; }
}
TString file_new_dir = internal_path ;
TFile * file_new = 0 ;
if ( CMP_RED_FILE != "" )
{
file_new = TFile::Open(CMP_RED_FILE) ;
if (file_new!=0)
{
std::cout<<"open "<<CMP_RED_FILE<<std::endl ;
if (file_new->cd(internal_path)==kTRUE)
{
std::cerr<<"cd "<<internal_path<<std::endl ;
file_new_dir = internal_path ;
file_new->cd() ;
}
else if (file_new->cd(old_internal_path)==kTRUE)
{
std::cerr<<"cd "<<old_internal_path<<std::endl ;
file_new_dir = old_internal_path ;
file_new->cd() ;
}
else
{
std::cerr<<"Failed move to: "<<internal_path<<" and "<<old_internal_path<<std::endl ;
file_new_dir = "" ;
}
}
else
{ std::cerr<<"Failed to open: "<<CMP_RED_FILE<<std::endl ; }
}
TCanvas * canvas ;
TH1 * histo_ref, * histo_new ;
TPaveStats * st_ref, * st_new ;
std::ofstream web_page(CMP_DIR+"/index.html") ;
web_page
<<"<!DOCTYPE HTML PUBLIC \"-//W3C//DTD HTML 3.2 Final//EN\">\n"
<<"<html>\n"
<<"<head>\n"
<<"<meta http-equiv=\"content-type\" content=\"text/html; charset=UTF-8\" />\n"
<<"<title>"<<CMP_TITLE<<"</title>\n"
<<"</head>\n"
<<"<h1><a href=\"../\"><img border=0 width=\"22\" height=\"22\" src=\"../../../../img/up.gif\" alt=\"Up\"/></a> "<<CMP_TITLE<<"</h1>\n" ;
if (file_ref==0)
{
web_page
<<"<p>In all plots below"
<<", there was no reference histograms to compare with"
<<", and the "<<CMP_RED_NAME<<" histograms are in red." ;
}
else
{
/* web_page
<<"<p>In all plots below"
<<", the "<<CMP_RED_NAME<<" histograms are in red"
<<", and the "<<CMP_BLUE_NAME<<" histograms are in blue." ;*/
web_page
<<"<p>In all plots below"
<<", the "<<CMP_RED_RELEASE<<" histograms are in red"
<<", and the "<<CMP_BLUE_RELEASE<<" histograms are in blue." ;
/* std::cout <<"<p>In all plots below "
<<", the "<<CMP_RED_RELEASE<<" histograms are in red"
<<", and the "<<CMP_BLUE_RELEASE<<" histograms are in blue." << std::endl ;*/
}
std::ifstream red_comment_file(CMP_RED_COMMENT) ;
std::string red_comment ;
std::getline(red_comment_file,red_comment) ;
red_comment_file.close() ;
std::ifstream blue_comment_file(CMP_BLUE_COMMENT) ;
std::string blue_comment ;
std::getline(blue_comment_file,blue_comment) ;
blue_comment_file.close() ;
web_page<<" "<<red_comment<<" "<<blue_comment
<<" Some more details"
<<": <a href=\"electronCompare.C\">script</a> used to make the plots"
<<", <a href=\""<<CMP_CONFIG<<"\">specification</a> of histograms"
<<", <a href=\"gifs/\">images</a> of histograms"
<<"." ;
web_page<<"</p>\n" ;
// canvas_name std::string => TString
TString canvas_name, histo_name, histo_full_path, gif_name, gif_path ;
TString short_histo_name ;
TString first_short_histo_name, first_histo_name ;
TString dl_short_histo_name, dl_histo_name ;
TString num_ref, denom_ref, num_full, denom_full ;
Int_t n_ele_charge ;
int scaled, err ;
int divide;
std::string cat, line, histo_path, num, denom ;
int eol ; // end of line
int eoc ; // enf of category
double rescale_factor = 0. ;
std::ifstream histo_file1(CMP_CONFIG) ;
web_page
<<"<br><table border=\"1\" cellpadding=\"5\" width=\"100%\">"
<<"<tr valign=\"top\"><td width=\"20%\">\n" ;
int cat_num = 0 ;
cat = "" ;
do
{
std::getline(histo_file1,cat) ;
} while (cat.find_first_not_of(" \t")==std::string::npos) ;
web_page<<"<b>"<<cat<<"</b><br><br>" ;
while (std::getline(histo_file1,line))
{
if (line.empty()) continue ;
std::size_t first = line.find_first_not_of(" \t") ;
if (first==std::string::npos) continue ;
if (line[first]=='#') continue ;
std::istringstream linestream(line) ;
divide = 0 ; num = denom = "" ;
linestream >> histo_path >> scaled >> err >> eol >> eoc >> divide >> num >> denom ;
histo_name = histo_path ;
Ssiz_t pos = histo_name.Last('/') ;
if (pos!=kNPOS) histo_name.Remove(0,pos+1) ;
short_histo_name = histo_name ;
short_histo_name.Remove(0,2) ;
if (short_histo_name.BeginsWith("ele_"))
{ short_histo_name.Remove(0,4) ; }
if (first_histo_name.IsNull())
{
first_short_histo_name = short_histo_name ;
first_histo_name = histo_name ;
dl_short_histo_name = short_histo_name ;
dl_histo_name = histo_name ;
}
else
{
TString common ;
TString first_specific ;
TString second_specific ;
CompareHistoNames(first_short_histo_name,short_histo_name,common,first_specific,second_specific) ;
if (!dl_histo_name.IsNull())
{
if (first_specific.IsNull())
{ web_page<<"<a href=\"#"<<first_histo_name<<"\">"<<first_short_histo_name<<"</a>" ; }
else
{ web_page<<common<<" | <a href=\"#"<<first_histo_name<<"\">"<<first_specific<<"</a>" ; }
dl_short_histo_name = dl_histo_name = "" ;
}
web_page<<" | <a href=\"#"<<histo_name<<"\">"<<second_specific<<"</a>" ;
}
if ((eol)||(eoc))
{
if (!dl_histo_name.IsNull())
{ web_page<<"<a href=\"#"<<dl_histo_name<<"\">"<<dl_short_histo_name<<"</a>" ; }
web_page<<"<br>\n" ;
first_histo_name = first_short_histo_name = "" ;
}
if (eoc)
{
cat_num++ ;
if ((cat_num%5)==0)
{ web_page<<"<br></td></tr>\n<tr valign=\"top\"><td width=\"20%\">" ; }
else
{ web_page<<"<br></td><td width=\"20%\">\n" ; }
cat = "" ;
do
{
std::getline(histo_file1,cat) ;
} while (cat.find_first_not_of(" \t")==std::string::npos) ;
web_page<<"<b>"<<cat<<"</b><br><br>" ;
}
}
web_page<<"<br></td></tr></table>\n" ;
histo_file1.close() ;
web_page<<"<br><br><table cellpadding=\"5\"><tr valign=\"top\"><td><a href=\""<<CMP_URL<<"/\"><img width=\"18\" height=\"18\" border=\"0\" align=\"middle\" src=\"../../../../img/up.gif\" alt=\"Top\"/></a></td><td>\n" ;
std::ifstream histo_file2(CMP_CONFIG) ;
n_ele_charge = 0 ;
cat = "" ;
do
{
std::getline(histo_file2,cat) ;
} while (cat.empty()) ;
while (std::getline(histo_file2,line))
{
if (line.empty()) continue ;
std::size_t first = line.find_first_not_of(" \t") ;
if (first==std::string::npos) continue ;
if (line[first]=='#') continue ;
std::istrstream linestream(line) ;
divide = 0 ; num = denom = "" ;
linestream >> histo_path >> scaled >> err >> eol >> eoc >> divide >> num >> denom ;
histo_name = histo_path.c_str() ;
histo_ref = 0 ;
histo_new = 0 ;
st_ref = 0 ;
st_new = 0 ;
Ssiz_t pos = histo_name.Last('/') ;
if (pos!=kNPOS) histo_name.Remove(0,pos+1) ;
short_histo_name = histo_name ;
short_histo_name.Remove(0,2) ;
gif_name = "gifs/" ; gif_name += histo_name ; gif_name += ".gif" ;
gif_path = CMP_DIR ; gif_path += "/" ; gif_path += gif_name ;
canvas_name = "c" ; canvas_name += histo_name ;
canvas = new TCanvas(canvas_name) ;
canvas->SetFillColor(10) ;
//std::cout << "canvas width " << canvas->GetWindowWidth() << std::endl; // 800 default
//std::cout << "canvas height " << canvas->GetWindowHeight() << std::endl; // 600 default
//std::cout << "canvas Y real " << canvas->GetYsizeReal() << std::endl;
//std::cout << "canvas Y user " << canvas->GetYsizeUser() << std::endl;
web_page<<"<a id=\""<<histo_name<<"\" name=\""<<short_histo_name<<"\"></a>" ;
// search histo_ref
if ( file_ref != 0 )
{
if (file_ref_dir.IsNull())
{ histo_full_path = histo_name ; }
else
{ histo_full_path = file_ref_dir ; histo_full_path += histo_path.c_str() ; }
histo_ref = (TH1 *)file_ref->Get(histo_full_path) ;
if (histo_ref!=0)
{
// renaming those histograms avoid very strange bugs because they
// have the same names as the ones loaded from the new file
histo_ref->SetName(TString(histo_ref->GetName())+"_ref") ;
}
else
{
web_page<<"No <b>"<<histo_path<<"</b> for "<<CMP_BLUE_NAME<<".<br>" ;
}
}
// search histo_new
histo_full_path = file_new_dir ; histo_full_path += histo_path.c_str() ;
histo_new = (TH1 *)file_new->Get(histo_full_path) ;
//std::cout << "size " << histo_new->GetSize() << std::endl ;
// special treatments
if ((scaled==1)&&(histo_new!=0)&&(histo_ref!=0)&&(histo_ref->GetEntries()!=0))
{
Int_t new_entries = histo_new->GetEntries() ;
if (new_entries==0) { new_entries = n_ele_charge ; }
if (new_entries==0)
{ std::cerr<<"DO NOT KNOW HOW TO RESCALE "<<histo_name<<std::endl ; }
else
{
// we want to reuse the rescale factor of the first histogram
// for all the subsequent histograms.
if (rescale_factor==0.) { rescale_factor = new_entries/histo_ref->GetEntries() ; }
histo_ref->Scale(rescale_factor) ;
}
}
if ((histo_new!=0)&&(histo_ref!=0)&&(histo_ref->GetMaximum()>histo_new->GetMaximum()))
{ histo_new->SetMaximum(histo_ref->GetMaximum()*1.1) ; }
if (histo_new==0)
{
web_page<<"No <b>"<<histo_path<<"</b> for "<<CMP_RED_NAME<<".<br>" ;
}
else
{
// catch n_ele_charge
if (histo_name=="h_ele_charge")
{ n_ele_charge = histo_new->GetEntries() ; }
// draw histo_new
TString newDrawOptions(err==1?"E1 P":"hist") ;
gErrorIgnoreLevel = kWarning ;
if (divide!=0)
{
num_full = file_new_dir ; num_full += num.c_str() ;
denom_full = file_new_dir ; denom_full += denom.c_str() ;
histo_new = DivideHistos(file_new,histo_new,num_full,denom_full) ;
}
histo_new->SetLineColor(kRed) ;
histo_new->SetMarkerColor(2) ;
histo_new->SetLineWidth(3) ;
RenderHisto(histo_new,canvas) ;
histo_new->Draw(newDrawOptions) ;
// std::cout << "SIZE : " << canvas->GetWw() << std::endl ; // 796 default
// std::cout << "SIZE : " << canvas->GetWh() << std::endl ; // 572 default
//canvas->Update() ;
//canvas->SetWindowSize(440, 600);
canvas->SetCanvasSize(960, 600);
canvas->Update() ;
st_new = (TPaveStats*)histo_new->FindObject("stats");
st_new->SetTextColor(kRed) ;
// draw histo_ref
if (histo_ref!=0)
{
if (divide!=0)
{
num_ref = num ;
denom_ref = denom ;
if (file_ref_dir.IsNull())
{
pos = num_ref.Last('/') ;
if (pos!=kNPOS) num_ref.Remove(0,pos+1) ;
pos = denom_ref.Last('/') ;
if (pos!=kNPOS) denom_ref.Remove(0,pos+1) ;
}
histo_ref = DivideHistos(file_ref,histo_ref,file_ref_dir+num_ref,file_ref_dir+denom_ref) ;
}
RenderHisto(histo_ref,canvas) ;
histo_ref->SetLineColor(kBlue) ;
histo_ref->SetLineWidth(3) ;
histo_ref->Draw("sames hist") ;
canvas->Update() ;
st_ref = (TPaveStats*)histo_ref->FindObject("stats");
st_ref->SetTextColor(kBlue) ;
Double_t y1 = st_ref->GetY1NDC() ;
Double_t y2 = st_ref->GetY2NDC() ;
st_ref->SetY1NDC(2*y1-y2) ;
st_ref->SetY2NDC(y1) ;
//Double_t x1 = st_ref->GetX1NDC() ;
//Double_t x2 = st_ref->GetX2NDC() ;
//std::cout << "position s x1 = " << x1 << std::endl ; // 0.78 par defaut
//std::cout << "position s x2 = " << x2 << std::endl ; // 0.98 par defaut
//std::cout << "position s y1 = " << y1 << std::endl ; // 0.755 ou 0.835 par defaut
//std::cout << "position s y2 = " << y2 << std::endl ; // 0.995 par defaut
}
// Redraws
newDrawOptions = "sames " ;
newDrawOptions += (err==1?"E1 P":"hist") ;
histo_new->Draw(newDrawOptions) ;
if (st_ref!=0) st_ref->Draw() ;
if (st_new!=0) st_new->Draw() ;
// eventual log scale
//if ( (log==1) && ( (histo_new->GetEntries()>0) || ( (histo_ref!=0) && (histo_ref->GetEntries()!=0) ) ) )
// { canvas->SetLogy(1) ; }
// ne pas oublier de decommenter les 4 lignes suivantes
std::cout<<histo_name
<<" has "<<histo_new->GetEffectiveEntries()<<" entries"
<<" of mean value "<<histo_new->GetMean()
<<std::endl ;
canvas->SaveAs(gif_path.Data()) ;
web_page<<"<a href=\""<<gif_name<<"\"><img border=\"0\" class=\"image\" width=\"440\" src=\""<<gif_name<<"\"></a><br>" ;
}
// else if ((file_ref!=0)&&(histo_ref!=0))
// {
// std::cout<<histo_path<<" NOT FOUND"<<std::endl ;
// web_page<<"<br>(no such histo)" ;
// canvas->SaveAs(gif_path.Data()) ;
// web_page<<"<a href=\""<<gif_name<<"\"><img border=\"0\" class=\"image\" width=\"440\" src=\""<<gif_name<<"\"></a><br>" ;
// }
if (eol)
{ web_page<<"</td></tr>\n<tr valign=\"top\"><td><a href=\""<<CMP_URL<<"/\"><img width=\"18\" height=\"18\" border=\"0\" align=\"middle\" src=\"../../../../img/up.gif\" alt=\"Top\"/></a></td><td>" ; }
else
{ web_page<<"</td><td>" ; }
if (eoc)
{
cat = "" ;
do
{
std::getline(histo_file2,cat) ;
} while (cat.empty()) ;
}
}
histo_file2.close() ;
web_page<<"</td></tr></table>\n" ;
web_page<<"\n</html>"<<std::endl ;
web_page.close() ;
}
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 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();
}
// 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){
// 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){
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();
}
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 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 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 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();
}
/////
// 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 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 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 setStyle(){
//Style//
TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");
// For the canvas:
tdrStyle->SetCanvasBorderMode(0);
tdrStyle->SetCanvasColor(kWhite);
tdrStyle->SetCanvasDefH(600); //Height of canvas
tdrStyle->SetCanvasDefW(600); //Width of canvas
tdrStyle->SetCanvasDefX(0); //POsition on screen
tdrStyle->SetCanvasDefY(0);
// For the Pad:
tdrStyle->SetPadBorderMode(0);
tdrStyle->SetPadColor(kWhite);
tdrStyle->SetPadGridX(false);
tdrStyle->SetPadGridY(false);
tdrStyle->SetGridColor(0);
tdrStyle->SetGridStyle(3);
tdrStyle->SetGridWidth(1);
// For the frame:
tdrStyle->SetFrameBorderMode(0);
tdrStyle->SetFrameBorderSize(1);
tdrStyle->SetFrameFillColor(0);
tdrStyle->SetFrameFillStyle(0);
tdrStyle->SetFrameLineColor(1);
tdrStyle->SetFrameLineStyle(1);
tdrStyle->SetFrameLineWidth(1);
tdrStyle->SetPadTopMargin(0.05);
tdrStyle->SetPadBottomMargin(0.13);
tdrStyle->SetPadLeftMargin(0.16);
tdrStyle->SetPadRightMargin(0.02);
// For the Global title:
tdrStyle->SetOptTitle(0);
tdrStyle->SetTitleFont(42);
tdrStyle->SetTitleColor(1);
tdrStyle->SetTitleTextColor(1);
tdrStyle->SetTitleFillColor(10);
tdrStyle->SetTitleFontSize(0.05);
// For the axis titles:
tdrStyle->SetTitleColor(1, "XYZ");
tdrStyle->SetTitleFont(42, "XYZ");
tdrStyle->SetTitleSize(0.06, "XYZ");
tdrStyle->SetTitleXOffset(0.9);
tdrStyle->SetTitleYOffset(0.9);
// For the axis labels:
tdrStyle->SetLabelColor(1, "XYZ");
tdrStyle->SetLabelFont(42, "XYZ");
tdrStyle->SetLabelOffset(0.007, "XYZ");
tdrStyle->SetLabelSize(0.05, "XYZ");
// For the axis:
tdrStyle->SetAxisColor(1, "XYZ");
//tdrStyle->SetStripDecimals(kTRUE);
tdrStyle->SetTickLength(0.03, "XYZ");
tdrStyle->SetNdivisions(510, "XYZ");
tdrStyle->SetPadTickX(1); // To get tick marks on the opposite side of the frame
tdrStyle->SetPadTickY(1);
tdrStyle->SetPaperSize(20.,20.);
tdrStyle->SetHatchesLineWidth(5);
tdrStyle->SetHatchesSpacing(0.05);
tdrStyle->SetOptStat(0);
tdrStyle->cd();
}
void 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");
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();
}
int main()
{
gROOT->Reset();
TStyle * Sty = new TStyle("MyStyle","MyStyle");
Sty->SetOptTitle(0);
Sty->SetOptStat(0);
Sty->SetPalette(8,0);
Sty->SetCanvasColor(10);
Sty->SetCanvasBorderMode(0);
Sty->SetFrameLineWidth(3);
Sty->SetFrameFillColor(10);
Sty->SetPadColor(10);
Sty->SetPadTickX(1);
Sty->SetPadTickY(1);
Sty->SetPadBottomMargin(.2);
Sty->SetPadLeftMargin(.2);
Sty->SetHistLineWidth(3);
Sty->SetHistLineColor(kRed);
Sty->SetFuncWidth(3);
Sty->SetFuncColor(kGreen);
Sty->SetLineWidth(3);
Sty->SetLabelSize(0.06,"xyz");
Sty->SetLabelOffset(0.02,"y");
Sty->SetLabelOffset(0.02,"x");
Sty->SetLabelColor(kBlack,"xyz");
Sty->SetTitleSize(0.06,"xyz");
Sty->SetTitleOffset(1.4,"y");
Sty->SetTitleOffset(1.3,"x");
Sty->SetTitleFillColor(10);
Sty->SetTitleTextColor(kBlack);
Sty->SetTickLength(.05,"xz");
Sty->SetTickLength(.025,"y");
Sty->SetNdivisions(5,"xyz");
Sty->SetEndErrorSize(0);
gROOT->SetStyle("MyStyle");
gROOT->ForceStyle();
line FermiN;
line N0s12;
line N0p32;
line N0p12;
line N0d52;
line N1s12;
line N0d32;
line N0f72;
line N0f52;
line N1p12;
line N1p32;
line N0g92;
line N1d52;
line NN0d52;
line NN0d32;
line NN1s12;
line NN0f72;
line NN1p12;
line NN1p32;
line NN0f52;
line FermiP;
line P0s12;
line P0p32;
line P0p12;
line P0d52;
line P1s12;
line P0d32;
line P0f72;
line P0f52;
line P1p12;
line P1p32;
line P0g92;
line P1d52;
line PP0d52;
line PP0d32;
line PP1s12;
line PP0f72;
line PP1p12;
line PP1p32;
line PP0f52;
float z,a,three,Efermi;
//TFile *f = new TFile("alev.root","RECREATE");
char stuff[100];
ostringstream outstring;
string name;
ifstream file;
int Nneut = 0;
float j,energy,RMS,occup,spect,Delta;
float error,RMSerror,spectError,DeltaError;
int l,N,colour,Efit,Rfit,Dfit,sfit;
// neutrons
for (int A=36;A<74;A+=2)
{
outstring.str("");
outstring << "nca" << A << ".inp";
name = outstring.str();
file.open(name.c_str());
if (file.fail())
{
cout << name << " not opened " << endl;
file.clear();
continue;
}
Nneut++;
file >> three >> z >> a >> Efermi;
file.close();
file.clear();
if (Efermi == 0.)
{
outstring.str("");
outstring << "nca" << A << ".fermi";
name = outstring.str();
file.open(name.c_str());
if (file.fail())
{
cout << name << " not opened " << endl;
file.clear();
continue;
}
file >> Efermi;
file.close();
file.clear();
}
cout << Efermi << endl;
FermiN.add((float)A,Efermi);
outstring.str("");
outstring << "nca" << A << ".level";
name = outstring.str();
file.open(name.c_str());
if (file.fail())
{
cout << name << " not opened " << endl;
file.clear();
continue;
}
file.getline(stuff,100);
for (;;)
{
file >> N >> j >> l >> energy >> RMS >> occup
>> spect >> Delta>> colour;
if (file.eof())break;
if (l == 0 && j == 0.5 && N ==0 ) N0s12.add((float)A,energy);
if (l == 0 && j == 0.5 && N ==1 ) N1s12.add((float)A,energy);
if (l == 1 && j == 0.5 && N ==0 ) N0p12.add((float)A,energy);
if (l == 1 && j == 1.5 && N ==0 ) N0p32.add((float)A,energy);
if (l == 2 && j == 2.5 && N ==0 ) N0d52.add((float)A,energy);
if (l == 2 && j == 1.5 && N ==0 ) N0d32.add((float)A,energy);
if (l == 3 && j == 3.5 && N ==0 ) N0f72.add((float)A,energy);
if (l == 3 && j == 2.5 && N ==0 ) N0f52.add((float)A,energy);
if (l == 1 && j == 0.5 && N ==1 ) N1p12.add((float)A,energy);
if (l == 1 && j == 1.5 && N ==1 ) N1p32.add((float)A,energy);
if (l == 4 && j == 4.5 && N ==0 ) N0g92.add((float)A,energy);
if (l == 2 && j == 2.5 && N ==1 ) N1d52.add((float)A,energy);
}
file.close();
file.clear();
outstring.str("");
outstring << "nca" << A << ".lev";
name = outstring.str();
file.open(name.c_str());
if (file.fail())
{
cout << name << " not opened " << endl;
file.clear();
continue;
}
file.getline(stuff,100);
file.getline(stuff,100);
for (;;)
{
file >> energy >> error >> N >> j >> l >> colour >>Efit>>
RMS >> RMSerror >> Rfit >> Delta>>DeltaError >> Dfit >>
spect>>spectError>>sfit;
if (file.eof())break;
//if (l == 0 && j == 0.5 && N ==0 ) NN0s12.add((float)A,energy);
if (l == 0 && j == 0.5 && N ==1 ) NN1s12.add((float)A,energy);
//if (l == 1 && j == 0.5 && N ==0 ) NN0p12.add((float)A,energy);
//if (l == 1 && j == 1.5 && N ==0 ) NN0p32.add((float)A,energy);
if (l == 2 && j == 2.5 && N ==0 ) NN0d52.add((float)A,energy);
if (l == 2 && j == 1.5 && N ==0 ) NN0d32.add((float)A,energy);
if (l == 3 && j == 3.5 && N ==0 ) NN0f72.add((float)A,energy);
if (l == 3 && j == 2.5 && N ==0 ) NN0f52.add((float)A,energy);
//if (l == 3 && j == 2.5 && N ==0 ) NN0f52.add((float)A,energy);
if (l == 1 && j == 0.5 && N ==1 ) NN1p12.add((float)A,energy);
if (l == 1 && j == 1.5 && N ==1 ) NN1p32.add((float)A,energy);
//if (l == 4 && j == 4.5 && N ==0 ) NN0g92.add((float)A,energy);
}
file.close();
file.clear();
}
