diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/art_gen_llp_nupi0_m200_20000.root b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/art_gen_llp_nupi0_m200_20000.root new file mode 100644 index 000000000..08e26ce80 Binary files /dev/null and b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/art_gen_llp_nupi0_m200_20000.root differ diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/hist_llp_nupi0_gen_m200_20000.root b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/hist_llp_nupi0_gen_m200_20000.root new file mode 100644 index 000000000..4aba1a55d Binary files /dev/null and b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/hist_llp_nupi0_gen_m200_20000.root differ diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/llp_config.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/llp_config.fcl new file mode 100644 index 000000000..ad820bb74 --- /dev/null +++ b/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 diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_llp_nupi0_m_200_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_llp_nupi0_m_200_TPC_sbnd.fcl new file mode 100644 index 000000000..3ecd0fe08 --- /dev/null +++ b/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" diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_llp_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_llp_sbnd.fcl new file mode 100644 index 000000000..81ae88332 --- /dev/null +++ b/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"