Merge pull request #126 from soleti/nofield

Geometry and configuration files for extracted position and field off

CampaignConfig/mdc2020_flateminus.cfg

@@ -6,7 +6,7 @@ submitter = mu2e
 outdir_fcl = /pnfs/mu2e/persistent/datasets/phy-etc/cnf/mu2e/
 outdir_dig_tape = /pnfs/mu2e/tape/phy-sim/dig/mu2e/
 outdir_dts_tape = /pnfs/mu2e/tape/phy-sim/dts/mu2e/
-outdir_mcs_tape = /pnfs/mu2e/tape/phy-sim/dts/mu2e/
+outdir_mcs_tape = /pnfs/mu2e/tape/phy-sim/mcs/mu2e/
 logdir_bck = /pnfs/mu2e/persistent/datasets/phy-etc/bck/mu2e/
 run_number = 1210
 njobs = 10

JobConfig/cosmic/ExtractedNoField.fcl

@@ -0,0 +1,57 @@
+# This configuration simulate cosmics with no field and extracted position
+# This is analogous to S1DSStops, so it needs to be included by another FCL
+# file that specifies the generator (e.g. physics.producers.generate : @local::Cosmic.generateCRY)
+# Stefano Roberto Soleti, 2021
+
+#include "Offline/fcl/standardServices.fcl"
+#include "Production/JobConfig/common/prolog.fcl"
+#include "Production/JobConfig/primary/prolog.fcl"
+#include "Production/JobConfig/cosmic/prolog.fcl"
+
+process_name : Primary
+
+services: @local::Services.Sim
+
+physics: {
+  producers : {
+    @table::Common.producers
+    @table::Primary.producers
+  }
+
+  filters : {
+    @table::Common.filters
+    @table::Primary.filters
+  }
+
+  analyzers: {
+    @table::Primary.analyzers
+  }
+
+  PrimaryPath :  [ @sequence::Common.generateSequence, @sequence::Common.g4Sequence, @sequence::Primary.PrimarySequence ]
+  EndPath : @local::Primary.EndPath
+  trigger_paths : [ PrimaryPath ]
+  end_paths : [ EndPath ]
+}
+
+outputs : { Output : @local::Primary.Output }
+
+outputs.Output.fileName : "sim.owner.CosmicNoFieldExtracted.version.sequencer.art"
+outputs.Output.outputCommands : [@sequence::Primary.KeptProducts, "keep mu2e::CosmicLivetime_*_*_*", "keep mu2e::GenParticles_*_*_*"]
+
+physics.producers.g4run.TrajectoryControl: @local::mu2eg4DefaultTrajectories
+physics.producers.g4run.Mu2eG4CommonCut: @local::Cosmic.Mu2eG4CommonCutCosmicNoField
+physics.producers.FindMCPrimary.PrimaryGenIds: ["CosmicCRY", "CosmicCORSIKA"]
+
+#
+# final configuration
+#
+#include "Production/JobConfig/common/epilog.fcl"
+#include "Production/JobConfig/primary/epilog.fcl"
+
+physics.end_paths : [ EndPath ]
+# no histogram output
+services.TFileService.fileName : "/dev/null"
+services.SeedService.baseSeed         :  773651
+
+# Need this to increase the world volume.
+services.GeometryService.inputFile: "Production/JobConfig/cosmic/geom_extracted_nofield.txt"

JobConfig/cosmic/ExtractedNoFieldCORSIKA.fcl

@@ -0,0 +1,20 @@
+#include "Production/JobConfig/cosmic/ExtractedNoField.fcl"
+
+source: {
+  module_type: FromCorsikaBinary
+  fileNames: ["/pnfs/mu2e/persistent/users/srsoleti/corsika/sim.srsoleti.corsika.v1.30001.csk"]
+  runNumber          : 1205
+  showerAreaExtension  : 10000
+  projectToTargetBox : true
+  targetBoxXmin: -10000
+  targetBoxXmax: 3000
+  targetBoxYmin: -5000
+  targetBoxYmax: 9880 # defined by CRV
+  targetBoxZmin: 17105 # defined by CRV
+  targetBoxZmax: 32220 # defined by Calo Z1
+  fluxConstant: 1.8e4
+}
+
+physics.producers.generate : @local::Cosmic.generateCORSIKA
+
+outputs.Output.fileName         : "sim.owner.CosmicCORSIKANoFieldExtracted.version.sequencer.art"

JobConfig/cosmic/ExtractedNoFieldCRY.fcl

@@ -0,0 +1,6 @@
+#include "Production/JobConfig/cosmic/ExtractedNoField.fcl"
+
+source.module_type: EmptyEvent
+physics.producers.generate : @local::Cosmic.generateCRY
+physics.producers.generate.inputFile    : "Offline/EventGenerator/defaultConfigs/extractedCRYconfig.txt"
+outputs.Output.fileName : "sim.owner.CosmicCRYNoFieldExtracted.version.sequencer.art"

JobConfig/cosmic/geom_extracted_nofield.txt

@@ -0,0 +1,10 @@
+#include "Offline/Mu2eG4/geom/geom_common_extracted.txt"
+vector<string> bfield.innerMaps = {
+  "BFieldMaps/Mau13/DS-nofield.txt",
+  "BFieldMaps/Mau13/PSMap.header",
+  "BFieldMaps/Mau13/TSuMap.header",
+  "BFieldMaps/Mau13/TSdMap.header",
+  "BFieldMaps/Mau13/PStoDumpAreaMap.header",
+  "BFieldMaps/Mau13/ProtonDumpAreaMap.header",
+  "BFieldMaps/Mau13/DSExtension.header"
+};

JobConfig/cosmic/prolog.fcl

@@ -144,7 +144,7 @@ Cosmic: {
   }
   generateCRY: {
     module_type: CRYEventGenerator
-    inputFile    : "EventGenerator/defaultConfigs/defaultCRYconfig.txt"
+    inputFile    : "Offline/EventGenerator/defaultConfigs/defaultCRYconfig.txt"
   }
   generateCORSIKA: {
     module_type  : CORSIKAEventGenerator
@@ -183,6 +183,17 @@ Cosmic: {
       }
     ]
   }
+  Mu2eG4CommonCutCosmicNoField: {
+    type: union
+    pars:
+    [
+      @local::cosmicKineticEnergyCutS1,
+      {
+        type: inVolume
+        pars: [ worldDirtBottom ]
+      }
+    ]
+  }
 
   Output : {
     module_type : RootOutput