Add support for validation

JobConfig/examples/create_mixin_files.sh

@@ -0,0 +1,22 @@
+# You must have previously setup mu2efiletools for this script to run
+dsconf=MDC2018a
+for mixin in \
+oot-CRV-cat-cat \
+neutron-CRV-cat \
+dio-CRV-cat \
+photon-CRV-cat \
+PS-CRV-cut-cat \
+TS-CRV-cat-cat \
+DS-CRV-cut-cat \
+oot-TrkCal-cat-cat \
+neutron-TrkCal-cat \
+dio-TrkCal-cat \
+photon-TrkCal-cat \
+DS-flash-TrkCal-cut-cat \
+proton-TrkCal \
+deuteron-TrkCal ; do
+  outfile="$mixin.txt"
+  dataset="sim.mu2e.$mixin.$dsconf.art"
+  echo "Writing files for $dataset to $outfile"
+  mu2eDatasetFileList $dataset > $outfile
+done

JobConfig/reco/CeRecoVal.fcl

@@ -0,0 +1,9 @@
+# example script for running reconstruction for Valdation of Ceminus and Ceplus
+#include "JobConfig/reco/mcdigis-trig_nocosmic.fcl"
+#include "JobConfig/reco/Validation_epilog.fcl"
+physics.RecoPath : [ cosmicTimeMap, @sequence::Reconstruction.DePath,
+  @sequence::Reconstruction.CrvRecoSequence,
+  @sequence::Reconstruction.RecoMCSequence
+]
+services.TFileService.fileName: "nts.owner.CeRecoVal.version.sequencer.root"
+outputs.Output.fileName: "mcs.owner.CeRecoVal.version.sequencer.art"

JobConfig/reco/CeSimRecoVal.fcl

@@ -0,0 +1,64 @@
+# example script for running reconstruction for Valdation of Ceminus and Ceplus
+# running simulation inline.  This doesn't run compression
+#include "fcl/minimalMessageService.fcl"
+#include "fcl/standardProducers.fcl"
+#include "fcl/standardServices.fcl"
+#include "JobConfig/primary/prolog.fcl"
+#include "JobConfig/reco/prolog.fcl"
+#include "Trigger/fcl/templates.fcl"
+
+source : { module_type :EmptyEvent }
+services: @local::Services.SimAndReco
+physics: {
+  producers : {
+    @table::Primary.producers
+    @table::Reconstruction.producers
+    generate : @local::CeEndpointGun
+    protonBunchIntensity : {
+      module_type: ProtonBunchIntensityFlat
+      mean: 0.0 # No PBI possible for pure signal! 
+      halfWidth : 1.0
+    } 
+  }
+
+  filters : { 
+    @table::Primary.filters
+    @table::Reconstruction.filters
+  }
+
+  analyzers: {
+    @table::Primary.analyzers
+    @table::Reconstruction.analyzers
+    @table::Trigger.analyzers
+  }
+
+  TriggerPath : [
+    @sequence::Primary.TriggerPath,
+    protonBunchIntensity,
+    @sequence::Reconstruction.DePath,
+    @sequence::Reconstruction.CrvRecoSequence,
+    @sequence::Reconstruction.RecoMCSequence
+  ]
+  EndPath : @local::Primary.EndPath
+}
+
+outputs : { 
+  Output : {
+     module_type : RootOutput
+     SelectEvents : [ TriggerPath ]
+     fileName    : @nil
+     # drop StepPoints from calorimeter
+     outputCommands : [ "drop *_*_*_*",
+      @sequence::Reconstruction.HighReco,
+      @sequence::Reconstruction.LowRecoAll,
+      @sequence::Reconstruction.DigiMC
+     ]
+  }
+}
+#include "JobConfig/primary/epilog.fcl"
+#include "JobConfig/reco/epilog.fcl"
+#include "JobConfig/reco/Validation_epilog.fcl"
+
+services.TFileService.fileName: "nts.owner.CeSimRecoVal.version.sequencer.root"
+outputs.Output.fileName: "mcs.owner.CeSimRecoVal.version.sequencer.art"
+

JobConfig/reco/MT_MDC2018.C

@@ -0,0 +1,14 @@
+//
+// Generate random misalignemnts for MDC2018.  Execute as:
+//  root -l JobConfig/reco/MT_MDC2018.C
+//
+#include "JobConfig/reco/MisalignTracker.C+"
+#include "math.h"
+{
+  double twist = 5; // 5 mm at the outer radius over the length of the tracker (~1 mrad)
+  double skew = 5; // 5mm displacement over the length of the tracker
+  double squeeze = 5; // 5mm over the length of the tracker
+  double anglesig = 0.0007; // 0.7 mrad, ~500 um at the outer edge
+  double possig = 0.5; // 500 um precision
+  MisalignTracker(twist,skew,squeeze,anglesig,possig,"MT_MDC2018.txt")
+}