1 |
loizides |
1.1 |
//--------------------------------------------------------------------------------------------------
|
2 |
loizides |
1.20 |
// $Id: ElectronIDMod.h,v 1.19 2009/06/02 05:30:44 loizides Exp $
|
3 |
loizides |
1.1 |
//
|
4 |
|
|
// ElectronIDMod
|
5 |
|
|
//
|
6 |
loizides |
1.4 |
// This module applies electron identification criteria and exports a pointer to a collection
|
7 |
loizides |
1.5 |
// of "good electrons" according to the specified identification scheme.
|
8 |
loizides |
1.1 |
//
|
9 |
loizides |
1.13 |
// See http://indico.cern.ch/contributionDisplay.py?contribId=1&confId=42251
|
10 |
|
|
//
|
11 |
loizides |
1.4 |
// Authors: S.Xie, C.Loizides
|
12 |
loizides |
1.1 |
//--------------------------------------------------------------------------------------------------
|
13 |
|
|
|
14 |
loizides |
1.4 |
#ifndef MITPHYSICS_MODS_ELECTRONIDMOD_H
|
15 |
|
|
#define MITPHYSICS_MODS_ELECTRONIDMOD_H
|
16 |
loizides |
1.1 |
|
17 |
|
|
#include "MitAna/TreeMod/interface/BaseMod.h"
|
18 |
loizides |
1.20 |
#include "MitAna/DataTree/interface/ElectronFwd.h"
|
19 |
|
|
#include "MitAna/DataTree/interface/VertexFwd.h"
|
20 |
|
|
#include "MitAna/DataTree/interface/DecayParticleFwd.h"
|
21 |
loizides |
1.1 |
|
22 |
|
|
namespace mithep
|
23 |
|
|
{
|
24 |
|
|
class ElectronIDMod : public BaseMod
|
25 |
|
|
{
|
26 |
|
|
public:
|
27 |
|
|
ElectronIDMod(const char *name="ElectronIDMod",
|
28 |
loizides |
1.4 |
const char *title="Electron identification module");
|
29 |
|
|
|
30 |
loizides |
1.13 |
Bool_t GetApplyConversionFilter() const { return fApplyConvFilter; }
|
31 |
ceballos |
1.14 |
Bool_t GetApplyD0Cut() const { return fApplyD0Cut; }
|
32 |
loizides |
1.8 |
Double_t GetCaloIsoCut() const { return fCaloIsolationCut; }
|
33 |
|
|
Double_t GetEcalJurIsoCut() const { return fEcalJuraIsoCut; }
|
34 |
|
|
const char *GetGoodName() const { return GetGoodElectronsName(); }
|
35 |
|
|
const char *GetGoodElectronsName() const { return fGoodElectronsName; }
|
36 |
|
|
Double_t GetHcalIsoCut() const { return fHcalIsolationCut; }
|
37 |
|
|
Double_t GetIDLikelihoodCut() const { return fIDLikelihoodCut; }
|
38 |
|
|
const char *GetIDType() const { return fElectronIDType; }
|
39 |
|
|
const char *GetInputName() const { return fElectronBranchName; }
|
40 |
|
|
const char *GetIsoType() const { return fElectronIsoType; }
|
41 |
|
|
const char *GetOutputName() const { return GetGoodElectronsName(); }
|
42 |
|
|
Double_t GetPtMin() const { return fElectronPtMin; }
|
43 |
|
|
Double_t GetTrackIsoCut() const { return fTrackIsolationCut; }
|
44 |
ceballos |
1.17 |
Bool_t GetReverseIsoCut() const { return fReverseIsoCut; }
|
45 |
ceballos |
1.18 |
Bool_t GetReverseD0Cut() const { return fReverseD0Cut; }
|
46 |
loizides |
1.13 |
void SetApplyConversionFilter(Bool_t b) { fApplyConvFilter = b; }
|
47 |
ceballos |
1.14 |
void SetApplyD0Cut(Bool_t b) { fApplyD0Cut = b; }
|
48 |
loizides |
1.8 |
void SetCaloIsoCut(Double_t cut) { fCaloIsolationCut = cut; }
|
49 |
|
|
void SetEcalJurIsoCut(Double_t cut) { fEcalJuraIsoCut = cut; }
|
50 |
|
|
void SetGoodName(const char *n) { SetGoodElectronsName(n); }
|
51 |
|
|
void SetGoodElectronsName(const char *n) { fGoodElectronsName = n; }
|
52 |
|
|
void SetHcalIsoCut(Double_t cut) { fHcalIsolationCut = cut; }
|
53 |
|
|
void SetIDLikelihoodCut(Double_t cut) { fIDLikelihoodCut = cut; }
|
54 |
|
|
void SetIDType(const char *type) { fElectronIDType = type; }
|
55 |
|
|
void SetInputName(const char *n) { fElectronBranchName = n; }
|
56 |
|
|
void SetIsoType(const char *type) { fElectronIsoType = type; }
|
57 |
|
|
void SetOutputName(const char *n) { SetGoodElectronsName(n); }
|
58 |
|
|
void SetPtMin(Double_t pt) { fElectronPtMin = pt; }
|
59 |
|
|
void SetTrackIsoCut(Double_t cut) { fTrackIsolationCut = cut; }
|
60 |
ceballos |
1.12 |
void SetD0Cut(Double_t cut) { fD0Cut = cut; }
|
61 |
ceballos |
1.17 |
void SetReverseIsoCut(Bool_t b) { fReverseIsoCut = b; }
|
62 |
ceballos |
1.18 |
void SetReverseD0Cut(Bool_t b) { fReverseD0Cut = b; }
|
63 |
loizides |
1.4 |
|
64 |
|
|
enum EElIdType {
|
65 |
|
|
kIdUndef = 0, //not defined
|
66 |
|
|
kTight, //"Tight"
|
67 |
|
|
kLoose, //"Loose"
|
68 |
|
|
kLikelihood, //"Likelihood"
|
69 |
loizides |
1.9 |
kNoId, //"NoId"
|
70 |
loizides |
1.4 |
kCustomId //"Custom"
|
71 |
|
|
};
|
72 |
|
|
enum EElIsoType {
|
73 |
|
|
kIsoUndef = 0, //not defined
|
74 |
|
|
kTrackCalo, //"TrackCalo"
|
75 |
|
|
kTrackJura, //"TrackJura"
|
76 |
|
|
kTrackJuraSliding, //"TrackJuraSliding"
|
77 |
|
|
kNoIso, //"NoIso"
|
78 |
|
|
kCustomIso //"Custom"
|
79 |
|
|
};
|
80 |
|
|
|
81 |
loizides |
1.1 |
protected:
|
82 |
ceballos |
1.12 |
void Process();
|
83 |
|
|
void SlaveBegin();
|
84 |
loizides |
1.7 |
|
85 |
ceballos |
1.12 |
TString fElectronBranchName; //name of electron collection (input)
|
86 |
|
|
TString fConversionBranchName; //name of electron collection (input)
|
87 |
loizides |
1.16 |
TString fGoodElectronsName; //name of exported "good electrons" col
|
88 |
ceballos |
1.12 |
TString fVertexName; //name of vertex collection
|
89 |
|
|
TString fElectronIDType; //type of electron ID we impose
|
90 |
|
|
TString fElectronIsoType; //type of electron Isolation we impose
|
91 |
|
|
Double_t fElectronPtMin; //min pt cut
|
92 |
|
|
Double_t fIDLikelihoodCut; //cut value for ID likelihood
|
93 |
|
|
Double_t fTrackIsolationCut; //cut value for track isolation
|
94 |
|
|
Double_t fCaloIsolationCut; //cut value for calo isolation
|
95 |
|
|
Double_t fEcalJuraIsoCut; //cut value for ecal jurassic isolation
|
96 |
|
|
Double_t fHcalIsolationCut; //cut value for hcal isolation
|
97 |
loizides |
1.13 |
Bool_t fApplyConvFilter; //whether remove conversions
|
98 |
loizides |
1.16 |
Bool_t fApplyD0Cut; //whether apply d0 cut
|
99 |
loizides |
1.13 |
Double_t fD0Cut; //max d0
|
100 |
loizides |
1.19 |
Bool_t fReverseIsoCut; //apply reversion iso cut (default=0)
|
101 |
|
|
Bool_t fReverseD0Cut; //apply reversion d0 cut (default=0)
|
102 |
loizides |
1.13 |
EElIdType fElIdType; //!identification scheme
|
103 |
|
|
EElIsoType fElIsoType; //!isolation scheme
|
104 |
ceballos |
1.12 |
const ElectronCol *fElectrons; //!electron collection
|
105 |
|
|
const DecayParticleCol *fConversions; //!conversion collection
|
106 |
loizides |
1.16 |
const VertexCol *fVertices; //!vertices branches
|
107 |
ceballos |
1.17 |
|
108 |
loizides |
1.8 |
ClassDef(ElectronIDMod, 1) // Electron identification module
|
109 |
loizides |
1.1 |
};
|
110 |
|
|
}
|
111 |
|
|
#endif
|