////////////////////////////////////////////////////////////////////
/// Specify inputs for the flasher generator.
/// This file stores the controls for setting the intensity, wavelength,
/// time, angular distribution and crosstalk for the flasher generator
/// as well as tables for the relevant distributions.
/// Contact: J Walker (john.walker@liverpool.ac.uk)
////////////////////////////////////////////////////////////////////
{
type: "FLASHER",
version: 1,
run_range: [0, 0],
pass: 0,
comment: "",
timestamp: "",

////////////////////////////////////////////////////////////////////
distribution_type: "uniform", // "uniform" to uniformly distribute number of photons between a min and max range, "SNO" to replicate the nhit distribution of flashers from SNO

intensity_max: 2000, // Maximum photons per flasher
intensity_min: 1000, // Minimum photons per flasher
pulse_mode: "fixed", // optional - "poisson" to poisson fluctuate number of photons, "fixed" for each pulse to have same intensity photons

// Parameters for replicating SNO flasher distribution:
//
// Nphotons vs Nhits
// Nphotons = exp((Nhits+[2])/[0])-[1]
blackbody_nphot_nhit_param_0: 2382.83
blackbody_nphot_nhit_param_1: 621844.0
blackbody_nphot_nhit_param_2: 31774.5
flat_nphot_nhit_param_0: 2683.29
flat_nphot_nhit_param_1: 799531.0
flat_nphot_nhit_param_2: 36456.7
helium_nphot_nhit_param_0: 2299.57
helium_nphot_nhit_param_1: 803503.0
helium_nphot_nhit_param_2: 31253.0
mono500_nphot_nhit_param_0: 2589.95
mono500_nphot_nhit_param_1: 234141.0
mono500_nphot_nhit_param_2: 32009.2
mono350_nphot_nhit_param_0: 3067.34
mono350_nphot_nhit_param_1: 621319.0
mono350_nphot_nhit_param_2: 40903.0
//
// Fit to SNO flasher distribution:
// 0<x<120: exp([0]+[1]*x)
// 120<x<200: exp([2]+[3]*x)
// 200<x<350: exp([4]+[5]*x)
// 350<x<525: exp([6]+[7]*x)
// 525<x<10000: exp([8]+[9]*x)
// Corresponding cumulative distribution:
// 0<x<120: exp([0]+[1]*x)/[1]
// 120<x<200: exp([2]+[3]*x)/[3]
// 200<x<350: exp([4]+[5]*x)/[5]
// 350<x<525: exp([6]+[7]*x)/[7]
// 525<x<10000: exp([8]+[9]*x)/[9]
sno_dist_param_0: -3.36661962526
sno_dist_param_1: -0.0356662745692
sno_dist_param_2: -5.36225747473
sno_dist_param_3: -0.0187205069641
sno_dist_param_4: -8.06595830939
sno_dist_param_5: -0.00539735674383
sno_dist_param_6: -10.0757754764
sno_dist_param_7: 0.000558351837771
sno_dist_param_8: -2.26162969064
sno_dist_param_9: -0.0139353768637
// Cumulative distribution boundaries:
sno_dist_bound_0: 0.954077928993
sno_dist_bound_1: 0.974651587293
sno_dist_bound_2: 0.985623320539
sno_dist_bound_3: 0.995030211377
// Maximum photons for SNO flasher distribution
sno_max_nphot: 300000

////////////////////////////////////////////////////////////////////
// wavelength settings:
wl_dist: "mono", // Choose "mono" for single wavelength, "LED" for LED distribution around 490nm, "flat" for flat distribution, "helium" for helium emission spectrum, and "blackbody" for blackbody radiation spectrum

mono_wl: 500.0, // Set to desired single wavelength value. To use SNO nhit distribution should be 350.0 or 500.0! Defaults to 500nm photon distribution otherwise.

led_dist_wl: [ 436.8, 439.3, 441.8, 444.3, 446.8, 449.3, 451.8, 454.3, 456.8, 459.3, 461.8, 464.3, 466.8, 469.3, 471.8, 474.3, 476.8, 479.3, 481.8, 484.3, 486.8, 489.3, 491.8, 494.3, 496.8, 499.3, 501.8, 504.3, 506.8, 509.3, 511.8, 514.3, 516.8, 519.3, 521.8, 524.3, 526.8, 529.3, 531.8, 534.3, 536.8, 539.3, 541.8, 544.3, 546.8, 549.3, 551.8, 554.3, 556.8, 559.3, 561.8, 564.3, 566.8, 569.3],
led_dist_wl_intensity: [ 0.0, 0.00000, 0.00149, 0.00295, 0.00564, 0.01037, 0.01837, 0.03131, 0.05137, 0.08113, 0.12334, 0.18048, 0.25422, 0.34468, 0.44984, 0.56510, 0.68333, 0.79536, 0.89112, 0.96103, 0.99764, 1.00000, 0.98480, 0.95537, 0.91301, 0.85952, 0.79710, 0.72820, 0.65533, 0.58097, 0.50736, 0.43648, 0.36990, 0.30881, 0.25396, 0.20574, 0.16420, 0.12908, 0.09997, 0.07627, 0.05732, 0.04243, 0.03095, 0.02223, 0.01573, 0.01097, 0.00753, 0.00510, 0.00340, 0.00223, 0.00144, 0.00000, 0.00000, 0.00000],

flat_dist_wl: [ 200.0, 210.0, 220.0, 230.0, 240.0, 250.0, 260.0, 270.0, 280.0, 290.0, 300.0, 310.0, 320.0, 330.0, 340.0, 350.0, 360.0, 370.0, 380.0, 390.0, 400.0, 410.0, 420.0, 430.0, 440.0, 450.0, 460.0, 470.0, 480.0, 490.0, 500.0, 510.0, 520.0, 530.0, 540.0, 550.0, 560.0, 570.0, 580.0, 590.0, 600.0, 610.0, 620.0, 630.0, 640.0, 650.0, 660.0, 670.0, 680.0, 690.0, 700.0, 710.0, 720.0, 730.0, 740.0, 750.0, 760.0, 770.0, 780.0, 790.0, 800.0],
flat_dist_wl_intensity: [ 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0],

helium_dist_wl: [ 388.0, 447.0, 501.0, 588.0, 668.0, 706.0, 728.0, 778.0],
helium_dist_wl_prob: [ 0.11, 0.06, 0.08, 0.28, 0.14, 0.21, 0.05, 0.07],

blackbody_temp:	6000, // Blackbody temperature in Kelvin


////////////////////////////////////////////////////////////////////
// time settings:
time_dist: "DEFAULT", // Set to "mono" for no time spread or set to "DEFAULT" to use the time distribution below

// time distribution (not used if "mono" specified above, Gaussian centred on 0, width = 1ns.)
dist_time: [-3.75, -3.25, -2.75, -2.25, -1.75, -1.25, -0.75, -0.25, 0.25,  0.75,  1.25,  1.75,  2.25,  2.75,  3.25,  3.75, ],
dist_time_intensity: [0.0002, 0.0010, 0.0046, 0.0159, 0.0432, 0.0913, 0.1506, 0.1933, 0.1933, 0.1506, 0.0913, 0.0432, 0.0159, 0.0046, 0.0010, 0.0002],


////////////////////////////////////////////////////////////////////
// angular settings:
angle_dist: "DEFAULT", // String - set to "iso" for isotropic distribution, "zero" for no dispersion, or "DEFAULT" to use the angular distribution below.

// specify unit "degrees" or "radians"
dist_angle_option: "degrees",
dist_angle: [0.00, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.00, 11.00, 12.00, 13.00, 14.00, 15.00, 16.00, 17.00, 18.00, 19.00, 20.00, 21.00, 22.00, 23.00, 24.00, 25.00, 26.00, 27.00, 28.00, 29.00, 30.00, 31.00, 32.00, 33.00, 34.00, 35.00, 36.00, 37.00, 38.00, 39.00, 40.00, 41.00, 42.00, 43.00, 44.00, 45.00, 46.00, 47.00, 48.00, 49.00, 50.00, 51.00, 52.00, 53.00, 54.00, 55.00, 56.00, 57.00, 58.00, 59.00, 60.00, 61.00, 62.00, 63.00, 64.00, 65.00, 66.00, 67.00, 68.00, 69.00, 70.00, 71.00, 72.00, 73.00, 74.00, 75.00, 76.00, 77.00, 78.00, 79.00, 80.00, ],
dist_angle_intensity: [1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., ],


////////////////////////////////////////////////////////////////////
// Cross talk setting
// Must leave off by default! Should not interfere with other simulations! Turn on in macro.
cross_talk_flag: 0, // Flag to turn on cross talk. Set to 0 (off) as default. Set to 1 to turn on.
}