PhotonIdentification/python/isolationCalculator_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from RecoEgamma.EgammaIsolationAlgos.egammaHBHERecHitThreshold_cff import egammaHBHERecHit
0004
0005 isolationSumsCalculator = cms.PSet(
0006     #required inputs
0007     ComponentName = cms.string('isolationSumsCalculator'),
0008
0009     barrelEcalRecHitCollection = cms.InputTag('ecalRecHit:EcalRecHitsEB'),
0010     endcapEcalRecHitCollection = cms.InputTag('ecalRecHit:EcalRecHitsEE'),
0011
0012     HBHERecHitCollection = egammaHBHERecHit.hbheRecHits,
0013     recHitEThresholdHB = egammaHBHERecHit.recHitEThresholdHB,
0014     recHitEThresholdHE = egammaHBHERecHit.recHitEThresholdHE,
0015     maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity,
0016
0017     # Photon will be marked as being near phi module boundary if
0018     #  it is closer than this.  Currently half a crystal.
0019     #  1 Ecal Crystal = 0.0174 radians = 1 degree
0020     modulePhiBoundary =   cms.double(0.0087),
0021     # Photon will be marked as being near an eta boundary if
0022     #  it is between the 0th and 1st element, or the 2nd and 3rd, or the 4th and 5th...
0023     moduleEtaBoundary = cms.vdouble(0.0, 0.02, 0.43, 0.46, 0.78, 0.81, 1.13, 1.15, 1.45, 1.58),
0024     #What collection of tracks do I use for Track Isolation?
0025     trackProducer = cms.InputTag("generalTracks"),
0026     #use beam spot for track isolation
0027     beamSpotProducer = cms.InputTag("offlineBeamSpot"),
0028     #switches, turn on quality cuts for various quantities.
0029     vetoClustered  = cms.bool(False),  #will remove clustered rechits from ecal iso sum
0030     useNumCrystals = cms.bool(True),  #will define the veto region by number of crystals in stead of geometrically
0031     #configuration of parameters for isolation
0032 #### BARREL
0033     #tracks
0034     isolationtrackThresholdA_Barrel = cms.double(0.0),
0035     TrackConeOuterRadiusA_Barrel    = cms.double(0.4),
0036     TrackConeInnerRadiusA_Barrel    = cms.double(0.04),
0037     isolationtrackEtaSliceA_Barrel  = cms.double(0.015),
0038     longImpactParameterA_Barrel     = cms.double(0.2),
0039     transImpactParameterA_Barrel    = cms.double(0.1),
0040 #
0041     isolationtrackThresholdB_Barrel = cms.double(0.0),
0042     TrackConeOuterRadiusB_Barrel    = cms.double(0.3),
0043     TrackConeInnerRadiusB_Barrel    = cms.double(0.04),
0044     isolationtrackEtaSliceB_Barrel  = cms.double(0.015),
0045     longImpactParameterB_Barrel     = cms.double(0.2),
0046     transImpactParameterB_Barrel    = cms.double(0.1),
0047     #Ecal rechits
0048     EcalRecHitInnerRadiusA_Barrel   = cms.double(3.5),
0049     EcalRecHitOuterRadiusA_Barrel   = cms.double(0.4),
0050     EcalRecHitEtaSliceA_Barrel      = cms.double(2.5),
0051     EcalRecHitThreshEA_Barrel       = cms.double(0.095),
0052     EcalRecHitThreshEtA_Barrel      = cms.double(0.0),
0053 #
0054     EcalRecHitInnerRadiusB_Barrel   = cms.double(3.5),
0055     EcalRecHitOuterRadiusB_Barrel   = cms.double(0.3),
0056     EcalRecHitEtaSliceB_Barrel      = cms.double(2.5),
0057     EcalRecHitThreshEB_Barrel       = cms.double(0.095),
0058     EcalRecHitThreshEtB_Barrel      = cms.double(0.0),
0059
0060     # hcal rechits
0061     HcalRecHitInnerRadiusA_Barrel = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15),
0062     HcalRecHitOuterRadiusA_Barrel = cms.vdouble(0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4),
0063
0064     HcalRecHitInnerRadiusB_Barrel = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15),
0065     HcalRecHitOuterRadiusB_Barrel = cms.vdouble(0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3),
0066
0067 #### ENDCAP
0068     #tracks
0069     isolationtrackThresholdA_Endcap  = cms.double(0.0),
0070     TrackConeOuterRadiusA_Endcap     = cms.double(0.4),
0071     TrackConeInnerRadiusA_Endcap     = cms.double(0.04),
0072     isolationtrackEtaSliceA_Endcap   = cms.double(0.015),
0073     longImpactParameterA_Endcap      = cms.double(0.2),
0074     transImpactParameterA_Endcap     = cms.double(0.1),
0075 ##
0076     isolationtrackThresholdB_Endcap  = cms.double(0.0),
0077     TrackConeOuterRadiusB_Endcap     = cms.double(0.3),
0078     TrackConeInnerRadiusB_Endcap     = cms.double(0.04),
0079     isolationtrackEtaSliceB_Endcap   = cms.double(0.015),
0080     longImpactParameterB_Endcap      = cms.double(0.2),
0081     transImpactParameterB_Endcap     = cms.double(0.1),
0082     #Ecal rechits
0083     EcalRecHitInnerRadiusA_Endcap    = cms.double(3.5),
0084     EcalRecHitOuterRadiusA_Endcap    = cms.double(0.4),
0085     EcalRecHitEtaSliceA_Endcap       = cms.double(2.5),
0086     EcalRecHitThreshEA_Endcap        = cms.double(0.0),
0087     EcalRecHitThreshEtA_Endcap       = cms.double(0.110),
0088 #
0089     EcalRecHitInnerRadiusB_Endcap    = cms.double(3.5),
0090     EcalRecHitOuterRadiusB_Endcap    = cms.double(0.3),
0091     EcalRecHitEtaSliceB_Endcap       = cms.double(2.5),
0092     EcalRecHitThreshEB_Endcap        = cms.double(0.0),
0093     EcalRecHitThreshEtB_Endcap       = cms.double(0.110),
0094
0095     #severityLevelCut = cms.int32(4),
0096     #severityRecHitThreshold = cms.double(5.0),
0097     #spikeIdString = cms.string('kSwissCrossBordersIncluded'),
0098     #spikeIdThreshold = cms.double(0.95),
0099
0100     # Hcal rechits
0101     HcalRecHitInnerRadiusA_Endcap = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15),
0102     HcalRecHitOuterRadiusA_Endcap = cms.vdouble(0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4),
0103
0104     HcalRecHitInnerRadiusB_Endcap = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15),
0105     HcalRecHitOuterRadiusB_Endcap = cms.vdouble(0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3),
0106
0107     #recHitFlagsToBeExcluded = cms.vstring(
0108     #    'kFaultyHardware',
0109     #    'kPoorCalib',
0110     #    ecalRecHitFlag_kSaturated,
0111     #    ecalRecHitFlag_kLeadingEdgeRecovered,
0112     #    ecalRecHitFlag_kNeighboursRecovered,
0113     #    'kTowerRecovered',
0114     #    'kDead'
0115     #),
0116
0117
0118
0119
0120 )
0121
0122