Skip to content

Commit b705e8a

Browse files
committed
Add configs to force decays of Xic for MC productions for triggered data
1 parent f470290 commit b705e8a

2 files changed

Lines changed: 99 additions & 0 deletions

File tree

Lines changed: 31 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,31 @@
1+
# Decay configuration taken from MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_with_decays_Mode2_pp_ref.cfg
2+
Init:showChangedParticleData = on
3+
4132:tau0 = 0.0455 # Xic0
4+
5+
# switch off all decay channels
6+
4232:onMode = off
7+
4132:onMode = off
8+
9+
## Xic decays
10+
### Ξc+ -> Ξ- π+ π+ (70%)
11+
4232:addChannel = 1 0.35000 0 3312 211 211 ### Ξc+ -> Ξ- π+ π+ 2.86%
12+
### Ξc+ -> Ξ*0 π+ (30%)
13+
4232:addChannel = 1 0.12500 0 3324 211
14+
15+
### Ξc+ -> Ξ*0 π+, Ξ*0 -> Ξ- π+
16+
4232:onIfMatch = 3324 211
17+
### Ξc+ -> Ξ- π+ π+
18+
4232:onIfMatch = 3312 211 211
19+
20+
## Xic0 decays
21+
### add Xic0 decays absent in PYTHIA8 decay table
22+
4132:oneChannel = 1 0.0143 0 3312 211 ### Xi_c()0 --> Xi- pi+ 0.01432524810113947
23+
24+
### Xic0 -> Xi- pi+
25+
4132:onIfMatch = 3312 211
26+
# Matching Exclusive Decay Channels:
27+
# -------------------------------------------------------
28+
# Xi_c()0 --> Xi- pi+ 0.01432524810113947
29+
# -> PDG Codes: 3312, 211
30+
# Xi_c()0 --> Xi- pi+ pi+ pi- 0.04775082700379833
31+
# -> PDG Codes: 3312, 211, 211, -211
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 xi_c0 anti_xi_c0 xi_c+ xi_c-

0 commit comments

Comments
 (0)