VHbbDataFormats/interface/VHbbEvent.h

#ifndef VHbbEvent__H 
#define VHbbEvent__H 

#include <TLorentzVector.h>
#include <TVector2.h>
#include <vector>

class VHbbEvent{
 public:


  class SimpleJet {
    public:
    SimpleJet(): flavour(-99), tche(-99), tchp(-99), jpb(-99), jp(-99), 
      ssvhe(-99), csv(-99), csvmva(-99), ntracks(-99), charge(-99),
      bestMCid(-99), bestMCmomid(-99){}
  public:
    double Pt() {return p4.Pt();}
    int flavour;
    float tche,tchp, jpb,jp , ssvhe, csv, csvmva;
    int ntracks;
    float charge;
    TLorentzVector p4;
    TLorentzVector chargedTracksFourMomentum;

    int bestMCid, bestMCmomid;
    // new
    TVector2 tVector;
  };


  class HardJet {
  public:
    HardJet(): constituents(-99){}
  public:
    int constituents;
    TLorentzVector p4;
    std::vector<TLorentzVector> subFourMomentum;
    std::vector<float> etaSub, phiSub;
  };


  class METInfo {
  public:
    METInfo(): sumEt(-99), metSig(-99), eLong(-99){}
  public:
    float sumEt, metSig, eLong;
    TLorentzVector p4;
  };

  class MuonInfo {
  public:
    MuonInfo(): charge(-99),tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfPhoIso(-99), pfNeuIso(-99),
      acop(-99), ipDb(-99), ipErrDb(-99), zPVPt(-99),zPVProb(-99), chi2(-99), globChi2(-99),
      cat(-99), nHits(-99), nPixelHits(-99), globNHits(-99),validMuStations(-99),
      mcId(-99), mcMomId(-99), mcgMomId(-99){}
  public:
    TLorentzVector p4;
    int charge;
    float tIso, eIso, hIso, pfChaIso,pfPhoIso,pfNeuIso,acop, ipDb, ipErrDb, zPVPt,zPVProb, chi2, globChi2;
    int cat, nHits, nPixelHits, globNHits, validMuStations;
    TLorentzVector mcFourMomentum;
    int mcId, mcMomId, mcgMomId;
  };

  class ElectronInfo {
  public:
    ElectronInfo() : scEta(-99), scPhi(-99), charge(-99), 
    tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfPhoIso(-99), pfNeuIso(-99), acop(-99),
      id95(-99),id85(-99),id80(-99),id70(-99),
      id95r(-99),id85r(-99),id80r(-99),id70r(-99),
      mcId(-99), mcMomId(-99), mcgMomId (-99){}
  public:
    TLorentzVector p4;
    float scEta, scPhi;
    int charge;
    float tIso, eIso, hIso, pfChaIso,pfPhoIso,pfNeuIso, acop;
    float  id95,id85,id80,id70,id95r, id85r,id80r, id70r;
    TLorentzVector mcFourMomentum;
    int mcId, mcMomId, mcgMomId;
  };

  class TauInfo{
  public:
    TauInfo()  : charge(-99), tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfPhoIso(-99), pfNeuIso(-99), acop(-99), 
    idbyIso(-99),idbyTrackIso(-99),idbyTaNCfrOnePercent(-99),
    idbyTaNCfrHalfPercent(-99), idbyTaNCfrQuarterPercent(-99), 
      idbyTaNCfrTenthPercent(-99), idbyTaNC(-99), mcId(-99), mcMomId(-99), mcgMomId(-99) {}
  public:
    TLorentzVector p4;
    int charge;
    float tIso, eIso, hIso,pfChaIso,pfPhoIso,pfNeuIso, acop;
    float  idbyIso,idbyTrackIso,idbyTaNCfrOnePercent,idbyTaNCfrHalfPercent, idbyTaNCfrQuarterPercent, idbyTaNCfrTenthPercent, idbyTaNC;
    TLorentzVector mcFourMomentum;
    int mcId, mcMomId, mcgMomId;
  };


  class DiMuonInfo  {
  public:
    TLorentzVector p4;
    MuonInfo daughter1, daughter2;
  };


  class DiElectronInfo {
  public:
    TLorentzVector p4;
    ElectronInfo daughter1, daughter2;
  };


 public:
  std::vector<SimpleJet> simpleJets;
  std::vector<SimpleJet> simpleJets2; //???
  std::vector<SimpleJet> subJets; //???
  std::vector<HardJet> hardJets;
  
 METInfo calomet;
  METInfo tcmet;
  METInfo pfmet;
  
  std::vector<MuonInfo> muInfo;
  std::vector<ElectronInfo> eleInfo;
  std::vector<TauInfo> tauInfo;
  
  std::vector<DiMuonInfo> diMuonInfo;
  std::vector<DiElectronInfo> diElectronInfo;
  
  //    HiggsCandidate dijetHiggs;
  //    HiggsCandidate fatjetHiggs;
  
  // to be decided.. one single or many V? 
  //  VCandidate vcand;
  //  WCandidate wele;
  //  WCandidate wmu;
  //WCandidate wtau;
  //  ZCandidate zmu;
  //  ZCandidate zele;
  //ZInvCandidate znunu;
  
  //
  // really needed????
  //

};
#endif
Revision:	1.7
Committed:	Wed Aug 10 08:17:27 2011 UTC (13 years, 8 months ago) by arizzi
Content type:	text/plain
Branch:	MAIN
Changes since 1.6:	+10 -9 lines
Log Message:	new isolation
#	User	Rev	Content
1	tboccali	1.1	#ifndef VHbbEvent__H
2			#define VHbbEvent__H
3
4			#include <TLorentzVector.h>
5			#include <TVector2.h>
6			#include <vector>
7
8			class VHbbEvent{
9			public:
10
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),
16			bestMCid(-99), bestMCmomid(-99){}
17			public:
18	tboccali	1.6	double Pt() {return p4.Pt();}
19	tboccali	1.1	int flavour;
20			float tche,tchp, jpb,jp , ssvhe, csv, csvmva;
21			int ntracks;
22			float charge;
23	tboccali	1.6	TLorentzVector p4;
24	tboccali	1.1	TLorentzVector chargedTracksFourMomentum;
25
26			int bestMCid, bestMCmomid;
27			// new
28			TVector2 tVector;
29			};
30
31
32			class HardJet {
33			public:
34			HardJet(): constituents(-99){}
35			public:
36			int constituents;
37	tboccali	1.6	TLorentzVector p4;
38	tboccali	1.1	std::vector<TLorentzVector> subFourMomentum;
39			std::vector<float> etaSub, phiSub;
40			};
41
42
43			class METInfo {
44			public:
45			METInfo(): sumEt(-99), metSig(-99), eLong(-99){}
46			public:
47			float sumEt, metSig, eLong;
48	tboccali	1.6	TLorentzVector p4;
49	tboccali	1.1	};
50
51			class MuonInfo {
52			public:
53	arizzi	1.7	MuonInfo(): charge(-99),tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfPhoIso(-99), pfNeuIso(-99),
54	tboccali	1.1	acop(-99), ipDb(-99), ipErrDb(-99), zPVPt(-99),zPVProb(-99), chi2(-99), globChi2(-99),
55			cat(-99), nHits(-99), nPixelHits(-99), globNHits(-99),validMuStations(-99),
56			mcId(-99), mcMomId(-99), mcgMomId(-99){}
57			public:
58	tboccali	1.6	TLorentzVector p4;
59	tboccali	1.1	int charge;
60	arizzi	1.7	float tIso, eIso, hIso, pfChaIso,pfPhoIso,pfNeuIso,acop, ipDb, ipErrDb, zPVPt,zPVProb, chi2, globChi2;
61	tboccali	1.1	int cat, nHits, nPixelHits, globNHits, validMuStations;
62			TLorentzVector mcFourMomentum;
63			int mcId, mcMomId, mcgMomId;
64			};
65
66			class ElectronInfo {
67			public:
68			ElectronInfo() : scEta(-99), scPhi(-99), charge(-99),
69	arizzi	1.7	tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfPhoIso(-99), pfNeuIso(-99), acop(-99),
70			id95(-99),id85(-99),id80(-99),id70(-99),
71			id95r(-99),id85r(-99),id80r(-99),id70r(-99),
72			mcId(-99), mcMomId(-99), mcgMomId (-99){}
73	tboccali	1.1	public:
74	tboccali	1.6	TLorentzVector p4;
75	tboccali	1.1	float scEta, scPhi;
76			int charge;
77	arizzi	1.7	float tIso, eIso, hIso, pfChaIso,pfPhoIso,pfNeuIso, acop;
78			float id95,id85,id80,id70,id95r, id85r,id80r, id70r;
79	tboccali	1.1	TLorentzVector mcFourMomentum;
80			int mcId, mcMomId, mcgMomId;
81			};
82
83			class TauInfo{
84			public:
85	arizzi	1.7	TauInfo() : charge(-99), tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfPhoIso(-99), pfNeuIso(-99), acop(-99),
86	tboccali	1.1	idbyIso(-99),idbyTrackIso(-99),idbyTaNCfrOnePercent(-99),
87			idbyTaNCfrHalfPercent(-99), idbyTaNCfrQuarterPercent(-99),
88			idbyTaNCfrTenthPercent(-99), idbyTaNC(-99), mcId(-99), mcMomId(-99), mcgMomId(-99) {}
89			public:
90	tboccali	1.6	TLorentzVector p4;
91	tboccali	1.1	int charge;
92	arizzi	1.7	float tIso, eIso, hIso,pfChaIso,pfPhoIso,pfNeuIso, acop;
93	tboccali	1.1	float idbyIso,idbyTrackIso,idbyTaNCfrOnePercent,idbyTaNCfrHalfPercent, idbyTaNCfrQuarterPercent, idbyTaNCfrTenthPercent, idbyTaNC;
94			TLorentzVector mcFourMomentum;
95			int mcId, mcMomId, mcgMomId;
96			};
97
98
99			class DiMuonInfo {
100			public:
101	tboccali	1.6	TLorentzVector p4;
102	tboccali	1.1	MuonInfo daughter1, daughter2;
103			};
104
105
106			class DiElectronInfo {
107			public:
108	tboccali	1.6	TLorentzVector p4;
109	tboccali	1.1	ElectronInfo daughter1, daughter2;
110			};
111
112
113
114			public:
115			std::vector<SimpleJet> simpleJets;
116			std::vector<SimpleJet> simpleJets2; //???
117			std::vector<SimpleJet> subJets; //???
118			std::vector<HardJet> hardJets;
119
120	tboccali	1.5	METInfo calomet;
121	tboccali	1.1	METInfo tcmet;
122			METInfo pfmet;
123
124			std::vector<MuonInfo> muInfo;
125			std::vector<ElectronInfo> eleInfo;
126			std::vector<TauInfo> tauInfo;
127
128			std::vector<DiMuonInfo> diMuonInfo;
129			std::vector<DiElectronInfo> diElectronInfo;
130
131			// HiggsCandidate dijetHiggs;
132			// HiggsCandidate fatjetHiggs;
133
134			// to be decided.. one single or many V?
135			// VCandidate vcand;
136			// WCandidate wele;
137			// WCandidate wmu;
138			//WCandidate wtau;
139			// ZCandidate zmu;
140			// ZCandidate zele;
141			//ZInvCandidate znunu;
142
143			//
144			// really needed????
145			//
146
147			};
148			#endif