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), isSemiLept(-99), isSemiLeptMCtruth(-99), SoftLeptpdgId(-99), SoftLeptIdlooseMu(-99), SoftLeptId95(-99), SoftLeptPt(-99), SoftLeptdR(-99), SoftLeptptRel(-99),  SoftLeptRelCombIso(-99),   tche(-99), tchp(-99), jpb(-99), jp(-99), 
      ssvhe(-99), csv(-99), csvmva(-99), csvivf(-99), cmva(-99),
      vtxMass(-99), vtx3dL(-99), vtx3deL(-99), vtxNTracks(-99), vtxProb(-99),
      ntracks(-99), charge(-99),jecunc(-99),
      chargedHadronEFraction(-99), neutralHadronEFraction(-99), chargedEmEFraction(-99), neutralEmEFraction(-99),nConstituents(-99), ptRaw(-99), ptLeadTrack(-99), jetArea(-99),
      bestMCid(-99), bestMCmomid(-99),
      puJetIdL(-99), puJetIdM(-99), puJetIdT(-99), puJetIdMva(-99) {}
  public:
    double Pt() {return p4.Pt();}
    int flavour, isSemiLept, isSemiLeptMCtruth, SoftLeptpdgId , SoftLeptIdlooseMu, SoftLeptId95;
    float SoftLeptPt, SoftLeptdR, SoftLeptptRel, SoftLeptRelCombIso; 
    float tche,tchp, jpb,jp , ssvhe, csv, csvmva,csvivf,cmva;
    // CSV vertex info
    float vtxMass, vtx3dL, vtx3deL, vtxNTracks, vtxProb;
    TVector3 vtxPosition;
    TLorentzVector  vtxP4;
    std::vector<unsigned int> vtxTrackIds, btagTrackIds, csvTrackIds;
    int csvNTracks, btagNTracks;
    int ntracks;
    float charge;
    TLorentzVector p4;
    TLorentzVector chargedTracksFourMomentum;

    // jec uncertainty

    float jecunc;
    // energy fractions as in 
    /*    //chargedHadronEnergyFraction (relative to uncorrected jet energy)
          float chargedHadronEnergyFraction() const {return chargedHadronEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
          /// neutralHadronEnergyFraction (relative to uncorrected jet energy)
          float neutralHadronEnergyFraction() const {return neutralHadronEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
          /// chargedEmEnergyFraction (relative to uncorrected jet energy)
          float chargedEmEnergyFraction()     const {return chargedEmEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
          /// neutralEmEnergyFraction (relative to uncorrected jet energy)
          float neutralEmEnergyFraction()     const {return neutralEmEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
    */

    float chargedHadronEFraction, neutralHadronEFraction, chargedEmEFraction, neutralEmEFraction;
    int nConstituents;
    float ptRaw,ptLeadTrack;  
    float jetArea;  

    //MC parton matching
    TLorentzVector bestMCp4,bestMCp4mom;
    int bestMCid, bestMCmomid;
    // new
    TVector2 tVector;
    //Scale factor
    float SF_CSVL;
    float SF_CSVM;
    float SF_CSVT; 
    float SF_CSVLerr;
    float SF_CSVMerr;
    float SF_CSVTerr; 
    // puJet Id
    float  puJetIdL, puJetIdM, puJetIdT, puJetIdMva;
    
  };


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

class RawJet {
public:
    RawJet(): Nconstituents(-99){}
public:
    int Nconstituents;
    std::vector<float> constituents_px;
    std::vector<float> constituents_py;
    std::vector<float> constituents_pz;
    std::vector<float> constituents_e;
    std::vector<float> constituents_pdgId;    
    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),pfChaPUIso(99), pfPhoIso(-99), pfNeuIso(-99),
      acop(-99), ipDb(-99), ipErrDb(-99), zPVPt(-99),zPVProb(-99), chi2(-99), globChi2(-99),
      cat(-99), nValidTracker(-99), nValidPixel(-99), nMatches(-99),nHits(-99), nPixelHits(-99), globNHits(-99),validMuStations(-99),emEnergy(-99), hadEnergy(-99),nValidLayers(-99),isPF(-99),
      mcId(-99), mcMomId(-99), mcgMomId(-99){}
  public:
    TLorentzVector p4;
    int charge;
    float tIso, eIso, hIso, pfChaIso,pfChaPUIso,  pfPhoIso,pfNeuIso,acop, ipDb, ipErrDb, zPVPt,zPVProb, chi2, globChi2;
    int cat, nValidTracker, nValidPixel, nMatches,nHits, nPixelHits, globNHits, validMuStations;
    float emEnergy, hadEnergy;
    TLorentzVector mcFourMomentum;
    int mcId, mcMomId, mcgMomId;
    int isPF,nValidLayers;
    std::vector< int > hltMatchedBits;
  };

  class ElectronInfo {
  public:
    ElectronInfo() : scEta(-99), scPhi(-99), charge(-99), 
      tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99),pfChaPUIso(-99), pfPhoIso(-99), pfNeuIso(-99), acop(-99),
      sihih(-99), Dphi(-99), Deta(-99), HoE(-99), convDist(-99), convDcot(-99), innerHits(-99),
      isEB(false),isEE(false),
      ipDb(-99), ipErrDb(-99),
      id95(-99),id85(-99),id80(-99),id70(-99),
      id95r(-99),id85r(-99),id80r(-99),id70r(-99),
      mcId(-99), mcMomId(-99), mcgMomId (-99),
      dxy(-99), dz(-99),pfPhoIsoDoubleCounted(-99),
     //MVA ELEID 2012 input vars (keep same names to make life easier in integrating with MVA evaluation code
      fMVAVar_EoP(-99), fMVAVar_HoE(-99), fMVAVar_IoEmIoP(-99), fMVAVar_PreShowerOverRaw(-99), fMVAVar_R9(-99), fMVAVar_d0(-99), fMVAVar_deta(-99), fMVAVar_detacalo(-99), fMVAVar_dphi(-99), fMVAVar_e1x5e5x5(-99), fMVAVar_eleEoPout(-99), fMVAVar_eta(-99), fMVAVar_etawidth(-99), fMVAVar_fbrem(-99), fMVAVar_gsfchi2(-99), fMVAVar_ip3d(-99), fMVAVar_kfchi2(-99), fMVAVar_kfhits(-99), fMVAVar_phiwidth(-99), fMVAVar_pt(-99), fMVAVar_see(-99), fMVAVar_spp(-99), mvaOut(-99),mvaOutTrig(-99)
{}
  public:
    TLorentzVector p4;
    float scEta, scPhi;
    int charge;
    float tIso, eIso, hIso, pfChaIso,pfChaPUIso, pfPhoIso,pfNeuIso, acop;
    float sihih, Dphi, Deta, HoE, convDist, convDcot;
    int innerHits;
    bool isEB,isEE;
    float       ipDb, ipErrDb, dxy,dz;
    float pfPhoIsoDoubleCounted;
    float  id95,id85,id80,id70,id95r, id85r,id80r, id70r;
    float fMVAVar_EoP, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_PreShowerOverRaw, fMVAVar_R9, fMVAVar_d0, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_e1x5e5x5, fMVAVar_eleEoPout, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_ip3d, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_phiwidth, fMVAVar_pt, fMVAVar_see, fMVAVar_spp, mvaOut,mvaOutTrig;
    TLorentzVector mcFourMomentum;
    int mcId, mcMomId, mcgMomId;
    std::vector< int > hltMatchedBits;
  };

  class TauInfo{
  public:
    TauInfo()  : charge(-99), tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfChaPUIso(-99),pfPhoIso(-99), pfNeuIso(-99), acop(-99), 
      byIsolation(-99),trackIsolation(-99),byTaNCfrOnePercent(-99),
      byTaNCfrHalfPercent(-99), byTaNCfrQuarterPercent(-99), 
      byTaNCfrTenthPercent(-99), byTaNC(-99), mcId(-99), mcMomId(-99), mcgMomId(-99),
      againstElectronLoose(-99),  againstElectronMedium(-99),  againstElectronTight(-99),  againstMuonLoose(-99),  againstMuonTight(-99),  
      byLooseIsolation(-99),  byMediumIsolation(-99),  byTightIsolation(-99),  byVLooseIsolation(-99),  decayModeFinding(-99),
      isolationPFChargedHadrCandsPtSum(-999999.), isolationPFGammaCandsEtSum(-999999.),leadPFChargedHadrCandPt(-999999.),
      NsignalPFChargedHadrCands (-99),NsignalPFGammaCands(-99), byLooseCombinedIsolationDeltaBetaCorr(-99), againstElectronMVA(-99) {}
  public:
    TLorentzVector p4;
    int charge;
    float tIso, eIso, hIso,pfChaIso,pfChaPUIso,pfPhoIso,pfNeuIso, acop;
    float  byIsolation,trackIsolation,byTaNCfrOnePercent,byTaNCfrHalfPercent, byTaNCfrQuarterPercent, byTaNCfrTenthPercent, byTaNC;
    int mcId, mcMomId, mcgMomId;
    float againstElectronLoose, againstElectronMedium, againstElectronTight, againstMuonLoose, againstMuonTight;
    float byLooseIsolation, byMediumIsolation, byTightIsolation, byVLooseIsolation, decayModeFinding;

    float isolationPFChargedHadrCandsPtSum,isolationPFGammaCandsEtSum,leadPFChargedHadrCandPt;
    int NsignalPFChargedHadrCands, NsignalPFGammaCands;
    float byLooseCombinedIsolationDeltaBetaCorr, againstElectronMVA;


    TLorentzVector mcFourMomentum;
  };


  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> simpleJets3;
  std::vector<SimpleJet> simpleJets4; //???
  std::vector<SimpleJet> subJets; //???
  std::vector<SimpleJet> filterJets; //???
  std::vector<HardJet> hardJets;
//  std::vector<RawJet> rawJets;
//  std::vector<RawJet> rawJetsGen;    
  
    std::vector<HardJet> CAmdft12;
    std::vector<SimpleJet> CAmdft12_subjets;
    std::vector<SimpleJet> CApr12_subjets;
    std::vector<SimpleJet> CAft12_subjets;    
    std::vector<RawJet> CA12wConstits;
    std::vector<RawJet> CA12wConstitsGen;    
  
  METInfo calomet;
  METInfo tcmet;
  METInfo pfmet;
  METInfo pfmetType1corr;
  METInfo pfmetType1p2corr;
  METInfo pfmetNoPUType1corr;
  METInfo pfmetNoPUType1p2corr;
  METInfo mht;
  METInfo metNoPU;
  METInfo metCh;
  std::vector<METInfo> metUncInfo;

  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.35.2.3
Committed:	Tue Nov 27 05:04:18 2012 UTC (12 years, 5 months ago) by ntran
Content type:	text/plain
Branch:	hbbsubstructDevPostHCP
CVS Tags:	hbbsubstructDev_11, hbbsubstructDev_10, hbbsubstructDev_9, hbbsubstructDev_8, hbbsubstructDev_7, hbbsubstructDev_6, hbbsubstructDev_5
Changes since 1.35.2.2:	+1 -1 lines
Log Message:	major update, including new proposal for FJ3 inputs to step 2
#	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	degrutto	1.24	SimpleJet(): flavour(-99), isSemiLept(-99), isSemiLeptMCtruth(-99), SoftLeptpdgId(-99), SoftLeptIdlooseMu(-99), SoftLeptId95(-99), SoftLeptPt(-99), SoftLeptdR(-99), SoftLeptptRel(-99), SoftLeptRelCombIso(-99), tche(-99), tchp(-99), jpb(-99), jp(-99),
15	arizzi	1.28	ssvhe(-99), csv(-99), csvmva(-99), csvivf(-99), cmva(-99),
16	tboccali	1.11	vtxMass(-99), vtx3dL(-99), vtx3deL(-99), vtxNTracks(-99), vtxProb(-99),
17			ntracks(-99), charge(-99),jecunc(-99),
18	malbouis	1.31	chargedHadronEFraction(-99), neutralHadronEFraction(-99), chargedEmEFraction(-99), neutralEmEFraction(-99),nConstituents(-99), ptRaw(-99), ptLeadTrack(-99), jetArea(-99),
19	degrutto	1.34	bestMCid(-99), bestMCmomid(-99),
20			puJetIdL(-99), puJetIdM(-99), puJetIdT(-99), puJetIdMva(-99) {}
21	tboccali	1.1	public:
22	tboccali	1.6	double Pt() {return p4.Pt();}
23	degrutto	1.24	int flavour, isSemiLept, isSemiLeptMCtruth, SoftLeptpdgId , SoftLeptIdlooseMu, SoftLeptId95;
24			float SoftLeptPt, SoftLeptdR, SoftLeptptRel, SoftLeptRelCombIso;
25	arizzi	1.28	float tche,tchp, jpb,jp , ssvhe, csv, csvmva,csvivf,cmva;
26	tboccali	1.11	// CSV vertex info
27			float vtxMass, vtx3dL, vtx3deL, vtxNTracks, vtxProb;
28	arizzi	1.29	TVector3 vtxPosition;
29			TLorentzVector vtxP4;
30	malbouis	1.31	std::vector<unsigned int> vtxTrackIds, btagTrackIds, csvTrackIds;
31			int csvNTracks, btagNTracks;
32	tboccali	1.1	int ntracks;
33			float charge;
34	tboccali	1.6	TLorentzVector p4;
35	tboccali	1.1	TLorentzVector chargedTracksFourMomentum;
36
37	tboccali	1.11	// jec uncertainty
38
39			float jecunc;
40			// energy fractions as in
41			/* //chargedHadronEnergyFraction (relative to uncorrected jet energy)
42			float chargedHadronEnergyFraction() const {return chargedHadronEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
43			/// neutralHadronEnergyFraction (relative to uncorrected jet energy)
44			float neutralHadronEnergyFraction() const {return neutralHadronEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
45			/// chargedEmEnergyFraction (relative to uncorrected jet energy)
46			float chargedEmEnergyFraction() const {return chargedEmEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
47			/// neutralEmEnergyFraction (relative to uncorrected jet energy)
48			float neutralEmEnergyFraction() const {return neutralEmEnergy()/((jecSetsAvailable() ? jecFactor(0) : 1.)*energy());}
49			*/
50
51			float chargedHadronEFraction, neutralHadronEFraction, chargedEmEFraction, neutralEmEFraction;
52			int nConstituents;
53	malbouis	1.31	float ptRaw,ptLeadTrack;
54			float jetArea;
55	tboccali	1.11
56	bortigno	1.10	//MC parton matching
57	tboccali	1.11	TLorentzVector bestMCp4,bestMCp4mom;
58	tboccali	1.1	int bestMCid, bestMCmomid;
59			// new
60			TVector2 tVector;
61	bortigno	1.10	//Scale factor
62	bortigno	1.9	float SF_CSVL;
63			float SF_CSVM;
64			float SF_CSVT;
65			float SF_CSVLerr;
66			float SF_CSVMerr;
67			float SF_CSVTerr;
68	degrutto	1.34	// puJet Id
69			float puJetIdL, puJetIdM, puJetIdT, puJetIdMva;
70	bortigno	1.9
71	tboccali	1.1	};
72
73
74			class HardJet {
75			public:
76			HardJet(): constituents(-99){}
77			public:
78			int constituents;
79	tboccali	1.6	TLorentzVector p4;
80	tboccali	1.1	std::vector<TLorentzVector> subFourMomentum;
81			std::vector<float> etaSub, phiSub;
82			};
83
84	ntran	1.35.2.1	class RawJet {
85			public:
86			RawJet(): Nconstituents(-99){}
87			public:
88			int Nconstituents;
89			std::vector<float> constituents_px;
90			std::vector<float> constituents_py;
91			std::vector<float> constituents_pz;
92			std::vector<float> constituents_e;
93			std::vector<float> constituents_pdgId;
94			TLorentzVector p4;
95			// std::vector<TLorentzVector> subFourMomentum;
96			// std::vector<float> etaSub, phiSub;
97			};
98	tboccali	1.1
99			class METInfo {
100			public:
101			METInfo(): sumEt(-99), metSig(-99), eLong(-99){}
102			public:
103			float sumEt, metSig, eLong;
104	tboccali	1.6	TLorentzVector p4;
105	tboccali	1.1	};
106
107			class MuonInfo {
108			public:
109	degrutto	1.13	MuonInfo(): charge(-99),tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99),pfChaPUIso(99), pfPhoIso(-99), pfNeuIso(-99),
110	tboccali	1.1	acop(-99), ipDb(-99), ipErrDb(-99), zPVPt(-99),zPVProb(-99), chi2(-99), globChi2(-99),
111	arizzi	1.33	cat(-99), nValidTracker(-99), nValidPixel(-99), nMatches(-99),nHits(-99), nPixelHits(-99), globNHits(-99),validMuStations(-99),emEnergy(-99), hadEnergy(-99),nValidLayers(-99),isPF(-99),
112	tboccali	1.1	mcId(-99), mcMomId(-99), mcgMomId(-99){}
113			public:
114	tboccali	1.6	TLorentzVector p4;
115	tboccali	1.1	int charge;
116	degrutto	1.13	float tIso, eIso, hIso, pfChaIso,pfChaPUIso, pfPhoIso,pfNeuIso,acop, ipDb, ipErrDb, zPVPt,zPVProb, chi2, globChi2;
117	tboccali	1.16	int cat, nValidTracker, nValidPixel, nMatches,nHits, nPixelHits, globNHits, validMuStations;
118	tboccali	1.19	float emEnergy, hadEnergy;
119	tboccali	1.1	TLorentzVector mcFourMomentum;
120			int mcId, mcMomId, mcgMomId;
121	arizzi	1.33	int isPF,nValidLayers;
122	bortigno	1.8	std::vector< int > hltMatchedBits;
123	tboccali	1.1	};
124
125			class ElectronInfo {
126			public:
127			ElectronInfo() : scEta(-99), scPhi(-99), charge(-99),
128	degrutto	1.13	tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99),pfChaPUIso(-99), pfPhoIso(-99), pfNeuIso(-99), acop(-99),
129	tboccali	1.17	sihih(-99), Dphi(-99), Deta(-99), HoE(-99), convDist(-99), convDcot(-99), innerHits(-99),
130			isEB(false),isEE(false),
131	tboccali	1.18	ipDb(-99), ipErrDb(-99),
132	arizzi	1.7	id95(-99),id85(-99),id80(-99),id70(-99),
133			id95r(-99),id85r(-99),id80r(-99),id70r(-99),
134	arizzi	1.32	mcId(-99), mcMomId(-99), mcgMomId (-99),
135	arizzi	1.35	dxy(-99), dz(-99),pfPhoIsoDoubleCounted(-99),
136	arizzi	1.32	//MVA ELEID 2012 input vars (keep same names to make life easier in integrating with MVA evaluation code
137	arizzi	1.33	fMVAVar_EoP(-99), fMVAVar_HoE(-99), fMVAVar_IoEmIoP(-99), fMVAVar_PreShowerOverRaw(-99), fMVAVar_R9(-99), fMVAVar_d0(-99), fMVAVar_deta(-99), fMVAVar_detacalo(-99), fMVAVar_dphi(-99), fMVAVar_e1x5e5x5(-99), fMVAVar_eleEoPout(-99), fMVAVar_eta(-99), fMVAVar_etawidth(-99), fMVAVar_fbrem(-99), fMVAVar_gsfchi2(-99), fMVAVar_ip3d(-99), fMVAVar_kfchi2(-99), fMVAVar_kfhits(-99), fMVAVar_phiwidth(-99), fMVAVar_pt(-99), fMVAVar_see(-99), fMVAVar_spp(-99), mvaOut(-99),mvaOutTrig(-99)
138	arizzi	1.32	{}
139	tboccali	1.1	public:
140	tboccali	1.6	TLorentzVector p4;
141	tboccali	1.1	float scEta, scPhi;
142			int charge;
143	degrutto	1.13	float tIso, eIso, hIso, pfChaIso,pfChaPUIso, pfPhoIso,pfNeuIso, acop;
144	tboccali	1.17	float sihih, Dphi, Deta, HoE, convDist, convDcot;
145			int innerHits;
146			bool isEB,isEE;
147	arizzi	1.32	float ipDb, ipErrDb, dxy,dz;
148	arizzi	1.35	float pfPhoIsoDoubleCounted;
149	arizzi	1.7	float id95,id85,id80,id70,id95r, id85r,id80r, id70r;
150	arizzi	1.33	float fMVAVar_EoP, fMVAVar_HoE, fMVAVar_IoEmIoP, fMVAVar_PreShowerOverRaw, fMVAVar_R9, fMVAVar_d0, fMVAVar_deta, fMVAVar_detacalo, fMVAVar_dphi, fMVAVar_e1x5e5x5, fMVAVar_eleEoPout, fMVAVar_eta, fMVAVar_etawidth, fMVAVar_fbrem, fMVAVar_gsfchi2, fMVAVar_ip3d, fMVAVar_kfchi2, fMVAVar_kfhits, fMVAVar_phiwidth, fMVAVar_pt, fMVAVar_see, fMVAVar_spp, mvaOut,mvaOutTrig;
151	tboccali	1.1	TLorentzVector mcFourMomentum;
152			int mcId, mcMomId, mcgMomId;
153	bortigno	1.8	std::vector< int > hltMatchedBits;
154	tboccali	1.1	};
155
156			class TauInfo{
157			public:
158	sethzenz	1.27	TauInfo() : charge(-99), tIso(-99), eIso(-99), hIso(-99),pfChaIso(-99), pfChaPUIso(-99),pfPhoIso(-99), pfNeuIso(-99), acop(-99),
159	sethzenz	1.26	byIsolation(-99),trackIsolation(-99),byTaNCfrOnePercent(-99),
160			byTaNCfrHalfPercent(-99), byTaNCfrQuarterPercent(-99),
161			byTaNCfrTenthPercent(-99), byTaNC(-99), mcId(-99), mcMomId(-99), mcgMomId(-99),
162			againstElectronLoose(-99), againstElectronMedium(-99), againstElectronTight(-99), againstMuonLoose(-99), againstMuonTight(-99),
163			byLooseIsolation(-99), byMediumIsolation(-99), byTightIsolation(-99), byVLooseIsolation(-99), decayModeFinding(-99),
164	sethzenz	1.23	isolationPFChargedHadrCandsPtSum(-999999.), isolationPFGammaCandsEtSum(-999999.),leadPFChargedHadrCandPt(-999999.),
165	sethzenz	1.26	NsignalPFChargedHadrCands (-99),NsignalPFGammaCands(-99), byLooseCombinedIsolationDeltaBetaCorr(-99), againstElectronMVA(-99) {}
166	tboccali	1.1	public:
167	tboccali	1.6	TLorentzVector p4;
168	tboccali	1.1	int charge;
169	sethzenz	1.27	float tIso, eIso, hIso,pfChaIso,pfChaPUIso,pfPhoIso,pfNeuIso, acop;
170	sethzenz	1.26	float byIsolation,trackIsolation,byTaNCfrOnePercent,byTaNCfrHalfPercent, byTaNCfrQuarterPercent, byTaNCfrTenthPercent, byTaNC;
171	dlopes	1.22	int mcId, mcMomId, mcgMomId;
172	sethzenz	1.26	float againstElectronLoose, againstElectronMedium, againstElectronTight, againstMuonLoose, againstMuonTight;
173			float byLooseIsolation, byMediumIsolation, byTightIsolation, byVLooseIsolation, decayModeFinding;
174	dlopes	1.22
175			float isolationPFChargedHadrCandsPtSum,isolationPFGammaCandsEtSum,leadPFChargedHadrCandPt;
176	sethzenz	1.23	int NsignalPFChargedHadrCands, NsignalPFGammaCands;
177	sethzenz	1.26	float byLooseCombinedIsolationDeltaBetaCorr, againstElectronMVA;
178
179
180	tboccali	1.1	TLorentzVector mcFourMomentum;
181			};
182
183
184			class DiMuonInfo {
185			public:
186	tboccali	1.6	TLorentzVector p4;
187	tboccali	1.1	MuonInfo daughter1, daughter2;
188			};
189
190
191			class DiElectronInfo {
192			public:
193	tboccali	1.6	TLorentzVector p4;
194	tboccali	1.1	ElectronInfo daughter1, daughter2;
195			};
196
197
198
199			public:
200			std::vector<SimpleJet> simpleJets;
201			std::vector<SimpleJet> simpleJets2; //???
202	tboccali	1.12	std::vector<SimpleJet> simpleJets3;
203			std::vector<SimpleJet> simpleJets4; //???
204	tboccali	1.1	std::vector<SimpleJet> subJets; //???
205	dlopes	1.22	std::vector<SimpleJet> filterJets; //???
206	tboccali	1.1	std::vector<HardJet> hardJets;
207	ntran	1.35.2.2	// std::vector<RawJet> rawJets;
208			// std::vector<RawJet> rawJetsGen;
209
210	ntran	1.35.2.3	std::vector<HardJet> CAmdft12;
211	ntran	1.35.2.2	std::vector<SimpleJet> CAmdft12_subjets;
212			std::vector<SimpleJet> CApr12_subjets;
213			std::vector<SimpleJet> CAft12_subjets;
214			std::vector<RawJet> CA12wConstits;
215			std::vector<RawJet> CA12wConstitsGen;
216	tboccali	1.1
217	degrutto	1.25	METInfo calomet;
218	tboccali	1.1	METInfo tcmet;
219			METInfo pfmet;
220	degrutto	1.25	METInfo pfmetType1corr;
221			METInfo pfmetType1p2corr;
222			METInfo pfmetNoPUType1corr;
223			METInfo pfmetNoPUType1p2corr;
224	tboccali	1.11	METInfo mht;
225	arizzi	1.14	METInfo metNoPU;
226	tboccali	1.20	METInfo metCh;
227	degrutto	1.25	std::vector<METInfo> metUncInfo;
228	tboccali	1.11
229	tboccali	1.1	std::vector<MuonInfo> muInfo;
230			std::vector<ElectronInfo> eleInfo;
231			std::vector<TauInfo> tauInfo;
232
233			std::vector<DiMuonInfo> diMuonInfo;
234			std::vector<DiElectronInfo> diElectronInfo;
235
236			// HiggsCandidate dijetHiggs;
237			// HiggsCandidate fatjetHiggs;
238
239			// to be decided.. one single or many V?
240			// VCandidate vcand;
241			// WCandidate wele;
242			// WCandidate wmu;
243			//WCandidate wtau;
244			// ZCandidate zmu;
245			// ZCandidate zele;
246			//ZInvCandidate znunu;
247
248			//
249			// really needed????
250			//
251
252			};
253			#endif