From 2c53ff092198be497321259461bf0935a84f54bb Mon Sep 17 00:00:00 2001 From: Jose Soto Date: Mon, 14 Oct 2024 02:36:22 +0200 Subject: [PATCH 1/2] Delete Rn from pdhd generators. Change name of prod_cosmics_1GeV to prod_cosmics, since there is no beam generator --- fcl/protodunehd/gen/prod_beam_cosmics_1GeV_protodunehd.fcl | 1 - ...cosmics_1GeV_protodunehd.fcl => prod_cosmics_protodunehd.fcl} | 1 - 2 files changed, 2 deletions(-) rename fcl/protodunehd/gen/{prod_cosmics_1GeV_protodunehd.fcl => prod_cosmics_protodunehd.fcl} (98%) diff --git a/fcl/protodunehd/gen/prod_beam_cosmics_1GeV_protodunehd.fcl b/fcl/protodunehd/gen/prod_beam_cosmics_1GeV_protodunehd.fcl index 950634ae..2654cbdb 100644 --- a/fcl/protodunehd/gen/prod_beam_cosmics_1GeV_protodunehd.fcl +++ b/fcl/protodunehd/gen/prod_beam_cosmics_1GeV_protodunehd.fcl @@ -41,7 +41,6 @@ physics: ar39: @local::protodunesp_39ar ar42: @local::protodunesp_42ar kr85: @local::protodunesp_85kr - rn222: @local::protodunesp_222rn } diff --git a/fcl/protodunehd/gen/prod_cosmics_1GeV_protodunehd.fcl b/fcl/protodunehd/gen/prod_cosmics_protodunehd.fcl similarity index 98% rename from fcl/protodunehd/gen/prod_cosmics_1GeV_protodunehd.fcl rename to fcl/protodunehd/gen/prod_cosmics_protodunehd.fcl index feba5f09..7b283a0f 100644 --- a/fcl/protodunehd/gen/prod_cosmics_1GeV_protodunehd.fcl +++ b/fcl/protodunehd/gen/prod_cosmics_protodunehd.fcl @@ -39,7 +39,6 @@ physics: ar39: @local::protodunesp_39ar ar42: @local::protodunesp_42ar kr85: @local::protodunesp_85kr - rn222: @local::protodunesp_222rn } From 208166c9bea24e5a1ecd8de43d2e12156591ea0e Mon Sep 17 00:00:00 2001 From: Jose Soto Date: Mon, 14 Oct 2024 02:40:42 +0200 Subject: [PATCH 2/2] Review of the pdhd light simulation workflow. Now only the semianalytical model will have a two separated ionandscint instances. New fhicl added for computable graph. --- .../g4/compgraph_g4_protodunehd.fcl | 71 +++++++++++++++++++ .../g4/compgraph_g4_protodunehd_stage2.fcl | 5 ++ .../semianalytical_g4_protodunehd_stage2.fcl | 6 ++ .../g4/standard_g4_protodunehd.fcl | 6 +- .../g4/standard_g4_protodunehd_stage2.fcl | 2 +- 5 files changed, 85 insertions(+), 5 deletions(-) create mode 100644 fcl/protodunehd/g4/compgraph_g4_protodunehd.fcl create mode 100644 fcl/protodunehd/g4/compgraph_g4_protodunehd_stage2.fcl create mode 100644 fcl/protodunehd/g4/semianalytical_g4_protodunehd_stage2.fcl diff --git a/fcl/protodunehd/g4/compgraph_g4_protodunehd.fcl b/fcl/protodunehd/g4/compgraph_g4_protodunehd.fcl new file mode 100644 index 00000000..4fa08b65 --- /dev/null +++ b/fcl/protodunehd/g4/compgraph_g4_protodunehd.fcl @@ -0,0 +1,71 @@ +#include "services_refactored_pdune.fcl" +#include "LArG4_dune.fcl" +#include "IonAndScint_dune.fcl" +#include "PDFastSim_dune.fcl" + +process_name: G4 + +services: +{ + TFileService: { fileName: "g4_protodunehd_hist.root" } + TimeTracker: {} + MemoryTracker: {} # default is one + RandomNumberGenerator: {} #ART native random number generator + message: @local::standard_info + + @table::protodunehd_refactored_simulation_services + + NuRandomService: @local::dune_prod_seedservice +} + + + +source: +{ + module_type: RootInput + maxEvents: 30000 + fileNames: ["gen_protodunehd.root"] +} + +physics: +{ + + producers: + { + + #retain largeant name for compatibility + largeant: @local::protodune_larg4 + IonAndScint: @local::protodunehd_ionandscint_correlated + PDFastSim: @local::protodune_hd_pdfastsim_ann_ar + rns: {module_type: "RandomNumberSaver"} + } + + analyzers: + { + + } + + simulate: [ rns, largeant, IonAndScint, PDFastSim, PDFastSimExternal ] + + stream1: [ out1 ] + + trigger_paths: [ simulate ] + end_paths: [ stream1 ] +} + +outputs: +{ + out1: + { + module_type: RootOutput + fileName: "%ifb_g4.root" + dataTier: "simulated" + #outputCommands: [ "keep *" ] + outputCommands: [ "keep *", "drop sim::SimEnergyDeposits_largeant_*_G4"] + #fastCloning: false #will fail if the split level is not the same as for the gen stage, so turn it off + compressionLevel: 1 #zlib argument (0-9) + #basketSize: 8192 #[Byte] buffer size at 8k + #splitLevel: 0 #reduces number of buffers + #treeMaxVirtualSize: 1 #[Byte] limits number of buffers/branch to 1 (default is 10) + } +} diff --git a/fcl/protodunehd/g4/compgraph_g4_protodunehd_stage2.fcl b/fcl/protodunehd/g4/compgraph_g4_protodunehd_stage2.fcl new file mode 100644 index 00000000..282698d5 --- /dev/null +++ b/fcl/protodunehd/g4/compgraph_g4_protodunehd_stage2.fcl @@ -0,0 +1,5 @@ +#include "compgraph_g4_protodunehd_stage2.fcl" + +process_name: G4Stage2 + +simulate: [ rns, IonAndScint, PDFastSim ] diff --git a/fcl/protodunehd/g4/semianalytical_g4_protodunehd_stage2.fcl b/fcl/protodunehd/g4/semianalytical_g4_protodunehd_stage2.fcl new file mode 100644 index 00000000..4e727f8e --- /dev/null +++ b/fcl/protodunehd/g4/semianalytical_g4_protodunehd_stage2.fcl @@ -0,0 +1,6 @@ +#include "semianalytical_g4_protodunehd.fcl" + +process_name: G4Stage2 + +physics.simulate: [ rns, IonAndScint, IonAndScintExternal, PDFastSim, PDFastSimExternal ] + diff --git a/fcl/protodunehd/g4/standard_g4_protodunehd.fcl b/fcl/protodunehd/g4/standard_g4_protodunehd.fcl index d5874b86..5a6a8fc1 100644 --- a/fcl/protodunehd/g4/standard_g4_protodunehd.fcl +++ b/fcl/protodunehd/g4/standard_g4_protodunehd.fcl @@ -36,10 +36,8 @@ physics: #retain largeant name for compatibility largeant: @local::protodune_larg4 - IonAndScint: @local::protodunehd_ionandscint - IonAndScintExternal: @local::protodunehd_ionandscint_external + IonAndScint: @local::protodunehd_ionandscint_correlated PDFastSim: @local::protodune_hd_pdfastsim_pvs - PDFastSimExternal: @local::protodune_hd_pdfastsim_pvs_external rns: {module_type: "RandomNumberSaver"} } @@ -48,7 +46,7 @@ physics: } - simulate: [ rns, largeant, IonAndScint, IonAndScintExternal, PDFastSim, PDFastSimExternal ] + simulate: [ rns, largeant, IonAndScint, PDFastSim ] stream1: [ out1 ] diff --git a/fcl/protodunehd/g4/standard_g4_protodunehd_stage2.fcl b/fcl/protodunehd/g4/standard_g4_protodunehd_stage2.fcl index 4b17f06e..2898464a 100644 --- a/fcl/protodunehd/g4/standard_g4_protodunehd_stage2.fcl +++ b/fcl/protodunehd/g4/standard_g4_protodunehd_stage2.fcl @@ -1,6 +1,6 @@ #include "standard_g4_protodunehd.fcl" process_name: G4Stage2 -physics.simulate: [ rns, IonAndScint, IonAndScintExternal, PDFastSim, PDFastSimExternal ] +physics.simulate: [ rns, IonAndScint, PDFastSim ] services.TFileService.fileName: "g4_protodunehd_hist.root" outputs.out1.fileName: "%ifb_g4_stage2.root"