test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
calotowermaker_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 calotowermaker = cms.EDProducer("CaloTowersCreator",
4  # Energy threshold for EB 5x5 crystal inclusion [GeV]
5  EBSumThreshold = cms.double(0.2),
6  # Weighting factor for HF short-fiber readouts
7  HF2Weight = cms.double(1.0),
8  # Weighting factor for EB
9  EBWeight = cms.double(1.0),
10  # Label of HFRecHitCollection to use
11  hfInput = cms.InputTag("hfreco"),
12  # Energy threshold for EE crystals-in-tower inclusion [GeV]
13  EESumThreshold = cms.double(0.45),
14  # Energy threshold for HO cell inclusion [GeV]
15  HOThreshold0 = cms.double(1.1),
16  HOThresholdPlus1 = cms.double(3.5),
17  HOThresholdMinus1 = cms.double(3.5),
18  HOThresholdPlus2 = cms.double(3.5),
19  HOThresholdMinus2 = cms.double(3.5),
20  HBGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
21  # Energy threshold for HB cell inclusion [GeV]
22  HBThreshold = cms.double(0.7),
23  EEWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
24  # Energy threshold for long-fiber HF readout inclusion [GeV]
25  HF1Threshold = cms.double(0.5),
26  HF2Weights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
27  HOWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
28  EEGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
29  # Weighting factor for HE 10-degree cells
30  HEDWeight = cms.double(1.0),
31  # Weighting factor for EE
32  EEWeight = cms.double(1.0),
33  # HO on/off flag for tower energy reconstruction
34  UseHO = cms.bool(True),
35  HBWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
36  # Weighting factor for HE 5-degree cells
37  HESWeight = cms.double(1.0),
38  # Weighting factor for HF long-fiber readouts
39  HF1Weight = cms.double(1.0),
40  HF2Grid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
41  HEDWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
42  HF1Grid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
43  EBWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
44  # Weighting factor for HO
45  HOWeight = cms.double(1.0),
46 
47  # Energy threshold for EB crystal inclusion [GeV]
48  EBThreshold = cms.double(0.07),
49  # Energy threshold for EE crystal inclusion [GeV]
50  EEThreshold = cms.double(0.3),
51  # Flags specifying if the above thresholds
52  # should be applied to Et (UseEtEXTreshold='True') or E ('False')
53  # Flags for use of symmetric thresholds: |e|>threshold
54  UseEtEBTreshold = cms.bool(False),
55  UseSymEBTreshold = cms.bool(True),
56  UseEtEETreshold = cms.bool(False),
57  UseSymEETreshold = cms.bool(True),
58 
59 
60  # Label of HBHERecHitCollection to use
61  hbheInput = cms.InputTag("hbhereco"),
62  # Global energy threshold on Hcal [GeV]
63  HcalThreshold = cms.double(-1000.0),
64  # Energy threshold for short-fiber HF readout inclusion [GeV]
65  HF2Threshold = cms.double(0.85),
66 
67  # Energy threshold for 5-degree (phi) HE cell inclusion [GeV]
68  HESThreshold = cms.double(0.8),
69  HF1Weights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
70  # Label of HORecHitCollection to use
71  hoInput = cms.InputTag("horeco"),
72  HESGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
73  #
74  HESWeights = cms.vdouble(1.0, 1.0, 1.0, 1.0, 1.0),
75  # Energy threshold for 10-degree (phi) HE cel inclusion [GeV]
76  HEDThreshold = cms.double(0.8),
77  # Global energy threshold on tower [GeV]
78  EcutTower = cms.double(-1000.0),
79  HEDGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
80  # Label of EcalRecHitCollections to use
81  ecalInputs = cms.VInputTag(cms.InputTag("ecalRecHit","EcalRecHitsEB"), cms.InputTag("ecalRecHit","EcalRecHitsEE")),
82  # Weighting factor for HB
83  HBWeight = cms.double(1.0),
84  HOGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
85  # Energy dependent weights and energy scale to be used
86  EBGrid = cms.vdouble(-1.0, 1.0, 10.0, 100.0, 1000.0),
87  #
88  #
89  # momentum assignment
90  # Method for momentum reconstruction
91  MomConstrMethod = cms.int32(1),
92  # Depth, fraction of the respective calorimeter [0,1]
93  MomHBDepth = cms.double(0.2),
94  MomHEDepth = cms.double(0.4),
95  MomEBDepth = cms.double(0.3),
96  MomEEDepth = cms.double(0.0),
97 
98 # parameters for handling of anomalous cells
99 #
100  # acceptable severity level
101  HcalAcceptSeverityLevel = cms.uint32(9),
102  #EcalAcceptSeverityLevel = cms.uint32(1),
103  EcalRecHitSeveritiesToBeExcluded = cms.vstring('kTime','kWeird','kBad'),
104 
105  # use of recovered hits
106  UseHcalRecoveredHits = cms.bool(True),
107  # The CaloTower code treats recovered cells as a separate category.
108  # The flag to use (or not use) them should be explicitly set
109  # regardless of the specified severity level in EcalAcceptSeverityLevel.
110  UseEcalRecoveredHits = cms.bool(False),
111 
112 
113 # NOTE: The following controls the creation of towers from
114 # rejected rechits.
115 # Always make sure that UseRejectedHitsOnly=false for
116 # normal reconstructions!!!
117 
118  UseRejectedHitsOnly = cms.bool(False),
119 
120 # Controls for hits to be included in the "bad" tower collection.
121 # They have no effect unless UseRejectedHitsOnly=true
122 # Hits passing the HcalAcceptSeverityLevel
123 # will be skipped as they are already in the default towers
124  HcalAcceptSeverityLevelForRejectedHit = cms.uint32(9999),
125 # List of ECAL problems that should be used in bad tower construction
126 # Note that these can only be of type already rejected in default
127 # reconstruction as specified in "EcalRecHitSeveritiesToBeExcluded"
128  EcalSeveritiesToBeUsedInBadTowers = cms.vstring(),
129 
130 
131 # The code also checks the settings of the flags for the default
132 # collection - if the recovered hits were used there, they
133 # will be skipped for the "bad" tower collection regardless of these settings
134  UseRejectedRecoveredHcalHits = cms.bool(True),
135  UseRejectedRecoveredEcalHits = cms.bool(False),
136 
137 
138 # flag to allow/disallow missing inputs
139  AllowMissingInputs = cms.bool(False),
140 
141 # specify hcal upgrade phase - 0, 1, 2
142  HcalPhase = cms.int32(0)
143 
144 )