100 likes | 175 Vues
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
E N D
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
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 64-51 64-52 64-49 64-48 72-50
New global fit to everything in SPICE py, toff, fSAM, asca, aabs, llhtot, … 2.50 34 0.3 1.0 1.0 114542.271848 113957.519150 5.847495 0.000032 1833665.290289 4133262.515452 2.50 34 0.3 1.0 1.0 114542.271848 2.50 32 0.3 1.0 1.0 113146.717127 112561.964429 5.847495 0.000032 1832326.020772 4133265.231722 2.50 36 0.3 1.0 1.0 116105.789244 115521.036547 5.847495 0.000032 1835341.083155 4133190.812949 2.50 32 0.25 1.0 1.0 113048.585029 112463.832332 5.847495 0.000032 1838294.821793 4134128.901970 2.50 32 0.35 1.0 1.0 113752.347034 113167.594336 5.847495 0.000032 1833671.750123 4133400.717904 2.4 32 0.25 1.0 1.0 112989.439125 112404.686428 5.847495 0.000032 1809424.967754 4131600.193859 2.6 32 0.25 1.0 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 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 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
Timing of oversized DOM MC 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: xR=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.zR*e.zR-1)*c)-D)*e.zR-h out-7: del=e.R-OMR Flashing 63-50 63-49 63-48 63-51 63-52 64-50 64-48 64-52
BAD multiprocessors (MPs) clist cudatest 0 1 2 3 4 5 cuda001 0 1 2 3 4 5 cuda002 0 1 2 3 4 5 cuda003 0 1 2 3 4 5 #badmps cuda001 3 22 cuda002 2 20 cuda002 4 10 Configured: xR=5 eff=0.95 sf=0.2 g=0.943 Loaded 12 angsens coefficients Loaded 6x170 dust layer points Loaded 16028 random multipliers Loaded 42 wavelenth points Loaded 171 ice layers Loaded 3540 DOMs (19x19) Processing f2k muons from stdin on device 2 Total GPU memory usage: 83053520 photons: 13762560 hits: 991 Error: TOT was a nan or an inf 1 times! Bad MP #20 photons: 13762560 hits: 393 photons: 13762560 hits: 570 photons: 13762560 hits: 501 photons: 13762560 hits: 832 photons: 13762560 hits: 717 CUDA Error: unspecified launch failure [dima@cuda002 gpu]$ cat mmc.1.f2k | BADMP=20 ./ppc 2 > /dev/null Configured: xR=5 eff=0.95 sf=0.2 g=0.943 Loaded 12 angsens coefficients Loaded 6x170 dust layer points Loaded 16028 random multipliers Loaded 42 wavelenth points Loaded 171 ice layers Loaded 3540 DOMs (19x19) Processing f2k muons from stdin on device 2 Not using MP #20 Total GPU memory usage: 83053520 photons: 13762560 hits: 871 … photons: 1813560 hits: 114 Device time: 31970.7 (in-kernel: 31725.6...31954.8) [ms] Total GPU memory usage: 83053520 photons: 13762560 hits: 938 Error: TOT was a nan or an inf 9 times! Bad MP #20 #20 #20 #20 photons: 13762560 hits: 442 photons: 13762560 hits: 627 CUDA Error: unspecified launch failure Disable 3 bad GPUs out of 24: 12.5% Disable 3 bad MPs out of 720: 0.4%! Failure rates: