1 | #ifndef PileUpJetID_h
|
---|
2 | #define PileUpJetID_h
|
---|
3 |
|
---|
4 | /** \class PileUpJetID
|
---|
5 | *
|
---|
6 | * CMS PileUp Jet ID Variables
|
---|
7 | *
|
---|
8 | * \author S. Zenz
|
---|
9 | *
|
---|
10 | */
|
---|
11 |
|
---|
12 | #include "classes/DelphesModule.h"
|
---|
13 |
|
---|
14 | #include <deque>
|
---|
15 |
|
---|
16 | class TObjArray;
|
---|
17 | class DelphesFormula;
|
---|
18 |
|
---|
19 | class PileUpJetID: public DelphesModule
|
---|
20 | {
|
---|
21 | public:
|
---|
22 | PileUpJetID();
|
---|
23 | ~PileUpJetID();
|
---|
24 |
|
---|
25 | void Init();
|
---|
26 | void Process();
|
---|
27 | void Finish();
|
---|
28 |
|
---|
29 | private:
|
---|
30 | Double_t fJetPTMin;
|
---|
31 | Double_t fParameterR;
|
---|
32 |
|
---|
33 | Double_t fMeanSqDeltaRMaxBarrel; // |eta| < 1.5
|
---|
34 | Double_t fBetaMinBarrel; // |eta| < 2.5
|
---|
35 | Double_t fMeanSqDeltaRMaxEndcap; // 1.5 < |eta| < 4.0
|
---|
36 | Double_t fBetaMinEndcap; // 1.5 < |eta| < 4.0
|
---|
37 | Double_t fMeanSqDeltaRMaxForward; // |eta| > 4.0
|
---|
38 |
|
---|
39 | Double_t fNeutralPTMin;
|
---|
40 | Double_t fJetPTMinForNeutrals;
|
---|
41 |
|
---|
42 | /*
|
---|
43 | JAY
|
---|
44 | ---
|
---|
45 |
|
---|
46 | |Eta|<1.5
|
---|
47 |
|
---|
48 | meanSqDeltaR betaStar SigEff BgdEff
|
---|
49 | 0.13 0.92 96% 8%
|
---|
50 | 0.13 0.95 97% 16%
|
---|
51 | 0.13 0.97 98% 27%
|
---|
52 |
|
---|
53 | |Eta|>1.5
|
---|
54 |
|
---|
55 | meanSqDeltaR betaStar SigEff BgdEff
|
---|
56 | 0.14 0.91 95% 15%
|
---|
57 | 0.14 0.94 97% 19%
|
---|
58 | 0.14 0.97 98% 29%
|
---|
59 |
|
---|
60 | BRYAN
|
---|
61 | -----
|
---|
62 |
|
---|
63 | Barrel (MeanSqDR, Beta, sig eff, bg eff):
|
---|
64 | 0.10, 0.08, 90%, 8%
|
---|
65 | 0.11, 0.12, 90%, 6%
|
---|
66 | 0.13, 0.16, 89%, 5%
|
---|
67 |
|
---|
68 | Endcap (MeanSqDR, Beta, sig eff, bg eff):
|
---|
69 | 0.07, 0.06, 89%, 4%
|
---|
70 | 0.08, 0.08, 92%, 6%
|
---|
71 | 0.09, 0.08, 95%, 10%
|
---|
72 | 0.10, 0.08, 97%, 13%
|
---|
73 |
|
---|
74 | SETH GUESSES FOR |eta| > 4.0
|
---|
75 | ----------------------------
|
---|
76 |
|
---|
77 | MeanSqDeltaR
|
---|
78 | 0.07
|
---|
79 | 0.10
|
---|
80 | 0.14
|
---|
81 | 0.2
|
---|
82 | */
|
---|
83 |
|
---|
84 | // If set to true, may have weird results for PFCHS
|
---|
85 | // If set to false, uses everything within dR < fParameterR even if in other jets &c.
|
---|
86 | // Results should be very similar for PF
|
---|
87 | Int_t fUseConstituents;
|
---|
88 |
|
---|
89 | Bool_t fAverageEachTower;
|
---|
90 |
|
---|
91 | TIterator *fItJetInputArray; //!
|
---|
92 |
|
---|
93 | const TObjArray *fJetInputArray; //!
|
---|
94 |
|
---|
95 | const TObjArray *fTrackInputArray; // SCZ
|
---|
96 | const TObjArray *fNeutralInputArray;
|
---|
97 |
|
---|
98 | TIterator *fItTrackInputArray; // SCZ
|
---|
99 | TIterator *fItNeutralInputArray; // SCZ
|
---|
100 |
|
---|
101 | TObjArray *fOutputArray; //!
|
---|
102 | TObjArray *fNeutralsInPassingJets; // SCZ
|
---|
103 |
|
---|
104 | ClassDef(PileUpJetID, 2)
|
---|
105 | };
|
---|
106 |
|
---|
107 | #endif
|
---|