Modified hist & beam macros to incorporate realistic LED theta-distribution from datasheet.

beam.mac

@@ -35,15 +35,16 @@
 /gps/ang/rot2 0 0 1
 /gps/ang/type beam1d
 
+
 #----- Set beam angular dispersion from histogram -----
-/gps/ang/surfnorm false
-/gps/ang/type iso
-/gps/ang/mintheta 0.0
-/gps/ang/maxtheta 90.0 deg
-
-#/control/execute hist.mac
-#/gps/ang/type user
-#/control/execute hist.mac
+#/gps/ang/surfnorm false
+#/gps/ang/type iso
+#/gps/ang/mintheta 0.0
+#/gps/ang/maxtheta 90.0 deg
+
+
+/gps/ang/type user
+/control/execute hist.mac
 
 #----- Set beam angular dispersion for fiber optic -----
 # N.A. is 0.22 so acceptance cone is ~25.4 degrees

hist.mac

-/gps/hist/type theta
-/gps/hist/point 0.000000 0.000000
-/gps/hist/point 0.008726 0.001677
-/gps/hist/point 0.026179 0.002438
-/gps/hist/point 0.043632 0.002966
-/gps/hist/point 0.061085 0.003407
-/gps/hist/point 0.078538 0.003771
-/gps/hist/point 0.095990 0.004110
-/gps/hist/point 0.113443 0.004415
-/gps/hist/point 0.130896 0.004724
-/gps/hist/point 0.148349 0.005010
-/gps/hist/point 0.165801 0.005293
-/gps/hist/point 0.183254 0.005592
-/gps/hist/point 0.200707 0.005862
-/gps/hist/point 0.218160 0.006130
-/gps/hist/point 0.235613 0.006412
-/gps/hist/point 0.253065 0.006682
-/gps/hist/point 0.270518 0.006950
-/gps/hist/point 0.287971 0.007234
-/gps/hist/point 0.305424 0.007515
-/gps/hist/point 0.322876 0.007781
-/gps/hist/point 0.340329 0.008063
-/gps/hist/point 0.357782 0.008323
-/gps/hist/point 0.375235 0.008598
-/gps/hist/point 0.392687 0.008907
-/gps/hist/point 0.410140 0.009175
-/gps/hist/point 0.427593 0.009454
-/gps/hist/point 0.445046 0.009775
-/gps/hist/point 0.462499 0.010074
-/gps/hist/point 0.479951 0.010385
-/gps/hist/point 0.497404 0.010707
-/gps/hist/point 0.514857 0.011049
-/gps/hist/point 0.532310 0.011366
-/gps/hist/point 0.549762 0.011748
-/gps/hist/point 0.567215 0.012158
-/gps/hist/point 0.584668 0.012550
-/gps/hist/point 0.602121 0.012953
-/gps/hist/point 0.619574 0.013384
-/gps/hist/point 0.637026 0.013836
-/gps/hist/point 0.654479 0.014333
-/gps/hist/point 0.671932 0.014900
-/gps/hist/point 0.689385 0.015575
-/gps/hist/point 0.706837 0.016303
-/gps/hist/point 0.724290 0.017249
-/gps/hist/point 0.741743 0.018648
-/gps/hist/point 0.759196 0.019849
-/gps/hist/point 0.776649 0.020756
-/gps/hist/point 0.794101 0.021327
-/gps/hist/point 0.811554 0.019696
-/gps/hist/point 0.829007 0.018203
-/gps/hist/point 0.846460 0.017120
-/gps/hist/point 0.863912 0.016291
-/gps/hist/point 0.881365 0.015641
-/gps/hist/point 0.898818 0.015079
-/gps/hist/point 0.916271 0.014629
-/gps/hist/point 0.933724 0.014212
-/gps/hist/point 0.951176 0.013883
-/gps/hist/point 0.968629 0.013545
-/gps/hist/point 0.986082 0.013284
-/gps/hist/point 1.003535 0.013057
-/gps/hist/point 1.020988 0.012847
-/gps/hist/point 1.038440 0.012644
-/gps/hist/point 1.055893 0.012443
-/gps/hist/point 1.073346 0.012321
-/gps/hist/point 1.090799 0.012159
-/gps/hist/point 1.108251 0.012022
-/gps/hist/point 1.125704 0.011891
-/gps/hist/point 1.143157 0.011784
-/gps/hist/point 1.160610 0.011684
-/gps/hist/point 1.178063 0.011600
-/gps/hist/point 1.195515 0.011499
-/gps/hist/point 1.212968 0.011423
-/gps/hist/point 1.230421 0.011358
-/gps/hist/point 1.247874 0.011273
-/gps/hist/point 1.265326 0.011216
-/gps/hist/point 1.282779 0.011145
-/gps/hist/point 1.300232 0.011114
-/gps/hist/point 1.317685 0.011054
-/gps/hist/point 1.335137 0.011028
-/gps/hist/point 1.352590 0.010988
-/gps/hist/point 1.370043 0.010938
-/gps/hist/point 1.387496 0.010915
-/gps/hist/point 1.404949 0.010911
-/gps/hist/point 1.422401 0.010852
-/gps/hist/point 1.439854 0.010845
-/gps/hist/point 1.457307 0.010834
-/gps/hist/point 1.474760 0.010807
-/gps/hist/point 1.492212 0.010819
-/gps/hist/point 1.509665 0.010799
-/gps/hist/point 1.527118 0.010811
-/gps/hist/point 1.544571 0.009312
-/gps/hist/point 1.562024 0.008725
-/gps/hist/point 1.570700 0.000000
+/gps/hist/type theta 
+/gps/hist/point	0	0.99453
+/gps/hist/point	0.06863	0.99446
+/gps/hist/point	0.10837	0.99169
+/gps/hist/point	0.15533	0.9889
+/gps/hist/point	0.19506	0.98883
+/gps/hist/point	0.22758	0.98608
+/gps/hist/point	0.26374	0.97791
+/gps/hist/point	0.29626	0.97515
+/gps/hist/point	0.32517	0.97239
+/gps/hist/point	0.35412	0.96424
+/gps/hist/point	0.39027	0.95877
+/gps/hist/point	0.42281	0.95331
+/gps/hist/point	0.45898	0.94514
+/gps/hist/point	0.49517	0.93426
+/gps/hist/point	0.52411	0.92611
+/gps/hist/point	0.5603	0.91523
+/gps/hist/point	0.59284	0.90977
+/gps/hist/point	0.62183	0.89621
+/gps/hist/point	0.64718	0.88535
+/gps/hist/point	0.67974	0.87719
+/gps/hist/point	0.71595	0.86361
+/gps/hist/point	0.74855	0.85004
+/gps/hist/point	0.78114	0.83648
+/gps/hist/point	0.81013	0.82291
+/gps/hist/point	0.84273	0.80934
+/gps/hist/point	0.86812	0.79308
+/gps/hist/point	0.90068	0.78492
+/gps/hist/point	0.92607	0.76866
+/gps/hist/point	0.95867	0.75509
+/gps/hist/point	0.98408	0.73613
+/gps/hist/point	1.00585	0.72257
+/gps/hist/point	1.03485	0.70631
+/gps/hist/point	1.05662	0.69276
+/gps/hist/point	1.08201	0.6765
+/gps/hist/point	1.111	0.66293
+/gps/hist/point	1.13643	0.64127
+/gps/hist/point	1.16185	0.62231
+/gps/hist/point	1.18732	0.59524
+/gps/hist/point	1.21633	0.57897
+/gps/hist/point	1.2417	0.56541
+/gps/hist/point	1.26352	0.54375
+/gps/hist/point	1.28898	0.51939
+/gps/hist/point	1.31804	0.49501
+/gps/hist/point	1.3399	0.46795
+/gps/hist/point	1.36897	0.44357
+/gps/hist/point	1.39444	0.4165
+/gps/hist/point	1.41269	0.38944
+/gps/hist/point	1.4309	0.36779
+/gps/hist/point	1.44554	0.34074
+/gps/hist/point	1.46373	0.32179
+/gps/hist/point	1.48196	0.29743
+/gps/hist/point	1.50017	0.27578
+/gps/hist/point	1.51122	0.24603
+/gps/hist/point	1.52584	0.22168
+/gps/hist/point	1.54048	0.19463
+/gps/hist/point	1.54794	0.16218
+/gps/hist/point	1.56262	0.12973
+/gps/hist/point	1.56644	0.09999
+/gps/hist/point	1.56666	0.07026
+/gps/hist/point	1.57052	0.03512
+/gps/hist/point	1.5708	0