MWDC track reconstruction
I have followed Ole's instruction and changed some of the parameters. I was
able to improve the track reconstruction efficiency from ~50% to ~75% (with respect to hits in
scintillators). The following
two animations show events without reconstructed tracks before and after
my modifications. With new parameters you can hardly find an event for which
a visual track could be imagined.
1.) 
2.) 
Here are the events with reconstructed tracks using new parameters:
3.) 
Here is a list of the changes that I made and corresponding analysis statistics.
*********** Before any change ********************
------------------------------------------------------------------------------
BLOCK: Decode
TriggerT3 (DL.evtypebits&8)==8 10000 292 (2.92%)
CoincT5 (DL.evtypebits&32)==32 10000 3597 (36%)
CoincT6 (DL.evtypebits&64)==64 10000 2450 (24.5%)
T1 (DL.evtypebits&2)==2 10000 3180 (31.8%)
Decode_master T1 10000 3180 (31.8%)
BLOCK: Tracking
uNRoads BB.mwdc.u.nroads>0 3180 2278 (71.6%)
uNgood BB.mwdc.u.ngood>0 3180 2042 (64.2%)
vNRoads BB.mwdc.v.nroads>0 3180 2267 (71.3%)
vNgood BB.mwdc.v.ngood>0 3180 2026 (63.7%)
xNRoads BB.mwdc.x.nroads>0 3180 2235 (70.3%)
xNgood BB.mwdc.x.ngood>0 3180 2032 (63.9%)
BBTrack BB.tr.n==0 3180 1747 (54.9%)
Tracking_master BBTrack 3180 1747 (54.9%)
BLOCK: Physics
HRSLGoldDp L.gold.dp>-0.045&&L.gold.dp<0.045 1747 68 (3.89%)
HRSLGoldTh L.gold.th>-0.06&&L.gold.th<0.06 1747 62 (3.55%)
HRSLGoldPh L.gold.ph>-0.03&&L.gold.ph<0.03 1747 63 (3.61%)
BBTpTrackMatch BB.tp.trHitIndex>-1 1747 1747 (100%)
TargetSizeCut ReactPt_L.z>-0.30&&ReactPt_L.z<0.30 1747 64 (3.66%)
*********** Change chi2_conflevel changed from 1e-5 --> 1e-3 ********************
------------------------------------------------------------------------------
BLOCK: Decode
TriggerT3 (DL.evtypebits&8)==8 10000 292 (2.92%)
CoincT5 (DL.evtypebits&32)==32 10000 3597 (36%)
CoincT6 (DL.evtypebits&64)==64 10000 2450 (24.5%)
T1 (DL.evtypebits&2)==2 10000 3180 (31.8%)
Decode_master T1 10000 3180 (31.8%)
BLOCK: Tracking
uNRoads BB.mwdc.u.nroads>0 3180 2278 (71.6%)
uNgood BB.mwdc.u.ngood>0 3180 1918 (60.3%)
vNRoads BB.mwdc.v.nroads>0 3180 2267 (71.3%)
vNgood BB.mwdc.v.ngood>0 3180 1906 (59.9%)
xNRoads BB.mwdc.x.nroads>0 3180 2235 (70.3%)
xNgood BB.mwdc.x.ngood>0 3180 1941 (61%)
BBTrack BB.tr.n==0 3180 1901 (59.8%)
Tracking_master BBTrack 3180 1901 (59.8%)
BLOCK: Physics
HRSLGoldDp L.gold.dp>-0.045&&L.gold.dp<0.045 1901 75 (3.95%)
HRSLGoldTh L.gold.th>-0.06&&L.gold.th<0.06 1901 69 (3.63%)
HRSLGoldPh L.gold.ph>-0.03&&L.gold.ph<0.03 1901 69 (3.63%)
BBTpTrackMatch BB.tp.trHitIndex>-1 1901 1901 (100%)
TargetSizeCut ReactPt_L.z>-0.30&&ReactPt_L.z<0.30 1901 70 (3.68%)
*********** Change chi2_conflevel 1e-8 ********************
------------------------------------------------------------------------------
BLOCK: Decode
TriggerT3 (DL.evtypebits&8)==8 10000 292 (2.92%)
CoincT5 (DL.evtypebits&32)==32 10000 3597 (36%)
CoincT6 (DL.evtypebits&64)==64 10000 2450 (24.5%)
T1 (DL.evtypebits&2)==2 10000 3180 (31.8%)
Decode_master T1 10000 3180 (31.8%)
BLOCK: Tracking
uNRoads BB.mwdc.u.nroads>0 3180 2278 (71.6%)
uNgood BB.mwdc.u.ngood>0 3180 2123 (66.8%)
vNRoads BB.mwdc.v.nroads>0 3180 2267 (71.3%)
vNgood BB.mwdc.v.ngood>0 3180 2096 (65.9%)
xNRoads BB.mwdc.x.nroads>0 3180 2235 (70.3%)
xNgood BB.mwdc.x.ngood>0 3180 2086 (65.6%)
BBTrack BB.tr.n==0 3180 1683 (52.9%)
Tracking_master BBTrack 3180 1683 (52.9%)
BLOCK: Physics
HRSLGoldDp L.gold.dp>-0.045&&L.gold.dp<0.045 1683 65 (3.86%)
HRSLGoldTh L.gold.th>-0.06&&L.gold.th<0.06 1683 59 (3.51%)
HRSLGoldPh L.gold.ph>-0.03&&L.gold.ph<0.03 1683 60 (3.57%)
BBTpTrackMatch BB.tp.trHitIndex>-1 1683 1683 (100%)
TargetSizeCut ReactPt_L.z>-0.30&&ReactPt_L.z<0.30 1683 61 (3.62%)
*********** Change chi2_conflevel 1e-3 && BB.mwdc.xp.res changed from 2.e-04
to 2.e-03 ********************
------------------------------------------------------------------------------
BLOCK: Decode
TriggerT3 (DL.evtypebits&8)==8 10000 292 (2.92%)
CoincT5 (DL.evtypebits&32)==32 10000 3597 (36%)
CoincT6 (DL.evtypebits&64)==64 10000 2450 (24.5%)
T1 (DL.evtypebits&2)==2 10000 3180 (31.8%)
Decode_master T1 10000 3180 (31.8%)
BLOCK: Tracking
uNRoads BB.mwdc.u.nroads>0 3180 2299 (72.3%)
uNgood BB.mwdc.u.ngood>0 3180 2290 (72%)
vNRoads BB.mwdc.v.nroads>0 3180 2288 (71.9%)
vNgood BB.mwdc.v.ngood>0 3180 2272 (71.4%)
xNRoads BB.mwdc.x.nroads>0 3180 2254 (70.9%)
xNgood BB.mwdc.x.ngood>0 3180 2240 (70.4%)
BBTrack BB.tr.n==0 3180 1229 (38.6%)
Tracking_master BBTrack 3180 1229 (38.6%)
BLOCK: Physics
HRSLGoldDp L.gold.dp>-0.045&&L.gold.dp<0.045 1229 46 (3.74%)
HRSLGoldTh L.gold.th>-0.06&&L.gold.th<0.06 1229 41 (3.34%)
HRSLGoldPh L.gold.ph>-0.03&&L.gold.ph<0.03 1229 43 (3.5%)
BBTpTrackMatch BB.tp.trHitIndex>-1 1229 1229 (100%)
TargetSizeCut ReactPt_L.z>-0.30&&ReactPt_L.z<0.30 1229 43 (3.5%)
*********** Changed 3d_chi2_conflevel 1e-9 --> 1e-4 ********************
------------------------------------------------------------------------------
BLOCK: Decode
TriggerT3 (DL.evtypebits&8)==8 10000 292 (2.92%)
CoincT5 (DL.evtypebits&32)==32 10000 3597 (36%)
CoincT6 (DL.evtypebits&64)==64 10000 2450 (24.5%)
T1 (DL.evtypebits&2)==2 10000 3180 (31.8%)
Decode_master T1 10000 3180 (31.8%)
BLOCK: Tracking
uNRoads BB.mwdc.u.nroads>0 3180 2299 (72.3%)
uNgood BB.mwdc.u.ngood>0 3180 2290 (72%)
vNRoads BB.mwdc.v.nroads>0 3180 2288 (71.9%)
vNgood BB.mwdc.v.ngood>0 3180 2272 (71.4%)
xNRoads BB.mwdc.x.nroads>0 3180 2254 (70.9%)
xNgood BB.mwdc.x.ngood>0 3180 2240 (70.4%)
BBTrack BB.tr.n==0 3180 2363 (74.3%)
Tracking_master BBTrack 3180 2363 (74.3%)
BLOCK: Physics
HRSLGoldDp L.gold.dp>-0.045&&L.gold.dp<0.045 2363 93 (3.94%)
HRSLGoldTh L.gold.th>-0.06&&L.gold.th<0.06 2363 91 (3.85%)
HRSLGoldPh L.gold.ph>-0.03&&L.gold.ph<0.03 2363 88 (3.72%)
BBTpTrackMatch BB.tp.trHitIndex>-1 2363 2363 (100%)
TargetSizeCut ReactPt_L.z>-0.30&&ReactPt_L.z<0.30 2363 87 (3.68%)
*********** 3d_chi2_conflevel 1e-9, changed 3d_matchcut from 2e-3 --> 1e-4
********************
------------------------------------------------------------------------------
BLOCK: Decode
TriggerT3 (DL.evtypebits&8)==8 10000 292 (2.92%)
CoincT5 (DL.evtypebits&32)==32 10000 3597 (36%)
CoincT6 (DL.evtypebits&64)==64 10000 2450 (24.5%)
T1 (DL.evtypebits&2)==2 10000 3180 (31.8%)
Decode_master T1 10000 3180 (31.8%)
BLOCK: Tracking
uNRoads BB.mwdc.u.nroads>0 3180 2299 (72.3%)
uNgood BB.mwdc.u.ngood>0 3180 2290 (72%)
vNRoads BB.mwdc.v.nroads>0 3180 2288 (71.9%)
vNgood BB.mwdc.v.ngood>0 3180 2272 (71.4%)
xNRoads BB.mwdc.x.nroads>0 3180 2254 (70.9%)
xNgood BB.mwdc.x.ngood>0 3180 2240 (70.4%)
BBTrack BB.tr.n==0 3180 1324 (41.6%)
Tracking_master BBTrack 3180 1324 (41.6%)
BLOCK: Physics
HRSLGoldDp L.gold.dp>-0.045&&L.gold.dp<0.045 1324 53 (4%)
HRSLGoldTh L.gold.th>-0.06&&L.gold.th<0.06 1324 50 (3.78%)
HRSLGoldPh L.gold.ph>-0.03&&L.gold.ph<0.03 1324 50 (3.78%)
BBTpTrackMatch BB.tp.trHitIndex>-1 1324 1324 (100%)
TargetSizeCut ReactPt_L.z>-0.30&&ReactPt_L.z<0.30 1324 50 (3.78%)
*********** Change chi2_conflevel 1e-3 && BB.mwdc.xp.res 2.e-03
&& BB.mwdc.maxmiss = 1********************
------------------------------------------------------------------------------
BLOCK: Decode
TriggerT3 (DL.evtypebits&8)==8 10000 292 (2.92%)
CoincT5 (DL.evtypebits&32)==32 10000 3597 (36%)
CoincT6 (DL.evtypebits&64)==64 10000 2450 (24.5%)
T1 (DL.evtypebits&2)==2 10000 3180 (31.8%)
Decode_master T1 10000 3180 (31.8%)
BLOCK: Tracking
uNRoads BB.mwdc.u.nroads>0 3180 2566 (80.7%)
uNgood BB.mwdc.u.ngood>0 3180 2557 (80.4%)
vNRoads BB.mwdc.v.nroads>0 3180 2538 (79.8%)
vNgood BB.mwdc.v.ngood>0 3180 2528 (79.5%)
xNRoads BB.mwdc.x.nroads>0 3180 2512 (79%)
xNgood BB.mwdc.x.ngood>0 3180 2499 (78.6%)
BBTrack BB.tr.n==0 3180 867 (27.3%)
Tracking_master BBTrack 3180 867 (27.3%)
BLOCK: Physics
HRSLGoldDp L.gold.dp>-0.045&&L.gold.dp<0.045 867 31 (3.58%)
HRSLGoldTh L.gold.th>-0.06&&L.gold.th<0.06 867 27 (3.11%)
HRSLGoldPh L.gold.ph>-0.03&&L.gold.ph<0.03 867 29 (3.34%)
BBTpTrackMatch BB.tp.trHitIndex>-1 867 867 (100%)
TargetSizeCut ReactPt_L.z>-0.30&&ReactPt_L.z<0.30 867 28 (3.23%)
*********** Change chi2_conflevel changed from 1e-3 to 5e-3 && BB.mwdc.xp.res 2.e-03
&& BB.mwdc.maxmiss = 1********************
--------------------------------------------------------------------------------------------------------------------------------------------------------
BLOCK: Decode
TriggerT3 (DL.evtypebits&8==8 10000 292 (2.92%)
CoincT5 (DL.evtypebits&32)==32 10000 3597 (36%)
CoincT6 (DL.evtypebits&64)==64 10000 2450 (24.5%)
T1 (DL.evtypebits&2)==2 10000 3180 (31.8%)
Decode_master T1 10000 3180 (31.8%)
BLOCK: Tracking
u1Hits BB.mwdc.u1.nhits>0 3180 2628 (82.6%)
u1pHits BB.mwdc.u1p.nhits>0 3180 2613 (82.2%)
v1Hits BB.mwdc.v1.nhits>0 3180 2596 (81.6%)
v1pHits BB.mwdc.v1p.nhits>0 3180 2627 (82.6%)
x1Hits BB.mwdc.x1.nhits>0 3180 2591 (81.5%)
x1pHits BB.mwdc.x1p.nhits>0 3180 2608 (82%)
u2Hits BB.mwdc.u2.nhits>0 3180 2746 (86.4%)
u2pHits BB.mwdc.u2p.nhits>0 3180 2791 (87.8%)
v2Hits BB.mwdc.v2.nhits>0 3180 2815 (88.5%)
v2pHits BB.mwdc.v2p.nhits>0 3180 2795 (87.9%)
x2Hits BB.mwdc.x2.nhits>0 3180 2789 (87.7%)
x2pHits BB.mwdc.x2p.nhits>0 3180 2796 (87.9%)
uProj (u1Hits&&u1pHits&&u2Hits) ||
(u1Hits&&u1pHits&&u2pHits)||
(u1pHits&&u2Hits&&u2pHits)||
(u1Hits&&u2Hits&&u2pHits) 3180 2564 (80.6%)
vProj (v1Hits&&v1pHits&&v2Hits) ||
(v1Hits&&v1pHits&&v2pHits)||
(v1pHits&&v2Hits&&v2pHits)||
(v1Hits&&v2Hits&&v2pHits) 3180 2584 (81.3%)
xProj (x1Hits&&x1pHits&&x2Hits) ||
(x1Hits&&x1pHits&&x2pHits)||
(x1pHits&&x2Hits&&x2pHits)||
(x1Hits&&x2Hits&&x2pHits) 3180 2570 (80.8%)
ThreeProjections uProj&&vProj&&xProj 3180 2429 (76.4%)
uNRoads BB.mwdc.u.nroads>0 3180 2566 (80.7%)
uNgood BB.mwdc.u.ngood>0 3180 2557 (80.4%)
vNRoads BB.mwdc.v.nroads>0 3180 2538 (79.8%)
vNgood BB.mwdc.v.ngood>0 3180 2527 (79.5%)
xNRoads BB.mwdc.x.nroads>0 3180 2512 (79%)
xNgood BB.mwdc.x.ngood>0 3180 2496 (78.5%)
YesBBTrack BB.tr.n>0 3180 2337 (73.5%)
BBTrack BB.tr.n==0 3180 843 (26.5%)
Tracking_master BBTrack 3180 843 (26.5%)
BLOCK: Physics
HRSLGoldDp L.gold.dp>-0.045&&L.gold.dp<0.045 843 29 (3.44%)
HRSLGoldTh L.gold.th>-0.06&&L.gold.th<0.06 843 25 (2.97%)
HRSLGoldPh L.gold.ph>-0.03&&L.gold.ph<0.03 843 27 (3.2%)
BBTpTrackMatch BB.tp.trHitIndex>-1 843 843 (100%)
TargetSizeCut ReactPt_L.z>-0.30&&ReactPt_L.z<0.30 843 26 (3.08%)
Here I need to emphasize, that we can not expect that all hits in scintillators
would have a track in chambers. There are many events, where only one chamber/wire-plane
was hit. For these events track reconstruction is impossible and they should
not be added to the bad statistics of the MWDCs. Therefore, I selected only
such events, for which I believe track reconstruction in the chambers should
be possible and tested how many of them have tracks. Using this procedure
I estimated, that the new track efficiency is approx. 96%.
However, we need to have in mind, that I was changing these parameters very
generously. This means that some bogus tracks could appear. This still needs to
be tested.
Last modified: 10/03/12
