HbbAnalysis/interface/Objects.hh

#ifndef HbbAnalysis_Objects_hh
#define HbbAnalysis_Objects_hh

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

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