From e7396de7442c6800879122cc0cc0722ea6c22751 Mon Sep 17 00:00:00 2001
From: Matthew Gignac <mgignac@slac.stanford.edu>
Date: Thu, 10 Jul 2025 06:39:06 -0700
Subject: [PATCH 1/3] adding PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
 for v8 MC production

---
 ...sRun2021MCRecon_KF_WithSpacing_pass3.lcsim | 207 ++++++++++++++++++
 1 file changed, 207 insertions(+)
 create mode 100644 steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim

diff --git a/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim b/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
new file mode 100644
index 000000000..e844eeb8c
--- /dev/null
+++ b/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
@@ -0,0 +1,207 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<lcsim xmlns:xs="http://www.w3.org/2001/XMLSchema-instance" xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/lcsim/1.0/lcsim.xsd">
+    <!-- 
+      @brief Steering file that will be used for the reconstruction of 2019 MC. 
+      @author <a href="mailto:omoreno1@ucsc.edu">Omar Moreno</a>
+    -->
+    <execute>
+        <driver name="EventMarkerDriver"/>
+        <!-- Ecal reconstruction drivers -->        
+        <driver name="EcalRawConverter" />
+        <driver name="EcalTimeCorrection"/>
+        <driver name="ReconClusterer" />
+        <driver name="CopyCluster" />
+        <driver name="HodoRunningPedestal"/>
+        <driver name="HodoRawConverter"/>
+        <!-- SVT hit reconstruction drivers -->
+        <driver name="RawTrackerHitSensorSetup"/>
+        <driver name="RawTrackerHitFitterDriver" />
+        <driver name="TrackerHitDriver"/>
+
+        <driver name="KalmanPatRecDriver"/> 
+        <driver name="TrackTruthMatching_KF" />
+	<driver name="ReconParticleDriver_Kalman" />
+	<!-- KFOutputDriver uses too much memory to run during production -->
+	<!--<driver name="KFOutputDriver"/>-->
+
+        <driver name="LCIOWriter"/>
+        <driver name="CleanupDriver"/>
+    </execute>    
+    <drivers>
+        <driver name="EventMarkerDriver" type="org.lcsim.job.EventMarkerDriver">
+            <eventInterval>1000</eventInterval>
+        </driver>
+        <!-- Ecal reconstruction drivers -->
+        <driver name="EcalRawConverter" type="org.hps.recon.ecal.EcalRawConverter2Driver"/>
+        <driver name="EcalTimeCorrection" type="org.hps.recon.ecal.EcalTimeCorrectionDriver"/>
+        <driver name="ReconClusterer" type="org.hps.recon.ecal.cluster.ReconClusterDriver">
+            <logLevel>WARNING</logLevel>
+            <outputClusterCollectionName>EcalClusters</outputClusterCollectionName>
+        </driver>
+        <driver name="CopyCluster" type="org.hps.recon.ecal.cluster.CopyClusterCollectionDriver">
+            <inputCollectionName>EcalClusters</inputCollectionName>
+            <outputCollectionName>EcalClustersCorr</outputCollectionName>
+        </driver>
+        <!-- Hodo reconstruction drivers -->
+        <driver name="HodoRunningPedestal"
+		type="org.hps.recon.ecal.HodoRunningPedestalDriver">
+            <inputCollectionName>HodoscopeReadoutHits</inputCollectionName>
+            <logLevel>CONFIG</logLevel>
+            <isMC>true</isMC>
+        </driver>
+        <driver name="HodoRawConverter" type="org.hps.recon.ecal.HodoRawConverterDriver">
+           <useRunningPedestal>true</useRunningPedestal>
+            <inputCollectionName>HodoscopeReadoutHits</inputCollectionName>
+           <tETAllChannels>8</tETAllChannels>
+           <logLevel>CONFIG</logLevel>
+            <isMC>true</isMC>
+        </driver>
+        <!-- SVT reconstruction drivers -->
+        <!-- Driver used to associate raw tracker hits to corresponding sensor. -->
+        <driver name="RawTrackerHitSensorSetup" type="org.lcsim.recon.tracking.digitization.sisim.config.RawTrackerHitSensorSetup">
+          <readoutCollections>SVTRawTrackerHits</readoutCollections>
+        </driver>
+        <!-- Fit the six raw samples and extract the amplitude and t0. --> 
+        <driver name="RawTrackerHitFitterDriver" type="org.hps.recon.tracking.RawTrackerHitFitterDriver">
+          <chiSqrThresh>.5</chiSqrThresh>
+          <doOldDT>1</doOldDT>
+	  <fitAlgorithm>Pileup</fitAlgorithm>
+          <fitTimeMinimizer>Migrad</fitTimeMinimizer>
+          <!--use this to correct for trigger time in MC instead of subtractTriggerTime-->   
+          <useTimestamps>true</useTimestamps>     
+          <!--offset to get times centered at 0 after timestamp correction-->                 
+	  <!-- set to 0 for unspaced; ~112 for spaced -->
+	  <tsCorrectionScale>120</tsCorrectionScale>
+          <!--correct for the SVT fit time offset...this should be on if <useTimingConditions> is turned on in readout-->        
+          <correctTimeOffset>true</correctTimeOffset>   
+          <!--per sensor shift...set false becasue it's not in readout sim-->       
+          <correctT0Shift>true</correctT0Shift>      
+          <!--use truth time for MC???  typically not used-->             
+          <useTruthTime>false</useTruthTime>    
+          <!--time of flight corrections-->              
+          <subtractTOF>true</subtractTOF>     
+          <!--set this false for MC, true for data-->              
+          <subtractTriggerTime>false</subtractTriggerTime>           
+          <!--per-strip timing correction from database...this should be on i f <useTimingConditions> is turned on in readout  -->  
+          <correctChanT0>false</correctChanT0>          
+          <isMC>true</isMC>
+          <debug>false</debug>
+        </driver>
+        <!-- 
+             Use the fitted raw tracker hits and create 1D clusters using a 
+             nearest neighbor algorithm. 
+        --> 
+        <driver name="TrackerHitDriver" type="org.hps.recon.tracking.DataTrackerHitDriver">
+             <neighborDeltaT>24.0</neighborDeltaT>
+            <neighborDeltaTSigma>3.0</neighborDeltaTSigma>
+            <saveMonsterEvents>false</saveMonsterEvents>
+            <thresholdMonsterEvents>400</thresholdMonsterEvents>
+            <clusterSeedThreshold>4.0</clusterSeedThreshold>
+            <doTimeError>1.0</doTimeError> 
+            <clusterNeighborThreshold>3.0</clusterNeighborThreshold>
+            <clusterThreshold>3.0</clusterThreshold> 
+            <doDeadFix>true</doDeadFix>
+            <doVSplit>true</doVSplit>
+            <debug>false</debug>
+        </driver>
+    
+        <driver name="ReconParticleDriver_Kalman" type="org.hps.recon.particle.HpsReconParticleDriver" > 
+            <ecalClusterCollectionName>EcalClustersCorr</ecalClusterCollectionName>
+            <trackCollectionNames>KalmanFullTracks</trackCollectionNames>
+            <matcherTrackCollectionName>KalmanFullTracks</matcherTrackCollectionName>
+            <trackClusterMatcherAlgo>TrackClusterMatcherMinDistance</trackClusterMatcherAlgo>
+            <unconstrainedV0CandidatesColName>UnconstrainedV0Candidates_KF</unconstrainedV0CandidatesColName>
+            <unconstrainedV0VerticesColName>UnconstrainedV0Vertices_KF</unconstrainedV0VerticesColName>
+            <beamConV0CandidatesColName>BeamspotConstrainedV0Candidates_KF</beamConV0CandidatesColName>
+            <beamConV0VerticesColName>BeamspotConstrainedV0Vertices_KF</beamConV0VerticesColName>
+            <targetConV0CandidatesColName>TargetConstrainedV0Candidates_KF</targetConV0CandidatesColName>
+            <targetConV0VerticesColName>TargetConstrainedV0Vertices_KF</targetConV0VerticesColName>
+            <finalStateParticlesColName>FinalStateParticles_KF</finalStateParticlesColName>
+            <otherElectronsColName>OtherElectrons_KF</otherElectronsColName>
+	    <includeUnmatchedTracksInFSP>true</includeUnmatchedTracksInFSP>		       
+            <beamPositionX> 0.0  </beamPositionX>
+            <beamSigmaX> 0.3 </beamSigmaX>
+            <beamPositionY> 0.0 </beamPositionY>
+            <beamSigmaY> 0.02 </beamSigmaY>
+            <beamPositionZ> -0.5 </beamPositionZ>
+            <trackClusterTimeOffset>25</trackClusterTimeOffset>
+            <maxMatchDt>40</maxMatchDt>
+            <useInternalVertexXYPositions>false</useInternalVertexXYPositions>
+            <minVertexChisqProb> 0.0 </minVertexChisqProb>
+            <maxVertexClusterDt>40.0</maxVertexClusterDt>
+	    <maxElectronP>7.0</maxElectronP>
+            <maxVertexP>7.0</maxVertexP>
+            <requireClustersForV0>false</requireClustersForV0>
+	    <useCorrectedClusterPositionsForMatching>false</useCorrectedClusterPositionsForMatching>
+            <applyClusterCorrections>true</applyClusterCorrections>
+            <useTrackPositionForClusterCorrection>true</useTrackPositionForClusterCorrection>  
+	    <isMC>true</isMC>
+	    <debug>false</debug>
+	    <makeMollerCols>true</makeMollerCols>
+            <unconstrainedMollerCandidatesColName>UnconstrainedMollerCandidates_KF</unconstrainedMollerCandidatesColName>
+            <unconstrainedMollerVerticesColName>UnconstrainedMollerVertices_KF</unconstrainedMollerVerticesColName>
+            <beamConMollerCandidatesColName>BeamspotConstrainedMollerCandidates_KF</beamConMollerCandidatesColName>
+            <beamConMollerVerticesColName>BeamspotConstrainedMollerVertices_KF</beamConMollerVerticesColName>
+            <targetConMollerCandidatesColName>TargetConstrainedMollerCandidates_KF</targetConMollerCandidatesColName>
+            <targetConMollerVerticesColName>TargetConstrainedMollerVertices_KF</targetConMollerVerticesColName>	    
+        </driver>         
+        <driver name="GBLRefitterDriver" type="org.hps.recon.tracking.gbl.GBLRefitterDriver"/>        
+        <driver name="KalmanPatRecDriver" type="org.hps.recon.tracking.kalman.KalmanPatRecDriver">
+	    <numPatRecIteration> 3 </numPatRecIteration>
+            <numKalmanIteration> 1 </numKalmanIteration>
+            <maxPtInverse> 8.757651 </maxPtInverse>
+            <maxD0> 38.0487 </maxD0>
+            <maxZ0> 3.98915 </maxZ0>
+            <maxChi2> 11.777395 </maxChi2>
+            <minHits> 0  </minHits>
+            <minStereo> 3  </minStereo>
+            <maxSharedHits> 3 </maxSharedHits>
+            <maxTimeRange> 39.95028 </maxTimeRange>
+            <maxTanLambda> 8.186345 </maxTanLambda>
+<!--            <maxResidual> 13.71568 </maxResidual>  -->
+            <maxChi2Inc> 13.52662 </maxChi2Inc>
+            <minChi2IncBad> 7.00678 </minChi2IncBad>
+<!--            <maxResidShare> 13.967129 </maxResidShare> -->
+            <maxChi2IncShare> 9.771546584 </maxChi2IncShare>
+            <mxChi2Vtx> 1.7652935 </mxChi2Vtx>
+            <numEvtPlots> 5 </numEvtPlots>
+            <doDebugPlots> false </doDebugPlots>
+            <siHitsLimit> 466 </siHitsLimit>
+            <seedCompThr> .725912 </seedCompThr>
+            <!--numStrategyIter1></numStrategyIter1 tthe mxChi2Vtx was 1.0-->
+            <beamPositionZ> 0.9 </beamPositionZ>
+            <beamSigmaZ> 0.02 </beamSigmaZ>
+            <beamPositionX> 0.0 </beamPositionX>
+            <beamSigmaX> 0.05 </beamSigmaX>
+            <beamPositionY> 0.0 </beamPositionY>
+            <beamSigmaY> 1.0 </beamSigmaY>
+            <lowPhThresh> 7.204329 </lowPhThresh>
+            <verbose> false </verbose>
+        </driver>
+        <driver name="TrackTruthMatching_KF" type="org.hps.analysis.MC.TrackToMCParticleRelationsDriver">
+            <trackCollectionName>KalmanFullTracks</trackCollectionName>
+            <kalmanTracks>true</kalmanTracks>
+            <debug>false</debug>
+        </driver>
+      
+        <driver name="LCIOWriter" type="org.lcsim.util.loop.LCIODriver">
+            <outputFilePath>${outputFile}.slcio</outputFilePath>
+        </driver>
+        <driver name="CleanupDriver" type="org.lcsim.recon.tracking.digitization.sisim.config.ReadoutCleanupDriver"/>
+
+	<driver name="KFOutputDriver" type="org.hps.recon.tracking.kalman.KFOutputDriver">
+          <outputPlotsFilename>${outputFile}.root</outputPlotsFilename>
+          <debug>false</debug>
+          <bsZ>0.9</bsZ>
+          <!--<trackCollectionName>GBLTracks</trackCollectionName> --> 
+          <trackCollectionName>KalmanFullTracks</trackCollectionName>
+          <minMom>0.1</minMom>
+          <maxMom>4.8</maxMom>
+          <chi2Cut>9999</chi2Cut>
+          <doKFresiduals>true</doKFresiduals>
+          <useParticles>true</useParticles>
+        </driver>
+
+
+      </drivers>
+</lcsim>

From b5f8a0908d23155f1974345af774d40b2ed56bc0 Mon Sep 17 00:00:00 2001
From: Matthew Gignac <mgignac@slac.stanford.edu>
Date: Thu, 10 Jul 2025 06:40:17 -0700
Subject: [PATCH 2/3] update description

---
 .../recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim     | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim b/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
index e844eeb8c..f1364cab6 100644
--- a/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
+++ b/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
@@ -1,8 +1,7 @@
 <?xml version="1.0" encoding="UTF-8"?>
 <lcsim xmlns:xs="http://www.w3.org/2001/XMLSchema-instance" xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/lcsim/1.0/lcsim.xsd">
     <!-- 
-      @brief Steering file that will be used for the reconstruction of 2019 MC. 
-      @author <a href="mailto:omoreno1@ucsc.edu">Omar Moreno</a>
+      @brief Steering file that will be used for the reconstruction of 2021 MC production with v8 alignment. 
     -->
     <execute>
         <driver name="EventMarkerDriver"/>

From 6427d03f2faeb50881a67a800a042a7ea3a6f17b Mon Sep 17 00:00:00 2001
From: Matthew Gignac <mgignac@slac.stanford.edu>
Date: Thu, 10 Jul 2025 08:12:28 -0700
Subject: [PATCH 3/3] also update KalmanPatRecDriver beam positions

---
 .../PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim      | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim b/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
index f1364cab6..137a1d612 100644
--- a/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
+++ b/steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass3.lcsim
@@ -168,12 +168,12 @@
             <siHitsLimit> 466 </siHitsLimit>
             <seedCompThr> .725912 </seedCompThr>
             <!--numStrategyIter1></numStrategyIter1 tthe mxChi2Vtx was 1.0-->
-            <beamPositionZ> 0.9 </beamPositionZ>
-            <beamSigmaZ> 0.02 </beamSigmaZ>
+            <beamPositionZ> -0.5 </beamPositionZ>
+            <beamSigmaZ> 0.5 </beamSigmaZ>
             <beamPositionX> 0.0 </beamPositionX>
-            <beamSigmaX> 0.05 </beamSigmaX>
+            <beamSigmaX> 0.3 </beamSigmaX>
             <beamPositionY> 0.0 </beamPositionY>
-            <beamSigmaY> 1.0 </beamSigmaY>
+            <beamSigmaY> 0.02 </beamSigmaY>
             <lowPhThresh> 7.204329 </lowPhThresh>
             <verbose> false </verbose>
         </driver>