int fgtRawSpectraQA( const Char_t *filenameIn = "testfile.daq",
const Char_t *filebaseOut = "testfile.rawQA",
Int_t date = 20120115,
Int_t time = 0,
Float_t pedRelThres = 0, // 0 = no cut, otherwise: # of sigma above pedestal for cut
Int_t nevents = 2000,
Int_t timebin = 2,
UInt_t statusMask,
Bool_t cutShortEvents = 1 ){
LoadLibs();
Int_t ierr = 0;
//
// START CONSTRUCTING THE CHAIN
//
cout << "Constructing the chain" << endl;
analysisChain = new StChain("eemcAnalysisChain");
//
// THE DATABASE
//
// note: DB is used to convert elec coords into geoIds in
// StEvent/StFgtStrip, and do determine the disc and quad given
// an rdo/arm/apv combination.
cout << "Constructing St_db_Maker" << endl;
TString dir0 = "MySQL:StarDb";
TString dir1 = "$STAR/StarDb";
St_db_Maker *dbMkr = new St_db_Maker( "dbMkr", dir0, dir1 );
dbMkr->SetDateTime(date,time);
cout << "Constructing StFgtDbMaker" << endl;
fgtDbMkr = new StFgtDbMaker( "fgtDbMkr" );
//
// NOW THE OTHER READERS AND MAKERS
//
cout << "Constructing the daq reader" << endl;
daqRdr = new StFgtRawDaqReader( "daqReader", filenameIn, "" );
daqRdr->setIsCosmic( 0 );
daqRdr->cutShortEvents( cutShortEvents );
// a2cMkr = new StFgtA2CMaker( "a2cMkr" );
// a2cMkr->setStatusMask( statusMask );
// a2cMkr->setAbsThres( -10000 ); // set to below -4096 to skip cut
// a2cMkr->setRelThres( pedRelThres ); // set to zero to skip cut
cout << "Constructing the QA Makers" << endl;
qaMkr = new StFgtQaRawOctAdc( "qaMkr", 2 );
qaMkr->setTimeBin( 2 );
// debug
analysisChain->ls(4);
cout << "Initializing" << endl;
ierr = analysisChain->Init();
if( ierr ){
cout << "Error initializing" << endl;
return;
};
if( nevents < 0 )
nevents = 1<<30; // a big number
cout << "max nevents = " << nevents << endl;
for( int i=0; i<nevents && !ierr; ++i ){
if( i+1 % 100 == 0 )
cout << "\ton event number " << i << endl;
//cout << "clear" << endl;
analysisChain->Clear();
//cout << "make" << endl;
ierr = analysisChain->Make();
};
//
// Calls the ::Finish() method on all makers
//
cout << "finish" << endl;
analysisChain->Finish();
cout << "Making plots" << endl;
gROOT->SetStyle("Plain");
gStyle->SetOptStat(0);
gStyle->SetPalette(1);
can = new TCanvas( "fgtRawQA", "fgtRawQA", 850, 1100);
can->Divide( 1, 4 );
// open output file
can->Print( (std::string(filebaseOut) + ".ps[").data() );
//std::vector< TH2F* >& hist = qaMkr->getHistVec();
cout << "First hist at " << qaMkr->getHist(0) << endl;
cout << "Finding max" << endl;
// get max
Float_t max = 0;
for( Int_t idx = 0; idx < 48; ++idx ){
if( qaMkr->getHist(idx) )
if( qaMkr->getHist(idx)->GetMaximum() > max )
max = qaMkr->getHist(idx)->GetMaximum();
};
max = 0.9*max;
cout << "Max is " << max << endl;
TH2F *dummy[4];
const Char_t dummyNames[4][20] = { "dummy1", "dummy2", "dummy3", "dummy4" };
Float_t xMin = qaMkr->getHist(0)->GetXaxis()->GetXmin();
Float_t xMax = qaMkr->getHist(0)->GetXaxis()->GetXmax();
Float_t yMin = qaMkr->getHist(0)->GetYaxis()->GetXmin();
Float_t yMax = qaMkr->getHist(0)->GetYaxis()->GetXmax();
for( Int_t i=0; i<4; ++i )
dummy[i] = new TH2F( dummyNames[i], "", 5, xMin, xMax, 1, yMin, yMax );
std::stringstream ss;
Int_t subpad = 1;
idx = 0;
for( Int_t rdo = 1; rdo < 3; ++rdo ){
for( Int_t arm = 0; arm < 6; ++arm ){
for( Int_t startIdx = 0; startIdx < 4; ++startIdx, ++idx ){
for( Int_t oct = 0; oct < 1; ++oct ){
cout << "rdo/arm/disc/oct = " << rdo << '/' << arm << '/' << startIdx << '/' << oct << endl;
if( qaMkr->getHist(idx)->GetEntries() > 0 ){
can->cd(subpad);
gPad->SetLeftMargin( 0.06 );
gPad->SetRightMargin( 0.05 );
gPad->SetBottomMargin( 0.11 );
qaMkr->getHist(idx)->SetMaximum( max );
dummy[subpad-1]->GetXaxis()->SetTitleOffset(0.9);
dummy[subpad-1]->GetXaxis()->SetTitleSize(0.06);
dummy[subpad-1]->GetXaxis()->SetLabelSize(0.06);
dummy[subpad-1]->GetYaxis()->SetTitleOffset(0.55);
dummy[subpad-1]->GetYaxis()->SetTitleSize(0.06);
dummy[subpad-1]->GetYaxis()->SetLabelSize(0.06);
gPad->SetGridx(0);
gPad->SetGridy(0);
for( Int_t i = 0; i<5; ++i ){
ss.str("");
ss.clear();
ss << "channels in APV "; // startArray[startIdx]+i;
dummy[subpad-1]->GetXaxis()->SetBinLabel( i+1, ss.str().data() );
};
dummy[subpad-1]->GetXaxis()->SetNdivisions(222,0);
dummy[subpad-1]->SetMinimum( 0 );
dummy[subpad-1]->SetMaximum( max );
dummy[subpad-1]->SetTitle( qaMkr->getHist(idx)->GetTitle() );
qaMkr->getHist(idx)->SetTitle("");
dummy[subpad-1]->GetYaxis()->SetTitle( qaMkr->getHist(idx)->GetYaxis()->GetTitle() );
qaMkr->getHist(idx)->GetYaxis()->SetTitle("");
dummy[subpad-1]->Draw("COLZ");
qaMkr->getHist(idx)->Draw("COLZ SAME");
gPad->Update();
gPad->Modified();
TPaveText *title = (TPaveText*)(gPad->GetPrimitive("title"));
if( title ){
title->SetX1NDC( 0.045 );
title->SetX2NDC( 0.55 );
title->SetY1NDC( 0.91 ) ;
title->SetY2NDC( 0.999 );
title->SetBorderSize(0);
title->SetTextAlign( 12 );
title->SetTextColor(kBlue);
title->Draw();
};
TPave *palette = (TPave*)(gPad->GetPrimitive("palette"));
if( palette ){
palette->SetX1NDC( 0.955 );
palette->SetX2NDC( 0.985 );
palette->Draw();
};
++subpad;
if( subpad == 5 ){
subpad = 1;
can->Print( (std::string(filebaseOut) + ".ps").data() );
};
};
};
};
};
};
if( subpad != 1 )
can->Print( (std::string(filebaseOut) + ".ps").data() );
can->Print( (std::string(filebaseOut) + ".ps]").data() );
gSystem->Exec(( std::string("ps2pdf -dAutoRotatePages=/None ") + filebaseOut + ".ps" ).data());
cerr << "\tall done" << endl;
return;
};
void drawPlots(canvasSet_t& cs,bool savePlots2file)
{
wCanvas_t *wc0 = cs.canvases[0];
unsigned npads = wc0->npadsx*wc0->npadsy;
unsigned npadsall = cs.ncanvases*npads;
if (!npads) {
if (gl_verbose) cout << "Nothing to draw, guess I'm done." << endl;
return; // no pads to draw on.
} else if (!wc0->pads.size()) {
/********************************************************
* CHECK MULTIPAD OPTION, ASSIGN HISTOS TO PADS/CANVASES
********************************************************/
if (wc0->multipads.size()) {
npadsall = assignHistos2Multipad(cs);
} else {
cerr << "npads>0, but no pad specs supplied, exiting." << endl;
return; // no pads to draw on.
}
} else if (cs.ncanvases>1) {
npadsall = assignPads2Canvases(cs);
} else {
npadsall = std::min(npadsall,(unsigned)wc0->pads.size());
}
wc0->c1->cd();
if (gl_verbose)
cout << "Drawing on " << npadsall << " pad(s)" << endl;
wLegend_t *wleg = NULL;
/***************************************************
* LOOP OVER PADS...
***************************************************/
//vector<vector<string> >::const_iterator it;
for (unsigned ipad = 0; ipad< npadsall; ipad++) {
if (gl_verbose) cout << "Drawing pad# " << ipad+1 << endl;
unsigned ipadc = ipad % npads;
unsigned cnum = (ipad / npads) + 1;
wCanvas_t *wc = cs.canvases[cnum-1];
if (!ipadc) {
if (cnum-1) { // first canvas already created
wc->c1 = new TCanvas(wc->title.c_str(),wc->title.c_str(),
wc->padxdim*wc->npadsx,
wc->padydim*wc->npadsy);
float left = wc->leftmargin;
float bot = wc->bottommargin;
float rhgt = 1-wc->rightmargin;
float top = 1-wc->topmargin;
wc->motherpad = new TPad("mother","",left,bot,rhgt,top);
wc->c1->SetFillColor(wc->fillcolor);
wc->motherpad->SetFillColor(wc->fillcolor);
wc->motherpad->Draw();
wc->motherpad->cd();
wc->motherpad->Divide(wc->npadsx,wc->npadsy); // , wc->padxmargin,wc->padymargin);
}
/***************************************************
* CHECK FOR LATEX OBJECTS ON THE CANVAS
***************************************************/
wc->c1->cd();
for (unsigned j=0; j<wc->latex_ids.size(); j++) {
string& lid = wc->latex_ids[j];
map<string,TLatex *>::const_iterator it = glmap_id2latex.find(lid);
if (it == glmap_id2latex.end()) {
cerr << "ERROR: latex id " << lid << " never defined in layout" << endl;
exit (-1);
}
TLatex *ltx = it->second;
ltx->Draw();
wc->c1->Update();
}
}
wPad_t *& wp = wc->pads[ipadc];
wp->vp = wc->motherpad->cd(ipadc+1);
if (!wp->histo_ids.size() &&
!wp->stack_ids.size() &&
!wp->graph_ids.size() &&
!wp->macro_ids.size()) {
cerr << "ERROR: pad #" << ipadc+1 << " has no ids defined for it";
cerr << ", continuing to the next" << endl;
continue;
}
#if 0
/***************************************************
* Draw the frame first:
* (Fix up frame since it can't be auto-scaled:)
***************************************************/
string& hid0 = wp->histo_ids[0];
map<string,wTH1 *>::const_iterator it = glmap_id2histo.find(hid0);
if (it == glmap_id2histo.end()) {
cerr << "ERROR: id0 " << hid0 << " never defined in layout" << endl;
return;
}
wTH1 *myHisto = it->second;
TH1 *h = myHisto->histo();
if (wp->hframe->histo()->GetXaxis()->GetXmin() <=
wp->hframe->histo()->GetXaxis()->GetXmax())
wp->hframe->histo()->GetXaxis()->SetRangeUser(h->GetXaxis()->GetXmin(),
h->GetXaxis()->GetXmax());
if (wp->hframe->histo()->GetYaxis()->GetXmin() <=
wp->hframe->histo()->GetYaxis()->GetXmax())
wp->hframe->histo()->GetYaxis()->SetRangeUser(h->GetYaxis()->GetXmin(),
h->GetYaxis()->GetXmax());
wp->hframe->SetStats(0);
//wp->hframe->Draw("AXIS");
#endif
/***************************************************
* Check for external macros to run on the pad
***************************************************/
for (size_t i=0; i<wp->macro_ids.size(); i++) {
map<string,string>::const_iterator it = glmap_objpath2id.find(wp->macro_ids[i]);
if (it != glmap_objpath2id.end()) {
string path = it->second;
int error;
gROOT->Macro(path.c_str(), &error, kTRUE); // update current pad
if (error) {
static const char *errorstr[] = {
"kNoError","kRecoverable","kDangerous","kFatal","kProcessing" };
cerr << "ERROR: error returned from macro: " << errorstr[error] << endl;
}
} else {
cerr << "ERROR: macro id " << wp->macro_ids[i];
cerr << " never defined in layout" << endl;
}
}
/***************************************************
* Check for existence of a legend, create it
***************************************************/
bool drawlegend = false;
if (wp->legid.size()) {
map<string,wLegend_t *>::const_iterator it=glmap_id2legend.find(wp->legid);
if (it != glmap_id2legend.end()) {
drawlegend = true;
wleg = it->second;
} else {
cerr << "ERROR: legend id " << wp->legid;
cerr << " never defined in layout" << endl;
}
} else {
// Maybe gPad already *has* a legend from macros...
TPave *testing = (TPave *)gPad->GetPrimitive("TPave");
if (testing &&
!strcmp(testing->IsA()->GetName(),"TLegend")) {
TLegend *pullTheOtherOne = (TLegend *)testing;
if (gl_verbose) cout << "Found legend from macro" << endl;
wleg = new wLegend_t();
wleg->leg = pullTheOtherOne;
drawlegend = true;
}
}
/***************************************************
* LOOP OVER STACKS DEFINED FOR PAD...
***************************************************/
if (wp->stack_ids.size()) {
wStack_t *ws=NULL;
for (unsigned j = 0; j < wp->stack_ids.size(); j++) {
string& sid = wp->stack_ids[j];
map<string,wStack_t *>::const_iterator it = glmap_id2stack.find(sid);
if (it == glmap_id2stack.end()) {
cerr << "ERROR: stack id " << sid << " never defined in layout" << endl;
exit (-1);
}
bool firstInPad = !j;
ws = it->second;
if (!ws) { cerr<< "find returned NULL stack pointer for " << sid << endl; continue; }
// Add the histos in the stack to any legend that exists
//
if (drawlegend) {
for (size_t i=0; i<ws->v_histos.size(); i++) {
wTH1 *wh = ws->v_histos[i];
wh->ApplySavedStyle();
if(wh->GetLegendEntry().size())
wh->Add2Legend(wleg->leg);
}
}
string drawopt("");
if (ws->sum->GetDrawOption().size()) {
drawopt = ws->sum->GetDrawOption();
cout << "drawopt stored with histo = " << drawopt << endl;
}
drawInPad(wp, ws, firstInPad, drawopt);
wp->vp->Update();
}
} // stack loop
/***************************************************
* LOOP OVER HISTOS DEFINED FOR PAD...
***************************************************/
for (unsigned j = 0; j < wp->histo_ids.size(); j++) {
string& hid = wp->histo_ids[j];
map<string,wTH1 *>::const_iterator it = glmap_id2histo.find(hid);
if (it == glmap_id2histo.end()) {
cerr << "ERROR: histo id " << hid << " never defined in layout" << endl;
exit (-1);
}
wTH1 *myHisto = it->second;
if (myHisto) {
bool firstInPad = !j && !wp->stack_ids.size();
if (gl_verbose) {
cout << "Drawing " << hid << " => ";
cout << myHisto->histo()->GetName() << endl;
cout << "firstInPad = " << firstInPad << endl;
}
drawInPad(wp,*myHisto,firstInPad);
myHisto->DrawFits("same");
if (drawlegend && myHisto->GetLegendEntry().size()) {
if (wleg->drawoption.size()) myHisto->SetDrawOption(wleg->drawoption);
myHisto->Add2Legend(wleg->leg);
}
if (myHisto->statsAreOn()) {
myHisto->DrawStats();
wp->vp->Update();
}
myHisto->ApplySavedStyle();
wp->vp->Update();
}
} // histos loop
/***************************************************
* LOOP OVER HISTOS DEFINED FOR ALTERNATE Y-AXIS
***************************************************/
Float_t rightmax=0.0,rightmin=0.0;
Float_t scale=0.0;
for (unsigned j = 0; j < wp->altyh_ids.size(); j++) {
string& hid = wp->altyh_ids[j];
map<string,wTH1 *>::const_iterator it = glmap_id2histo.find(hid);
if (it == glmap_id2histo.end()) {
cerr << "ERROR: histo id " << hid << " never defined in layout" << endl;
exit (-1);
}
wTH1 *myHisto = it->second;
TH1 *h = myHisto->histo();
if (!j) {
//scale second set of histos to the pad coordinates
rightmin = h->GetMinimum();
rightmax = 1.1*h->GetMaximum();
scale = gPad->GetUymax()/rightmax;
}
TH1 *scaled=(TH1 *)h->Clone(Form("%s_%d",h->GetName(),ipad));
scaled->Scale(scale);
scaled->Draw("same");
//draw an axis on the right side
TGaxis *axis = new TGaxis(gPad->GetUxmax(), gPad->GetUymin(),
gPad->GetUxmax(), gPad->GetUymax(),
rightmin,rightmax,505,"+L");
axis->Draw();
gPad->Update();
if (drawlegend && myHisto->GetLegendEntry().size()) {
if (wleg->drawoption.size()) myHisto->SetDrawOption(wleg->drawoption);
myHisto->Add2Legend(wleg->leg);
}
}
/***************************************************
* LOOP OVER GRAPHS DEFINED FOR PAD...
***************************************************/
#if 0
TMultiGraph *mg;
if (graph_ids.size())
mg = new TMultiGraph();
#endif
for( unsigned j = 0; j < wp->graph_ids.size(); j++ ) {
string& gid = wp->graph_ids[j];
wGraph_t *wg = findGraph(gid);
bool firstInPad = !j && !wp->histo_ids.size() && !wp->macro_ids.size();
if( firstInPad && wg->gr && wg->gr->IsA()==TGraph::Class() )
wg->drawopt += string("A"); // no histos drawn, need to draw the frame ourselves.
if( wg && wg->gr ) {
// "pre-draw" in order to define the plot elements
wg->gr->Draw(wg->drawopt.c_str());
if (firstInPad) {
// Now we can set the axis attributes and range:
wg->gr->GetXaxis()->ImportAttributes(wg->xax);
wg->gr->GetYaxis()->ImportAttributes(wg->yax);
cout << wg->xax->GetXmin() << " " << wg->xax->GetXmax() << endl;
if( wg->xax->GetXmax()>wg->xax->GetXmin() )
wg->gr->GetXaxis()->SetLimits(wg->xax->GetXmin(),wg->xax->GetXmax());
if( wg->yax->GetXmax()>wg->yax->GetXmin() )
wg->gr->GetYaxis()->SetRangeUser(wg->yax->GetXmin(),wg->yax->GetXmax());
}
// draw for good
drawInPad<TGraph>(wp,wg->gr,wg->drawopt.c_str(),firstInPad);
wp->vp->Update();
if( wg->fitfn )
wg->gr->Fit(wg->fitfn);
if( drawlegend && wg->leglabel.size() )
wleg->leg->AddEntry(wg->gr,wg->leglabel.c_str(),wg->legdrawopt.c_str());
}
if( wg && wg->gr2d ) {
drawInPad<TGraph2D>(wp,wg->gr2d,wg->drawopt.c_str(),firstInPad);
if (firstInPad) {
// Now we can set the axis attributes and range:
wg->gr2d->GetXaxis()->ImportAttributes(wg->xax);
wg->gr2d->GetYaxis()->ImportAttributes(wg->yax);
cout << wg->xax->GetXmin() << " " << wg->xax->GetXmax() << endl;
if( wg->xax->GetXmax()>wg->xax->GetXmin() )
wg->gr2d->GetXaxis()->SetLimits(wg->xax->GetXmin(),wg->xax->GetXmax());
if( wg->yax->GetXmax()>wg->yax->GetXmin() )
wg->gr2d->GetYaxis()->SetRangeUser(wg->yax->GetXmin(),wg->yax->GetXmax());
}
if (wg->contours) {
//cout << "setting contours "; wg->contours->Print();
wg->gr2d->GetHistogram()->SetContour(wg->contours->GetNoElements(),
wg->contours->GetMatrixArray());
wg->gr2d->SetLineStyle (wg->lstyle);
wg->gr2d->SetLineColor (wg->lcolor);
wg->gr2d->SetLineWidth (wg->lwidth);
}
wp->vp->Modified();
wp->vp->Update();
if( drawlegend && wg->leglabel.size() )
wleg->leg->AddEntry(wg->gr2d,wg->leglabel.c_str(),wg->legdrawopt.c_str());
}
} // graph loop
/***************************************************
* LOOP OVER LINES DEFINED FOR PAD...
***************************************************/
for( unsigned j = 0; j < wp->line_ids.size(); j++ ) {
string drawopt("L");
string& lid = wp->line_ids[j];
map<string,TLine *>::const_iterator it2 = glmap_id2line.find(lid);
if (it2 == glmap_id2line.end()) {
cerr << "ERROR: line id " << lid << " never defined in layout" << endl;
exit (-1);
}
TLine *line = it2->second;
if (!j && !wp->histo_ids.size() && !wp->macro_ids.size())
drawopt += string("A"); // no histos drawn, need to draw the frame ourselves.
if (line) {
drawInPad<TLine>(wp,line,drawopt.c_str());
//if (drawlegend)
//wleg->leg->AddEntry(line,lid.c_str(),"L");
}
}
/***************************************************
* LOOP OVER BOXES DEFINED FOR PAD...
***************************************************/
for (unsigned j = 0; j < wp->box_ids.size(); j++) {
string drawopt("L");
string& bid = wp->box_ids[j];
map<string,TBox *>::const_iterator it2 = glmap_id2box.find(bid);
if (it2 == glmap_id2box.end()) {
cerr << "ERROR: box id " << bid << " never defined in layout" << endl;
exit (-1);
}
TBox *box = it2->second;
if (box) {
drawInPad<TBox>(wp,box,drawopt.c_str());
}
}
/***************************************************
* Draw the legend
***************************************************/
if (drawlegend) {
wleg->leg->Draw("same");
wp->vp->Update();
}
/***************************************************
* Draw each latex/label object
***************************************************/
for (unsigned j=0; j<wp->latex_ids.size(); j++) {
string& lid = wp->latex_ids[j];
map<string,TLatex *>::const_iterator it2 = glmap_id2latex.find(lid);
if (it2 == glmap_id2latex.end()) {
cerr << "ERROR: latex id " << lid << " never defined in layout" << endl;
exit (-1);
}
if (gl_verbose) cout << "Drawing latex object " << lid << endl;
TLatex *ltx = it2->second;
ltx->Draw();
wp->vp->Update();
}
for (unsigned j = 0; j < wp->label_ids.size(); j++) {
string& lid = wp->label_ids[j];
map<string,wLabel_t *>::const_iterator it2 = glmap_id2label.find(lid);
if (it2 == glmap_id2label.end()) {
cerr << "ERROR: label id " << lid << " never defined in layout" << endl;
exit (-1);
}
if (gl_verbose) cout << "Drawing label object " << lid << endl;
wLabel_t *wlab = it2->second;
drawStandardText(wlab->text, wlab->x1ndc, wlab->y1ndc,-1,-1,wlab->textsize);
wp->vp->Update();
}
wc->c1->Update();
} // pad loop
//prdFixOverlay();
if (savePlots2file) {
wc0 = cs.canvases[0];
if (!wc0->savenamefmts.size()) // define a default
wc0->savenamefmts.push_back("%F_%C.png");
for (size_t i=0; i<cs.canvases.size(); i++) {
wCanvas_t *wc = cs.canvases[i];
wc->c1->cd();
for (size_t j=0; j<wc0->savenamefmts.size(); j++)
saveCanvas2File(wc,wc0->savenamefmts[j]);
}
}
} // drawPlots
void limit()
{
//=========Macro generated from canvas: limit/limit
//========= (Thu Apr 27 14:38:33 2017) by ROOT version6.02/05
TCanvas *limit = new TCanvas("limit", "limit",0,0,600,600);
gStyle->SetOptFit(1);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
limit->SetHighLightColor(2);
limit->Range(0,0,1,1);
limit->SetFillColor(0);
limit->SetBorderMode(0);
limit->SetBorderSize(2);
limit->SetLeftMargin(0.16);
limit->SetRightMargin(0.04);
limit->SetTopMargin(0.06);
limit->SetBottomMargin(0.12);
limit->SetFrameFillStyle(0);
limit->SetFrameBorderMode(0);
// ------------>Primitives in pad: pad
TPad *pad = new TPad("pad", "pad",0,0,1,1);
pad->Draw();
pad->cd();
pad->Range(-59.99999,-2.046818,1565,15.01);
pad->SetFillColor(0);
pad->SetBorderMode(0);
pad->SetLogy();
pad->SetBorderSize(2);
pad->SetTickx(1);
pad->SetTicky(1);
pad->SetLeftMargin(0.16);
pad->SetRightMargin(0.04);
pad->SetTopMargin(0.06);
pad->SetBottomMargin(0.12);
pad->SetFrameFillStyle(0);
pad->SetFrameBorderMode(0);
pad->SetFrameFillStyle(0);
pad->SetFrameBorderMode(0);
TH1F *tmp01 = new TH1F("tmp01","Graph",100,200,1500);
tmp01->SetMinimum(0.);
tmp01->SetMaximum(13.98659);
tmp01->SetStats(0);
tmp01->SetLineStyle(0);
tmp01->SetMarkerStyle(20);
tmp01->GetXaxis()->SetTitle("m_{H} (GeV)");
tmp01->GetXaxis()->SetLabelFont(42);
tmp01->GetXaxis()->SetLabelOffset(0.01);
tmp01->GetXaxis()->SetTitleSize(0.05);
tmp01->GetXaxis()->SetTickLength(0.02);
tmp01->GetXaxis()->SetTitleOffset(1.08);
tmp01->GetXaxis()->SetTitleFont(42);
tmp01->GetYaxis()->SetTitle("95% CL limit on #sigma/#sigma_{SM}");
tmp01->GetYaxis()->SetLabelFont(42);
tmp01->GetYaxis()->SetLabelOffset(0.007);
tmp01->GetYaxis()->SetTitleSize(0.05);
tmp01->GetYaxis()->SetTickLength(0.02);
tmp01->GetYaxis()->SetTitleOffset(1.56);
tmp01->GetYaxis()->SetTitleFont(42);
tmp01->GetZaxis()->SetLabelFont(42);
tmp01->GetZaxis()->SetLabelOffset(0.007);
tmp01->GetZaxis()->SetTitleSize(0.05);
tmp01->GetZaxis()->SetTickLength(0.02);
tmp01->GetZaxis()->SetTitleFont(42);
tmp01->Draw("AXIS");
TH1F *tmp02 = new TH1F("tmp02","Graph",100,200,1500);
tmp02->SetMinimum(0);
tmp02->SetMaximum(13.98659);
tmp02->SetStats(0);
tmp02->SetLineStyle(0);
tmp02->SetMarkerStyle(20);
tmp02->GetXaxis()->SetTitle("m_{H} (GeV)");
tmp02->GetXaxis()->SetLabelFont(42);
tmp02->GetXaxis()->SetLabelOffset(0.01);
tmp02->GetXaxis()->SetTitleSize(0.05);
tmp02->GetXaxis()->SetTickLength(0.02);
tmp02->GetXaxis()->SetTitleOffset(1.08);
tmp02->GetXaxis()->SetTitleFont(42);
tmp02->GetYaxis()->SetTitle("95% CL limit on #sigma/#sigma_{SM}");
tmp02->GetYaxis()->SetLabelFont(42);
tmp02->GetYaxis()->SetLabelOffset(0.007);
tmp02->GetYaxis()->SetTitleSize(0.05);
tmp02->GetYaxis()->SetTickLength(0.02);
tmp02->GetYaxis()->SetTitleOffset(1.56);
tmp02->GetYaxis()->SetTitleFont(42);
tmp02->GetZaxis()->SetLabelFont(42);
tmp02->GetZaxis()->SetLabelOffset(0.007);
tmp02->GetZaxis()->SetTitleSize(0.05);
tmp02->GetZaxis()->SetTickLength(0.02);
tmp02->GetZaxis()->SetTitleFont(42);
tmp02->Draw("AXIGSAME");
Double_t Graph0_fx3001[27] = {
200,
250,
300,
350,
400,
450,
500,
550,
600,
650,
700,
750,
800,
850,
900,
950,
1000,
1050,
1100,
1150,
1200,
1250,
1300,
1350,
1400,
1450,
1500};
Double_t Graph0_fy3001[27] = {
4.796875,
2.148438,
0.8398438,
0.3857422,
0.2001953,
0.1088867,
0.06591797,
0.04736328,
0.03564453,
0.02783203,
0.02294922,
0.01806641,
0.01513672,
0.01220703,
0.01123047,
0.009277344,
0.008300781,
0.007324219,
0.007324219,
0.006347656,
0.006347656,
0.005371094,
0.005371094,
0.004394531,
0.004394531,
0.004394531,
0.004394531};
Double_t Graph0_felx3001[27] = {
4.032346e-316,
0,
0,
0,
6.906313e-310,
4.032343e-316,
2.121996e-314,
4.032582e-316,
0,
4.032588e-316,
4.032343e-316,
0,
4.032584e-316,
1.508018e-316,
1.508219e-316,
4.032589e-316,
2.124327e-314,
4.032593e-316,
5.582942e-322,
6.906313e-310,
0,
4.032344e-316,
4.032583e-316,
4.032588e-316,
1.58101e-322,
4.032293e-316,
5.582942e-322};
Double_t Graph0_fely3001[27] = {
2.370331,
1.057434,
0.405159,
0.1845837,
0.09618759,
0.05231667,
0.03192902,
0.02331161,
0.01754379,
0.01435089,
0.01183319,
0.009315491,
0.008041382,
0.006484985,
0.005966187,
0.004928589,
0.00440979,
0.004119873,
0.004119873,
0.003570557,
0.003570557,
0.00302124,
0.00302124,
0.002471924,
0.002471924,
0.002471924,
0.002471924};
Double_t Graph0_fehx3001[27] = {
4.031643e-316,
0,
0,
0,
6.906313e-310,
4.032106e-316,
2.121996e-314,
4.032349e-316,
0,
4.032354e-316,
4.032106e-316,
0,
4.032351e-316,
1.508018e-316,
1.508219e-316,
4.032356e-316,
2.124327e-314,
4.03236e-316,
5.582942e-322,
6.906313e-310,
0,
4.032107e-316,
4.03235e-316,
4.032354e-316,
1.58101e-322,
4.032222e-316,
5.582942e-322};
Double_t Graph0_fehy3001[27] = {
5.693069,
2.47175,
0.8989364,
0.4081147,
0.2108151,
0.1170805,
0.07362081,
0.05438244,
0.04146756,
0.03243252,
0.02825768,
0.02361592,
0.01990321,
0.01752287,
0.01712826,
0.0153559,
0.01456667,
0.01341494,
0.01192601,
0.0124739,
0.01149977,
0.01055484,
0.01051202,
0.008670811,
0.008670811,
0.008635777,
0.008635777};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(27,Graph0_fx3001,Graph0_fy3001,Graph0_felx3001,Graph0_fehx3001,Graph0_fely3001,Graph0_fehy3001);
grae->SetName("Graph0");
grae->SetTitle("Graph");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ffff00");
grae->SetFillColor(ci);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph3001 = new TH1F("Graph_Graph3001","Graph",100,70,1630);
Graph_Graph3001->SetMinimum(0);
Graph_Graph3001->SetMaximum(11.53875);
Graph_Graph3001->SetDirectory(0);
Graph_Graph3001->SetStats(0);
Graph_Graph3001->SetLineStyle(0);
Graph_Graph3001->SetMarkerStyle(20);
Graph_Graph3001->GetXaxis()->SetLabelFont(42);
Graph_Graph3001->GetXaxis()->SetTitleSize(0.05);
Graph_Graph3001->GetXaxis()->SetTickLength(0.02);
Graph_Graph3001->GetXaxis()->SetTitleOffset(1.08);
Graph_Graph3001->GetXaxis()->SetTitleFont(42);
Graph_Graph3001->GetYaxis()->SetLabelFont(42);
Graph_Graph3001->GetYaxis()->SetLabelOffset(0.007);
Graph_Graph3001->GetYaxis()->SetTitleSize(0.05);
Graph_Graph3001->GetYaxis()->SetTickLength(0.02);
Graph_Graph3001->GetYaxis()->SetTitleOffset(1.56);
Graph_Graph3001->GetYaxis()->SetTitleFont(42);
Graph_Graph3001->GetZaxis()->SetLabelFont(42);
Graph_Graph3001->GetZaxis()->SetLabelOffset(0.007);
Graph_Graph3001->GetZaxis()->SetTitleSize(0.05);
Graph_Graph3001->GetZaxis()->SetTickLength(0.02);
Graph_Graph3001->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_Graph3001);
grae->Draw("3");
Double_t Graph1_fx3002[27] = {
200,
250,
300,
350,
400,
450,
500,
550,
600,
650,
700,
750,
800,
850,
900,
950,
1000,
1050,
1100,
1150,
1200,
1250,
1300,
1350,
1400,
1450,
1500};
Double_t Graph1_fy3002[27] = {
4.796875,
2.148438,
0.8398438,
0.3857422,
0.2001953,
0.1088867,
0.06591797,
0.04736328,
0.03564453,
0.02783203,
0.02294922,
0.01806641,
0.01513672,
0.01220703,
0.01123047,
0.009277344,
0.008300781,
0.007324219,
0.007324219,
0.006347656,
0.006347656,
0.005371094,
0.005371094,
0.004394531,
0.004394531,
0.004394531,
0.004394531};
Double_t Graph1_felx3002[27] = {
4.028644e-316,
2.608667e-321,
4.022558e-316,
0,
2.420922e-322,
4.028601e-316,
0,
2.11372e-314,
0,
4.940656e-323,
8.744962e-322,
0,
0,
4.940656e-323,
0,
6.906314e-310,
5.33769e+160,
0,
4.021573e-316,
6.906314e-310,
0,
1.188318e-312,
4.027998e-316,
0,
1.58101e-322,
4.027629e-316,
1.630417e-322};
Double_t Graph1_fely3002[27] = {
1.490716,
0.6650269,
0.2516418,
0.1146438,
0.05974151,
0.03228919,
0.01945674,
0.01420552,
0.01069075,
0.009417772,
0.007765532,
0.006113291,
0.005277157,
0.003647804,
0.00335598,
0.002772331,
0.002480507,
0.002317429,
0.002832413,
0.002231598,
0.002231598,
0.001888275,
0.001888275,
0.001544952,
0.001544952,
0.001544952,
0.002162933};
Double_t Graph1_fehx3002[27] = {
1.179415e-312,
2.608667e-321,
4.022558e-316,
0,
2.420922e-322,
4.028634e-316,
0,
2.113722e-314,
0,
4.940656e-323,
8.744962e-322,
0,
0,
4.940656e-323,
0,
6.906314e-310,
2.139459e+161,
0,
4.021573e-316,
6.906314e-310,
0,
1.188318e-312,
4.027998e-316,
0,
1.58101e-322,
4.027675e-316,
1.630417e-322};
Double_t Graph1_fehy3002[27] = {
2.351828,
1.053343,
0.3849801,
0.1768222,
0.09017253,
0.05078118,
0.03126747,
0.02303264,
0.01790216,
0.012869,
0.0113431,
0.009505779,
0.008205643,
0.007395977,
0.0068043,
0.006212621,
0.006088058,
0.005138258,
0.004204029,
0.005262822,
0.003643492,
0.004453157,
0.003768056,
0.004204029,
0.004204029,
0.003643492,
0.003643492};
grae = new TGraphAsymmErrors(27,Graph1_fx3002,Graph1_fy3002,Graph1_felx3002,Graph1_fehx3002,Graph1_fely3002,Graph1_fehy3002);
grae->SetName("Graph1");
grae->SetTitle("Graph");
ci = TColor::GetColor("#00ff00");
grae->SetFillColor(ci);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph3002 = new TH1F("Graph_Graph3002","Graph",100,-8.010917e+160,2.406781e+161);
Graph_Graph3002->SetMinimum(0);
Graph_Graph3002->SetMaximum(7.86335);
Graph_Graph3002->SetDirectory(0);
Graph_Graph3002->SetStats(0);
Graph_Graph3002->SetLineStyle(0);
Graph_Graph3002->SetMarkerStyle(20);
Graph_Graph3002->GetXaxis()->SetLabelFont(42);
Graph_Graph3002->GetXaxis()->SetTitleSize(0.05);
Graph_Graph3002->GetXaxis()->SetTickLength(0.02);
Graph_Graph3002->GetXaxis()->SetTitleOffset(1.08);
Graph_Graph3002->GetXaxis()->SetTitleFont(42);
Graph_Graph3002->GetYaxis()->SetLabelFont(42);
Graph_Graph3002->GetYaxis()->SetLabelOffset(0.007);
Graph_Graph3002->GetYaxis()->SetTitleSize(0.05);
Graph_Graph3002->GetYaxis()->SetTickLength(0.02);
Graph_Graph3002->GetYaxis()->SetTitleOffset(1.56);
Graph_Graph3002->GetYaxis()->SetTitleFont(42);
Graph_Graph3002->GetZaxis()->SetLabelFont(42);
Graph_Graph3002->GetZaxis()->SetLabelOffset(0.007);
Graph_Graph3002->GetZaxis()->SetTitleSize(0.05);
Graph_Graph3002->GetZaxis()->SetTickLength(0.02);
Graph_Graph3002->GetZaxis()->SetTitleFont(42);
grae->SetHistogram(Graph_Graph3002);
grae->Draw("3");
Double_t Graph2_fx1[27] = {
200,
250,
300,
350,
400,
450,
500,
550,
600,
650,
700,
750,
800,
850,
900,
950,
1000,
1050,
1100,
1150,
1200,
1250,
1300,
1350,
1400,
1450,
1500};
Double_t Graph2_fy1[27] = {
4.796875,
2.148438,
0.8398438,
0.3857422,
0.2001953,
0.1088867,
0.06591797,
0.04736328,
0.03564453,
0.02783203,
0.02294922,
0.01806641,
0.01513672,
0.01220703,
0.01123047,
0.009277344,
0.008300781,
0.007324219,
0.007324219,
0.006347656,
0.006347656,
0.005371094,
0.005371094,
0.004394531,
0.004394531,
0.004394531,
0.004394531};
TGraph *graph = new TGraph(27,Graph2_fx1,Graph2_fy1);
graph->SetName("Graph2");
graph->SetTitle("Graph");
graph->SetFillColor(1);
ci = TColor::GetColor("#ff0000");
graph->SetLineColor(ci);
graph->SetLineWidth(2);
graph->SetMarkerStyle(20);
TH1F *Graph_Graph1 = new TH1F("Graph_Graph1","Graph",100,70,1630);
Graph_Graph1->SetMinimum(0);
Graph_Graph1->SetMaximum(5.276123);
Graph_Graph1->SetDirectory(0);
Graph_Graph1->SetStats(0);
Graph_Graph1->SetLineStyle(0);
Graph_Graph1->SetMarkerStyle(20);
Graph_Graph1->GetXaxis()->SetLabelFont(42);
Graph_Graph1->GetXaxis()->SetTitleSize(0.05);
Graph_Graph1->GetXaxis()->SetTickLength(0.02);
Graph_Graph1->GetXaxis()->SetTitleOffset(1.08);
Graph_Graph1->GetXaxis()->SetTitleFont(42);
Graph_Graph1->GetYaxis()->SetLabelFont(42);
Graph_Graph1->GetYaxis()->SetLabelOffset(0.007);
Graph_Graph1->GetYaxis()->SetTitleSize(0.05);
Graph_Graph1->GetYaxis()->SetTickLength(0.02);
Graph_Graph1->GetYaxis()->SetTitleOffset(1.56);
Graph_Graph1->GetYaxis()->SetTitleFont(42);
Graph_Graph1->GetZaxis()->SetLabelFont(42);
Graph_Graph1->GetZaxis()->SetLabelOffset(0.007);
Graph_Graph1->GetZaxis()->SetTitleSize(0.05);
Graph_Graph1->GetZaxis()->SetTickLength(0.02);
Graph_Graph1->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph1);
graph->Draw("l");
Double_t Graph3_fx2[27] = {
200,
250,
300,
350,
400,
450,
500,
550,
600,
650,
700,
750,
800,
850,
900,
950,
1000,
1050,
1100,
1150,
1200,
1250,
1300,
1350,
1400,
1450,
1500};
Double_t Graph3_fy2[27] = {
4.79203,
2.152691,
0.8377939,
0.3861777,
0.1995838,
0.1090425,
0.0660598,
0.04759021,
0.03596779,
0.02751563,
0.02238263,
0.01815241,
0.01512282,
0.01228406,
0.01162986,
0.009836298,
0.008643669,
0.007866618,
0.007089044,
0.006702567,
0.005960677,
0.005744891,
0.005145786,
0.004956617,
0.004917011,
0.004476165,
0.004343825};
graph = new TGraph(27,Graph3_fx2,Graph3_fy2);
graph->SetName("Graph3");
graph->SetTitle("Graph");
graph->SetFillColor(1);
graph->SetLineWidth(2);
graph->SetMarkerStyle(20);
TH1F *Graph_Graph2 = new TH1F("Graph_Graph2","Graph",100,70,1630);
Graph_Graph2->SetMinimum(0);
Graph_Graph2->SetMaximum(5.270799);
Graph_Graph2->SetDirectory(0);
Graph_Graph2->SetStats(0);
Graph_Graph2->SetLineStyle(0);
Graph_Graph2->SetMarkerStyle(20);
Graph_Graph2->GetXaxis()->SetLabelFont(42);
Graph_Graph2->GetXaxis()->SetTitleSize(0.05);
Graph_Graph2->GetXaxis()->SetTickLength(0.02);
Graph_Graph2->GetXaxis()->SetTitleOffset(1.08);
Graph_Graph2->GetXaxis()->SetTitleFont(42);
Graph_Graph2->GetYaxis()->SetLabelFont(42);
Graph_Graph2->GetYaxis()->SetLabelOffset(0.007);
Graph_Graph2->GetYaxis()->SetTitleSize(0.05);
Graph_Graph2->GetYaxis()->SetTickLength(0.02);
Graph_Graph2->GetYaxis()->SetTitleOffset(1.56);
Graph_Graph2->GetYaxis()->SetTitleFont(42);
Graph_Graph2->GetZaxis()->SetLabelFont(42);
Graph_Graph2->GetZaxis()->SetLabelOffset(0.007);
Graph_Graph2->GetZaxis()->SetTitleSize(0.05);
Graph_Graph2->GetZaxis()->SetTickLength(0.02);
Graph_Graph2->GetZaxis()->SetTitleFont(42);
graph->SetHistogram(Graph_Graph2);
graph->Draw("pl");
TH1F *tmp0_copy3 = new TH1F("tmp0_copy3","Graph",100,200,1500);
tmp0_copy3->SetMinimum(0);
tmp0_copy3->SetMaximum(11.53875);
tmp0_copy3->SetDirectory(0);
tmp0_copy3->SetStats(0);
tmp0_copy3->SetLineStyle(0);
tmp0_copy3->SetMarkerStyle(20);
tmp0_copy3->GetXaxis()->SetLabelFont(42);
tmp0_copy3->GetXaxis()->SetTitleSize(0.05);
tmp0_copy3->GetXaxis()->SetTickLength(0.02);
tmp0_copy3->GetXaxis()->SetTitleOffset(1.08);
tmp0_copy3->GetXaxis()->SetTitleFont(42);
tmp0_copy3->GetYaxis()->SetLabelFont(42);
tmp0_copy3->GetYaxis()->SetLabelOffset(0.007);
tmp0_copy3->GetYaxis()->SetTitleSize(0.05);
tmp0_copy3->GetYaxis()->SetTickLength(0.02);
tmp0_copy3->GetYaxis()->SetTitleOffset(1.56);
tmp0_copy3->GetYaxis()->SetTitleFont(42);
tmp0_copy3->GetZaxis()->SetLabelFont(42);
tmp0_copy3->GetZaxis()->SetLabelOffset(0.007);
tmp0_copy3->GetZaxis()->SetTitleSize(0.05);
tmp0_copy3->GetZaxis()->SetTickLength(0.02);
tmp0_copy3->GetZaxis()->SetTitleFont(42);
tmp0_copy3->Draw("sameaxis");
TH1F *tmp0_copy4 = new TH1F("tmp0_copy4","Graph",100,200,1500);
tmp0_copy4->SetMinimum(0);
tmp0_copy4->SetMaximum(11.53875);
tmp0_copy4->SetDirectory(0);
tmp0_copy4->SetStats(0);
tmp0_copy4->SetLineStyle(0);
tmp0_copy4->SetMarkerStyle(20);
tmp0_copy4->GetXaxis()->SetLabelFont(42);
tmp0_copy4->GetXaxis()->SetTitleSize(0.05);
tmp0_copy4->GetXaxis()->SetTickLength(0.02);
tmp0_copy4->GetXaxis()->SetTitleOffset(1.08);
tmp0_copy4->GetXaxis()->SetTitleFont(42);
tmp0_copy4->GetYaxis()->SetLabelFont(42);
tmp0_copy4->GetYaxis()->SetLabelOffset(0.007);
tmp0_copy4->GetYaxis()->SetTitleSize(0.05);
tmp0_copy4->GetYaxis()->SetTickLength(0.02);
tmp0_copy4->GetYaxis()->SetTitleOffset(1.56);
tmp0_copy4->GetYaxis()->SetTitleFont(42);
tmp0_copy4->GetZaxis()->SetLabelFont(42);
tmp0_copy4->GetZaxis()->SetLabelOffset(0.007);
tmp0_copy4->GetZaxis()->SetTitleSize(0.05);
tmp0_copy4->GetZaxis()->SetTickLength(0.02);
tmp0_copy4->GetZaxis()->SetTitleFont(42);
tmp0_copy4->Draw("sameaxig");
TLegend *leg = new TLegend(0.495,0.825,0.945,0.925,NULL,"NBNDC");
leg->SetBorderSize(0);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
TLegendEntry *entry=leg->AddEntry("Graph3","Observed","LP");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph2","Expected","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph1","#pm1#sigma Expected","F");
ci = TColor::GetColor("#00ff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph0","#pm2#sigma Expected","F");
ci = TColor::GetColor("#ffff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
TPave *pave = new TPave(0.16,0.81,0.96,0.94,1,"brNDC");
pave->SetFillColor(0);
pave->Draw();
leg = new TLegend(0.495,0.825,0.945,0.925,NULL,"NBNDC");
leg->SetBorderSize(0);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
entry=leg->AddEntry("Graph3","Observed","LP");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph2","Expected","L");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph1","#pm1#sigma Expected","F");
ci = TColor::GetColor("#00ff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
entry=leg->AddEntry("Graph0","#pm2#sigma Expected","F");
ci = TColor::GetColor("#ffff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(42);
leg->Draw();
TLatex * tex = new TLatex(0.196,0.9113,"CMS");
tex->SetNDC();
tex->SetTextAlign(13);
tex->SetTextSize(0.048);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.196,0.8537,"Internal");
tex->SetNDC();
tex->SetTextAlign(13);
tex->SetTextFont(52);
tex->SetTextSize(0.03648);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.96,0.952,"");
tex->SetNDC();
tex->SetTextAlign(31);
tex->SetTextFont(42);
tex->SetTextSize(0.036);
tex->SetLineWidth(2);
tex->Draw();
tex = new TLatex(0.16,0.952,"");
tex->SetNDC();
tex->SetTextFont(42);
tex->SetTextSize(0.036);
tex->SetLineWidth(2);
tex->Draw();
pad->Modified();
limit->cd();
limit->Modified();
limit->cd();
limit->SetSelected(limit);
}
void plotsNoFitETau(bool mine=false,TString name="collMass",TString file="LFV_vbf_collMass_fakeRate_zjetsEmbed_newSignal.root", TString dir="preselectionMetFix_Nov2", TString dirInternal="vbfmutau",double blindA=100, double blindB=160, bool blind=false, TString Xaxis="M_{#mu,#tau}_{coll} [GeV]", TString Yaxis="Events / 20 GeV", TString Tex="#mu#tau_{h} 0 Jet", double xmin=0, double xmax=300, int rebinning=1, bool setLogY=false, double legx1=0.6, double legy1=0.9, double legx2=0.9, double legy2=0.5, double MAX=-1, double MIN=-1, double CORRFR=1., double scaleSignal=1.0, double ymin=-0.95, double ymax=0.95){
//gROOT->LoadMacro("tdrstyle.C");
setTDRStyle();
//gROOT->LoadMacro("CMS_lumi.C");
writeExtraText = false;
int iPeriod = 2;
// second parameter in example_plot is iPos, which drives the position of the CMS logo in the plot
int iPos=11;// : top-left, left-aligned
// iPos=33 : top-right, right-aligned
// iPos=22 : center, centered
// mode generally :
// iPos = 10*(alignement 1/2/3) + position (1/2/3 = left/center/right)
// example_plot( iPeriod, 11 ); // left-aligned
// example_plot( iPeriod, 33 ); // right-aligned
// example_plot( iPeriod, 0 ); // out of frame (in exceptional cases)
// example_plot( iPeriod, 11 ); // default: left-aligned
// example_plot( iPeriod, 22 ); // centered
// example_plot( iPeriod, 33 ); // right-aligned
double kForPlotting=1;
double branchingratioTauTau=0.063;
double branchingratioTauMu=0.1;
double Lumi=19717;
// Get Plots
// Get Plots
TFile *_file0= new TFile(dir+"/"+file);
TString find="DYrightarrowlljets";
if(!mine) find="zjetsother";
TH1F *hDY=_file0->Get(dirInternal+find); hDY->SetName("DY");
find="Diboson";
if(!mine) find="ww";
TH1F *hWW=_file0->Get(dirInternal+find); hWW->SetName("WW");
find="SingleTop";
if(!mine) find="singlet";
TH1F *hTOP=_file0->Get(dirInternal+find); hTOP->SetName("TOP");
find="tbart";
if(!mine) find="ttbar";
TH1F *hTT=_file0->Get(dirInternal+find); hTT->SetName("TT");
find="Zrightarrowtautauembedded";
if(!mine) find="ztautau";
TH1F *hZTauTau=_file0->Get(dirInternal+find); hZTauTau->SetName("ZTauTau");
find="Fakes";
if(!mine) find="fakes";
TH1F *hFAKES=_file0->Get(dirInternal+find); hFAKES->SetName("fakes");
hFAKES->Scale(CORRFR);
TH1F *hLFVVBF126=_file0->Get(dirInternal+"LFVvbfHiggsBR=1"); hLFVVBF126->SetName("LFVVBF126");
TH1F *hSMVBF126=_file0->Get(dirInternal+"SMvbfHiggs"); hSMVBF126->SetName("SMVBF126");
TH1F *hLFVGG126=_file0->Get(dirInternal+"LFVggHiggsBR=1"); hLFVGG126->SetName("LFVGG126");
TH1F *hSMGG126=_file0->Get(dirInternal+"SMHiggs"); hSMGG126->SetName("SMGG126");
find="Observed";
if(!mine) find="data_obs";
TH1F* hdata_obsNoErrFix=_file0->Get(dirInternal+find);
TH1F* hdata_obs = hdata_obsNoErrFix->Clone();
std::cout << hdata_obsNoErrFix->GetEntries() << std::endl;
for (i = 1; i <= hdata_obs->GetNbinsX(); i++){
if (hdata_obs->GetBinContent(i) == 0){
hdata_obs->SetBinContent(i,0.0);
hdata_obs->SetBinError(i,1.8);
}
}
hdata_obs->SetName("data_obs");
TH1F* hSMHIGGS=hSMGG126->Clone(); hSMHIGGS->SetName("HIGGSSM");
hSMHIGGS->Add(hSMVBF126);
// For the Control plots only we want signal to be scaled to 100% Br
hLFVVBF126->Scale(scaleSignal);
hLFVGG126->Scale(scaleSignal);
// From the fit
/* hDY->Scale(1.122280);
hWW->Scale(1.254172);
hTOP->Scale(1.348162);
hTT->Scale(1.082784);
hZTauTau->Scale(1.157596);
hFAKES->Scale(0.901000);
*/
// Daniel's Colors
/*
WGammaStar=kCyan
ZTauTauEmbedded=kOrange-4
ZLL_residual=kAzure+3
TTBar=40
SingleTop=kGreen-2
EWKDiBoson=kRed+2
WJets/QCD Multijets=kMagenta-10
SMHToTauTau=kMagenta
*/
hFAKES->SetFillColor(kMagenta-10); hFAKES->SetLineColor(kMagenta+4); hFAKES->SetLineWidth(1);
hZTauTau->SetFillColor(kOrange-4); hZTauTau->SetLineColor(kOrange+4); hZTauTau->SetLineWidth(1);
hDY->SetFillColor(kAzure+3); hDY->SetLineColor(kAzure+4); hDY->SetLineWidth(1);
hWW->SetFillColor(kAzure+3); hWW->SetLineColor(kAzure+3); hWW->SetLineWidth(1);
hTOP->SetFillColor(kGreen-2); hTOP->SetLineColor(kGreen+4); hTOP->SetLineWidth(1);
hTT->SetFillColor(kGreen-2); hTT->SetLineColor(kGreen-2); hTT->SetLineWidth(1);
hLFVGG126->SetLineColor(kBlue+1); hLFVGG126->SetLineWidth(3);
hLFVVBF126->SetLineColor(kBlue+1); hLFVVBF126->SetLineWidth(3); hLFVVBF126->SetLineStyle(kDashed);
hSMVBF126->SetLineColor(kMagenta); hSMVBF126->SetLineWidth(3); hSMVBF126->SetLineStyle(kDashed);
hSMGG126->SetLineColor(kMagenta); hSMGG126->SetLineWidth(3);
hSMHIGGS->SetFillColor(kMagenta); hSMHIGGS->SetLineColor(kMagenta+1); hSMHIGGS->SetLineWidth(1);
hdata_obs->SetMarkerSize(1); // Closer to Daniel's
// Rebin
hFAKES->Rebin(rebinning);
hZTauTau->Rebin(rebinning);
hDY->Rebin(rebinning);
hTOP->Rebin(rebinning);
hTT->Rebin(rebinning);
hWW->Rebin(rebinning);
hLFVGG126->Rebin(rebinning);
hLFVVBF126->Rebin(rebinning);
hSMVBF126->Rebin(rebinning);
hSMGG126->Rebin(rebinning);
hdata_obs->Rebin(rebinning);
hSMHIGGS->Rebin(rebinning);
// PLOT
TCanvas *c1 = new TCanvas("canvas_"+name);
TPad *Pad1= new TPad("pad1","",0,0.2,1,1); Pad1->Draw(); Pad1->cd();;
Pad1->SetLeftMargin(0.2147651);
Pad1->SetRightMargin(0.06543624);
Pad1->SetTopMargin(0.07);
Pad1->SetBottomMargin(0.04);
for (int i=0; i<hFAKES->GetNbinsX()+1; i++){
double content=hFAKES->GetBinContent(i);
double contentErr=0;//sqrt(content);//hFAKES->GetBinError(i);
double err=sqrt(0.3*0.3*content*content+contentErr*contentErr);
hFAKES->SetBinError(i,err);
}
TH1F* fullMC2=hZTauTau->Clone(); fullMC2->Add(hFAKES); fullMC2->Add(hTT); fullMC2->Add(hWW);
fullMC2->Add(hDY); fullMC2->Add(hTOP);
fullMC2->SetFillColorAlpha(kGray+2, 0.35);
//fullMC2->SetFillStyle(3002);
fullMC2->SetMarkerSize(0);
fullMC2->Draw("hist");
fullMC2->GetXaxis()->SetTitle("");
fullMC2->GetYaxis()->SetTitle(Yaxis);
fullMC2->GetXaxis()->SetRangeUser(xmin,xmax);
fullMC2->GetYaxis()->SetTitleOffset(1.2);
fullMC2->GetYaxis()->SetTitleSize(0.05);
fullMC2->GetXaxis()->SetNdivisions(505);
fullMC2->GetYaxis()->SetLabelSize(0.04);
fullMC2->GetXaxis()->SetLabelSize(0);
THStack* stack = new THStack("stack","");
stack->Add(hFAKES);
stack->Add(hWW);
stack->Add(hDY);
stack->Add(hTT);
stack->Add(hTOP);
stack->Add(hZTauTau);
stack->Add(hSMHIGGS);
cout<<"Yields"<<endl;
cout<<"VV " <<hWW->Integral()<<endl;
cout<<"TOP " <<hTOP->Integral()<<endl;
cout<<"TT " <<hTT->Integral()<<endl;
cout<<"DY " <<hDY->Integral()<<endl;
cout<<"ZTauTau " <<hZTauTau->Integral()<<endl;
cout<<"FAKES " <<hFAKES->Integral()<<endl;
cout<<"LFVGG126 " <<hLFVGG126->Integral()<<endl;
cout<<"LFVVBF126 " <<hLFVVBF126->Integral()<<endl;
stack->Draw("samehist");
fullMC2->Draw("sames,E2");
hLFVVBF126->Draw("sameshist");
hLFVGG126->Draw("sameshist");
double maxData=hdata_obs->GetMaximum();
double maxMC=stack->GetMaximum()*1.2;
double maxLFV=hLFVGG126->GetMaximum();
double maxLFV2=hLFVVBF126->GetMaximum();
double minMC=stack->GetMinimum();
if(maxData>maxMC) {maxMC=1.2*maxData;}
if(maxLFV>maxMC) {maxMC=1.2*maxLFV;}
if(maxLFV2>maxMC) {maxMC=1.2*maxLFV2;}
if(MAX!=-1) {maxMC=MAX;}
stack->SetMaximum(maxMC);
stack->GetYaxis()->SetRangeUser(minMC,maxMC);
fullMC2->SetMaximum(maxMC);
fullMC2->GetYaxis()->SetRangeUser(minMC,maxMC);
if(setLogY){fullMC2->SetMinimum(0.1); fullMC2->GetYaxis()->SetRangeUser(1,maxMC*1000);}
TH1F* histoDataUnblindedV4=hdata_obs->Clone(); histoDataUnblindedV4->SetName("data_unblinded");
if(!blind) hdata_obs->Draw("sames");
else {
int findBinA=hdata_obs->FindBin(blindA);
int findBinB=hdata_obs->FindBin(blindB);
for (int i=findBinA; i<findBinB; i++) hdata_obs->SetBinContent(i,-1000);
hdata_obs->Draw("sames");
TPave *pave = new TPave(blindA,minMC,blindB,maxMC,4,"br");
//pave->SetFillColorAlpha(kGray+3,0.35);
pave->SetFillColor(1);
pave->SetFillStyle(3003);
pave->SetDrawOption(0);
pave->SetBorderSize(0);
pave->Draw();
}
TLegend *leg = new TLegend(legx1,legy1,legx2,legy2,NULL,"brNDC");
leg->SetFillColor(0);
leg->SetBorderSize(0);
leg->SetFillStyle(0);
TLegendEntry *entry;
entry=leg->AddEntry(hdata_obs,"Data, e#tau_{h}","p");
entry=leg->AddEntry(fullMC2,"Bkgd. uncertainty","f");
entry=leg->AddEntry(hSMHIGGS,"SM Higgs","f");
eblindBy=leg->AddEntry(hZTauTau,"Z#rightarrow#tau#tau","f");
entry=leg->AddEntry(hTOP,"t#bar{t}, t, #bar{t}","f");
entry=leg->AddEntry(hDY,"Other","f");
entry=leg->AddEntry(hFAKES,"MisID'd #tau","f");
entry=leg->AddEntry(hLFVGG126,"LFV GG Higgs (B=100%)","l");
entry=leg->AddEntry(hLFVVBF126,"LFV VBF Higgs (B=100%)","l");
leg->Draw();
CMS_lumi( Pad1, iPeriod, iPos );
TLatex latex;
latex.SetNDC();
latex.SetTextAngle(0);
latex.SetTextColor(kBlack);
latex.SetTextFont(42);
latex.SetTextAlign(31);
latex.SetTextSize(.05);
latex.DrawLatex(0.39,0.8,Tex);
//cmsPrelim(Lumi);
Pad1->SetLogy(setLogY);
gPad->RedrawAxis();
c1->cd(); TPad *Pad2= new TPad("pad2","",0,0,1,0.23); Pad2->Draw(); Pad2->cd(); Pad2->SetGridy();
Pad2->SetLeftMargin(0.2147651);
Pad2->SetRightMargin(0.06543624);
Pad2->SetTopMargin(0.0);
Pad2->SetBottomMargin(0.37);
Pad2->SetFillStyle(0);
TH1F * Ratio=fullMC2->Clone(); Ratio->SetName("Ratio");
Ratio->Scale(-1);
Ratio->Add(hdata_obs,1);
Ratio->Divide(fullMC2);
Ratio->SetLineColor(kBlack);
Ratio->GetXaxis()->SetLabelFont(42);
Ratio->GetXaxis()->SetTitleFont(42);
Ratio->GetYaxis()->SetNdivisions(505);
Ratio->GetYaxis()->SetLabelFont(42);
Ratio->GetYaxis()->SetLabelSize(0.122);
Ratio->GetYaxis()->SetRangeUser(ymin,ymax);
Ratio->GetXaxis()->SetRangeUser(xmin,xmax);
Ratio->GetXaxis()->SetLabelSize(0.12);
Ratio->GetXaxis()->SetLabelFont(42);
Ratio->SetYTitle("#frac{Data-Bkgd}{Bkgd}");
Ratio->SetXTitle(Xaxis);
Ratio->GetXaxis()->SetNdivisions(505);
Ratio->GetYaxis()->CenterTitle(true);
Ratio->GetYaxis()->SetTitleOffset(0.4);
Ratio->GetYaxis()->SetTitleSize(0.11);
Ratio->GetXaxis()->SetTitleOffset(0.75);
Ratio->GetXaxis()->SetTitleSize(0.20);
Ratio->SetMarkerSize(1.);
TH1F* RatioError = Ratio->Clone(); RatioError->SetName("RatioError");
for (int i=0; i<RatioError->GetNbinsX()+1; i++){
double error=fullMC2->GetBinError(i)*hdata_obs->GetBinContent(i)/fullMC2->GetBinContent(i)/fullMC2->GetBinContent(i);
RatioError->SetBinContent(i,0);
RatioError->SetBinError(i,error);
}
//RatioError->SetFillStyle(3002);
RatioError->SetFillColorAlpha(kGray+2,0.35); RatioError->SetMarkerSize(0);
RatioError->Draw("E2");
//Ratio->Draw("samesE0"); //errorbarswithoutdatapoints
Ratio->Draw("sames");
if(!blind) Ratio->Draw("sames");
else {
int findBinA=Ratio->FindBin(blindA);
int findBinB=Ratio->FindBin(blindB);
for (int i=findBinA; i<findBinB; i++) Ratio->SetBinContent(i,-100);
Ratio->Draw("sames");
TPave *pave = new TPave(blindA,-1.0,blindB,1.0,4,"br");
pave->SetFillColor(1);
pave->SetFillStyle(3003);
//pave->SetFillColorAlpha(kGray+3,0.35);
pave->SetDrawOption(0);
pave->SetBorderSize(0);
pave->Draw();
}
if(!setLogY){
// c1->SaveAs(name+"ErrorBarsWithoutDataPoints.png");
// c1->SaveAs(name+"ErrorBarsWithoutDataPoints.pdf");
c1->SaveAs(name+".png");
c1->SaveAs(name+".pdf");
}
else {
c1->SaveAs(name+"_log.png");
c1->SaveAs(name+"_log.pdf");
}
}
void file(){
TCanvas *c1 = new TCanvas("c1","ROOT File description",200,10,700,550);
c1->Range(0,-0.25,21,14);
TPaveLabel *title = new TPaveLabel(5,12,15,13.7,c1->GetTitle());
title->SetFillColor(16);
title->Draw();
// horizonthal file layout
TPave *file = new TPave(1,8.5,20,11);
file->SetFillColor(11);
file->Draw();
TPave *fileh = new TPave(1,8.5,2.5,11);
fileh->SetFillColor(44);
fileh->Draw();
TPave *lrh = new TPave(2.5,8.5,3.3,11,1);
lrh->SetFillColor(33);
lrh->Draw();
lrh->DrawPave(6.9,8.5,7.7,11,1);
lrh->DrawPave(10.5,8.5,11.3,11,1);
lrh->DrawPave(14.5,8.5,15.3,11,1);
TLine *ldot = new TLine(1,8.5,0.5,6.5);
ldot->SetLineStyle(2);
ldot->Draw();
ldot->DrawLine(2.5, 8.5, 9.4, 6.5);
ldot->DrawLine(10.5, 8.5, 10, 6.5);
ldot->DrawLine(11.3, 8.5, 19.5, 6.5);
TLine *line = new TLine(2.6,11,2.6,11.5);
line->Draw();
line->DrawLine(2.6,11.5,7,11.5);
TArrow *arrow = new TArrow(7,11.5,7,11.1,0.01,"|>");
arrow->SetFillStyle(1001);
arrow->Draw();
line->DrawLine( 7, 8.5, 7, 8.0);
line->DrawLine( 7, 8.0, 10.6, 8);
arrow->DrawArrow( 10.6,8, 10.6, 8.4,0.01,"|>");
line->DrawLine( 10.6, 11, 10.6, 11.5);
line->DrawLine( 10.6, 11.5, 14.6, 11.5);
arrow->DrawArrow( 14.6,11.5, 14.6,11.1,0.01,"|>");
line->DrawLine( 14.6, 8.5, 14.6, 8.0);
line->DrawLine( 14.6, 8.0, 16, 8);
ldot->DrawLine(16, 8, 19, 8);
TText *vert = new TText(1.5,9.75,"File");
vert->SetTextAlign(21);
vert->SetTextAngle(90);
vert->SetTextSize(0.025);
vert->Draw();
vert->DrawText(2.0, 9.75,"Header");
vert->DrawText(2.9, 9.75,"Logical Record");
vert->DrawText(3.2, 9.75,"Header");
vert->DrawText(7.3, 9.75,"Logical Record");
vert->DrawText(7.6, 9.75,"Header");
vert->DrawText(10.9,9.75,"Logical Record");
vert->DrawText(11.2,9.75,"Header");
vert->DrawText(14.9,9.75,"Logical Record");
vert->DrawText(15.2,9.75,"Header");
TText *hori = new TText(4.75,10,"Object");
hori->SetTextAlign(22);
hori->SetTextSize(0.035);
hori->Draw();
hori->DrawText(4.75, 9.5,"Data");
hori->DrawText(9.2, 10,"Deleted");
hori->DrawText(9.2, 9.5,"Object");
line->DrawLine( 6.9, 8.5, 10.5, 11);
line->DrawLine( 6.9, 11, 10.5, 8.5);
TText *tbig = new TText(17,9.75,"............");
tbig->SetTextAlign(22);
tbig->SetTextSize(0.03);
tbig->Draw();
tbig->DrawText(2.6, 7, "fBEGIN");
tbig->DrawText(20., 7, "fEND");
arrow->DrawArrow( 2.6,7, 2.6,8.4,0.01,"|>");
arrow->DrawArrow( 20,7, 20,8.4,0.01,"|>");
//file header
TPaveText *header = new TPaveText(0.5,.2,9.4,6.5);
header->SetFillColor(44);
header->Draw();
TText *fh=header->AddText("File Header");
fh->SetTextAlign(22);
fh->SetTextSize(0.04);
header->SetTextSize(0.027);
header->SetTextAlign(12);
header->AddText(" ");
header->AddLine(0,0,0,0);
header->AddText("\"root\": Root File Identifier");
header->AddText("fVersion: File version identifier");
header->AddText("fBEGIN: Pointer to first data record");
header->AddText("fEND: Pointer to first free word at EOF");
header->AddText("fSeekFree: Pointer to FREE data record");
header->AddText("fNbytesFree: Number of bytes in FREE");
header->AddText("fNfree: Number of free data records");
header->AddText("fNbytesName: Number of bytes in name/title");
header->AddText("fUnits: Number of bytes for pointers");
header->AddText("fCompress: Compression level");
//logical record header
TPaveText *lrecord = new TPaveText(10,0.2,19.5,6.5);
lrecord->SetFillColor(33);
lrecord->Draw();
TText *tlrh=lrecord->AddText("Logical Record Header (TKEY)");
tlrh->SetTextAlign(22);
tlrh->SetTextSize(0.04);
lrecord->SetTextSize(0.027);
lrecord->SetTextAlign(12);
lrecord->AddText(" ");
lrecord->AddLine(0,0,0,0);
lrecord->AddText("fNbytes: Length of compressed object");
lrecord->AddText("fVersion: Key version identifier");
lrecord->AddText("fObjLen: Length of uncompressed object");
lrecord->AddText("fDatime: Date/Time when written to store");
lrecord->AddText("fKeylen: Number of bytes for the key");
lrecord->AddText("fCycle : Cycle number");
lrecord->AddText("fSeekKey: Pointer to object on file");
lrecord->AddText("fSeekPdir: Pointer to directory on file");
lrecord->AddText("fClassName: class name of the object");
lrecord->AddText("fName: name of the object");
lrecord->AddText("fTitle: title of the object");
c1->Update();
c1->Print("file.png");
}
void HistoData::drawResult(TH1 *Summary, bool Vertical, bool SetBinLabel) {
// add label to the summary if desired
if (SetBinLabel) {
Summary->GetXaxis()->SetBinLabel(bin,getRefHisto()->GetTitle());
//Summary->GetXaxis()->SetBinLabel(bin,name.c_str());
}
else
Summary->GetXaxis()->SetBinLabel(bin,name.c_str());
double minimum = Summary->GetMinimum();
// determine where to draw the result (score axis)
// 1: solid bar starts
// 2: solid bar ends, hatched bar starts
// 3: hatched bar ends
double score1 = minimum;
double score2 = (lowScore == 10. || lowScore < minimum) ? minimum : lowScore;
double score3 = (lowScore == 10.) ? 1 : ((highScore < minimum) ? minimum : highScore);
// determine where to draw the result (binning axis)
double binCenter = Summary->GetBinCenter(bin);
double binWidth = Summary->GetBinWidth(bin);
double bin1 = binCenter - binWidth / 3;
double bin2 = binCenter + binWidth / 3;
// set coordinates of display based on plot alignment
double solidX1, solidY1, solidX2, solidY2;
double hatchedX1, hatchedY1, hatchedX2, hatchedY2;
double axisX1, axisY1, axisX2, axisY2;
if (Vertical) {
solidX1 = bin1; solidX2 = bin2;
solidY1 = score1; solidY2 = score2;
hatchedX1 = bin1; hatchedX2 = bin2;
hatchedY1 = score2; hatchedY2 = score3;
axisX1 = bin1; axisX2 = bin2;
axisY1 = minimum; axisY2 = minimum;
} else {
solidX1 = score1; solidX2 = score2;
solidY1 = bin1; solidY2 = bin2;
hatchedX1 = score2; hatchedX2 = score3;
hatchedY1 = bin1; hatchedY2 = bin2;
axisX1 = minimum; axisX2 = minimum;
axisY1 = bin1; axisY2 = bin2;
}
// a solid bar is drawn from zero to the lowest score
if (lowScore > minimum && lowScore != 10.) {
TPave *solidBar = new TPave(solidX1,solidY1,solidX2,solidY2,1,"");
solidBar->SetBit(kCanDelete);
solidBar->SetLineColor(lineUseFillColor ? solidFillColor : solidLineColor);
solidBar->SetFillColor(solidFillColor);
solidBar->SetFillStyle(solidFillStyle);
solidBar->Draw();
}
// a hatched bar is drawn from the lowest score to the highest score
if ((lowScore != highScore && highScore > minimum) || lowScore == 10.) {
TPave *hatchedBar = new TPave(hatchedX1,hatchedY1,hatchedX2,hatchedY2,1,"");
hatchedBar->SetBit(kCanDelete);
hatchedBar->SetLineColor(lineUseFillColor ? shadedFillColor : shadedLineColor);
hatchedBar->SetFillColor(shadedFillColor);
hatchedBar->SetFillStyle(3004);
hatchedBar->Draw();
}
// paste a line over the base axis to fix border color artifacts
TLine *axisLine = new TLine(axisX1,axisY1,axisX2,axisY2);
if (Vertical) axisLine->SetLineColor(Summary->GetAxisColor("X"));
else axisLine->SetLineColor(Summary->GetAxisColor("Y"));
axisLine->SetBit(kCanDelete);
axisLine->Draw("SAME");
//paste a line before (in this proceeding that means after) Barrel and Endcap Plots
if ( name == "ERPt" || name == "BRPt" ) {
axisY1=axisY2=binCenter+binWidth/2;
axisX2=Summary->GetMaximum();
axisX1=Summary->GetMinimum()-200;
TLine *regionLine = new TLine(axisX1,axisY1,axisX2,axisY2);
regionLine->SetLineColor(Summary->GetAxisColor("X"));
regionLine->SetBit(kCanDelete);
regionLine->SetLineWidth(3);
regionLine->Draw();
}
}
