Difference between revisions of "Mattione Update 09172012"

From GlueXWiki
Jump to: navigation, search
(Cherenkov Detector Details)
(LHCb RICH)
 
(24 intermediate revisions by the same user not shown)
Line 13: Line 13:
 
* Separation between Pions & Kaons as a function of momentum.  
 
* Separation between Pions & Kaons as a function of momentum.  
 
** Theta dependence is present but not shown (not super large).  
 
** Theta dependence is present but not shown (not super large).  
** NOTE: Cut NOT performed when #sigmas < 3.0!!
 
 
[[Image:Mattione_Update_09172012_GasThreshold_PiKSeparation.png|left|450px]]
 
[[Image:Mattione_Update_09172012_GasThreshold_PiKSeparation.png|left|450px]]
 
<br style="clear:both;"/>
 
<br style="clear:both;"/>
Line 27: Line 26:
 
=== LHCb RICH ===
 
=== LHCb RICH ===
 
* BIG IMPORTANT NOTE: These studies were performed with the parameters of the LHCb RICH AS IS (at a ~$4 million cost)
 
* BIG IMPORTANT NOTE: These studies were performed with the parameters of the LHCb RICH AS IS (at a ~$4 million cost)
* From 0 to 20 degrees
+
* From 0 to ~11 degrees
 
* 99.3% K/Pi separation from 0.6 to 2.6 GeV/c
 
* 99.3% K/Pi separation from 0.6 to 2.6 GeV/c
 
* Perfect K/Pi separation above 2.6 GeV/c
 
* Perfect K/Pi separation above 2.6 GeV/c
Line 35: Line 34:
 
== Studies ==
 
== Studies ==
  
 +
* [[Mattione_Kaon_PID_Studies_phi3(1850) | Kaon PID Studies: &phi;<sub>3</sub>(1850)]]
 
* [[Mattione_Kaon_PID_Studies_Y2175 | Kaon PID Studies: Y(2175)]]
 
* [[Mattione_Kaon_PID_Studies_Y2175 | Kaon PID Studies: Y(2175)]]
* [[Mattione_Kaon_PID_Studies_eta1Prime(2300) | Kaon PID Studies: &eta;'(2300)]]
+
* [[Mattione_Kaon_PID_Studies_eta1Prime(2300) | Kaon PID Studies: &eta;<sub>1</sub>'(2300)]]
* [[Mattione_Kaon_PID_Studies_h2Prime(2600) | Kaon PID Studies: h2'(2600)]]
+
* [[Mattione_Kaon_PID_Studies_h2Prime(2600) | Kaon PID Studies: h<sub>2</sub>'(2600)]]
* [[Mattione_Kaon_PID_Studies_phi3(1850) | Kaon PID Studies: &phi;3(1850)]]
+
  
 
== Summary ==
 
== Summary ==
 +
 +
=== Analysis ===
 +
* Branching fractions and formulas used for determination of signal / background ratios: [[File:Mattione_Update_09172012_SignalToBackgroundCalc.pdf]]
 +
* Only the VertexP4 KinFit results (for both missing proton and no-missing) are shown here (results for no kinfit & p4 kinfit are in the reaction-specific sections)
 +
* Cuts: 8.5 < E<sub>&gamma;</sub> GeV < 9.0, PID confidence level (TOF, DC dE/dx) > 1%, Kinematic fit confidence level > 1%, +/- 1.5 &Gamma; (~80% of signal integral) cut on resonance peaks
 +
* "K: Inclusive" CKOV-Cut: 3-sigma cut on whether the particle COULD be a kaon
 +
* "K: 1 Exclusive" CKOV-Cut: 3-sigma cut on whether at least one kaon candidate WAS NOT a pion (require that the particle was definitely not a kaon at 99.7% certainty)
 +
** Note: the loss of events due to this was a lot smaller than I thought it would be.  However, this is because we get 4&sigma; separation between pions & kaons with the TOF up to p = 2.5 GeV/c (with an 80ps &sigma;) rather than up to 2 GeV/c as shown in the p-vs-theta plots in the proposal.
 +
[[Image:Mattione_Update_09172012_TOFPiKDeltaT.png|left|450px]]
 +
<br style="clear:both;"/>
 +
* Note: Improvements in signal / background ratio (No-CKOV) by using VertexP4 KinFit instead of just P4 KinFit (with proton missing):
 +
** &phi;<sub>3</sub>(1850): 2.0x, Y(2175): 18.38x, &eta;'<sub>1</sub>(2300): 8.5x, h'<sub>2</sub>(2600): 8.3x
 +
 +
=== Ratios as Function of Cuts ===
 +
* &phi;<sub>3</sub>(1850) VertexP4 KinFit Signal / Background Ratios as Function of Analysis Actions (for pre-CKOV ONLY!!)
 +
** NOTE: This isn't a fair signal / background ratio comparison because:
 +
*** I don't cut on the &phi;<sub>3</sub>(1850) mass until near the end
 +
*** In this pass of the analysis, I generated the bggen events with the entire photon energy spectrum (the next pass I'm only generating between 8.5 & 9.0 GeV)
 +
*** I didn't put a photon energy cut until the end. 
 +
{| class="wikitable"
 +
| [[Image:Mattione_Update_09172012_SigBackRatioVsCut.png|left|550px]]
 +
| [[Image:Mattione_Update_09172012_SigBackRatioVsCut_NoMissing.png|left|550px]]
 +
|}
 +
<br style="clear:both;"/>
 +
 +
=== Yield %'s ===
 +
* % of generated events remaining after all cuts
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="100"| Detector / PID Cut
 +
!width="100"| &phi;<sub>3</sub>(1850), (p)
 +
!width="100"| Y(2175), (p)
 +
!width="100"| &eta;'<sub>1</sub>(2300), (p)
 +
!width="100"| h'<sub>2</sub>(2600), (p)
 +
!width="100"| &phi;<sub>3</sub>(1850), p
 +
!width="100"| Y(2175), p
 +
!width="100"| &eta;'<sub>1</sub>(2300), p
 +
!width="100"| h'<sub>2</sub>(2600), p
 +
|-
 +
! No CKOV
 +
| 27.1151% || 7.50864% || 2.28122% || 2.99637% || 11.6797% || 2.64085% || 0.823292% || 0.98661%
 +
|-
 +
! He DC <br style="clear:both;"/> K: 1 Exclusive
 +
| 19.7511% || 5.92137% || 1.09893% || 2.05918% || 8.37895% || 2.06825% || 0.394798% || 0.698757%
 +
|-
 +
! Gas Threshold <br style="clear:both;"/> K: 1 Exclusive
 +
| 24.3453% || 7.32316% || 1.36843% || 2.31741% || 9.99005% || 2.57775% || 0.493155% || 0.785577%
 +
|-
 +
! HBRICH <br style="clear:both;"/> K: 1 Exclusive
 +
| 26.3715% || 7.18456% || 1.89643% || 2.61845% || 11.1512% || 2.53028% || 0.684736% || 0.909531%
 +
|-
 +
! LHCb RICH <br style="clear:both;"/> K: 1 Exclusive
 +
| 26.1134% || 7.2575% || 1.50837% || 2.43328% || 10.8679% || 2.55721% || 0.542613% || 0.829916%
 +
|}
 +
<br style="clear:both;"/>
 +
 +
=== Exclusive Kaon Cut Signal Yield Loss ===
 +
* % of pre-CKOV signal events LOST due to exclusive Kaon cut
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="100"| Detector / PID Cut
 +
!width="100"| &phi;<sub>3</sub>(1850), (p)
 +
!width="100"| Y(2175), (p)
 +
!width="100"| &eta;'<sub>1</sub>(2300), (p)
 +
!width="100"| h'<sub>2</sub>(2600), (p)
 +
!width="100"| &phi;<sub>3</sub>(1850), p
 +
!width="100"| Y(2175), p
 +
!width="100"| &eta;'<sub>1</sub>(2300), p
 +
!width="100"| h'<sub>2</sub>(2600), p
 +
|-
 +
! He DC <br style="clear:both;"/> K: 1 Exclusive
 +
| 27.1582% || 21.1393% || 51.8271% || 31.2778% || 28.2606% || 21.6826% || 52.0465% || 29.176%
 +
|-
 +
! Gas Threshold <br style="clear:both;"/> K: 1 Exclusive
 +
| 10.2149% || 2.47022% || 40.013% || 22.6596% || 14.4665% || 2.38949% || 40.0996% || 20.3761%
 +
|-
 +
! HBRICH <br style="clear:both;"/> K: 1 Exclusive
 +
| 2.74221% || 4.31607% || 16.8677% || 12.6127% || 4.52537% || 4.18701% || 16.8295% || 7.8125%
 +
|-
 +
! LHCb RICH <br style="clear:both;"/> K: 1 Exclusive
 +
| 3.69439% || 3.34465% || 33.8788% || 18.7925% || 6.95065% || 3.16726% || 34.0923% || 15.882%
 +
|}
 +
<br style="clear:both;"/>
 +
 +
=== Signal / Background Ratios ===
 +
* Signal / Background Ratios
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="100"| Detector / PID Cut
 +
!width="100"| &phi;<sub>3</sub>(1850), (p)
 +
!width="100"| Y(2175), (p)
 +
!width="100"| &eta;'<sub>1</sub>(2300), (p)
 +
!width="100"| h'<sub>2</sub>(2600), (p)
 +
!width="100"| &phi;<sub>3</sub>(1850), p
 +
!width="100"| Y(2175), p
 +
!width="100"| &eta;'<sub>1</sub>(2300), p
 +
!width="100"| h'<sub>2</sub>(2600), p
 +
|-
 +
! No CKOV
 +
| 0.0195118 || 0.0250697 || 0.00289624 || 0.0047081 || 1.07163 ||  ||  ||
 +
|-
 +
! He DC <br style="clear:both;"/> K: Inclusive
 +
| 0.0331994 || 0.0365619 || 0.00362402 || 0.00650644 || 1.13469 ||  ||  ||
 +
|-
 +
! Gas Threshold <br style="clear:both;"/> K: Inclusive
 +
| 0.0433289 || 0.0343859 || 0.00408972 || 0.00744195 || 1.15175 ||  ||  ||
 +
|-
 +
! HBRICH <br style="clear:both;"/> K: Inclusive
 +
| 0.0583342 || 0.0585217 || 0.00443097 || 0.0122539 || 1.15841 ||  ||  ||
 +
|-
 +
! LHCb RICH <br style="clear:both;"/> K: Inclusive
 +
| 0.0501617 || 0.0530719 || 0.00455084 || 0.0103902 || 1.15521 ||  ||  ||
 +
|-
 +
! He DC <br style="clear:both;"/> K: 1 Exclusive
 +
| 0.0840888 || 0.0658157 || 0.0069218 || 0.0119952 || 1.23877 ||  ||  ||
 +
|-
 +
! Gas Threshold <br style="clear:both;"/> K: 1 Exclusive
 +
| 0.115124 || 0.0813565 || 0.0080035 || 0.0123712 || 1.11751 ||  ||  ||
 +
|-
 +
! HBRICH <br style="clear:both;"/> K: 1 Exclusive
 +
| 0.154364 || 0.104682 || 0.0124215 || 0.0202096 || 1.14907 ||  ||  ||
 +
|-
 +
! LHCb RICH <br style="clear:both;"/> K: 1 Exclusive
 +
| 0.137432 || 0.105741 || 0.0112272 || 0.0167062 || 1.22991 ||  ||  ||
 +
|}
 +
<br style="clear:both;"/>
 +
 +
=== Signal / Background Ratios : Improvements Over No CKOV ===
 +
* Detector Signal / Background Ratio (StBR) DIVIDED BY No CKOV StBR
 +
{| border="1" cellpadding="2" align="left" valign="top" style="text-align:center;"
 +
!width="100"| Detector / PID Cut
 +
!width="100"| &phi;<sub>3</sub>(1850), (p)
 +
!width="100"| Y(2175), (p)
 +
!width="100"| &eta;'<sub>1</sub>(2300), (p)
 +
!width="100"| h'<sub>2</sub>(2600), (p)
 +
!width="100"| &phi;<sub>3</sub>(1850), p
 +
!width="100"| Y(2175), p
 +
!width="100"| &eta;'<sub>1</sub>(2300), p
 +
!width="100"| h'<sub>2</sub>(2600), p
 +
|-
 +
! He DC <br style="clear:both;"/> K: Inclusive
 +
| 1.7015 || 1.45841 || 1.25129 || 1.38197 || 1.05884 ||  ||  ||
 +
|-
 +
! Gas Threshold <br style="clear:both;"/> K: Inclusive
 +
| 2.22065 || 1.37161 || 1.41208 || 1.58067 || 1.07477 ||  ||  ||
 +
|-
 +
! HBRICH <br style="clear:both;"/> K: Inclusive
 +
| 2.98968 || 2.33436 || 1.5299 || 2.60273 || 1.08098 ||  ||  ||
 +
|-
 +
! LHCb RICH <br style="clear:both;"/> K: Inclusive
 +
| 2.57083 || 2.11697 || 1.57129 || 2.20688 || 1.07799 ||  ||  ||
 +
|-
 +
! He DC <br style="clear:both;"/> K: 1 Exclusive
 +
| 4.30963 || 2.6253 || 2.38993 || 2.54777 || 1.15596 ||  ||  ||
 +
|-
 +
! Gas Threshold <br style="clear:both;"/> K: 1 Exclusive
 +
| 5.90022 || 3.24521 || 2.76341 || 2.62763 || 1.04281 ||  ||  ||
 +
|-
 +
! HBRICH <br style="clear:both;"/> K: 1 Exclusive
 +
| 7.91128 || 4.17563 || 4.28883 || 4.29252 || 1.07227 ||  ||  ||
 +
|-
 +
! LHCb RICH <br style="clear:both;"/> K: 1 Exclusive
 +
| 7.04351 || 4.21788 || 3.87647 || 3.54839 || 1.1477 ||  ||  ||
 +
|}
 +
<br style="clear:both;"/>
 +
 +
<!--
 
* NOTE: Yields of &phi;<sub>3</sub>(1850) (p detected) (both signal and background) are somewhat inflated (it looks like the proton track uncertainties are too large, causing a spike at 0 in the pulls and 1 in the confidence level)
 
* NOTE: Yields of &phi;<sub>3</sub>(1850) (p detected) (both signal and background) are somewhat inflated (it looks like the proton track uncertainties are too large, causing a spike at 0 in the pulls and 1 in the confidence level)
 
=== genr8 Signal Events ===
 
=== genr8 Signal Events ===
Line 108: Line 271:
 
|}
 
|}
 
<br style="clear:both;"/>
 
<br style="clear:both;"/>
 +
-->

Latest revision as of 23:40, 4 October 2012

Cherenkov Detector Details

Helium Drift Chamber

  • From 1 to 18 degrees.
  • K/proton separation not shown but is included in the study.
  • Separation between Pions & Kaons as a function of momentum.
    • Theta dependence changes above 12 degrees but is not shown (not super large).
Mattione Update 09172012 HeDC PiKSeparation.png


Gas Threshold (Original Design)

  • From ~1 or ~0 to 10 degrees
  • No K/proton separation
  • Separation between Pions & Kaons as a function of momentum.
    • Theta dependence is present but not shown (not super large).
Mattione Update 09172012 GasThreshold PiKSeparation.png


HBRICH

  • From ~1 or ~0 to 20 degrees
  • K/proton separation not shown but is included in the study.
  • Separation between Pions & Kaons as a function of momentum.
    • Theta dependence is present but not shown (not super large).
Mattione Update 09172012 RICH PiKSeparation.png


LHCb RICH

  • BIG IMPORTANT NOTE: These studies were performed with the parameters of the LHCb RICH AS IS (at a ~$4 million cost)
  • From 0 to ~11 degrees
  • 99.3% K/Pi separation from 0.6 to 2.6 GeV/c
  • Perfect K/Pi separation above 2.6 GeV/c
  • 99.3% K/proton separation from 2.0 to 10.0 GeV/c
  • Perfect K/proton separation above 10.0 GeV/c

Studies

Summary

Analysis

  • Branching fractions and formulas used for determination of signal / background ratios: File:Mattione Update 09172012 SignalToBackgroundCalc.pdf
  • Only the VertexP4 KinFit results (for both missing proton and no-missing) are shown here (results for no kinfit & p4 kinfit are in the reaction-specific sections)
  • Cuts: 8.5 < Eγ GeV < 9.0, PID confidence level (TOF, DC dE/dx) > 1%, Kinematic fit confidence level > 1%, +/- 1.5 Γ (~80% of signal integral) cut on resonance peaks
  • "K: Inclusive" CKOV-Cut: 3-sigma cut on whether the particle COULD be a kaon
  • "K: 1 Exclusive" CKOV-Cut: 3-sigma cut on whether at least one kaon candidate WAS NOT a pion (require that the particle was definitely not a kaon at 99.7% certainty)
    • Note: the loss of events due to this was a lot smaller than I thought it would be. However, this is because we get 4σ separation between pions & kaons with the TOF up to p = 2.5 GeV/c (with an 80ps σ) rather than up to 2 GeV/c as shown in the p-vs-theta plots in the proposal.
Mattione Update 09172012 TOFPiKDeltaT.png


  • Note: Improvements in signal / background ratio (No-CKOV) by using VertexP4 KinFit instead of just P4 KinFit (with proton missing):
    • φ3(1850): 2.0x, Y(2175): 18.38x, η'1(2300): 8.5x, h'2(2600): 8.3x

Ratios as Function of Cuts

  • φ3(1850) VertexP4 KinFit Signal / Background Ratios as Function of Analysis Actions (for pre-CKOV ONLY!!)
    • NOTE: This isn't a fair signal / background ratio comparison because:
      • I don't cut on the φ3(1850) mass until near the end
      • In this pass of the analysis, I generated the bggen events with the entire photon energy spectrum (the next pass I'm only generating between 8.5 & 9.0 GeV)
      • I didn't put a photon energy cut until the end.
Mattione Update 09172012 SigBackRatioVsCut.png
Mattione Update 09172012 SigBackRatioVsCut NoMissing.png


Yield %'s

  •  % of generated events remaining after all cuts
Detector / PID Cut φ3(1850), (p) Y(2175), (p) η'1(2300), (p) h'2(2600), (p) φ3(1850), p Y(2175), p η'1(2300), p h'2(2600), p
No CKOV 27.1151% 7.50864% 2.28122% 2.99637% 11.6797% 2.64085% 0.823292% 0.98661%
He DC
K: 1 Exclusive
19.7511% 5.92137% 1.09893% 2.05918% 8.37895% 2.06825% 0.394798% 0.698757%
Gas Threshold
K: 1 Exclusive
24.3453% 7.32316% 1.36843% 2.31741% 9.99005% 2.57775% 0.493155% 0.785577%
HBRICH
K: 1 Exclusive
26.3715% 7.18456% 1.89643% 2.61845% 11.1512% 2.53028% 0.684736% 0.909531%
LHCb RICH
K: 1 Exclusive
26.1134% 7.2575% 1.50837% 2.43328% 10.8679% 2.55721% 0.542613% 0.829916%


Exclusive Kaon Cut Signal Yield Loss

  •  % of pre-CKOV signal events LOST due to exclusive Kaon cut
Detector / PID Cut φ3(1850), (p) Y(2175), (p) η'1(2300), (p) h'2(2600), (p) φ3(1850), p Y(2175), p η'1(2300), p h'2(2600), p
He DC
K: 1 Exclusive
27.1582% 21.1393% 51.8271% 31.2778% 28.2606% 21.6826% 52.0465% 29.176%
Gas Threshold
K: 1 Exclusive
10.2149% 2.47022% 40.013% 22.6596% 14.4665% 2.38949% 40.0996% 20.3761%
HBRICH
K: 1 Exclusive
2.74221% 4.31607% 16.8677% 12.6127% 4.52537% 4.18701% 16.8295% 7.8125%
LHCb RICH
K: 1 Exclusive
3.69439% 3.34465% 33.8788% 18.7925% 6.95065% 3.16726% 34.0923% 15.882%


Signal / Background Ratios

  • Signal / Background Ratios
Detector / PID Cut φ3(1850), (p) Y(2175), (p) η'1(2300), (p) h'2(2600), (p) φ3(1850), p Y(2175), p η'1(2300), p h'2(2600), p
No CKOV 0.0195118 0.0250697 0.00289624 0.0047081 1.07163
He DC
K: Inclusive
0.0331994 0.0365619 0.00362402 0.00650644 1.13469
Gas Threshold
K: Inclusive
0.0433289 0.0343859 0.00408972 0.00744195 1.15175
HBRICH
K: Inclusive
0.0583342 0.0585217 0.00443097 0.0122539 1.15841
LHCb RICH
K: Inclusive
0.0501617 0.0530719 0.00455084 0.0103902 1.15521
He DC
K: 1 Exclusive
0.0840888 0.0658157 0.0069218 0.0119952 1.23877
Gas Threshold
K: 1 Exclusive
0.115124 0.0813565 0.0080035 0.0123712 1.11751
HBRICH
K: 1 Exclusive
0.154364 0.104682 0.0124215 0.0202096 1.14907
LHCb RICH
K: 1 Exclusive
0.137432 0.105741 0.0112272 0.0167062 1.22991


Signal / Background Ratios : Improvements Over No CKOV

  • Detector Signal / Background Ratio (StBR) DIVIDED BY No CKOV StBR
Detector / PID Cut φ3(1850), (p) Y(2175), (p) η'1(2300), (p) h'2(2600), (p) φ3(1850), p Y(2175), p η'1(2300), p h'2(2600), p
He DC
K: Inclusive
1.7015 1.45841 1.25129 1.38197 1.05884
Gas Threshold
K: Inclusive
2.22065 1.37161 1.41208 1.58067 1.07477
HBRICH
K: Inclusive
2.98968 2.33436 1.5299 2.60273 1.08098
LHCb RICH
K: Inclusive
2.57083 2.11697 1.57129 2.20688 1.07799
He DC
K: 1 Exclusive
4.30963 2.6253 2.38993 2.54777 1.15596
Gas Threshold
K: 1 Exclusive
5.90022 3.24521 2.76341 2.62763 1.04281
HBRICH
K: 1 Exclusive
7.91128 4.17563 4.28883 4.29252 1.07227
LHCb RICH
K: 1 Exclusive
7.04351 4.21788 3.87647 3.54839 1.1477