Skip to content

Update track chi2 skimming cut #1098

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
mgignac merged 5 commits into from
Apr 24, 2025
Merged

Update track chi2 skimming cut #1098

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
ca499e6
also remove collections for skimmed output
mgignac Apr 3, 2025
f42222b
track chi2 cut divide by nDOF
mgignac Apr 22, 2025
953d965
update track chi2/nDOF cut to 30
mgignac Apr 22, 2025
3646c62
fix typo
mgignac Apr 22, 2025
e7d0b0e
fix typo again
mgignac Apr 22, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 3 additions & 3 deletions recon/src/main/java/org/hps/recon/skims/V0Skimmer.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -28,7 +28,7 @@ public class V0Skimmer extends Skimmer {
//private double _clusterTimingCut = 20.0; // only used if _tight is true
private double _posClusterEnergy = 0.2; //GeV
private double _v0Chi2Cut = 100.0;
private double _trackChi2Cut = 80.0;
private double _trackChi2Cut = 30.0;
private double _trackDtCut = 20.0; // the 2-track time difference
private double _trackPMax = 4.5; //GeV
private double _elePMax = 9999; //GeV
Expand Down Expand Up @@ -73,8 +73,8 @@ public boolean passSelection(EventHeader event){
if(_debug)System.out.println(this.getClass().getName()+":: failed nHitsMin "+electron.getTracks().get(0).getTrackerHits().size()+" "+positron.getTracks().get(0).getTrackerHits().size()+" nHitsMin = "+_nHitsMin);
continue;
}
if (electron.getTracks().get(0).getChi2() > _trackChi2Cut
|| positron.getTracks().get(0).getChi2() > _trackChi2Cut) {
if ((electron.getTracks().get(0).getChi2()/electron.getTracks().get(0).getNDF()) > _trackChi2Cut
|| (positron.getTracks().get(0).getChi2()/positron.getTracks().get(0).getNDF()) > _trackChi2Cut) {
if(_debug)System.out.println(this.getClass().getName()+":: failed track chi2");
continue;
}
Expand Down
2 changes: 1 addition & 1 deletion recon/src/main/resources/org/hps/recon/skims/v0skim_parameters_ver0.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
V0CandidateCollectionName UnconstrainedV0Candidates_KF
trackChi2Cut 80.0
trackChi2Cut 30.0
nHitsMin 9
v0Chi2Cut 30.0
trackDtCut 20.0
Expand Down
2 changes: 1 addition & 1 deletion recon/src/main/resources/org/hps/recon/skims/v0skim_parameters_ver1.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
V0CandidateCollectionName UnconstrainedV0Candidates_KF
trackChi2Cut 80.0
trackChi2Cut 30.0
nHitsMin 9
v0Chi2Cut 30.0
trackDtCut 20.0
Expand Down
1 change: 1 addition & 0 deletions ...-files/src/main/resources/org/hps/steering/recon/PhysicsRun2021_pass2_recon_skimmed.lcsim
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -176,6 +176,7 @@
<skimMoller>false</skimMoller>
<v0ParamFile>v0skim_parameters_ver1.txt</v0ParamFile>
<v0OutputFile>${outputFile}_v0skim.slcio</v0OutputFile>
<ignoreCollections>FPGAData HelicalTrackHitRelations HelicalTrackHits HelicalTrackMCRelations KFGBLStripClusterData KFGBLStripClusterDataRelations ReadoutTimestamps RotatedHelicalTrackHitRelations RotatedHelicalTrackHits RotatedHelicalTrackMCRelations SVTFittedRawTrackerHits SVTShapeFitParameters SVTTrueHitRelations StripClusterer_SiTrackerHitStrip1D SVTRawTrackerHits FADCGenericHits HodoReadoutHits HodoCalHits EcalReadoutHits EcalUncalHits HodoGenericClusters VTPBank EcalClusters</ignoreCollections>
</driver>
<driver name="CleanupDriver" type="org.lcsim.recon.tracking.digitization.sisim.config.ReadoutCleanupDriver"/>
<driver name="LCIOWriter" type="org.lcsim.util.loop.LCIODriver">
Expand Down
3 changes: 1 addition & 2 deletions ...c/main/resources/org/hps/steering/recon/PhysicsRun2021_pass2_recon_skimmed_dataqual.lcsim
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -186,9 +186,8 @@
<skimMoller>false</skimMoller>
<v0ParamFile>v0skim_parameters_ver1.txt</v0ParamFile>
<v0OutputFile>${outputFile}_v0skim.slcio</v0OutputFile>
<ignoreCollections>FPGAData HelicalTrackHitRelations HelicalTrackHits HelicalTrackMCRelations KFGBLStripClusterData KFGBLStripClusterDataRelations ReadoutTimestamps RotatedHelicalTrackHitRelations RotatedHelicalTrackHits RotatedHelicalTrackMCRelations SVTFittedRawTrackerHits SVTShapeFitParameters SVTTrueHitRelations StripClusterer_SiTrackerHitStrip1D SVTRawTrackerHits FADCGenericHits HodoReadoutHits HodoCalHits EcalReadoutHits EcalUncalHits HodoGenericClusters VTPBank EcalClusters</ignoreCollections>
</driver>


<driver name="EcalMonitoring" type="org.hps.analysis.dataquality.EcalMonitoring">
<triggerType>all</triggerType>
<clusterCollectionName>EcalClustersCorr</clusterCollectionName>
Expand Down