Fork me on GitHub

source: git/modules/DualReadoutCalorimeter.h@ dc883b4

Last change on this file since dc883b4 was a1b19ea, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago

first attempt at energy flow with Dual Readout

  • Property mode set to 100644
File size: 3.3 KB
Line 
1/*
2 * Delphes: a framework for fast simulation of a generic collider experiment
3 * Copyright (C) 2012-2014 Universite catholique de Louvain (UCL), Belgium
4 *
5 * This program is free software: you can redistribute it and/or modify
6 * it under the terms of the GNU General Public License as published by
7 * the Free Software Foundation, either version 3 of the License, or
8 * (at your option) any later version.
9 *
10 * This program is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 * GNU General Public License for more details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with this program. If not, see <http://www.gnu.org/licenses/>.
17 */
18
19#ifndef DualReadoutCalorimeter_h
20#define DualReadoutCalorimeter_h
21
22/** \class DualReadoutCalorimeter
23 *
24 * Fills DualReadoutCalorimeter towers, performs DualReadoutCalorimeter resolution smearing,
25 * and creates energy flow objects (tracks, photons, and neutral hadrons).
26 *
27 * \author M. Selvaggi - CERN
28 *
29 */
30
31#include "classes/DelphesModule.h"
32
33#include <map>
34#include <set>
35#include <vector>
36
37class TObjArray;
38class DelphesFormula;
39class Candidate;
40
41class DualReadoutCalorimeter: public DelphesModule
42{
43public:
44
45 DualReadoutCalorimeter();
46 ~DualReadoutCalorimeter();
47
48 void Init();
49 void Process();
50 void Finish();
51
52private:
53
54 typedef std::map< Long64_t, std::pair< Double_t, Double_t > > TFractionMap; //!
55 typedef std::map< Double_t, std::set< Double_t > > TBinMap; //!
56
57 Candidate *fTower;
58 Double_t fTowerEta, fTowerPhi, fTowerEdges[4];
59 Double_t fECalTowerEnergy, fHCalTowerEnergy;
60 Double_t fECalTrackEnergy, fHCalTrackEnergy;
61 Double_t fTrackEnergy;
62
63 Double_t fTimingEnergyMin;
64 Bool_t fElectronsFromTrack;
65
66 Int_t fTowerTrackHits, fTowerPhotonHits;
67
68 Double_t fECalEnergyMin;
69 Double_t fHCalEnergyMin;
70 Double_t fEnergyMin;
71
72 Double_t fECalEnergySignificanceMin;
73 Double_t fHCalEnergySignificanceMin;
74 Double_t fEnergySignificanceMin;
75
76 Double_t fECalTrackSigma;
77 Double_t fHCalTrackSigma;
78 Double_t fTrackSigma;
79
80 Bool_t fSmearTowerCenter;
81
82 TFractionMap fFractionMap; //!
83 TBinMap fBinMap; //!
84
85 std::vector < Double_t > fEtaBins;
86 std::vector < std::vector < Double_t >* > fPhiBins;
87
88 std::vector < Long64_t > fTowerHits;
89
90 std::vector < Double_t > fECalTowerFractions;
91 std::vector < Double_t > fHCalTowerFractions;
92
93 std::vector < Double_t > fECalTrackFractions;
94 std::vector < Double_t > fHCalTrackFractions;
95
96 DelphesFormula *fECalResolutionFormula; //!
97 DelphesFormula *fHCalResolutionFormula; //!
98
99 TIterator *fItParticleInputArray; //!
100 TIterator *fItTrackInputArray; //!
101
102 const TObjArray *fParticleInputArray; //!
103 const TObjArray *fTrackInputArray; //!
104
105 TObjArray *fTowerOutputArray; //!
106 TObjArray *fPhotonOutputArray; //!
107
108 TObjArray *fEFlowTrackOutputArray; //!
109 TObjArray *fEFlowPhotonOutputArray; //!
110 TObjArray *fEFlowNeutralHadronOutputArray; //!
111
112 TObjArray *fECalTowerTrackArray; //!
113 TIterator *fItECalTowerTrackArray; //!
114
115 TObjArray *fHCalTowerTrackArray; //!
116 TIterator *fItHCalTowerTrackArray; //!
117
118 TObjArray *fTowerTrackArray; //!
119 TIterator *fItTowerTrackArray; //!
120
121 void FinalizeTower();
122 Double_t LogNormal(Double_t mean, Double_t sigma);
123
124 ClassDef(DualReadoutCalorimeter, 1)
125};
126
127#endif
Note: See TracBrowser for help on using the repository browser.