/
githubmirror
/
cmssw
Обзор
Документация
Войти
/
githubmirror
/
cmssw
Код
Запросы
0
Пакеты
0
Релизы
0
Аналитика
Безопасность
master
RecoEgamma/PhotonIdentification/python/isolationCalculator_cfi.py
120 строк
5 KB
Afiq Anuar
remove HF and centralized hcal thresholds
08 июн 2021, 14:53
08 июн 2021, 14:53
4c513e7
Код
Авторство
О чём код?
import FWCore.ParameterSet.Config as cms from RecoEgamma.EgammaIsolationAlgos.egammaHBHERecHitThreshold_cff import egammaHBHERecHit isolationSumsCalculator = cms.PSet( #required inputs ComponentName = cms.string('isolationSumsCalculator'), barrelEcalRecHitCollection = cms.InputTag('ecalRecHit:EcalRecHitsEB'), endcapEcalRecHitCollection = cms.InputTag('ecalRecHit:EcalRecHitsEE'), HBHERecHitCollection = egammaHBHERecHit.hbheRecHits, recHitEThresholdHB = egammaHBHERecHit.recHitEThresholdHB, recHitEThresholdHE = egammaHBHERecHit.recHitEThresholdHE, maxHcalRecHitSeverity = egammaHBHERecHit.maxHcalRecHitSeverity, # Photon will be marked as being near phi module boundary if # it is closer than this. Currently half a crystal. # 1 Ecal Crystal = 0.0174 radians = 1 degree modulePhiBoundary = cms.double(0.0087), # Photon will be marked as being near an eta boundary if # it is between the 0th and 1st element, or the 2nd and 3rd, or the 4th and 5th... moduleEtaBoundary = cms.vdouble(0.0, 0.02, 0.43, 0.46, 0.78, 0.81, 1.13, 1.15, 1.45, 1.58), #What collection of tracks do I use for Track Isolation? trackProducer = cms.InputTag("generalTracks"), #use beam spot for track isolation beamSpotProducer = cms.InputTag("offlineBeamSpot"), #switches, turn on quality cuts for various quantities. vetoClustered = cms.bool(False), #will remove clustered rechits from ecal iso sum useNumCrystals = cms.bool(True), #will define the veto region by number of crystals in stead of geometrically #configuration of parameters for isolation #### BARREL #tracks isolationtrackThresholdA_Barrel = cms.double(0.0), TrackConeOuterRadiusA_Barrel = cms.double(0.4), TrackConeInnerRadiusA_Barrel = cms.double(0.04), isolationtrackEtaSliceA_Barrel = cms.double(0.015), longImpactParameterA_Barrel = cms.double(0.2), transImpactParameterA_Barrel = cms.double(0.1), # isolationtrackThresholdB_Barrel = cms.double(0.0), TrackConeOuterRadiusB_Barrel = cms.double(0.3), TrackConeInnerRadiusB_Barrel = cms.double(0.04), isolationtrackEtaSliceB_Barrel = cms.double(0.015), longImpactParameterB_Barrel = cms.double(0.2), transImpactParameterB_Barrel = cms.double(0.1), #Ecal rechits EcalRecHitInnerRadiusA_Barrel = cms.double(3.5), EcalRecHitOuterRadiusA_Barrel = cms.double(0.4), EcalRecHitEtaSliceA_Barrel = cms.double(2.5), EcalRecHitThreshEA_Barrel = cms.double(0.095), EcalRecHitThreshEtA_Barrel = cms.double(0.0), # EcalRecHitInnerRadiusB_Barrel = cms.double(3.5), EcalRecHitOuterRadiusB_Barrel = cms.double(0.3), EcalRecHitEtaSliceB_Barrel = cms.double(2.5), EcalRecHitThreshEB_Barrel = cms.double(0.095), EcalRecHitThreshEtB_Barrel = cms.double(0.0), # hcal rechits HcalRecHitInnerRadiusA_Barrel = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15), HcalRecHitOuterRadiusA_Barrel = cms.vdouble(0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4), HcalRecHitInnerRadiusB_Barrel = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15), HcalRecHitOuterRadiusB_Barrel = cms.vdouble(0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3), #### ENDCAP #tracks isolationtrackThresholdA_Endcap = cms.double(0.0), TrackConeOuterRadiusA_Endcap = cms.double(0.4), TrackConeInnerRadiusA_Endcap = cms.double(0.04), isolationtrackEtaSliceA_Endcap = cms.double(0.015), longImpactParameterA_Endcap = cms.double(0.2), transImpactParameterA_Endcap = cms.double(0.1), ## isolationtrackThresholdB_Endcap = cms.double(0.0), TrackConeOuterRadiusB_Endcap = cms.double(0.3), TrackConeInnerRadiusB_Endcap = cms.double(0.04), isolationtrackEtaSliceB_Endcap = cms.double(0.015), longImpactParameterB_Endcap = cms.double(0.2), transImpactParameterB_Endcap = cms.double(0.1), #Ecal rechits EcalRecHitInnerRadiusA_Endcap = cms.double(3.5), EcalRecHitOuterRadiusA_Endcap = cms.double(0.4), EcalRecHitEtaSliceA_Endcap = cms.double(2.5), EcalRecHitThreshEA_Endcap = cms.double(0.0), EcalRecHitThreshEtA_Endcap = cms.double(0.110), # EcalRecHitInnerRadiusB_Endcap = cms.double(3.5), EcalRecHitOuterRadiusB_Endcap = cms.double(0.3), EcalRecHitEtaSliceB_Endcap = cms.double(2.5), EcalRecHitThreshEB_Endcap = cms.double(0.0), EcalRecHitThreshEtB_Endcap = cms.double(0.110), #severityLevelCut = cms.int32(4), #severityRecHitThreshold = cms.double(5.0), #spikeIdString = cms.string('kSwissCrossBordersIncluded'), #spikeIdThreshold = cms.double(0.95), # Hcal rechits HcalRecHitInnerRadiusA_Endcap = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15), HcalRecHitOuterRadiusA_Endcap = cms.vdouble(0.4, 0.4, 0.4, 0.4, 0.4, 0.4, 0.4), HcalRecHitInnerRadiusB_Endcap = cms.vdouble(0.15, 0.15, 0.15, 0.15, 0.15, 0.15, 0.15), HcalRecHitOuterRadiusB_Endcap = cms.vdouble(0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3), #recHitFlagsToBeExcluded = cms.vstring( # 'kFaultyHardware', # 'kPoorCalib', # ecalRecHitFlag_kSaturated, # ecalRecHitFlag_kLeadingEdgeRecovered, # ecalRecHitFlag_kNeighboursRecovered, # 'kTowerRecovered', # 'kDead' #), )