Turn off vd steps by default but add a fcl file that demonstrates how to use it

Author: AndrewEdmonds11

fcl/README.md

@@ -22,6 +22,7 @@ where ```tier``` is the data tier of the input dataset, ```type``` is the type o
 | from_mcs-mockdata_separateTrkBranches.fcl | mock datasets | example on how to separate the tracks into separate branches again|
 | from_mcs-mockdata_selectorExample.fcl | mock datasets | example on how to use a selector to select certain types of tracks before putting them into the EventNtuple |
 | from_mcs-mixed_trkQualCompare.fcl | reconstructed mixed (i.e. primary+background hits) datasets | shows how to output result of more than one TrkQual |
+| from_mcs-primary_addVDSteps.fcl | reconstructed primary (i.e. no background hits) datasets | shows how to add the branch for virtual detector steps |
 | from_mcs-DeMCalib.fcl | reconstructed primary or mixed datasets | only writes one track per event |
 | from_mcs-OffSpill.fcl | off spill datasets | only contains ```CentralHelix``` tracks (i.e. field-on cosmics) |
 | from_dig-mockdata.fcl | mock datasets (digis) | runs reconstruction and creates EventNtuple in one job |

fcl/from_mcs-primary_addVDSteps.fcl

@@ -0,0 +1,3 @@
+#include "EventNtuple/fcl/from_mcs-primary.fcl"
+
+physics.analyzers.EventNtuple.StepPointMCTags : [ "compressRecoMCs:virtualdetector" ]

fcl/prolog.fcl

@@ -319,7 +319,7 @@ EventNtupleMaker : {
   ExtraMCStepCollectionTags : []
   SurfaceStepCollectionTag : "compressRecoMCs"
   FitType : LoopHelix
-  StepPointMCTags : [ "compressRecoMCs:virtualdetector" ]
+  StepPointMCTags : [ ]
 }
 # instance for processing trigger (digitization) output from simulation
 EventNtupleMakerTTMC: {

validation/test_fcls.sh

@@ -2,6 +2,10 @@
 # test_fcls.sh - runs fcl files to make sure that they complete successfully
 # Note: requires relevant filelists in a directory above this one
 #
+if [ $# -ne 1 ]; then
+    echo "Usage: . ./validation/test_fcls.sh MDC202X"
+    return
+fi
 vMDC=$1
 
 log_file="test_fcls.log"
@@ -137,6 +141,16 @@ else
     echo "FAIL"
 fi
 
+echo -n "from_mcs-primary_addVDSteps.fcl... "
+echo "mu2e -c fcl/from_mcs-primary_addVDSteps.fcl -S ../filelists/$primary_dataset.list --TFileName nts.ntuple.primaryVDSteps.root -n 100" >> ${log_file} 2>&1
+mu2e -c fcl/from_mcs-primary_addVDSteps.fcl -S ../filelists/$primary_dataset.list --TFileName nts.ntuple.primaryVDSteps.root -n 100 >> ${log_file} 2>&1
+if [ $? == 0 ]; then
+    echo "OK"
+else
+    echo "FAIL"
+fi
+
+
 echo -n "from_dig-mockdata.fcl... "
 echo "mu2e -c fcl/from_dig-mockdata.fcl -S ../filelists/${digi_dataset}.list --TFileName nts.ntuple.dig.root -n 100" >> ${log_file} 2>&1
 mu2e -c fcl/from_dig-mockdata.fcl -S ../filelists/${digi_dataset}.list --TFileName nts.ntuple.dig.root -n 100 >> ${log_file} 2>&1