Difference: MuCapPKTMuStop (1 vs. 7)

Revision 72011-04-21 - PeterWinter

Line: 1 to 1
 
META TOPICPARENT name="MuCapPKSomeNotes"
Changed:
<
<
-- PeterKammel - 28 Apr 2009

IsTrack old and new

Old

Procedure

  • find contiguous EL island (1 gap) with EH
  • define muon stop location y=ymax_EH, findz +/- 1 time bin, Bragg = EH in Bragg peak

Criteria

  • zmax_EL - zmax_EH <=2
  • zmin_EH - zmin_EL >=5
  • EL in outer fiducial, EBragg in inner fiducial

classification.png

New

Procedure (StudyPixelLists)

  • Fit one line to green island
  • If chi2 too large, fit two line segments with kink to green island
  • ImproveTrackFit
    • poor 2 line chi2, repeat with Hough transformation
    • relax head cut?
  • Fit of extra green pixels
    • all EL with zEL > zmin_EL are candidates
    • with Hough selection of features and outliers, somewhat difference in reconstruction book
    • fit of Hough selected points
    • pointing check
  • Combining of two tracks into one track, if pointing

Scattering

  • HasScatterEH
    • 2 Blue islands
  • HasExtraELScatter
    • 1 EH island, >=5 ExtraEL
      • basic line fit, repeated after Hough transform rejection of outliers
  • IsTrackTrackScatterer
  • new TwoTracksScatter, for crossed tracks, aux bit 31

TTpcPixelListTrigger

Criteria { "IsValid", "HasCoincStrips", "IsGoodStop", "IsBetterStop", "IsGoldenStop", "IsSpot", "HasEVH", "IsNotAlone", "OldIsTrack", "HeaddlminEHCut", "TaildlminEHCut", "ChisqDofCut", "HasEHScatter", "HasExtraELScatter", "HasMWPCHitInROI", "HasUnreliableIsTrack"}

Property Required Comment
HasCoincStrips Y  
IsValid Y  
IsGoodStop Y  
IsBetterStop Y  
IsGoldenStop    
HeaddlminEHCut Y  
TaildlminEHCut Y  
ChisqDofCut Y  
IsSpot    
IsNotAlone    
HasEVH    
OldIsTrack    
HasEHScatter    
HasExtraELScatter    
HasMWPCHitInRoi    
HasUnreliableTrack   ask Brendan

Additional classification bits found in fAuxBits, some main properties

  • discuss with Brendan
  • which info on scatter track?

We should have similar classification figure as for old.

Get event display going on my computer.

Can one improve track fitting

Certainly, but what is the motivation?

TDC info

  • better approximation of "step function"

tracks.png

  • z shift, PeterW, Dave

FADC info

  • Improve z info with upstream anodes
  • better time pickoff for mu-capture events

Description of TMuStop Class


This class inherits from TMuEntrance, which inherits from TObject.

Since this object is intended to be written to a Tree, the space
required for the member variables was considered:
1) In some cases a desired Float_t (4 bytes) is instead stored as a
UShort_t (2 bytes), first multiplying the Float_t by 100 before
casting to a UShort_t. Functions are provided to make this
conversion invisible to the user
(ex., void SetStopAnode(Float_t anode), Float_t GetStopAnode(),
etc.).
2) TPC tracking information is more compactly stored in a "pixel
space," which corresponds to TDC400 clock ticks and TPC wire
numbers. The value fTDiffY0 is the time difference in ns
between the "YPixel=0" column and fT (the muon entrance time).
In order to allow for fractional pixel locations, in many cases
the pixel-space Float_t value is multiplied by 100 and stored as
a UShort_t.
3) fTpcPixelList is naturally in pixel-space.
4) fOneLine and fTwoLine are in pixel space.


fStopFlags (UShort_t)
Bit Description
1 "IsValid." Since initially each next-nearest-neighbor-connected
island with an EH pixel gets a TMuStop object, and some islands
may be pieced together in case of a scatter, islands that are
incorporated into another will have this flag set to zero.

2 "HasCoincStrips." There are strips EL pixels coincident with the
stop drift time.

3 "IsGoodStop." All EL points are within "GoodBox" fiducial boundary.

4 "IsBetterStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within "BetterBox" fiducial boundary. This is the
standard fiducial cut for lifetime histograms.

5 "IsGoldenStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within a much tighter fiducial boundary, "GoldenBox."

6 "IsSpot." The extent of this island in Y and Z is small.
if(fAnodeMaxEL - fAnodeMaxEH < 3 && fAnodeMinEH - fAnodeMinEL < 3
&& fYPixelMaxEL - fYPixelMaxEH < 4 && fYPixelMinEH - fYPixelMinEL < 3
&& fAnodeMaxEH - fAnodeMinEH < 4 && fYPixelMaxEH - fYPixelMinEH < 5)

7 "HasEVH." At least 1 EVH pixel somewhere in track.

8 "IsNotAlone." There is another TMuStop in the same ROI as this track.

9 "OldIsTrack." A very minimal cut, it means there are not too many
EL pixels downstream of the stop (fAnodeMaxEL - fAnodeMaxEH
<= mutrack_parameters.head_maxEL_maxEH_diff), and there are
sufficient pixels leading up to any EH hit (fAnodeMinEH -
fAnodeMinEL >= mutrack_parameters.tail_minEH_minEL_diff).

10 "HeaddlminEHCut." There is an EH pixel sufficiently close to the
endpoint of the fit line. dl = distance from endpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_head = the minimum dl of the EH pixels wrt endpoint of line.
This flag is set to 1 if dlmin_head <= mutrack_parameters.head_dlminEH.

11 "TaildlminEHCut." There is no EH pixel too close to the startpoint of
the fit line. dl = distance from startpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_tail = the minimum dl of the EH pixels wrt startpoint of line.
This flag is set to 1 if dlmin_tail >= mutrack_parameters.tail_dlminEH.

12 "ChisqDofCut." The pixels fit well to the line (one or two segments).
This flag is set to 1 if
chisqDof <= mutrack_parameters.isTrackChisqDof, where chisqDof is
ReducedChisq of the TwoLineSeg if present, otherwise of the LineSeg.

13 "HasEHScatter." Two-line fit, and there is an EH pixel near scatter
point: d < mutrack_parameters.maxDistScatterC_EH, where d is the
distance in the yz-plane (pixel space) of an EH pixel to the
inflection point of a TTwoLineSeg fit.

14 "HasExtraELScatter." (set and used in MTA)!ExtraEL pixels form a line
pointing at the apparent stop location of this muon.

15 "HasMWPCHitInROI." (set and used in MTA) Another track in the same
ROI as this muon overlaps the MWPC part of the drift region.

16 "HasUnreliableIsTrack." (set and used in MTA) This flag indicates some
anomaly with the IsTrack conditions, in particular a possible veto by
extra-EL pixels, such that the OldIsTrack flag should instead be used.

Is that still relevant?

ExtraMuonInfo.h

00014 unsigned int fTrackQualityFlags;

00015 // Further inspection of tracks for which pixel lists are stored:

00016 // bit1: part of this island is too low (within MWPC space); can probably ignore this track

00017 // bit2: IsTrack && IsNotAlone && (chisq of fitline <= 0.6) && GetNExtraPixelsEL<8

00018 // No reason to regard this track as troublesome.

00019 //

00020 // bit3: IsTrack && IsNotAlone && (chisq of fitline <= 0.6) && GetNExtraPixelsEL>=8

00021 // && (the extra pixels point to the stop position)

00022 // This muon may have scattered outside the fiducial volume.

00023 // bit4: same as bit3 except: the extra pixels do not point to the stop position

00024 // The extra EL points may be unrelated, or pileup, or alvarez...

The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.

UInt_t TMuStop::fAuxBits 32 flags for use in MTA

Bit Where Defined Assignment
1 ImproveTrackFit Attempted to improve 2-line fit
2 ImproveTrackFit HasUnreliableIsTrack bit set based on bad ChisqDof
3 ImproveTrackFit Relaxed HeaddlminEHCut to < 8.0
4 ImproveTrackFit HasUnreliableIsTrack bit set based on deltheta>2.9
OR (deltheta > 2.0 AND length2*sin(deltheta) < 5.)
5 ImproveTrackFit Trouble re-finding island in TpcROI (no action)
6 FindScatters This muon is leg1 (IsValid) of a track-track scatter.
7 FindScatters This muon is leg2 of a track-track scatter
8 FindScatters ExtraEL line fit used FitRejectingOutliers w/ dmax = 3.
9 ScatterHistFidCut EH-scatter, theta2 > 0 (down)
10 ScatterHistFidCut EH-scatter, theta2 <= 0 (up)
11 ScatterHistFidCut EH-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
12 ScatterHistFidCut ExtraEL-scatter, theta2 > 0 (down)
13 ScatterHistFidCut ExtraEL-scatter, theta2 <= 0 (up)
14 ScatterHistFidCut ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
15 ScatterHistFidCut MWPC-scatter, theta2 > 0 (down)
16 ScatterHistFidCut MWPC-scatter, theta2 <= 0 (up)
17 ScatterHistFidCut MWPC-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
18 ScatterHistFidCut TrackTrack-scatter, theta2 > 0 (down)
19 ScatterHistFidCut TrackTrack-scatter, theta2 <= 0 (up)
20 ScatterHistFidCut TrackTrack-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
21 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 > 0 (down)
22 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 <= 0 (up)
23 ScatterHistFidCut EH- OR ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
24 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 > 0 (down)
25 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 <= 0 (up)
26 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, 0.7 <= theta2 <= 1.5,
27 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 > 0 (down)
28 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 <= 0 (up)
29 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, 0.7 <= theta2 <= 1.5,
30
31 FindScatters TwoTracksScatter
32 FindScatters ExELScatter points to stop, line needed to be reversed

>
>
Please go to https://muon.npl.washington.edu/twiki/bin/view/Main/MuCapPKTMuStop
 
META FILEATTACHMENT attachment="classification.png" attr="" comment="" date="1241107688" name="classification.png" path="classification.png" size="49053" stream="classification.png" tmpFilename="/usr/tmp/CGItemp35310" user="PeterKammel" version="1"
META FILEATTACHMENT attachment="tracks.png" attr="" comment="" date="1241124228" name="tracks.png" path="tracks.png" size="43337" stream="tracks.png" tmpFilename="/usr/tmp/CGItemp35362" user="PeterKammel" version="1"

Revision 62009-05-13 - PeterKammel

Line: 1 to 1
 
META TOPICPARENT name="MuCapPKSomeNotes"
-- PeterKammel - 28 Apr 2009

IsTrack old and new

Line: 13 to 13
 
  • zmin_EH - zmin_EL >=5
  • EL in outer fiducial, EBragg in inner fiducial
Changed:
<
<
classification.png
>
>
classification.png
 

New

Line: 45 to 45
 
Changed:
<
<
Criteria
 {                 "IsValid",                                "HasCoincStrips",                         "IsGoodStop",                             "IsBetterStop",                           "IsGoldenStop",                           "IsSpot",                                 "HasEVH",                                 "IsNotAlone",                             "OldIsTrack",                             "HeaddlminEHCut",                         "TaildlminEHCut",                         "ChisqDofCut",                            "HasEHScatter",                           "HasExtraELScatter",                      "HasMWPCHitInROI",                        "HasUnreliableIsTrack"} 
>
>
Criteria { "IsValid", "HasCoincStrips", "IsGoodStop", "IsBetterStop", "IsGoldenStop", "IsSpot", "HasEVH", "IsNotAlone", "OldIsTrack", "HeaddlminEHCut", "TaildlminEHCut", "ChisqDofCut", "HasEHScatter", "HasExtraELScatter", "HasMWPCHitInROI", "HasUnreliableIsTrack"}
 
Property Required Comment
HasCoincStrips Y  
Line: 81 to 80
 

TDC info

  • better approximation of "step function"
Changed:
<
<
tracks.png
>
>
tracks.png
 
  • z shift, PeterW, Dave
Line: 95 to 94
  Is that still relevant?

ExtraMuonInfo.h

Changed:
<
<
00014   unsigned int [[http://webusers.npl.illinois.edu/~kammel/html/classExtraMuonInfo.html#d9dd09466577bd31805ced4d35027328][fTrackQualityFlags]];
00015 // Further inspection of tracks for which pixel lists are stored:
00016 //  bit1: part of this island is too low (within MWPC space); can probably ignore this track
00017 //  bit2: IsTrack && !IsNotAlone && (chisq of fitline <= 0.6) && GetNExtraPixelsEL<8
00018 //        No reason to regard this track as troublesome.
00019 //
00020 //  bit3: IsTrack && !IsNotAlone && (chisq of fitline <= 0.6) && GetNExtraPixelsEL>=8
00021 //           && (the extra pixels point to the stop position)
00022 //        This muon may have scattered outside the fiducial volume.
00023 //  bit4: same as bit3 except: the extra pixels do not point to the stop position
00024 //        The extra EL points may be unrelated, or pileup, or alvarez...
>
>
00014 unsigned int fTrackQualityFlags;

00015 // Further inspection of tracks for which pixel lists are stored:

00016 // bit1: part of this island is too low (within MWPC space); can probably ignore this track

00017 // bit2: IsTrack && IsNotAlone && (chisq of fitline <= 0.6) && GetNExtraPixelsEL<8

00018 // No reason to regard this track as troublesome.

00019 //

00020 // bit3: IsTrack && IsNotAlone && (chisq of fitline <= 0.6) && GetNExtraPixelsEL>=8

00021 // && (the extra pixels point to the stop position)

00022 // This muon may have scattered outside the fiducial volume.

00023 // bit4: same as bit3 except: the extra pixels do not point to the stop position

00024 // The extra EL points may be unrelated, or pileup, or alvarez...

  The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.

Revision 52009-05-09 - PeterKammel

Line: 1 to 1
 
META TOPICPARENT name="MuCapPKSomeNotes"
-- PeterKammel - 28 Apr 2009

IsTrack old and new

Line: 46 to 46
 

Criteria

Added:
>
>
 {                 "IsValid",                                "HasCoincStrips",                         "IsGoodStop",                             "IsBetterStop",                           "IsGoldenStop",                           "IsSpot",                                 "HasEVH",                                 "IsNotAlone",                             "OldIsTrack",                             "HeaddlminEHCut",                         "TaildlminEHCut",                         "ChisqDofCut",                            "HasEHScatter",                           "HasExtraELScatter",                      "HasMWPCHitInROI",                        "HasUnreliableIsTrack"} 
 
Property Required Comment
HasCoincStrips Y  
Line: 90 to 91
 

Description of TMuStop Class

Changed:
<
<

This class inherits from TMuEntrance, which inherits from TObject.

Since this object is intended to be written to a Tree, the space
required for the member variables was considered:
1) In some cases a desired Float_t (4 bytes) is instead stored as a
UShort_t (2 bytes), first multiplying the Float_t by 100 before
casting to a UShort_t. Functions are provided to make this
conversion invisible to the user
(ex., void SetStopAnode(Float_t anode), Float_t GetStopAnode(),
etc.).
2) TPC tracking information is more compactly stored in a "pixel
space," which corresponds to TDC400 clock ticks and TPC wire
numbers. The value fTDiffY0 is the time difference in ns
between the "YPixel=0" column and fT (the muon entrance time).
In order to allow for fractional pixel locations, in many cases
the pixel-space Float_t value is multiplied by 100 and stored as
a UShort_t.
3) fTpcPixelList is naturally in pixel-space.
4) fOneLine and fTwoLine are in pixel space.


fStopFlags (UShort_t)
Bit Description
1 "IsValid." Since initially each next-nearest-neighbor-connected
island with an EH pixel gets a TMuStop object, and some islands
may be pieced together in case of a scatter, islands that are
incorporated into another will have this flag set to zero.

2 "HasCoincStrips." There are strips EL pixels coincident with the
stop drift time.

3 "IsGoodStop." All EL points are within "GoodBox" fiducial boundary.

4 "IsBetterStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within "BetterBox" fiducial boundary. This is the
standard fiducial cut for lifetime histograms.

5 "IsGoldenStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within a much tighter fiducial boundary, "GoldenBox."

6 "IsSpot." The extent of this island in Y and Z is small.
if(fAnodeMaxEL - fAnodeMaxEH < 3 && fAnodeMinEH - fAnodeMinEL < 3
&& fYPixelMaxEL - fYPixelMaxEH < 4 && fYPixelMinEH - fYPixelMinEL < 3
&& fAnodeMaxEH - fAnodeMinEH < 4 && fYPixelMaxEH - fYPixelMinEH < 5)

7 "HasEVH." At least 1 EVH pixel somewhere in track.

8 "IsNotAlone." There is another TMuStop in the same ROI as this track.

9 "OldIsTrack." A very minimal cut, it means there are not too many
EL pixels downstream of the stop (fAnodeMaxEL - fAnodeMaxEH
<= mutrack_parameters.head_maxEL_maxEH_diff), and there are
sufficient pixels leading up to any EH hit (fAnodeMinEH -
fAnodeMinEL >= mutrack_parameters.tail_minEH_minEL_diff).

10 "HeaddlminEHCut." There is an EH pixel sufficiently close to the
endpoint of the fit line. dl = distance from endpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_head = the minimum dl of the EH pixels wrt endpoint of line.
This flag is set to 1 if dlmin_head <= mutrack_parameters.head_dlminEH.

11 "TaildlminEHCut." There is no EH pixel too close to the startpoint of
the fit line. dl = distance from startpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_tail = the minimum dl of the EH pixels wrt startpoint of line.
This flag is set to 1 if dlmin_tail >= mutrack_parameters.tail_dlminEH.

12 "ChisqDofCut." The pixels fit well to the line (one or two segments).
This flag is set to 1 if
chisqDof <= mutrack_parameters.isTrackChisqDof, where chisqDof is
ReducedChisq of the TwoLineSeg if present, otherwise of the LineSeg.

13 "HasEHScatter." Two-line fit, and there is an EH pixel near scatter
point: d < mutrack_parameters.maxDistScatterC_EH, where d is the
distance in the yz-plane (pixel space) of an EH pixel to the
inflection point of a TTwoLineSeg fit.

14 "HasExtraELScatter." (set and used in MTA)!ExtraEL pixels form a line
pointing at the apparent stop location of this muon.

15 "HasMWPCHitInROI." (set and used in MTA) Another track in the same
ROI as this muon overlaps the MWPC part of the drift region.

16 "HasUnreliableIsTrack." (set and used in MTA) This flag indicates some
anomaly with the IsTrack conditions, in particular a possible veto by
extra-EL pixels, such that the OldIsTrack flag should instead be used.


The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.
>
>

This class inherits from TMuEntrance, which inherits from TObject.

Since this object is intended to be written to a Tree, the space
required for the member variables was considered:
1) In some cases a desired Float_t (4 bytes) is instead stored as a
UShort_t (2 bytes), first multiplying the Float_t by 100 before
casting to a UShort_t. Functions are provided to make this
conversion invisible to the user
(ex., void SetStopAnode(Float_t anode), Float_t GetStopAnode(),
etc.).
2) TPC tracking information is more compactly stored in a "pixel
space," which corresponds to TDC400 clock ticks and TPC wire
numbers. The value fTDiffY0 is the time difference in ns
between the "YPixel=0" column and fT (the muon entrance time).
In order to allow for fractional pixel locations, in many cases
the pixel-space Float_t value is multiplied by 100 and stored as
a UShort_t.
3) fTpcPixelList is naturally in pixel-space.
4) fOneLine and fTwoLine are in pixel space.


fStopFlags (UShort_t)
Bit Description
1 "IsValid." Since initially each next-nearest-neighbor-connected
island with an EH pixel gets a TMuStop object, and some islands
may be pieced together in case of a scatter, islands that are
incorporated into another will have this flag set to zero.

2 "HasCoincStrips." There are strips EL pixels coincident with the
stop drift time.

3 "IsGoodStop." All EL points are within "GoodBox" fiducial boundary.

4 "IsBetterStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within "BetterBox" fiducial boundary. This is the
standard fiducial cut for lifetime histograms.

5 "IsGoldenStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within a much tighter fiducial boundary, "GoldenBox."

6 "IsSpot." The extent of this island in Y and Z is small.
if(fAnodeMaxEL - fAnodeMaxEH < 3 && fAnodeMinEH - fAnodeMinEL < 3
&& fYPixelMaxEL - fYPixelMaxEH < 4 && fYPixelMinEH - fYPixelMinEL < 3
&& fAnodeMaxEH - fAnodeMinEH < 4 && fYPixelMaxEH - fYPixelMinEH < 5)

7 "HasEVH." At least 1 EVH pixel somewhere in track.

8 "IsNotAlone." There is another TMuStop in the same ROI as this track.

9 "OldIsTrack." A very minimal cut, it means there are not too many
EL pixels downstream of the stop (fAnodeMaxEL - fAnodeMaxEH
<= mutrack_parameters.head_maxEL_maxEH_diff), and there are
sufficient pixels leading up to any EH hit (fAnodeMinEH -
fAnodeMinEL >= mutrack_parameters.tail_minEH_minEL_diff).

10 "HeaddlminEHCut." There is an EH pixel sufficiently close to the
endpoint of the fit line. dl = distance from endpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_head = the minimum dl of the EH pixels wrt endpoint of line.
This flag is set to 1 if dlmin_head <= mutrack_parameters.head_dlminEH.

11 "TaildlminEHCut." There is no EH pixel too close to the startpoint of
the fit line. dl = distance from startpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_tail = the minimum dl of the EH pixels wrt startpoint of line.
This flag is set to 1 if dlmin_tail >= mutrack_parameters.tail_dlminEH.

12 "ChisqDofCut." The pixels fit well to the line (one or two segments).
This flag is set to 1 if
chisqDof <= mutrack_parameters.isTrackChisqDof, where chisqDof is
ReducedChisq of the TwoLineSeg if present, otherwise of the LineSeg.

13 "HasEHScatter." Two-line fit, and there is an EH pixel near scatter
point: d < mutrack_parameters.maxDistScatterC_EH, where d is the
distance in the yz-plane (pixel space) of an EH pixel to the
inflection point of a TTwoLineSeg fit.

14 "HasExtraELScatter." (set and used in MTA)!ExtraEL pixels form a line
pointing at the apparent stop location of this muon.

15 "HasMWPCHitInROI." (set and used in MTA) Another track in the same
ROI as this muon overlaps the MWPC part of the drift region.

16 "HasUnreliableIsTrack." (set and used in MTA) This flag indicates some
anomaly with the IsTrack conditions, in particular a possible veto by
extra-EL pixels, such that the OldIsTrack flag should instead be used.

Is that still relevant?

ExtraMuonInfo.h

00014   unsigned int [[http://webusers.npl.illinois.edu/~kammel/html/classExtraMuonInfo.html#d9dd09466577bd31805ced4d35027328][fTrackQualityFlags]];
00015 // Further inspection of tracks for which pixel lists are stored:
00016 //  bit1: part of this island is too low (within MWPC space); can probably ignore this track
00017 //  bit2: IsTrack && !IsNotAlone && (chisq of fitline <= 0.6) && GetNExtraPixelsEL<8
00018 //        No reason to regard this track as troublesome.
00019 //
00020 //  bit3: IsTrack && !IsNotAlone && (chisq of fitline <= 0.6) && GetNExtraPixelsEL>=8
00021 //           && (the extra pixels point to the stop position)
00022 //        This muon may have scattered outside the fiducial volume.
00023 //  bit4: same as bit3 except: the extra pixels do not point to the stop position
00024 //        The extra EL points may be unrelated, or pileup, or alvarez...

The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.

  UInt_t TMuStop::fAuxBits 32 flags for use in MTA

Bit Where Defined Assignment
1 ImproveTrackFit Attempted to improve 2-line fit
2 ImproveTrackFit HasUnreliableIsTrack bit set based on bad ChisqDof
3 ImproveTrackFit Relaxed HeaddlminEHCut to < 8.0
4 ImproveTrackFit HasUnreliableIsTrack bit set based on deltheta>2.9
OR (deltheta > 2.0 AND length2*sin(deltheta) < 5.)
5 ImproveTrackFit Trouble re-finding island in TpcROI (no action)
6 FindScatters This muon is leg1 (IsValid) of a track-track scatter.
7 FindScatters This muon is leg2 of a track-track scatter
8 FindScatters ExtraEL line fit used FitRejectingOutliers w/ dmax = 3.
9 ScatterHistFidCut EH-scatter, theta2 > 0 (down)
10 ScatterHistFidCut EH-scatter, theta2 <= 0 (up)
11 ScatterHistFidCut EH-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
12 ScatterHistFidCut ExtraEL-scatter, theta2 > 0 (down)
13 ScatterHistFidCut ExtraEL-scatter, theta2 <= 0 (up)
14 ScatterHistFidCut ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
15 ScatterHistFidCut MWPC-scatter, theta2 > 0 (down)
16 ScatterHistFidCut MWPC-scatter, theta2 <= 0 (up)
17 ScatterHistFidCut MWPC-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
18 ScatterHistFidCut TrackTrack-scatter, theta2 > 0 (down)
19 ScatterHistFidCut TrackTrack-scatter, theta2 <= 0 (up)
20 ScatterHistFidCut TrackTrack-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
21 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 > 0 (down)
22 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 <= 0 (up)
23 ScatterHistFidCut EH- OR ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
24 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 > 0 (down)
25 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 <= 0 (up)
26 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, 0.7 <= theta2 <= 1.5,
27 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 > 0 (down)
28 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 <= 0 (up)
29 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, 0.7 <= theta2 <= 1.5,
30
31 FindScatters TwoTracksScatter
32 FindScatters ExELScatter points to stop, line needed to be reversed

Revision 42009-04-30 - PeterKammel

Line: 1 to 1
 
META TOPICPARENT name="MuCapPKSomeNotes"
-- PeterKammel - 28 Apr 2009

IsTrack old and new

Line: 71 to 71
  We should have similar classification figure as for old.
Changed:
<
<
Read elogs again. Ask for more picture books or install on my computer.
>
>
Get event display going on my computer.

Can one improve track fitting

Certainly, but what is the motivation?

TDC info

  • better approximation of "step function"

tracks.png

  • z shift, PeterW, Dave

FADC info

  • Improve z info with upstream anodes
  • better time pickoff for mu-capture events
 

Description of TMuStop Class

Line: 80 to 95
 UInt_t TMuStop::fAuxBits 32 flags for use in MTA

Bit Where Defined Assignment
1 ImproveTrackFit Attempted to improve 2-line fit
2 ImproveTrackFit HasUnreliableIsTrack bit set based on bad ChisqDof
3 ImproveTrackFit Relaxed HeaddlminEHCut to < 8.0
4 ImproveTrackFit HasUnreliableIsTrack bit set based on deltheta>2.9
OR (deltheta > 2.0 AND length2*sin(deltheta) < 5.)
5 ImproveTrackFit Trouble re-finding island in TpcROI (no action)
6 FindScatters This muon is leg1 (IsValid) of a track-track scatter.
7 FindScatters This muon is leg2 of a track-track scatter
8 FindScatters ExtraEL line fit used FitRejectingOutliers w/ dmax = 3.
9 ScatterHistFidCut EH-scatter, theta2 > 0 (down)
10 ScatterHistFidCut EH-scatter, theta2 <= 0 (up)
11 ScatterHistFidCut EH-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
12 ScatterHistFidCut ExtraEL-scatter, theta2 > 0 (down)
13 ScatterHistFidCut ExtraEL-scatter, theta2 <= 0 (up)
14 ScatterHistFidCut ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
15 ScatterHistFidCut MWPC-scatter, theta2 > 0 (down)
16 ScatterHistFidCut MWPC-scatter, theta2 <= 0 (up)
17 ScatterHistFidCut MWPC-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
18 ScatterHistFidCut TrackTrack-scatter, theta2 > 0 (down)
19 ScatterHistFidCut TrackTrack-scatter, theta2 <= 0 (up)
20 ScatterHistFidCut TrackTrack-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
21 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 > 0 (down)
22 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 <= 0 (up)
23 ScatterHistFidCut EH- OR ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
24 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 > 0 (down)
25 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 <= 0 (up)
26 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, 0.7 <= theta2 <= 1.5,
27 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 > 0 (down)
28 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 <= 0 (up)
29 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, 0.7 <= theta2 <= 1.5,
30
31 FindScatters TwoTracksScatter
32 FindScatters ExELScatter points to stop, line needed to be reversed

META FILEATTACHMENT attachment="classification.png" attr="" comment="" date="1241107688" name="classification.png" path="classification.png" size="49053" stream="classification.png" tmpFilename="/usr/tmp/CGItemp35310" user="PeterKammel" version="1"
Added:
>
>
META FILEATTACHMENT attachment="tracks.png" attr="" comment="" date="1241124228" name="tracks.png" path="tracks.png" size="43337" stream="tracks.png" tmpFilename="/usr/tmp/CGItemp35362" user="PeterKammel" version="1"

Revision 32009-04-30 - PeterKammel

Line: 1 to 1
 
META TOPICPARENT name="MuCapPKSomeNotes"
-- PeterKammel - 28 Apr 2009

IsTrack old and new

Line: 13 to 13
 
  • zmin_EH - zmin_EL >=5
  • EL in outer fiducial, EBragg in inner fiducial
Added:
>
>
classification.png
 

New

Changed:
<
<
Procedure
>
>
Procedure (StudyPixelLists)
 
  • Fit one line to green island
  • If chi2 too large, fit two line segments with kink to green island
Added:
>
>
  • ImproveTrackFit
    • poor 2 line chi2, repeat with Hough transformation
    • relax head cut?
 
  • Fit of extra green pixels
    • all EL with zEL > zmin_EL are candidates
    • with Hough selection of features and outliers, somewhat difference in reconstruction book
Line: 25 to 30
 
    • pointing check
  • Combining of two tracks into one track, if pointing
Added:
>
>
Scattering
  • HasScatterEH
    • 2 Blue islands
  • HasExtraELScatter
    • 1 EH island, >=5 ExtraEL
      • basic line fit, repeated after Hough transform rejection of outliers
  • IsTrackTrackScatterer
  • new TwoTracksScatter, for crossed tracks, aux bit 31

TTpcPixelListTrigger

 Criteria

Property Required Comment
Line: 47 to 67
  Additional classification bits found in fAuxBits, some main properties
  • discuss with Brendan
Added:
>
>
  • which info on scatter track?

We should have similar classification figure as for old.

  Read elogs again. Ask for more picture books or install on my computer.
Changed:
<
<

Description of TMuStop Class

>
>

Description of TMuStop Class

 
Changed:
<
<

This class inherits from TMuEntrance, which inherits from TObject.

Since this object is intended to be written to a Tree, the space
required for the member variables was considered:
1) In some cases a desired Float_t (4 bytes) is instead stored as a
UShort_t (2 bytes), first multiplying the Float_t by 100 before
casting to a UShort_t. Functions are provided to make this
conversion invisible to the user
(ex., void SetStopAnode(Float_t anode), Float_t GetStopAnode(),
etc.).
2) TPC tracking information is more compactly stored in a "pixel
space," which corresponds to TDC400 clock ticks and TPC wire
numbers. The value fTDiffY0 is the time difference in ns
between the "YPixel=0" column and fT (the muon entrance time).
In order to allow for fractional pixel locations, in many cases
the pixel-space Float_t value is multiplied by 100 and stored as
a UShort_t.
3) fTpcPixelList is naturally in pixel-space.
4) fOneLine and fTwoLine are in pixel space.


fStopFlags (UShort_t)
Bit Description
1 "IsValid." Since initially each next-nearest-neighbor-connected
island with an EH pixel gets a TMuStop object, and some islands
may be pieced together in case of a scatter, islands that are
incorporated into another will have this flag set to zero.

2 "HasCoincStrips." There are strips EL pixels coincident with the
stop drift time.

3 "IsGoodStop." All EL points are within "GoodBox" fiducial boundary.

4 "IsBetterStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within "BetterBox" fiducial boundary. This is the
standard fiducial cut for lifetime histograms.

5 "IsGoldenStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within a much tighter fiducial boundary, "GoldenBox."

6 "IsSpot." The extent of this island in Y and Z is small.
if(fAnodeMaxEL - fAnodeMaxEH < 3 && fAnodeMinEH - fAnodeMinEL < 3
&& fYPixelMaxEL - fYPixelMaxEH < 4 && fYPixelMinEH - fYPixelMinEL < 3
&& fAnodeMaxEH - fAnodeMinEH < 4 && fYPixelMaxEH - fYPixelMinEH < 5)

7 "HasEVH." At least 1 EVH pixel somewhere in track.

8 "IsNotAlone." There is another TMuStop in the same ROI as this track.

9 "OldIsTrack." A very minimal cut, it means there are not too many
EL pixels downstream of the stop (fAnodeMaxEL - fAnodeMaxEH
<= mutrack_parameters.head_maxEL_maxEH_diff), and there are
sufficient pixels leading up to any EH hit (fAnodeMinEH -
fAnodeMinEL >= mutrack_parameters.tail_minEH_minEL_diff).

10 "HeaddlminEHCut." There is an EH pixel sufficiently close to the
endpoint of the fit line. dl = distance from endpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_head = the minimum dl of the EH pixels wrt endpoint of line.
This flag is set to 1 if dlmin_head <= mutrack_parameters.head_dlminEH.

11 "TaildlminEHCut." There is no EH pixel too close to the startpoint of
the fit line. dl = distance from startpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_tail = the minimum dl of the EH pixels wrt startpoint of line.
This flag is set to 1 if dlmin_tail >= mutrack_parameters.tail_dlminEH.

12 "ChisqDofCut." The pixels fit well to the line (one or two segments).
This flag is set to 1 if
chisqDof <= mutrack_parameters.isTrackChisqDof, where chisqDof is
ReducedChisq of the TwoLineSeg if present, otherwise of the LineSeg.

13 "HasEHScatter." Two-line fit, and there is an EH pixel near scatter
point: d < mutrack_parameters.maxDistScatterC_EH, where d is the
distance in the yz-plane (pixel space) of an EH pixel to the
inflection point of a TTwoLineSeg fit.

14 "HasExtraELScatter." (set and used in MTA) ExtraEL pixels form a line
pointing at the apparent stop location of this muon.

15 "HasMWPCHitInROI." (set and used in MTA) Another track in the same
ROI as this muon overlaps the MWPC part of the drift region.

16 "HasUnreliableIsTrack." (set and used in MTA) This flag indicates some
anomaly with the IsTrack conditions, in particular a possible veto by
extra-EL pixels, such that the OldIsTrack flag should instead be used.


The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.
>
>

This class inherits from TMuEntrance, which inherits from TObject.

Since this object is intended to be written to a Tree, the space
required for the member variables was considered:
1) In some cases a desired Float_t (4 bytes) is instead stored as a
UShort_t (2 bytes), first multiplying the Float_t by 100 before
casting to a UShort_t. Functions are provided to make this
conversion invisible to the user
(ex., void SetStopAnode(Float_t anode), Float_t GetStopAnode(),
etc.).
2) TPC tracking information is more compactly stored in a "pixel
space," which corresponds to TDC400 clock ticks and TPC wire
numbers. The value fTDiffY0 is the time difference in ns
between the "YPixel=0" column and fT (the muon entrance time).
In order to allow for fractional pixel locations, in many cases
the pixel-space Float_t value is multiplied by 100 and stored as
a UShort_t.
3) fTpcPixelList is naturally in pixel-space.
4) fOneLine and fTwoLine are in pixel space.


fStopFlags (UShort_t)
Bit Description
1 "IsValid." Since initially each next-nearest-neighbor-connected
island with an EH pixel gets a TMuStop object, and some islands
may be pieced together in case of a scatter, islands that are
incorporated into another will have this flag set to zero.

2 "HasCoincStrips." There are strips EL pixels coincident with the
stop drift time.

3 "IsGoodStop." All EL points are within "GoodBox" fiducial boundary.

4 "IsBetterStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within "BetterBox" fiducial boundary. This is the
standard fiducial cut for lifetime histograms.

5 "IsGoldenStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within a much tighter fiducial boundary, "GoldenBox."

6 "IsSpot." The extent of this island in Y and Z is small.
if(fAnodeMaxEL - fAnodeMaxEH < 3 && fAnodeMinEH - fAnodeMinEL < 3
&& fYPixelMaxEL - fYPixelMaxEH < 4 && fYPixelMinEH - fYPixelMinEL < 3
&& fAnodeMaxEH - fAnodeMinEH < 4 && fYPixelMaxEH - fYPixelMinEH < 5)

7 "HasEVH." At least 1 EVH pixel somewhere in track.

8 "IsNotAlone." There is another TMuStop in the same ROI as this track.

9 "OldIsTrack." A very minimal cut, it means there are not too many
EL pixels downstream of the stop (fAnodeMaxEL - fAnodeMaxEH
<= mutrack_parameters.head_maxEL_maxEH_diff), and there are
sufficient pixels leading up to any EH hit (fAnodeMinEH -
fAnodeMinEL >= mutrack_parameters.tail_minEH_minEL_diff).

10 "HeaddlminEHCut." There is an EH pixel sufficiently close to the
endpoint of the fit line. dl = distance from endpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_head = the minimum dl of the EH pixels wrt endpoint of line.
This flag is set to 1 if dlmin_head <= mutrack_parameters.head_dlminEH.

11 "TaildlminEHCut." There is no EH pixel too close to the startpoint of
the fit line. dl = distance from startpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_tail = the minimum dl of the EH pixels wrt startpoint of line.
This flag is set to 1 if dlmin_tail >= mutrack_parameters.tail_dlminEH.

12 "ChisqDofCut." The pixels fit well to the line (one or two segments).
This flag is set to 1 if
chisqDof <= mutrack_parameters.isTrackChisqDof, where chisqDof is
ReducedChisq of the TwoLineSeg if present, otherwise of the LineSeg.

13 "HasEHScatter." Two-line fit, and there is an EH pixel near scatter
point: d < mutrack_parameters.maxDistScatterC_EH, where d is the
distance in the yz-plane (pixel space) of an EH pixel to the
inflection point of a TTwoLineSeg fit.

14 "HasExtraELScatter." (set and used in MTA)!ExtraEL pixels form a line
pointing at the apparent stop location of this muon.

15 "HasMWPCHitInROI." (set and used in MTA) Another track in the same
ROI as this muon overlaps the MWPC part of the drift region.

16 "HasUnreliableIsTrack." (set and used in MTA) This flag indicates some
anomaly with the IsTrack conditions, in particular a possible veto by
extra-EL pixels, such that the OldIsTrack flag should instead be used.


The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.
  UInt_t TMuStop::fAuxBits 32 flags for use in MTA

Bit Where Defined Assignment
1 ImproveTrackFit Attempted to improve 2-line fit
2 ImproveTrackFit HasUnreliableIsTrack bit set based on bad ChisqDof
3 ImproveTrackFit Relaxed HeaddlminEHCut to < 8.0
4 ImproveTrackFit HasUnreliableIsTrack bit set based on deltheta>2.9
OR (deltheta > 2.0 AND length2*sin(deltheta) < 5.)
5 ImproveTrackFit Trouble re-finding island in TpcROI (no action)
6 FindScatters This muon is leg1 (IsValid) of a track-track scatter.
7 FindScatters This muon is leg2 of a track-track scatter
8 FindScatters ExtraEL line fit used FitRejectingOutliers w/ dmax = 3.
9 ScatterHistFidCut EH-scatter, theta2 > 0 (down)
10 ScatterHistFidCut EH-scatter, theta2 <= 0 (up)
11 ScatterHistFidCut EH-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
12 ScatterHistFidCut ExtraEL-scatter, theta2 > 0 (down)
13 ScatterHistFidCut ExtraEL-scatter, theta2 <= 0 (up)
14 ScatterHistFidCut ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
15 ScatterHistFidCut MWPC-scatter, theta2 > 0 (down)
16 ScatterHistFidCut MWPC-scatter, theta2 <= 0 (up)
17 ScatterHistFidCut MWPC-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
18 ScatterHistFidCut TrackTrack-scatter, theta2 > 0 (down)
19 ScatterHistFidCut TrackTrack-scatter, theta2 <= 0 (up)
20 ScatterHistFidCut TrackTrack-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
21 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 > 0 (down)
22 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 <= 0 (up)
23 ScatterHistFidCut EH- OR ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
24 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 > 0 (down)
25 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 <= 0 (up)
26 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, 0.7 <= theta2 <= 1.5,
27 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 > 0 (down)
28 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 <= 0 (up)
29 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, 0.7 <= theta2 <= 1.5,
30
31 FindScatters TwoTracksScatter
32 FindScatters ExELScatter points to stop, line needed to be reversed
Added:
>
>
META FILEATTACHMENT attachment="classification.png" attr="" comment="" date="1241107688" name="classification.png" path="classification.png" size="49053" stream="classification.png" tmpFilename="/usr/tmp/CGItemp35310" user="PeterKammel" version="1"

Revision 22009-04-30 - PeterKammel

Line: 1 to 1
 
META TOPICPARENT name="MuCapPKSomeNotes"
-- PeterKammel - 28 Apr 2009
Added:
>
>

IsTrack old and new

Old

 
Changed:
<
<
Description of TMuStop Class

This class inherits from TMuEntrance, which inherits from TObject.

Since this object is intended to be written to a Tree, the space
required for the member variables was considered:
1) In some cases a desired Float_t (4 bytes) is instead stored as a
UShort_t (2 bytes), first multiplying the Float_t by 100 before
casting to a UShort_t. Functions are provided to make this
conversion invisible to the user
(ex., void SetStopAnode(Float_t anode), Float_t GetStopAnode(),
etc.).
2) TPC tracking information is more compactly stored in a "pixel
space," which corresponds to TDC400 clock ticks and TPC wire
numbers. The value fTDiffY0 is the time difference in ns
between the "YPixel=0" column and fT (the muon entrance time).
In order to allow for fractional pixel locations, in many cases
the pixel-space Float_t value is multiplied by 100 and stored as
a UShort_t.
3) fTpcPixelList is naturally in pixel-space.
4) fOneLine and fTwoLine are in pixel space.


fStopFlags (UShort_t)
Bit Description
1 "IsValid." Since initially each next-nearest-neighbor-connected
island with an EH pixel gets a TMuStop object, and some islands
may be pieced together in case of a scatter, islands that are
incorporated into another will have this flag set to zero.

2 "HasCoincStrips." There are strips EL pixels coincident with the
stop drift time.

3 "IsGoodStop." All EL points are within "GoodBox" fiducial boundary.

4 "IsBetterStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within "BetterBox" fiducial boundary. This is the
standard fiducial cut for lifetime histograms.

5 "IsGoldenStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within a much tighter fiducial boundary, "GoldenBox."

6 "IsSpot." The extent of this island in Y and Z is small.
if(fAnodeMaxEL - fAnodeMaxEH < 3 && fAnodeMinEH - fAnodeMinEL < 3
&& fYPixelMaxEL - fYPixelMaxEH < 4 && fYPixelMinEH - fYPixelMinEL < 3
&& fAnodeMaxEH - fAnodeMinEH < 4 && fYPixelMaxEH - fYPixelMinEH < 5)

7 "HasEVH." At least 1 EVH pixel somewhere in track.

8 "IsNotAlone." There is another TMuStop in the same ROI as this track.

9 "OldIsTrack." A very minimal cut, it means there are not too many
EL pixels downstream of the stop (fAnodeMaxEL - fAnodeMaxEH
<= mutrack_parameters.head_maxEL_maxEH_diff), and there are
sufficient pixels leading up to any EH hit (fAnodeMinEH -
fAnodeMinEL >= mutrack_parameters.tail_minEH_minEL_diff).

10 "HeaddlminEHCut." There is an EH pixel sufficiently close to the
endpoint of the fit line. dl = distance from endpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_head = the minimum dl of the EH pixels wrt endpoint of line.
This flag is set to 1 if dlmin_head <= mutrack_parameters.head_dlminEH.

11 "TaildlminEHCut." There is no EH pixel too close to the startpoint of
the fit line. dl = distance from startpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_tail = the minimum dl of the EH pixels wrt startpoint of line.
This flag is set to 1 if dlmin_tail >= mutrack_parameters.tail_dlminEH.

12 "ChisqDofCut." The pixels fit well to the line (one or two segments).
This flag is set to 1 if
chisqDof <= mutrack_parameters.isTrackChisqDof, where chisqDof is
ReducedChisq of the TwoLineSeg if present, otherwise of the LineSeg.

13 "HasEHScatter." Two-line fit, and there is an EH pixel near scatter
point: d < mutrack_parameters.maxDistScatterC_EH, where d is the
distance in the yz-plane (pixel space) of an EH pixel to the
inflection point of a TTwoLineSeg fit.

14 "HasExtraELScatter." (set and used in MTA) ExtraEL pixels form a line
pointing at the apparent stop location of this muon.

15 "HasMWPCHitInROI." (set and used in MTA) Another track in the same
ROI as this muon overlaps the MWPC part of the drift region.

16 "HasUnreliableIsTrack." (set and used in MTA) This flag indicates some
anomaly with the IsTrack conditions, in particular a possible veto by
extra-EL pixels, such that the OldIsTrack flag should instead be used.


The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.
>
>
Procedure
  • find contiguous EL island (1 gap) with EH
  • define muon stop location y=ymax_EH, findz +/- 1 time bin, Bragg = EH in Bragg peak

Criteria

  • zmax_EL - zmax_EH <=2
  • zmin_EH - zmin_EL >=5
  • EL in outer fiducial, EBragg in inner fiducial

New

Procedure

  • Fit one line to green island
  • If chi2 too large, fit two line segments with kink to green island
  • Fit of extra green pixels
    • all EL with zEL > zmin_EL are candidates
    • with Hough selection of features and outliers, somewhat difference in reconstruction book
    • fit of Hough selected points
    • pointing check
  • Combining of two tracks into one track, if pointing

Criteria

Property Required Comment
HasCoincStrips Y  
IsValid Y  
IsGoodStop Y  
IsBetterStop Y  
IsGoldenStop    
HeaddlminEHCut Y  
TaildlminEHCut Y  
ChisqDofCut Y  
IsSpot    
IsNotAlone    
HasEVH    
OldIsTrack    
HasEHScatter    
HasExtraELScatter    
HasMWPCHitInRoi    
HasUnreliableTrack   ask Brendan

Additional classification bits found in fAuxBits, some main properties

  • discuss with Brendan

Read elogs again. Ask for more picture books or install on my computer.

Description of TMuStop Class


This class inherits from TMuEntrance, which inherits from TObject.

Since this object is intended to be written to a Tree, the space
required for the member variables was considered:
1) In some cases a desired Float_t (4 bytes) is instead stored as a
UShort_t (2 bytes), first multiplying the Float_t by 100 before
casting to a UShort_t. Functions are provided to make this
conversion invisible to the user
(ex., void SetStopAnode(Float_t anode), Float_t GetStopAnode(),
etc.).
2) TPC tracking information is more compactly stored in a "pixel
space," which corresponds to TDC400 clock ticks and TPC wire
numbers. The value fTDiffY0 is the time difference in ns
between the "YPixel=0" column and fT (the muon entrance time).
In order to allow for fractional pixel locations, in many cases
the pixel-space Float_t value is multiplied by 100 and stored as
a UShort_t.
3) fTpcPixelList is naturally in pixel-space.
4) fOneLine and fTwoLine are in pixel space.


fStopFlags (UShort_t)
Bit Description
1 "IsValid." Since initially each next-nearest-neighbor-connected
island with an EH pixel gets a TMuStop object, and some islands
may be pieced together in case of a scatter, islands that are
incorporated into another will have this flag set to zero.

2 "HasCoincStrips." There are strips EL pixels coincident with the
stop drift time.

3 "IsGoodStop." All EL points are within "GoodBox" fiducial boundary.

4 "IsBetterStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within "BetterBox" fiducial boundary. This is the
standard fiducial cut for lifetime histograms.

5 "IsGoldenStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within a much tighter fiducial boundary, "GoldenBox."

6 "IsSpot." The extent of this island in Y and Z is small.
if(fAnodeMaxEL - fAnodeMaxEH < 3 && fAnodeMinEH - fAnodeMinEL < 3
&& fYPixelMaxEL - fYPixelMaxEH < 4 && fYPixelMinEH - fYPixelMinEL < 3
&& fAnodeMaxEH - fAnodeMinEH < 4 && fYPixelMaxEH - fYPixelMinEH < 5)

7 "HasEVH." At least 1 EVH pixel somewhere in track.

8 "IsNotAlone." There is another TMuStop in the same ROI as this track.

9 "OldIsTrack." A very minimal cut, it means there are not too many
EL pixels downstream of the stop (fAnodeMaxEL - fAnodeMaxEH
<= mutrack_parameters.head_maxEL_maxEH_diff), and there are
sufficient pixels leading up to any EH hit (fAnodeMinEH -
fAnodeMinEL >= mutrack_parameters.tail_minEH_minEL_diff).

10 "HeaddlminEHCut." There is an EH pixel sufficiently close to the
endpoint of the fit line. dl = distance from endpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_head = the minimum dl of the EH pixels wrt endpoint of line.
This flag is set to 1 if dlmin_head <= mutrack_parameters.head_dlminEH.

11 "TaildlminEHCut." There is no EH pixel too close to the startpoint of
the fit line. dl = distance from startpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_tail = the minimum dl of the EH pixels wrt startpoint of line.
This flag is set to 1 if dlmin_tail >= mutrack_parameters.tail_dlminEH.

12 "ChisqDofCut." The pixels fit well to the line (one or two segments).
This flag is set to 1 if
chisqDof <= mutrack_parameters.isTrackChisqDof, where chisqDof is
ReducedChisq of the TwoLineSeg if present, otherwise of the LineSeg.

13 "HasEHScatter." Two-line fit, and there is an EH pixel near scatter
point: d < mutrack_parameters.maxDistScatterC_EH, where d is the
distance in the yz-plane (pixel space) of an EH pixel to the
inflection point of a TTwoLineSeg fit.

14 "HasExtraELScatter." (set and used in MTA) ExtraEL pixels form a line
pointing at the apparent stop location of this muon.

15 "HasMWPCHitInROI." (set and used in MTA) Another track in the same
ROI as this muon overlaps the MWPC part of the drift region.

16 "HasUnreliableIsTrack." (set and used in MTA) This flag indicates some
anomaly with the IsTrack conditions, in particular a possible veto by
extra-EL pixels, such that the OldIsTrack flag should instead be used.


The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.

  UInt_t TMuStop::fAuxBits 32 flags for use in MTA

Bit Where Defined Assignment
1 ImproveTrackFit Attempted to improve 2-line fit
2 ImproveTrackFit HasUnreliableIsTrack bit set based on bad ChisqDof
3 ImproveTrackFit Relaxed HeaddlminEHCut to < 8.0
4 ImproveTrackFit HasUnreliableIsTrack bit set based on deltheta>2.9
OR (deltheta > 2.0 AND length2*sin(deltheta) < 5.)
5 ImproveTrackFit Trouble re-finding island in TpcROI (no action)
6 FindScatters This muon is leg1 (IsValid) of a track-track scatter.
7 FindScatters This muon is leg2 of a track-track scatter
8 FindScatters ExtraEL line fit used FitRejectingOutliers w/ dmax = 3.
9 ScatterHistFidCut EH-scatter, theta2 > 0 (down)
10 ScatterHistFidCut EH-scatter, theta2 <= 0 (up)
11 ScatterHistFidCut EH-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
12 ScatterHistFidCut ExtraEL-scatter, theta2 > 0 (down)
13 ScatterHistFidCut ExtraEL-scatter, theta2 <= 0 (up)
14 ScatterHistFidCut ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
15 ScatterHistFidCut MWPC-scatter, theta2 > 0 (down)
16 ScatterHistFidCut MWPC-scatter, theta2 <= 0 (up)
17 ScatterHistFidCut MWPC-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
18 ScatterHistFidCut TrackTrack-scatter, theta2 > 0 (down)
19 ScatterHistFidCut TrackTrack-scatter, theta2 <= 0 (up)
20 ScatterHistFidCut TrackTrack-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
21 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 > 0 (down)
22 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 <= 0 (up)
23 ScatterHistFidCut EH- OR ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
24 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 > 0 (down)
25 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 <= 0 (up)
26 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, 0.7 <= theta2 <= 1.5,
27 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 > 0 (down)
28 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 <= 0 (up)
29 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, 0.7 <= theta2 <= 1.5,
30
31 FindScatters TwoTracksScatter
32 FindScatters ExELScatter points to stop, line needed to be reversed \ No newline at end of file

Revision 12009-04-28 - PeterKammel

Line: 1 to 1
Added:
>
>
META TOPICPARENT name="MuCapPKSomeNotes"
-- PeterKammel - 28 Apr 2009

Description of TMuStop Class

This class inherits from TMuEntrance, which inherits from TObject.

Since this object is intended to be written to a Tree, the space
required for the member variables was considered:
1) In some cases a desired Float_t (4 bytes) is instead stored as a
UShort_t (2 bytes), first multiplying the Float_t by 100 before
casting to a UShort_t. Functions are provided to make this
conversion invisible to the user
(ex., void SetStopAnode(Float_t anode), Float_t GetStopAnode(),
etc.).
2) TPC tracking information is more compactly stored in a "pixel
space," which corresponds to TDC400 clock ticks and TPC wire
numbers. The value fTDiffY0 is the time difference in ns
between the "YPixel=0" column and fT (the muon entrance time).
In order to allow for fractional pixel locations, in many cases
the pixel-space Float_t value is multiplied by 100 and stored as
a UShort_t.
3) fTpcPixelList is naturally in pixel-space.
4) fOneLine and fTwoLine are in pixel space.


fStopFlags (UShort_t)
Bit Description
1 "IsValid." Since initially each next-nearest-neighbor-connected
island with an EH pixel gets a TMuStop object, and some islands
may be pieced together in case of a scatter, islands that are
incorporated into another will have this flag set to zero.

2 "HasCoincStrips." There are strips EL pixels coincident with the
stop drift time.

3 "IsGoodStop." All EL points are within "GoodBox" fiducial boundary.

4 "IsBetterStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within "BetterBox" fiducial boundary. This is the
standard fiducial cut for lifetime histograms.

5 "IsGoldenStop." In addition to "GoodBox" EL cut (bit3), the stop
location is within a much tighter fiducial boundary, "GoldenBox."

6 "IsSpot." The extent of this island in Y and Z is small.
if(fAnodeMaxEL - fAnodeMaxEH < 3 && fAnodeMinEH - fAnodeMinEL < 3
&& fYPixelMaxEL - fYPixelMaxEH < 4 && fYPixelMinEH - fYPixelMinEL < 3
&& fAnodeMaxEH - fAnodeMinEH < 4 && fYPixelMaxEH - fYPixelMinEH < 5)

7 "HasEVH." At least 1 EVH pixel somewhere in track.

8 "IsNotAlone." There is another TMuStop in the same ROI as this track.

9 "OldIsTrack." A very minimal cut, it means there are not too many
EL pixels downstream of the stop (fAnodeMaxEL - fAnodeMaxEH
<= mutrack_parameters.head_maxEL_maxEH_diff), and there are
sufficient pixels leading up to any EH hit (fAnodeMinEH -
fAnodeMinEL >= mutrack_parameters.tail_minEH_minEL_diff).

10 "HeaddlminEHCut." There is an EH pixel sufficiently close to the
endpoint of the fit line. dl = distance from endpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_head = the minimum dl of the EH pixels wrt endpoint of line.
This flag is set to 1 if dlmin_head <= mutrack_parameters.head_dlminEH.

11 "TaildlminEHCut." There is no EH pixel too close to the startpoint of
the fit line. dl = distance from startpoint, projected onto
the line, of EH pixels that are not too far transversally
(d <= mutrack_parameters.stopMaxd) from the line.
dlmin_tail = the minimum dl of the EH pixels wrt startpoint of line.
This flag is set to 1 if dlmin_tail >= mutrack_parameters.tail_dlminEH.

12 "ChisqDofCut." The pixels fit well to the line (one or two segments).
This flag is set to 1 if
chisqDof <= mutrack_parameters.isTrackChisqDof, where chisqDof is
ReducedChisq of the TwoLineSeg if present, otherwise of the LineSeg.

13 "HasEHScatter." Two-line fit, and there is an EH pixel near scatter
point: d < mutrack_parameters.maxDistScatterC_EH, where d is the
distance in the yz-plane (pixel space) of an EH pixel to the
inflection point of a TTwoLineSeg fit.

14 "HasExtraELScatter." (set and used in MTA) ExtraEL pixels form a line
pointing at the apparent stop location of this muon.

15 "HasMWPCHitInROI." (set and used in MTA) Another track in the same
ROI as this muon overlaps the MWPC part of the drift region.

16 "HasUnreliableIsTrack." (set and used in MTA) This flag indicates some
anomaly with the IsTrack conditions, in particular a possible veto by
extra-EL pixels, such that the OldIsTrack flag should instead be used.


The class member TObject *fAuxObjPtr is intended to be optionally used as a
polymorphic object pointer in subsequent Mucap Tree Analysis (MTA) to point
to a user object, created and used in MTA, containing additional information
about the muon stop. In particular it is intended to facilitate scatter
identification. Also, there is UInt_t fAuxBits to be set and used in MTA.
fAuxObjPtr and fAuxBits are not streamed to the tree file.

UInt_t TMuStop::fAuxBits 32 flags for use in MTA

Bit Where Defined Assignment
1 ImproveTrackFit Attempted to improve 2-line fit
2 ImproveTrackFit HasUnreliableIsTrack bit set based on bad ChisqDof
3 ImproveTrackFit Relaxed HeaddlminEHCut to < 8.0
4 ImproveTrackFit HasUnreliableIsTrack bit set based on deltheta>2.9
OR (deltheta > 2.0 AND length2*sin(deltheta) < 5.)
5 ImproveTrackFit Trouble re-finding island in TpcROI (no action)
6 FindScatters This muon is leg1 (IsValid) of a track-track scatter.
7 FindScatters This muon is leg2 of a track-track scatter
8 FindScatters ExtraEL line fit used FitRejectingOutliers w/ dmax = 3.
9 ScatterHistFidCut EH-scatter, theta2 > 0 (down)
10 ScatterHistFidCut EH-scatter, theta2 <= 0 (up)
11 ScatterHistFidCut EH-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
12 ScatterHistFidCut ExtraEL-scatter, theta2 > 0 (down)
13 ScatterHistFidCut ExtraEL-scatter, theta2 <= 0 (up)
14 ScatterHistFidCut ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
15 ScatterHistFidCut MWPC-scatter, theta2 > 0 (down)
16 ScatterHistFidCut MWPC-scatter, theta2 <= 0 (up)
17 ScatterHistFidCut MWPC-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
18 ScatterHistFidCut TrackTrack-scatter, theta2 > 0 (down)
19 ScatterHistFidCut TrackTrack-scatter, theta2 <= 0 (up)
20 ScatterHistFidCut TrackTrack-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
21 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 > 0 (down)
22 ScatterHistFidCut EH- OR ExtraEL-scatter, theta2 <= 0 (up)
23 ScatterHistFidCut EH- OR ExtraEL-scatter, 0.7 <= theta2 <= 1.5, (down toward Al plate)
24 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 > 0 (down)
25 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, theta2 <= 0 (up)
26 ScatterHistFidCut EH- OR ExtraEL- OR MWPC-scatter, 0.7 <= theta2 <= 1.5,
27 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 > 0 (down)
28 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, theta2 <= 0 (up)
29 ScatterHistFidCut EH- OR ExtraEL- OR MWPC- OR TrackTrack-scatter, 0.7 <= theta2 <= 1.5,
30
31 FindScatters TwoTracksScatter
32 FindScatters ExELScatter points to stop, line needed to be reversed

 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback