CMS 3D CMS Logo

Functions
NeutronBGforMuons_cff Namespace Reference

Functions

def neutronBG (process)
 

Function Documentation

def NeutronBGforMuons_cff.neutronBG (   process)

Definition at line 3 of file NeutronBGforMuons_cff.py.

References mathSSE.return().

Referenced by NeutronBGforMuonsHP_cff.customise(), and NeutronBGforMuonsXS_cff.customise().

3 def neutronBG(process):
4 
5  # common fragment allowing to simulate neutron background in muon system
6 
7  if hasattr(process,'g4SimHits'):
8  # time window 10 second
9  TimeCut = cms.double(10000000000.0)
10  process.common_maximum_time.MaxTrackTime = TimeCut
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.FlagFluo = cms.bool(True)
16  process.g4SimHits.Physics.ThermalNeutrons = cms.untracked.bool(False)
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 = TimeCut
22  process.g4SimHits.StackingAction.DeadRegions = cms.vstring()
23  process.g4SimHits.StackingAction.GammaThreshold = cms.double(0.0)
24  # stepping action
25  process.g4SimHits.SteppingAction.MaxTrackTime = TimeCut
26  process.g4SimHits.SteppingAction.DeadRegions = cms.vstring()
27  # Russian roulette disabled
28  process.g4SimHits.StackingAction.RusRoGammaEnergyLimit = cms.double(0.0)
29  process.g4SimHits.StackingAction.RusRoNeutronEnergyLimit = cms.double(0.0)
30  # Calorimeter hits
31  process.g4SimHits.CaloSD.TmaxHit = TimeCut
32  process.g4SimHits.CaloSD.TmaxHits = cms.vdouble(10000000000,10000000000,10000000000,10000000000,10000000000)
33  # full simulation of HF
34  process.g4SimHits.HCalSD.UseShowerLibrary = cms.bool(False)
35  process.g4SimHits.HFShower.UseShowerLibrary = cms.bool(False)
36 
37  return(process)
38 
return((rh^lh)&mask)