Difference between revisions of "Spring 2017 Analysis Launch Cuts"
From GlueXWiki
(→Energy Cuts) |
(→Timing Cuts) |
||
(68 intermediate revisions by 7 users not shown) | |||
Line 1: | Line 1: | ||
− | + | '''These cuts also apply to 2018 and 2019 GlueX data''' | |
− | == | + | == PreSelect Cuts == |
+ | * For a given DChargedTrack, at least one hypothesis must have a hit in either: BCAL/FCAL/TOF/SC. | ||
+ | * For a given DNeutralShower, it must at least have a shower energy of 100 MeV/c. | ||
+ | * For neutral showers in the BCAL, they must have hits in at least 2 cells. | ||
== Timing Cuts == | == Timing Cuts == | ||
+ | * Timing cuts are applied only on the system with the best timing information available, in the order: BCAL/TOF/FCAL/SC. In other words, if a track has a hit in the TOF, then no cut is applied on the FCAL or SC time. | ||
+ | * Note that some hypotheses may not have a hit in any system and thus survive all of these cuts. | ||
+ | * Also, even if the SC is the best-available detector, no cut is placed on the timing if the track has been matched to more than one SC hit, UNLESS it's the only track in the combo. | ||
+ | ** This is because the SC timing is not really good enough for PID over such a short flight distance, so the cut is primarily used to reject tracks from other beam bunches or ghost tracks. | ||
+ | * Also, all kaons are cut that: | ||
+ | ** Do not have a hit in one of these timing detectors, AND | ||
+ | ** Do not have enough hits in the CDC to compute dE/dx | ||
+ | ** This is because their background is huge already, without letting these combos through. | ||
+ | * Cuts are re-placed after the kinematic fit, using updated momentum and vertex information. | ||
+ | * If a kinfit is not performed, then the timing cut at detached vertices is looser than otherwise: +/- 2ns WIDER for charged tracks, ~0.05ns for FCAL photons, ~0.5ns for BCAL photons | ||
+ | {| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;" | ||
+ | !width="50"| PID | ||
+ | !width="150"| BCAL/RF Δt (ns) | ||
+ | !width="150"| TOF/RF Δt (ns) | ||
+ | !width="150"| FCAL/RF Δt (ns) | ||
+ | !width="150"| SC/RF Δt (ns) | ||
+ | |- | ||
+ | ! γ | ||
+ | | ±1.5 || NA || ±2.5 || NA | ||
+ | |- | ||
+ | ! e<sup>+</sup> | ||
+ | | ±1.0 || ±0.5 || ±2.0 || ±2.5 | ||
+ | |- | ||
+ | ! e<sup>-</sup> | ||
+ | | ±1.0 || ±0.5 || ±2.0 || ±2.5 | ||
+ | |- | ||
+ | ! μ<sup>+</sup> | ||
+ | | ±1.0 || ±0.5 || ±2.0 || ±2.5 | ||
+ | |- | ||
+ | ! μ<sup>-</sup> | ||
+ | | ±1.0 || ±0.5 || ±2.0 || ±2.5 | ||
+ | |- | ||
+ | ! π<sup>+</sup> | ||
+ | | ±1.0 || ±0.5 || ±2.0 || ±2.5 | ||
+ | |- | ||
+ | ! π<sup>-</sup> | ||
+ | | ±1.0 || ±0.5 || ±2.0 || ±2.5 | ||
+ | |- | ||
+ | ! K<sup>+</sup> | ||
+ | | ±0.75 || ±0.3 || ±2.5 || ±2.5 | ||
+ | |- | ||
+ | ! K<sup>-</sup> | ||
+ | | ±0.75 || ±0.3 || ±2.5 || ±2.5 | ||
+ | |- | ||
+ | ! p | ||
+ | | ±1.0 || ±0.6 || ±2.0 || ±2.5 | ||
+ | |- | ||
+ | ! p̅ | ||
+ | | ±1.0 || ±0.6 || ±2.0 || ±2.5 | ||
+ | |} | ||
+ | <br style="clear:both;"/> | ||
− | == Energy Cuts (dE/dx | + | |
− | * | + | '''Note:''' The code where cuts are applied is in [https://github.com/JeffersonLab/halld_recon/blob/master/src/libraries/ANALYSIS/DSourceComboTimeHandler.cc halld_recon/src/libraries/ANALYSIS/DSourceComboTimeHandler.cc] and non-standard cuts can be applied by setting JANA configuration parameters like |
− | * | + | |
− | * | + | COMBO_TIMECUT:14_32=0.75 //Cut protons (14) in the FCAL (32) at +/- 0.75 |
− | + | ||
− | * | + | == Track Energy Loss Cuts == |
− | * | + | |
− | * | + | {| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;" |
+ | !width="50"| PID | ||
+ | !width="300"| CDC dE/dx (keV/cm) | ||
+ | !width="300"| FDC dE/dx (keV/cm) | ||
+ | !width="300"| SC dE/dx (MeV/cm) | ||
+ | !width="300"| TOF dE/dx (MeV/cm) | ||
+ | !width="200"| BCAL Shower_E/Track_p (c) | ||
+ | !width="200"| FCAL Shower_E/Track_p (c) | ||
+ | |- | ||
+ | ! p | ||
+ | | dE/dx > e^(-4.0*p + 2.25) + 1.0 || None || None || None || None || None | ||
+ | |- | ||
+ | ! p̅ | ||
+ | | dE/dx > e^(-4.0*p + 2.25) + 1.0 || None || None || None || None || None | ||
+ | |- | ||
+ | ! π<sup>+</sup> | ||
+ | | dE/dx < e^(-7.0*p + 3.0) + 6.2 || None || None || None || None || None | ||
+ | |- | ||
+ | ! π<sup>-</sup> | ||
+ | | dE/dx < e^(-7.0*p + 3.0) + 6.2 || None || None || None || None || None | ||
+ | |- | ||
+ | ! K<sup>+</sup> | ||
+ | | dE/dx < e^(-7.0*p + 3.0) + 6.2 || None || None || None || None || None | ||
+ | |- | ||
+ | ! K<sup>-</sup> | ||
+ | | dE/dx < e^(-7.0*p + 3.0) + 6.2 || None || None || None || None || None | ||
+ | |- | ||
+ | ! e<sup>+</sup> | ||
+ | | dE/dx < 5.5 || None || None || None || E/p > 0.7 || E/p > 0.7 | ||
+ | |- | ||
+ | ! e<sup>-</sup> | ||
+ | | dE/dx < 5.5 || None || None || None || E/p > 0.7 || E/p > 0.7 | ||
+ | |- | ||
+ | ! μ<sup>+</sup> | ||
+ | | dE/dx < e^(-7.0*p + 3.0) + 6.2 || None || None || None || E/p < 0.45 || E/p < 0.2 | ||
+ | |- | ||
+ | ! μ<sup>-</sup> | ||
+ | | dE/dx < e^(-7.0*p + 3.0) + 6.2 || None || None || None || E/p < 0.45 || E/p < 0.2 | ||
+ | |} | ||
+ | <br style="clear:both;"/> | ||
+ | |||
+ | == Mass Cuts == | ||
+ | * Note that if there are no particles missing, there is also a cut on the missing energy: -3.0 < Missing Energy (GeV) < 3.0 | ||
+ | |||
+ | * Detected Particles | ||
+ | {| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;" | ||
+ | !width="50"| PID | ||
+ | !width="200"| Invariant Mass (GeV/c<sup>2</sup>) | ||
+ | !width="200"| Missing Mass<sup>2</sup> (GeV/c<sup>2</sup>)<sup>2</sup> | ||
+ | |- | ||
+ | ! None | ||
+ | | NA || -0.1 < MM<sup>2</sup> < 0.1 | ||
+ | |- | ||
+ | ! γ | ||
+ | | NA || -0.1 < MM<sup>2</sup> < 0.1 | ||
+ | |- | ||
+ | ! e<sup>+</sup> | ||
+ | | NA || -1.0 < MM<sup>2</sup> < 1.0 | ||
+ | |- | ||
+ | ! e<sup>-</sup> | ||
+ | | NA || -1.0 < MM<sup>2</sup> < 1.0 | ||
+ | |- | ||
+ | ! π<sup>-</sup> | ||
+ | | None || -1.0 < MM<sup>2</sup> < 1.0 | ||
+ | |- | ||
+ | ! π<sup>+</sup> | ||
+ | | None || -1.0 < MM<sup>2</sup> < 1.0 | ||
+ | |- | ||
+ | ! K<sup>+</sup> | ||
+ | | None || -1.0 < MM<sup>2</sup> < 2.0 | ||
+ | |- | ||
+ | ! K<sup>-</sup> | ||
+ | | None || -1.0 < MM<sup>2</sup> < 2.0 | ||
+ | |- | ||
+ | ! p | ||
+ | | NA || -0.5 < MM<sup>2</sup> < 4.41 | ||
+ | |- | ||
+ | ! p̅ | ||
+ | | NA || -0.5 < MM<sup>2</sup> < 4.41 | ||
+ | |- | ||
+ | ! n | ||
+ | | None || -0.5 < MM<sup>2</sup> < 4.41 | ||
+ | |} | ||
+ | <br style="clear:both;"/> | ||
+ | |||
+ | * Decaying Mesons | ||
+ | {| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;" | ||
+ | !width="50"| PID | ||
+ | !width="200"| Invariant Mass (GeV/c<sup>2</sup>) | ||
+ | !width="200"| Missing Mass<sup>2</sup> (GeV/c<sup>2</sup>)<sup>2</sup> | ||
+ | |- | ||
+ | ! π<sup>0</sup> | ||
+ | | 0.08 < IM < 0.19 || -1.0 < MM<sup>2</sup> < 1.0 | ||
+ | |- | ||
+ | ! K<sub>S</sub> | ||
+ | | 0.3 < IM < 0.7 || -1.0 < MM<sup>2</sup> < 2.0 | ||
+ | |- | ||
+ | ! η | ||
+ | | 0.35 < IM < 0.75 || None | ||
+ | |- | ||
+ | ! ω | ||
+ | | 0.4 < IM < 1.2 || None | ||
+ | |- | ||
+ | ! η' | ||
+ | | 0.6 < IM < 1.3 || None | ||
+ | |- | ||
+ | ! φ | ||
+ | | 0.8 < IM < 1.2 || None | ||
+ | |- | ||
+ | ! D<sup>0</sup> | ||
+ | | 1.8 < IM < 1.92 || None | ||
+ | |- | ||
+ | ! J/ψ | ||
+ | | None || None | ||
+ | |} | ||
+ | <br style="clear:both;"/> | ||
+ | |||
+ | * Decaying Baryons | ||
+ | {| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;" | ||
+ | !width="50"| PID | ||
+ | !width="200"| Invariant Mass (GeV/c<sup>2</sup>) | ||
+ | !width="200"| Missing Mass<sup>2</sup> (GeV/c<sup>2</sup>)<sup>2</sup> | ||
+ | |- | ||
+ | ! Λ | ||
+ | | 1.0 < IM < 1.2 || None | ||
+ | |- | ||
+ | ! Λ̅ | ||
+ | | 1.0 < IM < 1.2 || None | ||
+ | |- | ||
+ | ! Σ<sup>0</sup> | ||
+ | | 1.1 < IM < 1.3 || None | ||
+ | |- | ||
+ | ! Σ<sup>+</sup> | ||
+ | | 1.1 < IM < 1.3 || None | ||
+ | |- | ||
+ | ! Σ<sup>-</sup> | ||
+ | | 1.1 < IM < 1.3 || None | ||
+ | |- | ||
+ | ! Ξ<sup>0</sup> | ||
+ | | 1.1 < IM < 1.5 || None | ||
+ | |- | ||
+ | ! Ξ<sup>-</sup> | ||
+ | | 1.1 < IM < 1.5 || None | ||
+ | |- | ||
+ | ! Ω<sup>-</sup> | ||
+ | | 1.32 < IM < 2.22 || 1.44 < MM<sup>2</sup> < 4.84 | ||
+ | |- | ||
+ | ! Λ<sub>c</sub> | ||
+ | | 2.0 < IM < 2.6 || None | ||
+ | |} | ||
+ | <br style="clear:both;"/> | ||
+ | |||
+ | == Kinematic Fit == | ||
+ | * If a kinematic fit has been performed, it must have converged. | ||
+ | |||
+ | == Fiducial Cuts == | ||
+ | * Implemented in August 2019, first used in Analysis Launch 2017-01 ver31, 2018-08 ver08, and 2018-08 ver02 | ||
+ | |||
+ | * Maximum Z in BCal: 393.0cm | ||
+ | * Maximum R in FCal: 105.5cm | ||
+ | * Block innermost 16 FCal channels around the beam hole | ||
+ | |||
+ | These cuts are hardcoded in [https://github.com/JeffersonLab/halld_recon/blob/master/src/libraries/PID/DNeutralShower_factory_PreSelect.cc halld_recon/src/libraries/PID/DNeutralShower_factory_PreSelect.cc]. | ||
+ | |||
+ | == Analysis Launches Using These Cuts== | ||
+ | * [https://halldweb.jlab.org/wiki-private/index.php/Spring_2016_Analysis_Launch Spring 2016]: ver09 -> ver17 (as of Aug.30, 2018) | ||
+ | * [https://halldweb.jlab.org/wiki-private/index.php/Spring_2017_Analysis_Launch Spring 2017]: ver05 -> ver16 (as of Aug.30, 2018) |
Latest revision as of 14:58, 7 June 2022
These cuts also apply to 2018 and 2019 GlueX data
Contents
PreSelect Cuts
- For a given DChargedTrack, at least one hypothesis must have a hit in either: BCAL/FCAL/TOF/SC.
- For a given DNeutralShower, it must at least have a shower energy of 100 MeV/c.
- For neutral showers in the BCAL, they must have hits in at least 2 cells.
Timing Cuts
- Timing cuts are applied only on the system with the best timing information available, in the order: BCAL/TOF/FCAL/SC. In other words, if a track has a hit in the TOF, then no cut is applied on the FCAL or SC time.
- Note that some hypotheses may not have a hit in any system and thus survive all of these cuts.
- Also, even if the SC is the best-available detector, no cut is placed on the timing if the track has been matched to more than one SC hit, UNLESS it's the only track in the combo.
- This is because the SC timing is not really good enough for PID over such a short flight distance, so the cut is primarily used to reject tracks from other beam bunches or ghost tracks.
- Also, all kaons are cut that:
- Do not have a hit in one of these timing detectors, AND
- Do not have enough hits in the CDC to compute dE/dx
- This is because their background is huge already, without letting these combos through.
- Cuts are re-placed after the kinematic fit, using updated momentum and vertex information.
- If a kinfit is not performed, then the timing cut at detached vertices is looser than otherwise: +/- 2ns WIDER for charged tracks, ~0.05ns for FCAL photons, ~0.5ns for BCAL photons
PID | BCAL/RF Δt (ns) | TOF/RF Δt (ns) | FCAL/RF Δt (ns) | SC/RF Δt (ns) |
---|---|---|---|---|
γ | ±1.5 | NA | ±2.5 | NA |
e+ | ±1.0 | ±0.5 | ±2.0 | ±2.5 |
e- | ±1.0 | ±0.5 | ±2.0 | ±2.5 |
μ+ | ±1.0 | ±0.5 | ±2.0 | ±2.5 |
μ- | ±1.0 | ±0.5 | ±2.0 | ±2.5 |
π+ | ±1.0 | ±0.5 | ±2.0 | ±2.5 |
π- | ±1.0 | ±0.5 | ±2.0 | ±2.5 |
K+ | ±0.75 | ±0.3 | ±2.5 | ±2.5 |
K- | ±0.75 | ±0.3 | ±2.5 | ±2.5 |
p | ±1.0 | ±0.6 | ±2.0 | ±2.5 |
p̅ | ±1.0 | ±0.6 | ±2.0 | ±2.5 |
Note: The code where cuts are applied is in halld_recon/src/libraries/ANALYSIS/DSourceComboTimeHandler.cc and non-standard cuts can be applied by setting JANA configuration parameters like
COMBO_TIMECUT:14_32=0.75 //Cut protons (14) in the FCAL (32) at +/- 0.75
Track Energy Loss Cuts
PID | CDC dE/dx (keV/cm) | FDC dE/dx (keV/cm) | SC dE/dx (MeV/cm) | TOF dE/dx (MeV/cm) | BCAL Shower_E/Track_p (c) | FCAL Shower_E/Track_p (c) |
---|---|---|---|---|---|---|
p | dE/dx > e^(-4.0*p + 2.25) + 1.0 | None | None | None | None | None |
p̅ | dE/dx > e^(-4.0*p + 2.25) + 1.0 | None | None | None | None | None |
π+ | dE/dx < e^(-7.0*p + 3.0) + 6.2 | None | None | None | None | None |
π- | dE/dx < e^(-7.0*p + 3.0) + 6.2 | None | None | None | None | None |
K+ | dE/dx < e^(-7.0*p + 3.0) + 6.2 | None | None | None | None | None |
K- | dE/dx < e^(-7.0*p + 3.0) + 6.2 | None | None | None | None | None |
e+ | dE/dx < 5.5 | None | None | None | E/p > 0.7 | E/p > 0.7 |
e- | dE/dx < 5.5 | None | None | None | E/p > 0.7 | E/p > 0.7 |
μ+ | dE/dx < e^(-7.0*p + 3.0) + 6.2 | None | None | None | E/p < 0.45 | E/p < 0.2 |
μ- | dE/dx < e^(-7.0*p + 3.0) + 6.2 | None | None | None | E/p < 0.45 | E/p < 0.2 |
Mass Cuts
- Note that if there are no particles missing, there is also a cut on the missing energy: -3.0 < Missing Energy (GeV) < 3.0
- Detected Particles
PID | Invariant Mass (GeV/c2) | Missing Mass2 (GeV/c2)2 |
---|---|---|
None | NA | -0.1 < MM2 < 0.1 |
γ | NA | -0.1 < MM2 < 0.1 |
e+ | NA | -1.0 < MM2 < 1.0 |
e- | NA | -1.0 < MM2 < 1.0 |
π- | None | -1.0 < MM2 < 1.0 |
π+ | None | -1.0 < MM2 < 1.0 |
K+ | None | -1.0 < MM2 < 2.0 |
K- | None | -1.0 < MM2 < 2.0 |
p | NA | -0.5 < MM2 < 4.41 |
p̅ | NA | -0.5 < MM2 < 4.41 |
n | None | -0.5 < MM2 < 4.41 |
- Decaying Mesons
PID | Invariant Mass (GeV/c2) | Missing Mass2 (GeV/c2)2 |
---|---|---|
π0 | 0.08 < IM < 0.19 | -1.0 < MM2 < 1.0 |
KS | 0.3 < IM < 0.7 | -1.0 < MM2 < 2.0 |
η | 0.35 < IM < 0.75 | None |
ω | 0.4 < IM < 1.2 | None |
η' | 0.6 < IM < 1.3 | None |
φ | 0.8 < IM < 1.2 | None |
D0 | 1.8 < IM < 1.92 | None |
J/ψ | None | None |
- Decaying Baryons
PID | Invariant Mass (GeV/c2) | Missing Mass2 (GeV/c2)2 |
---|---|---|
Λ | 1.0 < IM < 1.2 | None |
Λ̅ | 1.0 < IM < 1.2 | None |
Σ0 | 1.1 < IM < 1.3 | None |
Σ+ | 1.1 < IM < 1.3 | None |
Σ- | 1.1 < IM < 1.3 | None |
Ξ0 | 1.1 < IM < 1.5 | None |
Ξ- | 1.1 < IM < 1.5 | None |
Ω- | 1.32 < IM < 2.22 | 1.44 < MM2 < 4.84 |
Λc | 2.0 < IM < 2.6 | None |
Kinematic Fit
- If a kinematic fit has been performed, it must have converged.
Fiducial Cuts
- Implemented in August 2019, first used in Analysis Launch 2017-01 ver31, 2018-08 ver08, and 2018-08 ver02
- Maximum Z in BCal: 393.0cm
- Maximum R in FCal: 105.5cm
- Block innermost 16 FCal channels around the beam hole
These cuts are hardcoded in halld_recon/src/libraries/PID/DNeutralShower_factory_PreSelect.cc.
Analysis Launches Using These Cuts
- Spring 2016: ver09 -> ver17 (as of Aug.30, 2018)
- Spring 2017: ver05 -> ver16 (as of Aug.30, 2018)