bool PlotBundle::drawHistsAndGraphs(const PlotStyle &plot_style) const {
if (hasDrawables()) {
std::vector<DrawableDataObjectDrawOptionPair<TH1*> >::const_iterator hist_it;
for (hist_it = histograms.begin(); hist_it != histograms.end();
hist_it++) {
TH1* hist = hist_it->data_object;
std::string draw_option(hist_it->draw_option);
if (hist) {
if (hist->GetDimension() > 1) {
hist->GetZaxis()->SetLimits(plot_axis.z_axis_range.low,
plot_axis.z_axis_range.high);
hist->GetZaxis()->SetRangeUser(plot_axis.z_axis_range.low,
plot_axis.z_axis_range.high);
}
draw_option.append("SAME");
hist->Draw(draw_option.c_str());
}
}
std::vector<DrawableDataObjectDrawOptionPair<TGraph*> >::const_iterator graph_it;
for (graph_it = graphs.begin(); graph_it != graphs.end(); graph_it++) {
TGraph* graph = graph_it->data_object;
std::string draw_option(graph_it->draw_option);
if (graph) {
draw_option.append("SAME");
graph->Draw(draw_option.c_str());
}
}
std::vector<DrawableDataObjectDrawOptionPair<TGraph2D*> >::const_iterator graph2d_it;
for (graph2d_it = graphs2d.begin(); graph2d_it != graphs2d.end();
graph2d_it++) {
TGraph2D* graph = graph2d_it->data_object;
std::string draw_option(graph2d_it->draw_option);
if (graph) {
draw_option.append("SAME");
graph->Draw(draw_option.c_str());
}
}
return true;
} else {
std::cout
<< "Dude, you forgot to add a histogram or graph in the plot bundle..."
<< " Please add at least one and make sure it points to an existing object!"
<< std::endl;
return false;
}
}
void rf705_linearmorph()
{
// C r e a t e e n d p o i n t p d f s h a p e s
// ------------------------------------------------------
// Observable
RooRealVar x("x","x",-20,20) ;
// Lower end point shape: a Gaussian
RooRealVar g1mean("g1mean","g1mean",-10) ;
RooGaussian g1("g1","g1",x,g1mean,RooConst(2)) ;
// Upper end point shape: a Polynomial
RooPolynomial g2("g2","g2",x,RooArgSet(RooConst(-0.03),RooConst(-0.001))) ;
// C r e a t e i n t e r p o l a t i n g p d f
// -----------------------------------------------
// Create interpolation variable
RooRealVar alpha("alpha","alpha",0,1.0) ;
// Specify sampling density on observable and interpolation variable
x.setBins(1000,"cache") ;
alpha.setBins(50,"cache") ;
// Construct interpolating pdf in (x,a) represent g1(x) at a=a_min
// and g2(x) at a=a_max
RooIntegralMorph lmorph("lmorph","lmorph",g1,g2,x,alpha) ;
// P l o t i n t e r p o l a t i n g p d f a t v a r i o u s a l p h a
// -----------------------------------------------------------------------------
// Show end points as blue curves
RooPlot* frame1 = x.frame() ;
g1.plotOn(frame1) ;
g2.plotOn(frame1) ;
// Show interpolated shapes in red
alpha.setVal(0.125) ;
lmorph.plotOn(frame1,LineColor(kRed)) ;
alpha.setVal(0.25) ;
lmorph.plotOn(frame1,LineColor(kRed)) ;
alpha.setVal(0.375) ;
lmorph.plotOn(frame1,LineColor(kRed)) ;
alpha.setVal(0.50) ;
lmorph.plotOn(frame1,LineColor(kRed)) ;
alpha.setVal(0.625) ;
lmorph.plotOn(frame1,LineColor(kRed)) ;
alpha.setVal(0.75) ;
lmorph.plotOn(frame1,LineColor(kRed)) ;
alpha.setVal(0.875) ;
lmorph.plotOn(frame1,LineColor(kRed)) ;
alpha.setVal(0.95) ;
lmorph.plotOn(frame1,LineColor(kRed)) ;
// S h o w 2 D d i s t r i b u t i o n o f p d f ( x , a l p h a )
// -----------------------------------------------------------------------
// Create 2D histogram
TH1* hh = lmorph.createHistogram("hh",x,Binning(40),YVar(alpha,Binning(40))) ;
hh->SetLineColor(kBlue) ;
// F i t p d f t o d a t a s e t w i t h a l p h a = 0 . 8
// -----------------------------------------------------------------
// Generate a toy dataset at alpha = 0.8
alpha=0.8 ;
RooDataSet* data = lmorph.generate(x,1000) ;
// Fit pdf to toy data
lmorph.setCacheAlpha(kTRUE) ;
lmorph.fitTo(*data,Verbose(kTRUE)) ;
// Plot fitted pdf and data overlaid
RooPlot* frame2 = x.frame(Bins(100)) ;
data->plotOn(frame2) ;
lmorph.plotOn(frame2) ;
// S c a n - l o g ( L ) v s a l p h a
// -----------------------------------------
// Show scan -log(L) of dataset w.r.t alpha
RooPlot* frame3 = alpha.frame(Bins(100),Range(0.1,0.9)) ;
// Make 2D pdf of histogram
RooNLLVar nll("nll","nll",lmorph,*data) ;
nll.plotOn(frame3,ShiftToZero()) ;
lmorph.setCacheAlpha(kFALSE) ;
TCanvas* c = new TCanvas("rf705_linearmorph","rf705_linearmorph",800,800) ;
c->Divide(2,2) ;
c->cd(1) ; gPad->SetLeftMargin(0.15) ; frame1->GetYaxis()->SetTitleOffset(1.6) ; frame1->Draw() ;
c->cd(2) ; gPad->SetLeftMargin(0.20) ; hh->GetZaxis()->SetTitleOffset(2.5) ; hh->Draw("surf") ;
c->cd(3) ; gPad->SetLeftMargin(0.15) ; frame3->GetYaxis()->SetTitleOffset(1.4) ; frame3->Draw() ;
c->cd(4) ; gPad->SetLeftMargin(0.15) ; frame2->GetYaxis()->SetTitleOffset(1.4) ; frame2->Draw() ;
return ;
}
void drawCtauMass2DPlot(RooWorkspace& myws, // Local workspace
string outputDir, // Output directory
struct KinCuts cut, // Variable with current kinematic cuts
string plotLabel, // The label used to define the output file name
// Select the type of datasets to fit
string DSTAG, // Specifies the type of datasets: i.e, DATA, MCJPSINP, ...
bool isPbPb, // Define if it is PbPb (True) or PP (False)\
// Select the drawing options
map<string, double> binWidth={} // User-defined Location of the fit results
)
{
gStyle->SetOptStat(0);
if (DSTAG.find("_")!=std::string::npos) DSTAG.erase(DSTAG.find("_"));
double minRangeCtau = -0.5;
double maxRangeCtau = 2.0;
int nBinsCtau = min(int( round((maxRangeCtau - minRangeCtau)/binWidth["CTAU"]*2) ), 1000);
// myws.var("ctau")->setBin(nBinsCtau, Binning(nBinsCtau, minRangeCtau, maxRangeCtau));
double minRangeMass = cut.dMuon.M.Min;
double maxRangeMass = cut.dMuon.M.Max;
int nBinsMass = min(int( round((maxRangeMass - minRangeMass)/binWidth["MASS"]) ), 1000);
// myws.var("invMass")->setBin(nBinsCtau, Binning(nBinsCtau, minRangeCtau, maxRangeCtau));
string pdfTotName = Form("pdfCTAUMASS_Tot_%s", (isPbPb?"PbPb":"PP"));
TH1* hPDF = ((RooAbsReal*)myws.pdf(pdfTotName.c_str()))->createHistogram("PDF 2D",*myws.var("ctau"), Extended(kTRUE), Binning(nBinsCtau, minRangeCtau, maxRangeCtau), YVar(*myws.var("invMass"), Binning(nBinsMass, minRangeMass, maxRangeMass)));
string dsOSName = Form("dOS_%s_%s", DSTAG.c_str(), (isPbPb?"PbPb":"PP"));
TH1* hDATA = ((RooDataSet*)myws.data(dsOSName.c_str()))->createHistogram("DATA 2D",*myws.var("ctau"), Binning(nBinsCtau, minRangeCtau, maxRangeCtau), YVar(*myws.var("invMass"), Binning(nBinsMass, minRangeMass, maxRangeMass)));
// Create the main canvas
TCanvas *cFigPDF = new TCanvas(Form("cCtauMassPDF_%s", (isPbPb?"PbPb":"PP")), "cCtauMassPDF",2000,2000);
cFigPDF->cd();
hPDF->GetYaxis()->CenterTitle(kTRUE);
hPDF->GetYaxis()->SetTitleOffset(2.1);
hPDF->GetYaxis()->SetTitleSize(0.035);
hPDF->GetYaxis()->SetLabelSize(0.025);
hPDF->GetYaxis()->SetTitle("Mass [GeV/c]");
hPDF->GetXaxis()->CenterTitle(kTRUE);
hPDF->GetXaxis()->SetTitleOffset(2.1);
hPDF->GetXaxis()->SetTitleSize(0.035);
hPDF->GetXaxis()->SetLabelSize(0.025);
hPDF->GetXaxis()->SetTitle("#font[12]{l}_{J/#psi} (mm)");
hPDF->GetZaxis()->SetTitleOffset(2.0);
hPDF->GetZaxis()->SetTitleSize(0.035);
hPDF->GetZaxis()->SetLabelSize(0.025);
hPDF->GetYaxis()->SetRangeUser(minRangeMass, maxRangeMass);
hPDF->GetXaxis()->SetRangeUser(minRangeCtau, maxRangeCtau);
cFigPDF->SetLogz(kTRUE);
hPDF->Draw("LEGO2");
gSystem->mkdir(Form("%sctauMass/%s/plot/pdf2D/", outputDir.c_str(), DSTAG.c_str()), kTRUE);
cFigPDF->SaveAs(Form("%sctauMass/%s/plot/pdf2D/PLOT_%s_%s_%s%s_pt%.0f%.0f_rap%.0f%.0f_cent%d%d.pdf", outputDir.c_str(), DSTAG.c_str(), "CTAUMASSPDF", DSTAG.c_str(), (isPbPb?"PbPb":"PP"), plotLabel.c_str(), (cut.dMuon.Pt.Min*10.0), (cut.dMuon.Pt.Max*10.0), (cut.dMuon.AbsRap.Min*10.0), (cut.dMuon.AbsRap.Max*10.0), cut.Centrality.Start, cut.Centrality.End));
cFigPDF->Clear();
cFigPDF->Close();
// Create the main canvas
TCanvas *cFigDATA = new TCanvas(Form("cCtauMassPDF_%s", (isPbPb?"PbPb":"PP")), "cCtauMassPDF",2000,2000);
cFigDATA->cd();
hDATA->GetYaxis()->CenterTitle(kTRUE);
hDATA->GetYaxis()->SetTitleOffset(2.1);
hDATA->GetYaxis()->SetTitleSize(0.035);
hDATA->GetYaxis()->SetLabelSize(0.025);
hDATA->GetYaxis()->SetTitle("Mass [GeV/c]");
hDATA->GetXaxis()->CenterTitle(kTRUE);
hDATA->GetXaxis()->SetTitleOffset(2.1);
hDATA->GetXaxis()->SetTitleSize(0.035);
hDATA->GetXaxis()->SetLabelSize(0.025);
hDATA->GetXaxis()->SetTitle("#font[12]{l}_{J/#psi} (mm)");
hDATA->GetZaxis()->SetTitleOffset(2.0);
hDATA->GetZaxis()->SetTitleSize(0.035);
hDATA->GetZaxis()->SetLabelSize(0.025);
hDATA->GetYaxis()->SetRangeUser(minRangeMass, maxRangeMass);
hDATA->GetXaxis()->SetRangeUser(minRangeCtau, maxRangeCtau);
cFigDATA->SetLogz(kTRUE);
hDATA->Draw("LEGO2");
gSystem->mkdir(Form("%sctauMass/%s/plot/pdf2D/", outputDir.c_str(), DSTAG.c_str()), kTRUE);
cFigDATA->SaveAs(Form("%sctauMass/%s/plot/pdf2D/PLOT_%s_%s_%s%s_pt%.0f%.0f_rap%.0f%.0f_cent%d%d.pdf", outputDir.c_str(), DSTAG.c_str(), "CTAUMASSDATA", DSTAG.c_str(), (isPbPb?"PbPb":"PP"), plotLabel.c_str(), (cut.dMuon.Pt.Min*10.0), (cut.dMuon.Pt.Max*10.0), (cut.dMuon.AbsRap.Min*10.0), (cut.dMuon.AbsRap.Max*10.0), cut.Centrality.Start, cut.Centrality.End));
cFigDATA->Clear();
cFigDATA->Close();
delete hPDF;
delete hDATA;
};
void csv(TString input="tmva.csvoutput.txt", TString par1="par2", TString par2="par3", TString par3="", TString value="eventEffScaled_5") {
std::cout << "Usage:" << std::endl
<< ".x scripts/csv.C with default arguments" << std::endl
<< ".x scripts/csv.C(filename, par1, par2, value)" << std::endl
<< std::endl
<< " Optional arguments:" << std::endl
<< " filename path to CSV file" << std::endl
<< " par1 name of X-parameter branch" << std::endl
<< " par2 name of Y-parameter branch (if empty, efficiency is drawn as a function of par1)" << std::endl
<< " value name of result (efficiency) branch" << std::endl
<< std::endl;
TTree *tree = new TTree("data", "data");
tree->ReadFile(input);
gStyle->SetPalette(1);
gStyle->SetPadRightMargin(0.14);
TCanvas *canvas = new TCanvas("csvoutput", "CSV Output", 1200, 900);
tree->SetMarkerStyle(kFullDotMedium);
tree->SetMarkerColor(kRed);
if(par2.Length() > 0) {
//tree->Draw(Form("%s:%s", par2.Data(), par1.Data()));
if(par3.Length() > 0)
tree->Draw(Form("%s:%s:%s:%s", par1.Data(), par2.Data(), par3.Data(), value.Data()), "", "COLZ"); //, "", "Z");
else
tree->Draw(Form("%s:%s:%s", par2.Data(), par1.Data(), value.Data()), "", "COLZ"); //, "", "Z");
TH1 *histo = tree->GetHistogram();
if(!histo)
return;
histo->SetTitle(Form("%s with different classifier parameters", value.Data()));
histo->GetXaxis()->SetTitle(Form("Classifier parameter %s", par1.Data()));
histo->GetYaxis()->SetTitle(Form("Classifier parameter %s", par2.Data()));
if(par3.Length() > 0)
histo->GetZaxis()->SetTitle(Form("Classifier parameter %s", par3.Data()));
else
histo->GetZaxis()->SetTitle("");
if(par3.Length() == 0) {
float x = 0;
float y = 0;
float val = 0;
double maxVal = tree->GetMaximum(value);
double minVal = tree->GetMinimum(value);
tree->SetBranchAddress(par1, &x);
tree->SetBranchAddress(par2, &y);
tree->SetBranchAddress(value, &val);
TLatex l;
l.SetTextSize(0.03);
Long64_t nentries = tree->GetEntries();
for(Long64_t entry=0; entry < nentries; ++entry) {
tree->GetEntry(entry);
l.SetTextColor(textColor(val, maxVal, minVal));
l.DrawLatex(x, y, Form("%.3f", val*100));
}
}
}
else {
tree->Draw(Form("%s:%s", value.Data(), par1.Data()));
TH1 *histo = tree->GetHistogram();
if(!histo)
return;
histo->SetTitle(Form("%s with different classifier parameters", value.Data()));
histo->GetXaxis()->SetTitle(Form("Classifier parameter %s", par1.Data()));
histo->GetYaxis()->SetTitle(value);
}
}
void trigeff(){
gROOT->Reset();
//=========Macro generated from canvas: cEvtSel/EvtSel
//========= (Wed Apr 14 17:19:34 2010) by ROOT version5.22/00d
TCanvas *cEvtSel = new TCanvas("cEvtSel", "EvtSel",550,600);
/*
TCanvas *cEvtSel = new TCanvas("cEvtSel", "EvtSel",0,22,550,600);
gStyle->SetOptFit(1);
gStyle->SetOptStat(0);
cEvtSel->Range(-13.62658,-0.215625,63.58861,1.221875);
cEvtSel->SetFillColor(0);
cEvtSel->SetBorderMode(0);
cEvtSel->SetBorderSize(0);
cEvtSel->SetTickx(1);
cEvtSel->SetTicky(1);
cEvtSel->SetLeftMargin(0.17);
cEvtSel->SetRightMargin(0.04);
cEvtSel->SetTopMargin(0.05);
cEvtSel->SetBottomMargin(0.15);
cEvtSel->SetFrameLineColor(0);
cEvtSel->SetFrameBorderMode(0);
cEvtSel->SetFrameLineColor(0);
cEvtSel->SetFrameBorderMode(0);
*/
TH1 *hFrame1 = new TH2D("hFrame1","",1,-0.5,60.5,1,0,1.15);
hFrame1->SetDirectory(0);
hFrame1->SetStats(0);
hFrame1->GetXaxis()->SetTitle("Charged-particle multiplicity");
hFrame1->GetXaxis()->CenterTitle(true);
hFrame1->GetYaxis()->SetTitle("Selection efficiency");
hFrame1->GetYaxis()->CenterTitle(true);
hFrame1->Draw("");
hFrame1->GetXaxis()->SetNdivisions(312);
/*
TH1 *hFrame1 = new TH2D("hFrame1","",1,-0.5,60.5,1,0,1.15);
hFrame1->SetDirectory(0);
hFrame1->SetStats(0);
hFrame1->SetFillColor(1);
hFrame1->SetFillStyle(0);
hFrame1->SetLineStyle(0);
hFrame1->SetMarkerStyle(20);
hFrame1->SetMarkerSize(1.5);
hFrame1->GetXaxis()->SetTitle("Charged-particle multiplicity");
hFrame1->GetXaxis()->CenterTitle(true);
hFrame1->GetXaxis()->SetNdivisions(312);
hFrame1->GetXaxis()->SetLabelFont(42);
hFrame1->GetXaxis()->SetLabelOffset(0.01);
hFrame1->GetXaxis()->SetLabelSize(0.045);
hFrame1->GetXaxis()->SetTitleSize(0.055);
hFrame1->GetXaxis()->SetTitleFont(42);
hFrame1->GetYaxis()->SetTitle("Selection efficiency");
hFrame1->GetYaxis()->CenterTitle(true);
hFrame1->GetYaxis()->SetLabelFont(42);
hFrame1->GetYaxis()->SetLabelOffset(0.01);
hFrame1->GetYaxis()->SetLabelSize(0.045);
hFrame1->GetYaxis()->SetTitleSize(0.055);
hFrame1->GetYaxis()->SetTitleOffset(1.6);
hFrame1->GetYaxis()->SetTitleFont(42);
hFrame1->GetZaxis()->SetLabelFont(42);
hFrame1->GetZaxis()->SetLabelSize(0.045);
hFrame1->GetZaxis()->SetTitleFont(42);
hFrame1->Draw("");
*/
TH1 *trigeff = new TH1D("trigeff","trigeff",200,0,200);
trigeff->SetBinContent(1,0.007879584);
trigeff->SetBinContent(2,0.3524112);
trigeff->SetBinContent(3,0.5050529);
trigeff->SetBinContent(4,0.5978595);
trigeff->SetBinContent(5,0.6587795);
trigeff->SetBinContent(6,0.7034895);
trigeff->SetBinContent(7,0.7451056);
trigeff->SetBinContent(8,0.7911818);
trigeff->SetBinContent(9,0.8363224);
trigeff->SetBinContent(10,0.873251);
trigeff->SetBinContent(11,0.9012649);
trigeff->SetBinContent(12,0.9177872);
trigeff->SetBinContent(13,0.9395182);
trigeff->SetBinContent(14,0.9502278);
trigeff->SetBinContent(15,0.9600504);
trigeff->SetBinContent(16,0.9667429);
trigeff->SetBinContent(17,0.9722969);
trigeff->SetBinContent(18,0.977478);
trigeff->SetBinContent(19,0.9801504);
trigeff->SetBinContent(20,0.9840849);
trigeff->SetBinContent(21,0.9843521);
trigeff->SetBinContent(22,0.9856779);
trigeff->SetBinContent(23,0.9894354);
trigeff->SetBinContent(24,0.9892805);
trigeff->SetBinContent(25,0.9913111);
trigeff->SetBinContent(26,0.9934833);
trigeff->SetBinContent(27,0.9907814);
trigeff->SetBinContent(28,0.9928126);
trigeff->SetBinContent(29,0.990099);
trigeff->SetBinContent(30,0.9918119);
trigeff->SetBinContent(31,0.994504);
trigeff->SetBinContent(32,0.9938513);
trigeff->SetBinContent(33,0.993073);
trigeff->SetBinContent(34,0.9965539);
trigeff->SetBinContent(35,0.9955737);
trigeff->SetBinContent(36,0.9952005);
trigeff->SetBinContent(37,0.9940288);
trigeff->SetBinContent(38,0.9946929);
trigeff->SetBinContent(39,0.9973046);
trigeff->SetBinContent(40,0.9959936);
trigeff->SetBinContent(41,0.9964897);
trigeff->SetBinContent(42,0.9943966);
trigeff->SetBinContent(43,0.9940175);
trigeff->SetBinContent(44,0.9947941);
trigeff->SetBinContent(45,0.9960804);
trigeff->SetBinContent(46,0.9970282);
trigeff->SetBinContent(47,0.9963936);
trigeff->SetBinContent(48,0.9968187);
trigeff->SetBinContent(49,0.9956733);
trigeff->SetBinContent(50,0.996085);
trigeff->SetBinContent(51,0.9988466);
trigeff->SetBinContent(52,0.9982079);
trigeff->SetBinContent(53,0.99625);
trigeff->SetBinContent(54,0.9962802);
trigeff->SetBinContent(55,0.9980328);
trigeff->SetBinContent(56,0.9960396);
trigeff->SetBinContent(57,0.9993271);
trigeff->SetBinContent(58,0.9978261);
trigeff->SetBinContent(59,0.9950843);
trigeff->SetBinContent(60,0.9992459);
trigeff->SetBinContent(61,0.9977925);
trigeff->SetBinContent(62,0.9960599);
trigeff->SetBinContent(63,0.9952456);
trigeff->SetBinContent(64,0.9960661);
trigeff->SetBinContent(65,0.9965428);
trigeff->SetBinContent(66,0.995094);
trigeff->SetBinContent(67,0.9962441);
trigeff->SetBinContent(68,0.9990842);
trigeff->SetBinContent(69,0.9981413);
trigeff->SetBinContent(70,0.9981007);
trigeff->SetBinContent(71,0.9980676);
trigeff->SetBinContent(72,0.9990512);
trigeff->SetBinContent(73,0.9969636);
trigeff->SetBinContent(74,0.9959432);
trigeff->SetBinContent(75,0.9969008);
trigeff->SetBinContent(76,0.994929);
trigeff->SetBinContent(77,0.9978425);
trigeff->SetBinContent(78,0.9989396);
trigeff->SetBinContent(79,0.9955996);
trigeff->SetBinContent(80,0.9988975);
trigeff->SetBinContent(81,0.998829);
trigeff->SetBinContent(82,0.9960682);
trigeff->SetBinContent(83,0.9963986);
trigeff->SetBinContent(84,0.995098);
trigeff->SetBinContent(85,0.9962169);
trigeff->SetBinContent(86,0.9973615);
trigeff->SetBinContent(87,0.9987562);
trigeff->SetBinContent(88,0.9987212);
trigeff->SetBinContent(89,0.9987374);
trigeff->SetBinContent(90,0.9986264);
trigeff->SetBinContent(91,0.9972973);
trigeff->SetBinContent(92,1);
trigeff->SetBinContent(93,0.9959459);
trigeff->SetBinContent(94,0.9945504);
trigeff->SetBinContent(95,0.9985612);
trigeff->SetBinContent(96,1);
trigeff->SetBinContent(97,1);
trigeff->SetBinContent(98,0.9968454);
trigeff->SetBinContent(99,0.9983471);
trigeff->SetBinContent(100,0.9968);
trigeff->SetBinContent(101,0.9983871);
trigeff->SetBinContent(102,0.9966044);
trigeff->SetBinContent(103,0.9984);
trigeff->SetBinContent(104,1);
trigeff->SetBinContent(105,0.9946996);
trigeff->SetBinContent(106,0.9946333);
trigeff->SetBinContent(107,0.9981343);
trigeff->SetBinContent(108,1);
trigeff->SetBinContent(109,0.9962963);
trigeff->SetBinContent(110,0.9923225);
trigeff->SetBinContent(111,1);
trigeff->SetBinContent(112,0.9979633);
trigeff->SetBinContent(113,0.9978166);
trigeff->SetBinContent(114,1);
trigeff->SetBinContent(115,0.9980237);
trigeff->SetBinContent(116,0.9952494);
trigeff->SetBinContent(117,0.9979253);
trigeff->SetBinContent(118,1);
trigeff->SetBinContent(119,1);
trigeff->SetBinContent(120,0.9950125);
trigeff->SetBinContent(121,0.9975728);
trigeff->SetBinContent(122,0.9921466);
trigeff->SetBinContent(123,0.9974874);
trigeff->SetBinContent(124,1);
trigeff->SetBinContent(125,0.9974874);
trigeff->SetBinContent(126,0.9969419);
trigeff->SetBinContent(127,0.994723);
trigeff->SetBinContent(128,0.9940299);
trigeff->SetBinContent(129,1);
trigeff->SetBinContent(130,0.9966667);
trigeff->SetBinContent(131,1);
trigeff->SetBinContent(132,1);
trigeff->SetBinContent(133,0.9931507);
trigeff->SetBinContent(134,0.9963504);
trigeff->SetBinContent(135,1);
trigeff->SetBinContent(136,0.9959677);
trigeff->SetBinContent(137,1);
trigeff->SetBinContent(138,1);
trigeff->SetBinContent(139,1);
trigeff->SetBinContent(140,0.9959514);
trigeff->SetBinContent(141,1);
trigeff->SetBinContent(142,0.9955556);
trigeff->SetBinContent(143,1);
trigeff->SetBinContent(144,1);
trigeff->SetBinContent(145,1);
trigeff->SetBinContent(146,0.989418);
trigeff->SetBinContent(147,0.9940828);
trigeff->SetBinContent(148,1);
trigeff->SetBinContent(149,0.9942857);
trigeff->SetBinContent(150,1);
trigeff->SetBinContent(151,1);
trigeff->SetBinContent(152,1);
trigeff->SetBinContent(153,1);
trigeff->SetBinContent(154,1);
trigeff->SetBinContent(155,0.9923077);
trigeff->SetBinContent(156,1);
trigeff->SetBinContent(157,1);
trigeff->SetBinContent(158,1);
trigeff->SetBinContent(159,1);
trigeff->SetBinContent(160,0.9901961);
trigeff->SetBinContent(161,1);
trigeff->SetBinContent(162,1);
trigeff->SetBinContent(163,1);
trigeff->SetBinContent(164,1);
trigeff->SetBinContent(165,0.9767442);
trigeff->SetBinContent(166,1);
trigeff->SetBinContent(167,1);
trigeff->SetBinContent(168,1);
trigeff->SetBinContent(169,1);
trigeff->SetBinContent(170,1);
trigeff->SetBinContent(171,1);
trigeff->SetBinContent(172,1);
trigeff->SetBinContent(173,0.9807692);
trigeff->SetBinContent(174,1);
trigeff->SetBinContent(175,0.9850746);
trigeff->SetBinContent(176,1);
trigeff->SetBinContent(177,1);
trigeff->SetBinContent(178,1);
trigeff->SetBinContent(179,1);
trigeff->SetBinContent(180,1);
trigeff->SetBinContent(181,1);
trigeff->SetBinContent(182,1);
trigeff->SetBinContent(183,1);
trigeff->SetBinContent(184,1);
trigeff->SetBinContent(185,1);
trigeff->SetBinContent(186,1);
trigeff->SetBinContent(187,1);
trigeff->SetBinContent(188,1);
trigeff->SetBinContent(189,1);
trigeff->SetBinContent(190,1);
trigeff->SetBinContent(191,1);
trigeff->SetBinContent(192,1);
trigeff->SetBinContent(193,1);
trigeff->SetBinContent(194,1);
trigeff->SetBinContent(195,0.9333333);
trigeff->SetBinContent(196,1);
trigeff->SetBinContent(197,0.9230769);
trigeff->SetBinContent(198,1);
trigeff->SetBinContent(199,1);
trigeff->SetBinContent(200,1);
trigeff->SetBinContent(201,1);
trigeff->SetBinError(1,0.0006283812);
trigeff->SetBinError(2,0.005498682);
trigeff->SetBinError(3,0.0054516);
trigeff->SetBinError(4,0.005047444);
trigeff->SetBinError(5,0.004707652);
trigeff->SetBinError(6,0.004447388);
trigeff->SetBinError(7,0.004290336);
trigeff->SetBinError(8,0.003953715);
trigeff->SetBinError(9,0.003594607);
trigeff->SetBinError(10,0.003181503);
trigeff->SetBinError(11,0.002795854);
trigeff->SetBinError(12,0.002534088);
trigeff->SetBinError(13,0.002203234);
trigeff->SetBinError(14,0.002035759);
trigeff->SetBinError(15,0.001857668);
trigeff->SetBinError(16,0.001750359);
trigeff->SetBinError(17,0.001668633);
trigeff->SetBinError(18,0.001566698);
trigeff->SetBinError(19,0.001548763);
trigeff->SetBinError(20,0.001441238);
trigeff->SetBinError(21,0.001500853);
trigeff->SetBinError(22,0.001498831);
trigeff->SetBinError(23,0.001345494);
trigeff->SetBinError(24,0.001388057);
trigeff->SetBinError(25,0.001289631);
trigeff->SetBinError(26,0.001166616);
trigeff->SetBinError(27,0.001415891);
trigeff->SetBinError(28,0.0013075);
trigeff->SetBinError(29,0.001557715);
trigeff->SetBinError(30,0.001464586);
trigeff->SetBinError(31,0.001225565);
trigeff->SetBinError(32,0.001306868);
trigeff->SetBinError(33,0.001471707);
trigeff->SetBinError(34,0.001037251);
trigeff->SetBinError(35,0.001224911);
trigeff->SetBinError(36,0.001279626);
trigeff->SetBinError(37,0.001443898);
trigeff->SetBinError(38,0.0014146);
trigeff->SetBinError(39,0.001017398);
trigeff->SetBinError(40,0.001264398);
trigeff->SetBinError(41,0.001238902);
trigeff->SetBinError(42,0.001549757);
trigeff->SetBinError(43,0.00165428);
trigeff->SetBinError(44,0.001565537);
trigeff->SetBinError(45,0.001383086);
trigeff->SetBinError(46,0.001211415);
trigeff->SetBinError(47,0.001360627);
trigeff->SetBinError(48,0.001296707);
trigeff->SetBinError(49,0.001526394);
trigeff->SetBinError(50,0.001476827);
trigeff->SetBinError(51,0.0008151083);
trigeff->SetBinError(52,0.00103375);
trigeff->SetBinError(53,0.001528058);
trigeff->SetBinError(54,0.001515765);
trigeff->SetBinError(55,0.001134653);
trigeff->SetBinError(56,0.00161362);
trigeff->SetBinError(57,0.000672721);
trigeff->SetBinError(58,0.001253744);
trigeff->SetBinError(59,0.001853399);
trigeff->SetBinError(60,0.0007538634);
trigeff->SetBinError(61,0.001273096);
trigeff->SetBinError(62,0.001758596);
trigeff->SetBinError(63,0.001936339);
trigeff->SetBinError(64,0.001755834);
trigeff->SetBinError(65,0.001725618);
trigeff->SetBinError(66,0.001997934);
trigeff->SetBinError(67,0.001874404);
trigeff->SetBinError(68,0.0009153315);
trigeff->SetBinError(69,0.001313103);
trigeff->SetBinError(70,0.001341757);
trigeff->SetBinError(71,0.001365069);
trigeff->SetBinError(72,0.0009483164);
trigeff->SetBinError(73,0.001750424);
trigeff->SetBinError(74,0.002024279);
trigeff->SetBinError(75,0.001786534);
trigeff->SetBinError(76,0.00226206);
trigeff->SetBinError(77,0.001523934);
trigeff->SetBinError(78,0.001059883);
trigeff->SetBinError(79,0.002195374);
trigeff->SetBinError(80,0.001101928);
trigeff->SetBinError(81,0.001170274);
trigeff->SetBinError(82,0.002265586);
trigeff->SetBinError(83,0.002075545);
trigeff->SetBinError(84,0.002444966);
trigeff->SetBinError(85,0.00218004);
trigeff->SetBinError(86,0.001863254);
trigeff->SetBinError(87,0.001243007);
trigeff->SetBinError(88,0.001277954);
trigeff->SetBinError(89,0.001261829);
trigeff->SetBinError(90,0.001372683);
trigeff->SetBinError(91,0.001908515);
trigeff->SetBinError(93,0.00233586);
trigeff->SetBinError(94,0.002717361);
trigeff->SetBinError(95,0.001437813);
trigeff->SetBinError(98,0.0022271);
trigeff->SetBinError(99,0.001651526);
trigeff->SetBinError(100,0.002259118);
trigeff->SetBinError(101,0.001611602);
trigeff->SetBinError(102,0.002396962);
trigeff->SetBinError(103,0.001598719);
trigeff->SetBinError(105,0.00305204);
trigeff->SetBinError(106,0.003090155);
trigeff->SetBinError(107,0.00186393);
trigeff->SetBinError(109,0.00261406);
trigeff->SetBinError(110,0.003824007);
trigeff->SetBinError(112,0.002034585);
trigeff->SetBinError(113,0.002181021);
trigeff->SetBinError(115,0.001974331);
trigeff->SetBinError(116,0.003351189);
trigeff->SetBinError(117,0.002072536);
trigeff->SetBinError(120,0.003517911);
trigeff->SetBinError(121,0.002424237);
trigeff->SetBinError(122,0.004516325);
trigeff->SetBinError(123,0.002509404);
trigeff->SetBinError(125,0.002509404);
trigeff->SetBinError(126,0.003053424);
trigeff->SetBinError(127,0.003721576);
trigeff->SetBinError(128,0.004208913);
trigeff->SetBinError(130,0.003327773);
trigeff->SetBinError(133,0.004826582);
trigeff->SetBinError(134,0.003642969);
trigeff->SetBinError(136,0.00402412);
trigeff->SetBinError(140,0.004040379);
trigeff->SetBinError(142,0.004434557);
trigeff->SetBinError(146,0.007442916);
trigeff->SetBinError(147,0.005899627);
trigeff->SetBinError(149,0.005697936);
trigeff->SetBinError(155,0.007662665);
trigeff->SetBinError(160,0.009755745);
trigeff->SetBinError(165,0.01625201);
trigeff->SetBinError(173,0.01904496);
trigeff->SetBinError(175,0.01481357);
trigeff->SetBinError(195,0.06440612);
trigeff->SetBinError(197,0.0739053);
trigeff->SetEntries(385140);
trigeff->SetDirectory(0);
trigeff->SetFillColor(1);
trigeff->SetFillStyle(0);
trigeff->SetLineStyle(0);
trigeff->SetLineWidth(2);
trigeff->SetMarkerStyle(20);
trigeff->SetMarkerSize(1.5);
trigeff->GetXaxis()->SetTitle("M");
trigeff->GetXaxis()->SetLabelFont(42);
trigeff->GetXaxis()->SetLabelOffset(0.01);
trigeff->GetXaxis()->SetLabelSize(0.045);
trigeff->GetXaxis()->SetTitleSize(0.055);
trigeff->GetXaxis()->SetTitleFont(42);
trigeff->GetYaxis()->SetTitle("Selection efficiency");
trigeff->GetYaxis()->SetLabelFont(42);
trigeff->GetYaxis()->SetLabelOffset(0.01);
trigeff->GetYaxis()->SetLabelSize(0.045);
trigeff->GetYaxis()->SetTitleSize(0.055);
trigeff->GetYaxis()->SetTitleOffset(1.6);
trigeff->GetYaxis()->SetTitleFont(42);
trigeff->GetZaxis()->SetLabelFont(42);
trigeff->GetZaxis()->SetLabelSize(0.045);
trigeff->GetZaxis()->SetTitleFont(42);
trigeff->Draw("hist same");
TH1 *trigeff = new TH1D("trigeff","trigeff",200,0,200);
trigeff->SetBinContent(1,0.009007478);
trigeff->SetBinContent(2,0.2935561);
trigeff->SetBinContent(3,0.4702609);
trigeff->SetBinContent(4,0.6518799);
trigeff->SetBinContent(5,0.7428357);
trigeff->SetBinContent(6,0.8347119);
trigeff->SetBinContent(7,0.8777911);
trigeff->SetBinContent(8,0.904666);
trigeff->SetBinContent(9,0.9195271);
trigeff->SetBinContent(10,0.9292373);
trigeff->SetBinContent(11,0.9375605);
trigeff->SetBinContent(12,0.9489987);
trigeff->SetBinContent(13,0.9535653);
trigeff->SetBinContent(14,0.9493052);
trigeff->SetBinContent(15,0.9572732);
trigeff->SetBinContent(16,0.9589076);
trigeff->SetBinContent(17,0.9627003);
trigeff->SetBinContent(18,0.9629964);
trigeff->SetBinContent(19,0.9683917);
trigeff->SetBinContent(20,0.9674429);
trigeff->SetBinContent(21,0.9771153);
trigeff->SetBinContent(22,0.9762201);
trigeff->SetBinContent(23,0.9766491);
trigeff->SetBinContent(24,0.9828982);
trigeff->SetBinContent(25,0.9802921);
trigeff->SetBinContent(26,0.9826494);
trigeff->SetBinContent(27,0.9759428);
trigeff->SetBinContent(28,0.9843614);
trigeff->SetBinContent(29,0.9879452);
trigeff->SetBinContent(30,0.9855582);
trigeff->SetBinContent(31,0.9875656);
trigeff->SetBinContent(32,0.9883436);
trigeff->SetBinContent(33,0.9886389);
trigeff->SetBinContent(34,0.9900587);
trigeff->SetBinContent(35,0.9925793);
trigeff->SetBinContent(36,0.9926954);
trigeff->SetBinContent(37,0.9905261);
trigeff->SetBinContent(38,0.993372);
trigeff->SetBinContent(39,0.9899301);
trigeff->SetBinContent(40,0.9913687);
trigeff->SetBinContent(41,0.9920382);
trigeff->SetBinContent(42,0.9901091);
trigeff->SetBinContent(43,0.9941278);
trigeff->SetBinContent(44,0.9865421);
trigeff->SetBinContent(45,0.9902138);
trigeff->SetBinContent(46,0.9965451);
trigeff->SetBinContent(47,0.9961106);
trigeff->SetBinContent(48,0.9951392);
trigeff->SetBinContent(49,0.9951969);
trigeff->SetBinContent(50,0.9960957);
trigeff->SetBinContent(51,0.9981439);
trigeff->SetBinContent(52,0.9962282);
trigeff->SetBinContent(53,0.9973698);
trigeff->SetBinContent(54,0.9976959);
trigeff->SetBinContent(55,0.9933735);
trigeff->SetBinContent(56,0.9956897);
trigeff->SetBinContent(57,0.9953642);
trigeff->SetBinContent(58,0.9902303);
trigeff->SetBinContent(59,0.9985925);
trigeff->SetBinContent(60,0.9986178);
trigeff->SetBinContent(61,0.9991604);
trigeff->SetBinContent(62,0.9983897);
trigeff->SetBinContent(63,0.9940273);
trigeff->SetBinContent(64,0.9952786);
trigeff->SetBinContent(65,0.9977998);
trigeff->SetBinContent(66,0.9972015);
trigeff->SetBinContent(67,0.9908467);
trigeff->SetBinContent(68,1);
trigeff->SetBinContent(69,0.9988636);
trigeff->SetBinContent(70,0.9933155);
trigeff->SetBinContent(71,0.9986979);
trigeff->SetBinContent(72,0.9969651);
trigeff->SetBinContent(73,0.9971223);
trigeff->SetBinContent(74,0.9967213);
trigeff->SetBinContent(75,0.9982111);
trigeff->SetBinContent(76,0.9968504);
trigeff->SetBinContent(77,0.9967374);
trigeff->SetBinContent(78,0.9980237);
trigeff->SetBinContent(79,0.9977974);
trigeff->SetBinContent(80,1);
trigeff->SetBinContent(81,0.997669);
trigeff->SetBinContent(82,0.9977876);
trigeff->SetBinContent(83,0.9921671);
trigeff->SetBinContent(84,0.9942857);
trigeff->SetBinContent(85,1);
trigeff->SetBinContent(86,1);
trigeff->SetBinContent(87,0.9967532);
trigeff->SetBinContent(88,1);
trigeff->SetBinContent(89,0.9967427);
trigeff->SetBinContent(90,0.9961832);
trigeff->SetBinContent(91,1);
trigeff->SetBinContent(92,1);
trigeff->SetBinContent(93,1);
trigeff->SetBinContent(94,1);
trigeff->SetBinContent(95,1);
trigeff->SetBinContent(96,0.9940828);
trigeff->SetBinContent(97,1);
trigeff->SetBinContent(98,0.9927536);
trigeff->SetBinContent(99,1);
trigeff->SetBinContent(100,0.9896907);
trigeff->SetBinContent(101,0.990566);
trigeff->SetBinContent(102,0.988764);
trigeff->SetBinContent(103,1);
trigeff->SetBinContent(104,1);
trigeff->SetBinContent(105,1);
trigeff->SetBinContent(106,1);
trigeff->SetBinContent(107,1);
trigeff->SetBinContent(108,1);
trigeff->SetBinContent(109,1);
trigeff->SetBinContent(110,1);
trigeff->SetBinContent(111,1);
trigeff->SetBinContent(112,1);
trigeff->SetBinContent(113,1);
trigeff->SetBinContent(114,1);
trigeff->SetBinContent(115,1);
trigeff->SetBinContent(116,1);
trigeff->SetBinContent(117,1);
trigeff->SetBinContent(118,1);
trigeff->SetBinContent(119,1);
trigeff->SetBinContent(120,1);
trigeff->SetBinContent(121,1);
trigeff->SetBinContent(122,1);
trigeff->SetBinContent(123,1);
trigeff->SetBinContent(124,1);
trigeff->SetBinContent(125,1);
trigeff->SetBinContent(126,1);
trigeff->SetBinContent(127,1);
trigeff->SetBinContent(128,1);
trigeff->SetBinContent(129,1);
trigeff->SetBinContent(130,1);
trigeff->SetBinContent(131,1);
trigeff->SetBinContent(132,1);
trigeff->SetBinContent(133,1);
trigeff->SetBinContent(134,1);
trigeff->SetBinContent(135,1);
trigeff->SetBinContent(136,1);
trigeff->SetBinContent(137,1);
trigeff->SetBinContent(138,1);
trigeff->SetBinContent(139,1);
trigeff->SetBinContent(140,1);
trigeff->SetBinContent(141,1);
trigeff->SetBinContent(142,1);
trigeff->SetBinContent(143,1);
trigeff->SetBinContent(144,1);
trigeff->SetBinContent(145,1);
trigeff->SetBinContent(146,1);
trigeff->SetBinContent(147,1);
trigeff->SetBinContent(149,1);
trigeff->SetBinContent(151,1);
trigeff->SetBinContent(152,1);
trigeff->SetBinContent(154,1);
trigeff->SetBinContent(155,1);
trigeff->SetBinContent(157,1);
trigeff->SetBinContent(170,1);
trigeff->SetBinContent(173,1);
trigeff->SetBinContent(178,1);
trigeff->SetBinError(1,0.001231687);
trigeff->SetBinError(2,0.009082424);
trigeff->SetBinError(3,0.009308539);
trigeff->SetBinError(4,0.007834723);
trigeff->SetBinError(5,0.00598166);
trigeff->SetBinError(6,0.004397977);
trigeff->SetBinError(7,0.003418236);
trigeff->SetBinError(8,0.002747754);
trigeff->SetBinError(9,0.002375778);
trigeff->SetBinError(10,0.002196995);
trigeff->SetBinError(11,0.002011941);
trigeff->SetBinError(12,0.001831367);
trigeff->SetBinError(13,0.001784029);
trigeff->SetBinError(14,0.001916823);
trigeff->SetBinError(15,0.001817566);
trigeff->SetBinError(16,0.001819683);
trigeff->SetBinError(17,0.001818542);
trigeff->SetBinError(18,0.001890354);
trigeff->SetBinError(19,0.00177814);
trigeff->SetBinError(20,0.001864436);
trigeff->SetBinError(21,0.001663154);
trigeff->SetBinError(22,0.00170007);
trigeff->SetBinError(23,0.00173455);
trigeff->SetBinError(24,0.001481362);
trigeff->SetBinError(25,0.001655051);
trigeff->SetBinError(26,0.001603857);
trigeff->SetBinError(27,0.00195356);
trigeff->SetBinError(28,0.00163551);
trigeff->SetBinError(29,0.001474872);
trigeff->SetBinError(30,0.001623358);
trigeff->SetBinError(31,0.00154461);
trigeff->SetBinError(32,0.001534909);
trigeff->SetBinError(33,0.00158181);
trigeff->SetBinError(34,0.001491233);
trigeff->SetBinError(35,0.001286982);
trigeff->SetBinError(36,0.001328752);
trigeff->SetBinError(37,0.00152958);
trigeff->SetBinError(38,0.001348444);
trigeff->SetBinError(39,0.00166985);
trigeff->SetBinError(40,0.001624111);
trigeff->SetBinError(41,0.001586005);
trigeff->SetBinError(42,0.001827581);
trigeff->SetBinError(43,0.001420028);
trigeff->SetBinError(44,0.002227847);
trigeff->SetBinError(45,0.001874108);
trigeff->SetBinError(46,0.00114964);
trigeff->SetBinError(47,0.001293933);
trigeff->SetBinError(48,0.001462021);
trigeff->SetBinError(49,0.001515213);
trigeff->SetBinError(50,0.001377697);
trigeff->SetBinError(51,0.0009272125);
trigeff->SetBinError(52,0.001331017);
trigeff->SetBinError(53,0.001174711);
trigeff->SetBinError(54,0.001150746);
trigeff->SetBinError(55,0.001991336);
trigeff->SetBinError(56,0.001625642);
trigeff->SetBinError(57,0.001748087);
trigeff->SetBinError(58,0.00259828);
trigeff->SetBinError(59,0.0009945236);
trigeff->SetBinError(60,0.0009766661);
trigeff->SetBinError(61,0.000839278);
trigeff->SetBinError(62,0.001137741);
trigeff->SetBinError(63,0.002250715);
trigeff->SetBinError(64,0.0021065);
trigeff->SetBinError(65,0.001554078);
trigeff->SetBinError(66,0.001613457);
trigeff->SetBinError(67,0.003221342);
trigeff->SetBinError(69,0.001135718);
trigeff->SetBinError(70,0.002979388);
trigeff->SetBinError(71,0.001301235);
trigeff->SetBinError(72,0.00214274);
trigeff->SetBinError(73,0.00203191);
trigeff->SetBinError(74,0.002314579);
trigeff->SetBinError(75,0.001787308);
trigeff->SetBinError(76,0.002223598);
trigeff->SetBinError(77,0.00230327);
trigeff->SetBinError(78,0.001974331);
trigeff->SetBinError(79,0.002200216);
trigeff->SetBinError(81,0.002328284);
trigeff->SetBinError(82,0.002209941);
trigeff->SetBinError(83,0.00450458);
trigeff->SetBinError(84,0.004029049);
trigeff->SetBinError(87,0.003241478);
trigeff->SetBinError(89,0.00325202);
trigeff->SetBinError(90,0.003809503);
trigeff->SetBinError(96,0.005899627);
trigeff->SetBinError(98,0.007220074);
trigeff->SetBinError(100,0.010256);
trigeff->SetBinError(101,0.009389357);
trigeff->SetBinError(102,0.01117265);
trigeff->SetEntries(369636);
trigeff->SetDirectory(0);
trigeff->SetFillColor(1);
trigeff->SetFillStyle(0);
trigeff->SetLineStyle(2);
trigeff->SetLineWidth(3);
trigeff->SetMarkerStyle(20);
trigeff->GetXaxis()->SetTitle("M");
trigeff->GetXaxis()->SetLabelFont(42);
trigeff->GetXaxis()->SetLabelOffset(0.01);
trigeff->GetXaxis()->SetLabelSize(0.045);
trigeff->GetXaxis()->SetTitleSize(0.055);
trigeff->GetXaxis()->SetTitleFont(42);
trigeff->GetYaxis()->SetTitle("Selection efficiency");
trigeff->GetYaxis()->SetLabelFont(42);
trigeff->GetYaxis()->SetLabelOffset(0.01);
trigeff->GetYaxis()->SetLabelSize(0.045);
trigeff->GetYaxis()->SetTitleSize(0.055);
trigeff->GetYaxis()->SetTitleOffset(1.6);
trigeff->GetYaxis()->SetTitleFont(42);
trigeff->GetZaxis()->SetLabelFont(42);
trigeff->GetZaxis()->SetLabelSize(0.045);
trigeff->GetZaxis()->SetTitleFont(42);
trigeff->Draw("p same");
TLegend *leg = new TLegend(0.57,0.22,1,0.42,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(62);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(19);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("trigeff","PYTHIA 7 TeV","l");
entry->SetLineColor(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
//entry=leg->AddEntry("trigeff","(Atlas tune)","");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry=leg->AddEntry("trigeff","PHOJET 7 TeV","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
leg->Draw();
printFinalCanvases(cEvtSel,"trigeff");
/*
TLatex * tex = new TLatex(0.85,0.9,"CMS");
tex->SetNDC();
tex->SetTextSize(0.04);
tex->SetLineWidth(2);
tex->Draw();
cEvtSel->Modified();
cEvtSel->cd();
cEvtSel->SetSelected(cEvtSel);
*/
}
void diffrac(){
gROOT->Reset();
//gROOT->ProcessLine(".x rootlogon.C");
//=========Macro generated from canvas: cEvtSel/EvtSel
//========= (Wed Apr 14 17:19:02 2010) by ROOT version5.22/00d
TCanvas *cEvtSel = new TCanvas("cEvtSel", "EvtSel",550,600);
/*
TCanvas *cEvtSel = new TCanvas("cEvtSel", "EvtSel",0,22,550,600);
gStyle->SetOptFit(1);
gStyle->SetOptStat(0);
cEvtSel->Range(-32.99367,-0.009375,158.1456,0.053125);
cEvtSel->SetFillColor(0);
cEvtSel->SetBorderMode(0);
cEvtSel->SetBorderSize(0);
cEvtSel->SetTickx(1);
cEvtSel->SetTicky(1);
cEvtSel->SetLeftMargin(0.17);
cEvtSel->SetRightMargin(0.04);
cEvtSel->SetTopMargin(0.05);
cEvtSel->SetBottomMargin(0.15);
cEvtSel->SetFrameLineColor(0);
cEvtSel->SetFrameBorderMode(0);
cEvtSel->SetFrameLineColor(0);
cEvtSel->SetFrameBorderMode(0);
*/
//TH1 *hFrame2 = new TH2D("hFrame2","",1,-0.5,60.5,1,0,0.06);
TH1 *hFrame2 = new TH2D("hFrame2","",1,-0.5,150.5,1,0,0.05);
hFrame2->SetDirectory(0);
hFrame2->SetStats(0);
hFrame2->GetXaxis()->SetTitle("Charged-particle multiplicity");
hFrame2->GetXaxis()->CenterTitle(true);
hFrame2->GetYaxis()->SetTitle("Fraction of events");
hFrame2->GetYaxis()->CenterTitle(true);
hFrame2->GetYaxis()->SetTitleOffset(1.5);
hFrame2->Draw("");
hFrame2->GetXaxis()->SetNdivisions(312);
hFrame2->GetYaxis()->SetNdivisions(305);
/*
TH1 *hFrame2 = new TH2D("hFrame2","",1,-0.5,150.5,1,0,0.05);
hFrame2->SetDirectory(0);
hFrame2->SetStats(0);
hFrame2->SetFillColor(1);
hFrame2->SetFillStyle(0);
hFrame2->SetLineStyle(0);
hFrame2->SetMarkerStyle(20);
hFrame2->SetMarkerSize(1.5);
hFrame2->GetXaxis()->SetTitle("Charged-particle multiplicity");
hFrame2->GetXaxis()->CenterTitle(true);
hFrame2->GetXaxis()->SetNdivisions(312);
hFrame2->GetXaxis()->SetLabelFont(42);
hFrame2->GetXaxis()->SetLabelOffset(0.01);
hFrame2->GetXaxis()->SetLabelSize(0.045);
hFrame2->GetXaxis()->SetTitleSize(0.055);
hFrame2->GetXaxis()->SetTitleFont(42);
hFrame2->GetYaxis()->SetTitle("Fraction of events");
hFrame2->GetYaxis()->CenterTitle(true);
hFrame2->GetYaxis()->SetLabelFont(42);
hFrame2->GetYaxis()->SetLabelOffset(0.01);
hFrame2->GetYaxis()->SetLabelSize(0.045);
hFrame2->GetYaxis()->SetTitleSize(0.055);
hFrame2->GetYaxis()->SetTitleOffset(1.5);
hFrame2->GetYaxis()->SetTitleFont(42);
hFrame2->GetZaxis()->SetLabelFont(42);
hFrame2->GetZaxis()->SetLabelSize(0.045);
hFrame2->GetZaxis()->SetTitleFont(42);
hFrame2->Draw("");
*/
TH1 *diffrac = new TH1D("diffrac","diffrac",200,0,200);
diffrac->SetBinContent(1,0.0004697663);
diffrac->SetBinContent(2,0.008010118);
diffrac->SetBinContent(3,0.0127921);
diffrac->SetBinContent(4,0.01698988);
diffrac->SetBinContent(5,0.02012166);
diffrac->SetBinContent(6,0.022341);
diffrac->SetBinContent(7,0.02315105);
diffrac->SetBinContent(8,0.02518068);
diffrac->SetBinContent(9,0.02668032);
diffrac->SetBinContent(10,0.02875512);
diffrac->SetBinContent(11,0.03089617);
diffrac->SetBinContent(12,0.0324741);
diffrac->SetBinContent(13,0.03311853);
diffrac->SetBinContent(14,0.03265478);
diffrac->SetBinContent(15,0.03213081);
diffrac->SetBinContent(16,0.03054987);
diffrac->SetBinContent(17,0.0283245);
diffrac->SetBinContent(18,0.02640026);
diffrac->SetBinContent(19,0.02394001);
diffrac->SetBinContent(20,0.02234401);
diffrac->SetBinContent(21,0.02026921);
diffrac->SetBinContent(22,0.01865213);
diffrac->SetBinContent(23,0.01720369);
diffrac->SetBinContent(24,0.01639665);
diffrac->SetBinContent(25,0.01546013);
diffrac->SetBinContent(26,0.01423151);
diffrac->SetBinContent(27,0.01359311);
diffrac->SetBinContent(28,0.01247892);
diffrac->SetBinContent(29,0.01204529);
diffrac->SetBinContent(30,0.01130752);
diffrac->SetBinContent(31,0.01089798);
diffrac->SetBinContent(32,0.01070826);
diffrac->SetBinContent(33,0.009497711);
diffrac->SetBinContent(34,0.009579017);
diffrac->SetBinContent(35,0.008805107);
diffrac->SetBinContent(36,0.008741869);
diffrac->SetBinContent(37,0.008522043);
diffrac->SetBinContent(38,0.00790171);
diffrac->SetBinContent(39,0.007799325);
diffrac->SetBinContent(40,0.007486148);
diffrac->SetBinContent(41,0.006838714);
diffrac->SetBinContent(42,0.006947121);
diffrac->SetBinContent(43,0.006504457);
diffrac->SetBinContent(44,0.0063298);
diffrac->SetBinContent(45,0.006122019);
diffrac->SetBinContent(46,0.006061792);
diffrac->SetBinContent(47,0.005823898);
diffrac->SetBinContent(48,0.005661286);
diffrac->SetBinContent(49,0.005543845);
diffrac->SetBinContent(50,0.005363166);
diffrac->SetBinContent(51,0.005215611);
diffrac->SetBinContent(52,0.00503192);
diffrac->SetBinContent(53,0.004800048);
diffrac->SetBinContent(54,0.004839195);
diffrac->SetBinContent(55,0.004583233);
diffrac->SetBinContent(56,0.004544086);
diffrac->SetBinContent(57,0.004471814);
diffrac->SetBinContent(58,0.004146591);
diffrac->SetBinContent(59,0.004267044);
diffrac->SetBinContent(60,0.003990002);
diffrac->SetBinContent(61,0.004083353);
diffrac->SetBinContent(62,0.003806312);
diffrac->SetBinContent(63,0.003782221);
diffrac->SetBinContent(64,0.003812334);
diffrac->SetBinContent(65,0.003472055);
diffrac->SetBinContent(66,0.00366478);
diffrac->SetBinContent(67,0.003195013);
diffrac->SetBinContent(68,0.003285353);
diffrac->SetBinContent(69,0.00323416);
diffrac->SetBinContent(70,0.0031649);
diffrac->SetBinContent(71,0.003110696);
diffrac->SetBinContent(72,0.003170923);
diffrac->SetBinContent(73,0.002966153);
diffrac->SetBinContent(74,0.002957119);
diffrac->SetBinContent(75,0.002905926);
diffrac->SetBinContent(76,0.002954107);
diffrac->SetBinContent(77,0.002785473);
diffrac->SetBinContent(78,0.002836666);
diffrac->SetBinContent(79,0.002725247);
diffrac->SetBinContent(80,0.002728258);
diffrac->SetBinContent(81,0.002568658);
diffrac->SetBinContent(82,0.002288605);
diffrac->SetBinContent(83,0.002499398);
diffrac->SetBinContent(84,0.002445194);
diffrac->SetBinContent(85,0.002378945);
diffrac->SetBinContent(86,0.00227656);
diffrac->SetBinContent(87,0.002418092);
diffrac->SetBinContent(88,0.002351843);
diffrac->SetBinContent(89,0.002381956);
diffrac->SetBinContent(90,0.002189232);
diffrac->SetBinContent(91,0.002222356);
diffrac->SetBinContent(92,0.002198265);
diffrac->SetBinContent(93,0.002219345);
diffrac->SetBinContent(94,0.002198265);
diffrac->SetBinContent(95,0.002089858);
diffrac->SetBinContent(96,0.001996507);
diffrac->SetBinContent(97,0.001972416);
diffrac->SetBinContent(98,0.001903156);
diffrac->SetBinContent(99,0.001818839);
diffrac->SetBinContent(100,0.001876054);
diffrac->SetBinContent(101,0.001864009);
diffrac->SetBinContent(102,0.001767646);
diffrac->SetBinContent(103,0.001879065);
diffrac->SetBinContent(104,0.001894122);
diffrac->SetBinContent(105,0.001695375);
diffrac->SetBinContent(106,0.001674295);
diffrac->SetBinContent(107,0.001611058);
diffrac->SetBinContent(108,0.001556854);
diffrac->SetBinContent(109,0.001620092);
diffrac->SetBinContent(110,0.001556854);
diffrac->SetBinContent(111,0.001568899);
diffrac->SetBinContent(112,0.001475548);
diffrac->SetBinContent(113,0.001376174);
diffrac->SetBinContent(114,0.001448446);
diffrac->SetBinContent(115,0.001520718);
diffrac->SetBinContent(116,0.001261744);
diffrac->SetBinContent(117,0.001448446);
diffrac->SetBinContent(118,0.001282823);
diffrac->SetBinContent(119,0.001219586);
diffrac->SetBinContent(120,0.001201518);
diffrac->SetBinContent(121,0.001237654);
diffrac->SetBinContent(122,0.001141291);
diffrac->SetBinContent(123,0.001195495);
diffrac->SetBinContent(124,0.001150325);
diffrac->SetBinContent(125,0.001195495);
diffrac->SetBinContent(126,0.0009816912);
diffrac->SetBinContent(127,0.001135269);
diffrac->SetBinContent(128,0.00100277);
diffrac->SetBinContent(129,0.001005782);
diffrac->SetBinContent(130,0.0009003854);
diffrac->SetBinContent(131,0.0008823175);
diffrac->SetBinContent(132,0.0009576006);
diffrac->SetBinContent(133,0.0008732835);
diffrac->SetBinContent(134,0.0008220911);
diffrac->SetBinContent(135,0.000831125);
diffrac->SetBinContent(136,0.0007437967);
diffrac->SetBinContent(137,0.000662491);
diffrac->SetBinContent(138,0.0007046495);
diffrac->SetBinContent(139,0.0007287401);
diffrac->SetBinContent(140,0.0007407854);
diffrac->SetBinContent(141,0.0006263551);
diffrac->SetBinContent(142,0.0006745363);
diffrac->SetBinContent(143,0.000653457);
diffrac->SetBinContent(144,0.0006082872);
diffrac->SetBinContent(145,0.0005420381);
diffrac->SetBinContent(146,0.0005631173);
diffrac->SetBinContent(147,0.0005059022);
diffrac->SetBinContent(148,0.0005781739);
diffrac->SetBinContent(149,0.0005239701);
diffrac->SetBinContent(150,0.0004547097);
diffrac->SetBinContent(151,0.0004697663);
diffrac->SetBinContent(152,0.0004547097);
diffrac->SetBinContent(153,0.0004396531);
diffrac->SetBinContent(154,0.0004697663);
diffrac->SetBinContent(155,0.0003884606);
diffrac->SetBinContent(156,0.0003613587);
diffrac->SetBinContent(157,0.0003523247);
diffrac->SetBinContent(158,0.0003703927);
diffrac->SetBinContent(159,0.0003794266);
diffrac->SetBinContent(160,0.0003041436);
diffrac->SetBinContent(161,0.0003192002);
diffrac->SetBinContent(162,0.000280053);
diffrac->SetBinContent(163,0.0002770417);
diffrac->SetBinContent(164,0.0002228379);
diffrac->SetBinContent(165,0.0002529511);
diffrac->SetBinContent(166,0.0002469285);
diffrac->SetBinContent(167,0.0002318718);
diffrac->SetBinContent(168,0.0003041436);
diffrac->SetBinContent(169,0.0002198265);
diffrac->SetBinContent(170,0.0002499398);
diffrac->SetBinContent(171,0.0001565888);
diffrac->SetBinContent(172,0.0001656227);
diffrac->SetBinContent(173,0.0001535775);
diffrac->SetBinContent(174,0.0001626114);
diffrac->SetBinContent(175,0.0001987473);
diffrac->SetBinContent(176,0.0001445435);
diffrac->SetBinContent(177,0.0001174416);
diffrac->SetBinContent(178,0.0001445435);
diffrac->SetBinContent(179,0.0001415322);
diffrac->SetBinContent(180,0.0001445435);
diffrac->SetBinContent(181,0.0001264755);
diffrac->SetBinContent(182,6.62491e-05);
diffrac->SetBinContent(183,7.227174e-05);
diffrac->SetBinContent(184,0.0001053963);
diffrac->SetBinContent(185,9.033968e-05);
diffrac->SetBinContent(186,8.732835e-05);
diffrac->SetBinContent(187,6.022645e-05);
diffrac->SetBinContent(188,9.3351e-05);
diffrac->SetBinContent(189,8.732835e-05);
diffrac->SetBinContent(190,8.431703e-05);
diffrac->SetBinContent(191,9.033968e-05);
diffrac->SetBinContent(192,4.215852e-05);
diffrac->SetBinContent(193,6.62491e-05);
diffrac->SetBinContent(194,7.528306e-05);
diffrac->SetBinContent(195,4.215852e-05);
diffrac->SetBinContent(196,4.516984e-05);
diffrac->SetBinContent(197,3.613587e-05);
diffrac->SetBinContent(198,6.022645e-05);
diffrac->SetBinContent(199,4.516984e-05);
diffrac->SetBinContent(200,3.011323e-05);
diffrac->SetBinContent(201,0.000457721);
diffrac->SetBinError(1,3.761141e-05);
diffrac->SetBinError(2,0.0001553095);
diffrac->SetBinError(3,0.000196268);
diffrac->SetBinError(4,0.0002261902);
diffrac->SetBinError(5,0.0002461561);
diffrac->SetBinError(6,0.0002593761);
diffrac->SetBinError(7,0.0002640365);
diffrac->SetBinError(8,0.0002753673);
diffrac->SetBinError(9,0.0002834485);
diffrac->SetBinError(10,0.0002942634);
diffrac->SetBinError(11,0.0003050219);
diffrac->SetBinError(12,0.0003127139);
diffrac->SetBinError(13,0.0003158015);
diffrac->SetBinError(14,0.0003135827);
diffrac->SetBinError(15,0.0003110566);
diffrac->SetBinError(16,0.0003033076);
diffrac->SetBinError(17,0.0002920517);
diffrac->SetBinError(18,0.0002819569);
diffrac->SetBinError(19,0.0002684979);
diffrac->SetBinError(20,0.0002593936);
diffrac->SetBinError(21,0.0002470569);
diffrac->SetBinError(22,0.000236997);
diffrac->SetBinError(23,0.000227609);
diffrac->SetBinError(24,0.0002222062);
diffrac->SetBinError(25,0.0002157671);
diffrac->SetBinError(26,0.0002070161);
diffrac->SetBinError(27,0.0002023196);
diffrac->SetBinError(28,0.0001938506);
diffrac->SetBinError(29,0.0001904528);
diffrac->SetBinError(30,0.000184528);
diffrac->SetBinError(31,0.0001811555);
diffrac->SetBinError(32,0.0001795718);
diffrac->SetBinError(33,0.0001691173);
diffrac->SetBinError(34,0.0001698397);
diffrac->SetBinError(35,0.0001628343);
diffrac->SetBinError(36,0.0001622485);
diffrac->SetBinError(37,0.0001601956);
diffrac->SetBinError(38,0.000154255);
diffrac->SetBinError(39,0.0001532524);
diffrac->SetBinError(40,0.000150144);
diffrac->SetBinError(41,0.0001435046);
diffrac->SetBinError(42,0.0001446376);
diffrac->SetBinError(43,0.0001399536);
diffrac->SetBinError(44,0.0001380618);
diffrac->SetBinError(45,0.0001357769);
diffrac->SetBinError(46,0.0001351074);
diffrac->SetBinError(47,0.0001324297);
diffrac->SetBinError(48,0.0001305678);
diffrac->SetBinError(49,0.0001292064);
diffrac->SetBinError(50,0.0001270835);
diffrac->SetBinError(51,0.0001253231);
diffrac->SetBinError(52,0.0001230964);
diffrac->SetBinError(53,0.0001202268);
diffrac->SetBinError(54,0.0001207161);
diffrac->SetBinError(55,0.0001174802);
diffrac->SetBinError(56,0.0001169774);
diffrac->SetBinError(57,0.0001160434);
diffrac->SetBinError(58,0.000111744);
diffrac->SetBinError(59,0.0001133554);
diffrac->SetBinError(60,0.0001096138);
diffrac->SetBinError(61,0.0001108887);
diffrac->SetBinError(62,0.0001070609);
diffrac->SetBinError(63,0.0001067215);
diffrac->SetBinError(64,0.0001071455);
diffrac->SetBinError(65,0.000102252);
diffrac->SetBinError(66,0.0001050516);
diffrac->SetBinError(67,9.808779e-05);
diffrac->SetBinError(68,9.946486e-05);
diffrac->SetBinError(69,9.868688e-05);
diffrac->SetBinError(70,9.762446e-05);
diffrac->SetBinError(71,9.678486e-05);
diffrac->SetBinError(72,9.77173e-05);
diffrac->SetBinError(73,9.450948e-05);
diffrac->SetBinError(74,9.436545e-05);
diffrac->SetBinError(75,9.354508e-05);
diffrac->SetBinError(76,9.431739e-05);
diffrac->SetBinError(77,9.15858e-05);
diffrac->SetBinError(78,9.242357e-05);
diffrac->SetBinError(79,9.059027e-05);
diffrac->SetBinError(80,9.064031e-05);
diffrac->SetBinError(81,8.794918e-05);
diffrac->SetBinError(82,8.301643e-05);
diffrac->SetBinError(83,8.675536e-05);
diffrac->SetBinError(84,8.580948e-05);
diffrac->SetBinError(85,8.463906e-05);
diffrac->SetBinError(86,8.279768e-05);
diffrac->SetBinError(87,8.533261e-05);
diffrac->SetBinError(88,8.415556e-05);
diffrac->SetBinError(89,8.469261e-05);
diffrac->SetBinError(90,8.11941e-05);
diffrac->SetBinError(91,8.180606e-05);
diffrac->SetBinError(92,8.136146e-05);
diffrac->SetBinError(93,8.175061e-05);
diffrac->SetBinError(94,8.136146e-05);
diffrac->SetBinError(95,7.932992e-05);
diffrac->SetBinError(96,7.75379e-05);
diffrac->SetBinError(97,7.706868e-05);
diffrac->SetBinError(98,7.570348e-05);
diffrac->SetBinError(99,7.40075e-05);
diffrac->SetBinError(100,7.516251e-05);
diffrac->SetBinError(101,7.492083e-05);
diffrac->SetBinError(102,7.295857e-05);
diffrac->SetBinError(103,7.522281e-05);
diffrac->SetBinError(104,7.552359e-05);
diffrac->SetBinError(105,7.145152e-05);
diffrac->SetBinError(106,7.100594e-05);
diffrac->SetBinError(107,6.965209e-05);
diffrac->SetBinError(108,6.847035e-05);
diffrac->SetBinError(109,6.984711e-05);
diffrac->SetBinError(110,6.847035e-05);
diffrac->SetBinError(111,6.873472e-05);
diffrac->SetBinError(112,6.665847e-05);
diffrac->SetBinError(113,6.437472e-05);
diffrac->SetBinError(114,6.604346e-05);
diffrac->SetBinError(115,6.767106e-05);
diffrac->SetBinError(116,6.164024e-05);
diffrac->SetBinError(117,6.604346e-05);
diffrac->SetBinError(118,6.2153e-05);
diffrac->SetBinError(119,6.06017e-05);
diffrac->SetBinError(120,6.015112e-05);
diffrac->SetBinError(121,6.104895e-05);
diffrac->SetBinError(122,5.862419e-05);
diffrac->SetBinError(123,6.000018e-05);
diffrac->SetBinError(124,5.885576e-05);
diffrac->SetBinError(125,6.000018e-05);
diffrac->SetBinError(126,5.437084e-05);
diffrac->SetBinError(127,5.846931e-05);
diffrac->SetBinError(128,5.495148e-05);
diffrac->SetBinError(129,5.503393e-05);
diffrac->SetBinError(130,5.207063e-05);
diffrac->SetBinError(131,5.154554e-05);
diffrac->SetBinError(132,5.369957e-05);
diffrac->SetBinError(133,5.128098e-05);
diffrac->SetBinError(134,4.975521e-05);
diffrac->SetBinError(135,5.002785e-05);
diffrac->SetBinError(136,4.732665e-05);
diffrac->SetBinError(137,4.466513e-05);
diffrac->SetBinError(138,4.606438e-05);
diffrac->SetBinError(139,4.684519e-05);
diffrac->SetBinError(140,4.723075e-05);
diffrac->SetBinError(141,4.342991e-05);
diffrac->SetBinError(142,4.506935e-05);
diffrac->SetBinError(143,4.435955e-05);
diffrac->SetBinError(144,4.279894e-05);
diffrac->SetBinError(145,4.040113e-05);
diffrac->SetBinError(146,4.117922e-05);
diffrac->SetBinError(147,3.90312e-05);
diffrac->SetBinError(148,4.172611e-05);
diffrac->SetBinError(149,3.972207e-05);
diffrac->SetBinError(150,3.700375e-05);
diffrac->SetBinError(151,3.761141e-05);
diffrac->SetBinError(152,3.700375e-05);
diffrac->SetBinError(153,3.638595e-05);
diffrac->SetBinError(154,3.761141e-05);
diffrac->SetBinError(155,3.420205e-05);
diffrac->SetBinError(156,3.298739e-05);
diffrac->SetBinError(157,3.257243e-05);
diffrac->SetBinError(158,3.339718e-05);
diffrac->SetBinError(159,3.380201e-05);
diffrac->SetBinError(160,3.026342e-05);
diffrac->SetBinError(161,3.100346e-05);
diffrac->SetBinError(162,2.904014e-05);
diffrac->SetBinError(163,2.888359e-05);
diffrac->SetBinError(164,2.590438e-05);
diffrac->SetBinError(165,2.759923e-05);
diffrac->SetBinError(166,2.726869e-05);
diffrac->SetBinError(167,2.642425e-05);
diffrac->SetBinError(168,3.026342e-05);
diffrac->SetBinError(169,2.572875e-05);
diffrac->SetBinError(170,2.743445e-05);
diffrac->SetBinError(171,2.171496e-05);
diffrac->SetBinError(172,2.233257e-05);
diffrac->SetBinError(173,2.150514e-05);
diffrac->SetBinError(174,2.212861e-05);
diffrac->SetBinError(175,2.44641e-05);
diffrac->SetBinError(176,2.086305e-05);
diffrac->SetBinError(177,1.88057e-05);
diffrac->SetBinError(178,2.086305e-05);
diffrac->SetBinError(179,2.064459e-05);
diffrac->SetBinError(180,2.086305e-05);
diffrac->SetBinError(181,1.95156e-05);
diffrac->SetBinError(182,1.412435e-05);
diffrac->SetBinError(183,1.475241e-05);
diffrac->SetBinError(184,1.781522e-05);
diffrac->SetBinError(185,1.649369e-05);
diffrac->SetBinError(186,1.621647e-05);
diffrac->SetBinError(187,1.346704e-05);
diffrac->SetBinError(188,1.676633e-05);
diffrac->SetBinError(189,1.621647e-05);
diffrac->SetBinError(190,1.593442e-05);
diffrac->SetBinError(191,1.649369e-05);
diffrac->SetBinError(192,1.126734e-05);
diffrac->SetBinError(193,1.412435e-05);
diffrac->SetBinError(194,1.505661e-05);
diffrac->SetBinError(195,1.126734e-05);
diffrac->SetBinError(196,1.16628e-05);
diffrac->SetBinError(197,1.043153e-05);
diffrac->SetBinError(198,1.346704e-05);
diffrac->SetBinError(199,1.16628e-05);
diffrac->SetBinError(200,9.522638e-06);
diffrac->SetBinError(201,3.712608e-05);
diffrac->SetEntries(332232);
diffrac->SetDirectory(0);
diffrac->SetFillColor(1);
diffrac->SetFillStyle(0);
diffrac->SetLineStyle(0);
diffrac->SetLineWidth(2);
diffrac->SetMarkerStyle(20);
diffrac->SetMarkerSize(1.5);
diffrac->GetXaxis()->SetTitle("M");
diffrac->GetXaxis()->SetLabelFont(42);
diffrac->GetXaxis()->SetLabelOffset(0.01);
diffrac->GetXaxis()->SetLabelSize(0.045);
diffrac->GetXaxis()->SetTitleSize(0.055);
diffrac->GetXaxis()->SetTitleFont(42);
diffrac->GetYaxis()->SetTitle("Fraction of events");
diffrac->GetYaxis()->SetLabelFont(42);
diffrac->GetYaxis()->SetLabelOffset(0.01);
diffrac->GetYaxis()->SetLabelSize(0.045);
diffrac->GetYaxis()->SetTitleSize(0.055);
diffrac->GetYaxis()->SetTitleOffset(1.6);
diffrac->GetYaxis()->SetTitleFont(42);
diffrac->GetZaxis()->SetLabelFont(42);
diffrac->GetZaxis()->SetLabelSize(0.045);
diffrac->GetZaxis()->SetTitleFont(42);
diffrac->Draw("hist same");
TH1 *diffrac = new TH1D("diffrac","diffrac",200,0,200);
diffrac->SetBinContent(1,0.0001535956);
diffrac->SetBinContent(2,0.002138746);
diffrac->SetBinContent(3,0.003918136);
diffrac->SetBinContent(4,0.006984252);
diffrac->SetBinContent(5,0.01149359);
diffrac->SetBinContent(6,0.01725487);
diffrac->SetBinContent(7,0.02335522);
diffrac->SetBinContent(8,0.02994824);
diffrac->SetBinContent(9,0.03493575);
diffrac->SetBinContent(10,0.03668616);
diffrac->SetBinContent(11,0.03929439);
diffrac->SetBinContent(12,0.03968852);
diffrac->SetBinContent(13,0.03844526);
diffrac->SetBinContent(14,0.0360341);
diffrac->SetBinContent(15,0.03434745);
diffrac->SetBinContent(16,0.03306942);
diffrac->SetBinContent(17,0.03029311);
diffrac->SetBinContent(18,0.02782978);
diffrac->SetBinContent(19,0.02716903);
diffrac->SetBinContent(20,0.02540413);
diffrac->SetBinContent(21,0.02289154);
diffrac->SetBinContent(22,0.02272345);
diffrac->SetBinContent(23,0.02145412);
diffrac->SetBinContent(24,0.02181927);
diffrac->SetBinContent(25,0.02003698);
diffrac->SetBinContent(26,0.01887487);
diffrac->SetBinContent(27,0.01739977);
diffrac->SetBinContent(28,0.01641734);
diffrac->SetBinContent(29,0.01567544);
diffrac->SetBinContent(30,0.01542621);
diffrac->SetBinContent(31,0.01473069);
diffrac->SetBinContent(32,0.01400618);
diffrac->SetBinContent(33,0.01286146);
diffrac->SetBinContent(34,0.01269917);
diffrac->SetBinContent(35,0.0127919);
diffrac->SetBinContent(36,0.01181527);
diffrac->SetBinContent(37,0.01151387);
diffrac->SetBinContent(38,0.01042421);
diffrac->SetBinContent(39,0.01025613);
diffrac->SetBinContent(40,0.009320064);
diffrac->SetBinContent(41,0.009027363);
diffrac->SetBinContent(42,0.008412981);
diffrac->SetBinContent(43,0.00834053);
diffrac->SetBinContent(44,0.007647901);
diffrac->SetBinContent(45,0.007917418);
diffrac->SetBinContent(46,0.007523286);
diffrac->SetBinContent(47,0.006679959);
diffrac->SetBinContent(48,0.006526363);
diffrac->SetBinContent(49,0.006004718);
diffrac->SetBinContent(50,0.005914879);
diffrac->SetBinContent(51,0.006233662);
diffrac->SetBinContent(52,0.006123537);
diffrac->SetBinContent(53,0.005494665);
diffrac->SetBinContent(54,0.005019388);
diffrac->SetBinContent(55,0.004778851);
diffrac->SetBinContent(56,0.004686114);
diffrac->SetBinContent(57,0.004355739);
diffrac->SetBinContent(58,0.004112304);
diffrac->SetBinContent(59,0.004112304);
diffrac->SetBinContent(60,0.004187653);
diffrac->SetBinContent(61,0.003448656);
diffrac->SetBinContent(62,0.003593557);
diffrac->SetBinContent(63,0.003376205);
diffrac->SetBinContent(64,0.003054524);
diffrac->SetBinContent(65,0.002628513);
diffrac->SetBinContent(66,0.003097994);
diffrac->SetBinContent(67,0.002509694);
diffrac->SetBinContent(68,0.002495204);
diffrac->SetBinContent(69,0.002547368);
diffrac->SetBinContent(70,0.002153236);
diffrac->SetBinContent(71,0.002222789);
diffrac->SetBinContent(72,0.001904006);
diffrac->SetBinContent(73,0.002008335);
diffrac->SetBinContent(74,0.001762002);
diffrac->SetBinContent(75,0.001617101);
diffrac->SetBinContent(76,0.001834453);
diffrac->SetBinContent(77,0.001770696);
diffrac->SetBinContent(78,0.001463505);
diffrac->SetBinContent(79,0.001312808);
diffrac->SetBinContent(80,0.001193988);
diffrac->SetBinContent(81,0.001240357);
diffrac->SetBinContent(82,0.001307011);
diffrac->SetBinContent(83,0.001101251);
diffrac->SetBinContent(84,0.001008514);
diffrac->SetBinContent(85,0.0008838991);
diffrac->SetBinContent(86,0.0009969223);
diffrac->SetBinContent(87,0.0008896952);
diffrac->SetBinContent(88,0.001011412);
diffrac->SetBinContent(89,0.0008867972);
diffrac->SetBinContent(90,0.0007563858);
diffrac->SetBinContent(91,0.000681037);
diffrac->SetBinContent(92,0.0005911981);
diffrac->SetBinContent(93,0.0007360996);
diffrac->SetBinContent(94,0.0006259745);
diffrac->SetBinContent(95,0.0004955631);
diffrac->SetBinContent(96,0.000486869);
diffrac->SetBinContent(97,0.0004462966);
diffrac->SetBinContent(98,0.0003970301);
diffrac->SetBinContent(99,0.0003448656);
diffrac->SetBinContent(100,0.0002782109);
diffrac->SetBinContent(101,0.0003042931);
diffrac->SetBinContent(102,0.0002550266);
diffrac->SetBinContent(103,0.0003361715);
diffrac->SetBinContent(104,0.0001854739);
diffrac->SetBinContent(105,0.0003158853);
diffrac->SetBinContent(106,0.0002144542);
diffrac->SetBinContent(107,0.0001767798);
diffrac->SetBinContent(108,0.0001477995);
diffrac->SetBinContent(109,0.0001506976);
diffrac->SetBinContent(110,0.0001362074);
diffrac->SetBinContent(111,0.0002463325);
diffrac->SetBinContent(112,0.0001391054);
diffrac->SetBinContent(113,0.0002173522);
diffrac->SetBinContent(114,9.853302e-05);
diffrac->SetBinContent(115,0.0001420035);
diffrac->SetBinContent(116,0.0001159212);
diffrac->SetBinContent(117,8.983893e-05);
diffrac->SetBinContent(118,6.085863e-05);
diffrac->SetBinContent(119,6.375666e-05);
diffrac->SetBinContent(120,6.375666e-05);
diffrac->SetBinContent(121,5.79606e-05);
diffrac->SetBinContent(122,6.665469e-05);
diffrac->SetBinContent(123,4.926651e-05);
diffrac->SetBinContent(124,4.347045e-05);
diffrac->SetBinContent(125,4.636848e-05);
diffrac->SetBinContent(126,5.216454e-05);
diffrac->SetBinContent(127,3.477636e-05);
diffrac->SetBinContent(128,3.767439e-05);
diffrac->SetBinContent(129,2.318424e-05);
diffrac->SetBinContent(130,2.028621e-05);
diffrac->SetBinContent(131,3.187833e-05);
diffrac->SetBinContent(132,2.608227e-05);
diffrac->SetBinContent(133,3.477636e-05);
diffrac->SetBinContent(134,1.738818e-05);
diffrac->SetBinContent(135,1.738818e-05);
diffrac->SetBinContent(136,8.69409e-06);
diffrac->SetBinContent(137,1.449015e-05);
diffrac->SetBinContent(138,1.449015e-05);
diffrac->SetBinContent(139,2.89803e-06);
diffrac->SetBinContent(140,5.79606e-06);
diffrac->SetBinContent(141,8.69409e-06);
diffrac->SetBinContent(142,1.159212e-05);
diffrac->SetBinContent(143,2.89803e-06);
diffrac->SetBinContent(144,1.159212e-05);
diffrac->SetBinContent(145,2.89803e-06);
diffrac->SetBinContent(146,2.89803e-06);
diffrac->SetBinContent(147,5.79606e-06);
diffrac->SetBinContent(149,8.69409e-06);
diffrac->SetBinContent(151,5.79606e-06);
diffrac->SetBinContent(152,8.69409e-06);
diffrac->SetBinContent(154,2.89803e-06);
diffrac->SetBinContent(155,2.89803e-06);
diffrac->SetBinContent(157,5.79606e-06);
diffrac->SetBinContent(170,2.89803e-06);
diffrac->SetBinContent(173,5.79606e-06);
diffrac->SetBinContent(178,2.89803e-06);
diffrac->SetBinError(1,2.109798e-05);
diffrac->SetBinError(2,7.872833e-05);
diffrac->SetBinError(3,0.0001065593);
diffrac->SetBinError(4,0.0001422694);
diffrac->SetBinError(5,0.0001825069);
diffrac->SetBinError(6,0.0002236183);
diffrac->SetBinError(7,0.0002601617);
diffrac->SetBinError(8,0.000294603);
diffrac->SetBinError(9,0.00031819);
diffrac->SetBinError(10,0.0003260638);
diffrac->SetBinError(11,0.0003374556);
diffrac->SetBinError(12,0.0003391438);
diffrac->SetBinError(13,0.0003337896);
diffrac->SetBinError(14,0.0003231531);
diffrac->SetBinError(15,0.0003154995);
diffrac->SetBinError(16,0.0003095742);
diffrac->SetBinError(17,0.0002962943);
diffrac->SetBinError(18,0.0002839921);
diffrac->SetBinError(19,0.0002806005);
diffrac->SetBinError(20,0.0002713336);
diffrac->SetBinError(21,0.0002575662);
diffrac->SetBinError(22,0.0002566189);
diffrac->SetBinError(23,0.0002493485);
diffrac->SetBinError(24,0.0002514615);
diffrac->SetBinError(25,0.0002409725);
diffrac->SetBinError(26,0.0002338802);
diffrac->SetBinError(27,0.0002245552);
diffrac->SetBinError(28,0.0002181237);
diffrac->SetBinError(29,0.0002131382);
diffrac->SetBinError(30,0.0002114371);
diffrac->SetBinError(31,0.0002066155);
diffrac->SetBinError(32,0.0002014704);
diffrac->SetBinError(33,0.0001930619);
diffrac->SetBinError(34,0.0001918399);
diffrac->SetBinError(35,0.0001925391);
diffrac->SetBinError(36,0.0001850432);
diffrac->SetBinError(37,0.0001826679);
diffrac->SetBinError(38,0.0001738093);
diffrac->SetBinError(39,0.0001724023);
diffrac->SetBinError(40,0.0001643467);
diffrac->SetBinError(41,0.0001617454);
diffrac->SetBinError(42,0.0001561444);
diffrac->SetBinError(43,0.0001554706);
diffrac->SetBinError(44,0.0001488753);
diffrac->SetBinError(45,0.0001514758);
diffrac->SetBinError(46,0.0001476574);
diffrac->SetBinError(47,0.0001391356);
diffrac->SetBinError(48,0.0001375267);
diffrac->SetBinError(49,0.0001319161);
diffrac->SetBinError(50,0.0001309255);
diffrac->SetBinError(51,0.0001344074);
diffrac->SetBinError(52,0.0001332148);
diffrac->SetBinError(53,0.0001261892);
diffrac->SetBinError(54,0.0001206082);
diffrac->SetBinError(55,0.0001176829);
diffrac->SetBinError(56,0.0001165354);
diffrac->SetBinError(57,0.0001123524);
diffrac->SetBinError(58,0.0001091677);
diffrac->SetBinError(59,0.0001091677);
diffrac->SetBinError(60,0.0001101633);
diffrac->SetBinError(61,9.997153e-05);
diffrac->SetBinError(62,0.0001020502);
diffrac->SetBinError(63,9.891584e-05);
diffrac->SetBinError(64,9.40856e-05);
diffrac->SetBinError(65,8.727835e-05);
diffrac->SetBinError(66,9.475273e-05);
diffrac->SetBinError(67,8.528287e-05);
diffrac->SetBinError(68,8.503632e-05);
diffrac->SetBinError(69,8.59206e-05);
diffrac->SetBinError(70,7.899458e-05);
diffrac->SetBinError(71,8.026026e-05);
diffrac->SetBinError(72,7.428234e-05);
diffrac->SetBinError(73,7.629033e-05);
diffrac->SetBinError(74,7.145862e-05);
diffrac->SetBinError(75,6.845733e-05);
diffrac->SetBinError(76,7.291296e-05);
diffrac->SetBinError(77,7.16347e-05);
diffrac->SetBinError(78,6.512512e-05);
diffrac->SetBinError(79,6.168108e-05);
diffrac->SetBinError(80,5.882358e-05);
diffrac->SetBinError(81,5.995491e-05);
diffrac->SetBinError(82,6.154477e-05);
diffrac->SetBinError(83,5.6493e-05);
diffrac->SetBinError(84,5.406205e-05);
diffrac->SetBinError(85,5.061192e-05);
diffrac->SetBinError(86,5.375045e-05);
diffrac->SetBinError(87,5.077759e-05);
diffrac->SetBinError(88,5.413967e-05);
diffrac->SetBinError(89,5.069482e-05);
diffrac->SetBinError(90,4.681911e-05);
diffrac->SetBinError(91,4.442596e-05);
diffrac->SetBinError(92,4.139215e-05);
diffrac->SetBinError(93,4.6187e-05);
diffrac->SetBinError(94,4.259217e-05);
diffrac->SetBinError(95,3.789666e-05);
diffrac->SetBinError(96,3.756276e-05);
diffrac->SetBinError(97,3.596361e-05);
diffrac->SetBinError(98,3.392057e-05);
diffrac->SetBinError(99,3.161377e-05);
diffrac->SetBinError(100,2.839478e-05);
diffrac->SetBinError(101,2.969597e-05);
diffrac->SetBinError(102,2.718593e-05);
diffrac->SetBinError(103,3.121274e-05);
diffrac->SetBinError(104,2.318424e-05);
diffrac->SetBinError(105,3.025632e-05);
diffrac->SetBinError(106,2.49298e-05);
diffrac->SetBinError(107,2.263434e-05);
diffrac->SetBinError(108,2.069607e-05);
diffrac->SetBinError(109,2.089799e-05);
diffrac->SetBinError(110,1.986789e-05);
diffrac->SetBinError(111,2.671852e-05);
diffrac->SetBinError(112,2.007814e-05);
diffrac->SetBinError(113,2.509768e-05);
diffrac->SetBinError(114,1.689827e-05);
diffrac->SetBinError(115,2.028621e-05);
diffrac->SetBinError(116,1.832875e-05);
diffrac->SetBinError(117,1.613555e-05);
diffrac->SetBinError(118,1.328044e-05);
diffrac->SetBinError(119,1.359297e-05);
diffrac->SetBinError(120,1.359297e-05);
diffrac->SetBinError(121,1.296038e-05);
diffrac->SetBinError(122,1.389846e-05);
diffrac->SetBinError(123,1.194888e-05);
diffrac->SetBinError(124,1.122402e-05);
diffrac->SetBinError(125,1.159212e-05);
diffrac->SetBinError(126,1.22953e-05);
diffrac->SetBinError(127,1.003907e-05);
diffrac->SetBinError(128,1.0449e-05);
diffrac->SetBinError(129,8.196866e-06);
diffrac->SetBinError(130,7.667466e-06);
diffrac->SetBinError(131,9.611678e-06);
diffrac->SetBinError(132,8.69409e-06);
diffrac->SetBinError(133,1.003907e-05);
diffrac->SetBinError(134,7.098695e-06);
diffrac->SetBinError(135,7.098695e-06);
diffrac->SetBinError(136,5.019535e-06);
diffrac->SetBinError(137,6.480192e-06);
diffrac->SetBinError(138,6.480192e-06);
diffrac->SetBinError(139,2.89803e-06);
diffrac->SetBinError(140,4.098433e-06);
diffrac->SetBinError(141,5.019535e-06);
diffrac->SetBinError(142,5.79606e-06);
diffrac->SetBinError(143,2.89803e-06);
diffrac->SetBinError(144,5.79606e-06);
diffrac->SetBinError(145,2.89803e-06);
diffrac->SetBinError(146,2.89803e-06);
diffrac->SetBinError(147,4.098433e-06);
diffrac->SetBinError(149,5.019535e-06);
diffrac->SetBinError(151,4.098433e-06);
diffrac->SetBinError(152,5.019535e-06);
diffrac->SetBinError(154,2.89803e-06);
diffrac->SetBinError(155,2.89803e-06);
diffrac->SetBinError(157,4.098433e-06);
diffrac->SetBinError(170,2.89803e-06);
diffrac->SetBinError(173,4.098433e-06);
diffrac->SetBinError(178,2.89803e-06);
diffrac->SetEntries(345062);
diffrac->SetDirectory(0);
diffrac->SetFillColor(1);
diffrac->SetFillStyle(0);
diffrac->SetLineStyle(2);
diffrac->SetLineWidth(3);
diffrac->SetMarkerStyle(20);
diffrac->GetXaxis()->SetTitle("M");
diffrac->GetXaxis()->SetLabelFont(42);
diffrac->GetXaxis()->SetLabelOffset(0.01);
diffrac->GetXaxis()->SetLabelSize(0.045);
diffrac->GetXaxis()->SetTitleSize(0.055);
diffrac->GetXaxis()->SetTitleFont(42);
diffrac->GetYaxis()->SetTitle("Fraction of events");
diffrac->GetYaxis()->SetLabelFont(42);
diffrac->GetYaxis()->SetLabelOffset(0.01);
diffrac->GetYaxis()->SetLabelSize(0.045);
diffrac->GetYaxis()->SetTitleSize(0.055);
diffrac->GetYaxis()->SetTitleOffset(1.6);
diffrac->GetYaxis()->SetTitleFont(42);
diffrac->GetZaxis()->SetLabelFont(42);
diffrac->GetZaxis()->SetLabelSize(0.045);
diffrac->GetZaxis()->SetTitleFont(42);
diffrac->Draw("p same");
TLegend *leg = new TLegend(0.61,0.67,0.91,0.87,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(62);
leg->SetTextSize(0.035);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(19);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("diffrac","PYTHIA 7 TeV","l");
entry->SetLineColor(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
//entry=leg->AddEntry("diffrac","(Atlas tune)","");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry=leg->AddEntry("diffrac","PHOJET 7 TeV","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
leg->Draw();
printFinalCanvases(cEvtSel,"diffrac");
/*
TLatex * tex = new TLatex(0.85,0.9,"CMS");
tex->SetNDC();
tex->SetTextSize(0.04);
tex->SetLineWidth(2);
tex->Draw();
cEvtSel->Modified();
cEvtSel->cd();
cEvtSel->SetSelected(cEvtSel);
*/
}
void makePlots()
{
TFile * f2 = TFile::Open( "L1ITMuonBarrelPlots2.root");
TFile * f02 = TFile::Open( "L1ITMuonBarrelPlots0.2.root");
TFile * f005 = TFile::Open( "L1ITMuonBarrelPlots0.05.root");
TFile * f002 = TFile::Open( "L1ITMuonBarrelPlots0.02.root");
d2 = (TDirectory *) f2->Get( "L1ITMuPlotter");
d02 = (TDirectory *) f02->Get( "L1ITMuPlotter");
d005 = (TDirectory *) f005->Get( "L1ITMuPlotter");
d002 = (TDirectory *) f002->Get( "L1ITMuPlotter");
TString name;
//////////////////
/// Confirmed
drawConfirmed( "confirmed_st1" );
drawConfirmed( "confirmed_st2" );
drawConfirmed( "confirmed_st3" );
drawConfirmed( "confirmed_st4" );
TH1* conf = d005->Get( "confirmed_st1" )->Clone();
TH1* conf_2 = d005->Get( "confirmed_st2" );
TH1* conf_3 = d005->Get( "confirmed_st3" );
TH1* conf_4 = d005->Get( "confirmed_st4" );
conf->Add( conf_2 );
conf->Add( conf_3 );
conf->Add( conf_4 );
name = "confirmed_all";
TCanvas * confirmed = new TCanvas(name, name);
confirmed->cd();
// confirmed->SetLogy();
confirmed->SetGridy();
conf->SetLineWidth(2);
conf->Draw();
conf->Draw("sametext");
conf->SetTitle("confirmed in all stations");
conf->SetName("confirmed in all stations");
confirmed->SaveAs( name + "_dphi0.05.png" );
////////////////////
/// deltaPhiBin
name = "deltaPhiBin";
TCanvas * deltaPhiBin = new TCanvas(name, name);
deltaPhiBin->cd();
TH1* dphibin = d2->Get(name);
dphibin->SetLineWidth(2);
dphibin->Draw("text");
dphibin->Draw("same");
deltaPhiBin->SaveAs( name + ".png" );
////////////////////
/// deltaPhi
name = "deltaPhi";
TCanvas * deltaPhi = new TCanvas(name, name);
deltaPhi->cd();
TH1* dphi = d2->Get(name);
dphi->SetLineWidth(2);
dphi->Draw();
deltaPhi->SaveAs( name + ".png" );
////////////////////
/// deltaPhiDt
name = "deltaPhiDt";
TCanvas * deltaPhiDt = new TCanvas(name, name);
deltaPhiDt->cd();
TH1* dphiDt = d005->Get(name);
dphiDt->SetLineWidth(2);
dphiDt->Draw();
deltaPhiDt->SaveAs( name + ".png" );
/////////////////////////
/// rpcInHitsPerDtseg
name = "rpcInHitsPerDtseg";
TCanvas * rpcInHitsPerDtseg = new TCanvas(name, name);
rpcInHitsPerDtseg->SetGridy();
rpcInHitsPerDtseg->SetLogy();
rpcInHitsPerDtseg->cd();
TH1* rpcin = d005->Get(name);
rpcin->SetLineWidth(2);
rpcin->Draw();
rpcInHitsPerDtseg->SaveAs( name + "_dphi0.05.png" );
/////////////////////////
/// rpcOutHitsPerDtseg
name = "rpcOutHitsPerDtseg";
TCanvas * rpcOutHitsPerDtseg = new TCanvas(name, name);
rpcOutHitsPerDtseg->SetGridy();
rpcOutHitsPerDtseg->SetLogy();
rpcOutHitsPerDtseg->cd();
TH1* rpcout = d005->Get(name);
rpcout->SetLineWidth(2);
rpcout->Draw();
rpcOutHitsPerDtseg->SaveAs( name + "_dphi0.05.png" );
////////////////////
/// dtQuality_st1
name = "dtQuality_st1";
TCanvas * dtQuality_st1 = new TCanvas(name, name);
dtQuality_st1->SetGridx();
dtQuality_st1->SetGridy();
dtQuality_st1->cd();
// dtQuality_st1->SetLogz();
TH1* qual_st1 = d005->Get(name);
qual_st1->GetZaxis()->SetRangeUser( 0., 1. );
qual_st1->SetStats( 0 );
qual_st1->Draw("colztext");
dtQuality_st1->SaveAs( name + "_dphi0.05.png" );
////////////////////
/// dtQuality_st2
name = "dtQuality_st2";
TCanvas * dtQuality_st2 = new TCanvas(name, name);
dtQuality_st2->SetGridx();
dtQuality_st2->SetGridy();
dtQuality_st2->cd();
// dtQuality_st2->SetLogz();
TH1* qual_st2 = d005->Get(name);
qual_st2->GetZaxis()->SetRangeUser( 0., 1. );
qual_st2->SetStats( 0 );
qual_st2->Draw("colztext");
dtQuality_st2->SaveAs( name + "_dphi0.05.png" );
////////////////////
/// dtQualityNew
name = "dtQualityNew";
TCanvas * dtQualityNew = new TCanvas(name, name);
dtQualityNew->SetGridx();
dtQualityNew->SetGridy();
dtQualityNew->cd();
// dtQualityNew->SetLogz();
TH1* qualNew = d005->Get(name);
qualNew->GetZaxis()->SetRangeUser( 0., 0.5 );
qualNew->SetStats( 0 );
qualNew->Draw("colztext");
dtQualityNew->SaveAs( name + "_dphi0.05.png" );
}
