void go_plot_contours(
char *device, char *title, char *xvar, char *yvar, char *topline,
double **data, double xmin, double xmax, double ymin, double ymax,
double dx, double dy, long nx, long ny,
double *contour_level, long n_contours, long contour_label_interval,
long contour_label_offset, long layout[2], long ix, long iy,
char *shapes, int *pen, long flags, long pause_interval,
SHAPE_DATA *shape, long nshapes, unsigned long long tsetFlags, double xlableScale, double ylabelScale)
{
double map[4], average, spread;
double pmin, pmax, qmin, qmax;
double wpmin, wpmax=0, wqmin, wqmax=0;
long title_at_top=0;
char newLabel[50];
time_t timeValue;
if (!(flags&DEVICE_DEFINED))
change_term(device, strlen(device));
if ((ix == 0) && (iy == 0))
graphics_on();
if (layout[0] && layout[1]) {
wpmin = (1.0*ix+0)/layout[0];
wpmax = (1.0*ix+1)/layout[0];
wqmin = (layout[1]*1.0 - 1 - iy)/layout[1];
wqmax = (layout[1]*1.0 - 0 - iy)/layout[1];
set_wspace(wpmin, wpmax, wqmin, wqmax);
pmin = (wpmax - wpmin) * .15 + wpmin;
pmax = wpmax - (wpmax - wpmin) * .1;
qmin = (wqmax - wqmin) * .17 + wqmin;
qmax = wqmax - (wqmax - wqmin) * .08;
set_pspace(pmin, pmax, qmin, qmax);
}
if (contour_label_interval<0) {
get_pspace(&pmin, &pmax, &qmin, &qmax);
set_pspace(pmin, 0.9 * wpmax, qmin, 0.85 * wqmax);
}
if (flags&TITLE_AT_TOP) {
/* adjust pspace to make room at top for the title */
get_pspace(&pmin, &pmax, &qmin, &qmax);
spread = qmax-qmin;
qmax -= 0.04*spread;
qmin -= 0.04*spread;
/* qmax -= 0.04;
qmin -= 0.04;*/
set_pspace(pmin, pmax, qmin, qmax);
title_at_top = 1;
}
set_clipping(1, 1, 1);
if (flags&EQUAL_ASPECT1)
set_aspect(1.0L);
else if (flags&EQUAL_ASPECT_1)
set_aspect(-1.0L);
set_linetype(pen[0]);
get_mapping(map, map+1, map+2, map+3);
if (map[0]==map[1]) {
average = (xmax+xmin)/2;
spread = (xmax-xmin)/2;
map[0] = average - spread*1.05;
map[1] = average + spread*1.05;
}
if (map[2]==map[3]) {
average = (ymax+ymin)/2;
spread = (ymax-ymin)/2;
map[2] = average - spread*1.05;
map[3] = average + spread*1.05;
}
set_mapping(map[0], map[1], map[2], map[3]);
get_pspace(&pmin, &pmax, &qmin, &qmax);
get_wspace(&wpmin, &wpmax, &wqmin, &wqmax);
if (!(flags&NO_BORDER)) {
border();
if (!(flags&NO_SCALES)) {
if (tsetFlags&TICKSET_XTIME) {
timeValue = map[0];
if (timeValue!=DBL_MAX)
sprintf(newLabel, "Time starting %s", ctime(&timeValue));
else
sprintf(newLabel, "Undefined time values!");
delete_chars(newLabel, "\n");
makeTimeScales(0, 0.02, 0, map[2], map[2], newLabel,
0, 0.67*(qmin-wqmin)*(map[3]-map[2])/(qmax-qmin),
0, 1, 1);
} else {
make_scales_with_label(0, 0, 1, 0.0, 0.02, 0.0,
0.0, 1.0, 0, 0,
0, 0, 0, 0, xvar, 0, 0, xlableScale);
}
if (!(flags&NO_YSCALES)) {
if (tsetFlags&TICKSET_YTIME) {
timeValue = map[2];
if (timeValue!=DBL_MAX)
sprintf(newLabel, "Time starting %s", ctime(&timeValue));
else
sprintf(newLabel, "Undefined time values!");
delete_chars(newLabel, "\n");
makeTimeScales(1, 0.0125, 0, map[0], map[0], newLabel,
0, (pmin-wpmin)*(map[1]-map[0])/(pmax-pmin),
0, 1, 1);
} else {
make_scales_with_label(1, 0, 1, 0.0, 0.0125, 0.0,
0.0, 1.0, 0, 0,
0, 0, 0, 0, yvar, 0, 0, ylabelScale);
}
}
}
}
if (!(flags&NO_LABELS)) {
plotTitle(title, 1, title_at_top, 1.0, 0.0, 0, 0);
if (contour_label_interval>=0)
plotTitle(topline, 0, 0, 1.0, 0.0, 0, 0);
}
set_linetype(pen[1]);
draw_contours(data, xmin, xmax, ymin, ymax, nx, ny, contour_level, n_contours);
if (contour_label_interval!=0) {
set_linetype(pen[2]);
label_contours(data, nx, ny, contour_level, n_contours,
contour_label_interval, contour_label_offset,
xmin, (double)dx, ymin, (double)dy);
}
set_linetype(pen[0]);
if (nshapes)
PlotShapesData(shape, nshapes, xmin, xmax, ymin, ymax);
if (flags&DATE_STAMP)
time_date_stamp();
}
int main(int argc, char** argv){
std::cout << "Starting application" << std::endl;
TApplication* app = new TApplication("Analyse",&argc,argv);
//Set plot style for own root design
PlotStyle* s = new PlotStyle();
CanvasManager* cManager = new CanvasManager("combinedPlots");
//#######################
//# Creating Trees and Analyse objects
//#######################
std::vector<AnalysisWrapper*> wrapperVect;
wrapperVect.push_back(new AnalysisWrapper("deltaPhiGun/DeltaPhiGunPt100dPhi3_0/TrigTestDeltaPhiGun.root","DeltaPhiGunPt100dR3_0",10.,true));
// wrapperVect.push_back(new AnalysisWrapper("deltaPhiGun/DeltaPhiGunPt100dPhi0_3/TrigTestDeltaPhiGun.root","DeltaPhiGunPt100dR0_3",10.,true));
// wrapperVect.push_back(new AnalysisWrapper("deltaPhiGun/DeltaPhiGunPt100dPhi0_1/TrigTestDeltaPhiGun.root","DeltaPhiGunPt100dR0_1",10.,true));
// wrapperVect.push_back(new AnalysisWrapper("deltaPhiGun/DeltaPhiGunPt100dPhi0_05/TrigTestDeltaPhiGun.root","DeltaPhiGunPt100dR0_05",10.,true));
// wrapperVect.push_back(new AnalysisWrapper("deltaPhiGun/DeltaPhiGunPt100dPhi0_01/TrigTestDeltaPhiGun.root","DeltaPhiGunPt100dR0_01",10.,true));
// wrapperVect.push_back(new AnalysisWrapper("deltaPhiGun/DeltaPhiGunPt100dPhi0_005/TrigTestDeltaPhiGun.root","DeltaPhiGunPt100dR0_005",10.,true));
wrapperVect[0]->analyseBti();
wrapperVect[0]->analyseTraco();
wrapperVect[0]->analyseGenParticles();
wrapperVect[0]->savePlots();
//Create a vector with colors to access during scenario iterations
std::vector<int> colorVector;
colorVector.push_back(kBlack);
colorVector.push_back(kBlue);
colorVector.push_back(kViolet);
colorVector.push_back(kRed);
colorVector.push_back(kOrange);
colorVector.push_back(kGreen);
colorVector.push_back(kTeal);
colorVector.push_back(kCyan - 1);
//Create a vector with scenarios to access during scenario iterations
std::vector<TString> legendNameVector;
legendNameVector.push_back(TString("#Delta#Phi 3.0"));
legendNameVector.push_back(TString("#Delta#Phi 0.3"));
legendNameVector.push_back(TString("#Delta#Phi 0.1"));
legendNameVector.push_back(TString("#Delta#Phi 0.05"));
legendNameVector.push_back(TString("#Delta#Phi 0.01"));
legendNameVector.push_back(TString("#Delta#Phi 0.005"));
std::vector<std::vector<TH1*> > btiScenariosBestCase;
std::vector<std::vector<TH1*> > btiScenariosHtrg;
for(unsigned int i = 0 ; i < wrapperVect.size(); i++ ){
btiScenariosBestCase.push_back(wrapperVect[i]->analyseBtiTrigPerStatAndSlBestCase());
btiScenariosHtrg.push_back(wrapperVect[i]->analyseBtiTrigPerStatAndSlHtrg());
}
TCanvas* btiBestCaseCanvas = cManager->getDividedCanvas(1,1);
btiBestCaseCanvas->SetName("btiTrigsPerStatAndSlTogetherBestCase");
TCanvas* btiHtrgCanvas = cManager->getDividedCanvas(1,1);
btiHtrgCanvas->SetName("btiTrigsPerStatAndSlTogetherHtrg");
//Get the size of the first vector stored in the scenarios
for (int i = 0 ; i < btiScenariosBestCase[0].size()/8. ; i++){
int scenarioCounter = 0;
TLegend* lBestCase = new TLegend(0.5,.525,0.7,0.9);
TLegend* lBestCase2 = new TLegend(0.7,.525,0.9,0.9);
TLegend* lHtrg = new TLegend(0.5,.525,0.7,0.9);
TLegend* lHtrg2 = new TLegend(0.7,.525,0.9,0.9);
btiBestCaseCanvas->cd(i+1)->SetLogy();
//Build the frame
TString frameName = "BTI Trigger (Best Case) f#ddot{u}r Station ";
frameName += i+1;
frameName += ";# BTI Trigger;Relative H#ddot{a}ufigkeit";
TH2D* frame = new TH2D("frameBtiBestCase",frameName,1,-0.5,10.5,1,0.000001,1);
frame->SetStats(kFALSE);
frame->Draw();
for (std::vector<std::vector<TH1*> >::const_iterator scenarioIt = btiScenariosBestCase.begin();
scenarioIt != btiScenariosBestCase.end() ; scenarioIt++){
//Loop over the scenarios
scenarioIt->at(2*i)->SetLineColor(colorVector[scenarioCounter]);
scenarioIt->at(2*i)->Scale(1/scenarioIt->at(2*i)->Integral());
scenarioIt->at(2*i)->Draw("same");
scenarioIt->at(2*i+1)->SetLineColor(colorVector[scenarioCounter]);
scenarioIt->at(2*i+1)->SetLineStyle(2);
scenarioIt->at(2*i+1)->Scale(1/scenarioIt->at(2*i+1)->Integral());
scenarioIt->at(2*i+1)->Draw("same");
lBestCase->AddEntry(scenarioIt->at(2*i),(legendNameVector[scenarioCounter] + " SL 1").Data());
lBestCase2->AddEntry(scenarioIt->at(2*i+1),(legendNameVector[scenarioCounter] + " SL 3").Data());
scenarioCounter++;
}
lBestCase->Draw();
lBestCase2->Draw();
btiBestCaseCanvas->Update();
TString plotTitle("BTI Trigger (Best Case) f#ddot{u}r Station ");
plotTitle += i+1;
(new PlotTitle(plotTitle.Data()))->Draw();
scenarioCounter = 0;
btiHtrgCanvas->cd(i+1)->SetLogy();
//Build the frame
frameName = "BTI Trigger (HTRG) f#ddot{u}r Station ";
frameName += i+1;
frameName += ";# BTI Trigger;Relative H#ddot{a}ufigkeit";
frame = new TH2D("frameBtiHtrg",frameName,1,-0.5,25.5,1,0.000001,1);
frame->SetStats(kFALSE);
frame->Draw();
for (std::vector<std::vector<TH1*> >::const_iterator scenarioIt = btiScenariosHtrg.begin();
scenarioIt != btiScenariosHtrg.end() ; scenarioIt++){
//Loop over the scenarios
scenarioIt->at(2*i)->SetLineColor(colorVector[scenarioCounter]);
scenarioIt->at(2*i)->Scale(1/scenarioIt->at(2*i)->Integral());
scenarioIt->at(2*i)->Draw("same");
scenarioIt->at(2*i+1)->SetLineColor(colorVector[scenarioCounter]);
scenarioIt->at(2*i+1)->SetLineStyle(2);
scenarioIt->at(2*i+1)->Scale(1/scenarioIt->at(2*i+1)->Integral());
scenarioIt->at(2*i+1)->Draw("same");
lHtrg->AddEntry(scenarioIt->at(2*i),(legendNameVector[scenarioCounter] + " SL 1").Data());
lHtrg2->AddEntry(scenarioIt->at(2*i+1),(legendNameVector[scenarioCounter] + " SL 3").Data());
scenarioCounter++;
}
lHtrg->Draw();
lHtrg2->Draw();
btiHtrgCanvas->Update();
plotTitle = "BTI Trigger (HTRG) f#ddot{u}r Station ";
plotTitle += i+1;
(new PlotTitle(plotTitle.Data()))->Draw();
}
cManager->addCanvas("btiTrigsPerStatAndSlTogetherBestCase",btiBestCaseCanvas);
cManager->addCanvas("btiTrigsPerStatAndSlTogetherHtrg",btiHtrgCanvas);
//################################
//### plot TRACO triggers together
//################################
std::vector<std::vector<TH1*> > tracoScenariosBestCase;
std::vector<std::vector<TH1*> > tracoScenariosHtrg;
for(unsigned int i = 0 ; i < wrapperVect.size(); i++ ){
tracoScenariosBestCase.push_back(wrapperVect[i]->analyseTracoTrigPerStationBestCase());
tracoScenariosHtrg.push_back(wrapperVect[i]->analyseTracoTrigPerStationHtrig());
}
TCanvas* tracoBestCaseCanvas = cManager->getDividedCanvas(1,1);
tracoBestCaseCanvas->SetName("tracoTrigsPerStatTogetherBestCase");
TCanvas* tracoHtrgCanvas = cManager->getDividedCanvas(1,1);
tracoHtrgCanvas->SetName("tracoTrigsPerStatTogetherHtrg");
//Get the size of the first vector stored in the scenarios
for (int i = 0 ; i < tracoScenariosBestCase[0].size()/4. ; i++){
int scenarioCounter = 0;
TLegend* lBestCase = new TLegend(0.75,.45,0.90,0.90);
TLegend* lHtrg = new TLegend(0.75,.45,0.90,0.90);
tracoBestCaseCanvas->cd(i+1)->SetLogy();
//Build the frame
TString frameName = "TRACO Trigger (Best Case) f#ddot{u}r Station ";
frameName += i+1;
frameName += ";# TRACO Trigger;Relative H#ddot{a}ufigkeit";
TH2D* frame = new TH2D("frameTracoBestCase",frameName,1,-0.5,7.5,1,0.000001,1);
frame->SetStats(kFALSE);
frame->Draw();
for (std::vector<std::vector<TH1*> >::const_iterator scenarioIt = tracoScenariosBestCase.begin();
scenarioIt != tracoScenariosBestCase.end() ; scenarioIt++){
//Loop over the scenarios
scenarioIt->at(i)->SetLineColor(colorVector[scenarioCounter]);
scenarioIt->at(i)->Scale(1/scenarioIt->at(i)->Integral());
scenarioIt->at(i)->Draw("same");
lBestCase->AddEntry(scenarioIt->at(i),legendNameVector[scenarioCounter]);
scenarioCounter++;
}
lBestCase->Draw();
tracoBestCaseCanvas->Update();
TString plotTitle("TRACO Trigger (Best Case) f#ddot{u}r Station ");
plotTitle += i+1;
(new PlotTitle(plotTitle.Data()))->Draw();
scenarioCounter = 0;
tracoHtrgCanvas->cd(i+1)->SetLogy();
//Build the frame
frameName = "TRACO Trigger (HTRG) f#ddot{u}r Station ";
frameName += i+1;
frameName += ";# TRACO Trigger;Relative H#ddot{a}ufigkeit";
frame = new TH2D("frameTracoHtrg",frameName,1,-0.5,7.5,1,0.000001,1);
frame->SetStats(kFALSE);
frame->Draw();
for (std::vector<std::vector<TH1*> >::const_iterator scenarioIt = tracoScenariosHtrg.begin();
scenarioIt != tracoScenariosHtrg.end() ; scenarioIt++){
//Loop over the scenarios
scenarioIt->at(i)->SetLineColor(colorVector[scenarioCounter]);
scenarioIt->at(i)->Scale(1/scenarioIt->at(i)->Integral());
scenarioIt->at(i)->Draw("same");
lHtrg->AddEntry(scenarioIt->at(i),legendNameVector[scenarioCounter]);
scenarioCounter++;
}
lHtrg->Draw();
tracoHtrgCanvas->Update();
plotTitle = "TRACO Trigger (HTRG) f#ddot{u}r Station ";
plotTitle += i+1;
(new PlotTitle(plotTitle.Data()))->Draw();
}
cManager->addCanvas("tracoTrigsPerStatTogetherBestCase",tracoBestCaseCanvas);
cManager->addCanvas("tracoTrigsPerStatTogetherHtrg",tracoHtrgCanvas);
cManager->storePlots();
int scenarioCounter = 0;
std::cout << "N Entries in bti Best case Scenarios:" << std::endl;
for (std::vector<std::vector<TH1*> >::const_iterator scenarioIt = btiScenariosBestCase.begin();
scenarioIt != btiScenariosBestCase.end() ; scenarioIt++) {
std::cout << "\t" << legendNameVector[scenarioCounter] << ":\t\t" << scenarioIt->at(0)->GetEntries() << std::endl;
scenarioCounter++;
}
std::cout << std::endl;
scenarioCounter = 0;
std::cout << "N Entries in bti HTRG Scenarios:" << std::endl;
for (std::vector<std::vector<TH1*> >::const_iterator scenarioIt = btiScenariosHtrg.begin();
scenarioIt != btiScenariosHtrg.end() ; scenarioIt++) {
std::cout << "\t" << legendNameVector[scenarioCounter] << ":\t\t" << scenarioIt->at(0)->GetEntries() << std::endl;
scenarioCounter++;
}
std::cout << std::endl;
scenarioCounter = 0;
std::cout << "N Entries in traco HTRIG Scenarios:" << std::endl;
for (std::vector<std::vector<TH1*> >::const_iterator scenarioIt = tracoScenariosHtrg.begin();
scenarioIt != tracoScenariosHtrg.end() ; scenarioIt++) {
// std::cout << TString::Format("%s:\t%6d",legendNameVector[scenarioCounter].Data(),scenarioIt->at(0)->GetEntries()).Data() << std::endl;
std::cout << "\t" << legendNameVector[scenarioCounter] << ":\t\t" << scenarioIt->at(0)->GetEntries() << std::endl;
scenarioCounter++;
}
std::cout << std::endl;
scenarioCounter = 0;
std::cout << "N Entries in traco BestCase Scenarios:" << std::endl;
for (std::vector<std::vector<TH1*> >::const_iterator scenarioIt = tracoScenariosBestCase.begin();
scenarioIt != tracoScenariosBestCase.end() ; scenarioIt++) {
std::cout << "\t" << legendNameVector[scenarioCounter] << ":\t\t" << scenarioIt->at(0)->GetEntries() << std::endl;
scenarioCounter++;
}
std::cout << std::endl;
delete btiBestCaseCanvas;
delete btiHtrgCanvas;
delete tracoBestCaseCanvas;
delete tracoHtrgCanvas;
std::cout << "All done!\nExit via ctrl+c..." << std::endl;
// app->Run();
delete cManager;
return 0;
}