Skip to content

Commit 79b8494

Browse files
mpucciompuccio
committed
Add Omega_c decay table and G4 config to use external decayer
Co-authored with Fabrizio Grosa <fabrizio.grosa@cern.ch>
1 parent 52f7c9a commit 79b8494

File tree

2 files changed

+196
-0
lines changed

2 files changed

+196
-0
lines changed

‎MC/config/PWGHF/pythia8/decayer/force_hadronic_Omegac_forceOmegaChannels.cfg‎

Lines changed: 128 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,128 @@
1+
Init:showChangedParticleData = on
2+
3+
4332:mayDecay on # OmegaC
4+
4332:tau0 = 0.0803
5+
4332:onMode = off
6+
4332:addChannel = 1 0.5217000 0 211 3334
7+
4332:addChannel = 1 0.2306000 0 213 3334 211 111 22 22
8+
4332:addChannel = 1 0.0232000 0 111 211 3334 22 22
9+
4332:addChannel = 1 0.0107000 0 223 211 3334 211 -211 111 22 22
10+
4332:addChannel = 1 0.0368000 0 -13 14 3334
11+
4332:addChannel = 1 0.0025000 0 321 -311 3334 130
12+
4332:addChannel = 1 0.0078000 0 113 211 3334 211 -211
13+
4332:addChannel = 1 0.0100000 0 211 113 3334 211 -211
14+
4332:addChannel = 1 0.0018000 0 211 221 3334 211 -211 111 22 22
15+
4332:addChannel = 1 0.0294000 0 211 111 3334 22 22
16+
4332:addChannel = 1 0.0015000 0 211 111 111 3334 22 22 22 22
17+
4332:addChannel = 1 0.0377000 0 -11 12 3334
18+
4332:addChannel = 1 0.0001000 0 211 -211 211 111 3334 22 22
19+
4332:addChannel = 1 0.0024000 0 213 3334 211 111 22 11 -11
20+
4332:addChannel = 1 0.0026000 0 221 211 3334 211 -211 111 22 22
21+
4332:addChannel = 1 0.0008000 0 321 -311 3334 310 111 111 22 22 22 22
22+
4332:addChannel = 1 0.0160000 0 211 223 3334 211 -211 111 22 22
23+
4332:addChannel = 1 0.0017000 0 223 211 3334 111 22 22 22
24+
4332:addChannel = 1 0.0037000 0 211 221 3334 22 22
25+
4332:addChannel = 1 0.0114000 0 213 111 3334 211 111 22 22 22 22
26+
4332:addChannel = 1 0.0093000 0 111 213 3334 22 22 211 111 22 22
27+
4332:addChannel = 1 0.0033000 0 221 211 3334 22 22
28+
4332:addChannel = 1 0.0040000 0 221 211 3334 111 111 111 22 22 22 22 22 22
29+
4332:addChannel = 1 0.0002000 0 213 111 111 3334 211 111 22 22 22 22 22 22
30+
4332:addChannel = 1 0.0004000 0 223 211 3334 211 -211 111 22 11 -11
31+
4332:addChannel = 1 0.0009000 0 -311 321 3334 310 111 111 22 22 22 22
32+
4332:addChannel = 1 0.0001000 0 221 211 3334 111 111 111 22 22 22 11 -11 22 22
33+
4332:addChannel = 1 0.0022000 0 -311 321 3334 130
34+
4332:addChannel = 1 0.0017000 0 211 211 -211 3334
35+
4332:addChannel = 1 0.0034000 0 211 221 3334 111 111 111 22 22 22 22 22 22
36+
4332:addChannel = 1 0.0004000 0 211 213 -211 3334 211 111 22 22
37+
4332:addChannel = 1 0.0001000 0 211 111 113 3334 22 22 211 -211
38+
4332:addChannel = 1 0.0022000 0 -311 321 3334 310 211 -211
39+
4332:addChannel = 1 0.0005000 0 213 211 -211 3334 211 111 22 22
40+
4332:addChannel = 1 0.0010000 0 111 211 111 3334 22 22 22 22
41+
4332:addChannel = 1 0.0009000 0 321 3334 -311 310 111 111 22 22 22 22
42+
4332:addChannel = 1 0.0017000 0 321 -311 3334 310 211 -211
43+
4332:addChannel = 1 0.0007000 0 111 111 211 3334 22 22 22 22
44+
4332:addChannel = 1 0.0022000 0 321 3334 -311 310 211 -211
45+
4332:addChannel = 1 0.0030000 0 321 3334 -311 130
46+
4332:addChannel = 1 0.0004000 0 211 113 111 3334 211 -211 22 22
47+
4332:addChannel = 1 0.0003000 0 221 211 3334 211 -211 22
48+
4332:addChannel = 1 0.0009000 0 211 111 3334 22 11 -11
49+
4332:addChannel = 1 0.0002000 0 221 111 211 3334 211 -211 22 22 22
50+
4332:addChannel = 1 0.5217000 0 211 3334
51+
4332:addChannel = 1 0.2306000 0 213 3334 211 111 22 22
52+
4332:addChannel = 1 0.0232000 0 111 211 3334 22 22
53+
4332:addChannel = 1 0.0107000 0 223 211 3334 211 -211 111 22 22
54+
4332:addChannel = 1 0.0368000 0 -13 14 3334
55+
4332:addChannel = 1 0.0025000 0 321 -311 3334 130
56+
4332:addChannel = 1 0.0078000 0 113 211 3334 211 -211
57+
4332:addChannel = 1 0.0100000 0 211 113 3334 211 -211
58+
4332:addChannel = 1 0.0018000 0 211 221 3334 211 -211 111 22 22
59+
4332:addChannel = 1 0.0294000 0 211 111 3334 22 22
60+
4332:addChannel = 1 0.0015000 0 211 111 111 3334 22 22 22 22
61+
4332:addChannel = 1 0.0377000 0 -11 12 3334
62+
4332:addChannel = 1 0.0001000 0 211 -211 211 111 3334 22 22
63+
4332:addChannel = 1 0.0024000 0 213 3334 211 111 22 11 -11
64+
4332:addChannel = 1 0.0026000 0 221 211 3334 211 -211 111 22 22
65+
4332:addChannel = 1 0.0008000 0 321 -311 3334 310 111 111 22 22 22 22
66+
4332:addChannel = 1 0.0160000 0 211 223 3334 211 -211 111 22 22
67+
4332:addChannel = 1 0.0017000 0 223 211 3334 111 22 22 22
68+
4332:addChannel = 1 0.0037000 0 211 221 3334 22 22
69+
4332:addChannel = 1 0.0114000 0 213 111 3334 211 111 22 22 22 22
70+
4332:addChannel = 1 0.0093000 0 111 213 3334 22 22 211 111 22 22
71+
4332:addChannel = 1 0.0033000 0 221 211 3334 22 22
72+
4332:addChannel = 1 0.0040000 0 221 211 3334 111 111 111 22 22 22 22 22 22
73+
4332:addChannel = 1 0.0002000 0 213 111 111 3334 211 111 22 22 22 22 22 22
74+
4332:addChannel = 1 0.0004000 0 223 211 3334 211 -211 111 22 11 -11
75+
4332:addChannel = 1 0.0009000 0 -311 321 3334 310 111 111 22 22 22 22
76+
4332:addChannel = 1 0.0001000 0 221 211 3334 111 111 111 22 22 22 11 -11 22 22
77+
4332:addChannel = 1 0.0022000 0 -311 321 3334 130
78+
4332:addChannel = 1 0.0017000 0 211 211 -211 3334
79+
4332:addChannel = 1 0.0034000 0 211 221 3334 111 111 111 22 22 22 22 22 22
80+
4332:addChannel = 1 0.0022000 0 -311 321 3334 310 211 -211
81+
4332:addChannel = 1 0.0005000 0 213 211 -211 3334 211 111 22 22
82+
4332:addChannel = 1 0.0010000 0 111 211 111 3334 22 22 22 22
83+
4332:addChannel = 1 0.0009000 0 321 3334 -311 310 111 111 22 22 22 22
84+
4332:addChannel = 1 0.0017000 0 321 -311 3334 310 211 -211
85+
4332:addChannel = 1 0.0007000 0 111 111 211 3334 22 22 22 22
86+
4332:addChannel = 1 0.0022000 0 321 3334 -311 310 211 -211
87+
4332:addChannel = 1 0.0030000 0 321 3334 -311 130
88+
4332:addChannel = 1 0.0004000 0 211 113 111 3334 211 -211 22 22
89+
4332:addChannel = 1 0.0003000 0 221 211 3334 211 -211 22
90+
4332:addChannel = 1 0.0009000 0 211 111 3334 22 11 -11
91+
4332:addChannel = 1 0.0001000 0 213 111 3334 211 111 22 22 22 11 -11
92+
4332:addChannel = 1 0.0001000 0 321 -311 3334 310 -12 11 211
93+
4332:addChannel = 1 0.0003000 0 211 -211 213 3334 211 111 22 22
94+
4332:addChannel = 1 0.0017000 0 211 -211 211 3334
95+
4332:addChannel = 1 0.0002000 0 211 221 111 3334 22 22 22 22
96+
4332:addChannel = 1 0.0001000 0 111 213 111 3334 22 22 211 111 22 22 22 22
97+
4332:addChannel = 1 0.0008000 0 211 223 3334 111 22 22 22
98+
4332:addChannel = 1 0.0001000 0 211 -211 213 3334 211 111 22 11 -11
99+
4332:addChannel = 1 0.0002000 0 211 211 -213 3334 -211 111 22 22
100+
4332:addChannel = 1 0.0004000 0 211 221 3334 211 -211 22
101+
4332:addChannel = 1 0.0001000 0 211 221 3334 111 111 111 22 22 22 11 -11 22 22
102+
4332:addChannel = 1 0.0002000 0 213 -211 211 3334 211 111 22 22
103+
4332:addChannel = 1 0.0001000 0 211 223 3334 211 -211
104+
4332:addChannel = 1 0.0002000 0 211 221 111 3334 211 -211 111 22 22 22 22
105+
4332:addChannel = 1 0.0001000 0 211 111 221 3334 22 22 111 111 111 22 22 22 22 22 22
106+
4332:addChannel = 1 0.0001000 0 211 223 3334 211 -211 111 22 11 -11
107+
4332:addChannel = 1 0.0001000 0 221 111 211 3334 22 22 22 22
108+
4332:addChannel = 1 0.0001000 0 111 211 221 111 3334 22 22 22 22 22 22
109+
4332:addChannel = 1 0.0002000 0 221 211 3334 22 11 -11
110+
4332:addChannel = 1 0.0001000 0 213 111 3334 211 111 22 11 -11 22 22
111+
4332:addChannel = 1 0.0001000 0 111 221 211 3334 22 22 22 22
112+
4332:addChannel = 1 0.0001000 0 111 211 221 3334 22 22 111 111 111 22 22 22 22 22 22
113+
4332:addChannel = 1 0.0002000 0 213 3334 211 22
114+
4332:addChannel = 1 0.0001000 0 211 221 3334 111 111 111 22 22 22 22 22 11 -11
115+
4332:addChannel = 1 0.0001000 0 211 211 -211 221 3334 111 111 111 22 22 22 22 22 22
116+
4332:addChannel = 1 0.0001000 0 111 211 113 3334 22 22 211 -211
117+
4332:addChannel = 1 0.0002000 0 223 211 3334 211 -211
118+
4332:addChannel = 1 0.0002000 0 113 211 111 3334 211 -211 22 22
119+
4332:addChannel = 1 0.0003000 0 111 111 213 3334 22 22 22 22 211 111 22 22
120+
4332:addChannel = 1 0.0001000 0 211 -211 221 211 3334 111 111 111 22 22 22 22 22 22
121+
4332:addChannel = 1 0.0001000 0 111 213 3334 22 22 211 111 22 11 -11
122+
4332:addChannel = 1 0.0001000 0 223 211 3334 111 111 22 22 22 22 22
123+
4332:addChannel = 1 0.0001000 0 221 211 111 3334 111 111 111 22 22 22 22 22 11 -11 22 22
124+
4332:addChannel = 1 0.0001000 0 221 211 3334 111 111 111 22 11 -11 22 22 22 22
125+
4332:addChannel = 1 0.0001000 0 211 221 3334 22 11 -11
126+
4332:addChannel = 1 0.0001000 0 221 111 211 3334 111 111 111 22 22 22 22 22 22 22 22
127+
4332:addChannel = 1 0.0001000 0 211 113 3334 221 22 111 111 111 22 22 22 22 22 22
128+
4332:addChannel = 1 0.0001000 0 211 221 3334 111 111 111 22 11 -11 22 22 22 22

‎MC/config/PWGHF/pythia8/decayer/geant4_externaldecayer_omegac.in‎

Lines changed: 68 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,68 @@
1+
2+
/control/verbose 2
3+
/mcVerbose/all 1
4+
/mcVerbose/geometryManager 1
5+
/mcVerbose/opGeometryManager 1
6+
/mcTracking/loopVerbose 1
7+
/mcVerbose/composedPhysicsList 2
8+
/mcVerbose/runAction 2 # For looping thresholds control
9+
#/tracking/verbose 1
10+
#//control/cout/ignoreThreadsExcept 0
11+
12+
/mcPhysics/rangeCuts 0.001 mm
13+
/mcRegions/setRangePrecision 5
14+
/mcTracking/skipNeutrino true
15+
/mcDet/setIsMaxStepInLowDensityMaterials true
16+
/mcDet/setMaxStepInLowDensityMaterials 10 m
17+
/mcMagField/setConstDistance 1 mm
18+
/mcDet/setIsZeroMagField true
19+
/mcControl/useRootRandom true # couple G4 random seed to gRandom
20+
21+
# optical
22+
23+
/process/optical/verbose 0
24+
/process/optical/processActivation Scintillation 0
25+
/process/optical/processActivation OpWLS 0
26+
/process/optical/processActivation OpMieHG 0
27+
/process/optical/cerenkov/setTrackSecondariesFirst false
28+
/mcMagField/stepperType NystromRK4
29+
30+
# PAI for TRD
31+
# Geant4 VMC >= v3.2
32+
/mcPhysics/emModel/setEmModel PAI
33+
/mcPhysics/emModel/setRegions TRD_Gas-mix
34+
/mcPhysics/emModel/setParticles all
35+
/mcPrimaryGenerator/skipUnknownParticles true # don't crash when seeing unknown ion etc. (issue warning)
36+
37+
#
38+
# Precise Msc for EMCAL
39+
#
40+
# Geant4 VMC >= v3.2
41+
/mcPhysics/emModel/setEmModel SpecialUrbanMsc
42+
/mcPhysics/emModel/setRegions EMC_Lead$ EMC_Scintillator$
43+
/mcPhysics/emModel/setParticles e- e+
44+
45+
# combined transportation + Msc mode is currently broken for ALICE (Geant 10.2.0)
46+
/process/em/transportationWithMsc Disabled
47+
48+
#
49+
# Adding extra lines for fixing tracking bias
50+
#
51+
/mcMagField/setDeltaIntersection 1.0e-05 mm
52+
/mcMagField/setMinimumEpsilonStep 0.5e-05
53+
/mcMagField/setMaximumEpsilonStep 1.0e-05
54+
/mcMagField/printParameters
55+
56+
# Change default parameters for killing looping particles
57+
#
58+
/mcPhysics/useHighLooperThresholds
59+
/mcRun/setLooperThresholdImportantEnergy 100. MeV
60+
61+
# Define media with the INCLXX physics list; here basically in all ITS media
62+
#/mcVerbose/biasingConfigurationManager 3
63+
/mcPhysics/biasing/setModel inclxx
64+
/mcPhysics/biasing/setRegions ITS_AIR$ ITS_WATER$ ITS_COPPER$ ITS_KAPTON(POLYCH2)$ ITS_GLUE_IBFPC$ ITS_CERAMIC$ ITS_K13D2U2k$ ITS_K13D2U120$ ITS_F6151B05M$ ITS_M60J3K$ ITS_M55J6K$ ITS_FGS003$ ITS_CarbonFleece$ ITS_PEEKCF30$ ITS_GLUE$ ITS_ALUMINUM$ ITS_INOX304$ ALPIDE_METALSTACK$ ALPIDE_SI$
65+
/mcPhysics/biasing/setParticles proton neutron pi+ pi-
66+
67+
# external decayer
68+
/mcPhysics/setExtDecayerSelection omega_c0 anti_omega_c0

0 commit comments

Comments
 (0)