void spectrum(char *filename = "example.root")
{
TFile *infile = new TFile(filename);
TTree *data = (TTree*) infile->FindObjectAny("data");
TH1F *h[10];
TCanvas *c = new TCanvas("c","",960,480);
TPad *p = new TPad("p","",0,0,1,1);
p->Divide(4,2,0.01,0.01);
p->Draw();
double norm;
for (int i=0;i<8;i++)
{
p->cd(i+1);
p->GetPad(i+1)->SetLogy();
h[i] = new TH1F(Form("h%d",i),Form("Ch %d",i),60,-0.1,1.1);
// data->Draw(Form("-min[%d]>>h%d",i,i),"min[0]<min[1]&&min[0]<min[2]&&min[0]<min[3]&&min[0]<min[4]&&min[0]<min[5]&&min[0]<min[6]&&min[0]<min[7]");
data->Draw(Form("-min[%d]/256./0.86>>h%d",i,i));
if (i==0) norm = h[i]->GetEntries();
h[i]->SetTitle(Form("Ch %d",i));
h[i]->SetXTitle("Pulse height (Volt)");
h[i]->SetYTitle("#");
h[i]->SetMarkerSize(0.5);
h[i]->SetMarkerColor(2);
h[i]->Sumw2();
h[i]->Scale(1./norm);
h[i]->SetAxisRange(1e-7,10,"Y");
h[i]->Draw("e");
}
}
void RDK2AnalysisPlotter::makeEPGPlot(CoDet detType)
{
gROOT->cd();
TCanvas* theCanvas;
TPad* mainPad;
TPad* titlePad;
TPad* gEPad;
TPad* detPad;
TPad* gEPadSubs[3];
int numTitleLines=getTitleBoxLines( detType);
int canvasNumPixelsYPlot=600;
int canvasNumPixelsYPlots=2*canvasNumPixelsYPlot;
int canvasNumPixelsYTitle=40*numTitleLines;
int canvasNumPixelsY=canvasNumPixelsYPlots+canvasNumPixelsYTitle;
theCanvas = new TCanvas("EPGExpMCAnalysisComparisonPlot","EPGExpMCAnalysisComparisonPlot",10,10,1200,canvasNumPixelsY);
mainPad=(TPad*) theCanvas->GetPad(0);
double ylow,yhigh;
yhigh=1;
ylow=1.-canvasNumPixelsYTitle/(double)canvasNumPixelsY;
titlePad=new TPad("titlePad", "titlePad", 0., ylow, 1., yhigh, -1, 1, 1);
yhigh=ylow;
ylow=ylow-canvasNumPixelsYPlot/(double)canvasNumPixelsY;
gEPad=new TPad("gEPad", "gEPad", 0., ylow, 1., yhigh, -1, 1, 1);
yhigh=ylow;
ylow=0;
detPad=new TPad("detPad", "detPad", 0., ylow, 1., yhigh, -1, 1, 1);
titlePad->Draw();
gEPad->Draw();
detPad->Draw();
// titlePad=(TPad*) theCanvas->GetPad(1);
// gEPad=(TPad*) theCanvas->GetPad(2);
// detPad=(TPad*) theCanvas->GetPad(3);
// egTPad=(TPad*) theCanvas->GetPad(4);
titlePad->SetMargin(0.05,0.0,0.0,0.0);
gEPad->SetMargin(0.05,0.0,0.0,0.05);
if(numExp>1 && numMC>1)
{
gEPad->Divide(2);
gEPadSubs[0]=(TPad*) gEPad->GetPad(1);
TPad* tempPad;
tempPad=(TPad*) gEPad->GetPad(2);
tempPad->Divide(1,2);
gEPadSubs[1]=(TPad*) tempPad->GetPad(1);
gEPadSubs[2]=(TPad*) tempPad->GetPad(2);
TPad* allBasePads[5]={titlePad,gEPadSubs[0],gEPadSubs[1],gEPadSubs[2],detPad};
for (int i = 1;i< 5;i++)
{
allBasePads[i]->SetGrid(1,1);
allBasePads[i]->SetTickx(1);
allBasePads[i]->SetTicky(1);
}
gEPadSubs[2]->SetLogx();
gEPadSubs[0]->SetMargin(0.15,0.1,.1,0.1);
gEPadSubs[1]->SetMargin(0.1,0.1,.1,0.1);
gEPadSubs[2]->SetMargin(0.1,0.1,.1,0.1);
//Plot gE
gEPadSubs[0]->cd(); drawPlot(detType, PLOTVAR_GE,PLOT_COMP);
gEPadSubs[1]->cd(); drawPlot(detType, PLOTVAR_GE,PLOT_RESID);
gEPadSubs[2]->cd(); drawPlot(detType, PLOTVAR_GEVAR,PLOT_NORMRESID);
}
else
{
TPad* allBasePads[3]={titlePad,gEPad,detPad};
for (int i = 1;i< 3;i++)
{
allBasePads[i]->SetGrid(1,1);
allBasePads[i]->SetTickx(1);
allBasePads[i]->SetTicky(1);
}
gEPad->SetMargin(0.13,0.1,.1,0.1);
//Plot gE
gEPad->cd(); drawPlot(detType, PLOTVAR_GE,PLOT_COMP);
}
mainPad->SetFillColor(kGray);
///Make Title box
titlePad->cd();
titlePad->SetFillColor(kGray);
TPaveText* titleBox = makeTitleBox( detType);
titleBox->Draw();
detPad->SetMargin(0.1,0.05,.1,0.1);
//Plot det
detPad->cd(); drawPlot(detType, PLOTVAR_GE,PLOT_DETS);
TString coTypeString;
if(detType==DET_EP)
{
coTypeString="EP";
}
else if(detType==DET_EPG)
{
coTypeString="EPG";
}
else if(detType==DET_EPBG)
{
coTypeString="EPBG";
}
TString imagePath=GRAPHS_DIR;
imagePath+="layouts/PlotLayout";
if(numMC) imagePath+="_MC";
for (int i = 0;i< numMC;i++)
{
imagePath+=TString("_")+mc[i]->GetName();
}
if(numExp) imagePath+="_Exp";
for (int i = 0;i< numExp;i++)
{
imagePath+=TString("_")+exp[i]->GetName();
}
imagePath+="_"+coTypeString+"_EPGPlots.pdf";
theCanvas->SaveAs(imagePath);
delete theCanvas;
clearPlotHists();
}
void RDK2AnalysisPlotter::makeEPPlot(CoDet detType)
{
gROOT->cd();
TCanvas* theCanvas;
TPad* mainPad;
TPad* titlePad;
TPad* pTPad;
TPad* eEPad;
TPad* pEPad;
TPad* pTPadSubs[3];
TPad* eEPadSubs[3];
TPad* pEPadSubs[3];
int numTitleLines=getTitleBoxLines( detType);
int canvasNumPixelsYPlot=600;
int canvasNumPixelsYPlots=3*canvasNumPixelsYPlot;
int canvasNumPixelsYTitle=40*numTitleLines;
int canvasNumPixelsY=canvasNumPixelsYPlots+canvasNumPixelsYTitle;
theCanvas = new TCanvas("EPExpMCAnalysisComparisonPlot","EPExpMCAnalysisComparisonPlot",10,10,1200,canvasNumPixelsY);
mainPad=(TPad*) theCanvas->GetPad(0);
double ylow,yhigh;
yhigh=1;
ylow=1.-canvasNumPixelsYTitle/(double)canvasNumPixelsY;
titlePad=new TPad("titlePad", "titlePad", 0., ylow, 1., yhigh, -1, 1, 1);
yhigh=1.-canvasNumPixelsYTitle/(double)canvasNumPixelsY;
ylow=1.-canvasNumPixelsYTitle/(double)canvasNumPixelsY-canvasNumPixelsYPlot/(double)canvasNumPixelsY;
pTPad=new TPad("pTPad", "pTPad", 0., ylow, 1., yhigh, -1, 1, 1);
yhigh=1.-canvasNumPixelsYTitle/(double)canvasNumPixelsY-canvasNumPixelsYPlot/(double)canvasNumPixelsY;
ylow=1.-canvasNumPixelsYTitle/(double)canvasNumPixelsY- 2*canvasNumPixelsYPlot/(double)canvasNumPixelsY;
eEPad=new TPad("eEPad", "eEPad", 0., ylow, 1., yhigh, -1, 1, 1);
yhigh=1.-canvasNumPixelsYTitle/(double)canvasNumPixelsY-2*canvasNumPixelsYPlot/(double)canvasNumPixelsY;
ylow=0;
pEPad=new TPad("pEPad", "pEPad", 0., ylow, 1., yhigh, -1, 1, 1);
titlePad->Draw();
pTPad->Draw();
eEPad->Draw();
pEPad->Draw();
// titlePad=(TPad*) theCanvas->GetPad(1);
// pTPad=(TPad*) theCanvas->GetPad(2);
// eEPad=(TPad*) theCanvas->GetPad(3);
// pEPad=(TPad*) theCanvas->GetPad(4);
titlePad->SetMargin(0.05,0.0,0.0,0.0);
pTPad->SetMargin(0.05,0.0,0.0,0.05);
pEPad->SetMargin(0.05,0.0,0.0,0.05);
eEPad->SetMargin(0.05,0.0,0.0,0.05);
if(numExp+numMC >1)
{
pTPad->Divide(2);
eEPad->Divide(2);
pEPad->Divide(2);
pTPadSubs[0]=(TPad*) pTPad->GetPad(1);
eEPadSubs[0]=(TPad*) eEPad->GetPad(1);
pEPadSubs[0]=(TPad*) pEPad->GetPad(1);
TPad* tempPad;
tempPad=(TPad*) pTPad->GetPad(2);
tempPad->Divide(1,2);
pTPadSubs[1]=(TPad*) tempPad->GetPad(1);
pTPadSubs[2]=(TPad*) tempPad->GetPad(2);
tempPad=(TPad*) pEPad->GetPad(2);
tempPad->Divide(1,2);
pEPadSubs[1]=(TPad*) tempPad->GetPad(1);
pEPadSubs[2]=(TPad*) tempPad->GetPad(2);
tempPad=(TPad*) eEPad->GetPad(2);
tempPad->Divide(1,2);
eEPadSubs[1]=(TPad*) tempPad->GetPad(1);
eEPadSubs[2]=(TPad*) tempPad->GetPad(2);
TPad* allBasePads[10]={titlePad,pTPadSubs[0],pTPadSubs[1],pTPadSubs[2],eEPadSubs[0],eEPadSubs[1],eEPadSubs[2],pEPadSubs[0],pEPadSubs[1],pEPadSubs[2]};
for (int i = 1;i< 10;i++)
{
allBasePads[i]->SetGrid(1,1);
allBasePads[i]->SetTickx(1);
allBasePads[i]->SetTicky(1);
}
pTPadSubs[0]->SetMargin(0.13,0.1,.1,0.1);
pEPadSubs[0]->SetMargin(0.13,0.1,.1,0.1);
eEPadSubs[0]->SetMargin(0.13,0.1,.1,0.1);
pTPadSubs[1]->SetMargin(0.1,0.1,.1,0.1);
pEPadSubs[1]->SetMargin(0.1,0.1,.1,0.1);
eEPadSubs[1]->SetMargin(0.1,0.1,.1,0.1);
pTPadSubs[2]->SetMargin(0.1,0.1,.1,0.1);
pEPadSubs[2]->SetMargin(0.1,0.1,.1,0.1);
eEPadSubs[2]->SetMargin(0.1,0.1,.1,0.1);
//Plot pT
pTPadSubs[0]->cd(); drawPlot(detType, PLOTVAR_PT,PLOT_COMP);
pTPadSubs[1]->cd(); drawPlot(detType, PLOTVAR_PT,PLOT_RESID);
pTPadSubs[2]->cd(); drawPlot(detType, PLOTVAR_PT,PLOT_NORMRESID);
//Plot eE
eEPadSubs[0]->cd(); drawPlot(detType, PLOTVAR_EE,PLOT_COMP);
eEPadSubs[1]->cd(); drawPlot(detType, PLOTVAR_EE,PLOT_RESID);
eEPadSubs[2]->cd(); drawPlot(detType, PLOTVAR_EE,PLOT_NORMRESID);
//Plot pE
pEPadSubs[0]->cd(); drawPlot(detType, PLOTVAR_PE,PLOT_COMP);
pEPadSubs[1]->cd(); drawPlot(detType, PLOTVAR_PE,PLOT_RESID);
pEPadSubs[2]->cd(); drawPlot(detType, PLOTVAR_PE,PLOT_NORMRESID);
}
else
{
TPad* allBasePads[4]={titlePad,pTPad,eEPad,pEPad};
for (int i = 1;i< 4;i++)
{
allBasePads[i]->SetGrid(1,1);
allBasePads[i]->SetTickx(1);
allBasePads[i]->SetTicky(1);
}
pTPad->SetMargin(0.13,0.1,.1,0.1);
pEPad->SetMargin(0.13,0.1,.1,0.1);
eEPad->SetMargin(0.13,0.1,.1,0.1);
pTPad->cd(); drawPlot(detType, PLOTVAR_PT,PLOT_COMP);
eEPad->cd(); drawPlot(detType, PLOTVAR_EE,PLOT_COMP);
pEPad->cd(); drawPlot(detType, PLOTVAR_PE,PLOT_COMP);
}
mainPad->SetFillColor(kGray);
///Make Title box
titlePad->cd();
titlePad->SetFillColor(kGray);
TPaveText* titleBox = makeTitleBox( detType);
titleBox->Draw();
TString coTypeString;
if(detType==DET_EP)
{
coTypeString="EP";
}
else if(detType==DET_EPG)
{
coTypeString="EPG";
}
else if(detType==DET_EPBG)
{
coTypeString="EPBG";
}
TString imagePath=GRAPHS_DIR;
imagePath+="layouts/PlotLayout";
if(numMC>0) imagePath+="_MC";
for (int i = 0;i< numMC;i++)
{
imagePath+=TString("_")+mc[i]->GetName();
}
if(numExp>0) imagePath+="_Exp";
for (int i = 0;i< numExp;i++)
{
imagePath+=TString("_")+exp[i]->GetName();
}
imagePath+="_"+coTypeString+"_EPPlots.pdf";
theCanvas->SaveAs(imagePath);
delete theCanvas;
clearPlotHists();
}