Skip to content

Add Omega_c decay table and G4 config to use external decayer #1879

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
fgrosa merged 2 commits into from
Jan 20, 2025
Merged

Add Omega_c decay table and G4 config to use external decayer #1879

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
79b8494
Add Omega_c decay table and G4 config to use external decayer
mpuccio Jan 20, 2025
3be1e29
Constrain the decay table with PDG inputs
mpuccio Jan 20, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
48 changes: 48 additions & 0 deletions MC/config/PWGHF/pythia8/decayer/force_hadronic_Omegac_forceOmegaChannels.cfg
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,48 @@
Init:showChangedParticleData = on

4332:mayDecay on # OmegaC
4332:tau0 = 0.0803
4332:addChannel = 1 0.1062323 0 211 3334
4332:addChannel = 1 0.0531161 0 111 211 3334
4332:addChannel = 1 0.1381020 0 213 3334
4332:addChannel = 1 0.0318697 0 -211 211 211 3334
4332:addChannel = 1 0.2103399 0 -13 14 3334
4332:addChannel = 1 0.2103399 0 -11 12 3334
4332:addChannel = 1 0.1439810 0 22 22 111 211 213 3334
4332:addChannel = 1 0.0328422 0 22 22 111 211 3334
4332:addChannel = 1 0.0166708 0 -211 22 22 111 211 211 223 3334
4332:addChannel = 1 0.0048077 0 -311 130 321 3334
4332:addChannel = 1 0.0111139 0 -211 113 211 211 3334
4332:addChannel = 1 0.0027473 0 -211 22 22 111 211 211 221 3334
4332:addChannel = 1 0.0019980 0 22 22 22 22 111 111 211 3334
4332:addChannel = 1 0.0000624 0 -211 22 22 111 211 211 3334
4332:addChannel = 1 0.0014985 0 -11 11 22 111 211 213 3334
4332:addChannel = 1 0.0016234 0 -311 22 22 22 22 111 111 310 321 3334
4332:addChannel = 1 0.0015609 0 22 22 22 111 211 223 3334
4332:addChannel = 1 0.0043706 0 22 22 211 221 3334
4332:addChannel = 1 0.0129246 0 22 22 22 22 111 111 211 213 3334
4332:addChannel = 1 0.0046204 0 22 22 22 22 22 22 111 111 111 211 221 3334
4332:addChannel = 1 0.0003746 0 22 22 22 22 22 22 111 111 111 211 213 3334
4332:addChannel = 1 0.0003122 0 -211 -11 11 22 111 211 211 223 3334
4332:addChannel = 1 0.0003122 0 -11 11 22 22 22 22 22 111 111 111 211 221 3334
4332:addChannel = 1 0.0008741 0 -211 22 22 111 211 211 213 3334
4332:addChannel = 1 0.0004995 0 -211 22 22 111 113 211 211 3334
4332:addChannel = 1 0.0038087 0 -311 -211 211 310 321 3334
4332:addChannel = 1 0.0004371 0 -211 22 211 211 221 3334
4332:addChannel = 1 0.0005619 0 -11 11 22 111 211 3334
4332:addChannel = 1 0.0001249 0 -211 22 22 22 111 211 211 221 3334
4332:addChannel = 1 0.0001873 0 -11 11 22 22 22 111 111 211 213 3334
4332:addChannel = 1 0.0000624 0 -311 -12 11 211 310 321 3334
4332:addChannel = 1 0.0002498 0 22 22 22 22 111 211 221 3334
4332:addChannel = 1 0.0000624 0 -211 -11 11 22 111 211 211 213 3334
4332:addChannel = 1 0.0001249 0 -213 -211 22 22 111 211 211 3334
4332:addChannel = 1 0.0001873 0 -211 211 211 223 3334
4332:addChannel = 1 0.0001249 0 -211 22 22 22 22 111 111 211 211 221 3334
4332:addChannel = 1 0.0001873 0 22 22 22 22 22 22 22 22 111 111 111 111 211 221 3334
4332:addChannel = 1 0.0000624 0 22 22 22 22 22 22 111 111 211 221 3334
4332:addChannel = 1 0.0001873 0 -11 11 22 211 221 3334
4332:addChannel = 1 0.0001249 0 22 211 213 3334
4332:addChannel = 1 0.0001249 0 -211 22 22 22 22 22 22 111 111 111 211 211 221 3334
4332:addChannel = 1 0.0000624 0 22 22 22 22 22 111 111 211 223 3334
4332:addChannel = 1 0.0000624 0 -11 11 22 22 22 22 22 22 22 111 111 111 111 211 221 3334
4332:addChannel = 1 0.0000624 0 22 22 22 22 22 22 22 111 111 111 113 211 221 3334
68 changes: 68 additions & 0 deletions MC/config/PWGHF/pythia8/decayer/geant4_externaldecayer_omegac.in
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,68 @@

/control/verbose 2
/mcVerbose/all 1
/mcVerbose/geometryManager 1
/mcVerbose/opGeometryManager 1
/mcTracking/loopVerbose 1
/mcVerbose/composedPhysicsList 2
/mcVerbose/runAction 2 # For looping thresholds control
#/tracking/verbose 1
#//control/cout/ignoreThreadsExcept 0

/mcPhysics/rangeCuts 0.001 mm
/mcRegions/setRangePrecision 5
/mcTracking/skipNeutrino true
/mcDet/setIsMaxStepInLowDensityMaterials true
/mcDet/setMaxStepInLowDensityMaterials 10 m
/mcMagField/setConstDistance 1 mm
/mcDet/setIsZeroMagField true
/mcControl/useRootRandom true # couple G4 random seed to gRandom

# optical

/process/optical/verbose 0
/process/optical/processActivation Scintillation 0
/process/optical/processActivation OpWLS 0
/process/optical/processActivation OpMieHG 0
/process/optical/cerenkov/setTrackSecondariesFirst false
/mcMagField/stepperType NystromRK4

# PAI for TRD
# Geant4 VMC >= v3.2
/mcPhysics/emModel/setEmModel PAI
/mcPhysics/emModel/setRegions TRD_Gas-mix
/mcPhysics/emModel/setParticles all
/mcPrimaryGenerator/skipUnknownParticles true # don't crash when seeing unknown ion etc. (issue warning)

#
# Precise Msc for EMCAL
#
# Geant4 VMC >= v3.2
/mcPhysics/emModel/setEmModel SpecialUrbanMsc
/mcPhysics/emModel/setRegions EMC_Lead$ EMC_Scintillator$
/mcPhysics/emModel/setParticles e- e+

# combined transportation + Msc mode is currently broken for ALICE (Geant 10.2.0)
/process/em/transportationWithMsc Disabled

#
# Adding extra lines for fixing tracking bias
#
/mcMagField/setDeltaIntersection 1.0e-05 mm
/mcMagField/setMinimumEpsilonStep 0.5e-05
/mcMagField/setMaximumEpsilonStep 1.0e-05
/mcMagField/printParameters

# Change default parameters for killing looping particles
#
/mcPhysics/useHighLooperThresholds
/mcRun/setLooperThresholdImportantEnergy 100. MeV

# Define media with the INCLXX physics list; here basically in all ITS media
#/mcVerbose/biasingConfigurationManager 3
/mcPhysics/biasing/setModel inclxx
/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$
/mcPhysics/biasing/setParticles proton neutron pi+ pi-

# external decayer
/mcPhysics/setExtDecayerSelection omega_c0 anti_omega_c0
Loading