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 eta;
    float sigmaEtaEta;
    float sigmaIEtaIEta;
    float e1x5;
    float e2x5Max;
    float e5x5;
    float eOverP;
    SCVars():
      eta(0),
      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);
  double DeltaR(const TLorentzVector & v1, const TLorentzVector & v2);

  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 TLorentzVector & leg1, 
                        const TLorentzVector & leg2, 
                        const double mEx, 
                        const double mEy);
 
  double TransverseMass(const TLorentzVector & leg1, 
                        const double mEx, 
                        const double mEy);

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

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

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

}//namespace

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