UHE Simulation Efforts @ Chiba U. Motivation Motivation Java class structure Java class structure Provided Methods Provided Methods Future Plan Future.

UHE Simulation Efforts @ Chiba U.UHE Simulation Efforts @ Chiba U.

MotivationMotivationJava class structure Java class structure Provided MethodsProvided MethodsFuture PlanFuture Plan

Shigeru YoshidaShigeru YoshidaDept. of PhysicsDept. of Physics

CHIBA Univ.CHIBA Univ.

Laguna Beach 2003

Please integrate itTo YOUR library!!

Motivation.. Handle UHE(EeV ] eventsMotivation.. Handle UHE(EeV ] events

Include some minor interactionsInclude some minor interactions Update the CC/NC/Photo-Update the CC/NC/Photo-

Nuclear cross sectionsNuclear cross sections Capability to consider curvature Capability to consider curvature

of the polar ice layerof the polar ice layer

Laguna Beach 2003

e

e

e/

e/

Weak

Weak

Weak

Inco

min

gProducts

Weak

Weak

Weak

Cascades

Cascades

Decay DecayWeak

Pair/decayBremss Pair Pair PhotoNucl.

PhotoNucl.DecayPair Pair

PairBremssDecay

DecayWeak

Laguna Beach 2003

The Whole Structure

Java_lib

sources

classes

icecube UHE

particles

points

interactions decay

propagation

numRecipes

Laguna Beach 2003

particles/

Defines particles with mass, lifetime, and energy

flavor

doublet

0

0

1

1 2 3

e

e 0

+

10 GeV < E < 10 GeV6 11(Initial. Changeable)

Particle.java ParticleArray.java

Laguna Beach 2003

points/ParticlePoint.java

•Define/calculate the point vector to trace particles•Define/calculate the medium(Rock/Ice) in which particles travel

density, mass number, charge, . . .

In Future. . .

Xxxx.javaTake the configuration of the Antarctica Ice

Laguna Beach 2003

interactions/

Interactions.java (abstract)•Total/differential cross sections•Integral differential cross sections•Integral inelasticity term

ddy y

Bremsstrahlung.java KnockOnElectrons.java

NeutrinoCharge.java NeurtinoNeutral.java PairCreation.java

InteractionsMatrix.javaGenerate matrix of (logE), d/dy(logE)

Laguna Beach 2003

Laguna Beach 2003

decay/Decay.java

Calculate weak decay matrixWl

MuDecayMatrix.java TauDecayMatrix.java

Laguna Beach 2003

propagation/

PropagationMatrix.javaResolves the transport equation numerically

Build the transfer matrix for propagation overInfinitesimal distance dX[g/cm2], first

e

e

e/

e/

Weak

Weak

Weak

Inco

min

gProducts

Weak

Weak

Weak

Cascades

Cascades

Decay DecayWeak

Pair/decayBremss Pair Pair PhotoNucl.

PhotoNucl.DecayPair Pair

PairBremssDecay

DecayWeak

Laguna Beach 2003

Application Example of the PropagationMatrix

Laguna Beach 2003

Laguna Beach 2003

propagation/

PropagationMatrix.javaHas a function of switching on/off any interaction channel

bit6 bit5 bit4 bit3 bit2 bit1 bit0bit7

CCNCe+e-+-+-

Knock-onBremssPhoto-Nucl ->

Laguna Beach 2003

A factor of 2 enhancementDue to the heavier lepton production

Laguna Beach 2003

numRecipes/

Integration.javaNumerical Integration Methods

• Extended Trapezoidal rule• Romberg Integration

Interpolation.javaPolynomial Interpolation

SpecialFunctions.javaGaussian, Poisson, Gamma Function etc..

RandomDouble.java

Uniform DeviatesBased on “Mersenne Twister”Period 2^19937-1 623 dimensional equi-distribution

Future plan… Integration with YOUR Future plan… Integration with YOUR code!!code!!

Detector MC on the same base infrastructure.Detector MC on the same base infrastructure. Photon Propagator. Photonics? It needs to be faster and capable of LAEGER target volume for UHE as many events Photon Propagator. Photonics? It needs to be faster and capable of LAEGER target volume for UHE as many events OUTSIDEOUTSIDE the instrumental the instrumental

volume are likely to be triggered.volume are likely to be triggered. Events are horizontal/downward-going. Needs realistic mapping of the polar ice configuration in the simulation.Events are horizontal/downward-going. Needs realistic mapping of the polar ice configuration in the simulation. Reconstruction. UHE events are associated with many cascades which would cause sizable fluctuation.Reconstruction. UHE events are associated with many cascades which would cause sizable fluctuation. May need to run Geant4 for full simulation of UHE cascades and to be implemented in this library.May need to run Geant4 for full simulation of UHE cascades and to be implemented in this library. LET ‘S WORK TOGETHER!!LET ‘S WORK TOGETHER!!

Laguna Beach 2003