Down going Muon Analysis Charm Production in pN interactions

1
Down going Muon Analysis/ Charm Production in p-N
interactions

• Newt Ganugapati and
• Teresa Montaruli

2
Methodology

• Cosmic Ray muons are simulated using Corsika
  6.03 and QGSJET01 hadronic
• interaction model.
• Detector response with trigger conditions
  reconstruction (convoluted Pandel 64 iterations
• and early hits).

• Physics measurement of the true muon
  Intensity from experimental data will have to
  account for triggering reconstruction
  efficiencies.

• Method of diagonal unfolding from MC studies
• Finite angular resolution interbin
  correlations in zenith angle are not negligible.
  Adjust the bin size in zenith angle in such a way
  as to minimize them (leakage from one to other).
  Larger bin sizes when angular resolution is bad!
• Basic quality cuts to improve reconstruction
• Track Length, Chi Square, Difference of Odd and
  Even hits from a time ordered hit sample.

3
Muon flux vs zenith angle
2001

• www.amanda.wisc.edu/newt/reportfinal.pdf

4
SIMULATION DISCREPANCIES AND INTERACTION MODEL
Energy(GeV)
AMANDA data as well as L3 seem to prefer Lower
predictions as Sibyll NEW QGSJET02 should be
tested (J Kelley et al)
L3C data Source Ralph Engel
5
Angular resolution close to the vertical
6
Angular Resolution close to horizonzenith offset
reverses
Space Angle
Zenith
RecoSpace-TrueSpace(degrees)
RecoSpace-TrueSpace(degrees)
Before quality
Note the skewed nature of these distributions!
After quality
7
Zenith and Space Angle Resolution
After quality cuts Before QC
After quality cuts Before QC
The point is marked at the mean while the error
bar quantifies RMS value
8
Time Residuals and Convoluted Pandel
Patched Pandel uses a pdf for the likelihood
fit, which consists of 3 distinct functions (a
gaussian, a polynomial and the pandel function)
which are 'patched' together at 0ns residual
time and around 50ns time. Convoluted Pandel
convolution of Pandel and Gauss (analytical
solution)
ConvolutedPandel (SGR analysis) http//www.icecube
.wisc.edu/zornoza/sgr/Monte_Carlo
http//www.ifh.de/ackerman/2004/plots/Restime_03_
04.gif
9
Air Shower
Dmesons also produced (D0,D,Ds and other Kinds)
that decay to Prompt muons
During the production of these D mesons are
there any p / K also that are produced
(that decay into muons) ?If so statistically
what fraction of Momentum/energy do these p /
K take?
10
PYTHIA SIMULATION

• Make pp with PYTHIA then consider the momentum
  taken by the charm production with appropriate
  energy processes
• Subtract it to the Pbeam and use the Pbeam to
  make an event P-N with FRITIOF. This should
  account for the nucleons and production of
  accompanying p , K
• Estimate energy fractions and multiplicities in
  first interaction
• Better use FLUKA

11
1st interaction of p on NNumber of Particle and
Energy (Mean)
Eproton(GeV)
Eproton (GeV)
The RMS value is shown as an error bar!
Only events where there is charm (multiplicity
always 2)
12
COMPARISION OF FLUKA/PYTHIAFor 100TeV INCIDENT
PROTONmultiplicity distribution for Dcharged
There are normalization concerns so use FLUKA!
13
FLUKA dn/dx vs x(thanks G Battistoni)
PN at 100 TeV
The bump at high X is due to charmed baryons
since they are produced with quarks from the
projectile (diffractive) In these HE conditions
similar dn/dx for K e p. A leading baryon or
quasi-leading produced in the final state is made
by a projectile di-quark The di-quark takes on
average a large fraction of the projectile
momentum, hence it appears at large x . Eg LC
udc u-d come on average more from the proton
projectile and c must be taken from the sea. XsiC
usc or dsc, depending on charge. u-d come from
projectile
14
Charmed mesons