CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
isolationCalculator_cfi.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
3 from RecoEgamma.EgammaIsolationAlgos.egammaHBHERecHitThreshold_cff import egammaHBHERecHit
4 
5 isolationSumsCalculator = cms.PSet(
6  #required inputs
7  ComponentName = cms.string('isolationSumsCalculator'),
8 
9  barrelEcalRecHitCollection = cms.InputTag('ecalRecHit:EcalRecHitsEB'),
10  endcapEcalRecHitCollection = cms.InputTag('ecalRecHit:EcalRecHitsEE'),
11 
12  HBHERecHitCollection = egammaHBHERecHit.hbheRecHits,
13  recHitEThresholdHB = egammaHBHERecHit.recHitEThresholdHB,
14  recHitEThresholdHE = egammaHBHERecHit.recHitEThresholdHE,
15  maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity,
16 
17  # Photon will be marked as being near phi module boundary if
18  # it is closer than this. Currently half a crystal.
19  # 1 Ecal Crystal = 0.0174 radians = 1 degree
20  modulePhiBoundary = cms.double(0.0087),
21  # Photon will be marked as being near an eta boundary if
22  # it is between the 0th and 1st element, or the 2nd and 3rd, or the 4th and 5th...
23  moduleEtaBoundary = cms.vdouble(0.0, 0.02, 0.43, 0.46, 0.78, 0.81, 1.13, 1.15, 1.45, 1.58),
24  #What collection of tracks do I use for Track Isolation?
25  trackProducer = cms.InputTag("generalTracks"),
26  #use beam spot for track isolation
27  beamSpotProducer = cms.InputTag("offlineBeamSpot"),
28  #switches, turn on quality cuts for various quantities.
29  vetoClustered = cms.bool(False), #will remove clustered rechits from ecal iso sum
30  useNumCrystals = cms.bool(True), #will define the veto region by number of crystals in stead of geometrically
31  #configuration of parameters for isolation
32 #### BARREL
33  #tracks
34  isolationtrackThresholdA_Barrel = cms.double(0.0),
35  TrackConeOuterRadiusA_Barrel = cms.double(0.4),
36  TrackConeInnerRadiusA_Barrel = cms.double(0.04),
37  isolationtrackEtaSliceA_Barrel = cms.double(0.015),
38  longImpactParameterA_Barrel = cms.double(0.2),
39  transImpactParameterA_Barrel = cms.double(0.1),
40 #
41  isolationtrackThresholdB_Barrel = cms.double(0.0),
42  TrackConeOuterRadiusB_Barrel = cms.double(0.3),
43  TrackConeInnerRadiusB_Barrel = cms.double(0.04),
44  isolationtrackEtaSliceB_Barrel = cms.double(0.015),
45  longImpactParameterB_Barrel = cms.double(0.2),
46  transImpactParameterB_Barrel = cms.double(0.1),
47  #Ecal rechits
48  EcalRecHitInnerRadiusA_Barrel = cms.double(3.5),
49  EcalRecHitOuterRadiusA_Barrel = cms.double(0.4),
50  EcalRecHitEtaSliceA_Barrel = cms.double(2.5),
51  EcalRecHitThreshEA_Barrel = cms.double(0.095),
52  EcalRecHitThreshEtA_Barrel = cms.double(0.0),
53 #
54  EcalRecHitInnerRadiusB_Barrel = cms.double(3.5),
55  EcalRecHitOuterRadiusB_Barrel = cms.double(0.3),
56  EcalRecHitEtaSliceB_Barrel = cms.double(2.5),
57  EcalRecHitThreshEB_Barrel = cms.double(0.095),
58  EcalRecHitThreshEtB_Barrel = cms.double(0.0),
59 
60  # hcal rechits
61  HcalRecHitInnerRadiusA_Barrel = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15),
62  HcalRecHitOuterRadiusA_Barrel = cms.vdouble(0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4),
63 
64  HcalRecHitInnerRadiusB_Barrel = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15),
65  HcalRecHitOuterRadiusB_Barrel = cms.vdouble(0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3),
66 
67 #### ENDCAP
68  #tracks
69  isolationtrackThresholdA_Endcap = cms.double(0.0),
70  TrackConeOuterRadiusA_Endcap = cms.double(0.4),
71  TrackConeInnerRadiusA_Endcap = cms.double(0.04),
72  isolationtrackEtaSliceA_Endcap = cms.double(0.015),
73  longImpactParameterA_Endcap = cms.double(0.2),
74  transImpactParameterA_Endcap = cms.double(0.1),
75 ##
76  isolationtrackThresholdB_Endcap = cms.double(0.0),
77  TrackConeOuterRadiusB_Endcap = cms.double(0.3),
78  TrackConeInnerRadiusB_Endcap = cms.double(0.04),
79  isolationtrackEtaSliceB_Endcap = cms.double(0.015),
80  longImpactParameterB_Endcap = cms.double(0.2),
81  transImpactParameterB_Endcap = cms.double(0.1),
82  #Ecal rechits
83  EcalRecHitInnerRadiusA_Endcap = cms.double(3.5),
84  EcalRecHitOuterRadiusA_Endcap = cms.double(0.4),
85  EcalRecHitEtaSliceA_Endcap = cms.double(2.5),
86  EcalRecHitThreshEA_Endcap = cms.double(0.0),
87  EcalRecHitThreshEtA_Endcap = cms.double(0.110),
88 #
89  EcalRecHitInnerRadiusB_Endcap = cms.double(3.5),
90  EcalRecHitOuterRadiusB_Endcap = cms.double(0.3),
91  EcalRecHitEtaSliceB_Endcap = cms.double(2.5),
92  EcalRecHitThreshEB_Endcap = cms.double(0.0),
93  EcalRecHitThreshEtB_Endcap = cms.double(0.110),
94 
95  #severityLevelCut = cms.int32(4),
96  #severityRecHitThreshold = cms.double(5.0),
97  #spikeIdString = cms.string('kSwissCrossBordersIncluded'),
98  #spikeIdThreshold = cms.double(0.95),
99 
100  # Hcal rechits
101  HcalRecHitInnerRadiusA_Endcap = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15),
102  HcalRecHitOuterRadiusA_Endcap = cms.vdouble(0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4),
103 
104  HcalRecHitInnerRadiusB_Endcap = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15),
105  HcalRecHitOuterRadiusB_Endcap = cms.vdouble(0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3),
106 
107  #recHitFlagsToBeExcluded = cms.vstring(
108  # 'kFaultyHardware',
109  # 'kPoorCalib',
110  # ecalRecHitFlag_kSaturated,
111  # ecalRecHitFlag_kLeadingEdgeRecovered,
112  # ecalRecHitFlag_kNeighboursRecovered,
113  # 'kTowerRecovered',
114  # 'kDead'
115  #),
116 
117 
118 
119 
120 )
121 
122 
Generated for CMSSW Reference Manual by doxygen 1.8.5