| 1 |
#ifndef histo_h
|
| 2 |
#define histo_h
|
| 3 |
#include "triggerStudy.h"
|
| 4 |
//#include "plotRange.h"
|
| 5 |
//================================================================================= make_histograms
|
| 6 |
void susy::makeHistograms(TString suffix){
|
| 7 |
tools::enterFcn("make_histograms", 3);
|
| 8 |
|
| 9 |
// plotFakeRatio(myPlots["FakeRatio"]);
|
| 10 |
// myPlots["CutFlow"].addh1d(name_temp, title_temp, xtitle, ytitle, bin, low, up, value, weight);
|
| 11 |
//
|
| 12 |
// plot.addh2d("ScanYield"+suffix, "ScanYield"+suffix, "m_{0} [GeV/c^{2}]", "m_{1/2} [GeV/c^{2}]", "events", binScanM0, lowScanM0, upScanM0, binScanM12, lowScanM12, upScanM12, susyScanM0, susyScanM12, 1);
|
| 13 |
// TH1::SetDefaultSumw2(1);
|
| 14 |
|
| 15 |
tools::exitFcn("make_histograms", 3);
|
| 16 |
}
|
| 17 |
|
| 18 |
|
| 19 |
|
| 20 |
// //================================================================================= drawLine
|
| 21 |
// void susy::drawLine(Double_t x1, Double_t y1, Double_t x2, Double_t y2, Int_t w, Int_t c, Int_t st){
|
| 22 |
// tools::enterFcn("drawLine", 2);
|
| 23 |
// TLine* l = new TLine(x1, y1, x2, y2);
|
| 24 |
// l->SetLineWidth(w);
|
| 25 |
// l->SetLineColor(c);
|
| 26 |
// l->SetLineStyle(st);
|
| 27 |
// line.push_back(l);
|
| 28 |
// line.back()->Draw("SAME");
|
| 29 |
// tools::exitFcn("drawLine", 2);
|
| 30 |
// }
|
| 31 |
|
| 32 |
|
| 33 |
|
| 34 |
// //================================================================================= plotEventShape
|
| 35 |
// void susy::plotEventShape(plotMaker &plot){
|
| 36 |
// tools::enterFcn("plotEventShape", 3);
|
| 37 |
// plot.addh1d(it->first+"Aplanarity"+suffix, it->first+"Aplanarity"+suffix, it->first+" Aplanarity", "events", 50 , 0, 0.5, myShape.aplanarity());
|
| 38 |
// tools::exitFcn("plotEventShape", 3);
|
| 39 |
// }
|
| 40 |
|
| 41 |
// //================================================================================= saveHistograms
|
| 42 |
// void susy::saveHistograms(){
|
| 43 |
// tools::enterFcn("saveHistograms", 3);
|
| 44 |
// timer.Start("Save histograms");
|
| 45 |
// TString option = "RECREATE";
|
| 46 |
// for(map<TString, plotMaker>::iterator it=myPlots.begin(); it != myPlots.end(); ++it){
|
| 47 |
// it->second.savePlots(outputPath + outputName, option, it->first);
|
| 48 |
// option = "UPDATE";
|
| 49 |
// }
|
| 50 |
//
|
| 51 |
// //mergeAllOutputFiles(outputName, "objectsBefore objectsAfter diLepBefore diLepAfter kinematicsBefore kinematicsAfter RelIsoBefore RelIsoAfter cutFlow n1");
|
| 52 |
//
|
| 53 |
// //myPlots["CutFlow"].saveSingleFiles(myPlots["CutFlow"].h1d, 1,1,1, "output"); // 1,1,1 = *.root, *.eps, *.gif,
|
| 54 |
// //n1.savePlots(path + outputName, "UPDATE", "n1");
|
| 55 |
// //TPostScript *ps = new TPostScript("output.ps",112);
|
| 56 |
// //n1.h1d["ElCuts"]->Draw();
|
| 57 |
// //ps->Close();
|
| 58 |
// //gSystem->Exec("gs output.ps");
|
| 59 |
//
|
| 60 |
// TFile *outputFile = new TFile(outputPath + outputName, "UPDATE");
|
| 61 |
// outputData->Fill();
|
| 62 |
// outputData->Write();
|
| 63 |
// TLTree->Write();
|
| 64 |
// //TLTree->Close();
|
| 65 |
//
|
| 66 |
// outputFile->Close();
|
| 67 |
// delete outputFile;
|
| 68 |
// timer.Stop("Save histograms");
|
| 69 |
// tools::exitFcn("saveHistograms", 3);
|
| 70 |
// }
|
| 71 |
|
| 72 |
|
| 73 |
|
| 74 |
|
| 75 |
|
| 76 |
|
| 77 |
|
| 78 |
|
| 79 |
|
| 80 |
|
| 81 |
|
| 82 |
|
| 83 |
|
| 84 |
|
| 85 |
|
| 86 |
|
| 87 |
|
| 88 |
|
| 89 |
|
| 90 |
|
| 91 |
|
| 92 |
|
| 93 |
|
| 94 |
|
| 95 |
|
| 96 |
|
| 97 |
|
| 98 |
|
| 99 |
|
| 100 |
|
| 101 |
|
| 102 |
|
| 103 |
#endif
|
| 104 |
|
| 105 |
|
| 106 |
|
| 107 |
|
| 108 |
|
| 109 |
|
| 110 |
|
| 111 |
|
| 112 |
|
| 113 |
|
| 114 |
|
