Add mcsteps branch to ntuplehelper

Author: AndrewEdmonds11

helper/ntuplehelper.py

@@ -3,14 +3,14 @@
 class nthelper:
 
     single_object_branches = ['evtinfo', 'evtinfomc', 'hitcount', 'tcnt', 'crvsummary', 'crvsummarymc']
-    vector_object_branches = ['trk', 'trkmc', 'trkcalohit', 'trkcalohitmc', 'caloclusters', 'calohits', 'calorecodigis', 'calodigis', 'crvcoincs', 'crvcoincsmc', 'crvcoincsmcplane', 'trkqual', 'trkpid']
+    vector_object_branches = ['trk', 'trkmc', 'trkcalohit', 'trkcalohitmc', 'caloclusters', 'calohits', 'calorecodigis', 'calodigis', 'crvcoincs', 'crvcoincsmc', 'crvcoincsmcplane', 'trkqual', 'trkpid', 'mcsteps']
     vector_vector_object_branches = ['trksegs', 'trksegpars_lh', 'trksegpars_ch', 'trksegpars_kl', 'trkmcsim', 'trkhits', 'trkhitsmc', 'trkmats', 'trkhitcalibs', 'trkmcsci', 'trkmcssi', 'trksegsmc' ]
 
     evt_branches = ['evtinfo','evtinfomc','hitcount','tcnt']
     trk_branches = ['trk', 'trkmc', 'trkcalohit', 'trkcalohitmc', 'trkqual', 'trkpid']
     trksegs_branches = ['trksegs', 'trksegpars_lh', 'trksegpars_ch', 'trksegpars_kl', 'trksegsmc']
     straw_branches = ['trkhits', 'trkmats', 'trkhitsmc', 'trkhitcalibs']
-    mc_branches = ['trkmcsim']
+    mc_branches = ['trkmcsim', 'mcsteps']
     calo_branches = ['caloclusters', 'calohits', 'calorecodigis', 'calodigis']
     crv_branches = ['crvsummary','crvsummarymc','crvcoincs','crvcoincsmc','crvcoincsmcplane']
     deprecated_branches = ['trkmcsci','trkmcssi']
@@ -58,7 +58,8 @@ class nthelper:
                            "trkqual" : "MVAResultInfo",
                            "trkpid" : "MVAResultInfo",
                            "helices" : "HelixInfo",
-                           "trksegsmc" : "SurfaceStepInfo"
+                           "trksegsmc" : "SurfaceStepInfo",
+                           "mcsteps" : "MCStepInfo"
                           }
 
     #

inc/MCStepInfo.hh

@@ -14,8 +14,8 @@ namespace mu2e {
     float time = 0;  // time of this step WRT MC primary proton (ns)
     float de = 0; // energy deposit through this step (MeV)
     float dp = 0; // momentum magnitude change throw this step (MeV/c)
-    bool early = false;
-    bool late = false; // flag if this is the earliest or latest step
+    bool early = false; // flag if this is the earliest step
+    bool late = false; // flag if this is the latest step
     XYZVectorF mom; // particle momentum at the start of this step
     XYZVectorF pos;  // particle position at the start of this step
     void reset() {*this = MCStepInfo(); }