150 likes | 222 Views
Simulating the Solar Shadow. Allen I. Mincer NYU LANL 12/18/06. Stoermer Rule, as Modified by Henrici, adapted by AIM. Some More Calculation Details. Using Liouville’s theorem, isotropic C.R.s + B fields give isotropic C.R.s at Earth, unless absorbed, so
E N D
Simulating the Solar Shadow Allen I. Mincer NYU LANL 12/18/06
Some More Calculation Details • Using Liouville’s theorem, isotropic C.R.s + B fields give isotropic C.R.s at Earth, unless absorbed, so • Generate backwards going C.R.s from Milagro to Sun or Moon. Shadow if Sun or Moon is hit.
New Runs • Built in parametrization of MC of Milagro effective area and angular resolution from /data/links/G4MCData/Milagro/sets/20/proton/scattering_0.99_AbsL30m_ScatL50m/linkdirtrig/G4_proton_006769.trig.dat.gz • 1 year, take Sun/Moon position every 100 seconds. • For each position, select energy from E-2.7 spectrum starting at 0.1 TeV. • For each position generate 1.6K events centered on the vertices of a 40 x 40 grid ± 10 degrees around straight line to Sun/Moon : 0.1 degree steps in Θ, ΦcosΘ. • -> 5 x 108particles tracked per run. • Weight events by effective area. • Run varying fields: B Earth only, B Earth + Sun dipole parallel to Earth’s dipole, B Earth + Sun dipole perp. to Earth’s dipole. • Run with and without angular resolution
MC/Data Compare θ distribution Thanks Brian!
Shadowed Particle Energy Sun Shadowed Particle energy Moon Shadowed Particle energy
Sum over 1.2 degree radius bins Sun Tot B Moon Tot B Sun Rot B
Conclusions: • Solar large scale shadow shape is overwhelmingly determined by Milagro resolution. • Small scale (~1 degree) shape may have more sensitivity to B field • Shadow depth is determined to first order by whether particles strike sun. ie., strength of field. Next effect in MC is scattering out of 10 degrees. • How sensitive to real differences? • Next Step – Look at data, then depending on how rich data statistics are: • Build in realistic fields to test. • Build in better Milagro resolution