Difference between revisions of "Mattione Update 09042013"

From GlueXWiki
Jump to: navigation, search
(→‎Single π+ Track)
(→‎CDC Track Finding - Basics)
 
(89 intermediate revisions by the same user not shown)
Line 1: Line 1:
−
== CDC Track Finding ==
+
= CDC Track Finding - Basics =
 +
* Note: Some steps (and many details) are omitted entirely in this brief overview. 
  
−
== Single &pi;<sup>+</sup> Track ==
+
=== Build Super Layer Seeds ===
 +
* In each ring, link adjacent hits together (<span style="color:#0000FF">DCDCRingSeed</span>).
 +
* In each super layer, link each possible combination of adjacent <span style="color:#0000FF">DCDCRingSeed</span>'s together to form super layer seeds (<span style="color:#0000FF">DCDCSuperLayerSeed</span>).
 +
** <span style="color:red">NEW</span>: Do for both axial super layers and stereo super layers (previously just axial).
 +
* <span style="color:red">NEW</span>: Identify which <span style="color:#0000FF">DCDCSuperLayerSeed</span>'s contain potential spiral turns (where a spiraling track turns back inward/outward).
 +
{| class="wikitable"
 +
| [[Image:Mattione_Update_09042013_HDView1.png|thumb|middle|200px|b1pi Event]]
 +
| [[Image:Mattione_Update_08202013_SlowPiPlus_Event17_ReconHits.png|thumb|right|200px|Spiral Event]]
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== Build Track Circles ===
 +
* <span style="color:red">HEAVILY MODIFIED</span>: Link adjacent <span style="color:#0000FF">DCDCSuperLayerSeed</span>'s together to form tracks (<span style="color:#0000FF">DCDCTrackCircle</span>).
 +
** Previously, the axial <span style="color:#0000FF">DCDCSuperLayerSeed</span>'s were linked by progressively fitting track circles and selecting the next axial super layers that were close to the track fit.
 +
** Start in super layer 1 (axial) and expand outwards towards super layer 7 (axial) (was previously in reverse direction)
 +
** For each unique combination of axial super layers, create a new <span style="color:#0000FF">DCDCTrackCircle</span> object.
 +
** Keep track of which stereo <span style="color:#0000FF">DCDCSuperLayerSeed</span> were used when linking the <span style="color:#0000FF">DCDCTrackCircle</span> objects.  There may be (many) more than one for each super layer.
 +
 
 +
=== Fit Track Circles & Truncate: Axial ===
 +
[[Image:Mattione_Update_09042013_HDView2.png|thumb|right|200px|b1pi Event]]
 +
* Fit the hits in the axial super layers (plus the beamline) to a circle with <span style="color:#0000FF">DHelicalFit</span>.
 +
** Reject tracks for which half of the axial hits are not close to the circle fit.
 +
* <span style="color:red">NEW</span>: Truncate and re-fit track circles that have picked up erroneous super layer seeds in the outer super layers. 
 +
** This is used to save (for example) forward-going tracks that pick up bad hits from another track (e.g. a large-angle track near in &phi;)
 +
<br style="clear:both;"/>
 +
 
 +
=== Filter Track Circles: Axial ===
 +
[[Image:Mattione_Update_08202013_SlowPiPlus_Event17_Spiral.png|right|thumb|200px|Spiral Track]]
 +
* If two track circles share >= 50% of their hits, reject the <span style="color:#0000FF">DCDCTrackCircle</span> with the worse circle-fit weighted-chisq/ndf (weighted by #axial super layers squared)
 +
* <span style="color:red">NEW</span>: If two track circles are spiral-linked in their last axial super layer, have opposite charges, and have similar circle centers, reject the <span style="color:#0000FF">DCDCTrackCircle</span> with the worse circle-fit weighted-chisq/ndf.
 +
<br style="clear:both;"/>
 +
 
 +
=== Select Stereo Super Layer Seeds ===
 +
[[Image:Mattione_Update_09042013_HDView3.png|thumb|right|200px|b1pi Event]]
 +
* Previously, all stereo hits that intersected the circle and were nearby were lumped all together. Now, stereo super layer seeds are used instead.
 +
* Project all hits from all possible stereo super layers onto the circle.
 +
* For all possible combinations of stereo super layers:
 +
** Do a linear regression fit of the projected stereo hits to calculate &theta; & Vertex-Z
 +
** <span style="color:red">NEW</span>: Select the combination of stereo super layer seeds that have the best chisq/ndf from the regression fit.
 +
<br style="clear:both;"/>
 +
 
 +
=== Filter Track Circles: Stereo ===
 +
* <span style="color:red">NEW</span>: Find super layer seeds that are shared between <span style="color:#0000FF">DCDCTrackCircle</span>'s, and remove them from all <span style="color:#0000FF">DCDCTrackCircle</span>'s except the one with the best weighted chisq/ndf from the regression fit.
 +
** Like the axial truncation, this is used to save (for example) forward-going tracks that pick up bad hits from another track (e.g. a large-angle track near in &phi;).
 +
** However, do NOT remove the stereo super layer if it's the last one on the track. 
 +
* <span style="color:red">NEW</span>: If two track circles are spiral-linked in their last stereo super layer, have opposite charges, and have similar circle centers, reject the <span style="color:#0000FF">DCDCTrackCircle</span> with the worse circle-fit weighted-chisq/ndf.
 +
 
 +
=== Refinement ===
 +
[[Image:Mattione_Update_09042013_HDView4.png|thumb|right|200px|b1pi Event]]
 +
* Find single, unused hits in the next (outward) super layer and add them to track circles if nearby.
 +
** This primarily helps tracks that exit the CDC early (e.g. forward going). 
 +
* Redo circle fits, using both the new hits added just before, as well as the intersection points between the two stereo super layers (e.g. SL2 & SL3).
 +
* Re-project stereo hits onto the new circle fits.
 +
* Do a linear regression fit on a subset of the projected stereo hits to calculate &theta; & Vertex-Z
 +
** <span style="color:red">SEMI-NEW</span>: A subset is chosen to optimize the final fit results.  This is because there are regions/cases for which the projected stereo-hit-positions are bad.
 +
<br style="clear:both;"/>
 +
 
 +
=== Finalize ===
 +
[[Image:Mattione_Update_09042013_SpiralDeltaPhi.png|thumb|right|300px|Spiral Events (2nd peak due to wrong charge)]]
 +
* Calculate the position and momentum of the track at a point outside of the target and the start counter.
 +
** This is because the track finding does not take energy-loss into account.
 +
** Note that this skews the comparison between the <span style="color:#0000FF">DTrackCandidate</span> and <span style="color:#0000FF">DMCThrown</span> kinematics, since they are reported at different positions (r ~= 9.5cm, r = 0cm, respectively).
 +
* Create <span style="color:#0000FF">DTrackCandidate</span> objects.
 +
<br style="clear:both;"/>
 +
 
 +
<!--
 +
= CDC Track Finding - Details =
 +
=== Build Super Layer Seeds ===
 +
* In each ring, link adjacent hits together (DCDCRingSeed).
 +
** NEW: Do for both axial super layers and stereo super layers (previously just axial).
 +
* In each super layer, link each possible combination of adjacent DCDCRingSeed's together to form super layer seeds (DCDCSuperLayerSeed).
 +
** Note: There may be many combinations. For example, if two tracks cross in an "X" pattern, this will result in four DCDCSuperLayerSeed's.
 +
** NEW: Do for both axial super layers and stereo super layers (previously just axial).
 +
** NEW: Can skip up to 1 ring (in case the track did not deposit enough energy in one of the rings to register a hit). 
 +
** Require a minimum of two hits to form a DCDCSuperLayerSeed (unmatched hits will be added later).
 +
* NEW: In each super layer, identify regions with high hit density (> 60% in an area of ~24 straws) and reject all DCDCSuperLayerSeed's passing through those regions.
 +
** This occurs when a spiraling track loops around 6+ times before stopping.
 +
** This is necessary because in this new code, these events effectively freeze the program with millions of possible combinations of super layer seeds.
 +
* NEW: Identify which DCDCSuperLayerSeed's contain potential spiral turns (where a spiraling track turns back inward/outward)
 +
** Primary signature is many adjacent hits in the same ring.
 +
 
 +
=== Link Super Layer Seeds ===
 +
* HEAVILY MODIFIED: Link adjacent DCDCSuperLayerSeed's together to form tracks (DCDCTrackCircle).
 +
** Previously, the axial DCDCSuperLayerSeed's were linked by progressively fitting track circles and selecting the next axial super layers that were close to the track fit.
 +
** Note: The projection of the stereo straws onto the X-Y (axial) plane are lines. Therefore, adjacent super layers are determined by checking whether a point (axial) is near a line (stereo), or a line (stereo) is near a line (stereo). 
 +
** While the new method results in many more super layer combinations to disentangle, every possible combination is guaranteed to be considered.
 +
** NEW: There is an option to allow skipping of (at most) one super layer when linking: useful for when a high-voltage board is dead. 
 +
* Start in super layer 1 (axial) and expand outwards towards super layer 7 (axial) (was previously in reverse direction)
 +
* For each unique combination of axial super layers, create a new DCDCTrackCircle object.
 +
** Do not attempt seed link if the outer seed is definitely an outward-turning spiral, or if the inner seed is definitely an inward-turning spiral.
 +
** Keep track of which stereo DCDCSuperLayerSeed were used when linking the DCDCTrackCircle objects.  There may be (many) more than one for each super layer. 
 +
 
 +
=== Fit, Truncate, Filter Track Circles: Axial ===
 +
* Fit the hits in the axial super layers plus the beamline to a circle with DHelicalFit.
 +
* Truncate any re-fit track circles that have picked up erroneous super layer seeds in the outer super layers. 
 +
** For example, forward-going tracks that are near in &phi; to large-angle tracks: these track circles will often pick up hits that were on the large-angle track, causing that track to be rejected during filtering below. 
 +
*
 +
 
 +
=== Filter Track Circles: Axial ===
 +
* NEW: When filtering track circles, loop over them by comparing the worst circle-fits to the best ones (the order is key to reducing multiplicity).
 +
 
 +
=== Select Stereo Super Layer Seeds ===
 +
* Project all hits from all possible stereo super layers onto the circle
 +
** Ignore, but keep, hits that don't intersect the circle.
 +
 
 +
-->
 +
 
 +
= Single &pi;<sup>+</sup> Track =
 
=== Single &pi;<sup>+</sup> Reconstruction: Current (Trunk) Code, trackeff_hists plugin ===
 
=== Single &pi;<sup>+</sup> Reconstruction: Current (Trunk) Code, trackeff_hists plugin ===
 
* The below plots show the track reconstruction efficiency and resolution for the trunk using the single-track scripts plugin (trackeff_hists).   
 
* The below plots show the track reconstruction efficiency and resolution for the trunk using the single-track scripts plugin (trackeff_hists).   
Line 16: Line 124:
 
* This plugin compares matches all of the thrown and reconstructed particles in the event, and requires that their momentum be "close."
 
* This plugin compares matches all of the thrown and reconstructed particles in the event, and requires that their momentum be "close."
 
** "Close" for track candidates is defined as: &Delta;p/p < 20%, &Delta;&theta; < 40 degrees, &Delta;&phi; < 40 degrees, &Delta;Vertex-Z < 1 km
 
** "Close" for track candidates is defined as: &Delta;p/p < 20%, &Delta;&theta; < 40 degrees, &Delta;&phi; < 40 degrees, &Delta;Vertex-Z < 1 km
−
*** For thrown tracks with &theta; < 5 degrees, there is no cut on &Delta;&phi;
+
*** For thrown tracks with &theta; < 5 degrees, there is no cut on the candidate &Delta;&phi;
 
** "Close" for time-based tracks is defined as: &Delta;p/p < 10%, &Delta;&theta; < 15 degrees, &Delta;&phi; < 15 degrees, &Delta;Vertex-Z < 10 cm
 
** "Close" for time-based tracks is defined as: &Delta;p/p < 10%, &Delta;&theta; < 15 degrees, &Delta;&phi; < 15 degrees, &Delta;Vertex-Z < 10 cm
 
{| class="wikitable"
 
{| class="wikitable"
Line 29: Line 137:
 
* Note: The trackeffv2 plugin compares matches all of the thrown and reconstructed particles in the event, and requires that their momentum be "close."
 
* Note: The trackeffv2 plugin compares matches all of the thrown and reconstructed particles in the event, and requires that their momentum be "close."
 
** "Close" for track candidates is defined as: &Delta;p/p < 20%, &Delta;&theta; < 40 degrees, &Delta;&phi; < 40 degrees, &Delta;Vertex-Z < 1 km
 
** "Close" for track candidates is defined as: &Delta;p/p < 20%, &Delta;&theta; < 40 degrees, &Delta;&phi; < 40 degrees, &Delta;Vertex-Z < 1 km
−
*** For thrown tracks with &theta; < 5 degrees, there is no cut on &Delta;&phi;
+
*** For thrown tracks with &theta; < 5 degrees, there is no cut on the candidate &Delta;&phi;
 
** "Close" for time-based tracks is defined as: &Delta;p/p < 10%, &Delta;&theta; < 15 degrees, &Delta;&phi; < 15 degrees, &Delta;Vertex-Z < 10 cm
 
** "Close" for time-based tracks is defined as: &Delta;p/p < 10%, &Delta;&theta; < 15 degrees, &Delta;&phi; < 15 degrees, &Delta;Vertex-Z < 10 cm
 
{| class="wikitable"
 
{| class="wikitable"
 
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_Single_Spiral.png|thumb|left|500px|Spiral Code: Track Candidates Efficiency (Spiral Code)]]
 
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_Single_Spiral.png|thumb|left|500px|Spiral Code: Track Candidates Efficiency (Spiral Code)]]
 
| [[Image:Mattione_Update_09042013_Efficiency_TBT_Single_Spiral.png|thumb|left|500px|Spiral Code: Time-Based Track Efficiency (Spiral Code)]]
 
| [[Image:Mattione_Update_09042013_Efficiency_TBT_Single_Spiral.png|thumb|left|500px|Spiral Code: Time-Based Track Efficiency (Spiral Code)]]
 +
|}
 +
 +
* Selected Track Candidate 1-D Kinematics Comparisons
 +
{| class="wikitable"
 +
| [[Image:Mattione_Update_09042013_KinematicsV1_Candidates_Single.png|thumb|left|375px|Track Candidates (Forward Angle)]]
 +
| [[Image:Mattione_Update_09042013_KinematicsV2_Candidates_Single.png|thumb|left|375px|Track Candidates (Transition Region)]]
 +
| [[Image:Mattione_Update_09042013_KinematicsV3_Candidates_Single.png|thumb|left|375px|Track Candidates (Spirals)]]
 
|}
 
|}
  
Line 45: Line 160:
 
|}
 
|}
  
 +
<!--
 
== &gamma;p&rarr;&pi;<sup>+</sup>&pi;<sup>-</sup>p (&rho; Mesons) ==
 
== &gamma;p&rarr;&pi;<sup>+</sup>&pi;<sup>-</sup>p (&rho; Mesons) ==
 
* The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
 
* The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
Line 50: Line 166:
 
* Events were generated at beam energy of 9 GeV and a t-slope of 5.0
 
* Events were generated at beam energy of 9 GeV and a t-slope of 5.0
  
−
=== Protons ===
+
=== Proton ===
−
[[Image:Mattione_Update_09042013_Thrown_Proton_Rho.png|left|350px|Thrown Protons]]
+
[[Image:Mattione_Update_09042013_Thrown_Proton_Rho.png|thumb|left|350px|Thrown Protons]]
 
<br style="clear:both;"/>
 
<br style="clear:both;"/>
 
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
Line 71: Line 187:
  
 
=== &pi;<sup>+</sup> ===
 
=== &pi;<sup>+</sup> ===
−
[[Image:Mattione_Update_09042013_Thrown_PiPlus_Rho.png|left|350px|Thrown &pi;<sup>+</sup>'s]]
+
[[Image:Mattione_Update_09042013_Thrown_PiPlus_Rho.png|thumb|left|350px|Thrown &pi;<sup>+</sup>'s]]
 
<br style="clear:both;"/>
 
<br style="clear:both;"/>
 
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
Line 91: Line 207:
  
 
=== &pi;<sup>-</sup> ===
 
=== &pi;<sup>-</sup> ===
−
[[Image:Mattione_Update_09042013_Thrown_PiMinus_Rho.png|left|350px|Thrown &pi;<sup>-</sup>'s]]
+
[[Image:Mattione_Update_09042013_Thrown_PiMinus_Rho.png|thumb|left|350px|Thrown &pi;<sup>-</sup>'s]]
 
<br style="clear:both;"/>
 
<br style="clear:both;"/>
 
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
Line 109: Line 225:
 
|}
 
|}
 
<br style="clear:both;"/>
 
<br style="clear:both;"/>
 +
-->
  
−
== &gamma;p&rarr;&pi;<sup>+</sup>&pi;<sup>+</sup>&pi;<sup>-</sup>(n) ==
+
= &gamma;p&rarr;&pi;<sup>+</sup>&pi;<sup>+</sup>&pi;<sup>-</sup>(n) =
 +
* The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
 +
* The definition of closeness is the same as before (see the single-track section).
 +
* Events were generated at beam energy of 9 GeV and a t-slope of 5.0
  
−
== &gamma;p&rarr;&pi;<sup>+</sup>&pi;<sup>+</sup>&pi;<sup>-</sup>&pi;<sup>-</sup>&pi;<sup>0</sup>(p) (b1pi) ==
+
=== Summary ===
 +
* Overall track reconstruction improvement percentages for all tracks in this topology with momentum LESS than 3 GeV/c.
 +
** WARNING: This is not the overall gain for the topology, as many pions have larger than 3 GeV/c of momentum (their overall gain will be smaller). 
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Particle
 +
!width="150"| Track Candidates
 +
!width="150"| Time-Based Tracks
 +
|-
 +
! &pi;<sup>+</sup>
 +
| 10.3% || 5.38%
 +
|-
 +
! &pi;<sup>-</sup>
 +
| 8.08% || 4.47%
 +
|}
 +
<br style="clear:both;"/>
  
−
== bggen (pythia) ==
+
=== &pi;<sup>+</sup> ===
 +
* Thrown, Track Candidate 1-D Kinematics
 +
[[Image:Mattione_Update_09042013_Thrown_PiPlus_n3pi.png|thumb|left|350px|Thrown &pi;<sup>+</sup>'s]]
 +
[[Image:Mattione_Update_09042013_KinematicsV1_PiPlus_n3pi.png|thumb|middle|275px|&pi;<sup>+</sup> Track Candidate Kinematics]]
 +
<br style="clear:both;"/>
 +
* Efficiencies
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Stage
 +
!width="250"| Efficiency - Current Code
 +
!width="250"| Efficiency - Spiral Code
 +
!width="250"| Improvement
 +
!width="250"| Improvement - Zoomed
 +
|-
 +
! Track Candidates
 +
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_n3pi_PiPlus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_Candidates_n3pi_PiPlus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_Candidates_n3pi_PiPlus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_Candidates_n3pi_PiPlus.png|left|275px]]
 +
|-
 +
! Time-Based Tracks
 +
| [[Image:Mattione_Update_09042013_Efficiency_TimeBased_n3pi_PiPlus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_TimeBased_n3pi_PiPlus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_TimeBased_n3pi_PiPlus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_TimeBased_n3pi_PiPlus.png|left|275px]]
 +
|}
 +
<br style="clear:both;"/>
  
−
== &gamma;p&rarr;K<sup>+</sup>K<sup>+</sup>&Xi;<sup>-</sup> ==
+
= &gamma;p&rarr;&pi;<sup>+</sup>&pi;<sup>+</sup>&pi;<sup>-</sup>&pi;<sup>-</sup>&pi;<sup>0</sup>(p) (b1pi) =
 +
* The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
 +
* The definition of closeness is the same as before (see the single-track section).
 +
* Events were generated with a beam energy of 9 GeV and a t-slope of 5.0
 +
 
 +
=== Summary ===
 +
* Overall track reconstruction improvement percentages for all tracks in this topology with momentum LESS than 3 GeV/c.
 +
** WARNING: This is not the overall gain for the topology, as many pions have larger than 3 GeV/c of momentum (their overall gain will be smaller).
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Particle
 +
!width="150"| Track Candidates
 +
!width="150"| Time-Based Tracks
 +
|-
 +
! &pi;<sup>+</sup>
 +
| 7.93% || 5.21%
 +
|-
 +
! &pi;<sup>-</sup>
 +
| 12.5% || 6.23%
 +
|-
 +
! p
 +
| 23.7% || 22.6%
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== &pi;<sup>-</sup> ===
 +
* Thrown
 +
[[Image:Mattione_Update_09042013_Thrown_PiMinus_b1pi.png|thumb|left|275px|Thrown &pi;<sup>-</sup>'s]]
 +
<br style="clear:both;"/>
 +
* Efficiencies
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Stage
 +
!width="250"| Efficiency - Current Code
 +
!width="250"| Efficiency - Spiral Code
 +
!width="250"| Improvement
 +
!width="250"| Improvement - Zoomed
 +
|-
 +
! Track Candidates
 +
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_b1pi_PiMinus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_Candidates_b1pi_PiMinus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_Candidates_b1pi_PiMinus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_Candidates_b1pi_PiMinus.png|left|275px]]
 +
|-
 +
! Time-Based Tracks
 +
| [[Image:Mattione_Update_09042013_Efficiency_TimeBased_b1pi_PiMinus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_TimeBased_b1pi_PiMinus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_TimeBased_b1pi_PiMinus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_TimeBased_b1pi_PiMinus.png|left|275px]]
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== Proton ===
 +
* Thrown, Track Candidate 1-D Kinematics
 +
[[Image:Mattione_Update_09042013_Thrown_Proton_b1pi.png|thumb|left|275px|Thrown Proton's]]
 +
[[Image:Mattione_Update_09042013_KinematicsV1_Proton_b1pi.png|thumb|middle|275px|Track Candidates]]
 +
<br style="clear:both;"/>
 +
* Efficiencies
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Stage
 +
!width="250"| Efficiency - Current Code
 +
!width="250"| Efficiency - Spiral Code
 +
!width="250"| Improvement
 +
|-
 +
! Track Candidates
 +
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_b1pi_Proton_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_Candidates_b1pi_Proton_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_Candidates_b1pi_Proton.png|left|275px]]
 +
|-
 +
! Time-Based Tracks
 +
| [[Image:Mattione_Update_09042013_Efficiency_TimeBased_b1pi_Proton_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_TimeBased_b1pi_Proton_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_TimeBased_b1pi_Proton.png|left|275px]]
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
= &gamma;p&rarr;K<sup>+</sup>K<sup>+</sup>&Xi;<sup>-</sup> =
 +
* The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
 +
* The definition of closeness is the same as before (see the single-track section).
 +
* Events were generated with beam energies between 8.4 and 9 GeV
 +
 
 +
=== Summary ===
 +
* Overall track reconstruction improvement percentages for all tracks in this topology with momentum LESS than 3 GeV/c.
 +
** WARNING: This is not the overall gain for the topology, as some kaons have larger than 3 GeV/c of momentum (their overall gain will be slightly smaller).
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Particle
 +
!width="150"| Track Candidates
 +
!width="150"| Time-Based Tracks
 +
|-
 +
! K<sup>+</sup>
 +
| 21.7% || 14.1%
 +
|-
 +
! &pi;<sup>-</sup>
 +
| 35.1% || 28.7%
 +
|-
 +
! p
 +
| 12.8% || 16.7%
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== Sample Event ===
 +
 
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="200"| Reconstructed Hits
 +
!width="200"| Truth + Reconstructed Hits
 +
!width="200"| Current (09/03/13) Code Track Candidates
 +
!width="200"| Spiral Code Track Candidates
 +
|-
 +
| [[Image:Mattione_Update_09042013_Cascade_Event3_ReconHits.png|left|250px]] || [[Image:Mattione_Update_09042013_Cascade_Event3_TruthHits.png|left|250px]] || [[Image:Mattione_Update_09042013_Cascade_Event3_Current.png|left|250px]] || [[Image:Mattione_Update_09042013_Cascade_Event3_Spiral.png|left|250px]] ||
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== &pi;<sup>-</sup> ===
 +
* Thrown
 +
[[Image:Mattione_Update_09042013_Thrown_PiMinus_cascade.png|thumb|left|275px|Thrown &pi;<sup>-</sup>'s]]
 +
<br style="clear:both;"/>
 +
* Efficiencies
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Stage
 +
!width="250"| Efficiency - Current Code
 +
!width="250"| Efficiency - Spiral Code
 +
!width="250"| Improvement
 +
!width="250"| Improvement - Zoomed
 +
|-
 +
! Track Candidates
 +
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_cascade_PiMinus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_Candidates_cascade_PiMinus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_Candidates_cascade_PiMinus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_Candidates_cascade_PiMinus.png|left|275px]]
 +
|-
 +
! Time-Based Tracks
 +
| [[Image:Mattione_Update_09042013_Efficiency_TimeBased_cascade_PiMinus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_TimeBased_cascade_PiMinus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_TimeBased_cascade_PiMinus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_TimeBased_cascade_PiMinus.png|left|275px]]
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== K<sup>+</sup> ===
 +
* Thrown
 +
[[Image:Mattione_Update_09042013_Thrown_KPlus_cascade.png|thumb|left|275px|Thrown K<sup>+</sup>'s]]
 +
<br style="clear:both;"/>
 +
* Efficiencies
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Stage
 +
!width="250"| Efficiency - Current Code
 +
!width="250"| Efficiency - Spiral Code
 +
!width="250"| Improvement
 +
!width="250"| Improvement - Zoomed
 +
|-
 +
! Track Candidates
 +
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_cascade_KPlus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_Candidates_cascade_KPlus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_Candidates_cascade_KPlus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_Candidates_cascade_KPlus.png|left|275px]]
 +
|-
 +
! Time-Based Tracks
 +
| [[Image:Mattione_Update_09042013_Efficiency_TimeBased_cascade_KPlus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_TimeBased_cascade_KPlus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_TimeBased_cascade_KPlus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_TimeBased_cascade_KPlus.png|left|275px]]
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== Proton ===
 +
* Thrown
 +
[[Image:Mattione_Update_09042013_Thrown_Proton_cascade.png|thumb|left|275px|Thrown Proton's]]
 +
<br style="clear:both;"/>
 +
* Efficiencies
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Stage
 +
!width="250"| Efficiency - Current Code
 +
!width="250"| Efficiency - Spiral Code
 +
!width="250"| Improvement
 +
|-
 +
! Track Candidates
 +
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_cascade_Proton_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_Candidates_cascade_Proton_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_Candidates_cascade_Proton.png|left|275px]]
 +
|-
 +
! Time-Based Tracks
 +
| [[Image:Mattione_Update_09042013_Efficiency_TimeBased_cascade_Proton_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_TimeBased_cascade_Proton_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_TimeBased_cascade_Proton.png|left|275px]]
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
= bggen (pythia) =
 +
* The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
 +
* The definition of closeness is the same as before (see the single-track section).
 +
* Events were generated with beam energies between 8.4 and 9 GeV
 +
 
 +
=== Summary ===
 +
* Overall track reconstruction improvement percentages for all tracks in this topology with momentum LESS than 3 GeV/c.
 +
** WARNING: This is not the overall gain for the topology, as many pions have larger than 3 GeV/c of momentum (their overall gain will be smaller).
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Particle
 +
!width="150"| Track Candidates
 +
!width="150"| Time-Based Tracks
 +
|-
 +
! &pi;<sup>+</sup>
 +
| 16.7% || 11.3%
 +
|-
 +
! &pi;<sup>-</sup>
 +
| 21.8% || 12.7%
 +
|-
 +
! p
 +
| 12.4% || 6.92%
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== &pi;<sup>-</sup> ===
 +
* Thrown
 +
[[Image:Mattione_Update_09042013_Thrown_PiMinus_bggen.png|thumb|left|275px|Thrown &pi;<sup>-</sup>'s]]
 +
<br style="clear:both;"/>
 +
* Efficiencies
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Stage
 +
!width="250"| Efficiency - Current Code
 +
!width="250"| Efficiency - Spiral Code
 +
!width="250"| Improvement
 +
!width="250"| Improvement - Zoomed
 +
|-
 +
! Track Candidates
 +
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_bggen_PiMinus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_Candidates_bggen_PiMinus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_Candidates_bggen_PiMinus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_Candidates_bggen_PiMinus.png|left|275px]]
 +
|-
 +
! Time-Based Tracks
 +
| [[Image:Mattione_Update_09042013_Efficiency_TimeBased_bggen_PiMinus_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_TimeBased_bggen_PiMinus_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_TimeBased_bggen_PiMinus.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_TimeBased_bggen_PiMinus.png|left|275px]]
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
=== Proton ===
 +
* Thrown
 +
** Note this is in log-scale.
 +
[[Image:Mattione_Update_09042013_Thrown_Proton_bggen.png|thumb|left|275px|Thrown Proton's]]
 +
<br style="clear:both;"/>
 +
* Efficiencies
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="150"| Stage
 +
!width="250"| Efficiency - Current Code
 +
!width="250"| Efficiency - Spiral Code
 +
!width="250"| Improvement
 +
!width="250"| Improvement - Zoomed
 +
|-
 +
! Track Candidates
 +
| [[Image:Mattione_Update_09042013_Efficiency_Candidates_bggen_Proton_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_Candidates_bggen_Proton_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_Candidates_bggen_Proton.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_Candidates_bggen_Proton.png|left|275px]]
 +
|-
 +
! Time-Based Tracks
 +
| [[Image:Mattione_Update_09042013_Efficiency_TimeBased_bggen_Proton_Current.png|left|275px]] || [[Image:Mattione_Update_09042013_Efficiency_TimeBased_bggen_Proton_Spiral.png|left|275px]] || [[Image:Mattione_Update_09042013_EfficiencyDiff_TimeBased_bggen_Proton.png|left|275px]] ||
 +
[[Image:Mattione_Update_09042013_EfficiencyDiffZoomed_TimeBased_bggen_Proton.png|left|275px]]
 +
|}
 +
<br style="clear:both;"/>
 +
 
 +
= Hit Selector - bggen (pythia) Events =
 +
* New Code Only: Time-Based Track Efficiency improvements with hit selector disabled at all stages.
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="300"| Protons
 +
!width="300"| &pi;<sup>+</sup>'s
 +
!width="300"| &pi;<sup>-</sup>'s
 +
|-
 +
| [[Image:Mattione_Update_09042013_Efficiency_NoHitSelector_bggen_Proton.png|left|300px]] || [[Image:Mattione_Update_09042013_Efficiency_NoHitSelector_bggen_PiPlus.png|left|300px]] || [[Image:Mattione_Update_09042013_Efficiency_NoHitSelector_bggen_PiMinus.png|left|300px]]
 +
|}
 +
<br style="clear:both;"/>

Latest revision as of 13:26, 4 September 2013

CDC Track Finding - Basics

  • Note: Some steps (and many details) are omitted entirely in this brief overview.

Build Super Layer Seeds

  • In each ring, link adjacent hits together (DCDCRingSeed).
  • In each super layer, link each possible combination of adjacent DCDCRingSeed's together to form super layer seeds (DCDCSuperLayerSeed).
    • NEW: Do for both axial super layers and stereo super layers (previously just axial).
  • NEW: Identify which DCDCSuperLayerSeed's contain potential spiral turns (where a spiraling track turns back inward/outward).
b1pi Event
Spiral Event


Build Track Circles

  • HEAVILY MODIFIED: Link adjacent DCDCSuperLayerSeed's together to form tracks (DCDCTrackCircle).
    • Previously, the axial DCDCSuperLayerSeed's were linked by progressively fitting track circles and selecting the next axial super layers that were close to the track fit.
    • Start in super layer 1 (axial) and expand outwards towards super layer 7 (axial) (was previously in reverse direction)
    • For each unique combination of axial super layers, create a new DCDCTrackCircle object.
    • Keep track of which stereo DCDCSuperLayerSeed were used when linking the DCDCTrackCircle objects. There may be (many) more than one for each super layer.

Fit Track Circles & Truncate: Axial

b1pi Event
  • Fit the hits in the axial super layers (plus the beamline) to a circle with DHelicalFit.
    • Reject tracks for which half of the axial hits are not close to the circle fit.
  • NEW: Truncate and re-fit track circles that have picked up erroneous super layer seeds in the outer super layers.
    • This is used to save (for example) forward-going tracks that pick up bad hits from another track (e.g. a large-angle track near in φ)


Filter Track Circles: Axial

Spiral Track
  • If two track circles share >= 50% of their hits, reject the DCDCTrackCircle with the worse circle-fit weighted-chisq/ndf (weighted by #axial super layers squared)
  • NEW: If two track circles are spiral-linked in their last axial super layer, have opposite charges, and have similar circle centers, reject the DCDCTrackCircle with the worse circle-fit weighted-chisq/ndf.


Select Stereo Super Layer Seeds

b1pi Event
  • Previously, all stereo hits that intersected the circle and were nearby were lumped all together. Now, stereo super layer seeds are used instead.
  • Project all hits from all possible stereo super layers onto the circle.
  • For all possible combinations of stereo super layers:
    • Do a linear regression fit of the projected stereo hits to calculate θ & Vertex-Z
    • NEW: Select the combination of stereo super layer seeds that have the best chisq/ndf from the regression fit.


Filter Track Circles: Stereo

  • NEW: Find super layer seeds that are shared between DCDCTrackCircle's, and remove them from all DCDCTrackCircle's except the one with the best weighted chisq/ndf from the regression fit.
    • Like the axial truncation, this is used to save (for example) forward-going tracks that pick up bad hits from another track (e.g. a large-angle track near in φ).
    • However, do NOT remove the stereo super layer if it's the last one on the track.
  • NEW: If two track circles are spiral-linked in their last stereo super layer, have opposite charges, and have similar circle centers, reject the DCDCTrackCircle with the worse circle-fit weighted-chisq/ndf.

Refinement

b1pi Event
  • Find single, unused hits in the next (outward) super layer and add them to track circles if nearby.
    • This primarily helps tracks that exit the CDC early (e.g. forward going).
  • Redo circle fits, using both the new hits added just before, as well as the intersection points between the two stereo super layers (e.g. SL2 & SL3).
  • Re-project stereo hits onto the new circle fits.
  • Do a linear regression fit on a subset of the projected stereo hits to calculate θ & Vertex-Z
    • SEMI-NEW: A subset is chosen to optimize the final fit results. This is because there are regions/cases for which the projected stereo-hit-positions are bad.


Finalize

Spiral Events (2nd peak due to wrong charge)
  • Calculate the position and momentum of the track at a point outside of the target and the start counter.
    • This is because the track finding does not take energy-loss into account.
    • Note that this skews the comparison between the DTrackCandidate and DMCThrown kinematics, since they are reported at different positions (r ~= 9.5cm, r = 0cm, respectively).
  • Create DTrackCandidate objects.



Single π+ Track

Single π+ Reconstruction: Current (Trunk) Code, trackeff_hists plugin

  • The below plots show the track reconstruction efficiency and resolution for the trunk using the single-track scripts plugin (trackeff_hists).
  • However, this efficiency is extremely misleading:
    • If a track is reconstructed using any of the hits that were on the thrown track it is considered "found."
    • This is regardless of whether the reconstructed track momentum is anywhere close to the thrown values.
Current Code: Time-Based Tracks (trackeff_hists plugin)
Current Code: Momentum Resolution (trackeff_hists plugin)

Single π+ Reconstruction: Current (Trunk) Code, NEW trackeffv2 plugin

  • The below plots show the track candidate and time-based track reconstruction efficiency for the trunk using a new plugin (trackeff_hists).
  • This plugin compares matches all of the thrown and reconstructed particles in the event, and requires that their momentum be "close."
    • "Close" for track candidates is defined as: Δp/p < 20%, Δθ < 40 degrees, Δφ < 40 degrees, ΔVertex-Z < 1 km
      • For thrown tracks with θ < 5 degrees, there is no cut on the candidate Δφ
    • "Close" for time-based tracks is defined as: Δp/p < 10%, Δθ < 15 degrees, Δφ < 15 degrees, ΔVertex-Z < 10 cm
Current Code: Track Candidates (trackeffv2 plugin)
Current Code: Time-Based Tracks (trackeffv2 plugin)
  • Note that the degradation of the efficiency at forward angles matches the track resolution plot in the previous section.
    • If I expand the requirements for the "Close" window, the efficiency increases, and the chosen parameters were somewhat arbitrary.

Single π+ Reconstruction: NEW Spiral Code

  • The below plots show the track candidate and time-based track reconstruction efficiencies for the new spiral code.
  • Note: The trackeffv2 plugin compares matches all of the thrown and reconstructed particles in the event, and requires that their momentum be "close."
    • "Close" for track candidates is defined as: Δp/p < 20%, Δθ < 40 degrees, Δφ < 40 degrees, ΔVertex-Z < 1 km
      • For thrown tracks with θ < 5 degrees, there is no cut on the candidate Δφ
    • "Close" for time-based tracks is defined as: Δp/p < 10%, Δθ < 15 degrees, Δφ < 15 degrees, ΔVertex-Z < 10 cm
Spiral Code: Track Candidates Efficiency (Spiral Code)
Spiral Code: Time-Based Track Efficiency (Spiral Code)
  • Selected Track Candidate 1-D Kinematics Comparisons
Track Candidates (Forward Angle)
Track Candidates (Transition Region)
Track Candidates (Spirals)
  • The below plots show the improvements in the track candidate and time-based track reconstruction efficiencies with the new spiral code.
Track Candidate Efficiency Improvement
Zoomed Track Candidate Efficiency Improvement
Time-Based Track Efficiency Improvement
Zoomed Time-Based Track Efficiency Improvement


γp→π+π+π-(n)

  • The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
  • The definition of closeness is the same as before (see the single-track section).
  • Events were generated at beam energy of 9 GeV and a t-slope of 5.0

Summary

  • Overall track reconstruction improvement percentages for all tracks in this topology with momentum LESS than 3 GeV/c.
    • WARNING: This is not the overall gain for the topology, as many pions have larger than 3 GeV/c of momentum (their overall gain will be smaller).
Particle Track Candidates Time-Based Tracks
π+ 10.3% 5.38%
π- 8.08% 4.47%


π+

  • Thrown, Track Candidate 1-D Kinematics
Thrown π+'s
π+ Track Candidate Kinematics


  • Efficiencies
Stage Efficiency - Current Code Efficiency - Spiral Code Improvement Improvement - Zoomed
Track Candidates
Mattione Update 09042013 Efficiency Candidates n3pi PiPlus Current.png
Mattione Update 09042013 Efficiency Candidates n3pi PiPlus Spiral.png
Mattione Update 09042013 EfficiencyDiff Candidates n3pi PiPlus.png
Mattione Update 09042013 EfficiencyDiffZoomed Candidates n3pi PiPlus.png
Time-Based Tracks
Mattione Update 09042013 Efficiency TimeBased n3pi PiPlus Current.png
Mattione Update 09042013 Efficiency TimeBased n3pi PiPlus Spiral.png
Mattione Update 09042013 EfficiencyDiff TimeBased n3pi PiPlus.png
Mattione Update 09042013 EfficiencyDiffZoomed TimeBased n3pi PiPlus.png


γp→π+π+π-π-π0(p) (b1pi)

  • The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
  • The definition of closeness is the same as before (see the single-track section).
  • Events were generated with a beam energy of 9 GeV and a t-slope of 5.0

Summary

  • Overall track reconstruction improvement percentages for all tracks in this topology with momentum LESS than 3 GeV/c.
    • WARNING: This is not the overall gain for the topology, as many pions have larger than 3 GeV/c of momentum (their overall gain will be smaller).
Particle Track Candidates Time-Based Tracks
π+ 7.93% 5.21%
π- 12.5% 6.23%
p 23.7% 22.6%


π-

  • Thrown
Thrown π-'s


  • Efficiencies
Stage Efficiency - Current Code Efficiency - Spiral Code Improvement Improvement - Zoomed
Track Candidates
Mattione Update 09042013 Efficiency Candidates b1pi PiMinus Current.png
Mattione Update 09042013 Efficiency Candidates b1pi PiMinus Spiral.png
Mattione Update 09042013 EfficiencyDiff Candidates b1pi PiMinus.png
Mattione Update 09042013 EfficiencyDiffZoomed Candidates b1pi PiMinus.png
Time-Based Tracks
Mattione Update 09042013 Efficiency TimeBased b1pi PiMinus Current.png
Mattione Update 09042013 Efficiency TimeBased b1pi PiMinus Spiral.png
Mattione Update 09042013 EfficiencyDiff TimeBased b1pi PiMinus.png
Mattione Update 09042013 EfficiencyDiffZoomed TimeBased b1pi PiMinus.png


Proton

  • Thrown, Track Candidate 1-D Kinematics
Thrown Proton's
Track Candidates


  • Efficiencies
Stage Efficiency - Current Code Efficiency - Spiral Code Improvement
Track Candidates
Mattione Update 09042013 Efficiency Candidates b1pi Proton Current.png
Mattione Update 09042013 Efficiency Candidates b1pi Proton Spiral.png
Mattione Update 09042013 EfficiencyDiff Candidates b1pi Proton.png
Time-Based Tracks
Mattione Update 09042013 Efficiency TimeBased b1pi Proton Current.png
Mattione Update 09042013 Efficiency TimeBased b1pi Proton Spiral.png
Mattione Update 09042013 EfficiencyDiff TimeBased b1pi Proton.png


γp→K+K+Ξ-

  • The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
  • The definition of closeness is the same as before (see the single-track section).
  • Events were generated with beam energies between 8.4 and 9 GeV

Summary

  • Overall track reconstruction improvement percentages for all tracks in this topology with momentum LESS than 3 GeV/c.
    • WARNING: This is not the overall gain for the topology, as some kaons have larger than 3 GeV/c of momentum (their overall gain will be slightly smaller).
Particle Track Candidates Time-Based Tracks
K+ 21.7% 14.1%
π- 35.1% 28.7%
p 12.8% 16.7%


Sample Event

Reconstructed Hits Truth + Reconstructed Hits Current (09/03/13) Code Track Candidates Spiral Code Track Candidates
Mattione Update 09042013 Cascade Event3 ReconHits.png
Mattione Update 09042013 Cascade Event3 TruthHits.png
Mattione Update 09042013 Cascade Event3 Current.png
Mattione Update 09042013 Cascade Event3 Spiral.png


π-

  • Thrown
Thrown π-'s


  • Efficiencies
Stage Efficiency - Current Code Efficiency - Spiral Code Improvement Improvement - Zoomed
Track Candidates
Mattione Update 09042013 Efficiency Candidates cascade PiMinus Current.png
Mattione Update 09042013 Efficiency Candidates cascade PiMinus Spiral.png
Mattione Update 09042013 EfficiencyDiff Candidates cascade PiMinus.png
Mattione Update 09042013 EfficiencyDiffZoomed Candidates cascade PiMinus.png
Time-Based Tracks
Mattione Update 09042013 Efficiency TimeBased cascade PiMinus Current.png
Mattione Update 09042013 Efficiency TimeBased cascade PiMinus Spiral.png
Mattione Update 09042013 EfficiencyDiff TimeBased cascade PiMinus.png
Mattione Update 09042013 EfficiencyDiffZoomed TimeBased cascade PiMinus.png


K+

  • Thrown
Thrown K+'s


  • Efficiencies
Stage Efficiency - Current Code Efficiency - Spiral Code Improvement Improvement - Zoomed
Track Candidates
Mattione Update 09042013 Efficiency Candidates cascade KPlus Current.png
Mattione Update 09042013 Efficiency Candidates cascade KPlus Spiral.png
Mattione Update 09042013 EfficiencyDiff Candidates cascade KPlus.png
Mattione Update 09042013 EfficiencyDiffZoomed Candidates cascade KPlus.png
Time-Based Tracks
Mattione Update 09042013 Efficiency TimeBased cascade KPlus Current.png
Mattione Update 09042013 Efficiency TimeBased cascade KPlus Spiral.png
Mattione Update 09042013 EfficiencyDiff TimeBased cascade KPlus.png
Mattione Update 09042013 EfficiencyDiffZoomed TimeBased cascade KPlus.png


Proton

  • Thrown
Thrown Proton's


  • Efficiencies
Stage Efficiency - Current Code Efficiency - Spiral Code Improvement
Track Candidates
Mattione Update 09042013 Efficiency Candidates cascade Proton Current.png
Mattione Update 09042013 Efficiency Candidates cascade Proton Spiral.png
Mattione Update 09042013 EfficiencyDiff Candidates cascade Proton.png
Time-Based Tracks
Mattione Update 09042013 Efficiency TimeBased cascade Proton Current.png
Mattione Update 09042013 Efficiency TimeBased cascade Proton Spiral.png
Mattione Update 09042013 EfficiencyDiff TimeBased cascade Proton.png


bggen (pythia)

  • The below plots show the track candidate and time-based track reconstruction efficiencies for both track finding codes, and the improvement.
  • The definition of closeness is the same as before (see the single-track section).
  • Events were generated with beam energies between 8.4 and 9 GeV

Summary

  • Overall track reconstruction improvement percentages for all tracks in this topology with momentum LESS than 3 GeV/c.
    • WARNING: This is not the overall gain for the topology, as many pions have larger than 3 GeV/c of momentum (their overall gain will be smaller).
Particle Track Candidates Time-Based Tracks
π+ 16.7% 11.3%
π- 21.8% 12.7%
p 12.4% 6.92%


π-

  • Thrown
Thrown π-'s


  • Efficiencies
Stage Efficiency - Current Code Efficiency - Spiral Code Improvement Improvement - Zoomed
Track Candidates
Mattione Update 09042013 Efficiency Candidates bggen PiMinus Current.png
Mattione Update 09042013 Efficiency Candidates bggen PiMinus Spiral.png
Mattione Update 09042013 EfficiencyDiff Candidates bggen PiMinus.png
Mattione Update 09042013 EfficiencyDiffZoomed Candidates bggen PiMinus.png
Time-Based Tracks
Mattione Update 09042013 Efficiency TimeBased bggen PiMinus Current.png
Mattione Update 09042013 Efficiency TimeBased bggen PiMinus Spiral.png
Mattione Update 09042013 EfficiencyDiff TimeBased bggen PiMinus.png
Mattione Update 09042013 EfficiencyDiffZoomed TimeBased bggen PiMinus.png


Proton

  • Thrown
    • Note this is in log-scale.
Thrown Proton's


  • Efficiencies
Stage Efficiency - Current Code Efficiency - Spiral Code Improvement Improvement - Zoomed
Track Candidates
Mattione Update 09042013 Efficiency Candidates bggen Proton Current.png
Mattione Update 09042013 Efficiency Candidates bggen Proton Spiral.png
Mattione Update 09042013 EfficiencyDiff Candidates bggen Proton.png
Mattione Update 09042013 EfficiencyDiffZoomed Candidates bggen Proton.png
Time-Based Tracks
Mattione Update 09042013 Efficiency TimeBased bggen Proton Current.png
Mattione Update 09042013 Efficiency TimeBased bggen Proton Spiral.png
Mattione Update 09042013 EfficiencyDiff TimeBased bggen Proton.png
Mattione Update 09042013 EfficiencyDiffZoomed TimeBased bggen Proton.png


Hit Selector - bggen (pythia) Events

  • New Code Only: Time-Based Track Efficiency improvements with hit selector disabled at all stages.
Protons π+'s π-'s
Mattione Update 09042013 Efficiency NoHitSelector bggen Proton.png
Mattione Update 09042013 Efficiency NoHitSelector bggen PiPlus.png
Mattione Update 09042013 Efficiency NoHitSelector bggen PiMinus.png


Retrieved from "https://halldweb.jlab.org/wiki/index.php?title=Mattione_Update_09042013&oldid=50797"