Difference between revisions of "Mattione Update 10312011"
From GlueXWiki
(→1234 Segmentation BCAL Time-Smear Details) |
(→BCAL π- BCAL σt 334 vs 1234 Segmentation) |
||
(20 intermediate revisions by the same user not shown) | |||
Line 6: | Line 6: | ||
* The 1234 segmentation is not available in the HDGEANT simulation, so its BCAL shower time was smeared manually during reconstruction | * The 1234 segmentation is not available in the HDGEANT simulation, so its BCAL shower time was smeared manually during reconstruction | ||
** The BCAL hit-time smear was disabled in mc_smear | ** The BCAL hit-time smear was disabled in mc_smear | ||
− | * The shower time was smeared by: σ<sub>smear</sub> = sqrt(σ<sup>2</sup><sub> | + | * The shower time was smeared by: σ<sub>smear</sub> = sqrt(σ<sup>2</sup><sub>t_1234_Target</sub> - σ<sup>2</sup><sub>t_334_NoSmear</sub>), where σ<sub>t_1234_Target</sub> = σ<sub>t_334</sub>*Improvement_Ratio, Improvement_Ratio = σ<sub>t_1234</sub> / σ<sub>t_334</sub> |
** These time uncertainties were studied as a function of dE and θ | ** These time uncertainties were studied as a function of dE and θ | ||
− | ** The σ<sub>t_1234</sub> / σ<sub>t_334</sub> ratio was obtained from fitting the results of David Lawrence's studies of photons | + | ** The σ<sub>t_1234</sub> / σ<sub>t_334</sub> improvement ratio was obtained from fitting the results of David Lawrence's studies of photons |
** The σ<sup>2</sup><sub>t_334</sub> and σ<sup>2</sup><sub>t_334_NoSmear</sub> variances were determined separately for protons and pions, by comparing the reconstructed shower times to their true values. | ** The σ<sup>2</sup><sub>t_334</sub> and σ<sup>2</sup><sub>t_334_NoSmear</sub> variances were determined separately for protons and pions, by comparing the reconstructed shower times to their true values. | ||
=== Improvements from last week's study === | === Improvements from last week's study === | ||
* Looking at π(1600) -> b1π events with t-slope = 2.0 instead of π(2000) -> b1π events with t-slope = 5.0: Higher proton momentum in BCAL for studying separation | * Looking at π(1600) -> b1π events with t-slope = 2.0 instead of π(2000) -> b1π events with t-slope = 5.0: Higher proton momentum in BCAL for studying separation | ||
− | * The σ<sub>t_334_NoSmear</sub> resolution is no longer assumed to be zero: it is non-negligible | + | * The σ<sub>t_334_NoSmear</sub> resolution is no longer assumed to be zero: it is non-negligible for large energy deposits (presumably due to clustering?) |
* The σ<sub>t_334</sub> resolution is determined separately for protons and pions (was using the pion resolution for both previously) | * The σ<sub>t_334</sub> resolution is determined separately for protons and pions (was using the pion resolution for both previously) | ||
* Better fit to photon σ<sub>t_1234</sub> / σ<sub>t_334</sub> distributions due to extra data point at 50 degrees | * Better fit to photon σ<sub>t_1234</sub> / σ<sub>t_334</sub> distributions due to extra data point at 50 degrees | ||
* Low-θ, low-dE time resolutions better-simulated: instead of extrapolating fit, performing linear-interpolation between data-points. | * Low-θ, low-dE time resolutions better-simulated: instead of extrapolating fit, performing linear-interpolation between data-points. | ||
+ | * 2 (presumably canceling) bug fixes, final results on hold until confirmation (tomorrow possibly?). | ||
=== NOTES === | === NOTES === | ||
Line 37: | Line 38: | ||
* γ BCAL σ<sub>t_1234</sub> / σ<sub>t_334</sub> vs. dE (GeV): | * γ BCAL σ<sub>t_1234</sub> / σ<sub>t_334</sub> vs. dE (GeV): | ||
+ | ** NOTE: When applying this "Improvement Ratio," the maximum allowed value was set to 1 (shouldn't be worse for 1234 at large angles, right?) | ||
[[Image:Mattione_Update_10312011_BCALPhotonTimeResolutionRatio.gif|thumb|left|400px]] | [[Image:Mattione_Update_10312011_BCALPhotonTimeResolutionRatio.gif|thumb|left|400px]] | ||
<br style="clear:both;"/> | <br style="clear:both;"/> | ||
Line 49: | Line 51: | ||
<br style="clear:both;"/> | <br style="clear:both;"/> | ||
− | == | + | == BCAL π<sup>-</sup> & proton 334 Segmentation σ<sub>t</sub> Fits == |
* When the BCAL time-smear is enabled, the π<sup>-</sup> & proton time resolutions are similar. | * When the BCAL time-smear is enabled, the π<sup>-</sup> & proton time resolutions are similar. | ||
− | * When the BCAL time-smear is | + | * When the BCAL time-smear is disabled, the π<sup>-</sup> & proton time resolutions are similar, except at high shower energies. |
=== BCAL π<sup>-</sup> 334 Segmentation σ<sub>t</sub> Fits === | === BCAL π<sup>-</sup> 334 Segmentation σ<sub>t</sub> Fits === | ||
Line 58: | Line 60: | ||
*** Last week, the fit was extrapolated from high-θ over the low-θ range. | *** Last week, the fit was extrapolated from high-θ over the low-θ range. | ||
{| class="wikitable" | {| class="wikitable" | ||
− | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaT_Fits1.gif|thumb|left|300px]] | + | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaT_Fits1.gif|thumb|left|300px| σ<sub>t</sub> Linearly Interpolated]] |
| [[Image:Mattione_Update_10312011_BCALPiMinusSigmaT_Fits2.gif|thumb|left|300px]] | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaT_Fits2.gif|thumb|left|300px]] | ||
| [[Image:Mattione_Update_10312011_BCALPiMinusSigmaT_Fits3.gif|thumb|left|300px]] | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaT_Fits3.gif|thumb|left|300px]] | ||
Line 75: | Line 77: | ||
** NOTE: For the third and forth dE bins, the values were linearly interpolated between the points instead of fit to a function. | ** NOTE: For the third and forth dE bins, the values were linearly interpolated between the points instead of fit to a function. | ||
{| class="wikitable" | {| class="wikitable" | ||
− | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits1.gif|thumb|left|300px]] | + | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits1.gif|thumb|left|300px| Evaluated as σ<sub>t</sub> = 0]] |
− | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits2.gif|thumb|left|300px]] | + | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits2.gif|thumb|left|300px| Evaluated as σ<sub>t</sub> = 0]] |
− | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits3.gif|thumb|left|300px]] | + | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits3.gif|thumb|left|300px| σ<sub>t</sub> Linearly Interpolated]] |
|- | |- | ||
− | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits4.gif|thumb|left|300px]] | + | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits4.gif|thumb|left|300px| σ<sub>t</sub> Linearly Interpolated]] |
| [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits5.gif|thumb|left|300px]] | | [[Image:Mattione_Update_10312011_BCALPiMinusSigmaTNoTSmear_Fits5.gif|thumb|left|300px]] | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
|} | |} | ||
<br style="clear:both;"/> | <br style="clear:both;"/> | ||
Line 101: | Line 91: | ||
*** Last week, the π- fit was used, and was extrapolated from high-θ over the low-θ range. | *** Last week, the π- fit was used, and was extrapolated from high-θ over the low-θ range. | ||
{| class="wikitable" | {| class="wikitable" | ||
− | | [[Image:Mattione_Update_10312011_BCALProtonSigmaT_Fits1.gif|thumb|left|300px]] | + | | [[Image:Mattione_Update_10312011_BCALProtonSigmaT_Fits1.gif|thumb|left|300px| Linearly Interpolated]] |
| [[Image:Mattione_Update_10312011_BCALProtonSigmaT_Fits2.gif|thumb|left|300px]] | | [[Image:Mattione_Update_10312011_BCALProtonSigmaT_Fits2.gif|thumb|left|300px]] | ||
| [[Image:Mattione_Update_10312011_BCALProtonSigmaT_Fits3.gif|thumb|left|300px]] | | [[Image:Mattione_Update_10312011_BCALProtonSigmaT_Fits3.gif|thumb|left|300px]] | ||
Line 117: | Line 107: | ||
** NOTE: For the second dE bin, for θ <= 27.0, the values were linearly interpolated between the points. For θ > 27.0, the time resolution was set to zero. | ** NOTE: For the second dE bin, for θ <= 27.0, the values were linearly interpolated between the points. For θ > 27.0, the time resolution was set to zero. | ||
{| class="wikitable" | {| class="wikitable" | ||
− | | [[Image:Mattione_Update_10312011_BCALProtonSigmaTNoTSmear_Fits1.gif|thumb|left|300px]] | + | | [[Image:Mattione_Update_10312011_BCALProtonSigmaTNoTSmear_Fits1.gif|thumb|left|300px| Linearly Interpolated]] |
− | | [[Image:Mattione_Update_10312011_BCALProtonSigmaTNoTSmear_Fits2.gif|thumb|left|300px]] | + | | [[Image:Mattione_Update_10312011_BCALProtonSigmaTNoTSmear_Fits2.gif|thumb|left|300px| Linearly Interpolated]] |
| [[Image:Mattione_Update_10312011_BCALProtonSigmaTNoTSmear_Fits3.gif|thumb|left|300px]] | | [[Image:Mattione_Update_10312011_BCALProtonSigmaTNoTSmear_Fits3.gif|thumb|left|300px]] | ||
|- | |- | ||
Line 126: | Line 116: | ||
<br style="clear:both;"/> | <br style="clear:both;"/> | ||
+ | <!-- | ||
=== BCAL Proton BCAL σ<sub>t</sub> 334 vs 1234 Segmentation === | === BCAL Proton BCAL σ<sub>t</sub> 334 vs 1234 Segmentation === | ||
* BCAL proton σ<sub>t</sub> vs. θ in dE bins, 334 vs 1234 Segmentation: | * BCAL proton σ<sub>t</sub> vs. θ in dE bins, 334 vs 1234 Segmentation: | ||
Line 137: | Line 128: | ||
|} | |} | ||
<br style="clear:both;"/> | <br style="clear:both;"/> | ||
+ | --> | ||
+ | == π(2000) vs π(1600) b1pi Events, π<sup>+</sup> BCAL Δβ vs p == | ||
+ | Δβ = β<sub>p</sub> - β<sub>t</sub> = p/sqrt(p*p + m*m) - path/(c*TOF) | ||
+ | |||
+ | * b1pi Events, π<sup>+</sup> BCAL Δβ vs p, 334 Segmentation: | ||
+ | {| class="wikitable" | ||
+ | | [[Image:Mattione_Update_10312011_BCALPiPlusDeltaBetaVsP_334Segmentation.gif|thumb|left|600px|π(1600) t-slope = 2]] | ||
+ | | [[Image:Mattione_Update_10242011_BCALPiPlusDeltaBetaVsP_334Segmentation.gif|thumb|left|600px|π(2000) t-slope = 5]] | ||
+ | |} | ||
+ | <br style="clear:both;"/> | ||
+ | |||
+ | |||
+ | <!-- | ||
== π(1600) -> b1pi Events, π<sup>+</sup> BCAL Δβ vs p == | == π(1600) -> b1pi Events, π<sup>+</sup> BCAL Δβ vs p == | ||
Δβ = β<sub>p</sub> - β<sub>t</sub> = p/sqrt(p*p + m*m) - path/(c*TOF) | Δβ = β<sub>p</sub> - β<sub>t</sub> = p/sqrt(p*p + m*m) - path/(c*TOF) | ||
− | + | ||
− | + | ||
{| class="wikitable" | {| class="wikitable" | ||
| [[Image:Mattione_Update_10312011_BCALPiPlusDeltaBetaVsP_334Segmentation.gif|thumb|left|600px|334 Segmentation]] | | [[Image:Mattione_Update_10312011_BCALPiPlusDeltaBetaVsP_334Segmentation.gif|thumb|left|600px|334 Segmentation]] | ||
Line 161: | Line 164: | ||
| [[Image:Mattione_Update_10312011_BCALPiPlusDeltaBetaVsP1600MeV_1234Segmentation.gif|thumb|left|400px|p = 1.6 GeV/c, 1234 Segmentation]] | | [[Image:Mattione_Update_10312011_BCALPiPlusDeltaBetaVsP1600MeV_1234Segmentation.gif|thumb|left|400px|p = 1.6 GeV/c, 1234 Segmentation]] | ||
|} | |} | ||
+ | --> |
Latest revision as of 12:59, 7 November 2011
Contents
Procedure
- Reconstruct π(1600) -> b1π events with t-slope = 2.0 (was 5.0 for X(2000)) using both the 334 and 1234 segmentations
- Plot Δβ vs p for π+ candidates to see proton, π+ separation.
1234 Segmentation BCAL Time-Smear Details
- The 1234 segmentation is not available in the HDGEANT simulation, so its BCAL shower time was smeared manually during reconstruction
- The BCAL hit-time smear was disabled in mc_smear
- The shower time was smeared by: σsmear = sqrt(σ2t_1234_Target - σ2t_334_NoSmear), where σt_1234_Target = σt_334*Improvement_Ratio, Improvement_Ratio = σt_1234 / σt_334
- These time uncertainties were studied as a function of dE and θ
- The σt_1234 / σt_334 improvement ratio was obtained from fitting the results of David Lawrence's studies of photons
- The σ2t_334 and σ2t_334_NoSmear variances were determined separately for protons and pions, by comparing the reconstructed shower times to their true values.
Improvements from last week's study
- Looking at π(1600) -> b1π events with t-slope = 2.0 instead of π(2000) -> b1π events with t-slope = 5.0: Higher proton momentum in BCAL for studying separation
- The σt_334_NoSmear resolution is no longer assumed to be zero: it is non-negligible for large energy deposits (presumably due to clustering?)
- The σt_334 resolution is determined separately for protons and pions (was using the pion resolution for both previously)
- Better fit to photon σt_1234 / σt_334 distributions due to extra data point at 50 degrees
- Low-θ, low-dE time resolutions better-simulated: instead of extrapolating fit, performing linear-interpolation between data-points.
- 2 (presumably canceling) bug fixes, final results on hold until confirmation (tomorrow possibly?).
NOTES
- The photon σt_1234 / σt_334 is assumed to be applicable to protons and pions
David's BCAL Photon Time Resolution Studies
- γ BCAL σt vs. dE (GeV):
Time Resolution Segmentation Ratio
- γ BCAL σt_1234 / σt_334 vs. dE (GeV):
- NOTE: When applying this "Improvement Ratio," the maximum allowed value was set to 1 (shouldn't be worse for 1234 at large angles, right?)
- Fits to γ BCAL σt_1234 / σt_334 vs. θ in dE bins:
- Don't trust fits at low-θ, so linearly-interpolate between 12 and 16 degrees.
BCAL π- & proton 334 Segmentation σt Fits
- When the BCAL time-smear is enabled, the π- & proton time resolutions are similar.
- When the BCAL time-smear is disabled, the π- & proton time resolutions are similar, except at high shower energies.
BCAL π- 334 Segmentation σt Fits
- BCAL π- 334 Segmentation σt Fits vs. θ in dE bins:
- NOTE: For the first dE bin, the values were linearly interpolated between the points instead of fit to a function.
- Last week, the fit was extrapolated from high-θ over the low-θ range.
- NOTE: For the first dE bin, the values were linearly interpolated between the points instead of fit to a function.
BCAL π- 334 Segmentation σt Fits with BCAL Time Smear Disabled
- BCAL π- 334 Segmentation σt Fits with BCAL Time Smear Disabled vs. θ in dE bins:
- NOTE: These resolutions are negligible except at high dE.
- NOTE: The bad fit results were ignored.
- NOTE: For the first two dE bins, the time smear was set to zero.
- NOTE: For the third and forth dE bins, the values were linearly interpolated between the points instead of fit to a function.
BCAL Proton 334 Segmentation σt Fits
- BCAL proton 334 Segmentation σt Fits vs. θ in dE bins:
- NOTE: For the first dE bin, the values were linearly interpolated between the points instead of fit to a function.
- Last week, the π- fit was used, and was extrapolated from high-θ over the low-θ range.
- NOTE: For the first dE bin, the values were linearly interpolated between the points instead of fit to a function.
BCAL Proton 334 Segmentation σt Fits with BCAL Time Smear Disabled
- BCAL proton 334 Segmentation σt Fits vs. θ in dE bins:
- NOTE: These resolutions are negligible except at high dE.
- NOTE: The bad fit results were ignored.
- NOTE: For the first dE bin, the values were linearly interpolated between the points instead of fit to a function.
- NOTE: For the second dE bin, for θ <= 27.0, the values were linearly interpolated between the points. For θ > 27.0, the time resolution was set to zero.
π(2000) vs π(1600) b1pi Events, π+ BCAL Δβ vs p
Δβ = βp - βt = p/sqrt(p*p + m*m) - path/(c*TOF)
- b1pi Events, π+ BCAL Δβ vs p, 334 Segmentation: