Skip to content

Commit ba73282

Browse files
mpuccioalcaliva
mpuccio
authored and
alcaliva
committed
Add Omega_c decay table and G4 config to use external decayer (#1879)
* Add Omega_c decay table and G4 config to use external decayer Co-authored with Fabrizio Grosa <fabrizio.grosa@cern.ch> * Constrain the decay table with PDG inputs (cherry picked from commit 6c0023e)
1 parent 8a860ab commit ba73282

File tree

2 files changed

+116
-0
lines changed

2 files changed

+116
-0
lines changed

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

Lines changed: 48 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,48 @@
1+
Init:showChangedParticleData = on
2+
3+
4332:mayDecay on # OmegaC
4+
4332:tau0 = 0.0803
5+
4332:addChannel = 1 0.1062323 0 211 3334
6+
4332:addChannel = 1 0.0531161 0 111 211 3334
7+
4332:addChannel = 1 0.1381020 0 213 3334
8+
4332:addChannel = 1 0.0318697 0 -211 211 211 3334
9+
4332:addChannel = 1 0.2103399 0 -13 14 3334
10+
4332:addChannel = 1 0.2103399 0 -11 12 3334
11+
4332:addChannel = 1 0.1439810 0 22 22 111 211 213 3334
12+
4332:addChannel = 1 0.0328422 0 22 22 111 211 3334
13+
4332:addChannel = 1 0.0166708 0 -211 22 22 111 211 211 223 3334
14+
4332:addChannel = 1 0.0048077 0 -311 130 321 3334
15+
4332:addChannel = 1 0.0111139 0 -211 113 211 211 3334
16+
4332:addChannel = 1 0.0027473 0 -211 22 22 111 211 211 221 3334
17+
4332:addChannel = 1 0.0019980 0 22 22 22 22 111 111 211 3334
18+
4332:addChannel = 1 0.0000624 0 -211 22 22 111 211 211 3334
19+
4332:addChannel = 1 0.0014985 0 -11 11 22 111 211 213 3334
20+
4332:addChannel = 1 0.0016234 0 -311 22 22 22 22 111 111 310 321 3334
21+
4332:addChannel = 1 0.0015609 0 22 22 22 111 211 223 3334
22+
4332:addChannel = 1 0.0043706 0 22 22 211 221 3334
23+
4332:addChannel = 1 0.0129246 0 22 22 22 22 111 111 211 213 3334
24+
4332:addChannel = 1 0.0046204 0 22 22 22 22 22 22 111 111 111 211 221 3334
25+
4332:addChannel = 1 0.0003746 0 22 22 22 22 22 22 111 111 111 211 213 3334
26+
4332:addChannel = 1 0.0003122 0 -211 -11 11 22 111 211 211 223 3334
27+
4332:addChannel = 1 0.0003122 0 -11 11 22 22 22 22 22 111 111 111 211 221 3334
28+
4332:addChannel = 1 0.0008741 0 -211 22 22 111 211 211 213 3334
29+
4332:addChannel = 1 0.0004995 0 -211 22 22 111 113 211 211 3334
30+
4332:addChannel = 1 0.0038087 0 -311 -211 211 310 321 3334
31+
4332:addChannel = 1 0.0004371 0 -211 22 211 211 221 3334
32+
4332:addChannel = 1 0.0005619 0 -11 11 22 111 211 3334
33+
4332:addChannel = 1 0.0001249 0 -211 22 22 22 111 211 211 221 3334
34+
4332:addChannel = 1 0.0001873 0 -11 11 22 22 22 111 111 211 213 3334
35+
4332:addChannel = 1 0.0000624 0 -311 -12 11 211 310 321 3334
36+
4332:addChannel = 1 0.0002498 0 22 22 22 22 111 211 221 3334
37+
4332:addChannel = 1 0.0000624 0 -211 -11 11 22 111 211 211 213 3334
38+
4332:addChannel = 1 0.0001249 0 -213 -211 22 22 111 211 211 3334
39+
4332:addChannel = 1 0.0001873 0 -211 211 211 223 3334
40+
4332:addChannel = 1 0.0001249 0 -211 22 22 22 22 111 111 211 211 221 3334
41+
4332:addChannel = 1 0.0001873 0 22 22 22 22 22 22 22 22 111 111 111 111 211 221 3334
42+
4332:addChannel = 1 0.0000624 0 22 22 22 22 22 22 111 111 211 221 3334
43+
4332:addChannel = 1 0.0001873 0 -11 11 22 211 221 3334
44+
4332:addChannel = 1 0.0001249 0 22 211 213 3334
45+
4332:addChannel = 1 0.0001249 0 -211 22 22 22 22 22 22 111 111 111 211 211 221 3334
46+
4332:addChannel = 1 0.0000624 0 22 22 22 22 22 111 111 211 223 3334
47+
4332:addChannel = 1 0.0000624 0 -11 11 22 22 22 22 22 22 22 111 111 111 111 211 221 3334
48+
4332:addChannel = 1 0.0000624 0 22 22 22 22 22 22 22 111 111 111 113 211 221 3334

‎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)