1 / 16

SPICE Mie [mi:]

SPICE Mie [mi:]. Dmitry Chirkin, UW Madison. Updates to ppc and spice. PPC: Randomized the simulation based on system time (with us resolution) Added the implementation of the simple approximate Mie scattering function

albany
Download Presentation

SPICE Mie [mi:]

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. SPICE Mie [mi:] Dmitry Chirkin, UW Madison

  2. Updates to ppc and spice • PPC: • Randomized the simulation based on system time (with us resolution) • Added the implementation of the simple approximate Mie scattering function • New oversized DOM treatment (designed for minimum bias compared to oversize=1): • oversize only in direction perpendicular to the photon •  time needed to reach the nominal (non-oversized) DOM surface is added • re-use the photon after it hits a DOM and ensure the causality in the flasher simulation • Spice: • Fixed code determining the closest DOMs to the current layer (when using tilted ice) • Perform simultaneous global fit for py, time offset, scattering vs. absorption correlation coeff. • Optimize use of high-event flasher simulation: use 250-event simulation in the dust peak, 10 elsewhere. Eventually use 250-event simulation for the entire depth range. nominal DOM oversized DOM oversized ~ 5 times photon

  3. Timing of oversized DOM MC xR=1 default  64-48 Flashing 63-50  63-48 xR=1 default do not track back to detected DOM do not track after detection no ovesize delta correction! do not check causality del=(sqrtf(b*b+(1/(e.zR*e.zR-1)*c)-D)*e.zR-h del=e.R-OMR  64-52

  4. Simplified Mie Scattering Also known as the Liu scattering function Introduced by Jon Miller Single radius particles, described better as smaller angles by SAM

  5. New approximation to Mie fSAM

  6. Dependence on g=<cos(q)> and fSAM g=<cos(q)> fSAM 0.8 0 0.9 0 0.95 0 0.9 0.3 0.9 0.5 0.9 1.0 flashing 63-50  64-50

  7. New global fit to everything in SPICE 1. For some starting values, find best values of lsca ~ labs. 2. Find best values of py, toff, fSAM, asca, aabs, llhtot, … py photon yield factor toff global time offset (rising edge of the flasher pulse) fSAM fraction of SAM contribution to the scattering function asca scaling of scattering coefficient aabs scaling of absorption coefficient 3. Repeat until converged (~3 iterations) 4. Refine the fit with lsca and labs independent from each other Charge only Full likelihood with timing

  8. Verification with toy simulation Input table Simulated 60 x 250 events Reconstructed table with 10 event/flasher 250 event/flasher In the dust peak

  9. Correlation with dust logger With 10 events/flasher, 250 in dust peak With 250 events/flasher everywhere

  10. Plots for individual flashers SPICE Mie  AHA 

  11. Plots for CORSIKA/data  SPICE Mie   AHA 

  12. Plots from Anne (CORSIKA IC40)

  13. Plot from Jacob Feintzeig

  14. SPICE Mie: ice coefficients

  15. “unresolved” systematics: resolved Minimum is in the same place with both likelihoods!

  16. Conclusions • SPICE Mie is great •  fits timing perfectly! • SPICE paper is available (v. 0.01) •  please comment

More Related