Merge pull request #706 from sophiemiddleton/CaloCalibGun

Calo Source Calib + CaloCalibGun

CaloMC/fcl/RunCaloCalibGun.fcl

@@ -0,0 +1,95 @@
+// Driver FCL for the calocalibration (Reconstruct and truth)
+#include "fcl/minimalMessageService.fcl"
+#include "fcl/standardProducers.fcl"
+#include "fcl/standardServices.fcl"
+#include "CaloMC/fcl/prolog.fcl"
+#include "CaloReco/fcl/prolog.fcl"
+#include "CaloCluster/fcl/prolog.fcl"
+#include "CommonMC/fcl/prolog.fcl"
+#include "EventGenerator/fcl/prolog.fcl"
+#include "Production/JobConfig/digitize/prolog.fcl"
+process_name: CaloCalibExample
+
+source: { 
+   module_type: EmptyEvent
+   maxEvents : 1000
+}
+
+
+services : @local::Services.SimAndReco
+services.TFileService           : { fileName : "nts.owner.CaloCalib.version.sequencer.root" }
+services.SeedService.baseSeed         :  99
+services.SeedService.maxUniqueEngines :  20
+
+physics: {
+
+    producers: {
+
+      generate: {
+        module_type : CaloCalibGun
+        cosmin :  -1.0
+        cosmax :  1.0
+        phimin :  0.0
+        phimax : 2.0
+        tmin  :  500.
+        tmax : 1694.
+      }
+
+      g4run                : @local::g4run
+
+      protonTimeMap        : { module_type : GenerateProtonTimes }
+      muonTimeMap          : { module_type : GenerateMuonLife }
+      EWMProducer : { module_type : EventWindowMarkerProducer }
+      FindMCPrimary : { module_type : FindMCPrimary }//@local::CommonMC.FindMCPrimary
+      @table::EventGenerator.producers
+      @table::Digitize.producers
+      @table::CommonMC.producers
+      @table::CaloMC.producers
+      @table::CaloReco.producers
+      @table::CaloCluster.producers
+      @table::CaloMC.TruthProducers
+
+   }
+
+    filters: {
+      @table::Digitize.filters
+    }
+
+    analyzers: {
+
+      CaloExample: {
+        module_type                : CaloCalibAna
+        caloHitCollection     : CaloHitMaker
+        caloClusterCollection : CaloClusterMaker
+        caloHitTruth    : CaloHitTruthMatch 
+        vdCollection        : virtualdetector 
+        caloClusterTruth  : CaloClusterTruthMatch
+
+	    }
+
+    }
+
+    p1: [ generate, g4run, protonTimeMap, muonTimeMap, EWMProducer,
+          @sequence::CaloMC.StepSim,
+          @sequence::CommonMC.DigiSim,
+          @sequence::CaloMC.DigiSim,
+          @sequence::CaloReco.Reco,
+          @sequence::CaloCluster.Reco,	 
+	        @sequence::CaloMC.TruthMatch
+       ]
+
+    e1: [CaloExample]
+
+    trigger_paths: [p1]
+    end_paths:     [e1]
+}
+
+physics.producers.CaloHitTruthMatch.primaryParticle : "generate"
+physics.producers.FindMCPrimary.PrimaryProcess : "mu2eCaloCalib"
+
+physics.producers.EWMProducer.SpillType : 0
+physics.producers.muonTimeMap.applyToGenIds : [ "StoppedParticleReactionGun", "dioTail", "ExternalRMC", "InternalRMC",
+                 "MuCapProtonGenTool", "MuCapDeuteronGenTool", "DIOGenTool", "MuCapNeutronGenTool", "MuCapPhotonGenTool", "MuCapGammaRayGenTool" ]
+physics.producers.protonTimeMap.ignoredGenIds : ["unknown", "particleGun", "cosmicToy", "cosmicDYB", "cosmic", "CosmicCRY", "pbarFlat", "CosmicCORSIKA"]
+
+physics.producers.CaloShowerStepFromStepPt.physVolInfoInput : "g4run"

CommonMC/fcl/prolog.fcl

@@ -22,7 +22,6 @@ CommonMC: {
   FindMCPrimary : {
     module_type : FindMCPrimary
     debugLevel : 0
-    GenParticles : "compressDigiMCs"
     SimParticles : "compressDigiMCs"
     PrimaryGenIds : ["CeEndpoint", "CeLeadingLog", "dioTail",
         "cosmicDYB", "CosmicCRY", "CosmicCORSIKA",

EventGenerator/fcl/prolog.fcl

@@ -133,6 +133,16 @@ CaloTBGun : {
     module_type    : CaloTBGun
 }
 
+
+CaloCalibGun : {
+  module_type : CaloCalibGun
+  cosmin :  -1.0
+  cosmax :  1.0
+  phimin :  0.0
+  phimax : 2.0
+  tmin  :  500.
+  tmax : 1694.
+}
 # Save state of random numbers at the start of each event
 # the module is defined in art/Framework/Modules
 
@@ -146,7 +156,7 @@ EventGenerator : {
     pionStops : {
         treeName               : "stoppedPionDumper/stops"
         branchName             : "stops"
-        inputFiles             : ["mergedMuonStops/nts.mu2e.pion-DS-TGTstops.MDC2018a.001002_00000000.root" ]
+        inputFiles             : ["mergedMuonStops/nts.mu2e.pion-DS-TGTstops.MDC2018a.001002_00000000.root" ] # does this need removing?
         averageNumRecordsToUse : 0
         verbosityLevel         : 0
     }