2 import FWCore.ParameterSet.Config
as cms
10 print "Input process set to", inputProcess
11 process._Process__name=
"EmbeddedRECO"
12 process.TFileService = cms.Service(
"TFileService", fileName = cms.string(
"histo_embedded.root") )
15 outputModule = process.output
19 outputModule = getattr(process,str(getattr(process,
list(process.endpaths)[-1])))
23 print "Changing eventcontent to RAW+AODSIM + misc. "
24 outputModule.outputCommands = cms.untracked.vstring(
"drop *")
25 outputModule.outputCommands.extend(process.RAWEventContent.outputCommands )
26 outputModule.outputCommands.extend(process.AODSIMEventContent.outputCommands )
28 keepMC = cms.untracked.vstring(
"keep *_*_zMusExtracted_*",
29 "keep *_*_zmmCands_*",
30 "keep *_removedInputMuons_*_*",
31 "keep *_generator_*_*",
32 "keep *_PhotonIDProd_*_*",
34 "keep *_photonCore_*_*",
35 "keep *_genParticles_*_*",
36 "keep *_particleFlow_*_*",
37 "keep *_generator_*_*",
38 "keep *_tmfTracks_*_EmbeddedRECO",
39 "keep *_offlinePrimaryVertices_*_EmbeddedRECO",
40 "keep *_offlinePrimaryVerticesWithBS_*_EmbeddedRECO"
42 outputModule.outputCommands.extend(keepMC)
46 for item
in outputModule.outputCommands[:]:
47 if item ==
"drop *" and index != 0:
49 del outputModule.outputCommands[index]
53 hltProcessName =
"HLT"
56 hltProcessName = __HLT__
61 process.dimuonsHLTFilter.TriggerResultsTag.processName = hltProcessName
62 process.goodZToMuMuAtLeast1HLT.TrigTag.processName = hltProcessName
63 process.goodZToMuMuAtLeast1HLT.triggerEvent.processName = hltProcessName
64 process.hltTrigReport,HLTriggerResults.processName = hltProcessName
68 process.VtxSmeared = cms.EDProducer(
"FlatEvtVtxGenerator",
69 MaxZ = cms.double(0.0),
70 MaxX = cms.double(0.0),
71 MaxY = cms.double(0.0),
72 MinX = cms.double(0.0),
73 MinY = cms.double(0.0),
74 MinZ = cms.double(0.0),
75 TimeOffset = cms.double(0.0),
76 src = cms.InputTag(
"generator",
"unsmeared")
79 import FWCore.ParameterSet.VarParsing
as VarParsing
81 options.register (
'mdtau',
83 VarParsing.VarParsing.multiplicity.singleton,
84 VarParsing.VarParsing.varType.int,
85 "mdtau value for tauola")
87 options.register (
'transformationMode',
89 VarParsing.VarParsing.multiplicity.singleton,
90 VarParsing.VarParsing.varType.int,
91 "transformation mode. 0=mumu->mumu, 1=mumu->tautau")
93 options.register (
'minVisibleTransverseMomentum',
95 VarParsing.VarParsing.multiplicity.singleton,
96 VarParsing.VarParsing.varType.string,
97 "generator level cut on visible transverse momentum (typeN:pT,[...];[...])")
99 options.register (
'useJson',
101 VarParsing.VarParsing.multiplicity.singleton,
102 VarParsing.VarParsing.varType.int,
103 "should I enable json usage?")
105 options.register (
'overrideBeamSpot',
107 VarParsing.VarParsing.multiplicity.singleton,
108 VarParsing.VarParsing.varType.int,
109 "should I override beamspot in globaltag?")
118 if not hasattr(process,
"doNotParse"):
120 if hasattr(sys,
"argv") ==
True:
121 if not sys.argv[0].endswith(
'cmsDriver.py'):
122 options.parseArguments()
125 print "CL parsing disabled!"
129 print "Setting mdtau to ", options.mdtau
130 process.generator.ZTauTau.TauolaOptions.InputCards.mdtau = options.mdtau
131 process.newSource.ZTauTau.TauolaOptions.InputCards.mdtau = options.mdtau
132 process.generator.ParticleGun.ExternalDecays.Tauola.InputCards.mdtau = options.mdtau
133 process.newSource.ParticleGun.ExternalDecays.Tauola.InputCards.mdtau = options.mdtau
135 print "Setting minVisibleTransverseMomentum to ", options.minVisibleTransverseMomentum
136 process.newSource.ZTauTau.minVisibleTransverseMomentum = cms.untracked.string(options.minVisibleTransverseMomentum)
137 process.generator.ZTauTau.minVisibleTransverseMomentum = cms.untracked.string(options.minVisibleTransverseMomentum)
139 print "Setting transformationMode to ", options.transformationMode
140 process.generator.ZTauTau.transformationMode = cms.untracked.int32(options.transformationMode)
141 process.newSource.ZTauTau.transformationMode = cms.untracked.int32(options.transformationMode)
144 if hasattr(process,
"changeBSfromPS"):
145 changeBSfromPS =
True
147 if setFromCL
and options.overrideBeamSpot != 0
or changeBSfromPS :
148 print "options.overrideBeamSpot", options.overrideBeamSpot
150 bs = cms.string(
"BeamSpotObjects_2009_LumiBased_SigmaZ_v21_offline")
156 process.GlobalTag.toGet = cms.VPSet(
157 cms.PSet(record = cms.string(
"BeamSpotObjectsRcd"),
159 connect = cms.untracked.string(
"frontier://FrontierProd/CMS_COND_31X_BEAMSPOT")
162 print "BeamSpot in globaltag set to ", bs
164 print "BeamSpot in globaltag not changed"
166 if setFromCL
and options.useJson != 0:
167 print "Enabling json usage"
168 import PhysicsTools.PythonAnalysis.LumiList
as LumiList
169 import FWCore.ParameterSet.Types
as CfgTypes
171 process.source.lumisToProcess = CfgTypes.untracked(CfgTypes.VLuminosityBlockRange())
172 process.source.lumisToProcess.extend(myLumis)
174 process.generalTracksORG = process.generalTracks.clone()
176 process.generalTracks = cms.EDProducer(
"RecoTracksMixer",
177 trackCol1 = cms.InputTag(
"generalTracksORG",
"",
"EmbeddedRECO"),
178 trackCol2 = cms.InputTag(
"removedInputMuons",
"tracks")
181 for p
in process.paths:
182 pth = getattr(process,p)
183 if "generalTracks" in pth.moduleNames():
184 pth.replace(process.generalTracks, process.generalTracksORG*process.generalTracks)
189 process.trackerDrivenElectronSeedsORG = process.trackerDrivenElectronSeeds.clone()
190 process.trackerDrivenElectronSeedsORG.TkColList = cms.VInputTag(cms.InputTag(
"generalTracksORG"))
192 process.trackerDrivenElectronSeeds = cms.EDProducer(
"ElectronSeedTrackRefUpdater",
193 PreIdLabel = process.trackerDrivenElectronSeedsORG.PreIdLabel,
194 PreGsfLabel = process.trackerDrivenElectronSeedsORG.PreGsfLabel,
195 targetTracks = cms.InputTag(
"generalTracks"),
196 inSeeds = cms.InputTag(
"trackerDrivenElectronSeedsORG", process.trackerDrivenElectronSeedsORG.PreGsfLabel.value()),
197 inPreId = cms.InputTag(
"trackerDrivenElectronSeedsORG", process.trackerDrivenElectronSeedsORG.PreIdLabel.value()),
200 for p
in process.paths:
201 pth = getattr(process,p)
202 if "trackerDrivenElectronSeeds" in pth.moduleNames():
203 pth.replace(process.trackerDrivenElectronSeeds, process.trackerDrivenElectronSeedsORG*process.trackerDrivenElectronSeeds)
207 print "TODO: check if photon core hack helps"
209 process.trackerDrivenElectronSeedsMerged = cms.EDProducer("ElectronSeedTrackRefUpdaterAndMerger",
210 PreIdLabel = process.trackerDrivenElectronSeedsORG.PreIdLabel,
211 PreGsfLabel = process.trackerDrivenElectronSeedsORG.PreGsfLabel,
212 targetTracks = cms.InputTag("generalTracks"),
213 inSeeds1 = cms.InputTag("trackerDrivenElectronSeedsORG", process.trackerDrivenElectronSeeds.PreGsfLabel.value()),
214 inPreId1 = cms.InputTag("trackerDrivenElectronSeedsORG", process.trackerDrivenElectronSeeds.PreIdLabel.value()),
215 inSeeds2 = cms.InputTag("trackerDrivenElectronSeeds", process.trackerDrivenElectronSeeds.PreGsfLabel.value(), inputProcess ) ,
216 inPreId2 = cms.InputTag("trackerDrivenElectronSeeds", process.trackerDrivenElectronSeeds.PreIdLabel.value(), inputProcess)
219 process.electronMergedSeedsPhotonCoreHack = cms.EDProducer("ElectronSeedMerger",
220 EcalBasedSeeds = cms.InputTag("ecalDrivenElectronSeeds"),
221 TkBasedSeeds = cms.InputTag("trackerDrivenElectronSeedsMerged","SeedsForGsf")
224 process.photonCore.pixelSeedProducer = cms.string('electronMergedSeedsPhotonCoreHack')
227 for p in process.paths:
228 pth = getattr(process,p)
229 if "photonCore" in pth.moduleNames():
230 pth.replace(process.photonCore,
231 process.trackerDrivenElectronSeedsMerged * process.electronMergedSeedsPhotonCoreHack *process.photonCore)
237 process.electronGsfTracksORG = process.electronGsfTracks.clone()
238 process.electronGsfTracks = cms.EDProducer(
"GsfTrackMixer",
239 col1 = cms.InputTag(
"electronGsfTracksORG",
"",
"EmbeddedRECO"),
240 col2= cms.InputTag(
"electronGsfTracks",
"", inputProcess),
245 process.gsfConversionTrackProducer.TrackProducer = cms.string(
'electronGsfTracksORG')
247 for p
in process.paths:
248 pth = getattr(process,p)
249 if "electronGsfTracks" in pth.moduleNames():
250 pth.replace(process.electronGsfTracks, process.electronGsfTracksORG*process.electronGsfTracks)
253 process.electronMergedSeedsORG = process.electronMergedSeeds.clone()
254 process.electronMergedSeeds = cms.EDProducer("ElectronSeedsMixer",
255 col1 = cms.InputTag("electronMergedSeeds","", inputProcess),
256 col2 = cms.InputTag("electronMergedSeedsORG","","EmbeddedRECO")
258 for p in process.paths:
259 pth = getattr(process,p)
260 if "electronMergedSeeds" in pth.moduleNames():
261 pth.replace(process.electronMergedSeeds, process.electronMergedSeedsORG*process.electronMergedSeeds)
264 process.generalConversionTrackProducer.TrackProducer = cms.string(
'generalTracksORG')
267 process.uncleanedOnlyGeneralConversionTrackProducer.TrackProducer = cms.string(
'generalTracksORG')
270 process.gsfElectronsORG = process.gsfElectrons.clone()
271 process.gsfElectrons = cms.EDProducer(
"GSFElectronsMixer",
272 col1 = cms.InputTag(
"gsfElectronsORG"),
273 col2 = cms.InputTag(
"gsfElectrons",
"",inputProcess),
275 for p
in process.paths:
276 pth = getattr(process,p)
277 if "gsfElectrons" in pth.moduleNames():
278 pth.replace(process.gsfElectrons, process.gsfElectronsORG*process.gsfElectrons)
280 for p
in process.paths:
281 pth = getattr(process,p)
284 for mod
in pth.moduleNames():
285 if mod.find(
"dedx") != -1:
286 print "Removing", mod
287 module=getattr(process,mod)
291 clConfig = cms.PSet (
292 depsPlus = cms.InputTag(
"anaDeposits",
"plus" ),
293 depsMinus = cms.InputTag(
"anaDeposits",
"minus" ),
294 ZmumuCands = cms.InputTag(
"goldenZmumuCandidatesGe2IsoMuons")
297 cleanerConfigTMVA = clConfig.clone()
299 cleanerConfigTMVA.names = cms.vstring(
"H_Ecal_EcalBarrel",
"H_Ecal_EcalEndcap",
300 "H_Hcal_HcalBarrel",
"H_Hcal_HcalOuter",
"H_Hcal_HcalEndcap" )
303 cleanerConfigTMVA.H_Ecal_EcalBarrel = cms.vstring(
"MLP",
"ToyCalo_MLP_ANN_H_Ecal_EcalBarrel.weights.xml" )
304 cleanerConfigTMVA.H_Ecal_EcalEndcap = cms.vstring(
"MLP",
"ToyCalo_MLP_ANN_H_Ecal_EcalEndcap.weights.xml" )
305 cleanerConfigTMVA.H_Hcal_HcalBarrel = cms.vstring(
"MLP",
"ToyCalo_MLP_ANN_H_Hcal_HcalBarrel.weights.xml" )
306 cleanerConfigTMVA.H_Hcal_HcalOuter = cms.vstring(
"MLP",
"ToyCalo_MLP_ANN_H_Hcal_HcalOuter.weights.xml" )
307 cleanerConfigTMVA.H_Hcal_HcalEndcap = cms.vstring(
"MLP",
"ToyCalo_MLP_ANN_H_Hcal_HcalEndcap.weights.xml" )
310 cleanerConfigConst = clConfig.clone()
311 cleanerConfigConst.names = cms.vstring(
"H_Ecal_EcalBarrel",
"H_Ecal_EcalEndcap",
312 "H_Hcal_HcalBarrel",
"H_Hcal_HcalOuter",
"H_Hcal_HcalEndcap" )
313 cleanerConfigConst.H_Ecal_EcalBarrel = cms.double(0.26)
314 cleanerConfigConst.H_Ecal_EcalEndcap = cms.double(1)
315 cleanerConfigConst.H_Hcal_HcalBarrel = cms.double(2)
316 cleanerConfigConst.H_Hcal_HcalOuter = cms.double(1.16)
317 cleanerConfigConst.H_Hcal_HcalEndcap = cms.double(1.82)
319 cleanerConfigConst.H_Ecal_EcalBarrel = cms.double(0.2)
321 cleanerConfigConst.H_Ecal_EcalEndcap = cms.double(0.0)
322 cleanerConfigConst.H_Hcal_HcalBarrel = cms.double(6)
323 cleanerConfigConst.H_Hcal_HcalOuter = cms.double(6)
324 cleanerConfigConst.H_Hcal_HcalEndcap = cms.double(6)
327 process.castorrecoORG = process.castorreco.clone()
328 process.castorreco = cms.EDProducer(
"CastorRecHitMixer",
329 typeEnergyDepositMap = cms.string(
'absolute'),
330 srcEnergyDepositMapMuPlus = cms.InputTag(
"muonCaloEnergyDepositsAllCrossed",
"energyDepositsMuPlus"),
331 srcEnergyDepositMapMuMinus = cms.InputTag(
"muonCaloEnergyDepositsAllCrossed",
"energyDepositsMuMinus"),
332 todo = cms.VPSet(cms.PSet(
333 collection2 = cms.InputTag(
"castorrecoORG"),
334 collection1 = cms.InputTag(
"castorreco",
"",
"RECO")
337 for p
in process.paths:
338 pth = getattr(process,p)
339 if "castorreco" in pth.moduleNames():
340 pth.replace(process.castorreco, process.castorrecoORG*process.castorreco)
342 process.ecalPreshowerRecHitORG = process.ecalPreshowerRecHit.clone()
343 process.ecalPreshowerRecHit = cms.EDProducer(
"EcalRecHitMixer",
344 typeEnergyDepositMap = cms.string(
'absolute'),
345 srcEnergyDepositMapMuPlus = cms.InputTag(
"muonCaloEnergyDepositsAllCrossed",
"energyDepositsMuPlus"),
346 srcEnergyDepositMapMuMinus = cms.InputTag(
"muonCaloEnergyDepositsAllCrossed",
"energyDepositsMuMinus"),
347 todo = cms.VPSet(cms.PSet(
348 collection2 = cms.InputTag(
"ecalPreshowerRecHitORG",
"EcalRecHitsES"),
349 collection1 = cms.InputTag(
"ecalPreshowerRecHit",
"EcalRecHitsES",
"RECO")
352 for p
in process.paths:
353 pth = getattr(process,p)
354 if "ecalPreshowerRecHit" in pth.moduleNames():
355 pth.replace(process.ecalPreshowerRecHit, process.ecalPreshowerRecHitORG*process.ecalPreshowerRecHit)
359 process.ecalRecHitORG = process.ecalRecHit.clone()
360 process.ecalRecHit = cms.EDProducer(
"EcalRecHitMixer",
361 typeEnergyDepositMap = cms.string(
'absolute'),
362 srcEnergyDepositMapMuPlus = cms.InputTag(
"muonCaloEnergyDepositsByDistance",
"energyDepositsMuPlus"),
363 srcEnergyDepositMapMuMinus = cms.InputTag(
"muonCaloEnergyDepositsByDistance",
"energyDepositsMuMinus"),
364 todo = cms.VPSet(cms.PSet(
365 collection2 = cms.InputTag(
"ecalRecHitORG",
"EcalRecHitsEB"),
366 collection1 = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEB",
"RECO")
369 collection2 = cms.InputTag(
"ecalRecHitORG",
"EcalRecHitsEE"),
370 collection1 = cms.InputTag(
"ecalRecHit",
"EcalRecHitsEE",
"RECO")
373 for p
in process.paths:
374 pth = getattr(process,p)
375 if "ecalRecHit" in pth.moduleNames():
376 pth.replace(process.ecalRecHit, process.ecalRecHitORG*process.ecalRecHit)
380 process.hbherecoORG = process.hbhereco.clone()
381 process.hbhereco = cms.EDProducer(
"HBHERecHitMixer",
382 typeEnergyDepositMap = cms.string(
'absolute'),
383 srcEnergyDepositMapMuPlus = cms.InputTag(
"muonCaloEnergyDepositsByDistance",
"energyDepositsMuPlus"),
384 srcEnergyDepositMapMuMinus = cms.InputTag(
"muonCaloEnergyDepositsByDistance",
"energyDepositsMuMinus"),
385 todo = cms.VPSet(cms.PSet(
386 collection2 = cms.InputTag(
"hbherecoORG"),
387 collection1 = cms.InputTag(
"hbhereco",
"",
"RECO")
390 for p
in process.paths:
391 pth = getattr(process,p)
392 if "hbhereco" in pth.moduleNames():
393 pth.replace(process.hbhereco, process.hbherecoORG*process.hbhereco)
395 process.hfrecoORG = process.hfreco.clone()
396 process.hfreco = cms.EDProducer(
"HFRecHitMixer",
397 typeEnergyDepositMap = cms.string(
'absolute'),
398 srcEnergyDepositMapMuPlus = cms.InputTag(
"muonCaloEnergyDepositsAllCrossed",
"energyDepositsMuPlus"),
399 srcEnergyDepositMapMuMinus = cms.InputTag(
"muonCaloEnergyDepositsAllCrossed",
"energyDepositsMuMinus"),
402 collection1 = cms.InputTag(
"hfreco",
"", inputProcess),
403 collection2 = cms.InputTag(
"hfrecoORG")
407 for p
in process.paths:
408 pth = getattr(process, p)
409 if "hfreco" in pth.moduleNames():
410 pth.replace(process.hfreco, process.hfrecoORG*process.hfreco)
412 process.horecoORG = process.horeco.clone()
413 process.horeco = cms.EDProducer(
"HORecHitMixer",
414 typeEnergyDepositMap = cms.string(
'absolute'),
415 srcEnergyDepositMapMuPlus = cms.InputTag(
"muonCaloEnergyDepositsByDistance",
"energyDepositsMuPlus"),
416 srcEnergyDepositMapMuMinus = cms.InputTag(
"muonCaloEnergyDepositsByDistance",
"energyDepositsMuMinus"),
417 todo = cms.VPSet(cms.PSet(
418 collection2 = cms.InputTag(
"horecoORG"),
419 collection1 = cms.InputTag(
"horeco",
"",
"RECO")
422 for p
in process.paths:
423 pth = getattr(process,p)
424 if "horeco" in pth.moduleNames():
425 pth.replace(process.horeco, process.horecoORG*process.horeco)
428 l1ExtraCollections = [
429 [
"L1EmParticle",
"Isolated" ],
430 [
"L1EmParticle",
"NonIsolated" ],
431 [
"L1EtMissParticle",
"MET" ],
432 [
"L1EtMissParticle",
"MHT" ],
433 [
"L1JetParticle",
"Central" ],
434 [
"L1JetParticle",
"Forward" ],
435 [
"L1JetParticle",
"Tau" ],
436 [
"L1MuonParticle",
"" ]
438 l1extraParticleCollections = []
439 for l1ExtraCollection
in l1ExtraCollections:
440 inputType = l1ExtraCollection[0]
442 if inputType ==
"L1EmParticle":
443 pluginType =
"L1ExtraEmParticleMixerPlugin"
444 elif inputType ==
"L1EtMissParticle":
445 pluginType =
"L1ExtraMEtMixerPlugin"
446 elif inputType ==
"L1JetParticle":
447 pluginType =
"L1ExtraJetParticleMixerPlugin"
448 elif inputType ==
"L1MuonParticle":
449 pluginType =
"L1ExtraMuonParticleMixerPlugin"
451 raise ValueError(
"Invalid L1Extra type = %s !!" % inputType)
452 instanceLabel = l1ExtraCollection[1]
453 l1extraParticleCollections.append(cms.PSet(
454 pluginType = cms.string(pluginType),
455 instanceLabel = cms.string(instanceLabel)))
456 process.l1extraParticlesORG = process.l1extraParticles.clone()
457 process.l1extraParticles = cms.EDProducer(
'L1ExtraMixer',
458 src1 = cms.InputTag(
'l1extraParticles::HLT'),
459 src2 = cms.InputTag(
'l1extraParticlesORG'),
460 collections = cms.VPSet(l1extraParticleCollections)
462 for p
in process.paths:
463 pth = getattr(process,p)
464 if "l1extraParticles" in pth.moduleNames():
465 pth.replace(process.l1extraParticles, process.l1extraParticlesORG*process.l1extraParticles)
470 for s
in process.sequences:
471 seq = getattr(process,s)
472 seq.remove(process.offlineBeamSpot)
475 process.metreco.remove(process.BeamHaloId)
480 outputModule = process.output
484 outputModule = getattr(process,str(getattr(process,
list(process.endpaths)[-1])))
488 process.filterEmptyEv.src = cms.untracked.InputTag(
"generatorSmeared",
"",
"EmbeddedRECO")
491 process.schedule.remove(process.DQM_FEDIntegrity_v3)
495 process.RECOSIMoutput.outputCommands.extend([
'keep *_goldenZmumuCandidatesGe0IsoMuons_*_*'])
496 process.RECOSIMoutput.outputCommands.extend([
'keep *_goldenZmumuCandidatesGe1IsoMuons_*_*'])
497 process.RECOSIMoutput.outputCommands.extend([
'keep *_goldenZmumuCandidatesGe2IsoMuons_*_*'])
500 process.RECOSIMoutput.outputCommands.extend([
'keep *_*ORG_*_*'])
510 if hasattr(process,
"doZmumuSkim"):
511 print "Enabling Zmumu skim"
515 cmssw_ver = os.environ[
"CMSSW_VERSION"]
516 if cmssw_ver.find(
"CMSSW_4_2") != -1:
518 print "Using legacy version of Zmumu skim. Note, that muon isolation is disabled"
520 process.load(
"TauAnalysis/MCEmbeddingTools/ZmumuStandaloneSelectionLegacy_cff")
522 process.RandomNumberGeneratorService.dummy = cms.PSet(
523 initialSeed = cms.untracked.uint32(123456789),
524 engineName = cms.untracked.string(
'HepJamesRandom')
528 process.load(
"TauAnalysis/MCEmbeddingTools/ZmumuStandaloneSelection_cff")
529 process.load(
"TrackingTools/TransientTrack/TransientTrackBuilder_cfi")
532 for path
in process.paths:
533 getattr(process,path)._seq = process.goldenZmumuSelectionSequence * getattr(process,path)._seq
535 process.options = cms.untracked.PSet(
536 wantSummary = cms.untracked.bool(
True)
540 print "Zmumu skim not enabled"
542 print "# ######################################################################################"
543 print " Following parameters can be added before customize function "
544 print " call in order to controll process customization: "
546 print "process.doNotParse = cms.PSet() # disables CL parsing for crab compat"
547 print "process.doZmumuSkim = cms.PSet() # adds Zmumu skimming before embedding is run"
548 print "process.changeBSfromPS = cms.PSet() # overide bs"
549 print "# ######################################################################################"
553 process.castorreco.cleaningAlgo = cms.string("CaloCleanerNone")
554 process.ecalPreshowerRecHit.cleaningAlgo = cms.string("CaloCleanerNone")
555 process.ecalRecHit.cleaningAlgo = cms.string("CaloCleanerNone")
556 process.hbhereco.cleaningAlgo = cms.string("CaloCleanerNone")
557 process.hfreco.cleaningAlgo = cms.string("CaloCleanerNone")
558 process.horeco.cleaningAlgo = cms.string("CaloCleanerNone")
562 print "Cleaning set to all crossed!"
563 process.castorreco.cleaningAlgo = cms.string("CaloCleanerAllCrossed")
564 process.ecalPreshowerRecHit.cleaningAlgo = cms.string("CaloCleanerAllCrossed")
565 process.ecalRecHit.cleaningAlgo = cms.string("CaloCleanerAllCrossed")
566 process.hbhereco.cleaningAlgo = cms.string("CaloCleanerAllCrossed")
567 process.hfreco.cleaningAlgo = cms.string("CaloCleanerAllCrossed")
568 process.horeco.cleaningAlgo = cms.string("CaloCleanerAllCrossed")
return(e1-e2)*(e1-e2)+dp *dp
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger list("!*","!HLTx*"if it matches 2 triggers or more) will accept the event if all the matching triggers are FAIL.It will reject the event if any of the triggers are PASS or EXCEPTION(this matches the behavior of"!*"before the partial wildcard feature was incorporated).Triggers which are in the READY state are completely ignored.(READY should never be returned since the trigger paths have been run