| 1 |
bendavid |
1.1 |
//--------------------------------------------------------------------------------------------------
|
| 2 |
loizides |
1.2 |
// $Id: PFTau.h,v 1.1 2009/03/19 23:09:36 bendavid Exp $
|
| 3 |
bendavid |
1.1 |
//
|
| 4 |
|
|
// PFTau
|
| 5 |
|
|
//
|
| 6 |
loizides |
1.2 |
// This class holds information about reconstructed tau based on PFCandidates.
|
| 7 |
bendavid |
1.1 |
//
|
| 8 |
|
|
// Authors: J.Bendavid
|
| 9 |
|
|
//--------------------------------------------------------------------------------------------------
|
| 10 |
|
|
|
| 11 |
|
|
#ifndef MITANA_DATATREE_PFTAU_H
|
| 12 |
|
|
#define MITANA_DATATREE_PFTAU_H
|
| 13 |
|
|
|
| 14 |
|
|
#include "MitAna/DataTree/interface/Tau.h"
|
| 15 |
|
|
#include "MitAna/DataCont/interface/RefArray.h"
|
| 16 |
|
|
#include "MitAna/DataCont/interface/Ref.h"
|
| 17 |
|
|
#include "MitAna/DataTree/interface/PFCandidate.h"
|
| 18 |
|
|
#include "MitAna/DataTree/interface/BasicCluster.h"
|
| 19 |
|
|
#include "MitAna/DataTree/interface/PFJet.h"
|
| 20 |
|
|
|
| 21 |
|
|
namespace mithep
|
| 22 |
|
|
{
|
| 23 |
|
|
class PFTau : public Tau
|
| 24 |
|
|
{
|
| 25 |
|
|
public:
|
| 26 |
loizides |
1.2 |
PFTau() : fLeadPFCandSignD0Sig(0), fHCalTotalEOverP(0), fHCalMaxEOverP(0),
|
| 27 |
bendavid |
1.1 |
fHCal3x3EOverP(0), fIsoChargedHadronPtSum(0), fIsoGammaEtSum(0),
|
| 28 |
|
|
fMaxHCalPFClusterEt(0), fEMFraction(0), fECalStripSumEOverP(0),
|
| 29 |
loizides |
1.2 |
fBremRecoveryEOverP(0), fElectronPreIDOutput(0), fCaloCompatibility(0),
|
| 30 |
|
|
fSegmentCompatibility(0), fElectronPreIDDecision(kFALSE),
|
| 31 |
|
|
fMuonDecision(kFALSE) {}
|
| 32 |
bendavid |
1.1 |
|
| 33 |
|
|
PFTau(Double_t px, Double_t py, Double_t pz, Double_t e) :
|
| 34 |
|
|
Tau(px,py,pz,e),
|
| 35 |
loizides |
1.2 |
fLeadPFCandSignD0Sig(0), fHCalTotalEOverP(0), fHCalMaxEOverP(0),
|
| 36 |
bendavid |
1.1 |
fHCal3x3EOverP(0), fIsoChargedHadronPtSum(0), fIsoGammaEtSum(0),
|
| 37 |
|
|
fMaxHCalPFClusterEt(0), fEMFraction(0), fECalStripSumEOverP(0),
|
| 38 |
loizides |
1.2 |
fBremRecoveryEOverP(0), fElectronPreIDOutput(0), fCaloCompatibility(0),
|
| 39 |
|
|
fSegmentCompatibility(0), fElectronPreIDDecision(kFALSE),
|
| 40 |
|
|
fMuonDecision(kFALSE) {}
|
| 41 |
|
|
|
| 42 |
|
|
void AddIsoPFCand(const PFCandidate *p) { fIsoPFCands.Add(p); }
|
| 43 |
|
|
void AddSignalPFCand(const PFCandidate *p) { fSignalPFCands.Add(p); }
|
| 44 |
|
|
Double_t BremRecoveryEOverP() const { return fBremRecoveryEOverP; }
|
| 45 |
|
|
Double_t CaloCompatibility() const { return fCaloCompatibility; }
|
| 46 |
|
|
Double_t ECalStripSumEOverP() const { return fECalStripSumEOverP; }
|
| 47 |
|
|
Double_t EMFraction() const { return fEMFraction; }
|
| 48 |
|
|
const Track *ElectronTrack() const { return fElectronTrack.Obj(); }
|
| 49 |
|
|
Bool_t ElectronPreIDDecision() const { return fElectronPreIDDecision; }
|
| 50 |
|
|
Double_t ElectronPreIDOutput() const { return fElectronPreIDOutput; }
|
| 51 |
|
|
Double_t HCal3x3EOverP() const { return fHCal3x3EOverP; }
|
| 52 |
|
|
Double_t HCalMaxEOverP() const { return fHCalMaxEOverP; }
|
| 53 |
|
|
Double_t HCalTotalEOverP() const { return fHCalTotalEOverP; }
|
| 54 |
|
|
Double_t IsoChargedHadronPtSum() const { return fIsoChargedHadronPtSum; }
|
| 55 |
|
|
Double_t IsoGammaEtSum() const { return fIsoGammaEtSum; }
|
| 56 |
|
|
const PFCandidate *IsoPFCand(UInt_t i) const { return fIsoPFCands.At(i); }
|
| 57 |
|
|
const PFCandidate *LeadChargedHadronPFCand() const { return fLeadChargedHadPFCand.Obj(); }
|
| 58 |
|
|
const PFCandidate *LeadNeutralHadronPFCand() const { return fLeadNeutralPFCand.Obj(); }
|
| 59 |
|
|
const PFCandidate *LeadPFCand() const { return fLeadPFCand.Obj(); }
|
| 60 |
|
|
Double_t LeadPFCandSignD0Sig() const { return fLeadPFCandSignD0Sig; }
|
| 61 |
|
|
Double_t MaxHCalPFClusterEt() const { return fMaxHCalPFClusterEt; }
|
| 62 |
|
|
Bool_t MuonDecision() const { return fMuonDecision; }
|
| 63 |
|
|
UInt_t NIsoPFCandS() const { return fIsoPFCands.Entries(); }
|
| 64 |
|
|
UInt_t NSignalPFCands() const { return fSignalPFCands.GetEntries(); }
|
| 65 |
|
|
EObjType ObjType() const { return kPFTau; }
|
| 66 |
|
|
Double_t SegmentCompatibility() const { return fSegmentCompatibility; }
|
| 67 |
|
|
void SetBremRecoveryEOverP(Double_t x) { fBremRecoveryEOverP = x; }
|
| 68 |
|
|
void SetCaloCompatibility(Double_t x) { fCaloCompatibility = x; }
|
| 69 |
|
|
void SetECalStripSumEOverP(Double_t x) { fECalStripSumEOverP = x; }
|
| 70 |
|
|
void SetEMFraction(Double_t x) { fEMFraction = x; }
|
| 71 |
|
|
void SetElectronPreIDDecision(Bool_t b) { fElectronPreIDDecision = b; }
|
| 72 |
|
|
void SetElectronPreIDOutput(Double_t x) { fElectronPreIDOutput = x; }
|
| 73 |
|
|
void SetElectronTrack(const Track *t) { fElectronTrack = t; }
|
| 74 |
|
|
void SetHCal3x3EOverP(Double_t x) { fHCal3x3EOverP = x; }
|
| 75 |
|
|
void SetHCalMaxEOverP(Double_t x) { fHCalMaxEOverP = x; }
|
| 76 |
|
|
void SetHCalTotalEOverP(Double_t x) { fHCalTotalEOverP = x; }
|
| 77 |
|
|
void SetIsoChargedHadronPtSum(Double_t x) { fIsoChargedHadronPtSum = x; }
|
| 78 |
|
|
void SetIsoGammaEtSum(Double_t x) { fIsoGammaEtSum = x; }
|
| 79 |
|
|
void SetLeadChargedHadronPFCand(const PFCandidate *p)
|
| 80 |
|
|
{ fLeadChargedHadPFCand = p; }
|
| 81 |
|
|
void SetLeadNeutralPFCand(const PFCandidate *p) { fLeadNeutralPFCand = p; }
|
| 82 |
|
|
void SetLeadPFCand(const PFCandidate *p) { fLeadPFCand = p; }
|
| 83 |
|
|
void SetLeadPFCandSignD0Sig(Double_t x) { fLeadPFCandSignD0Sig = x; }
|
| 84 |
|
|
void SetMaxHCalPFClusterEt(Double_t x) { fMaxHCalPFClusterEt = x; }
|
| 85 |
|
|
void SetMuonDecision(Bool_t b) { fMuonDecision = b; }
|
| 86 |
|
|
void SetPFJet(const PFJet *j) { fPFJet = j; }
|
| 87 |
|
|
void SetSegmentCompatibility(Double_t x) { fSegmentCompatibility = x; }
|
| 88 |
|
|
const PFCandidate *SignalPFCand(UInt_t i) const { return fSignalPFCands.At(i); }
|
| 89 |
|
|
const PFJet *SourcePFJet() const { return fPFJet.Obj(); }
|
| 90 |
|
|
const Jet *SourceJet() const { return SourcePFJet(); }
|
| 91 |
bendavid |
1.1 |
|
| 92 |
|
|
protected:
|
| 93 |
loizides |
1.2 |
Double32_t fLeadPFCandSignD0Sig; //[0,0,14]signed lead track D0 significance
|
| 94 |
|
|
Double32_t fHCalTotalEOverP; //[0,0,14]total hcal e / lead ch had pfcand mom
|
| 95 |
|
|
Double32_t fHCalMaxEOverP; //[0,0,14]max hcal e / lead ch had pfcand. mom
|
| 96 |
|
|
Double32_t fHCal3x3EOverP; //[0,0,14]3x3 hcal e / lead ch hadron pfcand. mom
|
| 97 |
|
|
Double32_t fIsoChargedHadronPtSum; //[0,0,14]sum pt of sel. ch had pfcands in iso cone
|
| 98 |
|
|
Double32_t fIsoGammaEtSum; //[0,0,14]sum et of sel. photon pfcands in iso cone
|
| 99 |
|
|
Double32_t fMaxHCalPFClusterEt; //[0,0,14]et of largest et hcal pfcluster
|
| 100 |
|
|
Double32_t fEMFraction; //[0,0,14]em energy fraction
|
| 101 |
|
|
Double32_t fECalStripSumEOverP; //[0,0,14]simple brem recovery e / lead ch had mom
|
| 102 |
|
|
Double32_t fBremRecoveryEOverP; //[0,0,14]brem recovery E / lead charged hadron P
|
| 103 |
|
|
Double32_t fElectronPreIDOutput; //[0,0,14]pfel pre id bdt output to be an el
|
| 104 |
|
|
Double32_t fCaloCompatibility; //[0,0,14]calo comp. for this tau to be a muon
|
| 105 |
|
|
Double32_t fSegmentCompatibility; //[0,0,14]segment comp. for this tau to be a muon
|
| 106 |
|
|
Bool_t fElectronPreIDDecision; //pf electron pre id decision
|
| 107 |
|
|
Bool_t fMuonDecision; //pf muon id decision
|
| 108 |
|
|
Ref<PFCandidate> fLeadPFCand; //leading signal pf candidate (charged or neutral)
|
| 109 |
|
|
Ref<PFCandidate> fLeadChargedHadPFCand; //leading charged hadron signal pf candidate
|
| 110 |
|
|
Ref<PFCandidate> fLeadNeutralPFCand; //leading neutral signal pf candidate
|
| 111 |
|
|
Ref<PFJet> fPFJet; //original reconstructed pf jet
|
| 112 |
|
|
Ref<Track> fElectronTrack; //track corresponding to possible matching el cand.
|
| 113 |
|
|
RefArray<PFCandidate> fSignalPFCands; //selected pf candidates in signal cone
|
| 114 |
|
|
RefArray<PFCandidate> fIsoPFCands; //selected pf candidates in isolation annulus
|
| 115 |
bendavid |
1.1 |
|
| 116 |
|
|
ClassDef(PFTau, 1) // PFTau class
|
| 117 |
|
|
};
|
| 118 |
|
|
}
|
| 119 |
|
|
#endif
|
