CMSLLGAnalysisCode/LLGAnalysis.h at master · ssilvado/CMSLLGAnalysisCode · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
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
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
#include <iostream>
#include <fstream>
#include <TFile.h>
#include <TChain.h>
#include <TCanvas.h>
#include <TROOT.h>
#include <stdlib.h>
#include <TStyle.h>
#include <TH2D.h>
#include <TLegend.h>
#include <TColor.h>
#include <math.h>
#include <THnSparse.h>
#include <TChain.h>
#include <map>
#include <string>
#include <vector>
#include <TRandom3.h>
#include "THnSparse.h"
#include "TF1.h"
#include "TSystem.h"

using namespace std;

class LLGAnalysis {
    public:
        static LLGAnalysis* GetInstance( char* configFileName );
        ~LLGAnalysis() {}

        vector<double> CalculateVertex( vector<double> x, vector<double> y, vector<double> z, vector<double> weight, vector<int> charge, vector<double> distance, unsigned int &nConsidered, double &weightednConsidered, vector<double> &error );

        bool Init();
        void RunEventLoop( int nEventsMax = -1);
        void FinishRun();

    private:
        static LLGAnalysis* _instance;

        void makeHist( string nametitle, int nbinsx, double xmin, double xmax, int nbinsy, double ymin, double ymax, string xtitle, string ytitle, string ztitle, string drawOption = "", double xAxisOffset = 1., double yAxisOffset = 1.2, double zAxisOffset = 1. );
        void makeHist( string nametitle, int nbins, double xmin, double xmax, string xtitle, string ytitle, string drawOption = "", double xAxisOffset = 1., double yAxisOffset = 1.2 );
        void setStyle(double ytoff = 1.0, bool marker = true, double left_margin = 0.15);
        void MakeEfficiencyPlot( TH1D hpass, TH1D htotal, TCanvas *c, string triggerName = "");
        void FillEfficiencyHistograms();

        void SetupSignalRegion();
        void SignalRegionSelection();
        void SetupTTJetsCR();
        void TTJetsCRSelection();
        //void SetupSingleTopCR();
        //void SingleTopCRSelection();
        // INSERT YOUR SELECTION HERE


    private:
        LLGAnalysis() {}
        LLGAnalysis( char* configFileName );

        map<string, int>        _cutFlow;
        map<string, TH1D>       _histograms1D;
        map<string, TH2D>       _histograms2D;
        map<string, string>     _histograms1DDrawOptions;
        map<string, string>     _histograms2DDrawOptions;
        vector<string>          _inputFileNames;
        string                  _inputTreeName;
        TChain                  *_inputTree;
        string					_outputFileName;

        vector<double> *recoJet_pt;
        vector<double> *recoJet_phi;
        vector<double> *recoJet_eta;
        vector<double> *recoJet_btag_combinedInclusiveSecondaryVertexV2BJetTags;
        vector<double> *recoJet_btag_jetBProbabilityBJetTags;
        vector<double> *recoJet_btag_jetProbabilityBJetTags;
        vector<double> *recoJet_btag_trackCountingHighPurBJetTags;
        vector<double> *recoJet_btag_trackCountingHighEffBJetTags;
        vector<double> *muon_px;
        vector<double> *muon_py;
        vector<double> *muon_pz;
        vector<double> *muon_phi;
        vector<double> *muon_eta;
        vector<double> *electron_px;
        vector<double> *electron_py;
        vector<double> *electron_pz;
        vector<double> *electron_phi;
        vector<double> *electron_eta;
        vector<int> *triggerBits;
        vector<string> *triggerNames;
        vector<vector<double> > *recoJet_constVertex_x;
        vector<vector<double> > *recoJet_constVertex_y;
        vector<vector<double> > *recoJet_constVertex_z;
        vector<vector<double> > *recoJet_const_pt;
        vector<vector<double> > *recoJet_const_closestVertex_dxy;
        vector<vector<double> > *recoJet_const_closestVertex_dz;
        vector<vector<double> > *recoJet_const_closestVertex_d;
        vector<vector<int> > *recoJet_const_charge;

        vector<double> *vertex_x;
        vector<double> *vertex_y;
        vector<double> *vertex_z;
        vector<double> *secVertex_x;
        vector<double> *secVertex_y;
        vector<double> *secVertex_z;
        vector<double> *vertex_nTracks;
        vector<double> *vertex_pt;

        double met;


        double evtWeight;
        double JET_PT_CUT_SV;
        double JET_PT_CUT_PV;
        double JET_ETA_CUT;
        double MUON_PT_CUT;
        double ELECTRON_PT_CUT;
        double MET_CUT;
        std::string SELECTION;
        int N_JET_CUT;
        int N_BJET_CUT;

        // the total number of events per sample (take it from DAS!)
        double PROC_NTOT;
        // the total cross section for the process in pb
        double PROC_XSEC;
        // the target luminosity in fb-1
        double TARGET_LUMI;

        bool applyEventWeights;


        vector<string> _plotFormats;

};