GalProp

1
GalProp

• Initial functionality
• - Accurate propagation 2D(r,z,p) and 3D(x,y,z,p)
  model with isotropic diffusion (
  Dxxconst), reacceleration and convection.
• - Primary injection spectrum for electrons and
  ions with r-z source distribution
• - Possibility of stochastic primary sources (SNR
  explosion)
• - Reasonable r-z parametrization of the gas
  distribution (nHI, HII,H2)
• - Radial model of global galactic Bfield (r,z)
• - Very detailed description of charged particles
  energy losses(ion.,bremss, IC)
• - Detailed implimentation of interactions
  (pp,pN,HeN) and spallation
• of nuclei component up to Z26.
• - Iterative network of secondary particles
  abundances
• - Detailed gamma skymap descritpion ( IC based
  on realistic ISRF(l,b))
• - and more

Added to existing code - DM halo and DM source
functions from DarkSusy for pb,ep,e and gamma. -
anisotropic propagation Dxx(r,z) and nonuniform
convection V(r). - variable XCO(r) for the gas
distribution.
2
GalProp
Motivation - minuit minimization to the
existing experimental data p,pb,e,ep,
he,B/c,He3/He4,Be19/Be9,gamma etc which
require merging of simulation and analysis parts
and optimization of memory allocations. -
introduce detector elements (acceptance,
contaminations) for discovery plots. - improve
maintenance and developments(heliosphere prop.)
Restructure the code, add functionality Code
is written! Debugging-testing stage now. New
subpackages structure and readiness

Working close to not
yet
GpGamma
GpGalaxy
GpPropagation
GpAnalysis
GpCommon
5 weeks are needed to finish (with the
present speed)
(Also A.Nordt)
GpParticles
3
CMS LM9
Tight time schedule to April first results with
the full CMS simulation chain are needed (for
30fb-1)
MC
Analysis (2me or 2em)
PAX (MartinValery)
Generator level ( ntpl) PYTHIA,HERWIG,COMPHEP
etc CMKIN
FAMOS ROOT (Valery)
Fast simulation FAMOS
Full MC Geant4 (hits) OSCAR
ORCA ROOT (Valery) PAX (Anja,Martin)
Digitization(a la real) (digi) ORCA
Reconstruction (dst) ORCA
ROOT (Valery)
Have to synchronize analysis in each step! Agree
on root tree structure.
Full production chain is ready in ekpplus!
4
CMS LM9
Signal and Backgrounds cross sections (fb) signal
(2me.2em) chi2chi- 48 bkg
(2lX) susy 5000 ttbar
24000 ZW
390 ZZ 160
WW
3300 Zbbar 3000 Zttbar
1600 Wbbar
170
Part of backgrounds can be done with FAMOS (after
some verification with OSCAR)
Need NLO cross sections!
Need carefull settings of the datacards and
feedback from generator level analysis before
mass production!
PileUp of minimum bias on ORCA(FAMOS) level ?
HI 17.5 LI 3.5
5
AMS status
1. Analysis meeting CERN 10.02 (next
attempt...) - MC simulation model -
physics (CR, DM, Gamma sources)
2. TRD G4 model. M.Schmanau. Inroduce to the
AMS for fututre (now still G3 model)
3. Propagation in the heliosphere, solar
modulation... (Milano group)Mark HeroldValery.
Geophysical acceptances.
4. new MC mass production starts in february
revive Condor(or ekpplus?) (ValeryMarkCristia
n)
5. Discovery potential for AMS (Pamela. Glast)
update with GalProp (ValeryA.Nordt)
6
Antiprotons.
Problem Taking the DM halo complaint with the
EGRET data causes factor 8 too much antiprotons
as compare with BESS data. Most of
antiprotons are coming from the solar (2 kpc)
vicinity.
Tried solutions - anisotropic propagation D(r),
Vconvection(r,z) - relax a bit on angular fit to
EGRET data and modify DM halo. - reduce the
EGRET gamma excess by modifying the GalProp
model.
Possible - heliosphere acceptance, p and pb
propagation may help - cross sections? - some
compromise solution smaller EGRET excess by
modified GalProp model, propagation, optimized
DM halo, etc.