HbbAnalysis/interface/Objects.hh

#ifndef HbbAnalysis_Objects_hh
#define HbbAnalysis_Objects_hh

#include <vector>
#include <string>
#include <iostream>

#include "TLorentzVector.h"
#include "TMath.h"

namespace HbbAnalysis {

  struct PUVars {
    int bunchCrossing;
    unsigned int numInteractions;
    PUVars(){
      bunchCrossing = 0;
      numInteractions = 0;
    }
  };

  struct MCVars {
    unsigned int index;
    double E;
    double pT;
    double eta;
    double y;
    double phi;
    int pdgId;
    int status;

    MCVars():
      index(0),
      E(0),
      pT(0),
      eta(0),
      y(0),
      phi(0),
      pdgId(0),
      status(0)
    {;}
  };

  struct HLTVars {
    float pT;
    float eta;
    float phi;
    int id;
    float mass;
    unsigned int type;//1=EM, 2=muon, 3=jet, 0=unknown
    HLTVars():
      pT(0),
      eta(0),
      phi(0),
      id(0),
      mass(0),
      type(0)
    {}
  };

  struct L1Vars {
    //type:
    //0 = empty
    //1 = central jet
    //2 = tau jet
    //3 = forward jet
    unsigned int type;
    int bx;
    float pT;
    float ET;
    float eta;
    float phi; 
    L1Vars():
      type(0),
      bx(0),
      pT(0),
      ET(0),
      eta(0),
      phi(0)
    {}

  };

  struct GenVars {
    bool valid;
    double E;
    double pT;
    double eta;
    double y;
    double phi;
    int charge;
    int pdgId;
    int status;
    double mass;
    double vx;
    double vy;
    double vz;
    GenVars():
      valid(false),
      E(0),
      pT(0),
      eta(0),
      y(0),
      phi(0),
      charge(0),
      pdgId(0),
      status(0),
      mass(0),
      vx(0),
      vy(0),
      vz(0)
    {}

  };

  struct BaseVars {
    double E;
    double pT;
    double eta;
    double y;
    double phi;
    float charge;
    double vx;
    double vy;
    double vz;
    bool hltMatch;

    BaseVars():
      E(0),
      pT(0),
      eta(0),
      y(0),
      phi(0),
      charge(0),
      vx(0),
      vy(0),
      vz(0),
      hltMatch(true)
    {}
  };

  struct SCVars {
    float sigmaEtaEta;
    float sigmaIEtaIEta;
    float e1x5;
    float e2x5Max;
    float e5x5;
    float eOverP;
    SCVars():
      sigmaEtaEta(0),
      sigmaIEtaIEta(0),
      e1x5(0),
      e2x5Max(0),
      e5x5(0),
      eOverP(0)
    {}
  };

  struct EleIsoVars {
    float calo;
    float track;
    float ecal;
    float hcal;
    EleIsoVars():
      calo(0),
      track(0),
      ecal(0),
      hcal(0)
    {}
  };
  
  struct EleIDVars{
    unsigned short idAndIso;
    std::vector<std::pair<std::string,float> > electronIDs;
    float hOverE;
    float deltaPhiIn;
    float deltaEtaIn;
    bool ecalDrivenSeed;
    bool trackerDrivenSeed;
    float dB;
    EleIDVars():
      idAndIso(0),
      hOverE(0),
      deltaPhiIn(0),
      deltaEtaIn(0),
      ecalDrivenSeed(false),
      trackerDrivenSeed(false),
      dB(0)
    {
      electronIDs.clear();
    }
  };

  struct MuTrkVars {
    float dxy;
    float dz;
    float normalizedChi2;
    unsigned int muonHits;
    int charge;
    unsigned int trackerHits;
    unsigned int pixelHits;
    MuTrkVars():
      dxy(0),
      dz(0),
      normalizedChi2(0),
      muonHits(0),
      charge(0),
      trackerHits(0),
      pixelHits(0)
    {}
  };
  
  struct MuIsoVars {
    float sumPt;
    float emEt;
    float hadEt;
    float nTracks;
    float nJets;
    MuIsoVars():
      sumPt(0),
      emEt(0),
      hadEt(0),
      nTracks(0),
      nJets(0)
    {}
  };
  
  struct MuIDVars{
    unsigned short type;
    std::vector<unsigned short> ids;
    float caloCompat;
    float segCompat;
    unsigned int nChambers;
    unsigned int nMatchesLoose;
    unsigned int nMatchesMedium;
    unsigned int nMatchesTight;
    float dB;
    MuIDVars():
      type(0),
      caloCompat(0),
      segCompat(0),
      nChambers(0),
      nMatchesLoose(0),
      nMatchesMedium(0),
      nMatchesTight(0),
      dB(0)
    {
      ids.clear();
    }


  };

  
  struct TrkVars{
    double pT;
    double eta;
    double phi;
    double matchDist;
    double IPxy;
    double IPz;
    float signedSipt;
    TrkVars():
      pT(0),
      eta(0),
      phi(0),
      matchDist(0),
      IPxy(0),
      IPz(0),
      signedSipt(0)
    {}

  };
  

  struct CaloTauIsoVars{
    unsigned int nIsoTracks;
    unsigned int nSigTracks;
    float leadTrackHCAL3x3hitsEtSum;
    float leadTrackHCAL3x3hottesthitDEta;
    float signalTracksInvariantMass;
    float tracksInvariantMass; 
    float isolationTracksPtSum;
    float isolationECALhitsEtSum;
    float maximumHCALhitEt;
    CaloTauIsoVars():
      nIsoTracks(0),
      nSigTracks(0),
      leadTrackHCAL3x3hitsEtSum(0),
      leadTrackHCAL3x3hottesthitDEta(0),
      signalTracksInvariantMass(0),
      tracksInvariantMass(0), 
      isolationTracksPtSum(0),
      isolationECALhitsEtSum(0),
      maximumHCALhitEt(0)
    {}


  };

  struct CaloTauIDVars{
    float byIsolation;
    float byLeadingTrackFinding;
    float byLeadingTrackPtCut;
    CaloTauIDVars():
      byIsolation(0),
      byLeadingTrackFinding(0),
      byLeadingTrackPtCut(0)
    {}

  };

  struct PFTauIsoVars{
    unsigned int nSigCands;
    unsigned int nIsoCands;
    double maximumHCALPFClusterEt;
    //double emFraction;
    double hcalTotOverPLead;
    double hcalMaxOverPLead;
    double hcal3x3OverPLead;
    double ecalStripSumEOverPLead;
    double bremsRecoveryEOverPLead;
    PFTauIsoVars():
      nSigCands(0),
      nIsoCands(0),
      maximumHCALPFClusterEt(0),
                //emFraction(0),
      hcalTotOverPLead(0),
      hcalMaxOverPLead(0),
      hcal3x3OverPLead(0),
      ecalStripSumEOverPLead(0),
      bremsRecoveryEOverPLead(0)
    {}
  };

  //for hadr, neutr and gamma cands
  struct PFTauCandVars{
    unsigned int nSigCands;
    unsigned int nIsoCands;
    double isolationPtSum;
    PFTauCandVars():
      nSigCands(0),
      nIsoCands(0),
      isolationPtSum(0)
    {}
  };

  struct PFTauIDVars{
    float byLeadingTrackFinding;
    float byLeadingTrackPtCut;
    float byTrackIsolation;
    float byECALIsolation;
    float byIsolation;
    float againstElectron;
    float againstMuon;
    float byIsolationUsingLeadingPion;
    float byTaNC;
    float byTaNCfrHalfPercent;
    float byTaNCfrOnePercent;
    float byTaNCfrQuarterPercent;
    float byTaNCfrTenthPercent;
    float ecalIsolationUsingLeadingPion;
    float leadingPionPtCut;
    float trackIsolationUsingLeadingPion;
    PFTauIDVars():
      byLeadingTrackFinding(0),
      byLeadingTrackPtCut(0),
      byTrackIsolation(0),
      byECALIsolation(0),
      byIsolation(0),
      againstElectron(0),
      againstMuon(0),
      byIsolationUsingLeadingPion(0),
      byTaNC(0),
      byTaNCfrHalfPercent(0),
      byTaNCfrOnePercent(0),
      byTaNCfrQuarterPercent(0),
      byTaNCfrTenthPercent(0),
      ecalIsolationUsingLeadingPion(0),
      leadingPionPtCut(0),
      trackIsolationUsingLeadingPion(0)
    {}
  };

  struct PFTauEleIDVars{
    double pT;
    double eta;
    double phi;
    float output;
    float decision;

    PFTauEleIDVars():
      pT(0),
      eta(0),
      phi(0),
      output(0),
      decision(0)
    {}
  };
    
  struct PFTauMuIDVars{
    float caloCompat;
    float segCompat;
    float decision;
    PFTauMuIDVars():
      caloCompat(0),
      segCompat(0),
      decision(0)
    {}
  };

  struct JetVars{
    //1=had, 2=e, 3=mu for heavy quarks; 
    //2=uds, 3=g for light quarks
    unsigned int flavour;
    int partonFlavour;
    unsigned int nAssociatedTracks;
    double rawpT;
    float etaMean;
    float phiMean;
    float etaEtaMoment;
    float phiPhiMoment;
    float etaPhiMoment;
    //double rawEta;
    //double rawPhi;
    JetVars():
      flavour(0),
      partonFlavour(0),
      nAssociatedTracks(0),
      rawpT(0),
      etaMean(0),
      phiMean(0),
      etaEtaMoment(0),
      phiPhiMoment(0),
      etaPhiMoment(0)
    {}
  };

  struct CaloJetVars{
    double maxEInEmTowers;
    double maxEInHadTowers;
    double energyFractionHadronic;
    double emEnergyFraction;
    double hadEnergyInHB;
    double hadEnergyInHO;
    double hadEnergyInHE;
    double hadEnergyInHF;
    double emEnergyInEB;
    double emEnergyInEE;
    double emEnergyInHF;
    double towersArea;
    double n90;
    double n60;
    CaloJetVars():
      maxEInEmTowers(0),
      maxEInHadTowers(0),
      energyFractionHadronic(0),
      emEnergyFraction(0),
      hadEnergyInHB(0),
      hadEnergyInHO(0),
      hadEnergyInHE(0),
      hadEnergyInHF(0),
      emEnergyInEB(0),
      emEnergyInEE(0),
      emEnergyInHF(0),
      towersArea(0),
      n90(0),
      n60(0)
    {}
  };

  struct JPTJetVars{
    float zspCorrection;
    int elecMultiplicity;
    float calopT;
    float caloEta;
    //std::vector<TrkVars> pionsInCone;
    //std::vector<TrkVars> pionsCurledOut;
    //std::vector<TrkVars> pionsCurledIn;
    //std::vector<TrkVars> elecsInCone;
    //std::vector<TrkVars> elecsCurledOut;
    //std::vector<TrkVars> elecsCurledIn;
    //std::vector<TrkVars> muonsInCone;
    //std::vector<TrkVars> muonsCurledOut;
    //std::vector<TrkVars> muonsCurledIn;
    std::vector< std::pair<unsigned int,double> > particleStatusPt;
    JPTJetVars():
      zspCorrection(0),
      elecMultiplicity(0),
      calopT(0),
      caloEta(0)
    {
      particleStatusPt.clear();
    }
  };

  struct JPTPFJetVars{
    double chargedHadronEnergy;
    double chargedHadronEnergyFraction;
    double neutralHadronEnergy;
    double neutralHadronEnergyFraction;
    double chargedEmEnergy;
    double chargedEmEnergyFraction;
    double neutralEmEnergy;
    double neutralEmEnergyFraction;
    double chargedMultiplicity;
    double muonMultiplicity;
    JPTPFJetVars():
      chargedHadronEnergy(0),
      chargedHadronEnergyFraction(0),
      neutralHadronEnergy(0),
      neutralHadronEnergyFraction(0),
      chargedEmEnergy(0),
      chargedEmEnergyFraction(0),
      neutralEmEnergy(0),
      neutralEmEnergyFraction(0),
      chargedMultiplicity(0),
      muonMultiplicity(0)
    {}
  };

  struct PFJetVars{
    double chargedMuEnergy;
    double chargedMuEnergyFraction;
    double neutralMultiplicity;
    unsigned numberOfDaughters;
    double HFHadronEnergy;
    PFJetVars():
      chargedMuEnergy(0),
      chargedMuEnergyFraction(0),
      neutralMultiplicity(0),
      numberOfDaughters(0),
      HFHadronEnergy(0)
    {}
  };

  struct JetBtagVars{
    double cSV;
    double cSVMVA;
    double iPMVA;
    double bProba;
    double probability;
    double sSVHE;
    double sSVHP;
    double softElectronByPt;
    double softElectronByIP3d;
    double softMuon;
    double softMuonByPt;
    double softMuonByIP3d;
    double tCHE;
    double tCHP;
    JetBtagVars():
      cSV(0),
      cSVMVA(0),
      iPMVA(0),
      bProba(0),
      probability(0),
      sSVHE(0),
      sSVHP(0),
      softElectronByPt(0),
      softElectronByIP3d(0),
      softMuon(0),
      softMuonByPt(0),
      softMuonByIP3d(0),
      tCHE(0),
      tCHP(0)
    {}
  };


  struct JetIDVars{
    double fHPD; 
    double fRBX;
    int n90Hits;
    //double fSubDetector1;
    //double fSubDetector2;
    //double fSubDetector3;
    //double fSubDetector4;
    double restrictedEMF;
    int    nHCALTowers;
    int    nECALTowers;
    //double approximatefHPD;
    //double approximatefRBX;
    int hitsInN90;
    // muon hits id
    //int numberOfHits2RPC;
    //int numberOfHits3RPC;
    //int numberOfHitsRPC;
    JetIDVars():
      fHPD(0), 
      fRBX(0),
      n90Hits(0),
      restrictedEMF(0),
      nHCALTowers(0),
      nECALTowers(0),
      hitsInN90(0)
    {}
  };

  struct JetECorVars{
    std::vector<std::pair<std::string,double> > Levels;
    JetECorVars(){
      Levels.clear();
    }
  };


  struct MetVars{
    double mET;
    double mEx;
    double mEy;
    double sumET;
    double phi;
    double mEtSig;
    MetVars():
      mET(0),
      mEx(0),
      mEy(0),
      sumET(0),
      phi(0),
      mEtSig(0)
    {}

  };

  struct TriggerVars{
    std::string name;
    unsigned int index;
    bool accept;
    TriggerVars():
      name(""),
      index(0),
      accept(false)
    {}
  };


  struct VertexVars{
    std::vector<float> trackWeights;
    std::vector<float> trackpT;
    double chi2;
    double ndof;
    double x;
    double y;
    double z;
    double rho;
    double xError;
    double yError;
    double zError;
    double cov01;
    double cov02;
    double cov12;
    VertexVars():
      chi2(0),
      ndof(0),
      x(0),
      y(0),
      z(0),
      rho(0),
      xError(0),
      yError(0),
      zError(0),
      cov01(0),
      cov02(0),
      cov12(0)
    {
      trackWeights.clear();
      trackpT.clear();
    }
  };

  struct BeamSpotVars{
    double x0;
    double y0;
    double z0;
    double sigmaZ;
    double BeamWidthX;
    double BeamWidthY;
    double x0Error;
    double y0Error;
    double z0Error;
    double sigmaZ0Error;
    double BeamWidthXError;
    double BeamWidthYError;     
    //double dxdz;
    //double dydz;
    //double dxdzError;
    //double dydzError;
    //BeamType type();
    //enum BeamType { Unknown=-1, Fake=0, LHC=1, Tracker=2 };
    //double emittanceX() const { return emittanceX_; }
    //double emittanceY() const { return emittanceY_; }
    //double betaStar() const { return betaStar_; }
    BeamSpotVars():
      x0(0),
      y0(0),
      z0(0),
      sigmaZ(0),
      BeamWidthX(0),
      BeamWidthY(0),
      x0Error(0),
      y0Error(0),
      z0Error(0),
      sigmaZ0Error(0),
      BeamWidthXError(0),
      BeamWidthYError(0)
    {}

  };

  double DeltaPhi(const double phi1, const double phi2);
 
  template <class T1, class T2>
  double DeltaR(const T1 & v1, const T2 & v2){
    double dEta = v1.eta - v2.eta;
    double dPhi = fabs(v1.phi - v2.phi);
    if (dPhi > TMath::Pi()) dPhi = (2.0*TMath::Pi() - dPhi);
    return sqrt(dEta*dEta+dPhi*dPhi);
  }

  double SameSign(const BaseVars & v1, const BaseVars & v2);
 
  double OppSign(const BaseVars & v1, const BaseVars & v2);
 
  TLorentzVector FourMomentum(const BaseVars & v, const double scale=1) ;
 
  double TransverseMass(const BaseVars & leg1, 
                        const BaseVars & leg2, 
                        const double mEx, 
                        const double mEy);
 
  double TransverseMass(const BaseVars & leg1, 
                        const double mEx, 
                        const double mEy);

  TLorentzVector FourMomentumCDFmethod(const BaseVars & leg1, 
                                       const BaseVars & leg2, 
                                       double mEx, 
                                       double mEy);

  TLorentzVector FourMomentumCollinearApprox(const BaseVars & leg1,
                                             const BaseVars & leg2,
                                             double mEx, 
                                             double mEy);

  double EtaDetector(const BaseVars & v1);
  double EtaDetector(const GenVars & v1);

}//namespace

#endif
Revision:	1.26
Committed:	Sun Aug 7 09:51:36 2011 UTC (13 years, 8 months ago) by agilbert
Content type:	text/plain
Branch:	MAIN
Changes since 1.25:	+3 -1 lines
Log Message:	Change DR in FilterZbb
#	User	Rev	Content
1	amagnan	1.1	#ifndef HbbAnalysis_Objects_hh
2			#define HbbAnalysis_Objects_hh
3
4			#include <vector>
5			#include <string>
6	amagnan	1.2	#include <iostream>
7
8			#include "TLorentzVector.h"
9	agilbert	1.26	#include "TMath.h"
10	amagnan	1.1
11			namespace HbbAnalysis {
12	amagnan	1.12
13	agilbert	1.23	struct PUVars {
14			int bunchCrossing;
15			unsigned int numInteractions;
16			PUVars(){
17			bunchCrossing = 0;
18			numInteractions = 0;
19			}
20			};
21
22	amagnan	1.1	struct MCVars {
23			unsigned int index;
24	amagnan	1.2	double E;
25	amagnan	1.1	double pT;
26			double eta;
27	amagnan	1.9	double y;
28	amagnan	1.1	double phi;
29			int pdgId;
30			int status;
31	amagnan	1.24
32			MCVars():
33			index(0),
34			E(0),
35			pT(0),
36			eta(0),
37			y(0),
38			phi(0),
39			pdgId(0),
40			status(0)
41			{;}
42	amagnan	1.1	};
43
44	amagnan	1.11	struct HLTVars {
45			float pT;
46			float eta;
47			float phi;
48			int id;
49			float mass;
50	amagnan	1.12	unsigned int type;//1=EM, 2=muon, 3=jet, 0=unknown
51	amagnan	1.24	HLTVars():
52			pT(0),
53			eta(0),
54			phi(0),
55			id(0),
56			mass(0),
57			type(0)
58			{}
59	amagnan	1.11	};
60
61			struct L1Vars {
62			//type:
63			//0 = empty
64			//1 = central jet
65			//2 = tau jet
66			//3 = forward jet
67			unsigned int type;
68			int bx;
69			float pT;
70			float ET;
71			float eta;
72			float phi;
73	amagnan	1.24	L1Vars():
74			type(0),
75			bx(0),
76			pT(0),
77			ET(0),
78			eta(0),
79			phi(0)
80			{}
81
82	amagnan	1.11	};
83
84	amagnan	1.1	struct GenVars {
85			bool valid;
86	amagnan	1.2	double E;
87	amagnan	1.1	double pT;
88			double eta;
89	amagnan	1.9	double y;
90	amagnan	1.1	double phi;
91			int charge;
92			int pdgId;
93			int status;
94			double mass;
95			double vx;
96			double vy;
97			double vz;
98	amagnan	1.24	GenVars():
99			valid(false),
100			E(0),
101			pT(0),
102			eta(0),
103			y(0),
104			phi(0),
105			charge(0),
106			pdgId(0),
107			status(0),
108			mass(0),
109			vx(0),
110			vy(0),
111			vz(0)
112			{}
113
114	amagnan	1.1	};
115
116			struct BaseVars {
117	amagnan	1.2	double E;
118	amagnan	1.1	double pT;
119			double eta;
120	amagnan	1.9	double y;
121	amagnan	1.1	double phi;
122			float charge;
123			double vx;
124			double vy;
125			double vz;
126	agilbert	1.22	bool hltMatch;
127	amagnan	1.24
128			BaseVars():
129			E(0),
130			pT(0),
131			eta(0),
132			y(0),
133			phi(0),
134			charge(0),
135			vx(0),
136			vy(0),
137			vz(0),
138			hltMatch(true)
139			{}
140	amagnan	1.1	};
141
142			struct SCVars {
143			float sigmaEtaEta;
144			float sigmaIEtaIEta;
145			float e1x5;
146			float e2x5Max;
147			float e5x5;
148			float eOverP;
149	amagnan	1.24	SCVars():
150			sigmaEtaEta(0),
151			sigmaIEtaIEta(0),
152			e1x5(0),
153			e2x5Max(0),
154			e5x5(0),
155			eOverP(0)
156			{}
157	amagnan	1.1	};
158
159			struct EleIsoVars {
160			float calo;
161			float track;
162			float ecal;
163			float hcal;
164	amagnan	1.24	EleIsoVars():
165			calo(0),
166			track(0),
167			ecal(0),
168			hcal(0)
169			{}
170	amagnan	1.1	};
171
172			struct EleIDVars{
173	amagnan	1.17	unsigned short idAndIso;
174	amagnan	1.1	std::vector<std::pair<std::string,float> > electronIDs;
175			float hOverE;
176			float deltaPhiIn;
177			float deltaEtaIn;
178	amagnan	1.7	bool ecalDrivenSeed;
179			bool trackerDrivenSeed;
180	agilbert	1.22	float dB;
181	amagnan	1.24	EleIDVars():
182			idAndIso(0),
183			hOverE(0),
184			deltaPhiIn(0),
185			deltaEtaIn(0),
186			ecalDrivenSeed(false),
187			trackerDrivenSeed(false),
188			dB(0)
189			{
190			electronIDs.clear();
191			}
192	amagnan	1.1	};
193
194	amagnan	1.3	struct MuTrkVars {
195	amagnan	1.13	float dxy;
196			float dz;
197			float normalizedChi2;
198			unsigned int muonHits;
199	amagnan	1.14	int charge;
200	amagnan	1.13	unsigned int trackerHits;
201	amagnan	1.15	unsigned int pixelHits;
202	amagnan	1.24	MuTrkVars():
203			dxy(0),
204			dz(0),
205			normalizedChi2(0),
206			muonHits(0),
207			charge(0),
208			trackerHits(0),
209			pixelHits(0)
210			{}
211	amagnan	1.3	};
212
213	amagnan	1.1	struct MuIsoVars {
214			float sumPt;
215			float emEt;
216			float hadEt;
217			float nTracks;
218			float nJets;
219	amagnan	1.24	MuIsoVars():
220			sumPt(0),
221			emEt(0),
222			hadEt(0),
223			nTracks(0),
224			nJets(0)
225			{}
226	amagnan	1.1	};
227
228			struct MuIDVars{
229			unsigned short type;
230			std::vector<unsigned short> ids;
231			float caloCompat;
232			float segCompat;
233			unsigned int nChambers;
234			unsigned int nMatchesLoose;
235			unsigned int nMatchesMedium;
236			unsigned int nMatchesTight;
237	agilbert	1.22	float dB;
238	amagnan	1.24	MuIDVars():
239			type(0),
240			caloCompat(0),
241			segCompat(0),
242			nChambers(0),
243			nMatchesLoose(0),
244			nMatchesMedium(0),
245			nMatchesTight(0),
246			dB(0)
247			{
248			ids.clear();
249			}
250
251
252	amagnan	1.1	};
253
254	agilbert	1.20
255	amagnan	1.9	struct TrkVars{
256	amagnan	1.1	double pT;
257			double eta;
258			double phi;
259			double matchDist;
260			double IPxy;
261			double IPz;
262			float signedSipt;
263	amagnan	1.24	TrkVars():
264			pT(0),
265			eta(0),
266			phi(0),
267			matchDist(0),
268			IPxy(0),
269			IPz(0),
270			signedSipt(0)
271			{}
272
273	amagnan	1.1	};
274	agilbert	1.20
275	amagnan	1.1
276			struct CaloTauIsoVars{
277			unsigned int nIsoTracks;
278			unsigned int nSigTracks;
279			float leadTrackHCAL3x3hitsEtSum;
280			float leadTrackHCAL3x3hottesthitDEta;
281			float signalTracksInvariantMass;
282			float tracksInvariantMass;
283			float isolationTracksPtSum;
284			float isolationECALhitsEtSum;
285			float maximumHCALhitEt;
286	amagnan	1.24	CaloTauIsoVars():
287			nIsoTracks(0),
288			nSigTracks(0),
289			leadTrackHCAL3x3hitsEtSum(0),
290			leadTrackHCAL3x3hottesthitDEta(0),
291			signalTracksInvariantMass(0),
292			tracksInvariantMass(0),
293			isolationTracksPtSum(0),
294			isolationECALhitsEtSum(0),
295			maximumHCALhitEt(0)
296			{}
297
298
299
300	amagnan	1.1	};
301
302			struct CaloTauIDVars{
303			float byIsolation;
304			float byLeadingTrackFinding;
305			float byLeadingTrackPtCut;
306	amagnan	1.24	CaloTauIDVars():
307			byIsolation(0),
308			byLeadingTrackFinding(0),
309			byLeadingTrackPtCut(0)
310			{}
311
312	amagnan	1.1	};
313
314			struct PFTauIsoVars{
315			unsigned int nSigCands;
316			unsigned int nIsoCands;
317			double maximumHCALPFClusterEt;
318			//double emFraction;
319			double hcalTotOverPLead;
320			double hcalMaxOverPLead;
321			double hcal3x3OverPLead;
322			double ecalStripSumEOverPLead;
323			double bremsRecoveryEOverPLead;
324	amagnan	1.24	PFTauIsoVars():
325			nSigCands(0),
326			nIsoCands(0),
327			maximumHCALPFClusterEt(0),
328			//emFraction(0),
329			hcalTotOverPLead(0),
330			hcalMaxOverPLead(0),
331			hcal3x3OverPLead(0),
332			ecalStripSumEOverPLead(0),
333			bremsRecoveryEOverPLead(0)
334			{}
335	amagnan	1.1	};
336
337			//for hadr, neutr and gamma cands
338			struct PFTauCandVars{
339			unsigned int nSigCands;
340			unsigned int nIsoCands;
341			double isolationPtSum;
342	amagnan	1.24	PFTauCandVars():
343			nSigCands(0),
344			nIsoCands(0),
345			isolationPtSum(0)
346			{}
347	amagnan	1.1	};
348
349			struct PFTauIDVars{
350			float byLeadingTrackFinding;
351			float byLeadingTrackPtCut;
352			float byTrackIsolation;
353			float byECALIsolation;
354			float byIsolation;
355			float againstElectron;
356			float againstMuon;
357	amagnan	1.4	float byIsolationUsingLeadingPion;
358			float byTaNC;
359			float byTaNCfrHalfPercent;
360			float byTaNCfrOnePercent;
361			float byTaNCfrQuarterPercent;
362			float byTaNCfrTenthPercent;
363			float ecalIsolationUsingLeadingPion;
364			float leadingPionPtCut;
365			float trackIsolationUsingLeadingPion;
366	amagnan	1.24	PFTauIDVars():
367			byLeadingTrackFinding(0),
368			byLeadingTrackPtCut(0),
369			byTrackIsolation(0),
370			byECALIsolation(0),
371			byIsolation(0),
372			againstElectron(0),
373			againstMuon(0),
374			byIsolationUsingLeadingPion(0),
375			byTaNC(0),
376			byTaNCfrHalfPercent(0),
377			byTaNCfrOnePercent(0),
378			byTaNCfrQuarterPercent(0),
379			byTaNCfrTenthPercent(0),
380			ecalIsolationUsingLeadingPion(0),
381			leadingPionPtCut(0),
382			trackIsolationUsingLeadingPion(0)
383			{}
384	amagnan	1.1	};
385
386			struct PFTauEleIDVars{
387			double pT;
388			double eta;
389			double phi;
390			float output;
391			float decision;
392	amagnan	1.24
393			PFTauEleIDVars():
394			pT(0),
395			eta(0),
396			phi(0),
397			output(0),
398			decision(0)
399			{}
400	amagnan	1.1	};
401
402			struct PFTauMuIDVars{
403			float caloCompat;
404			float segCompat;
405			float decision;
406	amagnan	1.24	PFTauMuIDVars():
407			caloCompat(0),
408			segCompat(0),
409			decision(0)
410			{}
411	amagnan	1.1	};
412
413			struct JetVars{
414			//1=had, 2=e, 3=mu for heavy quarks;
415			//2=uds, 3=g for light quarks
416			unsigned int flavour;
417			int partonFlavour;
418			unsigned int nAssociatedTracks;
419			double rawpT;
420	amagnan	1.5	float etaMean;
421			float phiMean;
422			float etaEtaMoment;
423			float phiPhiMoment;
424			float etaPhiMoment;
425	amagnan	1.1	//double rawEta;
426			//double rawPhi;
427	amagnan	1.24	JetVars():
428			flavour(0),
429			partonFlavour(0),
430			nAssociatedTracks(0),
431			rawpT(0),
432			etaMean(0),
433			phiMean(0),
434			etaEtaMoment(0),
435			phiPhiMoment(0),
436			etaPhiMoment(0)
437			{}
438	amagnan	1.1	};
439
440			struct CaloJetVars{
441			double maxEInEmTowers;
442			double maxEInHadTowers;
443			double energyFractionHadronic;
444			double emEnergyFraction;
445			double hadEnergyInHB;
446			double hadEnergyInHO;
447			double hadEnergyInHE;
448			double hadEnergyInHF;
449			double emEnergyInEB;
450			double emEnergyInEE;
451			double emEnergyInHF;
452			double towersArea;
453			double n90;
454			double n60;
455	amagnan	1.24	CaloJetVars():
456			maxEInEmTowers(0),
457			maxEInHadTowers(0),
458			energyFractionHadronic(0),
459			emEnergyFraction(0),
460			hadEnergyInHB(0),
461			hadEnergyInHO(0),
462			hadEnergyInHE(0),
463			hadEnergyInHF(0),
464			emEnergyInEB(0),
465			emEnergyInEE(0),
466			emEnergyInHF(0),
467			towersArea(0),
468			n90(0),
469			n60(0)
470			{}
471	amagnan	1.1	};
472
473	amagnan	1.9	struct JPTJetVars{
474	amagnan	1.10	float zspCorrection;
475			int elecMultiplicity;
476	agilbert	1.20	float calopT;
477			float caloEta;
478			//std::vector<TrkVars> pionsInCone;
479			//std::vector<TrkVars> pionsCurledOut;
480			//std::vector<TrkVars> pionsCurledIn;
481			//std::vector<TrkVars> elecsInCone;
482			//std::vector<TrkVars> elecsCurledOut;
483			//std::vector<TrkVars> elecsCurledIn;
484			//std::vector<TrkVars> muonsInCone;
485			//std::vector<TrkVars> muonsCurledOut;
486			//std::vector<TrkVars> muonsCurledIn;
487			std::vector< std::pair<unsigned int,double> > particleStatusPt;
488	amagnan	1.24	JPTJetVars():
489			zspCorrection(0),
490			elecMultiplicity(0),
491			calopT(0),
492			caloEta(0)
493			{
494			particleStatusPt.clear();
495			}
496	amagnan	1.9	};
497
498			struct JPTPFJetVars{
499	amagnan	1.1	double chargedHadronEnergy;
500			double chargedHadronEnergyFraction;
501			double neutralHadronEnergy;
502			double neutralHadronEnergyFraction;
503			double chargedEmEnergy;
504			double chargedEmEnergyFraction;
505			double neutralEmEnergy;
506			double neutralEmEnergyFraction;
507			double chargedMultiplicity;
508	amagnan	1.9	double muonMultiplicity;
509	amagnan	1.24	JPTPFJetVars():
510			chargedHadronEnergy(0),
511			chargedHadronEnergyFraction(0),
512			neutralHadronEnergy(0),
513			neutralHadronEnergyFraction(0),
514			chargedEmEnergy(0),
515			chargedEmEnergyFraction(0),
516			neutralEmEnergy(0),
517			neutralEmEnergyFraction(0),
518			chargedMultiplicity(0),
519			muonMultiplicity(0)
520			{}
521	amagnan	1.9	};
522
523			struct PFJetVars{
524			double chargedMuEnergy;
525			double chargedMuEnergyFraction;
526	amagnan	1.1	double neutralMultiplicity;
527	agilbert	1.22	unsigned numberOfDaughters;
528			double HFHadronEnergy;
529	amagnan	1.24	PFJetVars():
530			chargedMuEnergy(0),
531			chargedMuEnergyFraction(0),
532			neutralMultiplicity(0),
533			numberOfDaughters(0),
534			HFHadronEnergy(0)
535			{}
536	amagnan	1.1	};
537
538			struct JetBtagVars{
539			double cSV;
540			double cSVMVA;
541			double iPMVA;
542			double bProba;
543			double probability;
544	amagnan	1.8	double sSVHE;
545			double sSVHP;
546	amagnan	1.4	double softElectronByPt;
547			double softElectronByIP3d;
548	amagnan	1.1	double softMuon;
549	amagnan	1.4	double softMuonByPt;
550			double softMuonByIP3d;
551	amagnan	1.1	double tCHE;
552			double tCHP;
553	amagnan	1.24	JetBtagVars():
554			cSV(0),
555			cSVMVA(0),
556			iPMVA(0),
557			bProba(0),
558			probability(0),
559			sSVHE(0),
560			sSVHP(0),
561			softElectronByPt(0),
562			softElectronByIP3d(0),
563			softMuon(0),
564			softMuonByPt(0),
565			softMuonByIP3d(0),
566			tCHE(0),
567			tCHP(0)
568			{}
569	amagnan	1.1	};
570
571	amagnan	1.4
572			struct JetIDVars{
573			double fHPD;
574			double fRBX;
575	amagnan	1.9	int n90Hits;
576	amagnan	1.4	//double fSubDetector1;
577			//double fSubDetector2;
578			//double fSubDetector3;
579			//double fSubDetector4;
580	amagnan	1.10	double restrictedEMF;
581			int nHCALTowers;
582			int nECALTowers;
583	amagnan	1.4	//double approximatefHPD;
584			//double approximatefRBX;
585	amagnan	1.9	int hitsInN90;
586	amagnan	1.4	// muon hits id
587			//int numberOfHits2RPC;
588			//int numberOfHits3RPC;
589			//int numberOfHitsRPC;
590	amagnan	1.24	JetIDVars():
591			fHPD(0),
592			fRBX(0),
593			n90Hits(0),
594			restrictedEMF(0),
595			nHCALTowers(0),
596			nECALTowers(0),
597			hitsInN90(0)
598			{}
599	amagnan	1.4	};
600
601	amagnan	1.18	struct JetECorVars{
602	amagnan	1.19	std::vector<std::pair<std::string,double> > Levels;
603	amagnan	1.24	JetECorVars(){
604			Levels.clear();
605			}
606	amagnan	1.18	};
607
608	amagnan	1.4
609	amagnan	1.1	struct MetVars{
610			double mET;
611			double mEx;
612			double mEy;
613			double sumET;
614			double phi;
615			double mEtSig;
616	amagnan	1.24	MetVars():
617			mET(0),
618			mEx(0),
619			mEy(0),
620			sumET(0),
621			phi(0),
622			mEtSig(0)
623			{}
624
625	amagnan	1.1	};
626
627			struct TriggerVars{
628			std::string name;
629			unsigned int index;
630			bool accept;
631	amagnan	1.24	TriggerVars():
632			name(""),
633			index(0),
634			accept(false)
635			{}
636	amagnan	1.1	};
637
638	amagnan	1.4
639			struct VertexVars{
640			std::vector<float> trackWeights;
641	agilbert	1.21	std::vector<float> trackpT;
642	amagnan	1.4	double chi2;
643			double ndof;
644			double x;
645			double y;
646			double z;
647	amagnan	1.16	double rho;
648	amagnan	1.4	double xError;
649			double yError;
650			double zError;
651			double cov01;
652			double cov02;
653			double cov12;
654	amagnan	1.24	VertexVars():
655			chi2(0),
656			ndof(0),
657			x(0),
658			y(0),
659			z(0),
660			rho(0),
661			xError(0),
662			yError(0),
663			zError(0),
664			cov01(0),
665			cov02(0),
666			cov12(0)
667			{
668			trackWeights.clear();
669			trackpT.clear();
670			}
671	amagnan	1.4	};
672
673	amagnan	1.14	struct BeamSpotVars{
674			double x0;
675			double y0;
676			double z0;
677			double sigmaZ;
678			double BeamWidthX;
679			double BeamWidthY;
680			double x0Error;
681			double y0Error;
682			double z0Error;
683			double sigmaZ0Error;
684			double BeamWidthXError;
685			double BeamWidthYError;
686			//double dxdz;
687			//double dydz;
688			//double dxdzError;
689			//double dydzError;
690			//BeamType type();
691			//enum BeamType { Unknown=-1, Fake=0, LHC=1, Tracker=2 };
692			//double emittanceX() const { return emittanceX_; }
693			//double emittanceY() const { return emittanceY_; }
694			//double betaStar() const { return betaStar_; }
695	amagnan	1.24	BeamSpotVars():
696			x0(0),
697			y0(0),
698			z0(0),
699			sigmaZ(0),
700			BeamWidthX(0),
701			BeamWidthY(0),
702			x0Error(0),
703			y0Error(0),
704			z0Error(0),
705			sigmaZ0Error(0),
706			BeamWidthXError(0),
707			BeamWidthYError(0)
708			{}
709
710	amagnan	1.14	};
711
712	amagnan	1.2	double DeltaPhi(const double phi1, const double phi2);
713
714	amagnan	1.24	template <class T1, class T2>
715			double DeltaR(const T1 & v1, const T2 & v2){
716			double dEta = v1.eta - v2.eta;
717	agilbert	1.26	double dPhi = fabs(v1.phi - v2.phi);
718			if (dPhi > TMath::Pi()) dPhi = (2.0*TMath::Pi() - dPhi);
719	amagnan	1.24	return sqrt(dEtadEta+dPhidPhi);
720			}
721
722	amagnan	1.2	double SameSign(const BaseVars & v1, const BaseVars & v2);
723
724			double OppSign(const BaseVars & v1, const BaseVars & v2);
725
726			TLorentzVector FourMomentum(const BaseVars & v, const double scale=1) ;
727
728			double TransverseMass(const BaseVars & leg1,
729			const BaseVars & leg2,
730			const double mEx,
731			const double mEy);
732
733			double TransverseMass(const BaseVars & leg1,
734			const double mEx,
735			const double mEy);
736
737			TLorentzVector FourMomentumCDFmethod(const BaseVars & leg1,
738			const BaseVars & leg2,
739			double mEx,
740			double mEy);
741
742			TLorentzVector FourMomentumCollinearApprox(const BaseVars & leg1,
743			const BaseVars & leg2,
744			double mEx,
745			double mEy);
746	amagnan	1.1
747	amagnan	1.9	double EtaDetector(const BaseVars & v1);
748			double EtaDetector(const GenVars & v1);
749
750	amagnan	1.2	}//namespace
751	amagnan	1.1
752			#endif