Skip to content

Commit 9a2cb54

Browse files
EBerzinEBerzin
committed
Adding smearing and killing constants, hit-based smearing
1 parent e55181b commit 9a2cb54

20 files changed

Lines changed: 2572 additions & 0 deletions

‎steering-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021MCRecon_KF_WithSpacing_pass5_physics_HitSmear_HitKill_14272.lcsim‎

Lines changed: 221 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,221 @@
1+
<?xml version="1.0" encoding="UTF-8"?>
2+
<lcsim xmlns:xs="http://www.w3.org/2001/XMLSchema-instance" xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/lcsim/1.0/lcsim.xsd">
3+
<!--
4+
Steering file for pass-5 2021 reconstruction on MC readout data
5+
created: 12/18/2025
6+
@author Matthew Gignac <mgignac@slac.stanford.edu>, Sarah Gaiser <sgaiser@slac.stanford.edu>
7+
-->
8+
<execute>
9+
10+
<driver name="EventMarkerDriver"/>
11+
12+
<!-- Ecal reconstruction drivers -->
13+
<driver name="EcalRawConverter" />
14+
<driver name="EcalTimeCorrection"/>
15+
<driver name="ReconClusterer" />
16+
<driver name="CopyCluster" />
17+
18+
<!-- Hodoscope drivers -->
19+
<driver name="HodoRunningPedestal"/>
20+
<driver name="HodoRawConverter"/>
21+
22+
<!-- SVT reconstruction drivers -->
23+
<driver name="RawTrackerHitSensorSetup"/>
24+
<driver name="RawTrackerHitFitterDriver" />
25+
<driver name="RawHitTimeSmearer"/>
26+
<driver name="TrackerHitDriver"/>
27+
<driver name="StripHitKiller"/>
28+
29+
30+
<driver name="StripHitNHitsSmearer"/>
31+
32+
33+
<driver name="KalmanPatRecDriver"/>
34+
<driver name="TrackTimeSmearer"/>
35+
<driver name="TrackTruthMatching_KF"/>
36+
<driver name="ReconParticleDriver_Kalman" />
37+
<driver name="SVTHitEfficiency"/>
38+
39+
40+
41+
<!-- Event filtering -->
42+
<driver name="LCIOWriter"/>
43+
<driver name="AidaSaveDriver"/>
44+
<driver name="CleanupDriver"/>
45+
46+
</execute>
47+
<drivers>
48+
<driver name="SVTHitEfficiency" type="org.hps.recon.tracking.kalman.SvtHitEfficiencyKalman">
49+
<maxLayer>14</maxLayer>
50+
<useTrkTimeCut>false</useTrkTimeCut>
51+
<!--<trkTimeCut>10.0</trkTimeCut> -->
52+
<!-- <trkTimeMean>15.0</trkTimeMean> -->
53+
<debug>false</debug>
54+
</driver>
55+
<driver name="StripHitKiller" type="org.hps.recon.tracking.StripHitKiller">
56+
<ratioFiles>pass5v9run14272_L1b_axial_hole.txt pass5v9run14272_L1b_stereo_hole.txt pass5v9run14272_L1t_axial_hole.txt pass5v9run14272_L1t_stereo_hole.txt pass5v9run14272_L2b_axial_hole.txt pass5v9run14272_L2b_stereo_hole.txt pass5v9run14272_L2t_axial_hole.txt pass5v9run14272_L2t_stereo_hole.txt</ratioFiles>
57+
<debug>false</debug>
58+
</driver>
59+
<driver name="RawHitTimeSmearer" type="org.hps.recon.tracking.RawHitTimeSmearer">
60+
<smearTimeFile>timeSmearing-2nsL4L7-4nsL1L3.txt</smearTimeFile>
61+
<debug>false</debug>
62+
</driver>
63+
<driver name="TrackTimeSmearer" type="org.hps.recon.tracking.TrackTimeSmearer">
64+
<debug>false</debug>
65+
<smearBottom>0.8</smearBottom>
66+
<smearTop>0.8</smearTop>
67+
</driver>
68+
<driver name="StripHitNHitsSmearer" type="org.hps.recon.tracking.StripHitNHitsSmearer">
69+
<smearPositionFile>/sdf/group/hps/users/eberzin/src/hps-java/tracking/src/main/resources/org/hps/recon/tracking/svtTimeAndPositionSmearing/posSmearing_pass5_v9.txt</smearPositionFile>
70+
<debug>false</debug>
71+
</driver>
72+
<driver name="EventMarkerDriver" type="org.lcsim.job.EventMarkerDriver">
73+
<eventInterval>1000</eventInterval>
74+
</driver>
75+
<driver name="HodoRunningPedestal" type="org.hps.recon.ecal.HodoRunningPedestalDriver">
76+
<inputCollectionName>HodoscopeReadoutHits</inputCollectionName>
77+
<logLevel>CONFIG</logLevel>
78+
<isMC>true</isMC>
79+
</driver>
80+
<driver name="HodoRawConverter" type="org.hps.recon.ecal.HodoRawConverterDriver">
81+
<inputCollectionName>HodoscopeReadoutHits</inputCollectionName>
82+
<tETAllChannels>8</tETAllChannels>
83+
<logLevel>CONFIG</logLevel>
84+
<isMC>true</isMC>
85+
</driver>
86+
87+
<!-- Ecal reconstruction drivers -->
88+
<driver name="EcalRawConverter" type="org.hps.recon.ecal.EcalRawConverter2Driver">
89+
</driver>
90+
<driver name="EcalTimeCorrection" type="org.hps.recon.ecal.EcalTimeCorrectionDriver">
91+
<additionalShift>13.3</additionalShift>
92+
</driver>
93+
<driver name="ReconClusterer" type="org.hps.recon.ecal.cluster.ReconClusterDriver">
94+
<logLevel>WARNING</logLevel>
95+
<outputClusterCollectionName>EcalClusters</outputClusterCollectionName>
96+
</driver>
97+
<driver name="CopyCluster" type="org.hps.recon.ecal.cluster.CopyClusterCollectionDriver">
98+
<inputCollectionName>EcalClusters</inputCollectionName>
99+
<outputCollectionName>EcalClustersCorr</outputCollectionName>
100+
</driver>
101+
<!-- SVT reconstruction drivers -->
102+
<driver name="RawTrackerHitSensorSetup" type="org.lcsim.recon.tracking.digitization.sisim.config.RawTrackerHitSensorSetup">
103+
<readoutCollections>SVTRawTrackerHits</readoutCollections>
104+
</driver>
105+
<driver name="RawTrackerHitFitterDriver" type="org.hps.recon.tracking.RawTrackerHitFitterDriver">
106+
<chiSqrThresh>.5</chiSqrThresh>
107+
<doOldDT>1</doOldDT>
108+
<fitAlgorithm>Pileup</fitAlgorithm>
109+
<fitTimeMinimizer>Migrad</fitTimeMinimizer>
110+
<useTimestamps>true</useTimestamps>
111+
<tsCorrectionScale>165</tsCorrectionScale>
112+
<correctTimeOffset>true</correctTimeOffset>
113+
<correctT0Shift>true</correctT0Shift>
114+
<useTruthTime>false</useTruthTime>
115+
<subtractTOF>true</subtractTOF>
116+
<subtractTriggerTime>false</subtractTriggerTime>
117+
<correctChanT0>false</correctChanT0>
118+
<isMC>true</isMC>
119+
<debug>false</debug>
120+
</driver>
121+
<driver name="TrackerHitDriver" type="org.hps.recon.tracking.DataTrackerHitDriver">
122+
<neighborDeltaT>24.0</neighborDeltaT>
123+
<neighborDeltaTSigma>3.0</neighborDeltaTSigma>
124+
<saveMonsterEvents>false</saveMonsterEvents>
125+
<thresholdMonsterEvents>400</thresholdMonsterEvents>
126+
<clusterSeedThreshold>4.0</clusterSeedThreshold>
127+
<doTimeError>1.0</doTimeError>
128+
<clusterNeighborThreshold>3.0</clusterNeighborThreshold>
129+
<clusterThreshold>3.0</clusterThreshold>
130+
<doDeadFix>true</doDeadFix>
131+
<doVSplit>true</doVSplit>
132+
<debug>false</debug>
133+
</driver>
134+
<driver name="KalmanPatRecDriver" type="org.hps.recon.tracking.kalman.KalmanPatRecDriver">
135+
<addResiduals>true</addResiduals>
136+
<addKinks>true</addKinks>
137+
<numPatRecIteration> 2 </numPatRecIteration>
138+
<numKalmanIteration> 1 </numKalmanIteration>
139+
<maxPtInverse> 8.881915381218574 </maxPtInverse>
140+
<maxD0> 39.249197341740356 </maxD0>
141+
<maxZ0> 8.37847994612359 </maxZ0>
142+
<maxChi2> 11.546843987796496 </maxChi2>
143+
<minHitsTopIter1> 9 </minHitsTopIter1>
144+
<minHitsBotIter1> 9 </minHitsBotIter1>
145+
<minHitsTopIter2> 8 </minHitsTopIter2>
146+
<minHitsBotIter2> 8 </minHitsBotIter2>
147+
<minStereo> 3 </minStereo>
148+
<maxSharedHits> 2 </maxSharedHits>
149+
<maxTimeRange> 40.0 </maxTimeRange>
150+
<maxTanLambda> 5.0 </maxTanLambda>
151+
<maxChi2Inc> 12.320066328390354 </maxChi2Inc>
152+
<minChi2IncBad> 9.206482863412027 </minChi2IncBad>
153+
<maxChi2IncShare> 5.862027198856136 </maxChi2IncShare>
154+
<mxChi2Vtx> 5.508828061070076 </mxChi2Vtx>
155+
<numEvtPlots> 5 </numEvtPlots>
156+
<doDebugPlots> false </doDebugPlots>
157+
<siHitsLimit> 400 </siHitsLimit>
158+
<seedCompThr> 0.3473319986601534 </seedCompThr>
159+
<useBeamPositionConditions>true</useBeamPositionConditions>
160+
<beamSigmaX>0.055</beamSigmaX>
161+
<beamSigmaY>0.045</beamSigmaY>
162+
<useFixedVertexZPosition>true</useFixedVertexZPosition>
163+
<beamPositionZ>-1.1</beamPositionZ>
164+
<lowPhThresh> 7.197594353612735 </lowPhThresh>
165+
<verbose> false </verbose>
166+
</driver>
167+
<driver name="TrackTruthMatching_KF" type="org.hps.analysis.MC.TrackToMCParticleRelationsDriver">
168+
<trackCollectionName>KalmanFullTracks</trackCollectionName>
169+
<kalmanTracks>true</kalmanTracks>
170+
<debug>false</debug>
171+
</driver>
172+
<driver name="ReconParticleDriver_Kalman" type="org.hps.recon.particle.HpsReconParticleDriver" >
173+
<ecalClusterCollectionName>EcalClustersCorr</ecalClusterCollectionName>
174+
<trackCollectionNames>KalmanFullTracks</trackCollectionNames>
175+
<matcherTrackCollectionName>KalmanFullTracks</matcherTrackCollectionName>
176+
<trackClusterMatcherAlgo>TrackClusterMatcherMinDistance</trackClusterMatcherAlgo>
177+
<unconstrainedV0CandidatesColName>UnconstrainedV0Candidates_KF</unconstrainedV0CandidatesColName>
178+
<unconstrainedV0VerticesColName>UnconstrainedV0Vertices_KF</unconstrainedV0VerticesColName>
179+
<beamConV0CandidatesColName>BeamspotConstrainedV0Candidates_KF</beamConV0CandidatesColName>
180+
<beamConV0VerticesColName>BeamspotConstrainedV0Vertices_KF</beamConV0VerticesColName>
181+
<targetConV0CandidatesColName>TargetConstrainedV0Candidates_KF</targetConV0CandidatesColName>
182+
<targetConV0VerticesColName>TargetConstrainedV0Vertices_KF</targetConV0VerticesColName>
183+
<finalStateParticlesColName>FinalStateParticles_KF</finalStateParticlesColName>
184+
<otherElectronsColName>OtherElectrons_KF</otherElectronsColName>
185+
<includeUnmatchedTracksInFSP>true</includeUnmatchedTracksInFSP>
186+
<useInternalVertexXYPositions>false</useInternalVertexXYPositions>
187+
<useBeamPositionConditions>true</useBeamPositionConditions>
188+
<requireClustersForV0>false</requireClustersForV0>
189+
<beamSigmaX>0.055</beamSigmaX>
190+
<beamSigmaY>0.045</beamSigmaY>
191+
<beamPositionZ>-1.1</beamPositionZ>
192+
<maxElectronP>7.0</maxElectronP>
193+
<maxVertexP>7.0</maxVertexP>
194+
<minVertexChisqProb>0.0</minVertexChisqProb>
195+
<maxVertexClusterDt>40.0</maxVertexClusterDt>
196+
<maxMatchDt>40</maxMatchDt>
197+
<trackClusterTimeOffset>40</trackClusterTimeOffset>
198+
<useCorrectedClusterPositionsForMatching>false</useCorrectedClusterPositionsForMatching>
199+
<applyClusterCorrections>true</applyClusterCorrections>
200+
<useTrackPositionForClusterCorrection>true</useTrackPositionForClusterCorrection>
201+
<debug>false</debug>
202+
<isMC>true</isMC>
203+
<makeMollerCols>true</makeMollerCols>
204+
<unconstrainedMollerCandidatesColName>UnconstrainedMollerCandidates_KF</unconstrainedMollerCandidatesColName>
205+
<unconstrainedMollerVerticesColName>UnconstrainedMollerVertices_KF</unconstrainedMollerVerticesColName>
206+
<beamConMollerCandidatesColName>BeamspotConstrainedMollerCandidates_KF</beamConMollerCandidatesColName>
207+
<beamConMollerVerticesColName>BeamspotConstrainedMollerVertices_KF</beamConMollerVerticesColName>
208+
<targetConMollerCandidatesColName>TargetConstrainedMollerCandidates_KF</targetConMollerCandidatesColName>
209+
<targetConMollerVerticesColName>TargetConstrainedMollerVertices_KF</targetConMollerVerticesColName>
210+
</driver>
211+
<driver name="AidaSaveDriver" type="org.lcsim.job.AidaSaveDriver">
212+
<outputFileName>${outputFile}_hit_eff.root</outputFileName>
213+
</driver>
214+
<driver name="LCIOWriter" type="org.lcsim.util.loop.LCIODriver">
215+
<outputFilePath>${outputFile}.slcio</outputFilePath>
216+
</driver>
217+
<driver name="CleanupDriver" type="org.lcsim.recon.tracking.digitization.sisim.config.ReadoutCleanupDriver"/>
218+
219+
</drivers>
220+
</lcsim>
221+

0 commit comments

Comments
 (0)