Updates to ppc and spice PPC LONG simulate

Updates to ppc and spice PPC: • LONG: simulate longitudinal cascade development • ANGW: smear cherenkov cone due to shower development • Made the code scalable with the number of GPU multiprocessors • The flasher simulation now uses the wavelength profile read from file wv. dat • Randomized the simulation based on system time (with us resolution) • Modified code to run CPU and GPU parts concurrently • Added option to disable a multiprocessor • Added the implementation of the simple approximate Mie scattering function • Added a configuration file "cfg. txt" • The time needed to reach the nominal (non-oversized) DOM surface is added • When DOM is oversized: re-use the photon after it hits a DOM and ensure the causality in the flasher simulation while keeping the bias low 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. Dima Chirkin, UW Madison

Photon angular profile from thesis of Christopher Wiebusch

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

New approximation to Mie

Dependence on g=<cos(q)> and f. SAM g=<cos(q)> f. SAM 0. 8 0 0. 95 0 0. 9 0. 3 0. 5 1. 0 flashing 63 -50 64 -51 64 -52 64 -49 64 -48 72 -50

Above/below the flasher

New global fit to everything in SPICE py, toff, f. SAM, asca, aabs, llhtot, … 2. 50 34 0. 3 1. 0 114542. 271848 113957. 519150 5. 847495 0. 000032 1833665. 290289 4133262. 515452 2. 50 34 0. 3 1. 0 114542. 271848 2. 50 32 0. 3 1. 0 113146. 717127 112561. 964429 5. 847495 0. 000032 1832326. 020772 4133265. 231722 2. 50 36 0. 3 1. 0 116105. 789244 115521. 036547 5. 847495 0. 000032 1835341. 083155 4133190. 812949 2. 50 32 0. 25 1. 0 113048. 585029 112463. 832332 5. 847495 0. 000032 1838294. 821793 4134128. 901970 2. 50 32 0. 35 1. 0 113752. 347034 113167. 594336 5. 847495 0. 000032 1833671. 750123 4133400. 717904 2. 4 32 0. 25 1. 0 112989. 439125 112404. 686428 5. 847495 0. 000032 1809424. 967754 4131600. 193859 2. 6 32 0. 25 1. 0 114643. 008467 114058. 255769 5. 847495 0. 000032 1867203. 704262 4136517. 482188 2. 4 32 0. 25 0. 98 1. 0 112581. 738041 111996. 983925 5. 847509 0. 000032 1828322. 503196 4133559. 858160 2. 4 32 0. 25 1. 02 1. 0 114492. 497808 113907. 746398 5. 847483 0. 000031 1790745. 309544 4129540. 675800 2. 4 32 0. 25 0. 98 113963. 336657 113378. 582876 5. 847507 0. 000030 1883860. 156599 4138858. 359864 2. 4 32 0. 25 0. 98 1. 02 112990. 183286 112405. 429204 5. 847511 0. 000029 1781714. 176100 4129307. 539358 2. 4 32 0. 25 0. 98 1. 0 112581. 738041 2. 4 30 0. 25 0. 98 1. 0 112066. 814058 111482. 059942 5. 847509 0. 000032 1828834. 382717 4133747. 256010 2. 4 34 0. 25 0. 98 1. 0 113223. 684205 112638. 930089 5. 847509 0. 000032 1829167. 340748 4133460. 368426 2. 4 30 0. 2 0. 98 1. 0 112664. 435687 112079. 681571 5. 847509 0. 000032 1830564. 675999 4134247. 526197 2. 4 30 0. 3 0. 98 1. 0 111916. 893671 111332. 139554 5. 847509 0. 000032 1825074. 226225 4133191. 418501 2. 3 30 0. 3 0. 98 1. 0 113352. 776400 112768. 022284 5. 847509 0. 000032 1799515. 691724 4131145. 935607 2. 5 30 0. 3 0. 98 1. 0 112451. 687974 111866. 933858 5. 847509 0. 000032 1856490. 845080 4135989. 607000 2. 4 30 0. 3 0. 96 1. 0 112879. 219373 112294. 466821 5. 847493 0. 000032 1850083. 084389 4135669. 806547 2. 4 30 0. 3 1 1. 0 112236. 475165 111651. 722467 5. 847495 0. 000032 1805746. 804534 4131542. 251044 2. 4 30 0. 3 0. 98 112905. 364261 112320. 610481 5. 847507 0. 000030 1879521. 070825 4138522. 367169 2. 4 30 0. 3 0. 98 1. 02 112324. 877135 111740. 123052 5. 847511 0. 000029 1778024. 480857 4129036. 281341 2. 4 30 0. 3 0. 98 1. 0 111916. 893671

Reconstructing ice from toy MC

Timing of oversized DOM MC Flashing 63 -50 63 -49 63 -48 63 -51 63 -52 64 -50 64 -48 64 -52 16 # over-R: DOM radius "oversize" scaling factor 0. 95 # overall DOM efficiency correction 0. 2 # 0=HG; 1=SAM 0. 943 # g=<cos(theta)> out-0: x. R=1 out-1: default out-2: do not track back to detected DOM out-3: do not track after detection out-4: no ovesize delta correction! out-5: do not check causality out-6: del=(sqrtf(b*b+(1/(e. z. R*e. z. R-1)*c)-D)*e. z. R-h out-7: del=e. R-OMR

BAD multiprocessors (MPs) clist cudatest cuda 001 cuda 002 cuda 003 Configured: x. R=5 eff=0. 95 sf=0. 2 g=0. 943 [dima@cuda 002 gpu]$ cat mmc. 1. f 2 k | BADMP=20. /ppc 2 > /dev/null 0 1 2 3 4 5 Loaded 12 angsens coefficients Configured: x. R=5 eff=0. 95 sf=0. 2 g=0. 943 0 1 2 3 4 5 Loaded 6 x 170 dust layer points Loaded 12 angsens coefficients 0 1 2 3 4 5 Loaded 16028 random multipliers Loaded layer points 0 1 6 x 170 2 3 4 5 dust Loaded 42 wavelenth points Loaded 16028 random multipliers Loaded 171 ice layers Loaded 42 wavelenth points #badmps Loaded 3540 DOMs (19 x 19) Loaded cuda 001 3 22171 ice layers Processing f 2 k muons from stdin on device 2 Loaded 3540 DOMs (19 x 19) cuda 002 2 20 Total GPU memory usage: 83053520 Processing from stdin on device 2 cuda 002 4 10 f 2 k muons photons: 13762560 hits: 991 Not using MP #20 Error: TOT was a nan or an inf 1 times! Bad MP #20 Total GPU memory usage: 83053520 photons: 13762560 hits: 393 photons: hits: 871 photons: 13762560 hits: 13762560 938 photons: 13762560 hits: 570 … Error: TOT was a nan or an inf 9 photons: times! Bad MP #20 13762560 hits: #20 501#20 photons: 1813560 hits: 114 photons: 13762560 hits: 442 photons: 13762560 hits: 832 photons: 13762560 hits: 627 photons: 13762560 hits: 717 Device time: 31970. 7 (in-kernel: 31725. 6. . . 31954. 8) [ms] CUDA Error: unspecified launch failure Failure rates: Disable 3 bad GPUs out of 24: 12. 5% Disable 3 bad MPs out of 720: 0. 4%!