MeVPrtl LLP (fcls and files)

sbndcode/JobConfigurations/standard/gen/MeVPrtl/llp_config.fcl

@@ -0,0 +1,71 @@
+#include "bnb_kaon_sbnd.fcl"
+BEGIN_PROLOG
+
+llpM: 0.265
+
+#FLUX CONFIGS
+kaon2llp: {
+  tool_type: Kaon2LLPFlux
+  M: @local::llpM # GeV/c^2
+  Target2Absorber: 5000 #cm decay horn length
+  MaxImportanceWeight: 1 #BNB gsimple flat weight at 1
+
+  Beam2DetectorRotation: @local::sbnd_bnb_beam2detector_rotation
+  BeamOrigin: @local::sbnd_bnb_beam_origin
+
+  // timing config
+  # Bunch sigma taken from 
+  # https://beamdocs.fnal.gov/AD/DocDB/0050/005000/001/bunchLength_1st_draft.pdf
+  # Bunch spacing taken from
+  # https://inspirehep.net/files/610a942fd8632bbbca2c8ad90da86670
+  SpillTimeConfig:  "evgb::EvtTimeFNALBeam  booster,  dtbucket=18.936, sigma=1.308"
+  GlobalTimeOffset: 0
+  Verbose: false
+}
+
+# RAY TRACE CONFIGS
+mcraytrace: @local::rethrow_ray_trace_box
+mcraytrace.ReferencePrtlMass: @local::llpM # GeV
+mcraytrace.ReferenceScndPDG: 13 # muon
+mcraytrace.NThrows: 250 
+
+wgtraytrace: @local::weighted_ray_trace_box
+wgtraytrace.ReferencePrtlMass: @local::llpM # GeV
+wgtraytrace.ReferenceScndPDG: 13 # muon
+
+mxdraytrace: @local::mixedweight_ray_trace_box
+mxdraytrace.ReferencePrtlMass: @local::llpM # GeV
+mxdraytrace.ReferenceScndPDG: 13 # muon
+
+
+#DECAY CONFIGS
+decay_llp: {
+  tool_type: LLPMakeDecay
+  Decays: ["nu_pi0"]
+  Majorana: true
+  DecayIsThreeBodyAnisotropic: false
+  ReferenceUE4: 0
+  ReferenceUT4: 0
+  ReferenceUM4: 0
+  ReferenceLLPMass: @local::llpM
+  ReferenceRayLength: 7.549834e2 # sqrt(57) = 7.549834 = diagonal length of detector
+  ReferenceRayDistance: 110e2 #110m = distance from beam target to detector
+  ReferenceLLPEnergyFromKaonEnergy: @local::kaon_energy
+  WidthDecays: ["mu_pi", "e_pi", "nu_mu_mu", "nu_e_e", "nu_nu_nu", "nu_pi0", "nu_eta", "nu_etap", "nu_rho0", "nu_mu_e"]
+  Verbose: false
+}
+
+#FULL GEN CONFIGS
+llp_gen: {
+  module_type: MeVPrtlGen
+  Deweight: false
+  Produce: true
+  AnaOutput: true
+  Verbose: false
+  MesonGen: @local::bnb_kaon
+  Flux: @local::kaon2llp
+  RayTrace: @local::mxdraytrace
+  Decay: @local::decay_llp
+}
+
+END_PROLOG

sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_llp_nupi0_m_200_TPC_sbnd.fcl

@@ -0,0 +1,13 @@
+#include "prodMeVPrtl_llp_sbnd.fcl"
+
+###----------Set mass of LLP [GeV]----------###
+llpM: 0.200
+physics.producers.generator.Flux.M: @local::llpM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::llpM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::llpM
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"

sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_llp_sbnd.fcl

@@ -0,0 +1,53 @@
+#include "llp_config.fcl"
+
+# service configuration
+#include "simulationservices_sbnd.fcl"
+
+
+process_name: Gen
+
+services:
+{
+  @table::sbnd_simulation_services
+  TFileService: { fileName: "hists_prodMeVPrtl_LLP_sbnd_%p-%tc.root" }
+  IFDH: {}
+}
+
+source: {
+  module_type: "EmptyEvent"
+  firstEvent: 1
+  firstRun: 1
+  timestampPlugin: {
+    plugin_type: "GeneratedEventTimestamp"
+  }
+}
+
+physics:
+{
+  producers:
+  {
+    rns:       { module_type: "RandomNumberSaver" }
+    generator: @local::llp_gen
+  }
+
+  simulate:  [ rns, generator ]
+  stream: [ out ]
+
+  trigger_paths: [ simulate ]
+  end_paths:     [ stream ] 
+
+} # physics
+
+outputs: 
+{
+  out:
+  {
+    module_type: RootOutput
+    fileName: "prodMeVPrtl_LLP_sbnd_%p-%tc.root" # default file name, can override from command line with -o or --output
+    dataTier: "simulated"
+    compressionLevel: 1
+    saveMemoryObjectThreshold: 0
+  }
+}
+
+services.NuRandomService.policy: "random"