test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Functions
customise_PPwithHI Namespace Reference

Functions

def addHIIsolationProducer
 
def customisePPrecoForPeripheralPbPb
 
def customisePPrecoforPPb
 
def modifyClusterLimits
 
def storeCaloTowersAOD
 

Function Documentation

def customise_PPwithHI.addHIIsolationProducer (   process)

Definition at line 4 of file customise_PPwithHI.py.

Referenced by customisePPrecoForPeripheralPbPb(), and customisePPrecoforPPb().

4 
5 def addHIIsolationProducer(process):
6 
7  process.load('Configuration.EventContent.EventContent_cff')
8 
9  # extend RecoEgammaFEVT content
10  process.RecoEgammaFEVT.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
11  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*'
12  ])
13 
14  # extend RecoEgammaRECO content
15  process.RECOEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
16  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
17  'keep recoCaloClusters_islandBasicClusters_*_*'
18  ])
19 
20  process.FEVTEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
21  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
22  'keep recoCaloClusters_islandBasicClusters_*_*'
23  ])
24  process.FEVTSIMEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
25  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
26  'keep recoCaloClusters_islandBasicClusters_*_*'
27  ])
28  # extend RecoEgammaRECO content
29  process.RAWRECOEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
30  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
31  'keep recoCaloClusters_islandBasicClusters_*_*'
32  ])
33 
34  process.RECOSIMEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
35  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
36  'keep recoCaloClusters_islandBasicClusters_*_*'
37  ])
38 
39  process.RAWRECOSIMHLTEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
40  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
41  'keep recoCaloClusters_islandBasicClusters_*_*'
42  ])
43 
44  process.RAWRECODEBUGHLTEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
45  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
46  'keep recoCaloClusters_islandBasicClusters_*_*'
47  ])
48 
49  process.FEVTHLTALLEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
50  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
51  'keep recoCaloClusters_islandBasicClusters_*_*'
52  ])
53 
54  process.FEVTDEBUGEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
55  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*',
56  'keep recoCaloClusters_islandBasicClusters_*_*'
57  ])
58 
59  # extend RecoEgammaAOD content
60  process.AODEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
61  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*'
62  ])
63 
64  process.AODSIMEventContent.outputCommands.extend(['keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerppGED_*_*',
65  'keep recoHIPhotonIsolationedmValueMap_photonIsolationHIProducerpp_*_*'
66  ])
67 
68  # add HI Photon isolation sequence to pp RECO
69  process.load('RecoHI.HiEgammaAlgos.photonIsolationHIProducer_cfi')
70  process.load('RecoEcal.EgammaClusterProducers.islandBasicClusters_cfi')
71 
72  process.photonIsolationHISequencePP = cms.Sequence(process.islandBasicClusters
73  * process.photonIsolationHIProducerpp
74  * process.photonIsolationHIProducerppGED)
75 
76  process.reconstruction *= process.photonIsolationHISequencePP
77 
78  return process
79 
80 
# modify cluster limits to run pp reconstruction on peripheral PbPb
def customise_PPwithHI.customisePPrecoForPeripheralPbPb (   process)

Definition at line 147 of file customise_PPwithHI.py.

References addHIIsolationProducer(), modifyClusterLimits(), and storeCaloTowersAOD().

149 
150  process=addHIIsolationProducer(process)
151  process=modifyClusterLimits(process)
152  process=storeCaloTowersAOD(process)
153 
154  return process
155 
def customise_PPwithHI.customisePPrecoforPPb (   process)

Definition at line 140 of file customise_PPwithHI.py.

References addHIIsolationProducer(), and storeCaloTowersAOD().

141 def customisePPrecoforPPb(process):
142 
143  process=addHIIsolationProducer(process)
144  process=storeCaloTowersAOD(process)
145 
146  return process
def customise_PPwithHI.modifyClusterLimits (   process)

Definition at line 81 of file customise_PPwithHI.py.

Referenced by customisePPrecoForPeripheralPbPb().

81 
82 def modifyClusterLimits(process):
83 
84  hiClusterCut = cms.string("strip < 400000 && pixel < 40000 && (strip < 60000 + 7.0*pixel) && (pixel < 8000 + 0.14*strip)")
85 
86  if hasattr(process,'initialStepSeedsPreSplitting'): process.initialStepSeedsPreSplitting.ClusterCheckPSet.cut = hiClusterCut
87  if hasattr(process,'initialStepSeeds'): process.initialStepSeeds.ClusterCheckPSet.cut = hiClusterCut
88  if hasattr(process,'lowPtTripletStepSeeds'): process.lowPtTripletStepSeeds.ClusterCheckPSet.cut = hiClusterCut
89  if hasattr(process,'globalSeedsFromTriplets'): process.globalSeedsFromTriplets.ClusterCheckPSet.cut = hiClusterCut
90  if hasattr(process,'detachedTripletStepSeeds'): process.detachedTripletStepSeeds.ClusterCheckPSet.cut = hiClusterCut
91  if hasattr(process,'mixedTripletStepSeedsA'): process.mixedTripletStepSeedsA.ClusterCheckPSet.cut = hiClusterCut
92  if hasattr(process,'mixedTripletStepSeedsB'): process.mixedTripletStepSeedsB.ClusterCheckPSet.cut = hiClusterCut
93  if hasattr(process,'globalMixedSeeds'): process.globalMixedSeeds.ClusterCheckPSet.cut = hiClusterCut
94  if hasattr(process,'pixelLessStepSeeds'): process.pixelLessStepSeeds.ClusterCheckPSet.cut = hiClusterCut
95  if hasattr(process,'globalPixelLessSeeds'): process.globalPixelLessSeeds.ClusterCheckPSet.cut = hiClusterCut
96  if hasattr(process,'globalPixelSeeds'): process.globalPixelSeeds.ClusterCheckPSet.cut = hiClusterCut
97  if hasattr(process,'pixelPairStepSeeds'): process.pixelPairStepSeeds.ClusterCheckPSet.cut = hiClusterCut
98  if hasattr(process,'globalSeedsFromPairsWithVertices'): process.globalSeedsFromPairsWithVertices.ClusterCheckPSet.cut = hiClusterCut
99  if hasattr(process,'tobTecStepSeedsPair'): process.tobTecStepSeedsPair.ClusterCheckPSet.cut = hiClusterCut
100  if hasattr(process,'tobTecStepSeedsTripl'): process.tobTecStepSeedsTripl.ClusterCheckPSet.cut = hiClusterCut
101  if hasattr(process,'pixelPairElectronSeeds'): process.pixelPairElectronSeeds.ClusterCheckPSet.cut = hiClusterCut
102  if hasattr(process,'regionalCosmicTrackerSeeds'): process.regionalCosmicTrackerSeeds.ClusterCheckPSet.cut = hiClusterCut
103  if hasattr(process,'stripPairElectronSeeds'): process.stripPairElectronSeeds.ClusterCheckPSet.cut = hiClusterCut
104  if hasattr(process,'photonConvTrajSeedFromSingleLeg'): process.photonConvTrajSeedFromSingleLeg.ClusterCheckPSet.cut = hiClusterCut
105  if hasattr(process,'photonConvTrajSeedFromQuadruplets'): process.photonConvTrajSeedFromQuadruplets.ClusterCheckPSet.cut = hiClusterCut
106  if hasattr(process,'tripletElectronSeeds'): process.tripletElectronSeeds.ClusterCheckPSet.cut = hiClusterCut
107  if hasattr(process,'jetCoreRegionalStepSeeds'): process.jetCoreRegionalStepSeeds.ClusterCheckPSet.cut = hiClusterCut
108 
109 
110  maxElement = cms.uint32(1000000)
111 
112  if hasattr(process,'initialStepSeedsPreSplitting'): process.initialStepSeedsPreSplitting.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = maxElement
113  if hasattr(process,'initialStepSeeds'): process.initialStepSeeds.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = maxElement
114  if hasattr(process,'lowPtTripletStepSeeds'): process.lowPtTripletStepSeeds.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = maxElement
115  if hasattr(process,'mixedTripletStepSeedsA'): process.mixedTripletStepSeedsA.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = maxElement
116  if hasattr(process,'mixedTripletStepSeedsB'): process.mixedTripletStepSeedsB.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = maxElement
117  if hasattr(process,'detachedTripletStepSeeds'): process.detachedTripletStepSeeds.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = maxElement
118  if hasattr(process,'pixelLessStepSeeds'): process.pixelLessStepSeeds.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = maxElement
119  if hasattr(process,'tobTecStepSeedsTripl'): process.tobTecStepSeedsTripl.OrderedHitsFactoryPSet.GeneratorPSet.maxElement = maxElement
120  if hasattr(process,'tobTecStepSeedsPair'): process.tobTecStepSeedsPair.OrderedHitsFactoryPSet.maxElement = maxElement
121  if hasattr(process,'pixelPairStepSeeds'): process.pixelPairStepSeeds.OrderedHitsFactoryPSet.maxElement = maxElement
122  if hasattr(process,'jetCoreRegionalStepSeeds'): process.jetCoreRegionalStepSeeds.OrderedHitsFactoryPSet.maxElement = maxElement
123 
124  return process
125 
126 
# Add caloTowers to AOD event content
def customise_PPwithHI.storeCaloTowersAOD (   process)

Definition at line 127 of file customise_PPwithHI.py.

Referenced by customisePPrecoForPeripheralPbPb(), and customisePPrecoforPPb().

128 def storeCaloTowersAOD(process):
129 
130  process.load('Configuration.EventContent.EventContent_cff')
131 
132  # extend AOD content
133  if hasattr(process,'AODoutput'):
134  process.AODoutput.outputCommands.extend(['keep *_towerMaker_*_*'])
135 
136  if hasattr(process,'AODSIMoutput'):
137  process.AODSIMoutput.outputCommands.extend(['keep *_towerMaker_*_*'])
138 
139  return process