def _initPR(self):
""" Describes the physics regions for LAr H6.
"""
from G4AtlasApps.SimFlags import simFlags
from AthenaCommon.DetFlags import DetFlags
if (simFlags.SimLayout.get_Value()=="tb_LArH6_2003"):
if (DetFlags.FCal_on()):
self.H6_lar_PhysReg=PyG4Atlas.PhysicsReg('FCAL')
self.H6_lar_PhysReg.add_Volumes('LArMgr::LAr::FCAL::LiquidArgonC')
self.H6_lar_PhysReg.add_Cuts('e-',0.03)
self.H6_lar_PhysReg.add_Cuts('e+',0.03)
self.H6_lar_PhysReg.add_Cuts('gamma',0.03)
self.H6_lar.add_PhysicsReg(self.H6_lar_PhysReg)
elif (simFlags.SimLayout.get_Value()=="tb_LArH6_2002"):
if (DetFlags.HEC_on()):
self.H6_lar_PhysReg=PyG4Atlas.PhysicsReg('HEC')
self.H6_lar_PhysReg.add_Volumes('LArMgr::LAr::HEC::LiquidArgon')
self.H6_lar_PhysReg.add_Cuts('e-',0.03)
self.H6_lar_PhysReg.add_Cuts('e+',0.03)
self.H6_lar_PhysReg.add_Cuts('gamma',0.03)
self.H6_lar.add_PhysicsReg(self.H6_lar_PhysReg)
elif (simFlags.SimLayout.get_Value()=="tb_LArH6EC_2002"):
if (DetFlags.em_on()):
self.H6_lar_PhysReg=PyG4Atlas.PhysicsReg('EMEC')
self.H6_lar_PhysReg.add_Volumes('LArMgr::LAr::EMEC::Mother')
self.H6_lar_PhysReg.add_Cuts('e-',0.03)
self.H6_lar_PhysReg.add_Cuts('e+',0.03)
self.H6_lar_PhysReg.add_Cuts('gamma',0.03)
self.H6_lar.add_PhysicsReg(self.H6_lar_PhysReg)
elif (simFlags.SimLayout.get_Value()=="tb_LArH6_2004"):
if (simFlags.LArTB_H6Hec.get_Value()):
self.H6_lar_PhysReg=PyG4Atlas.PhysicsReg('HEC')
self.H6_lar_PhysReg.add_Volumes('LArMgr::LAr::HEC::LiquidArgon')
self.H6_lar_PhysReg.add_Cuts('e-',0.03)
self.H6_lar_PhysReg.add_Cuts('e+',0.03)
self.H6_lar_PhysReg.add_Cuts('gamma',0.03)
self.H6_lar.add_PhysicsReg(self.H6_lar_PhysReg)
if (simFlags.LArTB_H6Emec.get_Value()):
self.H6_lar_PhysReg1=PyG4Atlas.PhysicsReg('EMEC')
self.H6_lar_PhysReg1.add_Volumes('LArMgr::LAr::EMEC::Mother')
self.H6_lar_PhysReg1.add_Cuts('e-',0.03)
self.H6_lar_PhysReg1.add_Cuts('e+',0.03)
self.H6_lar_PhysReg1.add_Cuts('gamma',0.03)
self.H6_lar.add_PhysicsReg(self.H6_lar_PhysReg1)
if (simFlags.LArTB_H6Fcal.get_Value()):
self.H6_lar_PhysReg2=PyG4Atlas.PhysicsReg('FCAL')
self.H6_lar_PhysReg2.add_Volumes('LArMgr::LAr::FCAL::LiquidArgonC')
self.H6_lar_PhysReg2.add_Cuts('e-',0.03)
self.H6_lar_PhysReg2.add_Cuts('e+',0.03)
self.H6_lar_PhysReg2.add_Cuts('gamma',0.03)
self.H6_lar.add_PhysicsReg(self.H6_lar_PhysReg2)
def getTB_RegionCreatorList():
regionCreatorList = []
from G4AtlasApps.SimFlags import simFlags
from AthenaCommon.DetFlags import DetFlags
if (simFlags.SimLayout.get_Value() == "tb_LArH6_2003"):
if (DetFlags.FCal_on()):
regionCreatorList += ['FCALPhysicsRegionTool']
elif (simFlags.SimLayout.get_Value() == "tb_LArH6_2002"):
if (DetFlags.HEC_on()):
regionCreatorList += ['HECPhysicsRegionTool']
elif (simFlags.SimLayout.get_Value() == "tb_LArH6EC_2002"):
if (DetFlags.em_on()):
regionCreatorList += ['EMECPhysicsRegionTool']
elif (simFlags.SimLayout.get_Value() == "tb_LArH6_2004"):
if (simFlags.LArTB_H6Hec.get_Value()):
regionCreatorList += ['HECPhysicsRegionTool']
if (simFlags.LArTB_H6Emec.get_Value()):
regionCreatorList += ['EMECPhysicsRegionTool']
if (simFlags.LArTB_H6Fcal.get_Value()):
regionCreatorList += ['FCALPhysicsRegionTool']
return regionCreatorList
