void cmsPrel(double intLumi=-1, bool wide = false, long int sqrtS = 7) {
TLatex *latex = new TLatex();
latex->SetNDC();
latex->SetTextSize(0.045);
double wideLeftOffset = 0.19; //for JMEStyle, was 0.15 by default;
latex->SetTextAlign(31); // align right
latex->DrawLatex(wide ? 0.98 : 0.95, 0.96, ((TString)"#sqrt{s} = "+ sqrtS)+ " TeV");
if (intLumi > 0.) {
latex->SetTextAlign(11); // align left
//latex->DrawLatex(wideLeftOffset,0.96,Form("CMS preliminary, %.2g nb^{-1}",intLumi));
latex->DrawLatex(wide ? 0.06 : wideLeftOffset, 0.96,
// Form(" ",intLumi));
// Form("CMS, %.2g pb^{-1}",intLumi));
//was Form("CMS preliminary, %.2g pb^{-1}",intLumi));
Form("CMS preliminary, %.2g fb^{-1}",intLumi));
}
else if (intLumi==0) { // simulation
latex->SetTextAlign(11); // align left
latex->DrawLatex(wide ? 0.06 : wideLeftOffset, 0.96, "CMS simulation (Su12 QCD)");
}
else {
latex->SetTextAlign(11); // align left
// latex->DrawLatex(wideLeftOffset,0.96,"CMS preliminary");
latex->DrawLatex(wide ? 0.06 : wideLeftOffset,0.96,"CMS work in progress");
}
} // cmsPrel
TCanvas *subsupscript()
{
TCanvas *S = new TCanvas("script","Subscripts and Superscripts",400,250);
TLatex Tl;
Tl.SetTextSize(0.09);
Tl.SetTextAlign(12);
TLatex Tt;
Tt.SetTextSize(0.09);
Tt.SetTextFont(82);
Tt.SetTextAlign(12);
// Draw First Column
float y, x1, x2;
float step = 0.3;
y = 0.80; x1 = 0.02; x2 = x1+0.3;
Tl.DrawLatex(x2, y, "x^{2y}") ; Tt.DrawText(x1, y, "x^{2y}");
y -= step ; Tl.DrawLatex(x2, y, "x_{2y}") ; Tt.DrawText(x1, y, "x_{2y}");
y -= step ; Tl.DrawLatex(x2, y, "x^{y^{2}}") ; Tt.DrawText(x1, y, "x^{y^{2}}");
// Draw Second Column
y = 0.80; x1 = 0.52; x2 = x1+0.3;
Tl.DrawLatex(x2, y, "x^{y_{1}}") ; Tt.DrawText(x1, y, "x^{y_{1}}");
y -= step ; Tl.DrawLatex(x2, y, "x^{y}_{1}") ; Tt.DrawText(x1, y, "x^{y}_{1}");
y -= step ; Tl.DrawLatex(x2, y, "x_{1}^{y}") ; Tt.DrawText(x1, y, "x_{1}^{y}");
return S;
}
void loglikdistrib(Int_t ntrials = 10000, Bool_t print = kFALSE)
{
// compute distribution of log likelihood value
TH1D * hmc = gStack[gPadNr][gOrder[gPadNr][0]];
TH1D * hdata = gStack[gPadNr][gMaxProcess-1];
Int_t nbins = hmc->GetNbinsX();
Double_t loglik = loglikelihood(hmc, hdata, 1, nbins);
TH1D * htest = new TH1D(*hdata);
TH1D * lldistrib = new TH1D("lldistrib", "log(Likelihood) distribution",
1000, loglik-200, loglik+200);
setopt(lldistrib);
for (Int_t n = 0; n < ntrials; n++) {
// generate poisson around theorie
for (Int_t i = 1; i <= nbins; i++) {
htest->SetBinContent(i, gRandom->Poisson(hmc->GetBinContent(i)));
}
lldistrib->Fill(loglikelihood(hmc, htest, 1, nbins));
}
TCanvas * llcanvas = new TCanvas("llcanvas", "Log(Likelihood) distribution",
40, 40, 800, 600);
setopt(llcanvas);
lldistrib->SetFillColor(kYellow);
lldistrib->Draw();
lldistrib->GetYaxis()->SetTitle("Anzahl Ereignisse");
lldistrib->GetXaxis()->SetTitle("-ln L");
// autozoom
Int_t lowbin = 1;
while (lldistrib->GetBinContent(lowbin) == 0)
lowbin++;
Int_t highbin = lldistrib->GetNbinsX();
while (lldistrib->GetBinContent(highbin) == 0)
highbin--;
lldistrib->SetAxisRange(lldistrib->GetBinLowEdge(lowbin),
lldistrib->GetBinLowEdge(highbin));
TH1D * hworse = (TH1D *) lldistrib->Clone();
for (Int_t nbin = 1; nbin < 501; nbin++) {
hworse->SetBinContent(nbin, 0);
}
hworse->SetFillColor(95);
hworse->Draw("same");
Double_t pvalue = lldistrib->Integral(501,1000) / lldistrib->Integral();
TLatex * tex = new TLatex(0.18, 0.96, Form("-ln L_{obs} = %5.2f", loglik));
tex->SetNDC();
tex->SetTextAlign(13);
tex->Draw();
tex = new TLatex(0.18, 0.86, Form("CL_{obs} = %.3f", pvalue));
tex->SetNDC();
tex->SetTextAlign(13);
tex->Draw();
TLine * l = new TLine(loglik, 0, loglik, lldistrib->GetMaximum());
l->SetLineWidth(3);
l->SetLineColor(kBlue);
l->Draw();
llcanvas->Modified();
llcanvas->Update();
if (print)
llcanvas->Print("lldistrib.pdf");
cd(gPadNr+1);
}
void xyplot()
{
TCanvas *c = new TCanvas("c","XY plot",200,10,700,500);
// Remove the frame
c->SetFillColor(kWhite);
c->SetFrameLineColor(kWhite);
c->SetFrameBorderMode(0);
// Define and draw a curve the frame
const Int_t n = 4;
Double_t x[n] = {-1, -3, -9, 3};
Double_t y[n] = {-1000, 900, 300, 300};
gr = new TGraph(n,x,y);
gr->SetTitle("XY plot");
gr->SetMinimum(-1080);
gr->SetMaximum(1080);
gr->SetLineColor(kRed);
gr->Draw("AC*");
// Remove the frame's axis
gr->GetHistogram()->GetYaxis()->SetTickLength(0);
gr->GetHistogram()->GetXaxis()->SetTickLength(0);
gr->GetHistogram()->GetYaxis()->SetLabelSize(0);
gr->GetHistogram()->GetXaxis()->SetLabelSize(0);
gr->GetHistogram()->GetXaxis()->SetAxisColor(0);
gr->GetHistogram()->GetYaxis()->SetAxisColor(0);
gPad->Update();
// Draw orthogonal axis system centered at (0,0).
// Draw the Y axis. Note the 4th label is erased with SetLabelAttributes
TGaxis *yaxis = new TGaxis(0, gPad->GetUymin(),
0, gPad->GetUymax(),
gPad->GetUymin(),gPad->GetUymax(),6,"+LN");
yaxis->ChangeLabel(4,-1,0.);
yaxis->Draw();
// Draw the Y-axis title.
TLatex *ytitle = new TLatex(-0.5,gPad->GetUymax(),"Y axis");
ytitle->Draw();
ytitle->SetTextSize(0.03);
ytitle->SetTextAngle(90.);
ytitle->SetTextAlign(31);
// Draw the X axis
TGaxis *xaxis = new TGaxis(gPad->GetUxmin(), 0,
gPad->GetUxmax(), 0,
gPad->GetUxmin(),gPad->GetUxmax(),510,"+L");
xaxis->Draw();
// Draw the X axis title.
TLatex *xtitle = new TLatex(gPad->GetUxmax(),-200.,"X axis");
xtitle->Draw();
xtitle->SetTextAlign(31);
xtitle->SetTextSize(0.03);
}
void TPTiming ()
{
TAxis * ax = TPMatchEmul2D->GetZaxis() ;
ax->SetRangeUser(-1,6) ;
TCanvas* canv = new TCanvas("canv", "canv") ;
canv->SetLogz(0) ;
gStyle->SetOptStat(10) ;
int color[10]= {1,10,3,4,5,6,7,8,9,2} ;
gStyle->SetPalette(7, color) ;
TPMatchEmul2D->GetXaxis()->SetTitle("Phi index");
TPMatchEmul2D->GetYaxis()->SetTitle("Eta index");
TPMatchEmul2D->Draw("colz") ;
TH2I * label = new TH2I("label", "",72, 1, 73, 38, -19, 19) ;
label->SetMarkerSize(0.6);
label->SetBit(kCanDelete);
for (int x=3 ; x<73 ; x+=4) {
for (int y=21; y<=37; y++) {
int towernb = 4*(y-21)+1 ;
label->SetBinContent(x-1, y, towernb) ; //EB+
label->SetBinContent(x, 40-y, towernb) ; //EB-
}
}
label->Draw("same text") ;
TLatex txt;
txt.SetTextSize(0.02);
TLine line;
line.SetLineColor(1) ;
line.SetLineStyle(1) ;
line.SetLineWidth(1) ;
TAxis* xAxis = TPMatchEmul2D->GetXaxis();
TAxis* yAxis = TPMatchEmul2D->GetYaxis();
// draw SM borders and numbers
float sm ;
for (int i=0; i<36 ; i++ ) {
if (i<18) {
sm = 4*i+3 ;
line.DrawLine(sm, 1, sm, 18) ;
txt.SetTextAlign(32);
txt.DrawText(sm-1+0.3, -17.7, Form("-%d",i+1));
}
else {
sm = 4*(i-18)+3 ;
line.DrawLine(sm, 0, sm, -17) ;
txt.SetTextAlign(12);
txt.DrawText(sm-2+0.3, 18.5, Form("+%d",i-17));
}
}
line.DrawLine(1, 0, 73, 0) ;
line.DrawLine(1, -17, 73, -17) ;
line.DrawLine(1, 1, 73, 1) ;
line.DrawLine(1, 18, 73, 18) ;
}
////CMS Preliminary label and luminosity + any aditional inline title
void cmsPrelim(const char * label)
{
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.04);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.90,0.96,"CMS 2012, #sqrt{s} = 8 TeV");
latex.SetTextAlign(11); // align left
latex.DrawLatex(0.18,0.96, label );
}
void topTitle(const char *title)
{
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.06);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.90,0.92, title);
latex.SetTextAlign(11); // align left
latex.DrawLatex(0.18,0.92,"Preliminary");
}
////CMS Preliminary label and luminosity
void cmsPrelim()
{
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.04);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.90,0.96,"#sqrt{s} = 8 TeV");
latex.SetTextAlign(11); // align left
latex.DrawLatex(0.18,0.96,"CMS Preliminary 2012");
}
void spectra_mw(const char *prefix = "pt/PhiNsigma_KTPCnsig30_STD2010_00_DEFAULT_00")
{
TGraphErrors *g = ReadMWGraph(Form("%s/mass", prefix), 1);
g->Draw("ACP");
TLatex latex;
latex.SetTextFont(42);
latex.SetTextAlign(12);
latex.SetTextSize(0.03);
latex.SetTextColor(kBlue+1);
latex.DrawLatex(1.2, 1.01877, "MC generated");
latex.SetTextColor(kGreen+1);
latex.DrawLatex(1.2, 1.01877-0.00015, "MC reconstructed");
latex.SetTextColor(kRed+1);
latex.DrawLatex(1.2, 1.01877-2*0.000155, "Real data (uncorrected)");
latex.SetTextColor(kBlack);
latex.DrawLatex(1.2, 1.01877-3*0.00016, "Real data (corrected)");
gPad->Print(Form("%s.pdf", g->GetName()));
gPad->Clear();
gPad->SetLogy(0);
TGraphErrors *g = ReadMWGraph(Form("%s/mass", prefix), 3);
g->Draw("AP");
latex.SetTextAlign(22);
latex.SetTextSize(0.04);
latex.SetTextColor(kBlack);
latex.DrawLatex(2.0, 1.0207, "pp @ #sqrt{s} = 2.76 TeV");
latex.DrawLatex(2.0, 1.0205, "TPC 3.0 #sigma");
latex.DrawLatex(2.0, 1.0202, "#phi #rightarrow K^{ +} K^{ -}");
latex.SetTextSize(0.035);
latex.DrawLatex(1.7, 1.0182, "uncertainties: stat. (bars) & syst. (boxes)");
gPad->Print(Form("%s.pdf", g->GetName()));
gPad->Clear();
TGraphErrors *g = ReadMWGraph(Form("%s/width", prefix), 2);
g->Draw("AP");
TLatex latex;
latex.SetTextAlign(22);
latex.SetTextSize(0.04);
latex.SetTextColor(kBlack);
latex.DrawLatex(1.7, 0.0090, "pp @ #sqrt{s} = 2.76 TeV");
latex.DrawLatex(1.7, 0.0085, "TPC 3.0 #sigma");
latex.DrawLatex(1.7, 0.0077, "#phi #rightarrow K^{ +} K^{ -}");
latex.SetTextSize(0.035);
latex.DrawLatex(1.7, 0.0007, "uncertainties: stat. (bars) & syst. (boxes)");
gPad->Print(Form("%s.pdf", g->GetName()));
// gPad->Clear();
// gPad->SetLogy();
// TGraphErrors *g = ReadPtGraphSys("pt/out.txt");
// if (!g) return;
// g->Draw("AP");
// gPad->Print(Form("%s.pdf", g->GetName()));
}
// ----------------------------------------------------------------------
void styles::setTitles(TH1 *h, TString xTitle, TString yTitle, TString Left, TString Right) {
if (0==h) {
cout << " Histogram not defined" << endl;
} else {
h->SetXTitle(xTitle); h->SetYTitle(yTitle);
TLatex label; label.SetNDC(kTRUE);
label.SetTextSize(0.06);
label.SetTextAlign(11);
label.DrawLatex(PadLeftMargin,1-PadTopMargin+0.02,Left);
label.SetTextAlign(31);
label.DrawLatex(1-PadRightMargin,1-PadTopMargin+0.02,Right);
}
}
void cmsPrelim( double intLumi ){ TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.04);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.9,0.96,Form("%.1f fb^{-1}, #sqrt{s} = 8 TeV",intLumi/1000));
// if (intLumi > 0.) {
// latex.SetTextAlign(31); // align right
// latex.DrawLatex(0.9,0.9,Form("#int #font[12]{L} dt = %.1f fb^{-1}",intLumi));
// }
latex.SetTextAlign(11); // align left
latex.DrawLatex(0.25,0.96,"CMS preliminary");
}
////CMS Preliminary label and lumu - Thanks to Kalanand
void Histograms::cmsPrelim()
{
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.03);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.90,0.96,"#sqrt{s} = 8 TeV, 2.0E34 cm^{-2}s^{-1}");
latex.SetTextAlign(11); // align left
latex.DrawLatex(0.18,0.96,"CMS preliminary 2012");
}
//____________________________________________________________________
void DrawTitlePage(const TCollection* c)
{
fBody->cd();
Double_t y = .7;
TLatex* ltx = new TLatex(.5, y, "AOD #rightarrow P(#it{N}_{ch} )");
ltx->SetTextSize(0.07);
ltx->SetTextFont(62);
ltx->SetTextAlign(22);
ltx->SetNDC();
ltx->Draw();
Bool_t mc = GetObject(c, "mcVertex", false) != 0;
if (mc) {
y -= 0.05;
TLatex* sub = new TLatex(.5, y, "(Simulation input)");
sub->SetTextSize(0.04);
sub->SetTextFont(42);
sub->SetTextAlign(22);
sub->SetNDC();
sub->Draw();
}
Double_t save = fParName->GetTextSize();
fParName->SetTextSize(0.03);
fParVal->SetTextSize(0.03);
y = .6;
UShort_t sys;
UShort_t sNN;
ULong_t trig;
Double_t minIpZ;
Double_t maxIpZ;
GetParameter(c, "sys", sys);
GetParameter(c, "sNN", sNN);
GetParameter(c, "trigger", trig);
GetParameter(c, "minIpZ", minIpZ);
GetParameter(c, "maxIpZ", maxIpZ);
TString tT; TriggerString(trig, tT); DrawParameter(y, "Trigger", tT);
TString tS; SysString(sys, tS); DrawParameter(y, "System", tS);
TString tE; SNNString(sNN, tE); DrawParameter(y, "#sqrt{s_{NN}}", tE);
DrawParameter(y, "IP_{z} range", Form("%+5.2fcm - %+5.2fcm",
minIpZ, maxIpZ));
PrintCanvas("Title page");
fParName->SetTextSize(save);
fParVal->SetTextSize(save);
}
void cmspre()
{
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.05);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.85,0.93,"#sqrt{s} = 8 TeV");
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.65,0.93,Form("#scale[0.5]{#lower[-0.15]{#it{#int}}}#it{L} dt = %0.1f#kern[0.2]{fb}^{-1}", 19.3));
latex.SetTextAlign(11); // align left
// latex.DrawLatex(0.15,0.93,"CMS, #sqrt{s} = 7 TeV");//preliminary 2011");
latex.DrawLatex(0.15,0.93,"CMS Preliminary");
}
void drawCanvas()
{
cPhase = new TCanvas("cPhase","cPhase",0,0,500,500);
cPhase->SetLeftMargin(0.2);
cPhase->SetBottomMargin(0.2);
cPhase->Range(-0.2,-0.2,1.1,1.1);
TLatex *latex = new TLatex(1.0,-0.1,"#mu_{B}");
latex->SetTextAlign(31);
latex->Draw();
latex = new TLatex(-0.1,1.0,"T");
latex->SetTextAlign(31);
latex->SetTextAngle(90);
latex->Draw();
}
//----------------------------------------------------------------------
/// show the jets contained in filename (as produced by
/// ClusterSequence::print_jets_for_root()), with an optional label
TCanvas * showjets (const char * filename, const char * label = 0) {
// display the various 2-d drawing options
gROOT->Reset();
// set up canvas
TCanvas * lego = new TCanvas("lego","lego options",400,50,800,600);
lego->SetTheta(30.0);
lego->SetPhi(20.0);
// orientation used for plots in subtraction paper
//lego->SetTheta(62.15);
//lego->SetPhi(9.15);
////vector<double> col
set_default_colours(lego);
TPaveLabel pl;
JetHist * jets = new JetHist(filename);
jets->stack.Draw("lego1"); // cyl does not work with 5.16
if (label != 0) {
Float_t x1=0.63, y1=0.875, x2=0.95, y2=0.925;
pl.DrawPaveLabel(x1,y1,x2,y2,label,"brNDC");
} else if (jets->comment() != "") {
Float_t x1=0.15, y1=0.875, x2=0.95, y2=0.925;
pl.DrawPaveLabel(x1,y1,x2,y2,jets->comment().c_str(),"brNDC");
}
// normal histogram labels not working, so draw them by hand
TLatex l;
l.SetTextAlign(22);
l.SetTextSize(0.05);
//l.DrawLatex(0.0,0.85,"anti-k_{t}, R=1");
l.SetTextSize(0.04);
l.DrawLatex(0.20,-0.98,"y");
l.SetTextAlign(32);
l.DrawLatex(-0.7,0.8,"p_{t} [GeV]");
l.DrawLatex(-0.6,-0.78,"#phi");
// do not delete jets -- otherwise you lose everything!;
return lego;
///
lego->Update();
}
void myBoxText(Double_t x, Double_t y,Double_t boxsize,Int_t mcolor,char *text)
{
Double_t tsize=0.06;
TLatex l; l.SetTextAlign(12); //l.SetTextSize(tsize);
l.SetNDC();
l.DrawLatex(x,y,text);
Double_t y1=y-0.25*tsize;
Double_t y2=y+0.25*tsize;
Double_t x2=x-0.3*tsize;
Double_t x1=x2-boxsize;
printf("x1= %f x2= %f y1= %f y2= %f \n",x1,x2,y1,y2);
TPave *mbox= new TPave(x1,y1,x2,y2,0,"NDC");
mbox->SetFillColor(mcolor);
mbox->SetFillStyle(1001);
mbox->Draw();
TLine mline;
mline.SetLineWidth(4);
mline.SetLineColor(1);
mline.SetLineStyle(1);
Double_t y_new=(y1+y2)/2.;
mline.DrawLineNDC(x1,y_new,x2,y_new);
}
void Draw_vn( TH1D * h_vn_pt, TString MBorDtrig = "MBtrig", TString EPorSP = "EP", TString vnname = "v2", TString Ytitle = "v_{2}", int cent_low = 0, int cent_high = 100, double ptlow = 1.0, double pthigh = 35.0, TString fitoption = "poly3bkg")
{
TCanvas * cfg_vn = new TCanvas(Form("cfg_%s_%s", MBorDtrig.Data(), vnname.Data()));
h_vn_pt->GetXaxis()->SetRangeUser(ptlow+0.05, pthigh-0.05);
h_vn_pt->GetYaxis()->SetRangeUser(-0.2, 0.4);
h_vn_pt->GetXaxis()->SetTitle("D^{0} p_{T} (GeV/c)");
h_vn_pt->GetYaxis()->SetTitle(Ytitle);
h_vn_pt->GetYaxis()->SetTitleSize(0.05);
h_vn_pt->GetXaxis()->SetTitleSize(0.05);
if( EPorSP == "EP" )
{
h_vn_pt->SetMarkerColor(6.0);
h_vn_pt->SetLineColor(6.0);
h_vn_pt->SetMarkerStyle(22);
h_vn_pt->SetMarkerSize(1.0);
}
else if( EPorSP == "SP" )
{
h_vn_pt->SetMarkerColor(4.0);
h_vn_pt->SetLineColor(4.0);
h_vn_pt->SetMarkerStyle(21);
h_vn_pt->SetMarkerSize(1.3);
}
h_vn_pt->Draw();
TLatex Tl;
Tl.SetNDC();
Tl.SetTextAlign(12);
Tl.SetTextSize(0.05);
Tl.SetTextFont(42);
Tl.DrawLatex(0.125,0.965, "#font[61]{CMS} #scale[0.8]{Preliminary}");
Tl.DrawLatex(0.57,0.965, "#scale[0.8]{PbPb #sqrt{s_{NN}} = 5.02 TeV}");
TLatex* tex;
tex = new TLatex(0.20,0.83,"|y| < 1.0");
tex->SetNDC();
tex->SetTextFont(42);
tex->SetTextSize(0.04);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.20,0.78,Form("Cent. %d-%d%%", cent_low, cent_high));
tex->SetNDC();
tex->SetTextFont(42);
tex->SetTextSize(0.04);
tex->SetLineWidth(2);
tex->Draw();
TF1 * fun = new TF1("fun", "0.0", 0, 100);
fun->SetLineColor(1.0);
fun->SetLineStyle(3);
fun->SetLineWidth(1);
fun->Draw("same");
cfg_vn->SaveAs(Form("Plots_vn/cfg_vnvsmass_%s_%s_cent%dto%d_%s.pdf", MBorDtrig.Data(), vnname.Data(), cent_low, cent_high, fitoption.Data()));
}
void PrintItLog(TPad *pad, TString title)//char *title)
{
TLatex *latex = new TLatex();
latex->SetTextFont( 42);
latex->SetTextSize(0.05);
// Get the most recent changes
pad->Update();
double xmin = pad->GetUxmin();
double xmax = pad->GetUxmax();
//double ymin = pad->GetUymin();
double ymax = pad->GetUymax();
//cout << "xmin = " << xmin << endl;
//cout << "xmax = " << xmax << endl;
//cout << "ymin = " << ymin << endl;
//cout << "ymax = " << ymax << endl;
//cout << "TMath::Power(10,ymax) = " << TMath::Power(10,ymax) << endl;
double xpos = xmin + 0.50*(xmax - xmin);
double ypos = 1.5*TMath::Power(10,ymax);// + TMath::*(ymax - ymin);
latex->SetTextAlign(22);
latex->DrawLatex(xpos,ypos,title);
}
////CMS Preliminary label and lumu -- upper left corner
void cmsPrelim2()
{
const float LUMINOSITY = IntLUMI;
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.04);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.90,0.96,"#sqrt{s} = 13 TeV");
if (LUMINOSITY > 0.) {
latex.SetTextAlign(11); // align left
latex.DrawLatex(0.21,0.85,Form("#int #font[12]{L} dt = %.1f fb^{-1}", LUMINOSITY));
}
latex.SetTextAlign(11); // align left
latex.DrawLatex(0.18,0.96,"CMS preliminary");
}
void DrawTLatex(Double_t x,Double_t y, Double_t tsize, Short_t align, TString text){
TLatex* tl = new TLatex(x, y, text);
tl->SetNDC (true);
tl->SetTextAlign( align);
tl->SetTextFont ( 42);
tl->SetTextSize ( tsize);
tl->Draw("same");
}
void drawtext(double x, double y, int f, char *s)
{
TLatex *t = new TLatex(x,y,Form("#font[41]{%d :} %s",f,s));
t->SetTextFont(f);
t->SetTextAlign(12);
t->SetTextSize(0.048);
t->Draw();
}
void cmspre(double intlumifbinv)
{
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.04);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.85,0.93,"#sqrt{s} = 8 TeV");
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.65,0.93,Form("#int #font[12]{L} dt = %.1f fb^{-1}", (float)intlumifbinv));
// latex.DrawLatex(0.65,0.93,Form("2012B", (float)intlumifbinv));
latex.SetTextAlign(11); // align left
// latex.DrawLatex(0.15,0.93,"CMS, #sqrt{s} = 7 TeV");//preliminary 2011");
latex.DrawLatex(0.15,0.96,"CMS preliminary");
}
void comparecfback(TString file1, TString file2, TString histname1, TString histname2, TString histname3,TString sys) {
TFile* ifile1_ = new TFile(Form("%s",file1.Data()), "read");
TFile* ifile2_ = new TFile(Form("%s",file2.Data()), "read");
TH1D* hist1;
TH1D* hist2;
TH1D* hist3;
// if ( (strcmp(sys,"V0LL") == 0) || (strcmp(sys,"V0ALAL") == 0) || (strcmp(sys,"V0LAL") == 0) ) {
// hist1 = (TH1D*)ifile1_->Get(Form("Numcqinv%stpcM%d%s",sys,mult,kT));
// hist2 = (TH1D*)ifile2_->Get(Form("Numcqinv%stpcM%d%s",sys,mult,kT));
// }
// else {
hist1 = (TH1D*)ifile1_->Get(Form("%s",histname1.Data()));
hist2 = (TH1D*)ifile2_->Get(Form("%s",histname2.Data()));
hist3 = (TH1D*)ifile2_->Get(Form("%s",histname3.Data()));
// }
TCanvas *myCan = new TCanvas("asd","asd");
myCan->Draw();
myCan->cd();
hist1->SetMarkerColor(kBlack);
hist2->SetMarkerColor(kRed);
hist3->SetMarkerColor(kBlue);
hist1->SetMarkerStyle(20);
hist2->SetMarkerStyle(20);
hist3->SetMarkerStyle(20);
hist1->SetMarkerSize(1.6);
hist2->SetMarkerSize(1.6);
hist3->SetMarkerSize(1.6);
hist1->SetTitle("");
hist1->GetXaxis()->SetRangeUser(0.,1.0);
hist2->GetXaxis()->SetRangeUser(0.,1.0);
hist3->GetXaxis()->SetRangeUser(0.,1.0);
gStyle->SetOptStat(0);
//hist1->Scale(1./1.026);//1.026, 1.05
hist1->SetMaximum(1.015);
hist1->SetMinimum(0.9425);
hist1->GetXaxis()->SetTitleSize(0.05);
hist1->GetYaxis()->SetTitleSize(0.05);
hist1->Draw("p");
hist2->Draw("psame");
hist3->Draw("psame");
TLatex Tl;
Tl.SetTextAlign(23);
Tl.SetTextSize(0.08);
Tl.SetTextColor(kBlack);
Tl.DrawLatex(0.5,0.995,"p#bar{p}");
Tl.SetTextColor(kRed);
Tl.DrawLatex(0.5,0.980,"p#bar{#Lambda}");
Tl.SetTextColor(kBlue);
Tl.DrawLatex(0.5,0.965,"#bar{p}#Lambda");
myCan->SaveAs(Form("comp%s.eps",sys.Data()));
}
void topTitle(const char *title)
{
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.11);
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.88, 0.70, title);
}
void cmsPrel(double intLumi){
TLatex latex;
latex.SetNDC();
latex.SetTextSize(0.04);
latex.SetTextAlign(31); // align right
//latex.DrawLatex(0.98,0.965,"#font[22]{#sqrt{s} = 7 TeV}");
latex.DrawLatex(0.98,0.965,"#sqrt{s} = 7 TeV");
if (intLumi > 0.) {
latex.SetTextAlign(31); // align right
latex.DrawLatex(0.98,0.88,Form("#int #font[12]{L}dt = %.1f nb^{-1}",intLumi));
}
latex.SetTextAlign(11); // align left
//latex.DrawLatex(0.02,0.965,"#font[22]{CMS preliminary 2010}");
latex.DrawLatex(0.02,0.965,"CMS preliminary 2010");
//latex.DrawLatex(0.20,0.20,"-0.4 < #eta < 0.4");
return;
}
void overlay_plots(const string& fFile0, const string& fFile1, const double fYmax, const string& fPlot, const string& fLegendEnt1, const string& fLegendEnt2, const string& fName) {
TProfile *p[2];
TFile file0(fFile0.c_str());
TDirectoryFile *subDir = (TDirectoryFile*)file0.Get("offsetAnalysis");
p[0] = (TProfile*)subDir->Get(fPlot.c_str());
TFile file1(fFile1.c_str());
subDir = (TDirectoryFile*)file1.Get("offsetAnalysis");
p[1] = (TProfile*)subDir->Get(fPlot.c_str());
p[0]->SetTitleOffset(1.5,"Y");
p[0]->GetXaxis()->SetTitleSize(0.04);
p[0]->GetYaxis()->SetTitleSize(0.04);
double ymax = ((p[0]->GetMaximum())>(p[1]->GetMaximum())) ? p[0]->GetMaximum() : p[1]->GetMaximum();
p[0]->GetYaxis()->SetRangeUser(0.,fYmax);
TCanvas *c = new TCanvas("c","",800,800);
c->cd();
p[0]->SetLineWidth(3);
p[0]->SetLineColor(kRed);
p[0]->SetFillColor(kRed);
p[0]->Draw("hist");
p[1]->SetLineWidth(3);
p[1]->SetLineColor(kBlack);
p[1]->SetMarkerStyle(20);
p[1]->SetMarkerColor(kBlack);
p[1]->Draw("sames");
TLegend *legend = new TLegend(.57,.77,.9,.87);
legend->SetBorderSize(1);
legend->SetFillColor(0);
// legend->SetFillStyle(0);
legend->SetMargin(0.12);
legend->AddEntry(p[0],fLegendEnt1.c_str(),"l");
legend->AddEntry(p[1],fLegendEnt2.c_str(),"l");
legend->Draw();
TLatex l;
l.SetTextAlign(12);
l.SetTextSize(0.04);
l.SetTextFont(62);
l.SetNDC();
l.DrawLatex(0.15,0.85,"CMS 2009 Preliminary");
string fileName = fName;
c->SetGridy();
c->SaveAs(fileName.c_str());
delete legend;
delete c;
}
void cmsPrelim(double intLumi, bool prelim)
{
TLatex latex;
latex.SetNDC();
latex.SetTextFont(62);
if(prelim) latex.SetTextSize(0.04);
else latex.SetTextSize(0.045);
latex.SetTextAlign(11); // align left
//if(prelim) latex.DrawLatex(0.13,0.92,"CMS Preliminary");
//else latex.DrawLatex(0.13,0.92,"CMS");
if(prelim) latex.DrawLatex(0.19,0.92,"CMS Preliminary");
else latex.DrawLatex(0.19,0.92,"CMS");
latex.SetTextAlign(31); // align right
//latex.DrawLatex(0.89, 0.92, Form("#sqrt{s} = 7 TeV, L_{int} = %4.2f fb^{-1}", intLumi));
latex.DrawLatex(0.92, 0.92, Form("#sqrt{s} = 8 TeV, L_{int} = %4.1f fb^{-1}", intLumi));
}
void
LEPStyle::cmsPrelim( double intLumi )
{
TLatex latex;
latex.SetNDC();
latex.SetTextSize(_scale*0.055);
latex.SetTextAlign(31); // align right
latex.DrawLatex(1-_d/_W,_txtY,_txtRight);
if (intLumi > 0.) {
latex.SetTextAlign(11); // align right
latex.DrawLatex(2.5*_g/_W,0.80,Form("#int #font[12]{L} dt = %.1f pb^{-1}",intLumi));
}
latex.SetTextAlign(11); // align left
latex.DrawLatex(_g/_W,_txtY,_txtLeft);
// latex.SetTextAlign(01);
// latex.DrawLatex(0.5,_txtY,_txtCenter);
}
void add_plot_label( char* label, double x, double y, double size = 0.05, int color = 1, int font = 62, int align = 22 ){
TLatex *latex = new TLatex( x, y, label );
latex->SetNDC();
latex->SetTextSize(size);
latex->SetTextColor(color);
latex->SetTextFont(font);
latex->SetTextAlign(align);
latex->Draw();
}
