Difference between revisions of "Mattione Update 04252012"
From GlueXWiki
(→Spiral Fitting - Procedure Overview) |
|||
(19 intermediate revisions by the same user not shown) | |||
Line 1: | Line 1: | ||
− | == | + | == Summary == |
+ | * Promising start to single-track event pion spiral reconstruction, MUCH work still to do. | ||
+ | == Spiral Fitting - Reconstructed Track Candidate Viewer Comparison == | ||
+ | |||
+ | * Notes: | ||
+ | ** All events have one π<sup>+</sup> thrown at p = 150 MeV/c, θ = 90 degrees, φ = random. | ||
+ | ** Current Code = Code as of 04/08/12 | ||
+ | ** Current Code # Reconstructed Track Candidates in events below: 4, 6, 5, 12 | ||
+ | ** Spiral Code # Reconstructed Track Candidates in events below: 2, 1, 0, 1 | ||
+ | ** Event 1: Spiral Code reconstructs two tracks due to turn near beamline, will be difficult to reduce to one track. | ||
+ | ** Event 2: Spiral Code reconstructs both spiral arms (from the beamline to the edge) and then links them together to form one track | ||
+ | ** Event 3: Not sure yet why spiral code fails to reconstruct anything. | ||
+ | ** Event 294: Current code ignores stereo hits while trying to link axial hits together by phi: explosion in track candidates | ||
+ | ** Event 294: Spiral code has memory explosion (leak?) while attempting to link axial and stereo hits together, need to improve somehow... | ||
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;" | {| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;" | ||
!width="100"| Event # | !width="100"| Event # | ||
!width="200"| Hits (CDC & Truth) | !width="200"| Hits (CDC & Truth) | ||
− | !width="200"| Current (04/08/12) Code | + | !width="200"| Current (04/08/12) Code Track Candidates |
− | !width="200"| Spiral Code | + | !width="200"| Spiral Code Track Candidates |
|- | |- | ||
! Event 1 | ! Event 1 | ||
Line 20: | Line 33: | ||
|} | |} | ||
<br style="clear:both;"/> | <br style="clear:both;"/> | ||
+ | |||
+ | == Spiral Fitting - Histograms == | ||
+ | * # Reconstructed Tracks / Event: | ||
+ | ** NOTE: Hit selectors disabled for wire based & time based fits: just use hits from track candidates (and ignore hits after spiral turn back inwards (for now anyway)) | ||
+ | ** Not sure why so few wire-based & time-based tracks reconstructed, possibly due to bad θ reconstruction. | ||
+ | {| class="wikitable" | ||
+ | | [[Image:Mattione_Update_04252012_SlowPiPlus_NumCandidates.gif|thumb|left|350px|Track Candidates]] | ||
+ | | [[Image:Mattione_Update_04252012_SlowPiPlus_NumWireBased.gif|thumb|left|350px|Wire Based Tracks]] | ||
+ | | [[Image:Mattione_Update_04252012_SlowPiPlus_NumTimeBased.gif|thumb|left|350px|Time Based Tracks]] | ||
+ | |} | ||
+ | <br style="clear:both;"/> | ||
+ | |||
+ | * Track Candidates: Track Parameter Differences (BEST Reconstructed - Thrown) | ||
+ | ** Best Reconstructed = Track Candidate that is the closest match to the thrown track. | ||
+ | ** NOTE: Δφ = 20 not real, due to artifact: reconstructed track candidate parameters are reported at a point along the track (~)just outside the start counter, whereas thrown parameters reported at thrown vertex (beamline) | ||
+ | ** Need to work on θ reconstruction. | ||
+ | {| class="wikitable" | ||
+ | | [[Image:Mattione_Update_04252012_SlowPiPlus_DeltaP.gif|thumb|left|350px|p<sub>Reconstructed</sub> - p<sub>Thrown</sub>]] | ||
+ | | [[Image:Mattione_Update_04252012_SlowPiPlus_DeltaTheta.gif|thumb|left|350px|θ<sub>Reconstructed</sub> - θ<sub>Thrown</sub>]] | ||
+ | | [[Image:Mattione_Update_04252012_SlowPiPlus_DeltaPhi.gif|thumb|left|350px|φ<sub>Reconstructed</sub> - φ<sub>Thrown</sub>]] | ||
+ | |} | ||
+ | |||
+ | == Spiral Fitting - Procedure Overview == | ||
+ | * Link adjacent hits together within each superlayer | ||
+ | * Link adjacent superlayer hits group together (all possible combinations) | ||
+ | * Fit the axial superlayer hits (and the beamline) to a circle (for each unique axial group combination) | ||
+ | * Reject bad axial group combinations (bad combo if one of its axial groups was already used in a fit that got a lower chisq/ndf) | ||
+ | * Compare the different stereo superlayer hit group combinations to the circle fit, choose the best combination (best = combo with lowest standard deviation of reconstructed z-vertex hits) | ||
+ | * Calculate track parameters (p, theta, phi, x, y, z) | ||
+ | * Link spiral arms together (look for rings with many consecutive hits shared between spiral arms) |
Latest revision as of 11:23, 25 April 2012
Contents
Summary
- Promising start to single-track event pion spiral reconstruction, MUCH work still to do.
Spiral Fitting - Reconstructed Track Candidate Viewer Comparison
- Notes:
- All events have one π+ thrown at p = 150 MeV/c, θ = 90 degrees, φ = random.
- Current Code = Code as of 04/08/12
- Current Code # Reconstructed Track Candidates in events below: 4, 6, 5, 12
- Spiral Code # Reconstructed Track Candidates in events below: 2, 1, 0, 1
- Event 1: Spiral Code reconstructs two tracks due to turn near beamline, will be difficult to reduce to one track.
- Event 2: Spiral Code reconstructs both spiral arms (from the beamline to the edge) and then links them together to form one track
- Event 3: Not sure yet why spiral code fails to reconstruct anything.
- Event 294: Current code ignores stereo hits while trying to link axial hits together by phi: explosion in track candidates
- Event 294: Spiral code has memory explosion (leak?) while attempting to link axial and stereo hits together, need to improve somehow...
Event # | Hits (CDC & Truth) | Current (04/08/12) Code Track Candidates | Spiral Code Track Candidates |
---|---|---|---|
Event 1 | |||
Event 2 | |||
Event 3 | |||
Event 294 |
Spiral Fitting - Histograms
- # Reconstructed Tracks / Event:
- NOTE: Hit selectors disabled for wire based & time based fits: just use hits from track candidates (and ignore hits after spiral turn back inwards (for now anyway))
- Not sure why so few wire-based & time-based tracks reconstructed, possibly due to bad θ reconstruction.
- Track Candidates: Track Parameter Differences (BEST Reconstructed - Thrown)
- Best Reconstructed = Track Candidate that is the closest match to the thrown track.
- NOTE: Δφ = 20 not real, due to artifact: reconstructed track candidate parameters are reported at a point along the track (~)just outside the start counter, whereas thrown parameters reported at thrown vertex (beamline)
- Need to work on θ reconstruction.
Spiral Fitting - Procedure Overview
- Link adjacent hits together within each superlayer
- Link adjacent superlayer hits group together (all possible combinations)
- Fit the axial superlayer hits (and the beamline) to a circle (for each unique axial group combination)
- Reject bad axial group combinations (bad combo if one of its axial groups was already used in a fit that got a lower chisq/ndf)
- Compare the different stereo superlayer hit group combinations to the circle fit, choose the best combination (best = combo with lowest standard deviation of reconstructed z-vertex hits)
- Calculate track parameters (p, theta, phi, x, y, z)
- Link spiral arms together (look for rings with many consecutive hits shared between spiral arms)