CMS 3D CMS Logo

isolationCalculator_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 
4 isolationSumsCalculator = cms.PSet(
5  #required inputs
6  ComponentName = cms.string('isolationSumsCalculator'),
7 
8  barrelEcalRecHitCollection = cms.InputTag('ecalRecHit:EcalRecHitsEB'),
9  endcapEcalRecHitCollection = cms.InputTag('ecalRecHit:EcalRecHitsEE'),
10  HcalRecHitCollection = cms.InputTag('towerMaker'),
11 
12  # Photon will be marked as being near phi module boundary if
13  # it is closer than this. Currently half a crystal.
14  # 1 Ecal Crystal = 0.0174 radians = 1 degree
15  modulePhiBoundary = cms.double(0.0087),
16  # Photon will be marked as being near an eta boundary if
17  # it is between the 0th and 1st element, or the 2nd and 3rd, or the 4th and 5th...
18  moduleEtaBoundary = cms.vdouble(0.0, 0.02, 0.43, 0.46, 0.78, 0.81, 1.13, 1.15, 1.45, 1.58),
19  #What collection of tracks do I use for Track Isolation?
20  trackProducer = cms.InputTag("generalTracks"),
21  #use beam spot for track isolation
22  beamSpotProducer = cms.InputTag("offlineBeamSpot"),
23  #switches, turn on quality cuts for various quantities.
24  vetoClustered = cms.bool(False), #will remove clustered rechits from ecal iso sum
25  useNumCrystals = cms.bool(True), #will define the veto region by number of crystals in stead of geometrically
26  #configuration of parameters for isolation
27 #### BARREL
28  #tracks
29  isolationtrackThresholdA_Barrel = cms.double(0.0),
30  TrackConeOuterRadiusA_Barrel = cms.double(0.4),
31  TrackConeInnerRadiusA_Barrel = cms.double(0.04),
32  isolationtrackEtaSliceA_Barrel = cms.double(0.015),
33  longImpactParameterA_Barrel = cms.double(0.2),
34  transImpactParameterA_Barrel = cms.double(0.1),
35 #
36  isolationtrackThresholdB_Barrel = cms.double(0.0),
37  TrackConeOuterRadiusB_Barrel = cms.double(0.3),
38  TrackConeInnerRadiusB_Barrel = cms.double(0.04),
39  isolationtrackEtaSliceB_Barrel = cms.double(0.015),
40  longImpactParameterB_Barrel = cms.double(0.2),
41  transImpactParameterB_Barrel = cms.double(0.1),
42  #Ecal rechits
43  EcalRecHitInnerRadiusA_Barrel = cms.double(3.5),
44  EcalRecHitOuterRadiusA_Barrel = cms.double(0.4),
45  EcalRecHitEtaSliceA_Barrel = cms.double(2.5),
46  EcalRecHitThreshEA_Barrel = cms.double(0.095),
47  EcalRecHitThreshEtA_Barrel = cms.double(0.0),
48 #
49  EcalRecHitInnerRadiusB_Barrel = cms.double(3.5),
50  EcalRecHitOuterRadiusB_Barrel = cms.double(0.3),
51  EcalRecHitEtaSliceB_Barrel = cms.double(2.5),
52  EcalRecHitThreshEB_Barrel = cms.double(0.095),
53  EcalRecHitThreshEtB_Barrel = cms.double(0.0),
54  #Hcal towers
55  HcalTowerInnerRadiusA_Barrel = cms.double(0.15),
56  HcalTowerOuterRadiusA_Barrel = cms.double(0.4),
57  HcalTowerThreshEA_Barrel = cms.double(0.0),
58  HcalDepth1TowerInnerRadiusA_Barrel = cms.double(0.15),
59  HcalDepth1TowerOuterRadiusA_Barrel = cms.double(0.4),
60  HcalDepth1TowerThreshEA_Barrel = cms.double(0.0),
61  HcalDepth2TowerInnerRadiusA_Barrel = cms.double(0.15),
62  HcalDepth2TowerOuterRadiusA_Barrel = cms.double(0.4),
63  HcalDepth2TowerThreshEA_Barrel = cms.double(0.0),
64 #
65  HcalTowerInnerRadiusB_Barrel = cms.double(0.15),
66  HcalTowerOuterRadiusB_Barrel = cms.double(0.3),
67  HcalTowerThreshEB_Barrel = cms.double(0.0),
68  HcalDepth1TowerInnerRadiusB_Barrel = cms.double(0.15),
69  HcalDepth1TowerOuterRadiusB_Barrel = cms.double(0.3),
70  HcalDepth1TowerThreshEB_Barrel = cms.double(0.0),
71  HcalDepth2TowerInnerRadiusB_Barrel = cms.double(0.15),
72  HcalDepth2TowerOuterRadiusB_Barrel = cms.double(0.3),
73  HcalDepth2TowerThreshEB_Barrel = cms.double(0.0),
74 
75 #### ENDCAP
76  #tracks
77  isolationtrackThresholdA_Endcap = cms.double(0.0),
78  TrackConeOuterRadiusA_Endcap = cms.double(0.4),
79  TrackConeInnerRadiusA_Endcap = cms.double(0.04),
80  isolationtrackEtaSliceA_Endcap = cms.double(0.015),
81  longImpactParameterA_Endcap = cms.double(0.2),
82  transImpactParameterA_Endcap = cms.double(0.1),
83 
84  isolationtrackThresholdB_Endcap = cms.double(0.0),
85  TrackConeOuterRadiusB_Endcap = cms.double(0.3),
86  TrackConeInnerRadiusB_Endcap = cms.double(0.04),
87  isolationtrackEtaSliceB_Endcap = cms.double(0.015),
88  longImpactParameterB_Endcap = cms.double(0.2),
89  transImpactParameterB_Endcap = cms.double(0.1),
90  #Ecal rechits
91  EcalRecHitInnerRadiusA_Endcap = cms.double(3.5),
92  EcalRecHitOuterRadiusA_Endcap = cms.double(0.4),
93  EcalRecHitEtaSliceA_Endcap = cms.double(2.5),
94  EcalRecHitThreshEA_Endcap = cms.double(0.0),
95  EcalRecHitThreshEtA_Endcap = cms.double(0.110),
96 #
97  EcalRecHitInnerRadiusB_Endcap = cms.double(3.5),
98  EcalRecHitOuterRadiusB_Endcap = cms.double(0.3),
99  EcalRecHitEtaSliceB_Endcap = cms.double(2.5),
100  EcalRecHitThreshEB_Endcap = cms.double(0.0),
101  EcalRecHitThreshEtB_Endcap = cms.double(0.110),
102 
103  #severityLevelCut = cms.int32(4),
104  #severityRecHitThreshold = cms.double(5.0),
105  #spikeIdString = cms.string('kSwissCrossBordersIncluded'),
106  #spikeIdThreshold = cms.double(0.95),
107 
108  #Hcal towers
109  HcalTowerInnerRadiusA_Endcap = cms.double(0.15),
110  HcalTowerOuterRadiusA_Endcap = cms.double(0.4),
111  HcalTowerThreshEA_Endcap = cms.double(0.0),
112  HcalDepth1TowerInnerRadiusA_Endcap = cms.double(0.15),
113  HcalDepth1TowerOuterRadiusA_Endcap = cms.double(0.4),
114  HcalDepth1TowerThreshEA_Endcap = cms.double(0.0),
115  HcalDepth2TowerInnerRadiusA_Endcap = cms.double(0.15),
116  HcalDepth2TowerOuterRadiusA_Endcap = cms.double(0.4),
117  HcalDepth2TowerThreshEA_Endcap = cms.double(0.0),
118 #
119  HcalTowerInnerRadiusB_Endcap = cms.double(0.15),
120  HcalTowerOuterRadiusB_Endcap = cms.double(0.3),
121  HcalTowerThreshEB_Endcap = cms.double(0.0),
122  HcalDepth1TowerInnerRadiusB_Endcap = cms.double(0.15),
123  HcalDepth1TowerOuterRadiusB_Endcap = cms.double(0.3),
124  HcalDepth1TowerThreshEB_Endcap = cms.double(0.0),
125  HcalDepth2TowerInnerRadiusB_Endcap = cms.double(0.15),
126  HcalDepth2TowerOuterRadiusB_Endcap = cms.double(0.3),
127  HcalDepth2TowerThreshEB_Endcap = cms.double(0.0),
128 
129  #recHitFlagsToBeExcluded = cms.vstring(
130  # 'kFaultyHardware',
131  # 'kPoorCalib',
132  # ecalRecHitFlag_kSaturated,
133  # ecalRecHitFlag_kLeadingEdgeRecovered,
134  # ecalRecHitFlag_kNeighboursRecovered,
135  # 'kTowerRecovered',
136  # 'kDead'
137  #),
138 
139 
140 
141 
142 )
143 
144