10def load_files_for_quirks_scenario(flags):
11 from ExtraParticles.PDGHelpers import getPDGTABLE
12 if getPDGTABLE('PDGTABLE.MeV'):
13 quirk_firststring=1e-6
14 quirk_maxboost=1e-2
15 quirk_maxmerge=1e-6
16 quirk_debugenabled=False
17 quirk_debugdist=1000
18 quirk_debugsteps=1000
19 simdict = flags.Input.SpecialConfiguration
20 quirk_mass = float(simdict["MASS"])
21 quirk_charge = float(simdict["CHARGE"])
22 quirk_pdgid = int(simdict["PDGID"])
23 quirk_stringforce = float(simdict["STRINGFORCE"])
24
25 f = open('PDGTABLE.MeV', 'a')
26 f.write("M%8d %.8E +0.0E+00 -0.0E+00 Quirk +\n" % (quirk_pdgid, quirk_mass))
27 f.write("W%8d 0.E+00 +0.0E+00 -0.0E+00 Quirk +\n" % quirk_pdgid)
28 f.close()
29
30 if flags.Common.ProductionStep == ProductionStep.Simulation:
31 f = open('quirks_setup.txt', 'w')
32 for x in [quirk_mass, quirk_charge, quirk_pdgid, quirk_stringforce, quirk_firststring, quirk_maxboost, quirk_maxmerge, quirk_maxmerge]:
33 f.write(repr(x) + "\n")
34 if quirk_debugenabled:
35 f.write("1\n")
36 f.write(repr(quirk_debugdist) + "\n")
37 f.write(repr(quirk_debugsteps) + "\n")
38 else:
39 f.write("0\n")
40 f.close()
41 del quirk_firststring, quirk_maxboost, quirk_maxmerge, quirk_debugenabled, quirk_debugdist, quirk_debugsteps, doG4SimConfig, simdict, f, quirk_mass, quirk_charge, quirk_pdgid, quirk_stringforce
42
43