26 |
|
#include "MitAna/DataTree/interface/SuperClusterCol.h" |
27 |
|
#include "MitCommon/MathTools/interface/MathUtils.h" |
28 |
|
|
29 |
+ |
#include "MitAna/DataTree/interface/PFCandidateCol.h" |
30 |
+ |
|
31 |
|
class TRandom3; |
32 |
|
|
33 |
|
namespace mithep { |
71 |
|
}; |
72 |
|
|
73 |
|
enum eScaleCats { |
74 |
< |
kEBlowEtaGoldCenter = 0, |
73 |
< |
kEBlowEtaGoldGap, |
74 |
> |
kEBlowEtaGold = 0, |
75 |
|
kEBlowEtaBad, |
76 |
|
kEBhighEtaGold, |
77 |
|
kEBhighEtaBad, |
80 |
|
kEEhighEtaGold, |
81 |
|
kEEhighEtaBad |
82 |
|
}; |
83 |
< |
|
83 |
> |
|
84 |
> |
enum ShowerShapeScales { |
85 |
> |
kNoShowerShapeScaling = 0, |
86 |
> |
k2011ShowerShape, |
87 |
> |
k2012ShowerShape |
88 |
> |
}; |
89 |
|
|
90 |
|
static eScaleCats EScaleCat(const Photon *p); |
91 |
|
|
96 |
|
static void ScalePhotonR9(Photon *p, Double_t scale); |
97 |
|
static void ScalePhotonError(Photon *p, Double_t scale); |
98 |
|
|
99 |
+ |
static void ScalePhotonShowerShapes(Photon *p, ShowerShapeScales scale); |
100 |
+ |
|
101 |
+ |
|
102 |
|
static Bool_t PassSinglePhotonPresel(const Photon *p,const ElectronCol *els, const DecayParticleCol *conversions, const BaseVertex *bs, const TrackCol* trackCol, const Vertex *vtx, double rho, Bool_t applyElectronVeto = kTRUE, Bool_t invertElectronVeto = kFALSE); |
103 |
+ |
static Bool_t PassSinglePhotonPreselPFISO(const Photon *p,const ElectronCol *els, const DecayParticleCol *conversions, const BaseVertex *bs, const TrackCol* trackCol,const Vertex *vtx, double rho, const PFCandidateCol *fPFCands, Bool_t applyElectronVeto = kTRUE, Bool_t invertElectronVeto = kFALSE); |
104 |
+ |
static Bool_t PassSinglePhotonPreselPFISONoEcal(const Photon *p,const ElectronCol *els, const DecayParticleCol *conversions, const BaseVertex *bs, const TrackCol* trackCol,const Vertex *vtx, double rho, const PFCandidateCol *fPFCands, Bool_t applyElectronVeto = kTRUE, Bool_t invertElectronVeto = kFALSE); |
105 |
+ |
static Bool_t PassSinglePhotonPreselPFISONoEcalNoPFChargedIso(const Photon *p,const ElectronCol *els, const DecayParticleCol *conversions, const BaseVertex *bs, const TrackCol* trackCol,const Vertex *vtx, double rho, const PFCandidateCol *fPFCands, Bool_t applyElectronVeto = kTRUE, Bool_t invertElectronVeto = kFALSE); |
106 |
+ |
static Bool_t PassSinglePhotonPreselPFISO_NoTrigger(const Photon *p,const ElectronCol *els, const DecayParticleCol *conversions, const BaseVertex *bs, const TrackCol* trackCol,const Vertex *vtx, double rho, const PFCandidateCol *fPFCands, Bool_t applyElectronVeto = kTRUE, Bool_t invertElectronVeto = kFALSE); |
107 |
|
static Bool_t PassConversionId(const Photon *p, const DecayParticle *c); |
108 |
|
static Bool_t PassElectronVeto(const Photon *p, const ElectronCol *els); |
109 |
|
static Double_t ElectronVetoCiC(const Photon *p, const ElectronCol *els); |
125 |
|
|
126 |
|
static const DecayParticle *MatchedCiCConversion(const Photon *p, const DecayParticleCol *conversions, |
127 |
|
Double_t dPhiMin=0.1, Double_t dEtaMin=0.1,Double_t dRMin=0.1, |
128 |
< |
bool print = false); |
128 |
> |
bool print = false, |
129 |
> |
int* numLegs = NULL, int* convIdx = NULL); // for debugging |
130 |
|
|
131 |
|
|
132 |
|
static const Electron *MatchedElectron(const Photon *p, const ElectronCol *els); |
133 |
|
static const Photon *MatchedPhoton(const Electron *e, const PhotonCol *phs); |
134 |
|
static const SuperCluster *MatchedSC(const SuperCluster *psc, const SuperClusterCol *scs, Double_t drMin=0.3); |
135 |
< |
|
135 |
> |
|
136 |
> |
static const SuperCluster *MatchedPFSC(const SuperCluster *psc, const PhotonCol *pfphos, const ElectronCol *eles, Double_t drMin=0.1); |
137 |
> |
|
138 |
|
static bool PassCiCSelection(const Photon* ph, |
139 |
|
const Vertex* vtx, |
140 |
|
const TrackCol* trackCol, |
144 |
|
bool applyEleVeto = true, |
145 |
|
bool print = false, float* kin=NULL); |
146 |
|
|
147 |
< |
|
147 |
> |
static bool PassCiCPFIsoSelection(const Photon* ph, |
148 |
> |
const Vertex* vtx, |
149 |
> |
const PFCandidateCol* pfCol, |
150 |
> |
const VertexCol* vtxCol, |
151 |
> |
double rho, double ptmin, |
152 |
> |
std::vector<double>* kin = NULL); |
153 |
> |
|
154 |
> |
static bool PassCiCPFIsoSelectionWithEleVeto(const Photon* ph, |
155 |
> |
const ElectronCol *els, |
156 |
> |
const DecayParticleCol *conversions, const BaseVertex *bs, |
157 |
> |
const Vertex* vtx, |
158 |
> |
const PFCandidateCol* pfCol, |
159 |
> |
const VertexCol* vtxCol, |
160 |
> |
double rho, double ptmin, |
161 |
> |
Bool_t applyElectronVeto, Bool_t invertElectronVeto, |
162 |
> |
std::vector<double>* kin= NULL // store variables for debugging... |
163 |
> |
);// add for mono photon |
164 |
> |
|
165 |
> |
static bool PassVgamma2011Selection(const Photon* ph, double rho); |
166 |
> |
|
167 |
|
static const MCParticle *MatchMC(const Particle *ph, const MCParticleCol *c, Bool_t matchElectrons = kFALSE); |
168 |
|
ClassDef(PhotonTools, 0) // Muon tools |
169 |
|
}; |