void highlight1()
{
TCanvas *ch = new TCanvas("ch", "ch", 0, 0, 700, 500);
const Int_t n = 500;
Double_t x[n], y[n];
TH1F *h[n];
for (Int_t i = 0; i < n; i++) {
h[i] = new TH1F(TString::Format("h_%02d", i), "", 100, -3.0, 3.0);
// in practice gaus need reset parameters
h[i]->FillRandom("gaus", 1000);
h[i]->Fit("gaus", "Q");
h[i]->SetMaximum(250); // for n > 200
x[i] = i;
y[i] = h[i]->GetFunction("gaus")->GetParameter(2);
}
TGraph *g = new TGraph(n, x, y);
g->SetMarkerStyle(6);
for (Int_t i = 0; i < n; i++) g->AddHighlight(i, h[i]);
g->Draw("AP");
g->SetHighlight();
TPad *ph = new TPad("ph", "ph", 0.3, 0.4, 1.0, 1.0);
ph->SetFillColor(kBlue - 10);
ph->Draw();
ph->cd();
TText *info = new TText(0.5, 0.5, "please move the mouse over the graph");
info->SetTextAlign(22);
info->Draw();
ch->cd();
TGraph::SetHighlightPad(ph);
}
//--------------------------------------------------------------------------------------------------
void overlayFrame(TString text, bool align)
{
// Overlay a linear frame from user coordinates (0 - 1, 0 - 1) and put the frame text
// Create new transparent pad for the text
TPad *transPad = new TPad("transPad","Transparent Pad",0,0,1,1);
transPad->SetFillStyle(4000);
transPad->Draw();
transPad->cd();
// Overlay the text in a well defined frame
TText *plotText = new TText();
plotText->SetTextColor(kBlue);
plotText->SetTextSize(0.04);
plotText->SetNDC();
// Draw text at top right
if (align) {
plotText->SetTextColor(kBlack);
plotText->SetTextAlign(33);
plotText->DrawText(0.92,0.95,text.Data());
}
// Draw text at bottom left
else
plotText->DrawText(0.01,0.01,text.Data());
return;
}
void greyscale()
{
TCanvas *c = new TCanvas("grey", "Grey Scale", 500, 500);
c->SetBorderMode(0);
Int_t n = 200; // tunable parameter
Float_t n1 = 1./n;
for (int i = 0; i < n; i++) {
for (int j = 0; j < n; j++) {
TBox *b = new TBox(n1*j, n1*(n-1-i), n1*(j+1), n1*(n-i));
Float_t grey = Float_t(i*n+j)/(n*n);
b->SetFillColor(TColor::GetColor(grey, grey, grey));
b->Draw();
}
}
TPad *p = new TPad("p","p",0.3, 0.3, 0.7,0.7);
const char *guibackground = gEnv->GetValue("Gui.BackgroundColor", "");
p->SetFillColor(TColor::GetColor(guibackground));
p->Draw();
p->cd();
TText *t = new TText(0.5, 0.5, "GUI Background Color");
t->SetTextAlign(22);
t->SetTextSize(.09);
t->Draw();
c->SetEditable(kFALSE);
}
void PaintOverflow(TH1 *h)
{
// This function paint the histogram h with an extra bin for overflows
char* name = h->GetName();
char* title = h->GetTitle();
Int_t nx = h->GetNbinsX()+1;
Double_t x1 = h->GetBinLowEdge(1);
Double_t bw = h->GetBinWidth(nx);
Double_t x2 = h->GetBinLowEdge(nx)+bw;
// Book a temporary histogram having ab extra bin for overflows
TH1F *htmp = new TH1F(name, title, nx, x1, x2);
// Fill the new hitogram including the extra bin for overflows
for (Int_t i=1; i<=nx; i++) {
htmp->Fill(htmp->GetBinCenter(i), h->GetBinContent(i));
}
// Fill the underflows
htmp->Fill(x1-1, h->GetBinContent(0));
// Restore the number of entries
htmp->SetEntries(h->GetEntries());
// Draw the temporary histogram
htmp->Draw();
TText *t = new TText(x2-bw/2,h->GetBinContent(nx),"Overflow");
t->SetTextAngle(90);
t->SetTextAlign(12);
t->SetTextSize(0.03);;
t->Draw();
}
void table(Float_t x1, Float_t x2, Float_t yrange, TText *t,
const char **symbol, Bool_t octal)
{
Int_t i;
Int_t n = 0;
for (i=0;i<1000;i++) {
if (!strcmp(symbol[i],"END")) break;
n++;
}
Float_t y1 = 2.5;
Float_t y2 = yrange - 0.5;
Float_t dx = (x2-x1)/5;
Float_t dy = (y2 - 1 -y1)/(n+1);
Float_t y = y2 - 1 - 0.7*dy;
Float_t xc0 = x1 + 0.5*dx;
Float_t xc1 = xc0 + dx;
Float_t xc2 = xc1 + dx;
Float_t xc3 = xc2 + dx;
Float_t xc4 = xc3 + dx;
TLine *line = new TLine();
line->DrawLine(x1,y1,x1,y2);
line->DrawLine(x1,y1,x2,y1);
line->DrawLine(x1,y2,x2,y2);
line->DrawLine(x2,y1,x2,y2);
line->DrawLine(x1,y2-1,x2,y2-1);
line->DrawLine(x1+ dx,y1,x1+ dx,y2);
line->DrawLine(x1+2*dx,y1,x1+2*dx,y2);
line->DrawLine(x1+3*dx,y1,x1+3*dx,y2);
line->DrawLine(x1+4*dx,y1,x1+4*dx,y2);
TText *tit = new TText(0,0,"a");
tit->SetTextSize(0.015);
tit->SetTextFont(72);
tit->SetTextAlign(22);
tit->DrawText(xc0,y2-0.6,"Input");
tit->DrawText(xc1,y2-0.6,"Roman");
tit->DrawText(xc2,y2-0.6,"Greek");
tit->DrawText(xc3,y2-0.6,"Special");
tit->DrawText(xc4,y2-0.6,"Zapf");
char text[12];
for (i=0;i<n;i++) {
if (octal) {
unsigned char value = *symbol[i];
sprintf(text,"@\\ %3o",value);
} else {
strcpy(text,symbol[i]);
}
t->DrawText(xc0,y,text);
sprintf(text,"%s",symbol[i]);
t->DrawText(xc1,y,text);
sprintf(text,"`%s",symbol[i]);
t->DrawText(xc2,y,text);
sprintf(text,"'%s",symbol[i]);
t->DrawText(xc3,y,text);
sprintf(text,"~%s",symbol[i]);
t->DrawText(xc4,y,text);
y -= dy;
}
}
void SetupTowerDisplay(TH2F *hist)
{
hist->SetStats(kFALSE);
hist->SetXTitle("ieta");
hist->SetYTitle("iphi");
hist->GetXaxis()->CenterTitle();
hist->GetYaxis()->CenterTitle();
hist->GetXaxis()->SetNdivisions(65);
hist->GetXaxis()->SetLabelColor(0);
hist->GetXaxis()->SetTickLength(.78);
hist->GetXaxis()->SetTitleOffset(0.95);
hist->GetYaxis()->SetNdivisions(72);
hist->GetYaxis()->SetLabelColor(0);
hist->GetYaxis()->SetTitleOffset(0.85);
TText *yLabel = new TText();
TLine *pLine = new TLine();
pLine->SetLineStyle(1);
pLine->SetLineColor(1);
pLine->SetLineWidth(1);
yLabel->SetTextAlign(22);
yLabel->SetTextSize(0.015);
char phi_num[3];
char eta_num[3];
TText *xLabel = new TText();
xLabel->SetTextSize(0.015);
xLabel->SetTextAlign(22);
for (Int_t i=1; i<73; ++i)
{
sprintf(phi_num,"%d",i);
if(TMath::Abs(i%2)==1) {yLabel->DrawText(-33,0.5+i,phi_num);}
else {yLabel->DrawText(-34.5,0.5+i,phi_num);}
pLine->DrawLine(-32,i,33,i);
}
for (Int_t i=-32; i<33;++i)
{
sprintf(eta_num,"%d",i);
if(TMath::Abs(i%2)==0) {xLabel->DrawText(0.5+i,-0.5,eta_num);}
else {xLabel->DrawText(0.5+i,-2,eta_num);}
pLine->DrawLine(i,1,i,72);
}
}
//______________________________________________________________________________
void DrawPeriod(int runmin, int runmax, int run1, int run2, double ymin, double ymax, const char* label)
{
if ( run1 < runmin || run1 > runmax || run2 < runmin || run2 > runmax ) return;
TBox* b = new TBox(run1,ymin,run2,ymax);
b->SetFillColor(5);
b->Draw();
TText* text = new TText((run1+run2)/2.0,ymax*0.6,label);
text->SetTextAlign(22);
text->SetTextSize(0.02);
text->Draw();
}
void fill_with_text(TGraph *real, TGraph *down, TGraph *up, TVirtualPad *can, int scantype, std::string scanx = "") {
can->cd();
float xpos_of_text = 0.22;
TLegend* this_leg = new TLegend(xpos_of_text,0.6,0.38,0.75);
//TLegend* this_leg = new TLegend(xpos_of_text,0.55,0.45,0.75,"n_{jets} #geq 3"); // this was the style of the paper.
this_leg->SetFillColor(0);
this_leg->SetBorderSize(0);
this_leg->SetTextSize(0.035);
//this_leg->SetTextSize(0.04); // paper style.
if(scantype==PlottingSetup::SMS||scantype==PlottingSetup::GMSB) {
//this_leg->AddEntry(real,"#sigma^{prod} = #sigma^{NLO-QCD}" , "l");
//this_leg->AddEntry(up,"#sigma^{prod} = 3 #times #sigma^{NLO-QCD}" , "l");
//this_leg->AddEntry(down,"#sigma^{prod} = 1/3 #times #sigma^{NLO-QCD}" , "l");
this_leg->AddEntry(real,"#sigma^{NLO-QCD}" , "l");
this_leg->AddEntry(up,"3 #times #sigma^{NLO-QCD}" , "l");
this_leg->AddEntry(down,"1/3 #times #sigma^{NLO-QCD}" , "l");
} else {
write_warning(__FUNCTION__,"Not implemented yet for mSUGRA");
}
this_leg->Draw();
TText *title = write_text(xpos_of_text+0.005,0.52,"JZB");
title->SetTextSize(0.04);
title->SetTextAlign(13);
title->SetTextFont(62);
title->Draw();
write_SMS_text( scantype, scanx, xpos_of_text );
// //string legT5zz="pp #rightarrow #tilde{g} #tilde{g}, #tilde{g} #rightarrow 2j + #chi^{0}_{1}, #chi^{0}_{1} #rightarrow Z + #tilde{G}";
// string legT5zz="pp #rightarrow #tilde{g} #tilde{g}, #tilde{g} #rightarrow 2j + #chi^{0}_{2}, #chi^{0}_{2} #rightarrow Z #chi^{0}_{1}";
// string legT5zzl2="m(#tilde{q}) >> m(#tilde{g}), x = "+scanx;
// string legT5zzl3=" GMSB";
// TText *title = write_text(xpos_of_text,0.85,legT5zz);
// title->SetTextAlign(11);
// title->SetTextSize(0.035);
// if(scantype!=PlottingSetup::mSUGRA) title->Draw("same");
// TText *title2 = write_text(xpos_of_text,0.79,legT5zzl2);
// title2->SetTextAlign(11);
// title2->SetTextSize(0.035);
// if(scantype!=PlottingSetup::mSUGRA) title2->Draw("same");
// TText *title3 = write_text(0.40,0.79,legT5zzl3);
// title3->SetTextAlign(11);
// title3->SetTextSize(0.035);
// title3->SetTextColor(kRed);
//// if(scantype==PlottingSetup::GMSB) title3->Draw("same");
DrawPrelim();
draw_diagonal_xchange( scantype, scanx );
}
//______________________________________________________________________________
void AddText(TPaveText *pave, const char *datamember, Int_t value, const char *comment)
{
char line[128];
for (Int_t i=0; i<128; i++) line[i] = ' ';
memcpy(&line[0], datamember, strlen(datamember));
line[10] = '=';
char number[20];
sprintf(number, "%5i", value);
memcpy(&line[12], number, strlen(number));
line[26] = '=';
line[27] = '>';
sprintf(&line[30], "%s",comment);
TText *text = pave->AddText(line);
text->SetTextAlign(12);
}
//______________________________________________________________________________
void AddText(TPaveText *pave, TObject *pf, Int_t iaxis)
{
char line[128];
TGeoPatternFinder *finder = (TGeoPatternFinder*)pf;
if (!pave || !pf) return;
for (Int_t i=0; i<128; i++) line[i] = ' ';
TGeoVolume *volume = finder->GetVolume();
TGeoShape *sh = volume->GetShape();
sprintf(line, "Division of %s on axis %d (%s)", volume->GetName(), iaxis,sh->GetAxisName(iaxis));
TText *text = pave->AddText(line);
text->SetTextColor(3);
text->SetTextAlign(12);
AddText(pave, "fNdiv",finder->GetNdiv(),"number of divisions");
AddText(pave, "fStart",finder->GetStart(),"start divisioning position");
AddText(pave, "fStep",finder->GetStep(),"division step");
}
TCanvas* pHitSpecDet(const std::string& s)
{
TCanvas* c = new TCanvas( ("HitSpecDet"+s).c_str(), ("HitSpecDet"+s).c_str(), -2);
c->SetLogx();
TH2D* h = (TH2D*)gROOT->FindObject(s.c_str());
h->GetXaxis()->SetRange(4,32);
std::stringstream str;
str<<h->GetTitle();
TText t;
t.SetTextColor(4);
t.SetTextAlign(11);
t.SetTextSize(0.035);
h->DrawCopy("box");
t.DrawTextNDC(0.17,0.2, str.str().c_str());
return c;
}
void write_SMS_text(int scantype, std::string& scanx, float xpos = 0.22 ) {
string legT5zz="pp #rightarrow #tilde{g} #tilde{g}, #tilde{g} #rightarrow 2j + #chi^{0}_{2}, #chi^{0}_{2} #rightarrow Z #chi^{0}_{1}";
if ( scantype==PlottingSetup::GMSB) legT5zz="pp #rightarrow #tilde{g} #tilde{g}, #tilde{g} #rightarrow 2j + #chi^{0}_{1}, #chi^{0}_{1} #rightarrow Z #tilde{G}";
string legT5zzl2="m(#tilde{q}) >> m(#tilde{g})";
if (scantype==PlottingSetup::SMS) legT5zzl2 += ", x = "+scanx;
string legT5zzl3=" GMSB";
TText *title = write_text(xpos,0.85,legT5zz);
title->SetTextAlign(11);
title->SetTextSize(0.035);
if(scantype!=PlottingSetup::mSUGRA) title->Draw("same");
TText *title2 = write_text(xpos,0.79,legT5zzl2);
title2->SetTextAlign(11);
title2->SetTextSize(0.035);
if(scantype!=PlottingSetup::mSUGRA) title2->Draw("same");
TText *title3 = write_text(0.40,0.79,legT5zzl3);
title3->SetTextAlign(11);
title3->SetTextSize(0.035);
title3->SetTextColor(kRed);
// if(scantype==PlottingSetup::GMSB) title3->Draw("same");
}
void plot_example(){
THStack *hs = new THStack("hs","");
TH1F *h1 = new TH1F("h1","test hstack",10,-4,4);
h1->FillRandom("gaus",20000);
h1->SetFillColor(kRed);
hs->Add(h1);
TH1F *h2 = new TH1F("h2","test hstack",10,-4,4);
h2->FillRandom("gaus",15000);
h2->SetFillColor(kBlue);
hs->Add(h2);
TH1F *h3 = new TH1F("h3","test hstack",10,-4,4);
h3->FillRandom("gaus",10000);
h3->SetFillColor(kGreen);
hs->Add(h3);
TCanvas *cs = new TCanvas("cs","cs",10,10,700,900);
TText T; T.SetTextFont(42); T.SetTextAlign(21);
cs->Divide(2,2);
cs->cd(1); hs->Draw(); T.DrawTextNDC(.5,.95,"Default drawing option");
cs->cd(2); hs->Draw("nostack"); T.DrawTextNDC(.5,.95,"Option \"nostack\"");
cs->cd(3); hs->Draw("nostackb"); T.DrawTextNDC(.5,.95,"Option \"nostackb\"");
cs->cd(4); hs->Draw("lego1"); T.DrawTextNDC(.5,.95,"Option \"lego1\"");
}
void hlHisto4()
{
TCanvas *c1 = new TCanvas("c1", "", 0, 0, 600, 400);
TF1 *f1 = new TF1("f1", "x*gaus(0) + [3]*abs(sin(x)/x)", -50.0, 50.0);
f1->SetParameters(20.0, 4.0, 1.0, 20.0);
TH1F *h1 = new TH1F("h1", "Test random numbers", 200, -50.0, 50.0);
h1->FillRandom("f1", 100000);
h1->Draw();
h1->Fit(f1, "Q");
gStyle->SetGridColor(kGray);
c1->SetGrid();
TText *info = new TText(0.0, h1->GetMaximum()*0.7, "please move the mouse over the frame");
info->SetTextSize(0.04);
info->SetTextAlign(22);
info->SetTextColor(kRed-1);
info->SetBit(kCannotPick);
info->Draw();
c1->Update();
h1->SetHighlight();
c1->HighlightConnect("HighlightZoom(TVirtualPad*,TObject*,Int_t,Int_t)");
}
void pstable()
{
// Display all possible types of ROOT/Postscript characters
const char *symbol1[] =
{"A","B","C","D","E","F","G","H","I","J","K","L","M","N",
"O","P","Q","R","S","T","U","V","W","X","Y","Z",
"0","1","2","3","4","5","6","7","8","9",
".",",","+","-","*","/","=","(",")","{","}","END"};
const char *symbol2[] =
{"a","b","c","d","e","f","g","h","i","j","k","l","m","n",
"o","p","q","r","s","t","u","v","w","x","y","z",
":",";","@","\\","_","|","%",
"@'","<",">","[","]","\42","@\43","@\136",
"@\77","@\41","@&","$","@\176"," ","END"};
const char *symbol3[] =
{"\241","\242","\243","\244","\245","\246","\247","\250",
"\251","\252","\253","\254","\255","\256","\257","\260",
"\261","\262","\263","\264","\265","\266","\267","\270",
"\271","\272","\273","\274","\275","\276","\277","\300",
"\301","\302","\303","\304","\305","\306","\307","\310",
"\311","\312","\313","\314","\315","\316","\317","END"};
const char *symbol4[] =
{"\321","\322","\323","\324","\325","\326","\327","\330",
"\331","\332","\333","\334","\335","\336","\337","\340",
"\341","\342","\343","\344","\345","\346","\347","\340",
"\351","\352","\353","\354","\355","\356","\357","\360",
"\361","\362","\363","\364","\365","\366","\367","\370",
"\371","\372","\373","\374","\375","\376","\377","END"};
const char *symbol5[] =
{"\177","\200","\201","\202","\203","\204","\205","\206",
"\207","\210","\211","\212","\213","\214","\215","\216",
"\217","\220","\221","\222","\223","\224","\225","\226",
"\227","\230","\231","\232","\233","\234","\235","\236",
"\237","\240","END"};
Float_t xrange = 18;
Float_t yrange = 25;
Int_t w = 650;
Int_t h = w*yrange/xrange;
TCanvas *c1 = new TCanvas("c1","c1",200,10,w,h);
c1->Range(0,0,xrange,yrange);
TText *t = new TText(0,0,"a");
t->SetTextSize(0.02);
t->SetTextFont(62);
t->SetTextAlign(22);
table(0.5,0.5*xrange-0.5,yrange,t,symbol1,0);
table(0.5*xrange+0.5,xrange-0.5,yrange,t,symbol2,0);
TText *tlabel = new TText(0,0,"a");
tlabel->SetTextFont(72);
tlabel->SetTextSize(0.018);
tlabel->SetTextAlign(22);
tlabel->DrawText(0.5*xrange,1.3,
"Input characters are standard keyboard characters");
c1->Modified();
c1->Update();
c1->Print("pstable1.ps");
TCanvas *c2 = new TCanvas("c2","c2",220,20,w,h);
c2->Range(0,0,xrange,yrange);
table(0.5,0.5*xrange-0.5,yrange,t,symbol3,1);
table(0.5*xrange+0.5,xrange-0.5,yrange,t,symbol4,1);
tlabel->DrawText(0.5*xrange,1.3,
"Input characters using backslash and octal numbers");
c2->Modified();
c2->Update();
c2->Print("pstable2.ps");
TCanvas *c3 = new TCanvas("c3","c3",240,20,w,h);
c3->Range(0,0,xrange,yrange);
table(0.5,0.5*xrange-0.5,yrange,t,symbol5,1);
tlabel->DrawText(0.5*xrange,1.3,
"Input characters using backslash and octal numbers");
c3->Modified();
c3->Update();
c3->Print("pstable3.ps");
}
void analyze() {
TCanvas *c1 = new TCanvas("c1","Analyze.mac",620,790);
c1->Range(-1,0,19,30);
TPaveLabel *pl1 = new TPaveLabel(0,27,3.5,29,"Analyze");
pl1->SetFillColor(42);
pl1->Draw();
TPaveText *pt1 = new TPaveText(0,22.8,4,25.2);
TText *t1 = pt1->AddText("Parenthesis matching");
TText *t2 = pt1->AddText("Remove unnecessary");
TText *t2a = pt1->AddText("parenthesis");
pt1->Draw();
TPaveText *pt2 = new TPaveText(6,23,10,25);
TText *t3 = pt2->AddText("break of");
TText *t4 = pt2->AddText("Analyze");
pt2->Draw();
TPaveText *pt3 = new TPaveText(0,19,4,21);
t4=pt3->AddText("look for simple");
TText *t5 = pt3->AddText("operators");
pt3->Draw();
TPaveText *pt4 = new TPaveText(0,15,4,17);
TText *t6 = pt4->AddText("look for an already");
TText *t7 = pt4->AddText("defined expression");
pt4->Draw();
TPaveText *pt5 = new TPaveText(0,11,4,13);
TText *t8 = pt5->AddText("look for usual");
TText *t9 = pt5->AddText("functions :cos sin ..");
pt5->Draw();
TPaveText *pt6 = new TPaveText(0,7,4,9);
TText *t10 = pt6->AddText("look for a");
TText *t11 = pt6->AddText("numeric value");
pt6->Draw();
TPaveText *pt7 = new TPaveText(6,18.5,10,21.5);
TText *t12 = pt7->AddText("Analyze left and");
TText *t13 = pt7->AddText("right part of");
TText *t14 = pt7->AddText("the expression");
pt7->Draw();
TPaveText *pt8 = new TPaveText(6,15,10,17);
TText *t15 = pt8->AddText("Replace expression");
pt8->Draw();
TPaveText *pt9 = new TPaveText(6,11,10,13);
TText *t16 = pt9->AddText("Analyze");
pt9->SetFillColor(42);
pt9->Draw();
TPaveText *pt10 = new TPaveText(6,7,10,9);
TText *t17 = pt10->AddText("Error");
TText *t18 = pt10->AddText("Break of Analyze");
pt10->Draw();
TPaveText *pt11 = new TPaveText(14,22,17,24);
pt11->SetFillColor(42);
TText *t19 = pt11->AddText("Analyze");
TText *t19a = pt11->AddText("Left");
pt11->Draw();
TPaveText *pt12 = new TPaveText(14,19,17,21);
pt12->SetFillColor(42);
TText *t20 = pt12->AddText("Analyze");
TText *t20a = pt12->AddText("Right");
pt12->Draw();
TPaveText *pt13 = new TPaveText(14,15,18,18);
TText *t21 = pt13->AddText("StackNumber++");
TText *t22 = pt13->AddText("operator[StackNumber]");
TText *t23 = pt13->AddText("= operator found");
pt13->Draw();
TPaveText *pt14 = new TPaveText(12,10.8,17,13.2);
TText *t24 = pt14->AddText("StackNumber++");
TText *t25 = pt14->AddText("operator[StackNumber]");
TText *t26 = pt14->AddText("= function found");
pt14->Draw();
TPaveText *pt15 = new TPaveText(6,7,10,9);
TText *t27 = pt15->AddText("Error");
TText *t28 = pt15->AddText("break of Analyze");
pt15->Draw();
TPaveText *pt16 = new TPaveText(0,2,7,5);
TText *t29 = pt16->AddText("StackNumber++");
TText *t30 = pt16->AddText("operator[StackNumber] = 0");
TText *t31 = pt16->AddText("value[StackNumber] = value found");
pt16->Draw();
TArrow *ar = new TArrow(2,27,2,25.4,0.012,"|>");
ar->SetFillColor(1);
ar->Draw();
ar->DrawArrow(2,22.8,2,21.2,0.012,"|>");
ar->DrawArrow(2,19,2,17.2,0.012,"|>");
ar->DrawArrow(2,15,2,13.2,0.012,"|>");
ar->DrawArrow(2,11,2, 9.2,0.012,"|>");
ar->DrawArrow(2, 7,2, 5.2,0.012,"|>");
ar->DrawArrow(4,24,6,24,0.012,"|>");
ar->DrawArrow(4,20,6,20,0.012,"|>");
ar->DrawArrow(4,16,6,16,0.012,"|>");
ar->DrawArrow(4,12,6,12,0.012,"|>");
ar->DrawArrow(4, 8,6, 8,0.012,"|>");
ar->DrawArrow(10,20,14,20,0.012,"|>");
ar->DrawArrow(12,23,14,23,0.012,"|>");
ar->DrawArrow(12,16.5,14,16.5,0.012,"|>");
ar->DrawArrow(10,12,12,12,0.012,"|>");
TText *ta = new TText(2.2,22.2,"err = 0");
ta->SetTextFont(71);
ta->SetTextSize(0.015);
ta->SetTextColor(4);
ta->SetTextAlign(12);
ta->Draw();
ta->DrawText(2.2,18.2,"not found");
ta->DrawText(2.2,6.2,"found");
TText *tb = new TText(4.2,24.1,"err != 0");
tb->SetTextFont(71);
tb->SetTextSize(0.015);
tb->SetTextColor(4);
tb->SetTextAlign(11);
tb->Draw();
tb->DrawText(4.2,20.1,"found");
tb->DrawText(4.2,16.1,"found");
tb->DrawText(4.2,12.1,"found");
tb->DrawText(4.2, 8.1,"not found");
TLine *l1 = new TLine(12,16.5,12,23);
l1->Draw();
}
void process_syst_plot(TH2F *rhisto,string saveto,int scantype, std::string scanx = "") {
TH2F *histo = prep_histo(rhisto,scantype); // this is to be independent of the style used at creation time
float rightmargin=gStyle->GetPadRightMargin();
gStyle->SetPadRightMargin(0.20);
TString name = rhisto->GetName();
if(name.Contains("Nevents")) gStyle->SetPadRightMargin(0.22);
TCanvas *can = new TCanvas("syst_plot","Systematics Plot");
set_range(histo,scantype,true);
histo->GetZaxis()->CenterTitle();
gStyle->SetStripDecimals(false);
histo->GetXaxis()->SetDecimals(true);
if(name.Contains("efficiency")) {
histo->GetZaxis()->SetTitle("A #times #varepsilon (#geq 1 Z(ll))");
histo->GetZaxis()->CenterTitle(0);
//histo->GetZaxis()->SetRangeUser(0.0,0.15);
}
if(name.Contains("Nevents")) {
histo->GetZaxis()->SetTitle("N(events)");
histo->GetZaxis()->SetTitleOffset(histo->GetZaxis()->GetTitleOffset()+0.4);
}
if(name.Contains("sysjes")) {
histo->GetZaxis()->SetTitle("Jet Energy Scale");
histo->GetZaxis()->SetRangeUser(0.0,0.2);
}
if(name.Contains("sysjsu")) {
histo->GetZaxis()->SetTitle("JZB Scale Uncertainty");
histo->GetZaxis()->SetRangeUser(0.0,0.5);
}
if(name.Contains("sysresmap")) {
histo->GetZaxis()->SetTitle("Resulution");
histo->GetZaxis()->SetRangeUser(0.0,0.5);
}
if(name.Contains("sysstatmap")) {
histo->GetZaxis()->SetTitle("Statistical Error");
histo->GetZaxis()->SetRangeUser(0.0,0.01);
}
if(name.Contains("systotmap")) {
histo->GetZaxis()->SetTitle("All Systematic Errors");
histo->GetZaxis()->SetRangeUser(0.0,0.5);
}
histo->Draw("COLZ");
DrawPrelim();
float xpos_of_text = 0.22;
if(name.Contains("_noscefficiencymap")) {
// Add some more decorations on these publication plots
write_SMS_text( scantype, scanx, xpos_of_text );
// Decode name and get JZB cut
TPRegexp pat("\\d+$");
size_t index = name.Index(pat,0);
string cut = string("JZB > ")+(name(index,name.Length()-index).Data())+" GeV";
//string cut = string("#splitline{JZB > ")+(name(index,name.Length()-index).Data())+" GeV}{n_{jets} #geq 3}"; //paper style
TText *text = write_text(xpos_of_text,0.73,cut);
text->SetTextAlign(11);
text->SetTextSize(0.035);
text->Draw();
draw_diagonal_xchange( scantype, scanx );
}
CompleteSave(can,(saveto+(string)histo->GetName()));
gStyle->SetPadRightMargin(rightmargin);
delete can;
}
void EMCDistribution(TString gain = "CALIB", bool log_scale = false)
{
TText *t;
TCanvas *c1 = new TCanvas(
"EMCDistribution_" + gain + TString(log_scale ? "_Log" : "") + cuts,
"EMCDistribution_" + gain + TString(log_scale ? "_Log" : "") + cuts, 1800,
1000);
c1->Divide(8, 8, 0., 0.01);
int idx = 1;
TPad *p;
for (int iphi = 8 - 1; iphi >= 0; iphi--)
{
for (int ieta = 0; ieta < 8; ieta++)
{
p = (TPad *) c1->cd(idx++);
c1->Update();
p->SetLogy();
p->SetGridx(0);
p->SetGridy(0);
TString hname = Form("hEnergy_ieta%d_iphi%d", ieta, iphi) + TString(log_scale ? "_Log" : "");
TH1 *h = NULL;
if (log_scale)
h = new TH1F(hname,
Form(";Calibrated Tower Energy Sum (GeV);Count / bin"), 300,
5e-3, 3096);
else
// h = new TH1F(hname,
// Form(";Calibrated Tower Energy Sum (GeV);Count / bin"), 196,
// 1900, 2096);
h = new TH1F(hname,
Form(";Calibrated Tower Energy Sum (GeV);Count / bin"), 596,
-96, 500);
h->SetLineWidth(0);
h->SetLineColor(kBlue + 3);
h->SetFillColor(kBlue + 3);
h->GetXaxis()->SetTitleSize(.09);
h->GetXaxis()->SetLabelSize(.08);
h->GetYaxis()->SetLabelSize(.08);
if (log_scale)
QAHistManagerDef::useLogBins(h->GetXaxis());
T->Draw(
"TOWER_" + gain + "_CEMC[].get_energy_power_law_exp()>>" + hname,
Form(
"TOWER_%s_CEMC[].get_bineta()==%d && TOWER_%s_CEMC[].get_binphi()==%d",
gain.Data(), ieta, gain.Data(), iphi),
"");
TText *t = new TText(.9, .9, Form("Col%d Row%d", ieta, iphi));
t->SetTextAlign(33);
t->SetTextSize(.15);
t->SetNDC();
t->Draw();
// return;
}
}
SaveCanvas(c1,
TString(_file0->GetName()) + TString("_DrawPrototype3EMCalTower_") + TString(c1->GetName()), false);
}
PigsGUI::PigsGUI(const TGWindow *p) : TGMainFrame(p, fGUIsizeX, fGUIsizeY) {
// Creates the GUI
if(fVerbose) std::cout<<__PRETTY_FUNCTION__ << std::endl;
daq = 0; storage = 0; ev = 0; // Initialize local variables
year = month = day = hour = min = sec = 0;
fAcqThread = 0;
keepAcquiring = kFALSE;
useIntegration = kTRUE;
const int32_t fHistColors[] = { kMagenta+1, kGreen+1, kBlue+1, kRed+1 };
fAboutMsg = (char*)
"\n"
"\n"
" _____ _____ ______ _______\n"
" |_____] | | ____ |______\n"
" | __|__ |_____| ______|\n"
"\n"
"\n"
" *** Position Indicating Gamma Sensor ***\n"
" * CAEN DT-5781 Data Acquisition System *\n"
" Four Channel Version\n"
"\n"
" by Ondrej Chvala <*****@*****.**>\n"
" version 0.097, July 2015\n"
" https://github.com/ondrejch/DAQ-DT5781\n"
" GNU/GPL";
int32_t i = 0; // helper variable
for (i=0; i<4; i++) {
fScaleFactor[i] = 1.0;
fNormAvgH[i] = 0;
}
fIntegralMin = 1;
fIntegralMax = 16384;
// *** Main GUI window ***
fMainGUIFrame = new TGMainFrame(gClient->GetRoot(),10,10,kMainFrame | kVerticalFrame);
#include "fpigsicon.xpm"
TImage *tmpicon = TImage::Create();
tmpicon->SetImageBuffer((char**)fpigsicon_xpm, TImage::kXpm);
gVirtualX->SetIconPixmap(fMainGUIFrame->GetId(),tmpicon->GetPixmap());
delete tmpicon;
fMainGUIFrame->SetName("fMainGUIFrame");
fMainGUIFrame->SetWindowName("F-PIGS"); // GUI window name
fMainGUIFrame->SetLayoutBroken(kTRUE);
ufont = gClient->GetFont("-*-*-bold-r-*-*-16-*-*-*-*-*-*-*");
// ufont = gClient->GetFont("-*-helvetica-medium-r-normal-*-14-*-*-*-*-*-iso8859-1");
// ufont = gClient->GetFont("-urw-nimbus sans l-bold-r-normal--0-0-0-0-p-0-iso8859-1");
valTitle.fMask = kGCForeground | kGCBackground | kGCFillStyle | kGCFont | kGCGraphicsExposures;
gClient->GetColorByName("#0000FF",valTitle.fForeground);
gClient->GetColorByName("#e0e0e0",valTitle.fBackground);
valTitle.fFillStyle = kFillSolid;
valTitle.fFont = ufont->GetFontHandle();
valTitle.fGraphicsExposures = kFALSE;
uGC = gClient->GetGC(&valTitle, kTRUE);
fMainTitle = new TGLabel(fMainGUIFrame,"Four-channel Position Identifying Gamma Sensor (F-PIGS)",
uGC->GetGC(),ufont->GetFontStruct());
fMainTitle->SetTextJustify(36);
fMainTitle->SetMargins(0,0,0,0);
fMainTitle->SetWrapLength(-1);
fMainGUIFrame->AddFrame(fMainTitle, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
fMainTitle->MoveResize(0,0,fGUIsizeX-4,32);
fMainGUIFrame->Connect("CloseWindow()", "PigsGUI", this, "~PigsGUI()"); // call class destructor on alt+f4
fMainGUIFrame->DontCallClose();
// Buttons for main GUI
fStartDAQ = new TGTextButton(fMainGUIFrame, "Start DAQ"); // start DAQ
fStartDAQ->SetTextJustify(36);
fStartDAQ->SetMargins(0,0,0,0);
fStartDAQ->Resize(90,25);
fMainGUIFrame->AddFrame(fStartDAQ, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
fStartDAQ->MoveResize(50,fGUIsizeY-30,90,25);
gClient->GetColorByName("green", fColor);
fStartDAQ->ChangeBackground(fColor);
fStartDAQ->SetState(kButtonDisabled);
fStartDAQ->Connect("Clicked()","PigsGUI",this,"RunAcquisition()");
fStopDAQ = new TGTextButton(fMainGUIFrame, "Stop DAQ"); // stop DAQ
fStopDAQ->SetTextJustify(36);
fStopDAQ->SetMargins(0,0,0,0);
fStopDAQ->Resize(90,25);
fMainGUIFrame->AddFrame(fStopDAQ, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
fStopDAQ->MoveResize(fGUIsizeX-50-90,fGUIsizeY-30,90,25);
gClient->GetColorByName("red", fColor);
fStopDAQ->ChangeBackground(fColor);
fStopDAQ->SetState(kButtonDisabled);
fStopDAQ->Connect("Clicked()","PigsGUI",this,"StopAcquisition()");
fExitDAQ = new TGTextButton(fMainGUIFrame, "Exit DAQ"); // exit DAQ
fExitDAQ->SetTextJustify(36);
fExitDAQ->SetMargins(0,0,0,0);
fExitDAQ->Resize(90,25);
fMainGUIFrame->AddFrame(fExitDAQ, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
fExitDAQ->MoveResize(fGUIsizeX/2-45,fGUIsizeY-30,90,25);
fExitDAQ->Connect("Clicked()","PigsGUI",this,"~PigsGUI()");
// *** Tab widget ****
fTabHolder = new TGTab(fMainGUIFrame,fGUIsizeX-4,fGUIsizeX-4);//,uGC->GetGC());
// *** Container of "CurrentHistogram" ***
fCurHistFrame = fTabHolder->AddTab("CurrentHistogram");
fCurHistFrame->SetLayoutManager(new TGVerticalLayout(fCurHistFrame));
// embedded canvas
fLatestHistoCanvas = new TRootEmbeddedCanvas("CurrentHEC",fCurHistFrame,fGUIsizeX-10,fGUIsizeY-140);
Int_t wfLatestHistoCanvas = fLatestHistoCanvas->GetCanvasWindowId();
cCurrHCanvas = new TCanvas("cCurrHCanvas", 10, 10, wfLatestHistoCanvas);
fLatestHistoCanvas->AdoptCanvas(cCurrHCanvas);
cCurrHCanvas->Divide(2,2);
fCurHistFrame->AddFrame(fLatestHistoCanvas, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
fHCurrHProgressBar = new TGHProgressBar(fCurHistFrame,fGUIsizeX-5);
fHCurrHProgressBar->SetFillType(TGProgressBar::kBlockFill);
fHCurrHProgressBar->ChangeOptions(kSunkenFrame | kDoubleBorder | kOwnBackground);
// will reflect user color changes
gClient->GetColorByName("#ffffff",fColor);
fHCurrHProgressBar->SetBackgroundColor(fColor);
fHCurrHProgressBar->SetPosition(1);
fCurHistFrame->AddFrame(fHCurrHProgressBar, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
// *** Container of "History" ***
fTabHisto = fTabHolder->AddTab("History");
fTabHisto->SetLayoutManager(new TGVerticalLayout(fTabHisto));
// embedded canvas
fLastMeas = new TRootEmbeddedCanvas("HistoryHEC",fTabHisto,fGUIsizeX-10,fGUIsizeY-110);
Int_t wfLastNspectra = fLastMeas->GetCanvasWindowId();
cLastMeas = new TCanvas("cLastMeas", 10, 10, wfLastNspectra);
fLastMeas->AdoptCanvas(cLastMeas);
fMG = new TMultiGraph("fMG","");
for (i=0; i<4; i++) {
fGraph[i] = new TGraph();
fGraph[i]->SetName(Form("gCh%d",i));
fGraph[i]->SetDrawOption("AP");
fGraph[i]->SetMarkerColor(fHistColors[i]);
fGraph[i]->SetMarkerStyle(21);
fGraph[i]->SetMarkerSize(2.0);
fGraph[i]->SetLineWidth(0.5);
fGraph[i]->SetLineColor(i+12);
fGraph[i]->SetFillStyle(0);
fMG->Add(fGraph[i]);
}
fTabHisto->AddFrame(fLastMeas, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
// *** Container of "Average" ***
fTabSum = fTabHolder->AddTab("Average");
fTabSum->SetLayoutManager(new TGVerticalLayout(fTabSum));
// embedded canvas
fSumSpectra = new TRootEmbeddedCanvas("SumHEC",fTabSum,fGUIsizeX-10,fGUIsizeY-110);
Int_t wfSumSpectra = fSumSpectra->GetCanvasWindowId();
cSumSpectra = new TCanvas("cSumSpectra", 10, 10, wfSumSpectra);
fSumSpectra->AdoptCanvas(cSumSpectra);
cSumSpectra->Divide(2,2);
for (i=1; i<5; i++) {
cSumSpectra->GetPad(i)->SetLogx();
cSumSpectra->GetPad(i)->SetLogy();
}
fTabSum->AddFrame(fSumSpectra, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
// *** Container of "Arrow" ***
fTabArrow = fTabHolder->AddTab("Arrow");
fTabArrow->SetLayoutManager(new TGVerticalLayout(fTabArrow));
fArrowECanvas = new TRootEmbeddedCanvas("ArrowHEC",fTabArrow,fGUIsizeX-10,fGUIsizeY-110);
Int_t wfArrowECanvas = fArrowECanvas->GetCanvasWindowId();
cArrowCanvas = new TCanvas("cArrowCanvas", 5, 5, wfArrowECanvas);
// Draw compass circle
TEllipse *el1 = new TEllipse(0.5,0.5,0.48,0.48);
el1->SetFillColor(14);
el1->SetFillStyle(1001);
el1->SetLineColor(1);
el1->SetLineWidth(6);
el1->Draw();
// Add bearing labels
TText *north = new TText(0.5,0.9,"N");
north->SetTextColor(2);
north->SetTextSize(0.1);
north->SetTextAlign(12);
north->SetTextAlign(21);
north->Draw();
TText *south = new TText(0.5,0.04,"S");
south->SetTextSize(0.1);
south->SetTextAlign(12);
south->SetTextAlign(21);
south->Draw();
TText *east = new TText(0.92,0.5,"E");
east->SetTextSize(0.1);
east->SetTextAlign(12);
east->SetTextAlign(21);
east->Draw();
TText *west = new TText(0.08,0.5,"W");
west->SetTextSize(0.1);
west->SetTextAlign(12);
west->SetTextAlign(21);
west->Draw();
// Draw initial arrow pointing North
ar1 = new TArrow(0.5,0.3,0.5,0.7,0.3,"|>");
ar1->SetAngle(30);
ar1->SetLineWidth(5);
ar1->SetFillColor(4);
ar1->Draw();
fArrowECanvas->AdoptCanvas(cArrowCanvas);
fTabArrow->AddFrame(fArrowECanvas, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
// *** Container of "Config" ***
fTabConfig = fTabHolder->AddTab("Config");
fTabConfig->SetLayoutManager(new TGVerticalLayout(fTabConfig));
// Acquisition time settings
fControlFrame = new TGGroupFrame(fTabConfig, "Acquisition time [sec]");
fControlFrame->SetTitlePos(TGGroupFrame::kCenter);
// Acquisition time entry
fAcqTimeFrame = new TGCompositeFrame(fControlFrame, 200, 10, kHorizontalFrame);
fAcqTimeSlider = new TGHSlider(fAcqTimeFrame,300,kSlider1 | kScaleBoth,-1);
//fAcqTimeSlider = new TGHSlider(fControlFrame,300,kSlider1 | kScaleBoth,-1);
fAcqTimeSlider->Connect("PositionChanged(Int_t)", "PigsGUI", this, "SetAcquisitionLoopTimeSlider()");
//fAcqTimeSlider->Connect("PositionChanged(Int_t)", "TGLabel", fAcqTimeLabel, "SetText(Int_t)");
fAcqTimeSlider->SetRange(1,6000); // time in 100 ms increments => [0.1 - 600 sec]
fAcqTimeSlider->SetPosition(fDefaultAcqTime*10); // 10 second acquire time by default
fAcqTimeFrame->AddFrame(fAcqTimeSlider, new TGLayoutHints(kLHintsNormal, 5, 5, 5, 5));
//fAcqTimeEntry = new TGNumberEntry(fControlFrame, (Double_t) fDefaultAcqTime ,5,-1, TGNumberFormat::kNESRealOne, TGNumberFormat::kNEAPositive,TGNumberFormat::kNELLimitMinMax, 0.1, 600);
fAcqTimeEntry = new TGNumberEntry(fAcqTimeFrame, (Double_t) fDefaultAcqTime ,5,-1, TGNumberFormat::kNESRealOne, TGNumberFormat::kNEAPositive,TGNumberFormat::kNELLimitMinMax, 0.1, 600);
fAcqTimeEntry->GetNumberEntry()->SetToolTipText("Time for one DAQ loop in seconds.");
fAcqTimeEntry->GetNumberEntry()->Connect("TextChanged(char*)", "PigsGUI", this,
"SetAcquisitionLoopTimeNumberEntry()");
fAcqTimeEntry->GetNumberEntry()->Connect("ReturnPressed()", "PigsGUI", this,
"SetAcquisitionLoopTimeNumberEntry()");
fAcqTimeFrame->AddFrame(fAcqTimeEntry, new TGLayoutHints(kLHintsNormal, 5, 5, 5, 5));
//fControlFrame->AddFrame(fAcqTimeEntry, new TGLayoutHints(kLHintsNormal, 5, 5, 5, 5));
// Disable acquisition time adjustment until DAQ initialized
fAcqTimeEntry->SetState(kFALSE);
fAcqTimeSlider->SetState(kFALSE);
/*
fAcqTimeLabelFrame = new TGCompositeFrame(fTabConfig, 100, 10, kHorizontalFrame);
fAcqTimeLabelText = new TGLabel(fAcqTimeLabelFrame,"Acquire time:",
uGC->GetGC(),ufont->GetFontStruct());
fAcqTimeLabelText->SetTextJustify(kTextLeft);
fAcqTimeLabelText->SetWrapLength(-1);
fAcqTimeLabel = new TGLabel(fAcqTimeLabelFrame," ",
uGC->GetGC(),ufont->GetFontStruct());
fAcqTimeLabel->SetTextColor(0x0066ff);
fAcqTimeLabel->SetTextJustify(kTextCenterX);
fAcqTimeLabel->SetWrapLength(-1);
fAcqTimeLabel->SetMinWidth(3);
fAcqTimeLabelFrame->AddFrame(fAcqTimeLabelText, new TGLayoutHints(kLHintsLeft, 10, 5, 10, 10));
fAcqTimeLabelFrame->AddFrame(fAcqTimeLabel, new TGLayoutHints(kLHintsRight, 0, 10, 10, 10 ));
*/
fControlFrame->AddFrame(fAcqTimeFrame, new TGLayoutHints(kLHintsTop, 10, 10, 5, 5));
// fControlFrame->AddFrame(fAcqTimeLabelFrame, new TGLayoutHints(kLHintsBottom, 10, 10, 5, 5));
fTabConfig->AddFrame(fControlFrame, new TGLayoutHints(kLHintsNormal, 10, 10, 10, 10));
// Scale Factor setting
fScalerFrame = new TGGroupFrame(fTabConfig, "Channel gain compensation");
fScalerFrame->SetTitlePos(TGGroupFrame::kCenter);
for (i=0; i<4; i++){
fScalerInput[i] = new PigsScalerInput(fScalerFrame, Form("ch %d scaling", i));
fScalerInput[i]->GetEntry()->Connect("TextChanged(char*)", "PigsGUI", this,
Form("SetGainScalerCh%d()",i));
fScalerInput[i]->GetEntry()->SetToolTipText(
"Channel gain is a multiplicative factor used in detector response calculation.");
fScalerFrame->AddFrame(fScalerInput[i], new TGLayoutHints(kLHintsNormal, 0, 0, 2, 2));
}
fTabConfig->AddFrame(fScalerFrame, new TGLayoutHints(kLHintsNormal, 10, 10, 10, 10));
// Integration Limits
fIntLimFrame = new TGGroupFrame(fTabConfig, "ADC window for integration");
fIntLimFrame->SetTitlePos(TGGroupFrame::kCenter);
fUseIntegration = new TGCheckButton(fIntLimFrame, "Energy integration On/Off");
fUseIntegration->SetOn(kFALSE); // Start with regular counts
fUseIntegration->SetToolTipText("If enabled, the detector response is calculated by integrating "
"the energy deposited in ADC bins within the limits specidied below.\n"
"If disabled, the hit count is used as a detector response.");
fUseIntegration->Connect("Toggled(Bool_t)", "PigsGUI", this, "ToggleUseIntegration()");
fIntLimFrame->AddFrame(fUseIntegration, new TGLayoutHints(kLHintsNormal, 0, 0, 2, 2));
fIntLimInputMin = new PigsIntLimInput(fIntLimFrame, "Lower limit");
fIntLimInputMin->GetEntry()->SetIntNumber(fIntegralMin);
fIntLimInputMin->GetEntry()->Connect("TextChanged(char*)", "PigsGUI", this, "SetIntegralLimitMin()");
fIntLimFrame->AddFrame(fIntLimInputMin, new TGLayoutHints(kLHintsNormal, 0, 0, 2, 2));
fIntLimInputMax = new PigsIntLimInput(fIntLimFrame, "Upper limit");
fIntLimInputMax->GetEntry()->SetIntNumber(fIntegralMax);
fIntLimInputMax->GetEntry()->Connect("TextChanged(char*)", "PigsGUI", this, "SetIntegralLimitMax()");
fIntLimFrame->AddFrame(fIntLimInputMax, new TGLayoutHints(kLHintsNormal, 0, 0, 2, 2));
fTabConfig->AddFrame(fIntLimFrame, new TGLayoutHints(kLHintsNormal, 10, 10, 10, 10));
this->ToggleUseIntegration(); // Updates the integration entries' status
// *** Container of "DT5781" ***
fTabDT5781 = fTabHolder->AddTab("DT5781");
fTabDT5781->SetLayoutManager(new TGVerticalLayout(fTabDT5781));
gClient->GetColorByName("white", fColor);
fDTinfo = new TGTextView(fTabDT5781,fGUIsizeX-120,fGUIsizeY-150,"DAQ not initialized.",kSunkenFrame,fColor);
fTabDT5781->AddFrame(fDTinfo, new TGLayoutHints(kLHintsNormal));
fDTinfo->MoveResize(10,50,fGUIsizeX-120,fGUIsizeY-150);
fInitDAQ = new TGTextButton(fTabDT5781, "Init DAQ"); // button InitDAQ
fInitDAQ->SetTextJustify(36);
fInitDAQ->SetMargins(0,0,0,0);
fInitDAQ->Resize(90,25);
fMainGUIFrame->AddFrame(fInitDAQ, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
gClient->GetColorByName("light blue", fColor);
fInitDAQ->ChangeBackground(fColor);
fInitDAQ->MoveResize(fGUIsizeX-50-60,50,90,25);
fInitDAQ->Connect("Clicked()","PigsGUI",this,"InitDAQ()");
fDisconnectDAQ = new TGTextButton(fTabDT5781, "Disconnect DAQ"); // buttons DisconnectDAQ
fDisconnectDAQ->SetTextJustify(36);
fDisconnectDAQ->SetMargins(0,0,0,0);
fDisconnectDAQ->Resize(90,25);
fMainGUIFrame->AddFrame(fDisconnectDAQ, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
gClient->GetColorByName("light blue", fColor);
fDisconnectDAQ->MoveResize(fGUIsizeX-50-60,150,90,25);
fDisconnectDAQ->ChangeBackground(fColor);
fDisconnectDAQ->Connect("Clicked()","PigsGUI",this,"DisconnectDAQ()");
// *** Container of "About" ***
fTabAbout = fTabHolder->AddTab("About");
fTabAbout->SetLayoutManager(new TGVerticalLayout(fTabAbout));
ufont = gClient->GetFont("-*-fixed-medium-r-*-*-15-*-*-*-*-*-*-*");
fAboutText = new TGTextView(fTabAbout,1,1,"SPIGS",kSunkenFrame);
fAboutText->SetFont(ufont->GetFontStruct());
fTabAbout->AddFrame(fAboutText, new TGLayoutHints(kLHintsNormal));
fAboutText->LoadBuffer(fAboutMsg);
//-------------------------------------------------------------------------
// Change to the starting tab
fTabHolder->SetTab("DT5781");
// Display the GUI
fTabHolder->Resize(fTabHolder->GetDefaultSize());
fMainGUIFrame->AddFrame(fTabHolder, new TGLayoutHints(kLHintsLeft | kLHintsTop,2,2,2,2));
fTabHolder->MoveResize(0,32,fGUIsizeX-2,fGUIsizeY-80);
fMainGUIFrame->SetMWMHints(kMWMDecorAll, kMWMFuncAll, kMWMInputModeless);
fMainGUIFrame->MapSubwindows();
fMainGUIFrame->Resize(fMainGUIFrame->GetDefaultSize());
fMainGUIFrame->MapWindow();
fMainGUIFrame->Resize(fGUIsizeX,fGUIsizeY);
static const int32_t tmpw = 410; // Constants for About window placement
static const int32_t tmph = 260;
fAboutText->MoveResize((fGUIsizeX-tmpw)/2,(fGUIsizeY-tmph)/3,tmpw,tmph);
}
void event(){
TCanvas *c1 = new TCanvas("c1","ROOT Event description",700,500);
c1->Range(0,0,14,15.5);
TPaveText *event = new TPaveText(1,13,3,15);
event->SetFillColor(11);
event->Draw();
event->AddText("Event");
TLine *line = new TLine(1.1,13,1.1,1.5);
line->SetLineWidth(2);
line->Draw();
line->DrawLine(1.3,13,1.3,3.5);
line->DrawLine(1.5,13,1.5,5.5);
line->DrawLine(1.7,13,1.7,7.5);
line->DrawLine(1.9,13,1.9,9.5);
line->DrawLine(2.1,13,2.1,11.5);
TArrow *arrow = new TArrow(1.1,1.5,3.9,1.5,0.02,"|>");
arrow->SetFillStyle(1001);
arrow->SetFillColor(1);
arrow->Draw();
arrow->DrawArrow(1.3,3.5,3.9,3.5,0.02,"|>");
arrow->DrawArrow(1.5,5.5,3.9,5.5,0.02,"|>");
arrow->DrawArrow(1.7,7.5,3.9,7.5,0.02,"|>");
arrow->DrawArrow(1.9,9.5,3.9,9.5,0.02,"|>");
arrow->DrawArrow(2.1,11.5,3.9,11.5,0.02,"|>");
TPaveText *p1 = new TPaveText(4,1,11,2);
p1->SetTextAlign(12);
p1->SetFillColor(42);
p1->AddText("1 Mbyte");
p1->Draw();
TPaveText *p2 = new TPaveText(4,3,10,4);
p2->SetTextAlign(12);
p2->SetFillColor(42);
p2->AddText("100 Kbytes");
p2->Draw();
TPaveText *p3 = new TPaveText(4,5,9,6);
p3->SetTextAlign(12);
p3->SetFillColor(42);
p3->AddText("10 Kbytes");
p3->Draw();
TPaveText *p4 = new TPaveText(4,7,8,8);
p4->SetTextAlign(12);
p4->SetFillColor(42);
p4->AddText("1 Kbytes");
p4->Draw();
TPaveText *p5 = new TPaveText(4,9,7,10);
p5->SetTextAlign(12);
p5->SetFillColor(42);
p5->AddText("100 bytes");
p5->Draw();
TPaveText *p6 = new TPaveText(4,11,6,12);
p6->SetTextAlign(12);
p6->SetFillColor(42);
p6->AddText("10 bytes");
p6->Draw();
TText *text = new TText();
text->SetTextAlign(12);
text->SetTextSize(0.04);
text->SetTextFont(72);
text->DrawText(6.2,11.5,"Header:Event_flag");
text->DrawText(7.2,9.5,"Trigger_Info");
text->DrawText(8.2,7.5,"Muon_Detector: TOF");
text->DrawText(9.2,5.5,"Calorimeters");
text->DrawText(10.2,3.5,"Forward_Detectors");
text->DrawText(11.2,1.5,"TPCs");
}
void plotpurity()
{
Bool_t dolines = kFALSE;
//Bool_t dolines = kTRUE;
TString lines = "";
TString xn = "";
//TString xn = "_onlyQCD";
//TString xn = "_onlyGJ";
//TString xn = "_denominator";
if(dolines){lines="/lines";}
std::vector<TString> rebins;
rebins.push_back("");
// rebins.push_back("rebin1");
// rebins.push_back("rebin2");
// rebins.push_back("rebin3");
// rebins.push_back("rebin4");
// rebins.push_back("rebin5");
std::vector<TString> ptranges;
// ptranges.push_back("90to120");
// ptranges.push_back("120to155");
// ptranges.push_back("155to175");
ptranges.push_back("175to190");
ptranges.push_back("190to250");
ptranges.push_back("250to400");
ptranges.push_back("400to700");
ptranges.push_back("700to1000");
ptranges.push_back("175to1000");
//ptranges.push_back("250to1000");
//ptranges.push_back("400to1000");
std::vector<TString> isovars;
//isovars.push_back("wchiso");
isovars.push_back("chiso");
// isovars.push_back("sieieF5x5");
// isovars.push_back("sieipF5x5");
// isovars.push_back("sipipF5x5");
// isovars.push_back("rho");
// isovars.push_back("nVtx");
// isovars.push_back("phoet");
// isovars.push_back("pfMET");
std::vector<TString> cuts;
cuts.push_back("a_idnc");
cuts.push_back("b_idnc_mL30");
cuts.push_back("c_idnc_trig");
cuts.push_back("d_idnc_mL30_trig");
cuts.push_back("e_idnc_t175");
cuts.push_back("f_idnc_mL30_t175");
cuts.push_back("g_idnc_t250");
cuts.push_back("h_idnc_mL30_t250");
cuts.push_back("i_idnc_mL30_allt");
// cuts.push_back("oldj");
TString inpath = "/afs/hep.wisc.edu/cms/tperry/LoneG_slc6_491_CMSSW_7_4_14/src/LoneGamma/qcdStudy/gitignore/Lire/analyzed";
TString outpath = "/afs/hep.wisc.edu/cms/tperry/LoneG_slc6_491_CMSSW_7_4_14/src/LoneGamma/qcdStudy/gitignore/Lire/plots";
TString wwwpath = "/afs/hep.wisc.edu/home/tperry/www/MonoPhoton/qcdPlots/Lire/purity";
TString inname_GJ = inpath+"/purity_GJets_Merged.root";
TString inname_QCD = inpath+"/purity_QCD_Merged.root";
// TString inname_GJ = inpath+"/purity_mrg4bins_GJets.root";
// TString inname_QCD = inpath+"/purity_mrg4bins_QCD.root";
//TString inname_GJ = inpath+"/purity_mrg3bins_GJets.root";
//TString inname_QCD = inpath+"/purity_mrg3bins_QCD.root";
TFile *infile_GJ = new TFile(inname_GJ);
TFile *infile_QCD = new TFile(inname_QCD);
ofstream log;
ofstream log_latex;
log.open (outpath+"/Purity_log"+xn+".txt");
log_latex.open (outpath+"/Purity_log_latex"+xn+".txt");
Int_t fillcolor = 3;
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
gStyle->SetFrameLineWidth(3);
gStyle->SetLineWidth(2);
TCanvas* canvas = new TCanvas("canvas","canvas",900,100,500,500);
gStyle->SetOptStat(0);
gPad->SetTickx();
gPad->SetTicky();
gStyle->SetLineWidth(3);
TText* title = new TText(1,1,"") ;
title->SetTextSize(0.07);
title->SetTextColor(kBlack);
title->SetTextAlign(13);
title->SetTextFont(62);
TText* extra = new TText(1,1,"") ;
extra->SetTextSize(0.05);
extra->SetTextColor(kBlack);
extra->SetTextAlign(13);
extra->SetTextFont(52);
TText* lumi = new TText(1,1,"") ;
lumi->SetTextSize(0.05);
lumi->SetTextColor(kBlack);
lumi->SetTextAlign(31);
lumi->SetTextFont(42);
int c_sig = 4;
int c_bkg = 2;
int c_den = 8;
int ls_cmb = 1;
int ls_gj = 2;
int ls_qcd = 3;
//for(std::vector<TString>::iterator it = ptranges.begin(); it != ptranges.end(); ++it) {
for(unsigned i=0; i<ptranges.size(); i++) {
for(unsigned j=0; j<isovars.size(); j++) {
for(unsigned k=0; k<cuts.size(); k++) {
for(unsigned l=0; l<rebins.size(); l++) {
TString ptrange = ptranges.at(i);
TString isovar = isovars.at(j);
TString cut = cuts.at(k);
TString rebin = rebins.at(l);
TString cutname;
if(cut=="idnc"){cutname="ID";}
if(cut=="idnc_mL30"){cutname="ID+MET";}
if(cut=="idnc_trig"){cutname="ID+Trigger";}
if(cut=="idnc_mL30_trig"){cutname="ID+MET+Trigger";}
TString bwidth;
if(rebin=="rebin1"){bwidth="width=0.5";}
if(rebin=="rebin2"){bwidth="width=1.0";}
if(rebin=="rebin3"){bwidth="width=1.5";}
if(rebin=="rebin4"){bwidth="width=2.0";}
if(rebin=="rebin5"){bwidth="width=2.5";}
std::cout<<boost::format("%8s %15s %7s %7s ") % ptrange % cut % isovar % rebin<<std::endl;
//Signal purity
TH1F* h_Nsig_GJ = (TH1F*)infile_GJ->Get("h_Nsig_"+isovar+"_"+ptrange+"_"+cut);
TH1F* h_Nsig_QCD = (TH1F*)infile_QCD->Get("h_Nsig_"+isovar+"_"+ptrange+"_"+cut);
h_Nsig_GJ->Rebin(3); h_Nsig_QCD->Rebin(3); //h_Nsig_GJ->Rebin(l+1); h_Nsig_QCD->Rebin(l+1);
//TH1F* h_Nsig = (TH1F*)h_Nsig_QCD->Clone("h_Nsig");
TH1F* h_Nsig = (TH1F*)h_Nsig_GJ->Clone("h_Nsig");
h_Nsig->Add(h_Nsig_QCD);
h_Nsig->SetLineColor(c_sig); h_Nsig_GJ->SetLineColor(c_sig); h_Nsig_QCD->SetLineColor(c_sig);
h_Nsig->SetLineStyle(ls_cmb); h_Nsig_GJ->SetLineStyle(ls_gj); h_Nsig_QCD->SetLineStyle(ls_qcd);
h_Nsig->SetLineWidth(2); h_Nsig_GJ->SetLineWidth(2); h_Nsig_QCD->SetLineWidth(2);
//Background purity
TH1F* h_Nbkg_GJ = (TH1F*)infile_GJ->Get("h_Nbkg_"+isovar+"_"+ptrange+"_"+cut);
TH1F* h_Nbkg_QCD = (TH1F*)infile_QCD->Get("h_Nbkg_"+isovar+"_"+ptrange+"_"+cut);
h_Nbkg_GJ->Rebin(3); h_Nbkg_QCD->Rebin(3);
//TH1F* h_Nbkg = (TH1F*)h_Nbkg_QCD->Clone("h_Nbkg");
TH1F* h_Nbkg = (TH1F*)h_Nbkg_GJ->Clone("h_Nbkg");
h_Nbkg->Add(h_Nbkg_QCD);
h_Nbkg->SetLineColor(c_bkg); h_Nbkg_GJ->SetLineColor(c_bkg); h_Nbkg_QCD->SetLineColor(c_bkg);
h_Nbkg->SetLineStyle(ls_cmb); h_Nbkg_GJ->SetLineStyle(ls_gj); h_Nbkg_QCD->SetLineStyle(ls_qcd);
h_Nbkg->SetLineWidth(2); h_Nbkg_GJ->SetLineWidth(2); h_Nbkg_QCD->SetLineWidth(2);
//Denominator
TH1F* h_Deno_GJ = (TH1F*)h_Nsig_GJ->Clone("h_Deno_GJ");
h_Deno_GJ->Add(h_Nbkg_GJ);
TH1F* h_Deno_QCD = (TH1F*)h_Nsig_QCD->Clone("h_Deno_QCD");
h_Deno_QCD->Add(h_Nbkg_QCD);
// TH1F* h_Deno_GJ = (TH1F*)infile_GJ->Get("h_Deno_"+isovar+"_"+ptrange+"_"+cut);
// TH1F* h_Deno_QCD = (TH1F*)infile_QCD->Get("h_Deno_"+isovar+"_"+ptrange+"_"+cut);
// h_Deno_GJ->Rebin(3); h_Deno_QCD->Rebin(3);
// //TH1F* h_Deno = (TH1F*)h_Deno_QCD->Clone("h_Deno");
TH1F* h_Deno = (TH1F*)h_Deno_GJ->Clone("h_Deno");
h_Deno->Add(h_Deno_QCD);
h_Deno->SetLineColor(c_den); h_Deno_GJ->SetLineColor(c_den); h_Deno_QCD->SetLineColor(c_den);
h_Deno->SetLineStyle(ls_cmb); h_Deno_GJ->SetLineStyle(ls_gj); h_Deno_QCD->SetLineStyle(ls_qcd);
h_Deno->SetLineWidth(1); h_Deno_GJ->SetLineWidth(1); h_Deno_QCD->SetLineWidth(1);
// should probably add these as last bin..
h_Nsig->ClearUnderflowAndOverflow();
h_Nbkg->ClearUnderflowAndOverflow();
h_Deno->ClearUnderflowAndOverflow();
//Define purity
TGraphAsymmErrors *purity_signal = new TGraphAsymmErrors(h_Nsig,h_Deno,"n");
purity_signal->SetLineWidth(2);
purity_signal->SetLineColor(4);
purity_signal->SetMarkerColor(4);
purity_signal->SetMarkerStyle(20);
TGraphAsymmErrors *purity_background = new TGraphAsymmErrors(h_Nbkg,h_Deno,"n");
purity_background->SetLineWidth(2);
purity_background->SetLineColor(2);
purity_background->SetMarkerColor(2);
purity_background->SetMarkerStyle(24);
//-----------------------------------------------------
// Calculations
log<<boost::format("%8s %15s %7s %7s ") % ptrange % cut % isovar % rebin ;
// Start lower bound calculation
Int_t nbins;
nbins = purity_background->GetN();
log<<boost::format("(%2i bins) \n") % nbins ;
std::vector<Double_t> thresholds, xsatthresh, ysatthresh, xsatqt, avgpur;
std::vector<Int_t> batthresh, batqt;
thresholds.clear(); // purity thresholds
xsatthresh.clear(); // x st. f(x) > thresh
ysatthresh.clear(); // y = f(x)x at thresh
batthresh.clear(); // bin corresponding to x
batqt.clear(); // bin at upper threshold
xsatqt.clear(); // x value at upper thresh
avgpur.clear(); // average purity over range
thresholds.push_back(0.50);
thresholds.push_back(0.55);
thresholds.push_back(0.60);
thresholds.push_back(0.65);
thresholds.push_back(0.70);
Int_t nthreshs;
nthreshs = thresholds.size();
// find purity lower bounds = x st. f(x)>thresh
Double_t pointx, pointy;
for(int m=0; m<nbins; ++m){
purity_background->GetPoint(m, pointx, pointy);
//log<<boost::format(" %2i %2.2f %2.2f \n") % i % pointx % pointy;
for(int n=0; n<nthreshs; ++n){
if(pointy > thresholds.at(n) && ysatthresh.size()==n ){
ysatthresh.push_back(pointy);
xsatthresh.push_back(pointx);
batthresh.push_back(h_Nsig->FindBin(pointx));
}
}
}
// for thresholds never reached, fill with content of xmax
purity_background->GetPoint(nbins-1,pointx,pointy);
Int_t nrnpassthresh = 0;
nrnpassthresh = nthreshs - ysatthresh.size();
for(int n=0; n<nrnpassthresh; ++n){
ysatthresh.push_back(pointy);
xsatthresh.push_back(pointx);
batthresh.push_back(h_Nsig->FindBin(pointx));
}
//// dump lower bound threshold info
//for(int n=0; n<nthreshs; ++n){
// log<<boost::format(" Threshold: %.2f %2i (%2.2f, %0.2f)\n")
// % thresholds.at(n) % batthresh.at(n) % xsatthresh.at(n) % ysatthresh.at(n) ;
//}
log<<boost::format(" Threshold bin (xx.x,yy.y) TotalSize above|below\n");
// Start upper bound calculation
Double_t fullsize = h_Nbkg->Integral(0,-1);
for(int n=0; n<nthreshs; ++n){
Double_t sizedown = h_Nbkg->Integral(0,batthresh.at(n)-1);
Double_t sizeup = h_Nbkg->Integral(batthresh.at(n),-1);
log<<boost::format(" %7.2f %3i (%2.1f,%1.2f) %4.1f %6.1f|%6.1f \n")
% thresholds.at(n) % batthresh.at(n) % xsatthresh.at(n) % ysatthresh.at(n)
% fullsize % sizeup % sizedown ;
}
// upper bound calculation
// keep adding contents of bins starting at threshbin
// until sum > fraction (1/4) of full integral
Double_t goalsize = fullsize / 4.;
for(int n=0; n<nthreshs; ++n){
Double_t thissize = 0.;
int thisbin = batthresh.at(n);
while( thisbin < nbins+1 ){
thissize += h_Nbkg->GetBinContent(thisbin);
if(thissize > goalsize){
break;
}
thisbin++;
}
batqt.push_back(thisbin);
xsatqt.push_back(h_Nbkg->GetBinCenter(thisbin));
}
log<<" Target Eventcount: "<<goalsize<<" = fullsize/4 = "<<fullsize<<"/4 \n";
log<<boost::format(" Threshold blo bhi (xlo,xhi) Int(blo,bhi) Int(blo,bhi+1) Int(blo,bhi-1)\n");
// find average purity over range
// must weight each purity bin by nr. events in bin
for(int n=0; n<nthreshs; ++n){
int nbins = batqt.at(n) - batthresh.at(n);
//std::cout<<nbins<<" "<<batthresh.at(n)<<" "<<batqt.at(n)<<std::endl;
Double_t sumpurs, sumweights;
sumpurs = 0.;
sumweights = 0.;
for(int o=0; o<nbins+1; ++o){
int thisbin = batthresh.at(n) + o;
int thispoint = thisbin - 1; // bin nr = point nr + 1
purity_background->GetPoint(thispoint, pointx,pointy);
sumpurs += pointy*h_Nbkg->GetBinContent(thisbin);
sumweights += h_Nbkg->GetBinContent(thisbin);
//std::cout<<" "<<thispoint<<"|"<<pointy<<"|"<<h_Nbkg->GetBinContent(thisbin)<<" "<<std::endl;;
}
//std::cout<<std::endl;
//std::cout<<sumpurs/sumweights<<std::endl<<std::endl;;
avgpur.push_back( sumpurs/sumweights );
}
log_latex<<"\\begin{tabular}{r|r|r|r|r|r|r|r}\n";
log_latex<<boost::format("\\multicolumn{2}{c|}{%8s} & \\multicolumn{2}{c|}{%15s} & \\multicolumn{2}{c|}{%7s} & \\multicolumn{2}{c|}{%7s} \\\\ \\hline \\hline \n") % ptrange % cutname % isovar % bwidth ;
log_latex<<boost::format(" %10s & %5s & %5s & %10s & %5s & %8s & %15s & %10s \\\\ \\hline \n")
% "Threshold"
% "x lo" % "x hi"
% "avg purity"
% "err \\%"
% "evts. tot."
% "evts. > x lo"
% "evts in rng." ;
///log_latex<<"\\begin{tabular}{r|r|r|r|r}\n";
///log_latex<<boost::format("%8s & \\multicolumn{2}{c}{%15s} & %7s & %7s \\\\ \n \\hline \\hline \n") % ptrange % cutname % isovar % bwidth ;
///log_latex<<boost::format(" %10s & %5s & %5s & %10s & %5s \\\\ \n \\hline \n")
/// % "Threshold"
/// % "x lo" % "x hi"
/// % "avg purity"
/// % "err \\%" ;
for(int n=0; n<nthreshs; ++n){
log<<boost::format(" %7.2f %3i %3i (%2.1f,%1.2f) %10.2f %10.2f %10.2f \n")
% thresholds.at(n) % batthresh.at(n) %batqt.at(n)
% xsatthresh.at(n) % xsatqt.at(n)
% h_Nbkg->Integral(batthresh.at(n),batqt.at(n))
% h_Nbkg->Integral(batthresh.at(n),batqt.at(n)+1)
% h_Nbkg->Integral(batthresh.at(n),batqt.at(n)-1) ;
Double_t nrevents, errevents, pcterr;
nrevents = h_Nbkg->IntegralAndError(batthresh.at(n),batqt.at(n),errevents);
pcterr = errevents/nrevents;
log_latex<<boost::format(" %10.2f & %5.2f & %5.2f & %10.4f & %5.2f & %8.1f & %15.1f & %10.1f \\\\ \n")
% thresholds.at(n)
% h_Nbkg->GetBinLowEdge(h_Nbkg->FindBin(xsatthresh.at(n))) % h_Nbkg->GetBinLowEdge(h_Nbkg->FindBin(xsatqt.at(n) ))
//% xsatthresh.at(n) % xsatqt.at(n)
% avgpur.at(n)
% (pcterr*100)
% fullsize
% h_Nbkg->Integral(batthresh.at(n),-1)
% h_Nbkg->Integral(batthresh.at(n),batqt.at(n));
///log_latex<<boost::format(" %10.2f & %5.2f & %5.2f & %7.3f & %5.2f \\\\ \n")
/// % thresholds.at(n)
/// % xsatthresh.at(n) % xsatqt.at(n)
/// % avgpur.at(n)
/// % (pcterr*100) ;
}
log_latex<<"\\hline \n \\end{tabular}\n";
//-----------------------------------------------------
// all set, now to start thinking about plotting
//Make lower bound lines (purity)
Double_t frac = 1./nthreshs;
TLine *lgb_lo0 = new TLine(xsatthresh.at(0), 0*frac, xsatthresh.at(0), 1*frac);
TLine *lgb_lo1 = new TLine(xsatthresh.at(1), 1*frac, xsatthresh.at(1), 2*frac);
TLine *lgb_lo2 = new TLine(xsatthresh.at(2), 2*frac, xsatthresh.at(2), 3*frac);
TLine *lgb_lo3 = new TLine(xsatthresh.at(3), 3*frac, xsatthresh.at(3), 4*frac);
TLine *lgb_lo4 = new TLine(xsatthresh.at(4), 4*frac, xsatthresh.at(4), 5*frac);
TLine *lgc_lo0 = new TLine(xsatthresh.at(0), 0*frac, xsatthresh.at(0), 1*frac);
TLine *lgc_lo1 = new TLine(xsatthresh.at(1), 1*frac, xsatthresh.at(1), 2*frac);
TLine *lgc_lo2 = new TLine(xsatthresh.at(2), 2*frac, xsatthresh.at(2), 3*frac);
TLine *lgc_lo3 = new TLine(xsatthresh.at(3), 3*frac, xsatthresh.at(3), 4*frac);
TLine *lgc_lo4 = new TLine(xsatthresh.at(4), 4*frac, xsatthresh.at(4), 5*frac);
lgb_lo0->SetLineColor(1);
lgb_lo1->SetLineColor(1);
lgb_lo2->SetLineColor(1);
lgb_lo3->SetLineColor(1);
lgb_lo4->SetLineColor(1);
lgc_lo0->SetLineColor(2);
lgc_lo1->SetLineColor(3);
lgc_lo2->SetLineColor(4);
lgc_lo3->SetLineColor(5);
lgc_lo4->SetLineColor(6);
lgb_lo0->SetLineWidth(4);
lgb_lo1->SetLineWidth(4);
lgb_lo2->SetLineWidth(4);
lgb_lo3->SetLineWidth(4);
lgb_lo4->SetLineWidth(4);
lgc_lo0->SetLineWidth(2);
lgc_lo1->SetLineWidth(2);
lgc_lo2->SetLineWidth(2);
lgc_lo3->SetLineWidth(2);
lgc_lo4->SetLineWidth(2);
//Make upper bound lines (purity)
TLine *lgb_hi0 = new TLine(xsatqt.at(0), 0*frac, xsatqt.at(0), 1*frac);
TLine *lgb_hi1 = new TLine(xsatqt.at(1), 1*frac, xsatqt.at(1), 2*frac);
TLine *lgb_hi2 = new TLine(xsatqt.at(2), 2*frac, xsatqt.at(2), 3*frac);
TLine *lgb_hi3 = new TLine(xsatqt.at(3), 3*frac, xsatqt.at(3), 4*frac);
TLine *lgb_hi4 = new TLine(xsatqt.at(4), 4*frac, xsatqt.at(4), 5*frac);
TLine *lgc_hi0 = new TLine(xsatqt.at(0), 0*frac, xsatqt.at(0), 1*frac);
TLine *lgc_hi1 = new TLine(xsatqt.at(1), 1*frac, xsatqt.at(1), 2*frac);
TLine *lgc_hi2 = new TLine(xsatqt.at(2), 2*frac, xsatqt.at(2), 3*frac);
TLine *lgc_hi3 = new TLine(xsatqt.at(3), 3*frac, xsatqt.at(3), 4*frac);
TLine *lgc_hi4 = new TLine(xsatqt.at(4), 4*frac, xsatqt.at(4), 5*frac);
lgb_hi0->SetLineColor(1);
lgb_hi1->SetLineColor(1);
lgb_hi2->SetLineColor(1);
lgb_hi3->SetLineColor(1);
lgb_hi4->SetLineColor(1);
lgc_hi0->SetLineColor(2);
lgc_hi1->SetLineColor(3);
lgc_hi2->SetLineColor(4);
lgc_hi3->SetLineColor(5);
lgc_hi4->SetLineColor(6);
lgb_hi0->SetLineWidth(4);
lgb_hi1->SetLineWidth(4);
lgb_hi2->SetLineWidth(4);
lgb_hi3->SetLineWidth(4);
lgb_hi4->SetLineWidth(4);
lgc_hi0->SetLineWidth(2);
lgc_hi1->SetLineWidth(2);
lgc_hi2->SetLineWidth(2);
lgc_hi3->SetLineWidth(2);
lgc_hi4->SetLineWidth(2);
lgc_hi0->SetLineStyle(2);
lgc_hi1->SetLineStyle(2);
lgc_hi2->SetLineStyle(2);
lgc_hi3->SetLineStyle(2);
lgc_hi4->SetLineStyle(2);
//Make avg purity horizontal lines
TLine *lgb_pur0 = new TLine(xsatthresh.at(0), avgpur.at(0), xsatqt.at(0), avgpur.at(0));
TLine *lgb_pur1 = new TLine(xsatthresh.at(1), avgpur.at(1), xsatqt.at(1), avgpur.at(1));
TLine *lgb_pur2 = new TLine(xsatthresh.at(2), avgpur.at(2), xsatqt.at(2), avgpur.at(2));
TLine *lgb_pur3 = new TLine(xsatthresh.at(3), avgpur.at(3), xsatqt.at(3), avgpur.at(3));
TLine *lgb_pur4 = new TLine(xsatthresh.at(4), avgpur.at(4), xsatqt.at(4), avgpur.at(4));
TLine *lgc_pur0 = new TLine(xsatthresh.at(0), avgpur.at(0), xsatqt.at(0), avgpur.at(0));
TLine *lgc_pur1 = new TLine(xsatthresh.at(1), avgpur.at(1), xsatqt.at(1), avgpur.at(1));
TLine *lgc_pur2 = new TLine(xsatthresh.at(2), avgpur.at(2), xsatqt.at(2), avgpur.at(2));
TLine *lgc_pur3 = new TLine(xsatthresh.at(3), avgpur.at(3), xsatqt.at(3), avgpur.at(3));
TLine *lgc_pur4 = new TLine(xsatthresh.at(4), avgpur.at(4), xsatqt.at(4), avgpur.at(4));
lgb_pur0->SetLineColor(1);
lgb_pur1->SetLineColor(1);
lgb_pur2->SetLineColor(1);
lgb_pur3->SetLineColor(1);
lgb_pur4->SetLineColor(1);
//lgb_pur0->SetLineColor(2);
//lgb_pur1->SetLineColor(3);
//lgb_pur2->SetLineColor(4);
//lgb_pur3->SetLineColor(5);
//lgb_pur4->SetLineColor(6);
lgc_pur0->SetLineColor(2);
lgc_pur1->SetLineColor(3);
lgc_pur2->SetLineColor(4);
lgc_pur3->SetLineColor(5);
lgc_pur4->SetLineColor(6);
lgb_pur0->SetLineWidth(4);
lgb_pur1->SetLineWidth(4);
lgb_pur2->SetLineWidth(4);
lgb_pur3->SetLineWidth(4);
lgb_pur4->SetLineWidth(4);
lgc_pur0->SetLineWidth(2);
lgc_pur1->SetLineWidth(2);
lgc_pur2->SetLineWidth(2);
lgc_pur3->SetLineWidth(2);
lgc_pur4->SetLineWidth(2);
//lgc_pur0->SetLineStyle(2);
//lgc_pur1->SetLineStyle(2);
//lgc_pur2->SetLineStyle(2);
//lgc_pur3->SetLineStyle(2);
//lgc_pur4->SetLineStyle(2);
//-----------------------------------------------------
Double_t themax = 0 ;
themax = std::max(themax, h_Nsig->GetMaximum());
themax = std::max(themax, h_Nbkg->GetMaximum());
themax = std::max(themax, h_Deno->GetMaximum());
//Make lower bound lines (hist)
Double_t logmax;
Double_t logmin;
logmax = std::log(themax);
logmin = std::log(1);
Double_t spacing = (logmax-logmin)/nthreshs;
//std::cout<<spacing<<std::endl;
//Double_t fifmax = themax/5.;
TLine *lhb_lo0 = new TLine(xsatthresh.at(0), exp(logmin+0*spacing), xsatthresh.at(0), exp(logmin+1*spacing));
TLine *lhb_lo1 = new TLine(xsatthresh.at(1), exp(logmin+1*spacing), xsatthresh.at(1), exp(logmin+2*spacing));
TLine *lhb_lo2 = new TLine(xsatthresh.at(2), exp(logmin+2*spacing), xsatthresh.at(2), exp(logmin+3*spacing));
TLine *lhb_lo3 = new TLine(xsatthresh.at(3), exp(logmin+3*spacing), xsatthresh.at(3), exp(logmin+4*spacing));
TLine *lhb_lo4 = new TLine(xsatthresh.at(4), exp(logmin+4*spacing), xsatthresh.at(4), exp(logmin+5*spacing));
TLine *lhc_lo0 = new TLine(xsatthresh.at(0), exp(logmin+0*spacing), xsatthresh.at(0), exp(logmin+1*spacing));
TLine *lhc_lo1 = new TLine(xsatthresh.at(1), exp(logmin+1*spacing), xsatthresh.at(1), exp(logmin+2*spacing));
TLine *lhc_lo2 = new TLine(xsatthresh.at(2), exp(logmin+2*spacing), xsatthresh.at(2), exp(logmin+3*spacing));
TLine *lhc_lo3 = new TLine(xsatthresh.at(3), exp(logmin+3*spacing), xsatthresh.at(3), exp(logmin+4*spacing));
TLine *lhc_lo4 = new TLine(xsatthresh.at(4), exp(logmin+4*spacing), xsatthresh.at(4), exp(logmin+5*spacing));
lhb_lo0->SetLineColor(1);
lhb_lo1->SetLineColor(1);
lhb_lo2->SetLineColor(1);
lhb_lo3->SetLineColor(1);
lhb_lo4->SetLineColor(1);
lhc_lo0->SetLineColor(2);
lhc_lo1->SetLineColor(3);
lhc_lo2->SetLineColor(4);
lhc_lo3->SetLineColor(5);
lhc_lo4->SetLineColor(6);
lhb_lo0->SetLineWidth(4);
lhb_lo1->SetLineWidth(4);
lhb_lo2->SetLineWidth(4);
lhb_lo3->SetLineWidth(4);
lhb_lo4->SetLineWidth(4);
lhc_lo0->SetLineWidth(2);
lhc_lo1->SetLineWidth(2);
lhc_lo2->SetLineWidth(2);
lhc_lo3->SetLineWidth(2);
lhc_lo4->SetLineWidth(2);
//Make upper bound lines (hist)
TLine *lhb_hi0 = new TLine(xsatqt.at(0), exp(logmin+0*spacing), xsatqt.at(0), exp(logmin+1*spacing));
TLine *lhb_hi1 = new TLine(xsatqt.at(1), exp(logmin+1*spacing), xsatqt.at(1), exp(logmin+2*spacing));
TLine *lhb_hi2 = new TLine(xsatqt.at(2), exp(logmin+2*spacing), xsatqt.at(2), exp(logmin+3*spacing));
TLine *lhb_hi3 = new TLine(xsatqt.at(3), exp(logmin+3*spacing), xsatqt.at(3), exp(logmin+4*spacing));
TLine *lhb_hi4 = new TLine(xsatqt.at(4), exp(logmin+4*spacing), xsatqt.at(4), exp(logmin+5*spacing));
TLine *lhc_hi0 = new TLine(xsatqt.at(0), exp(logmin+0*spacing), xsatqt.at(0), exp(logmin+1*spacing));
TLine *lhc_hi1 = new TLine(xsatqt.at(1), exp(logmin+1*spacing), xsatqt.at(1), exp(logmin+2*spacing));
TLine *lhc_hi2 = new TLine(xsatqt.at(2), exp(logmin+2*spacing), xsatqt.at(2), exp(logmin+3*spacing));
TLine *lhc_hi3 = new TLine(xsatqt.at(3), exp(logmin+3*spacing), xsatqt.at(3), exp(logmin+4*spacing));
TLine *lhc_hi4 = new TLine(xsatqt.at(4), exp(logmin+4*spacing), xsatqt.at(4), exp(logmin+5*spacing));
lhb_hi0->SetLineColor(1);
lhb_hi1->SetLineColor(1);
lhb_hi2->SetLineColor(1);
lhb_hi3->SetLineColor(1);
lhb_hi4->SetLineColor(1);
lhc_hi0->SetLineColor(2);
lhc_hi1->SetLineColor(3);
lhc_hi2->SetLineColor(4);
lhc_hi3->SetLineColor(5);
lhc_hi4->SetLineColor(6);
lhb_hi0->SetLineWidth(4);
lhb_hi1->SetLineWidth(4);
lhb_hi2->SetLineWidth(4);
lhb_hi3->SetLineWidth(4);
lhb_hi4->SetLineWidth(4);
lhc_hi0->SetLineWidth(2);
lhc_hi1->SetLineWidth(2);
lhc_hi2->SetLineWidth(2);
lhc_hi3->SetLineWidth(2);
lhc_hi4->SetLineWidth(2);
lhc_hi0->SetLineStyle(2);
lhc_hi1->SetLineStyle(2);
lhc_hi2->SetLineStyle(2);
lhc_hi3->SetLineStyle(2);
lhc_hi4->SetLineStyle(2);
//-----------------------------------------------------
//-----------------------------------------------------
//Draw Histograms
TLegend *meg = new TLegend(0.4,0.6,0.88,0.88);
//meg->SetTextSize(0.04) ;
meg->SetFillColor(0); meg->SetShadowColor(0);meg->SetBorderSize(2);
meg->SetTextFont(22.);
meg->AddEntry(h_Nsig,"Signal","l");
meg->AddEntry(h_Nsig_GJ,"Signal (GJ MC)","l");
meg->AddEntry(h_Nsig_QCD,"Signal (QCD MC)","l");
meg->AddEntry(h_Nbkg,"Background","l");
meg->AddEntry(h_Nbkg_GJ,"Background (GJ MC)","l");
meg->AddEntry(h_Nbkg_QCD,"Background (QCD MC)","l");
meg->AddEntry(h_Deno,"Denominator","l");
meg->AddEntry(h_Deno_GJ,"Denominator (GJ MC)","l");
meg->AddEntry(h_Deno_QCD,"Denominator (QCD MC)","l");
TString leg_h0; leg_h0.Form("Purity > %0.2f (%2.1f-%2.1f)", thresholds.at(0), xsatthresh.at(0), xsatqt.at(0));
TString leg_h1; leg_h1.Form("Purity > %0.2f (%2.1f-%2.1f)", thresholds.at(1), xsatthresh.at(1), xsatqt.at(1));
TString leg_h2; leg_h2.Form("Purity > %0.2f (%2.1f-%2.1f)", thresholds.at(2), xsatthresh.at(2), xsatqt.at(2));
TString leg_h3; leg_h3.Form("Purity > %0.2f (%2.1f-%2.1f)", thresholds.at(3), xsatthresh.at(3), xsatqt.at(3));
TString leg_h4; leg_h4.Form("Purity > %0.2f (%2.1f-%2.1f)", thresholds.at(4), xsatthresh.at(4), xsatqt.at(4));
TLegend *neg = new TLegend(0.15,0.5,0.4,0.75);
//meg->SetTextSize(0.04) ;
neg->SetFillColor(0); neg->SetShadowColor(0);neg->SetBorderSize(2);
neg->SetTextFont(22.);
neg->AddEntry(lhc_lo0,leg_h0,"l");
neg->AddEntry(lhc_lo1,leg_h1,"l");
neg->AddEntry(lhc_lo2,leg_h2,"l");
neg->AddEntry(lhc_lo3,leg_h3,"l");
neg->AddEntry(lhc_lo4,leg_h4,"l");
h_Nsig->SetTitle(cut);
h_Nsig->SetXTitle(h_Nsig->GetTitle());
h_Nsig->SetYTitle("Events ");
h_Nsig->SetMaximum(1000*themax);
h_Nsig->SetMinimum(1);
h_Nsig->Draw("hist");
h_Nsig_GJ->Draw("hist,same");
h_Nsig_QCD->Draw("hist,same");
h_Nbkg->Draw("hist,same");
h_Nbkg_GJ->Draw("hist,same");
h_Nbkg_QCD->Draw("hist,same");
h_Deno->Draw("hist,same");
h_Deno_GJ->Draw("hist,same");
h_Deno_QCD->Draw("hist,same");
meg->SetHeader("photon p_{T}: "+ptrange);
meg->Draw();
if(dolines){
neg->Draw();
}
title->DrawTextNDC(0.17,0.87,"CMS");
extra->DrawTextNDC(0.17,0.81,"Preliminary");
lumi->DrawTextNDC(0.9,0.91,"2.32 /fb (13 TeV)");
if(dolines){
lhb_lo0->Draw();
lhb_lo1->Draw();
lhb_lo2->Draw();
lhb_lo3->Draw();
lhb_lo4->Draw();
lhb_hi0->Draw();
lhb_hi1->Draw();
lhb_hi2->Draw();
lhb_hi3->Draw();
lhb_hi4->Draw();
lhc_lo0->Draw();
lhc_lo1->Draw();
lhc_lo2->Draw();
lhc_lo3->Draw();
lhc_lo4->Draw();
lhc_hi0->Draw();
lhc_hi1->Draw();
lhc_hi2->Draw();
lhc_hi3->Draw();
lhc_hi4->Draw();
}
gPad->SetLogy();
canvas->SaveAs(outpath+"/Histos_"+isovar+"_"+ptrange+"_"+cut+rebin+xn+".pdf");
canvas->SaveAs(wwwpath+"/pdf/"+ptrange+"/"+isovar+"/"+cut+lines+"/Histos_"+isovar+"_"+ptrange+"_"+cut+rebin+xn+".pdf");
canvas->SaveAs(wwwpath+"/png/"+ptrange+"/"+isovar+"/"+cut+lines+"/Histos_"+isovar+"_"+ptrange+"_"+cut+rebin+xn+".png");
canvas->Clear();
gPad->SetLogy(kFALSE);
//-----------------------------------------------------
//-----------------------------------------------------
//Draw Purity
TLegend *leg = new TLegend(0.5,0.75,0.88,0.88);
leg->SetTextSize(0.04) ;
leg->SetFillColor(0); leg->SetShadowColor(0);leg->SetBorderSize(0);
leg->SetTextFont(22.);
leg->AddEntry(purity_signal,"Signal","lep");
leg->AddEntry(purity_background,"Background","lep");
float xmin ; //= 0.;
float xmax ; //= 25.;
xmin = h_Nsig->GetXaxis()->GetXmin();
xmax = h_Nsig->GetXaxis()->GetXmax();
TLine *lowline = new TLine(xmin,0,xmax,0);
TLine *hiline = new TLine(xmin,1,xmax,1);
lowline->SetLineWidth(1); lowline->SetLineColor(1); lowline->SetLineStyle(3);
hiline->SetLineWidth(1); hiline->SetLineColor(1); hiline->SetLineStyle(3);
TH1F *hframe = canvas->DrawFrame(xmin,-0.05,xmax,1.4,"");
lowline->Draw();
hiline->Draw();
title->DrawTextNDC(0.17,0.87,"CMS");
extra->DrawTextNDC(0.17,0.81,"Preliminary");
lumi->DrawTextNDC(0.9,0.91,"2.32 /fb (13 TeV)");
hframe->SetTitle(cut);
hframe->SetXTitle(isovar);
hframe->SetYTitle("Purity");
// draw points
purity_signal->Draw("P");
purity_background->Draw("P");
leg->SetHeader("photon p_{T}: "+ptrange);
leg->Draw();
//canvas->SaveAs(outpath+"/Purity_"+isovar+"_"+ptrange+"_"+cut+"_"+rebin+xn+".pdf");
if(dolines){
lgb_lo0->Draw();
lgb_lo1->Draw();
lgb_lo2->Draw();
lgb_lo3->Draw();
lgb_lo4->Draw();
lgb_hi0->Draw();
lgb_hi1->Draw();
lgb_hi2->Draw();
lgb_hi3->Draw();
lgb_hi4->Draw();
lgc_lo0->Draw();
lgc_lo1->Draw();
lgc_lo2->Draw();
lgc_lo3->Draw();
lgc_lo4->Draw();
lgc_hi0->Draw();
lgc_hi1->Draw();
lgc_hi2->Draw();
lgc_hi3->Draw();
lgc_hi4->Draw();
lgb_pur0->Draw();
lgb_pur1->Draw();
lgb_pur2->Draw();
lgb_pur3->Draw();
lgb_pur4->Draw();
lgc_pur0->Draw();
lgc_pur1->Draw();
lgc_pur2->Draw();
lgc_pur3->Draw();
lgc_pur4->Draw();
}
canvas->SaveAs(outpath+"/Purity_"+isovar+"_"+ptrange+"_"+cut+rebin+xn+".pdf");
canvas->SaveAs(wwwpath+"/pdf/"+ptrange+"/"+isovar+"/"+cut+lines+"/Purity_"+isovar+"_"+ptrange+"_"+cut+rebin+xn+".pdf");
canvas->SaveAs(wwwpath+"/png/"+ptrange+"/"+isovar+"/"+cut+lines+"/Purity_"+isovar+"_"+ptrange+"_"+cut+rebin+xn+".png");
canvas->Clear();
//-----------------------------------------------------
// clean your room!
h_Nsig->Delete();
h_Nsig_GJ->Delete();
h_Nsig_QCD->Delete();
h_Nbkg->Delete();
h_Nbkg_GJ->Delete();
h_Nbkg_QCD->Delete();
h_Deno->Delete();
h_Deno_GJ->Delete();
h_Deno_QCD->Delete();
}
log<<"--------------------------------------\n";
log_latex<<"--------------------------------------\n";
}
}
}
log.close();
log_latex.close();
}
void EMCDistribution_ADC(bool log_scale = true)
{
TString gain = "RAW";
TText *t;
TCanvas *c1 = new TCanvas(
"EMCDistribution_ADC_" + gain + TString(log_scale ? "_Log" : "") + cuts,
"EMCDistribution_ADC_" + gain + TString(log_scale ? "_Log" : "") + cuts,
1800, 1000);
c1->Divide(8, 8, 0., 0.01);
int idx = 1;
TPad *p;
for (int iphi = 8 - 1; iphi >= 0; iphi--)
{
for (int ieta = 0; ieta < 8; ieta++)
{
p = (TPad *) c1->cd(idx++);
c1->Update();
if (log_scale)
{
p->SetLogz();
}
p->SetGridx(0);
p->SetGridy(0);
TString hname = Form("hEnergy_ieta%d_iphi%d", ieta, iphi) + TString(log_scale ? "_Log" : "");
TH1 *h = NULL;
if (log_scale)
h = new TH2F(hname,
Form(";Calibrated Tower Energy Sum (GeV);Count / bin"), 24, -.5,
23.5,
// 128+64, 0, 3096);
4098, -1, 4097);
// else
// h = new TH2F(hname,
// Form(";Calibrated Tower Energy Sum (GeV);Count / bin"), 100,
// -.050, .5,128,0,2048);
h->SetLineWidth(0);
h->SetLineColor(kBlue + 3);
h->SetFillColor(kBlue + 3);
h->GetXaxis()->SetTitleSize(.09);
h->GetXaxis()->SetLabelSize(.08);
h->GetYaxis()->SetLabelSize(.08);
h->GetYaxis()->SetRangeUser(0, 4096);
// if (log_scale)
// QAHistManagerDef::useLogBins(h->GetYaxis());
TString sdraw = "TOWER_" + gain + "_CEMC[].signal_samples[]:fmod(Iteration$,24)>>" + hname;
TString scut =
Form(
"TOWER_%s_CEMC[].get_bineta()==%d && TOWER_%s_CEMC[].get_binphi()==%d",
gain.Data(), ieta, gain.Data(), iphi);
cout << "T->Draw(\"" << sdraw << "\",\"" << scut << "\");" << endl;
T->Draw(sdraw, scut, "colz");
TText *t = new TText(.9, .9, Form("Col%d Row%d", ieta, iphi));
t->SetTextAlign(33);
t->SetTextSize(.15);
t->SetNDC();
t->Draw();
// return;
}
}
SaveCanvas(c1,
TString(_file0->GetName()) + TString("_DrawPrototype3EMCalTower_") + TString(c1->GetName()), false);
}
// all users - please change the name of this file to lhcbStyle.C
// Commits to lhcbdocs svn of .C files are not allowed
void lhcbStyle(){
// define names for colours
Int_t black = 1;
Int_t red = 2;
Int_t green = 3;
Int_t blue = 4;
Int_t yellow = 5;
Int_t magenta= 6;
Int_t cyan = 7;
Int_t purple = 9;
////////////////////////////////////////////////////////////////////
// PURPOSE:
//
// This macro defines a standard style for (black-and-white)
// "publication quality" LHCb ROOT plots.
//
// USAGE:
//
// Include the lines
// gROOT->ProcessLine(".L lhcbstyle.C");
// lhcbStyle();
// at the beginning of your root macro.
//
// Example usage is given in myPlot.C
//
// COMMENTS:
//
// Font:
//
// The font is chosen to be 132, this is Times New Roman (like the text of
// your document) with precision 2.
//
// "Landscape histograms":
//
// The style here is designed for more or less square plots.
// For longer histograms, or canvas with many pads, adjustements are needed.
// For instance, for a canvas with 1x5 histograms:
// TCanvas* c1 = new TCanvas("c1", "L0 muons", 600, 800);
// c1->Divide(1,5);
// Adaptions like the following will be needed:
// gStyle->SetTickLength(0.05,"x");
// gStyle->SetTickLength(0.01,"y");
// gStyle->SetLabelSize(0.15,"x");
// gStyle->SetLabelSize(0.1,"y");
// gStyle->SetStatW(0.15);
// gStyle->SetStatH(0.5);
//
// Authors: Thomas Schietinger, Andrew Powell, Chris Parkes, Niels Tuning
// Maintained by Editorial board member (currently Niels)
///////////////////////////////////////////////////////////////////
// Use times new roman, precision 2
Int_t lhcbFont = 132; // Old LHCb style: 62;
// Line thickness
Double_t lhcbWidth = 2.00; // Old LHCb style: 3.00;
// Text size
Double_t lhcbTSize = 0.06;
// use plain black on white colors
gROOT->SetStyle("Plain");
TStyle *lhcbStyle= new TStyle("lhcbStyle","LHCb plots style");
//lhcbStyle->SetErrorX(0); // don't suppress the error bar along X
lhcbStyle->SetFillColor(1);
lhcbStyle->SetFillStyle(1001); // solid
lhcbStyle->SetFrameFillColor(0);
lhcbStyle->SetFrameBorderMode(0);
lhcbStyle->SetPadBorderMode(0);
lhcbStyle->SetPadColor(0);
lhcbStyle->SetCanvasBorderMode(0);
lhcbStyle->SetCanvasColor(0);
lhcbStyle->SetStatColor(0);
lhcbStyle->SetLegendBorderSize(0);
// If you want the usual gradient palette (blue -> red)
lhcbStyle->SetPalette(1);
// If you want colors that correspond to gray scale in black and white:
int colors[8] = {0,5,7,3,6,2,4,1};
lhcbStyle->SetPalette(8,colors);
// set the paper & margin sizes
lhcbStyle->SetPaperSize(20,26);
lhcbStyle->SetPadTopMargin(0.05);
lhcbStyle->SetPadRightMargin(0.05); // increase for colz plots
lhcbStyle->SetPadBottomMargin(0.16);
lhcbStyle->SetPadLeftMargin(0.14);
// use large fonts
lhcbStyle->SetTextFont(lhcbFont);
lhcbStyle->SetTextSize(lhcbTSize);
lhcbStyle->SetLabelFont(lhcbFont,"x");
lhcbStyle->SetLabelFont(lhcbFont,"y");
lhcbStyle->SetLabelFont(lhcbFont,"z");
lhcbStyle->SetLabelSize(lhcbTSize,"x");
lhcbStyle->SetLabelSize(lhcbTSize,"y");
lhcbStyle->SetLabelSize(lhcbTSize,"z");
lhcbStyle->SetTitleFont(lhcbFont);
lhcbStyle->SetTitleFont(lhcbFont,"x");
lhcbStyle->SetTitleFont(lhcbFont,"y");
lhcbStyle->SetTitleFont(lhcbFont,"z");
lhcbStyle->SetTitleSize(1.2*lhcbTSize,"x");
lhcbStyle->SetTitleSize(1.2*lhcbTSize,"y");
lhcbStyle->SetTitleSize(1.2*lhcbTSize,"z");
// use medium bold lines and thick markers
lhcbStyle->SetLineWidth(lhcbWidth);
lhcbStyle->SetFrameLineWidth(lhcbWidth);
lhcbStyle->SetHistLineWidth(lhcbWidth);
lhcbStyle->SetFuncWidth(lhcbWidth);
lhcbStyle->SetGridWidth(lhcbWidth);
lhcbStyle->SetLineStyleString(2,"[12 12]"); // postscript dashes
lhcbStyle->SetMarkerStyle(20);
lhcbStyle->SetMarkerSize(1.0);
// label offsets
lhcbStyle->SetLabelOffset(0.010,"X");
lhcbStyle->SetLabelOffset(0.010,"Y");
// by default, do not display histogram decorations:
//lhcbStyle->SetOptStat(0);
lhcbStyle->SetOptStat("emr"); // show only nent -e , mean - m , rms -r
// full opts at http://root.cern.ch/root/html/TStyle.html#TStyle:SetOptStat
lhcbStyle->SetStatFormat("6.3g"); // specified as c printf options
lhcbStyle->SetOptTitle(0);
lhcbStyle->SetOptFit(0);
//lhcbStyle->SetOptFit(1011); // order is probability, Chi2, errors, parameters
//titles
lhcbStyle->SetTitleOffset(0.95,"X");
lhcbStyle->SetTitleOffset(0.95,"Y");
lhcbStyle->SetTitleOffset(1.2,"Z");
lhcbStyle->SetTitleFillColor(0);
lhcbStyle->SetTitleStyle(0);
lhcbStyle->SetTitleBorderSize(0);
lhcbStyle->SetTitleFont(lhcbFont,"title");
lhcbStyle->SetTitleX(0.0);
lhcbStyle->SetTitleY(1.0);
lhcbStyle->SetTitleW(1.0);
lhcbStyle->SetTitleH(0.05);
// look of the statistics box:
lhcbStyle->SetStatBorderSize(0);
lhcbStyle->SetStatFont(lhcbFont);
lhcbStyle->SetStatFontSize(0.05);
lhcbStyle->SetStatX(0.9);
lhcbStyle->SetStatY(0.9);
lhcbStyle->SetStatW(0.25);
lhcbStyle->SetStatH(0.15);
// put tick marks on top and RHS of plots
lhcbStyle->SetPadTickX(1);
lhcbStyle->SetPadTickY(1);
// histogram divisions: only 5 in x to avoid label overlaps
lhcbStyle->SetNdivisions(505,"x");
lhcbStyle->SetNdivisions(510,"y");
gROOT->SetStyle("lhcbStyle");
gROOT->ForceStyle();
/*
// add LHCb label
lhcbName = new TPaveText(gStyle->GetPadLeftMargin() + 0.05,
0.87 - gStyle->GetPadTopMargin(),
gStyle->treeGetPadLeftMargin() + 0.20,
0.95 - gStyle->GetPadTopMargin(),
"BRNDC");
lhcbName->AddText("LHCb");
lhcbName->SetFillColor(0);
lhcbName->SetTextAlign(12);
lhcbName->SetBorderSize(0);
*/
TText *lhcbLabel = new TText();
lhcbLabel->SetTextFont(lhcbFont);
lhcbLabel->SetTextColor(1);
lhcbLabel->SetTextSize(lhcbTSize);
lhcbLabel->SetTextAlign(12);
TLatex *lhcbLatex = new TLatex();
lhcbLatex->SetTextFont(lhcbFont);
lhcbLatex->SetTextColor(1);
lhcbLatex->SetTextSize(lhcbTSize);
lhcbLatex->SetTextAlign(12);
//std::cout << "-------------------------" << std::endl;
//std::cout << "Set LHCb Style - Feb 2012" << std::endl;
//std::cout << "-------------------------" << std::endl;
}
void crossfeeds_nondiag(TString title,
TString bkgfile,
TString epsfile,
TString txtfile,
Double_t alpha_,
Double_t mass_,
Double_t n_,
Double_t sigma_
)
{
RooRealVar mbc("mbc", "m_{BC}", 1.83, 1.89, "GeV");
RooRealVar ebeam("ebeam", "Ebeam", 0., 100., "GeV");
RooRealVar chg("chg", "Charge", -2, 2);
RooCategory passed("passed", "Event should be used for plot");
passed.defineType("yes", 1);
passed.defineType("no", 0);
RooRealVar arg_cutoff ("arg_cutoff", "Argus cutoff", 1.8865, 1.885, 1.8875,"GeV");
RooRealVar arg_slope ("arg_slope", "Argus slope", -13, -100, 40);
RooRealVar mbc_float ("mbc_float", "Floating D mass", mass_, "GeV");
RooRealVar sigma ("sigma", "CB width", sigma_, "GeV");
RooRealVar alpha("alpha", "CB shape cutoff", alpha_);
RooRealVar n("n", "CB tail parameter", n_);
RooCBShape cb_float ("cb_float", "Floating Crystal Barrel", mbc, mbc_float, sigma, alpha, n);
RooArgusBG argus("argus", "Argus BG", mbc, arg_cutoff, arg_slope);
RooRealVar yld("yield", "D yield", 0, -30, 100000);
RooRealVar bkg("bkg", "Background", 20, 0, 40000);
// Build pdf
RooAddPdf sumpdf_float("sumpdf_float", "Generic D sum pdf", RooArgList(cb_float, argus),
RooArgList(yld, bkg));
RooDataSet* dset = RooDataSet::read(bkgfile, RooArgList(mbc, ebeam, passed), "", "");
RooPlot* xframe = mbc.frame();
RooDataSet* dset2 = dset->reduce("passed==1");
dset2->plotOn(xframe);
// RooFitResult* rv = sumpdf_float.fitTo(*dset2, Extended(kTRUE), Save(kTRUE),
// Hesse(kTRUE), Verbose(kTRUE));
RooFitResult* rv = sumpdf_float.fitTo(*dset2, "ermh");
sumpdf_float.paramOn(xframe, dset2);
if ((yld.getVal() < 0) && (-yld.getVal()/bkg.getVal() > 0.5)){
yld.setVal(0);
bkg.setVal(1);
}
sumpdf_float.plotOn(xframe);
sumpdf_float.plotOn(xframe, Components(RooArgSet(argus)),
LineColor(kRed), LineStyle(kDashed));
TCanvas* c1 = new TCanvas("c1","Canvas", 2);
xframe->SetTitleOffset(2.2, "Y");
xframe->SetTitleOffset(1.1, "X");
xframe->SetTitle(title);
c1->SetLeftMargin(0.17);
xframe->Draw();
if ( rv && rv->covQual() != 3){
// fit has failed
TText *txt = new TText();
txt->SetTextSize(.08);
txt->SetTextAlign(22);
txt->SetTextAngle(30);
txt->DrawTextNDC(0.5, 0.5, "FAILED");
}
c1->Update();
c1->Print(epsfile);
c1->Clear();
FILE* table = fopen(txtfile.Data(), "w+");
fprintf(table, "Name\t|| Value\t|| Error\n");
// fprintf(table, "yldsigma\t| %.10f\t| \n", yld.getVal()/yld.getError());
fprintf(table, "entries\t| %.10f\t| \n", dset->numEntries());
fprintf(table, "yld\t| %.10f\t| %.10f\n", yld.getVal(), yld.getError());
// fprintf(table, "ratio\t| %.10f\t| \n", yld.getVal()/dset->numEntries());
// fprintf(table, "ratioerr\t| %.10f\t| \n", yld.getError()/dset->numEntries());
fclose(table);
cout << "Saved output as: " << txtfile << endl;
rv->Delete();
}
/*#include "TTree.h"
#include "TCanvas.h"
#include "TGraph.h"
#include "TMultiGraph.h"
#include "TRint.h"
#include <time.h>
#include <stdio.h>
const short MaxNN = 6;
*/
void make_scatterplot(TString data_file, TString plot_name, TString header)
{
//TApplication program = new TRInt();
//pull in data
TTree *t = new TTree();
t->ReadFile(data_file);
t->SetName("t");
//TString name = "run2.png";
TCanvas *BG = new TCanvas("c1", "Read Velocity on Local Disk for CMS3 Files", 1920, 1080);
BG->cd();
TPad *c = new TPad("MainPad", "My main pad", 0, 0, 1, 0.9);
c->Draw();
c->Divide(2,2);
//In first slot have Time vs. Buffer for 1 Concurrent Read
c->cd(1);
int n1 = t->Draw("VelocityMBps:BufferSize", "ConcurrentReads==1", "goff");
TGraph *ghist1 = new TGraph(n1, t->GetV2(), t->GetV1());
ghist1->SetName("ghist1");
ghist1->SetMarkerStyle(3);
ghist1->SetMarkerColor(1);
ghist1->SetTitle("Single File Read");
ghist1->GetXaxis()->SetTitle("Buffer Size (Bytes)");
ghist1->GetYaxis()->SetTitle("Read Velocity (MB/s)");
ghist1->Draw("ap");
//In second slot have Time vs. Buffer for 3 Concurrent Reads
c->cd(2);
int n2 = t->Draw("VelocityMBps:BufferSize", "ConcurrentReads==3", "goff");
TGraph *ghist2 = new TGraph(n2, t->GetV2(), t->GetV1());
Double_t DP2x[n2], DP2y[n2];
ghist2->SetName("ghist2");
ghist2->SetMarkerStyle(3);
ghist2->SetMarkerColor(1);
ghist2->SetTitle("3 Concurrent Reads");
ghist2->GetXaxis()->SetTitle("Buffer Size (Bytes)");
ghist2->GetYaxis()->SetTitle("Read Velocity (MB/s)");
ghist2->Draw("ap");
//In third slot have Time vs. Buffer for 3 Concurrent Reads
c->cd(3);
int n3 = t->Draw("VelocityMBps:BufferSize", "ConcurrentReads==6", "goff");
TGraph *ghist3 = new TGraph(n3, t->GetV2(), t->GetV1());
Double_t DP3x[n3], DP3y[n3];
ghist3->SetName("ghist3");
ghist3->SetMarkerStyle(3);
ghist3->SetMarkerColor(1);
ghist3->SetTitle("6 Concurrent Reads");
ghist3->GetXaxis()->SetTitle("Buffer Size (Bytes)");
ghist3->GetYaxis()->SetTitle("Read Velocity (MB/s)");
ghist3->Draw("ap");
//In fourth slot have Time vs. Buffer for 10 Concurrent Reads
c->cd(4);
int n4 = t->Draw("VelocityMBps:BufferSize", "ConcurrentReads==10", "goff");
TGraph *ghist4 = new TGraph(n4, t->GetV2(), t->GetV1());
ghist4->SetMarkerStyle(3);
ghist4->SetName("ghist4");
ghist4->SetMarkerColor(1);
ghist4->SetTitle("10 Concurrent Reads");
ghist4->GetXaxis()->SetTitle("Buffer Size (Bytes)");
ghist4->GetYaxis()->SetTitle("Read Velocity (MB/s)");
ghist4->Draw("ap");
//Draw to screen
//Initialize Canvas
c->cd(0);
c->Draw();
BG->cd();
TText *title = new TText(.5,.95, header);
title->SetTextAlign(22);
title->Draw();
gDirectory->Add(ghist1);
gDirectory->Add(ghist2);
gDirectory->Add(ghist3);
gDirectory->Add(ghist4);
gDirectory->Add(t);
gPad->SaveAs(plot_name);
}
// -----------------------------------------------------------------------------
//
TCanvas* createPlot( TDirectory* output_file,
std::string output_canvas,
//std::vector< std::vector<std::string> > input_files,
std::vector<std::string> input_files,
std::vector<std::string> input_histos,
std::vector<std::string> input_legend,
std::vector<int> marker_style,
std::vector<int> marker_colour,
std::vector<float> marker_size,
std::vector<float> lumis,
double lumi,
int rebin,
bool norm,
bool log,
double min = -1.,
double max = -1.,
int integral = -1 )
{
// SetSomeStyles();
// Check integral
if ( integral >= input_files.size() ) { integral = -1; }
// Loop through histogram names
std::vector<TH1D*> his;
for ( uint ihis = 0; ihis < input_histos.size(); ++ihis ) {
if ( integral < 0 || ihis <= integral ) {
his.push_back( (TH1D*)getHisto( input_files[ihis], input_histos[ihis], "QcdBkgdEst", rebin ) );
} else {
his.back()->Add( (TH1D*)getHisto( input_files[ihis], input_histos[ihis], "QcdBkgdEst", rebin ) );
}
}
// Create legend
TLegend* legend = new TLegend( 0.65, 0.52-0.03*his.size(), 0.9, 0.52, NULL, "brNDC" );
legend->SetFillColor(0);
legend->SetLineColor(0);
legend->SetTextAlign(31);
legend->SetTextSize(0.025);
// Create canvas
TCanvas* aCanvas = createCanvas( output_canvas, output_file, log );
TPaveText* stats = new TPaveText( 0.65, 0.54, 0.9, 0.54+0.03*his.size(), "brNDC" );
stats->SetFillColor(0);
stats->SetLineColor(0);
TLatex* prelim = new TLatex( 0.15, 0.96, "#scale[0.8]{CMS preliminary 2011}" );
prelim->SetTextSize(0.03);
prelim->SetNDC();
std::stringstream ssl; ssl << "#scale[0.8]{#int L dt = " << lumi/1000. << " fb^{-1}, #sqrt{s} = 7 TeV}";
TLatex* lumistxt = new TLatex( 0.95, 0.96, ssl.str().c_str() );
lumistxt->SetTextSize(0.03);
lumistxt->SetNDC();
lumistxt->SetTextAlign(31);
// For Ted
bool print_ted = true;
std::stringstream ted;
if ( print_ted ) { ted << "self._htMeans = ("; }
// Loop through histograms
double aMax = 0.;
double aMin = 1.e12;
for ( uint ihis = 0; ihis < his.size(); ++ihis ) {
if ( !his[ihis] ) { continue; }
// Line colour and fill
his[ihis]->Scale(lumis[ihis]/100.);
his[ihis]->SetMarkerStyle(marker_style[ihis]);
his[ihis]->SetMarkerColor(marker_colour[ihis]);
his[ihis]->SetMarkerSize(marker_size[ihis]);
his[ihis]->SetLineColor(marker_colour[ihis]);
his[ihis]->SetLineStyle(0);
his[ihis]->SetLineWidth(0);
// Populate legend
if ( input_legend.size() > ihis ) { legend->AddEntry( his[ihis], input_legend[ihis].c_str(), "ep" ); }
else { legend->AddEntry( his[ihis], input_histos[ihis].c_str(), "ep" ); }
// Populate stats box
std::stringstream ss;
//ss << "Mean=" << int(his[ihis]->GetMean()*100.)/100. << ", RMS=" << int(his[ihis]->GetRMS()*100.)/100.;
ss << "Entries=" << his[ihis]->GetEntries();
TText* text = stats->AddText(ss.str().c_str());
text->SetTextAlign(11);
text->SetTextSize(0.025);
text->SetTextColor(marker_colour[ihis]);
// For Ted
if ( print_ted ) ted << std::setw(9) << std::scientific << std::setprecision(3) << his[ihis]->GetMean() << ", ";
// Calc min/max number of entries
if ( his[ihis]->GetMaximum() > aMax ) { aMax = his[ihis]->GetMaximum(); }
if ( his[ihis]->GetMinimum(1.e-12) < aMin ) { aMin = his[ihis]->GetMinimum(1.e-12); }
}
// For Ted
if ( print_ted ) {
ted << ")" << std::endl;
std::cout << ted.str() << std::endl;
}
if ( !his.empty() ) {
if ( his[0] ) his[0]->GetYaxis()->SetTitleOffset(1.43);
if ( his[0] ) his[0]->GetYaxis()->SetTitleSize(0.06);
if ( his[0] ) his[0]->GetXaxis()->SetTitleSize(0.06);
if ( his[0] ) his[0]->GetXaxis()->SetTitleOffset(0.9);
}
// First histo be drawn
bool first = true;
for ( uint ihis = 0; ihis < his.size(); ++ihis ) {
if ( !his[ihis] ) { continue; }
//his[ihis]->GetYaxis()->SetTitle("a.u.");
if ( log ) {
his[ihis]->SetMaximum( aMax * 10. );
his[ihis]->SetMinimum( aMin * 0.1 );
} else {
his[ihis]->SetMaximum( aMax * 1.1 );
his[ihis]->SetMinimum( aMin * 0.9 );
}
if ( min > 0. ) his[ihis]->SetMinimum( min );
if ( max > 0. ) his[ihis]->SetMaximum( max );
if ( norm ) {
std::string options = "";
if ( first ) { options = "Ehist"; first = false; }
else { options = "psame"; }
if ( his[ihis]->GetEntries() > 0. ) { his[ihis]->DrawNormalized(options.c_str()); }
} else {
std::string options = "";
if ( first ) { options = "Ehist"; first = false; }
else { options = "psame"; }
his[ihis]->Draw(options.c_str());
}
} // Loop through histos
output_file->cd();
legend->Draw("same");
stats->Draw("same");
prelim->Draw("same");
lumistxt->Draw("same");
aCanvas->Modified();
//aCanvas->SaveAs( std::string(output_canvas+".png").c_str() );
aCanvas->SaveAs( std::string(output_canvas+".pdf").c_str() );
aCanvas->SaveAs( std::string(output_canvas+".C").c_str() );
aCanvas->Write();
return aCanvas;
}
void archi() {
TCanvas *c1 = new TCanvas("c1","Dictionary Architecture",20,10,750,930);
c1->SetBorderSize(0);
c1->Range(0,0,20.5,26);
TPaveLabel *title = new TPaveLabel(4,24,16,25.5,c1->GetTitle());
title->SetFillColor(46);
title->Draw();
TPavesText *dll = new TPavesText(0.5,19,4.5,23,5,"tr");
dll->SetFillColor(39);
dll->SetTextSize(0.023);
dll->AddText(" ");
dll->AddText("Dynamically");
dll->AddText("Linked");
dll->AddText("Libraries");
dll->Draw();
TPaveLabel *dlltitle = new TPaveLabel(1.5,22.6,3.5,23.3,"DLLs");
dlltitle->SetFillColor(28);
dlltitle->Draw();
TPavesText *cpp = new TPavesText(5.5,19,9.5,23,5,"tr");
cpp->SetTextSize(0.023);
cpp->AddText(" ");
cpp->AddText("Commented");
cpp->AddText("Header");
cpp->AddText("Files");
cpp->Draw();
TPaveLabel *cpptitle = new TPaveLabel(6.5,22.6,8.5,23.3,"C++");
cpptitle->SetFillColor(28);
cpptitle->Draw();
TPavesText *odl = new TPavesText(10.5,19,14.5,23,5,"tr");
odl->SetTextSize(0.023);
odl->AddText(" ");
odl->AddText("Objects");
odl->AddText("Description");
odl->AddText("Files");
odl->Draw();
TPaveLabel *odltitle = new TPaveLabel(11.5,22.6,13.5,23.3,"ODL");
odltitle->SetFillColor(28);
odltitle->Draw();
TPavesText *idl = new TPavesText(15.5,19,19.5,23,5,"tr");
idl->SetTextSize(0.023);
idl->AddText(" ");
idl->AddText("Interface");
idl->AddText("Definition");
idl->AddText("Language");
idl->Draw();
TPaveLabel *idltitle = new TPaveLabel(16.5,22.6,18.5,23.3,"IDL");
idltitle->SetFillColor(28);
idltitle->Draw();
TWbox *p1 = new TWbox(7.8,10,13.2,17,11,12,1);
p1->Draw();
TText *pro1 = new TText(10.5,15.8,"Process 1");
pro1->SetTextAlign(21);
pro1->SetTextSize(0.03);
pro1->Draw();
TPaveText *p1dict = new TPaveText(8.8,13.8,12.2,15.6);
p1dict->SetTextSize(0.023);
p1dict->AddText("Dictionary");
p1dict->AddText("in memory");
p1dict->Draw();
TPavesText *p1object = new TPavesText(8.6,10.6,12.1,13.0,5,"tr");
p1object->SetTextSize(0.023);
p1object->AddText("Objects");
p1object->AddText("in memory");
p1object->Draw();
TWbox *p2 = new TWbox(15.5,10,20,17,11,12,1);
p2->Draw();
TText *pro2 = new TText(17.75,15.8,"Process 2");
pro2->SetTextAlign(21);
pro2->SetTextSize(0.03);
pro2->Draw();
TPaveText *p2dict = new TPaveText(16,13.8,19.5,15.6);
p2dict->SetTextSize(0.023);
p2dict->AddText("Dictionary");
p2dict->AddText("in memory");
p2dict->Draw();
TPavesText *p2object = new TPavesText(16.25,10.6,19.25,13.0,5,"tr");
p2object->SetTextSize(0.023);
p2object->AddText("Objects");
p2object->AddText("in memory");
p2object->Draw();
TWbox *stub1 = new TWbox(12.9,11.5,13.6,15.5,49,3,1);
stub1->Draw();
TText *tstub1 = new TText(13.25,13.5,"Stub1");
tstub1->SetTextSize(0.025);
tstub1->SetTextAlign(22);
tstub1->SetTextAngle(90);
tstub1->Draw();
TWbox *stub2 = new TWbox(15.1,11.5,15.8,15.5,49,3,1);
stub2->Draw();
TText *tstub2 = new TText(15.45,13.5,"Stub2");
tstub2->SetTextSize(0.025);
tstub2->SetTextAlign(22);
tstub2->SetTextAngle(-90);
tstub2->Draw();
TArrow *ar1 = new TArrow();
ar1->SetLineWidth(6);
ar1->SetLineColor(1);
ar1->SetFillStyle(1001);
ar1->SetFillColor(1);
ar1->DrawArrow(13.5,14,15,14,0.012,"|>");
ar1->DrawArrow(15.1,13,13.51,13,0.012,"|>");
TPaveText *cint = new TPaveText(1.0,15.0,8.0,17.5);
cint->SetFillColor(39);
cint->SetBorderSize(1);
cint->SetTextSize(0.023);
cint->AddText("C++ Interpreter");
cint->AddText("and program builder");
cint->Draw();
TPaveText *command = new TPaveText(2.5,13.4,8.0,14.5);
command->SetTextSize(0.023);
command->SetFillColor(39);
command->SetBorderSize(1);
command->AddText("Command Thread");
command->Draw();
TPavesText *view = new TPavesText(1.0,9.5,7.7,12.6,3,"tr");
view->SetFillColor(39);
view->SetBorderSize(2);
view->SetTextSize(0.023);
view->AddText("Viewer Thread(s)");
view->AddText("Picking");
view->AddText("Context Menus");
view->AddText("Inspector/Browser");
view->Draw();
TPavesText *web = new TPavesText(0.5,5,6,8.5,5,"tr");
web->SetTextSize(0.023);
web->AddText(" ");
web->AddText("generated");
web->AddText("automatically");
web->AddText("from dictionary");
web->AddText("and source files");
web->Draw();
TPaveLabel *webtitle = new TPaveLabel(1.5,8.1,5.0,8.8,"HTML Files");
webtitle->SetFillColor(28);
webtitle->Draw();
TPavesText *printed = new TPavesText(0.5,1.0,6,4,5,"tr");
printed->SetTextSize(0.023);
printed->AddText(" ");
printed->AddText("generated");
printed->AddText("automatically");
printed->AddText("from HTML files");
printed->Draw();
TPaveLabel *printedtitle = new TPaveLabel(1.5,3.6,5.0,4.3,"Printed Docs");
printedtitle->SetFillColor(28);
printedtitle->Draw();
TBox *box1 = new TBox(0.2,9.2,14.25,17.8);
box1->SetFillStyle(0);
box1->SetLineStyle(2);
box1->Draw();
TBox *box2 = new TBox(10.2,18.7,20.2,23.6);
box2->SetFillStyle(0);
box2->SetLineStyle(3);
box2->Draw();
ar1->DrawArrow(2.5,17.5,2.5,18.9,0.012,"|>");
ar1->DrawArrow(5.5,9.2,5.5,8.7,0.012,"|>");
ar1->DrawArrow(5.5,5,5.5,4.2,0.012,"|>");
ar1->DrawArrow(8.5,9.2,8.5,8.2,0.012,"|>");
ar1->DrawArrow(9.5,8.1,9.5,9.0,0.012,"|>");
ar1->DrawArrow(6.5,19,6.5,17.6,0.012,"|>");
ar1->DrawArrow(8.5,19,8.5,17.1,0.012,"|>");
ar1->DrawArrow(11.5,19,11.5,17.1,0.012,"|>");
TPaveLabel *ootitle = new TPaveLabel(10.5,7.8,17,8.8,"Objects Data Base");
ootitle->SetFillColor(28);
ootitle->Draw();
TPad *pio = new TPad("pio","pio",0.37,0.02,0.95,0.31,49);
pio->Range(0,0,12,8);
pio->Draw();
pio->cd();
TPavesText *raw = new TPavesText(0.5,1,2.5,6,7,"tr");
raw->Draw();
TPavesText *dst1 = new TPavesText(4,1,5,3,7,"tr");
dst1->Draw();
TPavesText *dst2 = new TPavesText(6,1,7,3,7,"tr");
dst2->Draw();
TPavesText *dst3 = new TPavesText(4,4,5,6,7,"tr");
dst3->Draw();
TPavesText *dst4 = new TPavesText(6,4,7,6,7,"tr");
dst4->Draw();
Float_t xlow = 8.5;
Float_t ylow = 1;
Float_t dx = 0.5;
Float_t dy = 0.5;
for (Int_t j=1;j<9;j++) {
Float_t y0 = ylow + (j-1)*0.7;
Float_t y1 = y0 + dy;
for (Int_t i=1;i<5;i++) {
Float_t x0 = xlow +(i-1)*0.6;
Float_t x1 = x0 + dx;
TPavesText *anal = new TPavesText(x0,y0,x1,y1,7,"tr");
anal->Draw();
}
}
TText *daq = new TText();
daq->SetTextSize(0.07);
daq->SetTextAlign(22);
daq->DrawText(1.5,7.3,"DAQ");
daq->DrawText(6,7.3,"DST");
daq->DrawText(10.,7.3,"Physics Analysis");
daq->DrawText(1.5,0.7,"Events");
daq->DrawText(1.5,0.3,"Containers");
daq->DrawText(6,0.7,"Tracks/Hits");
daq->DrawText(6,0.3,"Containers");
daq->DrawText(10.,0.7,"Attributes");
daq->DrawText(10.,0.3,"Containers");
c1->cd();
}
void Drawmethodcomp(){
const int nDil = 8;
const double centDil[nDil+1] = {1.0,0.7,0.6,0.5,0.4,0.3,0.2,0.1,0};
double centDilbin[nDil];
for(int i=0;i<nDil;i++)
centDilbin[i]=0.29+i+N-1-nDil;
const double NpartDil[nDil] = {8.75,30.51,53.30,86.23,130.06,187.35,261.49,355.45};
const double NpartDilerr[nDil] = {1.13,3.02,3.95,4.35,4.60,4.44,3.96,2.83};
const double NcollDil[nDil] = {8.01,38.86,86.85,175.76,326.06,563.21,926.79,1484.49};
const double NcollDilerr[nDil] = {1.41,6.41,12.48,21.13,34.27,52.66,81.37,120.0};
TGraphErrors *graphNpartDil = new TGraphErrors(nDil,centDilbin,NpartDil,0,NpartDilerr);
TGraphErrors *graphNcollDil = new TGraphErrors(nDil,centDilbin,NcollDil,0,NcollDilerr);
TCanvas *c1 = new TCanvas("c1","c1",1,1,550,460);
c1->SetLogy();
c1->SetFillColor(10);
c1->SetFrameFillColor(0);
c1->SetFrameBorderSize(0);
c1->SetFrameBorderMode(0);
c1->SetLeftMargin(0.15);
c1->SetBottomMargin(0.15);
c1->SetTopMargin(0.02);
c1->SetRightMargin(0.02);
gStyle->SetOptStat(0);
c1->SetTicks(-1);
N=N-1;
TString str="Ncoll";
TH1D* hist = new TH1D("","",N,0,N);
hist->GetXaxis()->SetNdivisions(502);
if(method==0)
hist->SetXTitle("Centrality");
else
hist->SetXTitle("HF #Sigma E_{T} |#eta|>3");
hist->SetYTitle(Form("<%s> and systematic errors",str.Data()));
hist->SetMinimum(1);
hist->SetMaximum(3999.99);
hist->GetXaxis()->CenterTitle(0);
hist->GetYaxis()->CenterTitle(1);
hist->GetYaxis()->SetTitleOffset(1.1);
hist->GetXaxis()->SetTitleOffset(1.1);
hist->GetXaxis()->SetTitleSize(0.056);
hist->GetYaxis()->SetTitleSize(0.056);
hist->GetXaxis()->SetLabelSize(0.05);
hist->GetYaxis()->SetLabelSize(0.05);
hist->GetXaxis()->SetLabelOffset(99);
hist->Draw();
TFile *f=TFile::Open(outG);
TGraphErrors* graph = (TGraphErrors*)f->Get(Form("std/%s_graph",str.Data()));
// TGraphErrors* Gri055_graph = (TGraphErrors*)f->Get(Form("Gri055/%s_graph",str.Data()));
// TGraphErrors* Gri101_graph = (TGraphErrors*)f->Get(Form("Gri101/%s_graph",str.Data()));
TVectorD *centbin = (TVectorD*)f->Get(Form("std/G0/centbin"));
TVectorD *kpoint = (TVectorD*)f->Get(Form("std/G0/kpoint"));
TFile *f1=TFile::Open(Form("../../OneComp/double_side/%s",outG.Data()));
TGraphErrors* graph1 = (TGraphErrors*)f1->Get(Form("std/%s_graph",str.Data()));
// TGraphErrors* Gri055_graph = (TGraphErrors*)f->Get(Form("Gri055/%s_graph",str.Data()));
// TGraphErrors* Gri101_graph = (TGraphErrors*)f->Get(Form("Gri101/%s_graph",str.Data()));
TVectorD *centbin1 = (TVectorD*)f1->Get(Form("std/G0/centbin"));
TVectorD *kpoint1 = (TVectorD*)f1->Get(Form("std/G0/kpoint"));
graph->SetTitle("g1");
graph->SetMarkerStyle(20);
graph->SetMarkerColor(1);
graph->SetLineColor(1);
graph->SetLineWidth(2);
graph->SetMarkerSize(1.2);
graph->Draw("Psameez");
graph1->SetMarkerStyle(24);
graph1->SetMarkerColor(4);
graph1->SetLineColor(4);
graph1->SetLineWidth(2);
graph1->SetMarkerSize(1.2);
graph1->Draw("Psameez");
/*
Gri055_graph->SetTitle("g2");
Gri055_graph->SetMarkerStyle(33);
Gri055_graph->SetMarkerColor(2);
Gri055_graph->SetLineColor(2);
Gri055_graph->SetLineWidth(2);
Gri055_graph->SetMarkerSize(1.2);
Gri055_graph->Draw("Psameez");
Gri101_graph->SetTitle("g3");
Gri101_graph->SetMarkerStyle(34);
Gri101_graph->SetMarkerColor(4);
Gri101_graph->SetLineColor(4);
Gri101_graph->SetLineWidth(2);
Gri101_graph->SetMarkerSize(1.2);
Gri101_graph->Draw("Psameez");
*/
graphNpartDil->SetMarkerSize(1.2);
graphNpartDil->SetLineColor(2);
graphNpartDil->SetLineWidth(2);
graphNpartDil->SetMarkerStyle(33);
graphNpartDil->SetMarkerColor(2);
graphNcollDil->SetMarkerSize(1.2);
graphNcollDil->SetLineColor(2);
graphNcollDil->SetLineWidth(2);
graphNcollDil->SetMarkerStyle(33);
graphNcollDil->SetMarkerColor(2);
if(str=="Npart") graphNpartDil->Draw("Psame");
if(str=="Ncoll") graphNcollDil->Draw("Psame");
std::vector<TString> label(N);
for(int i=0;i<N;i++)
if(method==0)label[i] = Form("%.2f-%.2f%%",(*centbin)[i]*100,(*centbin)[i+1]*100);
else label[i] = Form("%.2f-%.2f",(*kpoint)[i],(*kpoint)[i+1]);
TLatex *tex1= new TLatex(0.2,0.9,"CMS Preliminary PbPb #sqrt{s_{NN}} = 2.76 TeV");
tex1->SetNDC();
tex1->SetTextColor(1);
tex1->SetTextFont(42);
tex1->SetTextSize(0.05);
tex1->Draw();
double y = gPad->GetUymin();
// - 0.2*h->GetYaxis()->GetBinWidth(1);
TText t;
t.SetTextAngle(45);
t.SetTextSize(0.03);
t.SetTextAlign(33);
for (int i=0;i<N;i++) {
double x = hist->GetXaxis()->GetBinCenter(i+1);
t.DrawText(x,y,label[i]);
}
TLegend *leg0 = new TLegend(0.18,0.70,0.50,0.85);
leg0->SetFillColor(10);
leg0->SetBorderSize(0);
leg0->SetTextFont(42);
leg0->SetTextSize(0.047);
leg0->AddEntry(graph,"From Ancestor fitting","p");
leg0->AddEntry(graph1,"From One Comp fitting","p");
// leg0->AddEntry(Gri055_graph,"Gribov #Omega=0.55","p");
// leg0->AddEntry(Gri101_graph,"Gribov #Omega=1.01","p");
if(str=="Npart") leg0->AddEntry(graphNpartDil,"Npart from run I","p");
if(str=="Ncoll") leg0->AddEntry(graphNcollDil,"Ncoll from run I","p");
leg0->Draw();
c1->SaveAs(Form("%sGri.png",str.Data()));
c1->SaveAs(Form("%sGri.pdf",str.Data()));
TCanvas *c2 = new TCanvas("c2","c2",1,1,550,460);
c2->SetFillColor(10);
c2->SetFrameFillColor(0);
c2->SetFrameBorderSize(0);
c2->SetFrameBorderMode(0);
c2->SetLeftMargin(0.15);
c2->SetBottomMargin(0.15);
c2->SetTopMargin(0.02);
c2->SetRightMargin(0.02);
c2->SetTicks(-1);
c2->cd();
TGraphErrors* fdvf1 = (TGraphErrors*)graph->Clone("fdvf1");
TGraphErrors* f1dvDil = (TGraphErrors*)graph->Clone("f1dvDil");
TGraphErrors* fdvDil = (TGraphErrors*)graph->Clone("fdvDil");
TGraphErrors* Dilsys = (TGraphErrors*)graph->Clone("Dilsys");
TGraphErrors* Dilsys2 = (TGraphErrors*)graph->Clone("Dilsys2");
for(int ip = 0;ip<fdvf1->GetN();ip++){
double x = graph->GetX()[ip];
double ey = graph->GetEY()[ip];
double y = graph->GetY()[ip];
double ey1 = graph1->GetEY()[ip];
double y1 = graph1->GetY()[ip];
double eyNcoll = graphNcollDil->GetEY()[ip];
double yNcoll = graphNcollDil->GetY()[ip];
double eyNpart = graphNpartDil->GetEY()[ip];
double yNpart = graphNpartDil->GetY()[ip];
if(y1!=0 && y!=0){
fdvf1->SetPoint(ip,x,y/y1);
fdvf1->SetPointError(ip,0,y/y1*sqrt((ey/y)**2+(ey1/y1)**2));
if(str=="Ncoll"){
fdvDil->SetPoint(ip,x,y/yNcoll);
double yerr = y/yNcoll*sqrt((ey/y)**2+(eyNcoll/yNcoll)**2);
fdvDil->SetPointError(ip,0,0);
f1dvDil->SetPoint(ip,x,y1/yNcoll);
double y1err = y1/yNcoll*sqrt((ey1/y1)**2+(eyNcoll/yNcoll)**2);
f1dvDil->SetPointError(ip,0,0);
Dilsys->SetPoint(ip,x,1.);
Dilsys->SetPointError(ip,0,TMath::Max(fabs(y1/yNcoll-1),fabs(y/yNcoll-1)));
Dilsys2->SetPoint(ip,x,1.);
Dilsys2->SetPointError(ip,0,TMath::Max(yerr,y1err));
}
else if(str=="Npart"){
fdvDil->SetPoint(ip,x,y/yNpart);
double yerr = y/yNpart*sqrt((ey/y)**2+(eyNcoll/yNpart)**2);
fdvDil->SetPointError(ip,0,0);
f1dvDil->SetPoint(ip,x,y1/yNpart);
double y1err = y1/yNpart*sqrt((ey1/y1)**2+(eyNcoll/yNcoll)**2);
f1dvDil->SetPointError(ip,0,0);
Dilsys->SetPoint(ip,x,1.);
Dilsys->SetPointError(ip,0,TMath::Max(fabs(y1/yNpart-1),fabs(y/yNpart-1)));
Dilsys2->SetPoint(ip,x,1.);
Dilsys2->SetPointError(ip,0,TMath::Max(yerr,y1err));
}
}
}
hist->SetMinimum(0.5);
hist->SetMaximum(1.5);
hist->SetYTitle(Form("<%s> ratio",str.Data()));
hist->Draw();
fdvf1->SetMarkerStyle(20);
fdvf1->SetMarkerColor(1);
fdvf1->SetLineColor(1);
fdvf1->SetLineWidth(2);
fdvf1->SetMarkerSize(1.4);
//fdvf1->Draw("Psameez");
f1dvDil->SetMarkerStyle(24);
f1dvDil->SetMarkerColor(2);
f1dvDil->SetLineColor(2);
f1dvDil->SetLineWidth(2);
f1dvDil->SetMarkerSize(1.4);
fdvDil->SetMarkerStyle(20);
fdvDil->SetMarkerColor(4);
fdvDil->SetLineColor(4);
fdvDil->SetLineWidth(2);
fdvDil->SetMarkerSize(1.4);
Dilsys->SetFillColor(kYellow+2);
Dilsys2->SetFillColor(kYellow+1);
//Dilsys2->Draw("e3same");
Dilsys->Draw("e3same");
fdvDil->Draw("Psameez");
f1dvDil->Draw("Psameez");
tex1->Draw();
TLine *l = new TLine(0,1,N,1);
l->SetLineStyle(2);
l->SetLineWidth(3);
l->Draw("same");
TLegend *leg1 = new TLegend(0.38,0.72,0.75,0.86);
leg1->SetFillColor(10);
leg1->SetBorderSize(0);
leg1->SetTextFont(42);
leg1->SetTextSize(0.04);
//leg1->AddEntry(fdvf1,"Ancestor fitting/One Comp fitting","p");
leg1->AddEntry(f1dvDil,"One Component fitting/Run I fitting","p");
leg1->AddEntry(fdvDil,"Two Component fitting/Run I fitting","p");
leg1->AddEntry(Dilsys,"RunI fitting systematics","f");
// leg1->AddEntry(Dilsys2,"RunI fitting ratio systematics","f");
double y = gPad->GetUymin()+0.5;
for (int i=0;i<N;i++) {
double x = hist->GetXaxis()->GetBinCenter(i+1);
t.DrawText(x,y,label[i]);
}
leg1->Draw();
c2->SaveAs(Form("%sGriratio.png",str.Data()));
c2->SaveAs(Form("%sGriratio.pdf",str.Data()));
}
void tree() {
//
// This macro displays the Tree data structures
//Author: Rene Brun
TCanvas* c1 = new TCanvas("c1","Tree Data Structure",200,10,750,940);
c1->Range(0,-0.1,1,1.15);
gBenchmark->Start("tree");
Int_t branchcolor = 26;
Int_t leafcolor = 30;
Int_t basketcolor = 42;
Int_t offsetcolor = 43;
TPaveLabel* title = new TPaveLabel(.3,1.05,.8,1.13,c1->GetTitle());
title->SetFillColor(16);
title->Draw();
TPaveText* treePave = new TPaveText(.01,.75,.15,1.00);
treePave->SetFillColor(18);
treePave->SetTextAlign(12);
TText *tnt = treePave->AddText("Tree");
tnt->SetTextAlign(22);
tnt->SetTextSize(0.030);
treePave->AddText("fScanField");
treePave->AddText("fMaxEventLoop");
treePave->AddText("fMaxVirtualSize");
treePave->AddText("fEntries");
treePave->AddText("fDimension");
treePave->AddText("fSelectedRows");
treePave->Draw();
TPavesText* farm = new TPavesText(.01,1.02,.15,1.1,9,"tr");
TText *tfarm = farm->AddText("CHAIN");
tfarm->SetTextSize(0.024);
farm->AddText("Collection");
farm->AddText("of Trees");
farm->Draw();
TLine* llink = new TLine(.15,.92,.80,.92);
llink->SetLineWidth(2);
llink->SetLineColor(1);
llink->Draw();
llink->DrawLine(.21,.87,.21,.275);
llink->DrawLine(.23,.87,.23,.375);
llink->DrawLine(.25,.87,.25,.775);
llink->DrawLine(.41,.25,.41,-.025);
llink->DrawLine(.43,.25,.43,.075);
llink->DrawLine(.45,.25,.45,.175);
TPaveLabel* branch0 = new TPaveLabel(.20,.87,.35,.97,"Branch 0");
branch0->SetTextSize(0.35);
branch0->SetFillColor(branchcolor);
branch0->Draw();
TPaveLabel* branch1 = new TPaveLabel(.40,.87,.55,.97,"Branch 1");
branch1->SetTextSize(0.35);
branch1->SetFillColor(branchcolor);
branch1->Draw();
TPaveLabel* branch2 = new TPaveLabel(.60,.87,.75,.97,"Branch 2");
branch2->SetTextSize(0.35);
branch2->SetFillColor(branchcolor);
branch2->Draw();
TPaveLabel* branch3 = new TPaveLabel(.80,.87,.95,.97,"Branch 3");
branch3->SetTextSize(0.35);
branch3->SetFillColor(branchcolor);
branch3->Draw();
TPaveLabel* leaf0 = new TPaveLabel(.4,.75,.5,.8,"Leaf 0");
leaf0->SetFillColor(leafcolor);
leaf0->Draw();
TPaveLabel* leaf1 = new TPaveLabel(.6,.75,.7,.8,"Leaf 1");
leaf1->SetFillColor(leafcolor);
leaf1->Draw();
TPaveLabel* leaf2 = new TPaveLabel(.8,.75,.9,.8,"Leaf 2");
leaf2->SetFillColor(leafcolor);
leaf2->Draw();
TPaveText* firstevent = new TPaveText(.4,.35,.9,.4);
firstevent->AddText("First event of each basket");
firstevent->AddText("Array of fMaxBaskets Integers");
firstevent->SetFillColor(basketcolor);
firstevent->Draw();
TPaveLabel* basket0 = new TPaveLabel(.4,.25,.5,.3,"Basket 0");
basket0->SetFillColor(basketcolor);
basket0->Draw();
TPaveLabel* basket1 = new TPaveLabel(.6,.25,.7,.3,"Basket 1");
basket1->SetFillColor(basketcolor);
basket1->Draw();
TPaveLabel* basket2 = new TPaveLabel(.8,.25,.9,.3,"Basket 2");
basket2->SetFillColor(basketcolor);
basket2->Draw();
TPaveText* offset = new TPaveText(.55,.15,.9,.2);
offset->AddText("Offset of events in fBuffer");
offset->AddText("Array of fEventOffsetLen Integers");
offset->AddText("(if variable length structure)");
offset->SetFillColor(offsetcolor);
offset->Draw();
TPaveText* buffer = new TPaveText(.55,.05,.9,.1);
buffer->AddText("Basket buffer");
buffer->AddText("Array of fBasketSize chars");
buffer->SetFillColor(offsetcolor);
buffer->Draw();
TPaveText* zipbuffer = new TPaveText(.55,-.05,.75,.0);
zipbuffer->AddText("Basket compressed buffer");
zipbuffer->AddText("(if compression)");
zipbuffer->SetFillColor(offsetcolor);
zipbuffer->Draw();
TArrow* ar1 = new TArrow();
ar1->SetLineWidth(2);
ar1->SetLineColor(1);
ar1->SetFillStyle(1001);
ar1->SetFillColor(1);
ar1->DrawArrow(.21,.275,.39,.275,0.015,"|>");
ar1->DrawArrow(.23,.375,.39,.375,0.015,"|>");
ar1->DrawArrow(.25,.775,.39,.775,0.015,"|>");
ar1->DrawArrow(.50,.775,.59,.775,0.015,"|>");
ar1->DrawArrow(.70,.775,.79,.775,0.015,"|>");
ar1->DrawArrow(.50,.275,.59,.275,0.015,"|>");
ar1->DrawArrow(.70,.275,.79,.275,0.015,"|>");
ar1->DrawArrow(.45,.175,.54,.175,0.015,"|>");
ar1->DrawArrow(.43,.075,.54,.075,0.015,"|>");
ar1->DrawArrow(.41,-.025,.54,-.025,0.015,"|>");
TLine* ldot = new TLine(.95,.92,.99,.92);
ldot->SetLineStyle(3);
ldot->Draw();
ldot->DrawLine(.9,.775,.99,.775);
ldot->DrawLine(.9,.275,.99,.275);
ldot->DrawLine(.55,.05,.55,0);
ldot->DrawLine(.9,.05,.75,0);
TText* pname = new TText(.46,.21,"fEventOffset");
pname->SetTextFont(72);
pname->SetTextSize(0.018);
pname->Draw();
pname->DrawText(.44,.11,"fBuffer");
pname->DrawText(.42,.01,"fZipBuffer");
pname->DrawText(.26,.81,"fLeaves = TObjArray of TLeaf");
pname->DrawText(.24,.40,"fBasketEvent");
pname->DrawText(.22,.31,"fBaskets = TObjArray of TBasket");
pname->DrawText(.20,1.0,"fBranches = TObjArray of TBranch");
TPaveText* ntleaf = new TPaveText(0.30,.42,.62,.7);
ntleaf->SetTextSize(0.014);
ntleaf->SetFillColor(leafcolor);
ntleaf->SetTextAlign(12);
ntleaf->AddText("fLen: number of fixed elements");
ntleaf->AddText("fLenType: number of bytes of data type");
ntleaf->AddText("fOffset: relative to Leaf0-fAddress");
ntleaf->AddText("fNbytesIO: number of bytes used for I/O");
ntleaf->AddText("fIsPointer: True if pointer");
ntleaf->AddText("fIsRange: True if leaf has a range");
ntleaf->AddText("fIsUnsigned: True if unsigned");
ntleaf->AddText("*fLeafCount: points to Leaf counter");
ntleaf->AddText(" ");
ntleaf->AddLine(0,0,0,0);
ntleaf->AddText("fName = Leaf name");
ntleaf->AddText("fTitle = Leaf type (see Type codes)");
ntleaf->Draw();
TPaveText* type = new TPaveText(.65,.42,.95,.7);
type->SetTextAlign(12);
type->SetFillColor(leafcolor);
type->AddText(" ");
type->AddText("C : a character string");
type->AddText("B : an 8 bit signed integer");
type->AddText("b : an 8 bit unsigned integer");
type->AddText("S : a 16 bit signed short integer");
type->AddText("s : a 16 bit unsigned short integer");
type->AddText("I : a 32 bit signed integer");
type->AddText("i : a 32 bit unsigned integer");
type->AddText("F : a 32 bit floating point");
type->AddText("D : a 64 bit floating point");
type->AddText("TXXXX : a class name TXXXX");
type->Draw();
TPaveLabel* typecode = new TPaveLabel(.7,.68,.9,.72,"fType codes");
typecode->SetFillColor(leafcolor);
typecode->Draw();
ldot->DrawLine(.4,.75,.30,.7);
ldot->DrawLine(.5,.75,.62,.7);
TPaveText* ntbasket = new TPaveText(0.02,-0.07,0.35,.25);
ntbasket->SetFillColor(basketcolor);
ntbasket->SetTextSize(0.014);
ntbasket->SetTextAlign(12);
ntbasket->AddText("fNbytes: Size of compressed Basket");
ntbasket->AddText("fObjLen: Size of uncompressed Basket");
ntbasket->AddText("fDatime: Date/Time when written to store");
ntbasket->AddText("fKeylen: Number of bytes for the key");
ntbasket->AddText("fCycle : Cycle number");
ntbasket->AddText("fSeekKey: Pointer to Basket on file");
ntbasket->AddText("fSeekPdir: Pointer to directory on file");
ntbasket->AddText("fClassName: 'TBasket'");
ntbasket->AddText("fName: Branch name");
ntbasket->AddText("fTitle: TreePave name");
ntbasket->AddText(" ");
ntbasket->AddLine(0,0,0,0);
ntbasket->AddText("fNevBuf: Number of events in Basket");
ntbasket->AddText("fLast: pointer to last used byte in Basket");
ntbasket->Draw();
ldot->DrawLine(.4,.3,0.02,0.25);
ldot->DrawLine(.5,.25,0.35,-.07);
ldot->DrawLine(.5,.3,0.35,0.25);
TPaveText* ntbranch = new TPaveText(0.02,0.40,0.18,0.68);
ntbranch->SetFillColor(branchcolor);
ntbranch->SetTextSize(0.015);
ntbranch->SetTextAlign(12);
ntbranch->AddText("fBasketSize");
ntbranch->AddText("fEventOffsetLen");
ntbranch->AddText("fMaxBaskets");
ntbranch->AddText("fEntries");
ntbranch->AddText("fAddress of Leaf0");
ntbranch->AddText(" ");
ntbranch->AddLine(0,0,0,0);
ntbranch->AddText("fName: Branchname");
ntbranch->AddText("fTitle: leaflist");
ntbranch->Draw();
ldot->DrawLine(.2,.97,.02,.68);
ldot->DrawLine(.35,.97,.18,.68);
ldot->DrawLine(.35,.87,.18,.40);
TPavesText* basketstore = new TPavesText(.8,-0.088,0.952,-0.0035,7,"tr");
basketstore->SetFillColor(28);
basketstore->AddText("Baskets");
basketstore->AddText("Stores");
basketstore->Draw();
c1->Update();
gBenchmark->Show("treePave");
}
void EMCDistribution_PeakSample_Fast(bool full_gain = false)
{
const TString gain = "RAW";
TString hname = "EMCDistribution_" + gain + TString(full_gain ? "_FullGain" : "") + cuts;
TH2 *h2 = NULL;
{
if (full_gain)
{
h2 = new TH2F(hname,
Form(";Calibrated Tower Energy Sum (ADC);Count / bin"), 100,
.05 * 100, 25 * 100, 64, -.5, 63.5);
QAHistManagerDef::useLogBins(h2->GetXaxis());
}
else
{
h2 = new TH2F(hname,
Form(";Calibrated Tower Energy Sum (ADC);Count / bin"), 260,
-.2 * 100, 5 * 100, 64, -.5, 63.5);
}
T->Draw(
"TOWER_" + gain + "_CEMC[].get_bineta() + 8* TOWER_" + gain + "_CEMC[].get_binphi():(TOWER_RAW_CEMC[].signal_samples[10] - TOWER_RAW_CEMC[].signal_samples[0])*(-1)>>" + hname, "", "goff");
}
TText *t;
TCanvas *c1 = new TCanvas(
"EMCDistribution_PeakSample_Fast_" + TString(full_gain ? "_FullGain" : "") + cuts,
"EMCDistribution_PeakSample_Fast_" + TString(full_gain ? "_FullGain" : "") + cuts, 1800, 950);
c1->Divide(8, 8, 0., 0.01);
int idx = 1;
TPad *p;
for (int iphi = 8 - 1; iphi >= 0; iphi--)
{
for (int ieta = 0; ieta < 8; ieta++)
{
p = (TPad *) c1->cd(idx++);
c1->Update();
p->SetLogy();
if (full_gain)
{
p->SetLogx();
}
p->SetGridx(0);
p->SetGridy(0);
TString hname = Form("hEnergy_ieta%d_iphi%d", ieta, iphi) + TString(full_gain ? "_FullGain" : "");
TH1 *h = h2->ProjectionX(hname, ieta + 8 * iphi + 1,
ieta + 8 * iphi + 1); // axis bin number is encoded as ieta+8*iphi+1
h->SetLineWidth(0);
h->SetLineColor(kBlue + 3);
h->SetFillColor(kBlue + 3);
h->GetXaxis()->SetTitleSize(.09);
h->GetXaxis()->SetLabelSize(.08);
h->GetYaxis()->SetLabelSize(.08);
h->Draw();
if (full_gain)
h->Fit("x*gaus", "M");
else
h->Fit("landau", "M");
double peak = -1;
TF1 *fit = ((TF1 *) (h->GetListOfFunctions()->At(0)));
if (fit)
{
fit->SetLineColor(kRed);
peak = fit->GetParameter(1);
}
cout << Form("Finished <Col%d Row%d> = %.1f", ieta, iphi, peak)
<< endl;
TText *t = new TText(.9, .9,
Form("<Col%d Row%d> = %.1f", ieta, iphi, peak));
t->SetTextAlign(33);
t->SetTextSize(.15);
t->SetNDC();
t->Draw();
}
}
SaveCanvas(c1,
TString(_file0->GetName()) + TString("_DrawPrototype3EMCalTower_") + TString(c1->GetName()), false);
}
