ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/VHbbAnalysis/VHbbDataFormats/interface/VHbbEvent.h
(Generate patch)

Comparing UserCode/VHbbAnalysis/VHbbDataFormats/interface/VHbbEvent.h (file contents):
Revision 1.2 by tboccali, Mon Jul 18 12:59:51 2011 UTC vs.
Revision 1.23 by sethzenz, Wed Feb 22 16:46:54 2012 UTC

# Line 8 | Line 8
8   class VHbbEvent{
9   public:
10  
11 –  class ParticleMCInfo {
12 –  public:
13 –    ParticleMCInfo(): status(-99), momid(-99), gmomid(-99), charge(-1){}
14 –  public:
15 –    int status;
16 –    int momid;
17 –    int gmomid;
18 –    int charge;
19 –    TLorentzVector fourMomentum;
20 –    //    int ndau;
21 –    std::vector<int> dauid;
22 –    std::vector<TLorentzVector> dauFourMomentum;    
23 –  };
11  
12    class SimpleJet {
13      public:
14 <    SimpleJet(): flavour(-99), tche(-99), tchp(-99), jpb(-99), jp(-99),
15 <      ssvhe(-99), csv(-99), csvmva(-99), ntracks(-99), charge(-99),
14 >    SimpleJet(): flavour(-99), isSemiLept(-99), isSemiLeptMCtruth(-99),  tche(-99), tchp(-99), jpb(-99), jp(-99),
15 >      ssvhe(-99), csv(-99), csvmva(-99),
16 >      vtxMass(-99), vtx3dL(-99), vtx3deL(-99), vtxNTracks(-99), vtxProb(-99),
17 >      ntracks(-99), charge(-99),jecunc(-99),
18 >      chargedHadronEFraction(-99), neutralHadronEFraction(-99), chargedEmEFraction(-99), neutralEmEFraction(-99),nConstituents(-99),
19        bestMCid(-99), bestMCmomid(-99){}
20    public:
21 <    int flavour;
21 >    double Pt() {return p4.Pt();}
22 >    int flavour, isSemiLept, isSemiLeptMCtruth;
23      float tche,tchp, jpb,jp , ssvhe, csv, csvmva;
24 +    // CSV vertex info
25 +    float vtxMass, vtx3dL, vtx3deL, vtxNTracks, vtxProb;
26 +    std::vector<unsigned int> vtxTrackIds;
27      int ntracks;
28      float charge;
29 <    TLorentzVector fourMomentum;
29 >    TLorentzVector p4;
30      TLorentzVector chargedTracksFourMomentum;
31  
32 +    // jec uncertainty
33 +
34 +    float jecunc;
35 +    // energy fractions as in
36 +    /*    //chargedHadronEnergyFraction (relative to uncorrected jet energy)
37 +          float chargedHadronEnergyFraction() const {return chargedHadronEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
38 +          /// neutralHadronEnergyFraction (relative to uncorrected jet energy)
39 +          float neutralHadronEnergyFraction() const {return neutralHadronEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
40 +          /// chargedEmEnergyFraction (relative to uncorrected jet energy)
41 +          float chargedEmEnergyFraction()     const {return chargedEmEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
42 +          /// neutralEmEnergyFraction (relative to uncorrected jet energy)
43 +          float neutralEmEnergyFraction()     const {return neutralEmEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
44 +    */
45 +
46 +    float chargedHadronEFraction, neutralHadronEFraction, chargedEmEFraction, neutralEmEFraction;
47 +    int nConstituents;
48 +
49 +    //MC parton matching
50 +    TLorentzVector bestMCp4,bestMCp4mom;
51      int bestMCid, bestMCmomid;
52      // new
53      TVector2 tVector;
54 +    //Scale factor
55 +    float SF_CSVL;
56 +    float SF_CSVM;
57 +    float SF_CSVT;
58 +    float SF_CSVLerr;
59 +    float SF_CSVMerr;
60 +    float SF_CSVTerr;
61 +    
62    };
63  
64  
# Line 46 | Line 67 | class VHbbEvent{
67      HardJet(): constituents(-99){}
68    public:
69      int constituents;
70 <    TLorentzVector fourMomentum;
70 >    TLorentzVector p4;
71      std::vector<TLorentzVector> subFourMomentum;
72      std::vector<float> etaSub, phiSub;
73    };
# Line 57 | Line 78 | class VHbbEvent{
78      METInfo(): sumEt(-99), metSig(-99), eLong(-99){}
79    public:
80      float sumEt, metSig, eLong;
81 <    TLorentzVector fourMomentum;
81 >    TLorentzVector p4;
82    };
83  
84    class MuonInfo {
85    public:
86 <    MuonInfo(): charge(-99),tIso(-99), eIso(-99), hIso(-99),
86 >    MuonInfo(): charge(-99),tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99),pfChaPUIso(99), pfPhoIso(-99), pfNeuIso(-99),
87        acop(-99), ipDb(-99), ipErrDb(-99), zPVPt(-99),zPVProb(-99), chi2(-99), globChi2(-99),
88 <      cat(-99), nHits(-99), nPixelHits(-99), globNHits(-99),validMuStations(-99),
88 >      cat(-99), nValidTracker(-99), nValidPixel(-99), nMatches(-99),nHits(-99), nPixelHits(-99), globNHits(-99),validMuStations(-99),emEnergy(-99), hadEnergy(-99),
89        mcId(-99), mcMomId(-99), mcgMomId(-99){}
90    public:
91 <    TLorentzVector fourMomentum;
91 >    TLorentzVector p4;
92      int charge;
93 <    float tIso, eIso, hIso, acop, ipDb, ipErrDb, zPVPt,zPVProb, chi2, globChi2;
94 <    int cat, nHits, nPixelHits, globNHits, validMuStations;
93 >    float tIso, eIso, hIso, pfChaIso,pfChaPUIso,  pfPhoIso,pfNeuIso,acop, ipDb, ipErrDb, zPVPt,zPVProb, chi2, globChi2;
94 >    int cat, nValidTracker, nValidPixel, nMatches,nHits, nPixelHits, globNHits, validMuStations;
95 >    float emEnergy, hadEnergy;
96      TLorentzVector mcFourMomentum;
97      int mcId, mcMomId, mcgMomId;
98 +    std::vector< int > hltMatchedBits;
99    };
100  
101    class ElectronInfo {
102    public:
103      ElectronInfo() : scEta(-99), scPhi(-99), charge(-99),
104 <    tIso(-99), eIso(-99), hIso(-99),
105 <      acop(-99), id95(-99),id85(-99),id70(-99),id95r(-99),
106 <      id70r(-99), id85r(-99),mcId(-99), mcMomId(-99), mcgMomId (-99){}
104 >      tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99),pfChaPUIso(-99), pfPhoIso(-99), pfNeuIso(-99), acop(-99),
105 >      sihih(-99), Dphi(-99), Deta(-99), HoE(-99), convDist(-99), convDcot(-99), innerHits(-99),
106 >      isEB(false),isEE(false),
107 >      ipDb(-99), ipErrDb(-99),
108 >      id95(-99),id85(-99),id80(-99),id70(-99),
109 >      id95r(-99),id85r(-99),id80r(-99),id70r(-99),
110 >      mcId(-99), mcMomId(-99), mcgMomId (-99){}
111    public:
112 <    TLorentzVector fourMomentum;
112 >    TLorentzVector p4;
113      float scEta, scPhi;
114      int charge;
115 <    float tIso, eIso, hIso, acop;
116 <    float  id95,id85,id70,id95r, id70r, id85r;
115 >    float tIso, eIso, hIso, pfChaIso,pfChaPUIso, pfPhoIso,pfNeuIso, acop;
116 >    float sihih, Dphi, Deta, HoE, convDist, convDcot;
117 >    int innerHits;
118 >    bool isEB,isEE;
119 >    float       ipDb, ipErrDb;
120 >    float  id95,id85,id80,id70,id95r, id85r,id80r, id70r;
121      TLorentzVector mcFourMomentum;
122      int mcId, mcMomId, mcgMomId;
123 +    std::vector< int > hltMatchedBits;
124    };
125  
126    class TauInfo{
127    public:
128 <    TauInfo()  : charge(-99), tIso(-99), eIso(-99), hIso(-99), acop(-99),
128 >    TauInfo()  : charge(-99), tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfPhoIso(-99), pfNeuIso(-99), acop(-99),
129      idbyIso(-99),idbyTrackIso(-99),idbyTaNCfrOnePercent(-99),
130      idbyTaNCfrHalfPercent(-99), idbyTaNCfrQuarterPercent(-99),
131 <      idbyTaNCfrTenthPercent(-99), idbyTaNC(-99), mcId(-99), mcMomId(-99), mcgMomId(-99) {}
131 >      idbyTaNCfrTenthPercent(-99), idbyTaNC(-99), mcId(-99), mcMomId(-99), mcgMomId(-99),
132 >      idagainstElectronLoose(-99),  idagainstElectronMedium(-99),  idagainstElectronTight(-99),  idagainstMuonLoose(-99),  idagainstMuonTight(-99),  
133 >      idbyLooseIsolation(-99),  idbyMediumIsolation(-99),  idbyTightIsolation(-99),  idbyVLooseIsolation(-99),  iddecayModeFinding(-99),
134 >      isolationPFChargedHadrCandsPtSum(-999999.), isolationPFGammaCandsEtSum(-999999.),leadPFChargedHadrCandPt(-999999.),
135 >      NsignalPFChargedHadrCands (-99),NsignalPFGammaCands(-99) {}
136    public:
137 <    TLorentzVector fourMomentum;
137 >    TLorentzVector p4;
138      int charge;
139 <    float tIso, eIso, hIso, acop;
139 >    float tIso, eIso, hIso,pfChaIso,pfPhoIso,pfNeuIso, acop;
140      float  idbyIso,idbyTrackIso,idbyTaNCfrOnePercent,idbyTaNCfrHalfPercent, idbyTaNCfrQuarterPercent, idbyTaNCfrTenthPercent, idbyTaNC;
105 –    TLorentzVector mcFourMomentum;
141      int mcId, mcMomId, mcgMomId;
142 +    float idagainstElectronLoose, idagainstElectronMedium, idagainstElectronTight, idagainstMuonLoose, idagainstMuonTight;
143 +    float idbyLooseIsolation, idbyMediumIsolation, idbyTightIsolation, idbyVLooseIsolation, iddecayModeFinding;
144 +
145 +    float isolationPFChargedHadrCandsPtSum,isolationPFGammaCandsEtSum,leadPFChargedHadrCandPt;
146 +    int NsignalPFChargedHadrCands, NsignalPFGammaCands;
147 +    TLorentzVector mcFourMomentum;
148    };
149  
150  
151    class DiMuonInfo  {
152    public:
153 <    TLorentzVector fourMomentum;
153 >    TLorentzVector p4;
154      MuonInfo daughter1, daughter2;
155    };
156  
157  
158    class DiElectronInfo {
159    public:
160 <    TLorentzVector fourMomentum;
160 >    TLorentzVector p4;
161      ElectronInfo daughter1, daughter2;
162    };
163  
164  
124 –  class TriggerInfo {
125 –  public:
126 –    TriggerInfo() :  triggerMu9(-99),
127 –      triggerIsoMu9(-99),
128 –      triggerIsoMu13_3(-99),
129 –      triggerMu11(-99),
130 –      triggerDoubleMu3(-99),
131 –      triggerDoubleMu3_2(-99),
132 –      triggerMu15(-99),
133 –      triggerMu15_1(-99),
134 –      triggerDoubleElec10(-99),
135 –      triggerDoubleElec15_1(-99),
136 –      triggerDoubleElec17_1(-99),
137 –      triggerMet100_1(-99),
138 –      triggerSingleEle1(-99),
139 –      triggerSingleEle2(-99),
140 –      triggerSingleEle3(-99),
141 –      triggerSingleEle4(-99),
142 –      triggerBtagMu1(-99),
143 –      triggerBtagMu2(-99),
144 –      triggerBtagMu0(-99),
145 –      triggerBtagMu11(-99),
146 –      triggerBtagMuJet1(-99),
147 –      triggerBtagMuJet2(-99),
148 –      triggerBtagMuJet3(-99),
149 –      triggerBtagMuJet4(-99),
150 –      triggerIsoMu15(-99),
151 –      triggerIsoMu17v5(-99),
152 –      triggerIsoMu17v6(-99) {
153 –      for (unsigned int i=0; i< 500; ++i){
154 –        flag[i]= -99;
155 –      }
156 –    }
157 –  public:
158 –    int flag[500];
159 –    int triggerMu9,
160 –      triggerIsoMu9,
161 –      triggerIsoMu13_3,
162 –      triggerMu11,
163 –      triggerDoubleMu3,
164 –      triggerDoubleMu3_2,
165 –      triggerMu15,
166 –      triggerMu15_1,
167 –      triggerDoubleElec10,
168 –      triggerDoubleElec15_1,
169 –      triggerDoubleElec17_1,
170 –      triggerMet100_1,
171 –      triggerSingleEle1,
172 –      triggerSingleEle2,
173 –      triggerSingleEle3,
174 –      triggerSingleEle4,
175 –      triggerBtagMu1,
176 –      triggerBtagMu2,
177 –      triggerBtagMu0,
178 –      triggerBtagMu11,
179 –      triggerBtagMuJet1,
180 –      triggerBtagMuJet2,
181 –      triggerBtagMuJet3,
182 –      triggerBtagMuJet4,
183 –      triggerIsoMu15,
184 –      triggerIsoMu17v5,
185 –      triggerIsoMu17v6;
186 –  };
187 –
188 –  class PrimaryVertexInfo {
189 –  public:
190 –    TVector3 firstPVInPT2,firstPVInProb;
191 –  };
192 –
193 –  class PUInfo{
194 –  public:
195 –    PUInfo(): rho(-99) {}
196 –  public:
197 –    float rho;
198 –  };
199 –
165  
166   public:
202 –  ParticleMCInfo mcH;
203 –  ParticleMCInfo mcW;
204 –  ParticleMCInfo mcZ;
167    std::vector<SimpleJet> simpleJets;
168    std::vector<SimpleJet> simpleJets2; //???
169 +  std::vector<SimpleJet> simpleJets3;
170 +  std::vector<SimpleJet> simpleJets4; //???
171    std::vector<SimpleJet> subJets; //???
172 +  std::vector<SimpleJet> filterJets; //???
173    std::vector<HardJet> hardJets;
174    
175 <  ParticleMCInfo mcBbar;
211 <  ParticleMCInfo mcB;
212 <  ParticleMCInfo mcC;    
213 <  
214 <  PUInfo puInfo;
215 <  
216 <  METInfo calomet;
175 > METInfo calomet;
176    METInfo tcmet;
177    METInfo pfmet;
178 <  
178 >  METInfo mht;
179 >  METInfo metNoPU;
180 >  METInfo metCh;
181 >
182    std::vector<MuonInfo> muInfo;
183    std::vector<ElectronInfo> eleInfo;
184    std::vector<TauInfo> tauInfo;
185    
224 –  TriggerInfo triggerInfo;
225 –  
226 –  PrimaryVertexInfo pvInfo;
227 –  
186    std::vector<DiMuonInfo> diMuonInfo;
187    std::vector<DiElectronInfo> diElectronInfo;
188    

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines