Updates to ppc and spice
This presentation is the property of its rightful owner.
Sponsored Links
1 / 10

Updates to ppc and spice PowerPoint PPT Presentation


  • 25 Views
  • Uploaded on
  • Presentation posted in: General

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

Download Presentation

Updates to ppc and spice

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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Updates to ppc and spice

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

Photon angular profile

from thesis of Christopher Wiebusch


Simplified mie scattering

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

New approximation to Mie


Dependence on g cos q and f sam

Dependence on g=<cos(q)> and fSAM

g=<cos(q)> fSAM

0.80

0.90

0.950

0.90.3

0.90.5

0.91.0

flashing 63-50

64-50

64-51

64-52

64-49

64-48

72-50


Above below the flasher

Above/below the flasher


New global fit to everything in spice

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


Reconstructing ice from toy mc

Reconstructing ice from toy MC


Timing of oversized dom mc

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

BAD multiprocessors (MPs)

clist

cudatest0 1 2 3 4 5

cuda0010 1 2 3 4 5

cuda0020 1 2 3 4 5

cuda0030 1 2 3 4 5

#badmps

cuda0013 22

cuda0022 20

cuda0024 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

[[email protected] 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:


  • Login