From 723e7c97f9c00cf68d0da76f092722ef57fcb107 Mon Sep 17 00:00:00 2001
From: quli <quli@cern.ch>
Date: Tue, 22 Sep 2020 04:26:13 +0200
Subject: [PATCH] preInclude file of fractionally charged particles
---
.../share/specialConfig/preInclude.fcp.py | 86 +++++++++++++++++++
1 file changed, 86 insertions(+)
create mode 100755 Simulation/SimulationJobOptions/share/specialConfig/preInclude.fcp.py
diff --git a/Simulation/SimulationJobOptions/share/specialConfig/preInclude.fcp.py b/Simulation/SimulationJobOptions/share/specialConfig/preInclude.fcp.py
new file mode 100755
index 00000000000..da3df03bdf8
--- /dev/null
+++ b/Simulation/SimulationJobOptions/share/specialConfig/preInclude.fcp.py
@@ -0,0 +1,86 @@
+#########################################################################
+# preInclude.fcp.py - Wendy Taylor 23 Jun 2011 #
+#########################################################################
+
+def load_files_for_fcp_scenario(MASS, CHARGE, X, Y):
+ import os, shutil, sys
+
+ print X, Y
+ CODE=int(20000000)+int(X)*1000+int(Y)*10
+ print CODE
+
+ ALINE1="M {code} {intmass}.E+03 +0.0E+00 -0.0E+00 fcp +".format(code=CODE,intmass=int(MASS))
+ ALINE2="W {code} 0.E+00 +0.0E+00 -0.0E+00 fcp +".format(code=CODE)
+ BLINE1="{code} {intmass}.00 {fcharge} 0.0 # fcp".format(code=CODE,intmass=int(MASS), fcharge=float(CHARGE))
+ BLINE2="-{code} {intmass}.00 -{fcharge} 0.0 # fcpBar".format(code=CODE,intmass=int(MASS), fcharge=float(CHARGE))
+
+ pdgmod = os.path.isfile('PDGTABLE.MeV')
+ if pdgmod is True:
+ os.remove('PDGTABLE.MeV')
+ os.system('get_files -data PDGTABLE.MeV')
+ f=open('PDGTABLE.MeV','a')
+ f.writelines(str(ALINE1))
+ f.writelines('\n')
+ f.writelines(str(ALINE2))
+ f.writelines('\n')
+ f.close()
+ partmod = os.path.isfile('particles.txt')
+ if partmod is True:
+ os.remove('particles.txt')
+ f=open('particles.txt','w')
+ f.writelines(str(BLINE1))
+ f.writelines('\n')
+ f.writelines(str(BLINE2))
+ f.writelines('\n')
+ f.close()
+
+ del ALINE1
+ del ALINE2
+ del BLINE1
+ del BLINE2
+
+doG4SimConfig = True
+from AthenaCommon.AthenaCommonFlags import athenaCommonFlags
+import PyUtils.AthFile as af
+try:
+ f = af.fopen(athenaCommonFlags.FilesInput()[0])
+
+ if "StreamHITS" in f.infos["stream_names"]:
+ from Digitization.DigitizationFlags import digitizationFlags
+ simdict = digitizationFlags.specialConfiguration.get_Value()
+ doG4SimConfig = False
+ else:
+ from G4AtlasApps.SimFlags import simFlags
+ if not "InteractingPDGCodes" in simFlags.specialConfiguration.get_Value():
+ assert "CHARGE" in simFlags.specialConfiguration.get_Value()
+ assert "X" in simFlags.specialConfiguration.get_Value()
+ assert "Y" in simFlags.specialConfiguration.get_Value()
+ CODE=int(20000000)+int(simFlags.specialConfiguration.get_Value()["X"])*1000+int(simFlags.specialConfiguration.get_Value()["Y"])*10
+ simFlags.specialConfiguration.get_Value()['InteractingPDGCodes'] = str([CODE,-1*CODE])
+ simdict = simFlags.specialConfiguration.get_Value()
+except:
+ from G4AtlasApps.SimFlags import simFlags
+ simdict = simFlags.specialConfiguration.get_Value()
+
+assert "MASS" in simdict
+assert "CHARGE" in simdict
+assert "X" in simdict
+assert "Y" in simdict
+load_files_for_fcp_scenario(simdict["MASS"], simdict["CHARGE"], simdict["X"], simdict["Y"])
+
+if doG4SimConfig:
+ from G4AtlasApps import AtlasG4Eng
+ AtlasG4Eng.G4Eng.log.info("Unlocking simFlags.EquationOfMotion to reset the value for Monopole simulation.")
+ from G4AtlasApps.SimFlags import simFlags
+ # FIXME ideally would include this file early enough, so that the unlocking is not required
+ #simFlags.EquationOfMotion.unlock()
+ #simFlags.EquationOfMotion.set_On()
+ #simFlags.EquationOfMotion.set_Value_and_Lock("G4mplEqMagElectricField")#"MonopoleEquationOfMotion")
+ #simFlags.G4Stepper.set_Value_and_Lock('ClassicalRK4')
+ simFlags.PhysicsOptions += ["MonopolePhysicsTool"]
+ # add monopole-specific configuration for looper killer
+ simFlags.OptionalUserActionList.addAction('G4UA::MonopoleLooperKillerTool')
+ # add default HIP killer
+ simFlags.OptionalUserActionList.addAction('G4UA::HIPKillerTool')
+
+del doG4SimConfig, simdict
--
GitLab