CMS 3D CMS Logo

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

Functions

def customise
 

Function Documentation

def NeutronBGforMuonsXS_cff.customise (   process)

Definition at line 3 of file NeutronBGforMuonsXS_cff.py.

3 
4 def customise(process):
5 
6  # fragment allowing to simulate neutron background in muon system
7 
8  if hasattr(process,'g4SimHits'):
9  # time window 100 millisecond
10  process.common_maximum_time.MaxTrackTime = cms.double(100000000.0)
11  process.common_maximum_time.DeadRegions = cms.vstring()
12  # Physics List XS
13  process.g4SimHits.Physics.type = cms.string('SimG4Core/Physics/FTFP_BERT_XS_EML')
14  process.g4SimHits.Physics.CutsOnProton = cms.untracked.bool(False)
15  process.g4SimHits.Physics.FlagMuNucl = cms.bool(True)
16  process.g4SimHits.Physics.FlagFluo = cms.bool(True)
17  # Eta cut
18  process.g4SimHits.Generator.MinEtaCut = cms.double(-7.0)
19  process.g4SimHits.Generator.MaxEtaCut = cms.double(7.0)
20  # stacking action
21  process.g4SimHits.StackingAction.MaxTrackTime = cms.double(100000000.0)
22  process.g4SimHits.StackingAction.DeadRegions = cms.vstring()
23  process.g4SimHits.StackingAction.KillHeavy = cms.bool(True)
24  process.g4SimHits.StackingAction.IonThreshold = cms.double(0.001)
25  process.g4SimHits.StackingAction.ProtonThreshold = cms.double(0.001)
26  process.g4SimHits.StackingAction.NeutronThreshold = cms.double(0.0)
27  process.g4SimHits.StackingAction.GammaThreshold = cms.double(0.0005)
28  # stepping action
29  process.g4SimHits.SteppingAction.MaxTrackTime = cms.double(100000000.0)
30  process.g4SimHits.SteppingAction.DeadRegions = cms.vstring()
31  # Russian roulette disabled
32  process.g4SimHits.StackingAction.RusRoGammaEnergyLimit = cms.double(0.0)
33  process.g4SimHits.StackingAction.RusRoNeutronEnergyLimit = cms.double(0.0)
34  # full simulation of HF
35  process.g4SimHits.HFShower.UseHFGflash = cms.bool(False)
36 
37  return(process)
return(e1-e2)*(e1-e2)+dp *dp