test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HCAL.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 def customise(process):
3 
4 # extend the particle gun acceptance
5 
6  process.generator.AddAntiParticle = cms.bool(False)
7 
8 # no magnetic field
9 
10  process.g4SimHits.UseMagneticField = cms.bool(False)
11  process.UniformMagneticFieldESProducer = cms.ESProducer("UniformMagneticFieldESProducer",
12  ZFieldInTesla = cms.double(0.0)
13  )
14 
15  process.prefer("UniformMagneticFieldESProducer")
16 
17 # modify the content
18 
19  process.output.outputCommands.append("keep *_simHcalUnsuppressedDigis_*_*")
20 
21 # user schedule: use only calorimeters digitization and local reconstruction
22 
23  del process.schedule[:]
24 
25  process.schedule.append(process.generation_step)
26  process.schedule.append(process.simulation_step)
27 
28  delattr(process,"hbhereco")
29  process.hbhereco = process.hbheprereco.clone()
30  process.hcalLocalRecoSequence.replace(process.hbheprereco,process.hbhereco)
31  process.hbhereco.digiLabel = cms.InputTag("simHcalUnsuppressedDigis")
32  process.horeco.digiLabel = cms.InputTag("simHcalUnsuppressedDigis")
33  process.hfreco.digiLabel = cms.InputTag("simHcalUnsuppressedDigis")
34 
35  process.local_digireco = cms.Path(process.mix * process.hcalDigiSequence * process.hbhereco * process.hfreco * process.horeco )
36 
37  process.schedule.append(process.local_digireco)
38 
39  process.load("Validation/Configuration/hcalSimValid_cff")
40 
41  process.AllRecHitsValidation.ecalselector = cms.untracked.string('no')
42 
43  process.local_validation = cms.Path(process.hcalSimHitStudy+process.hcalDigisValidationSequence+process.hcalRecHitsValidationSequence)
44  process.schedule.append(process.local_validation)
45 
46  process.schedule.append(process.endjob_step)
47  process.schedule.append(process.out_step)
48 
49  return(process)
def customise
Definition: HCAL.py:2
return(e1-e2)*(e1-e2)+dp *dp