1 |
loizides |
1.1 |
//--------------------------------------------------------------------------------------------------
|
2 |
|
|
// $Id: $
|
3 |
|
|
//
|
4 |
|
|
// JetPlusIsoTrackSelMod
|
5 |
|
|
//
|
6 |
|
|
// This module selects events containing one lepton and one isolated track
|
7 |
|
|
// Both gsfTracks and tracker tracks are considered.
|
8 |
|
|
//
|
9 |
|
|
// Authors: S. Xie
|
10 |
|
|
//--------------------------------------------------------------------------------------------------
|
11 |
|
|
|
12 |
|
|
#ifndef MITPHYSICS_SELMODS_JETPLUSISOTRACKSELMOD_H
|
13 |
|
|
#define MITPHYSICS_SELMODS_JETPLUSISOTRACKSELMOD_H
|
14 |
|
|
|
15 |
|
|
#include "MitAna/TreeMod/interface/BaseSelMod.h"
|
16 |
|
|
#include <TH1D.h>
|
17 |
|
|
|
18 |
|
|
namespace mithep
|
19 |
|
|
{
|
20 |
|
|
class JetPlusIsoTrackSelMod : public BaseSelMod
|
21 |
|
|
{
|
22 |
|
|
public:
|
23 |
|
|
JetPlusIsoTrackSelMod(const char *name="JetPlusIsoTrackSelMod",
|
24 |
|
|
const char *title="Jet plus isolated track selection module");
|
25 |
|
|
~JetPlusIsoTrackSelMod() {}
|
26 |
|
|
|
27 |
|
|
const char *GetJetColName() const { return fJetColName; }
|
28 |
|
|
const char *GetTrackerTrackColName() const { return fTrackerTrackColName; }
|
29 |
|
|
const char *GetGsfTrackColName() const { return fGsfTrackColName; }
|
30 |
|
|
Double_t GetJetEtaMin() const { return fJetEtaMin; }
|
31 |
|
|
Double_t GetJetEtaMax() const { return fJetEtaMax; }
|
32 |
|
|
Double_t GetJetPtMin() const { return fJetPtMin; }
|
33 |
|
|
Double_t GetJetPtMax() const { return fJetPtMax; }
|
34 |
|
|
Double_t GetTrackEtaMin() const { return fTrackEtaMin; }
|
35 |
|
|
Double_t GetTrackEtaMax() const { return fTrackEtaMax; }
|
36 |
|
|
Double_t GetTrackPtMin() const { return fTrackPtMin; }
|
37 |
|
|
Double_t GetTrackPtMax() const { return fTrackPtMax; }
|
38 |
|
|
void SetJetColName(const char *n) { fJetColName=n; }
|
39 |
|
|
void SetTrackerTrackColName(const char *n) { fTrackerTrackColName=n; }
|
40 |
|
|
void SetGsfTrackColName(const char *n) { fGsfTrackColName=n; }
|
41 |
|
|
void SetJetEtaMin(Double_t e) { fJetEtaMin = e; }
|
42 |
|
|
void SetJetEtaMax(Double_t e) { fJetEtaMax = e; }
|
43 |
|
|
void SetJetPtMin(Double_t pt) { fJetPtMin = pt; }
|
44 |
|
|
void SetJetPtMax(Double_t pt) { fJetPtMax = pt; }
|
45 |
|
|
void SetTrackEtaMin(Double_t e) { fTrackEtaMin = e; }
|
46 |
|
|
void SetTrackEtaMax(Double_t e) { fTrackEtaMax = e; }
|
47 |
|
|
void SetTrackPtMin(Double_t pt) { fTrackPtMin = pt; }
|
48 |
|
|
void SetTrackPtMax(Double_t pt) { fTrackPtMax = pt; }
|
49 |
|
|
|
50 |
|
|
protected:
|
51 |
|
|
void Process();
|
52 |
|
|
void SlaveBegin();
|
53 |
|
|
|
54 |
|
|
TString fJetColName; //name of input lepton collection
|
55 |
|
|
TString fTrackerTrackColName; //name of input lepton collection
|
56 |
|
|
TString fGsfTrackColName; //name of input lepton collection
|
57 |
|
|
Double_t fJetPtMin; //minimum pt required (def = 0 GeV)
|
58 |
|
|
Double_t fJetPtMax; //maximum pt required (def = 5000 GeV)
|
59 |
|
|
Double_t fJetEtaMin; //minimum eta required (def = -10)
|
60 |
|
|
Double_t fJetEtaMax; //maximum eta required (def = +10)
|
61 |
|
|
Double_t fTrackPtMin; //minimum pt required (def = 0 GeV)
|
62 |
|
|
Double_t fTrackPtMax; //maximum pt required (def = 5000 GeV)
|
63 |
|
|
Double_t fTrackEtaMin; //minimum eta required (def = -10)
|
64 |
|
|
Double_t fTrackEtaMax; //maximum eta required (def = +10)
|
65 |
|
|
const JetCol *fJetCol; //!pointer to collection
|
66 |
|
|
const TrackCol *fTrackerTrackCol; //!pointer to collection
|
67 |
|
|
const TrackCol *fGsfTrackCol; //!pointer to collection
|
68 |
|
|
TH1D *fNAccCounters; //!acceptance histogram
|
69 |
|
|
|
70 |
|
|
ClassDef(JetPlusIsoTrackSelMod,1) // Jet plus isolated track selection module
|
71 |
|
|
};
|
72 |
|
|
}
|
73 |
|
|
#endif
|