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.3 by tboccali, Mon Jul 18 13:14:05 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  
165  
166   public:
126  ParticleMCInfo mcH;
127  ParticleMCInfo mcW;
128  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;
135 <  ParticleMCInfo mcB;
136 <  ParticleMCInfo mcC;    
137 <    
138 <  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;

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines