Updates to ppc and spice

Document Sample
Updates to ppc and spice Powered By Docstoc
					            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 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
Above/below the flasher
 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
        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)>
Flashing 63-50
                   out-0: xR=1
63-49              out-1: default
                   out-2: do not track back to detected DOM
63-48              out-3: do not track after detection
                   out-4: no ovesize delta correction!
63-51              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
63-52

64-50

64-48

64-52
        BAD multiprocessors (MPs)
 clist                          Configured: xR=5 eff=0.95 sf=0.2 g=0.943
 cudatest     [dima@cuda002 gpu]$ cat mmc.1.f2k | BADMP=20 ./ppc 2 > /dev/null
                   0 1 2 3 4 5 Loaded 12 angsens coefficients
 cuda001      Configured: xR=5 eff=0.95 sf=0.2 g=0.943
                   0 1 2 3 4 5 Loaded 6x170 dust layer points
 cuda002      Loaded 123angsens coefficients
                   0 1 2 4 5 Loaded 16028 random multipliers
 cuda003           0 1 2 3 4 dust layer points
              Loaded 6x1705
                                Loaded 42 wavelenth points
              Loaded 16028 random multipliers
                                Loaded 171 ice layers
  #badmps     Loaded 42 wavelenth points
                                Loaded 3540 DOMs (19x19)
  cuda001     Loaded 171 ice layers
                   3 22         Processing f2k muons from stdin on device 2
  cuda002     Loaded 3540 DOMs (19x19)
                   2 20         Total GPU memory usage: 83053520
  cuda002     Processing f2k muons from stdin on device 2
                   4 10         photons: 13762560 hits: 991
              Not using MP #20
                                Error: TOT was a nan or an inf 1 times! Bad MP #20
                    GPU memory usage: 83053520
              Totalusage: 83053520
Total GPU memory                photons: 13762560 hits: 393
              photons: 13762560 hits: 871
photons: 13762560 hits: 938 photons: 13762560 hits: 570
Error: TOT was… nan or an inf 9 photons: 13762560 #20 #20 #20
               a                 times! Bad MP #20 hits: 501
              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]
                                failure
CUDA Error: unspecified launch CUDA Error: unspecified launch failure
                 Disable 3 bad GPUs out of 24: 12.5%
Failure rates:
                 Disable 3 bad MPs out of 720: 0.4%!

				
DOCUMENT INFO
Shared By:
Categories:
Tags:
Stats:
views:0
posted:5/21/2013
language:Unknown
pages:10