![Page 1: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/1.jpg)
Olaf HartmannStefan-Meyer-Institute for Subatomic Physics
Austrian Academy of SciencesVienna, Austria
Hadron Induced ReactionsBackground Calculations with
FLUKA Workshop, Ericeira/Portugal, October 4th-8th, 2010
![Page 2: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/2.jpg)
Overview
● Motivation – Physics Case– PANDA@FAIR
● Implementation with FLUKA– ROOT output files (FairRoot, PandaRoot)
● Comparison with other event generators● Discussion, Outlook
– Porting Geometry
– Detector Radiation Map
![Page 3: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/3.jpg)
Antiprotons at FAIR-HESR
• Antiproton ring HESR • Injection at 3.5 GeV/c• 1.5 – 15 GeV/c (0.83 – 14 GeV)• 1010 to 1011 stored particles• bunched beam• high resolution mode
• 2·1031cm-2s-1 (1010 pbar)• σ
p/p ≤ 2·10-5
• p ≤ 9 GeV/c, e− cooling• high luminosity mode
• 2·1032cm-2s-1 (1011 pbar)• σ
p/p ~ 10-4
• p ≤ 15 GeV/c, stochastic cooling
FAIR officially startedon monday this week!
![Page 4: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/4.jpg)
Internal Target (Pellets, Cluster Jet, ...)
www-panda.gsi.de
SolenoidDipole
![Page 5: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/5.jpg)
PANDA Physics Program
● Hadron Spectroscopy– Charmonium states, Open Charm
– Exotics (Glueballs, Hybrids)
– Baryon-Antibaryon
● Nucleon Structure (γγ)
● In-Medium Effects● S=-1, S=-2 Hypernuclei
– keV γ spectroscopy
pbar - p
pbar - A
![Page 6: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/6.jpg)
Implementation
● Starting point: using FLUKA as “event generator” for generic pbar-p/pbar-A events
● Interface: root files (example from fluka homepage)
● mgdraw.f → BXDRAW– Register boundary-crossing tracks which
emerge from the target– Collect them in a root tree (ascii file)– Read output into PandaRoot
![Page 7: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/7.jpg)
FLUKA input
Antiproton beam Δp/p = 10-4, FWHM = 3 mm
![Page 8: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/8.jpg)
FLUKA input
![Page 9: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/9.jpg)
Simple Geometry
Hydrogen target:
Gas volumed = 2 cmh = 20 cm
Nuclear targets:
“wire piece”d = 500 µmh = 20 cm
![Page 10: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/10.jpg)
Magnetic Field
B=(0,0,2T)
![Page 11: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/11.jpg)
USRBIN
r(z) beamparticles
r(z) charged particles
![Page 12: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/12.jpg)
MGDRAW
boundary crossingfrom the target to theoutside
energy > 0
at least oneinteraction
pid conversion toPDG codes
copy to arrays,later filled intoROOT Tree
![Page 13: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/13.jpg)
The ROOT File
![Page 14: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/14.jpg)
Encountered Problems
● In mgdraw.fJTRACK has to be >0. In rare cases, fragments with JTRACK<0 occur – and cause a crash when accessing AM(JTRACK).
● Particles curling in the magnetic fieldThe same particle crosses the region boundary many times → particle array overflow, weird output written to log file. Solution: magnetic field outside switched off
● LTRACK>1 – do not count beam particles which did not interact (majority)
![Page 15: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/15.jpg)
Comparison with other Generators
● Generic background from pbar-p reactions:Dual Parton Model (DPM)
– Elastic and inelastic scattering
– Fast calculation of large no. of events
● pbar-A reactions: UrQMD + SMMSMM – Statistical Multifragmentation Model
– SMM works on baryons from UrQMD
– Calculation takes long (10 kevents~6 h)
● Note: 1:1 comparison difficult – no “target” in DPM/UrQMD
![Page 16: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/16.jpg)
Comparison: pbar A
Used target materials: 12C, 40Ca, (63)CuBeam momenta: 3.0 GeV/c (Ξ+Ξ-), 4.1 GeV/c (J/Ψ), 6.2 GeV/c (Ψ')
normalization:no. of events
Peaks aroundBeam momenta:
Only in Flukaat 3.0/4.1 GeV/c
![Page 17: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/17.jpg)
Different normalization:number of nucleons (n+p)
Agreement for .3<p<1.0
Striking difference:3.0 GeV/c
![Page 18: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/18.jpg)
UrQMD:- no antibaryons- few electrons- π0 not decayed
PDG Code distributions
![Page 19: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/19.jpg)
Best agreement: 4.1 GeV/c pbar on 40Ca
![Page 20: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/20.jpg)
Only in Fluka:Antiprotons(i.a. scattered beam particles)
Peaks: beamEnergy.Missing for 6.2 GeV/c + CuIm
porta
nt fe
atur
e in
Fluka!
![Page 21: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/21.jpg)
pbar p
Nucleons with p≈0 in DPM Shifts for small momenta
![Page 22: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/22.jpg)
some phase spacedistributions ...
Protons
![Page 23: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/23.jpg)
![Page 24: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/24.jpg)
![Page 25: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/25.jpg)
Angular distributions
scattering peak inFluka much broaderand shifted to smallerangles
![Page 26: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/26.jpg)
Charged Pions
DPM hasmuch higheryield at smallangles
![Page 27: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/27.jpg)
Discussion and Outlook
● Fluka works as “event generator”● No target effects in DPM and UrQMD+SMM● More realistic beam-target geometries in Fluka
– Example: restgas particles
● For thin targets too many events are without any interaction (slows down the simulation)
– Biasing?
● For multiuser purposes: slim input files as interface? Code distribution? (Licenses)
● Particle transport has to be done in Geant3/4
![Page 28: Hadron Induced Reactions - FLUKA€¦ · Hadron Spectroscopy – Charmonium states, Open Charm – Exotics (Glueballs, Hybrids) – Baryon-Antibaryon Nucleon Structure (γγ) In-Medium](https://reader034.vdocument.in/reader034/viewer/2022051918/600abde821b70c67302de5be/html5/thumbnails/28.jpg)
Outlook
● Next step: make the Panda geometry available in Fluka
– In a flexible manner (frequent changes)
– Radiation map of the detector (useful tools like fluscw.f, ...)
● Fluka integration into PandaRoot? (VMC)● Revisit steering parameters for simulations with
Fluka● Later on, comparison/adaption to real data