Difference between revisions of "TD Ddoverd 2008"

From New IAC Wiki
Jump to navigation Jump to search
(Created page with ' =9/5/08= == SIDIS Analysis== a.) Cross-Section comparison i.) Calculate absolute cross section for <math>\cos(\theta_{pi})</math> = 0.5, 1 < Q^2 < 4 GeV^2 , W = 1.45 +/- 0.…')
 
 
(One intermediate revision by one other user not shown)
Line 1: Line 1:
 +
[[Analysis]]<br>
 +
[[Software]]<br>
  
 
=9/5/08=
 
=9/5/08=
Line 516: Line 518:
 
| h3-h2 || [[Image:phi_angle_cm_helicity_difference_for_h3_h2_negative_PT.gif|300px]] || [[Image:phi_angle_cm_helicity_difference_for_h3_h2_positive_PT.gif|300px]]  
 
| h3-h2 || [[Image:phi_angle_cm_helicity_difference_for_h3_h2_negative_PT.gif|300px]] || [[Image:phi_angle_cm_helicity_difference_for_h3_h2_positive_PT.gif|300px]]  
 
|}<br>
 
|}<br>
 +
 +
 +
 +
[[DeltaDoverD_Progress]]

Latest revision as of 04:52, 8 December 2012

Analysis
Software

9/5/08

SIDIS Analysis

a.) Cross-Section comparison

i.) Calculate absolute cross section for [math]\cos(\theta_{pi})[/math] = 0.5, 1 < Q^2 < 4 GeV^2 , W = 1.45 +/- 0.2 GeV

[math]{\sigma} = \frac{1}{ L_{int}}dN = \frac{35}{3.3 \times 10^{34}} \times cm^2 = 10.61 \times10^{-34} \times 10^{24} \times [barn] = 1.061 \times 10^{-3} [\mu barn][/math]

Luminosity Calculation

Phi angle in CM Frame vs Relative Rate cos theta 0-4 0-6 W 1-45.jpg


ii.) Plot [math]\phi_{diff}^{CM}[/math]-vs-[math]cos{\theta}_{\pi}^{LAB}[/math] when [math]cos{\theta}_{e}^{CM}[/math] = -0.3 also plot [math]\phi_{diff}^{CM}[/math]-vs- [math]cos{\theta}_{e}^{LAB}[/math]

As one can see from histograms of the pion_theta_angle_vs_phi_angle_in_CM_Frame and the electron_theta_angle_vs_phi_angle_in_CM_Frame, the pion and electron acceptance in the region of [math]\phi_{diff}^{CM}=180[/math] is nearly zero(significantly low).

electron sector pion_theta_angle_vs_phi_angle_in_CM_Frame_after_cuts(EC_inner>0.06, EC_tot/p>0.2, nphe>2.5, [math]0.9\lt M_x\lt 1.1[/math], [math]1.1\lt W\lt 1.5[/math] ) electron_theta_angle_vs_phi_angle_in_CM_Frame_after_cuts(EC_inner>0.06, EC_tot/p>0.2, nphe>2.5, [math]0.9\lt M_x\lt 1.1[/math], [math]1.1\lt W\lt 1.5[/math] )
1 Pion theta angle vs phi angle in cm frame after cuts e sector 1.gif Electron theta angle vs phi angle in cm frame after cuts e sector 1.gif
2 Pion theta angle vs phi angle in cm frame after cuts e sector 2.gif Electron theta angle vs phi angle in cm frame after cuts e sector 2.gif
3 Pion theta angle vs phi angle in cm frame after cuts e sector 3.gif Electron theta angle vs phi angle in cm frame after cuts e sector 3.gif
4 Pion theta angle vs phi angle in cm frame after cuts e sector 4.gif Electron theta angle vs phi angle in cm frame after cuts e sector 4.gif
5 Pion theta angle vs phi angle in cm frame after cuts e sector 5.gif Electron theta angle vs phi angle in cm frame after cuts e sector 5.gif
6 Pion theta angle vs phi angle in cm frame after cuts e sector 6.gif Electron theta angle vs phi angle in cm frame after cuts e sector 6.gif


Pion theta vs phi angle cm Theta pion lab 0 32 3D.gifPhi angle cm Theta pion lab 0 32 3D.gif

[math]Q^2[/math]_vs_[math]\phi_{diff}^{CM}[/math] plot shows that the [math]Q^2[/math] cut should not make much difference on [math]\phi_{diff}^{CM}[/math] plot. The cut around [math]1 GeV^2 \lt Q^2\lt 2 GeV^2[/math] should reduce the number of pions around 0 and 360 of phi angle.

Q sqrd vs phi angle in cm frame after cuts all Q.gif

9/19/08

SIDIS Analysis

a.) Cross-Section comparison

i.) Calculate absolute cross section for [math]\cos(\theta_{pi})[/math] = 0.5, 1 < Q^2 < 4 GeV^2 , W = 1.45 +/- 0.2 GeV

[math]\frac{d \sigma}{d \Omega^*_{\pi}} = \frac{1.061 \times 10^{-3} [\mu barn]}{0.88} = 0.0012 \ne 2.4[/math]?


Documentation


The five-fold differential cross section for single pion production is equal to the following:

[math]\frac{\partial^5 \sigma}{\partial E_f \partial \Omega_e \partial {\Omega_{\pi}}^*} = \frac{1}{2 \pi} \Sigma \frac{1}{L_{int} A_{cc} \epsilon_{CC} \Delta W \Delta Q^2 \Delta cos {\theta_{\pi}}^* {\phi_{\pi}}^*} \frac{d(W, Q^2)}{d(E_f, cos \theta_e)}[/math]

The Jacobian term can be given as:

[math]\frac{d(W, Q^2)}{d(E_f, cos \theta_e)} = \frac{2 M_p E_i E_f}{W}[/math]


[math]\frac{\partial^5 \sigma}{\partial E_f \partial \Omega_e \partial {\Omega_{\pi}}^*} = \Gamma_{v}\times \frac{d^2 \sigma}{d {\Omega_{\pi}}^*}[/math]


where [math]\Gamma_{v}[/math] is the virtual photon flux and can be written as

[math]\Gamma_{v} = \frac{\alpha}{2 \pi^2 Q^2} \frac{(W^2 - {M_p}^2)E_f}{2M_p E_e} \frac{1}{1-\epsilon}[/math]
[math]\epsilon = (1 + 2(1+\frac{{\nu}^2}{Q^2})tan^2 \frac{\theta_e}{2})^{-1}[/math]

The reason of having low cross section might be binning

In paper they use the following number of bins:



Variable Num. Bin Range Bin Size
W 27 1.15 - 1.7 GeV 20 MeV
[math]Q^2[/math] 7 1.1 - 5.0 [math]GeV^2[/math] variable
[math]cos {\theta_{\pi}}^*[/math] 10 -1.0 - 1.0 0.2
[math]{\phi_{\pi}}^*[/math] 24 -180. - 180 15


Number of bins in my case were 360.


Number of bins equal is ~24

Phi angle in cm frame vs cross section actual one cuts on MissingMass W cos theta 0-5 nocut onQsqr.jpg

When i applied [math]1.6\lt Q^2\lt 1.84[/math] cut, the number of entries were reduced to 152.

In paper they use liquid-hydrogen unpolarized target(we have NH3 and ND3 polarized). The polarized electron beam current is 8 nA(in our case it is about 6 nA). and luminosity is different too. It also depends on statistics.

ii.) What is smallest angular coverage of EC (8 degrees?) Determine phi region where there is no acceptance ( where should we stop plotting data)?

iii.) Asymmetry Calculation

a.) Beam Asymmetry Plot

In order to plot Beam Asymmetry we need to plot for the first time the histogram of the invariant_mass_vs_cross_section. Then determine the structure function fitting the cross section and calculating the beam asymmetry function which is given as :

[math]A_{TL^'} = \sqrt{2\epsilon (1 - \epsilon)} \sigma_{LT^'} \times sin{\phi_{\pi}}^*[/math]

for given cos(theta) of the pion in cm, Q^2 and W invariant mass.

Choose kinematics ( a single theta and phi point )to max our stats and comparison to paper Histogram the following

  1. Number of e-pi coincidence events, number of FC counts
  2. Number of e-pi coincidence events/FC counts
  3. Number of e-pi coincidence events/FC count/Pt

for groups of runs with h_e,P_t = ++, +-, -+, --

table with 4 columns of h_e, P_t and 3 rows of the above histograms

9/26/08

SIDIS Analysis

1.) make semi-inclusive spectrum:

a.) h>0 Pt>0 b.) h > 0 pt<0 c.)h<0 pt>0 d.)h<0 pt<0

2.) ad up opposite target polarization histograms

3.) subtract h> 0 and h<0


Electron beam helicity asymmetry vs phi angle cm from paper.gif


helflag = 1 =>

helflag = 2 =>

helflag = 3 =>

helflag = 4 =>


Helcode 1 phi angle cm frame opposite target polarizations are added W 1-4 costhetapionCM 0-1.gifHelcode 2 phi angle cm frame opposite target polarizations are added W 1-4 costhetapionCM 0-1.gif Helcode 3 phi angle cm frame opposite target polarizations are added W 1-4 costhetapionCM 0-1.gifHelcode 4 phi angle cm frame opposite target polarizations are added W 1-4 costhetapionCM 0-1.gif


1 4substructed helcodes phi angle cm frame.gif3 2substructed helcodes phi angle cm frame.gif
1 4sum helcodes phi angle cm frame.gif3 2sum helcodes phi angle cm frame.gif


1 4substructed helcodes phi angle cm frame xmgrace.jpg3 2substructed helcodes phi angle cm frame xmgrace.jpg

Phi_angle_CM vs Asymmetry

Compared data(there is no table for asymmetry values like cross section):

1 4substructed helcodes phi angle cm frame xmgrace and results from.jpg3 2substructed helcodes phi angle cm frame xmgrace and results from.jpg

11/21/08

0.) insert run summarry table, Pb, Pt, PB*Pt, Btorus

1.) make semi-inclusive spectrum:

a.) h>0 Pt>0

b.) h > 0 pt<0

c.)h<0 pt>0

d.)h<0 pt<0

2.) helicity difference plots for Pt>0 and Pt<0

3.) Asym plots for Pt>0 and Pt< 0

4.) unpolarized target asymmetry

Rebin asymmetry hisograms and combine the two into a total asymmetry histogram

Run Summary Table

The table below should have run ranges with 
[math]P_b*P_t \pm \Delta (P_b *P_t)[/math] according 
to the elastic asymmetry measurements stored at
http://www.jlab.org/Hall-B/secure/eg1/EG2000/eg1b_analysis_progress.htm.


http://www.jlab.org/Hall-B/secure/eg1/EG2000/josh/pbpt/

5.73 in
pos = 0.45743 +/- 0.04269
neg = -0.38816 +/- 0.04556


5.73 out
pos = 0.46604 +/- 0.03496
neg = -0.50684 +/- 0.03578

Start Run Number Beam Polarization(Pb) Target Polarization(Pt) from dstdump file Pb*Pt Beam Torus Beam Energy (MeV)
26998 0 -0.69 0.46604 -2250 5736
26996 -0.71 -0.69 0.46604 -2250 5736
26995 -0.71 -0.69 0.46604 -2250 5736
26994 -0.71 -0.73 0.46604 -2250 5736
26993 -0.71 -0.72 0.46604 -2250 5736
26992 -0.71 -0.73 0.46604 -2250 5736
26991 -0.71 -0.75 0.46604 -2250 5736
26990 -0.71 -0.76 0.46604 -2250 5736
26989 0.71 -0.67 -0.38816 2250 5736
26988 0.71 -0.68 -0.38816 2250 5736
26987 0.71 -0.67 -0.38816 2250 5736
26986 0.71 -0.68 -0.38816 2250 5736
26985 0.71 -0.68 -0.38816 2250 5736
26984 0.71 -0.68 -0.38816 2250 5736
26983 0.71 -0.68 -0.38816 2250 5736
26981 0.71 -0.68 -0.38816 2250 5736
26980 0.71 -0.7 -0.38816 2250 5736
26979 0.71 -0.71 -0.38816 2250 5736
26966 0.71 0.74 0.45743 2250 5736
26965 0.71 0.74 0.45743 2250 5736
26964 0.71 0.74 0.45743 2250 5736
26963 0.71 0.74 0.45743 2250 5736
26962 0.71 0.74 0.45743 2250 5736
26961 0.71 0.75 0.45743 2250 5736
26960 0.71 0.76 0.45743 2250 5736
26959 0.71 0.69 0.45743 2250 5736
26958 0.71 0.69 0.45743 2250 5736
26957 0.71 0.69 0.45743 2250 5736
26956 0.71 0.69 0.45743 2250 5736
26955 0.71 0.69 0.45743 2250 5736
26954 0.71 0.69 0.45743 2250 5736
26953 0.71 0.69 0.45743 2250 5736
26952 0.71 0.7 0.45743 2250 5736
26951 0.71 0.7 0.45743 2250 5736
26948 0.71 0.7 0.45743 2250 5736
26947 0.71 0.7 0.45743 2250 5736
26946 0.71 0.7 0.45743 2250 5736
26945 0.71 0.7 0.45743 2250 5736
26943 0.7 0.7 0.45743 2250 5736
26942 0.7 0.7 0.45743 2250 5736
26941 0.7 0.7 0.45743 2250 5736
26940 0.7 0.7 0.45743 2250 5736
26939 0.7 0.7 0.45743 2250 5736
26938 0.7 0.7 0.45743 2250 5736
26937 0.7 0.7 0.45743 2250 5736
26934 0.7 0.7 0.45743 2250 5736
26933 0.7 0.71 0.45743 2250 5736
26932 0.7 0.71 0.45743 2250 5736
26931 0.7 0.71 0.45743 2250 5736
26930 0.7 0.71 0.45743 2250 5736
26929 0.7 0.71 0.45743 2250 5736
26928 0.7 0.72 0.45743 2250 5736
26927 0.7 0.72 0.45743 2250 5736
26926 0.7 0.73 0.45743 2250 5736
26925 0.7 0.73 0.45743 2250 5736
27074 0.71 -0.81 -0.38816 2250 5736
27075 0.71 -0.77 -0.38816 2250 5736
27076 0.71 -0.74 -0.38816 2250 5736
27077 0.71 -0.74 -0.38816 2250 5736
27078 071 -0.75 -0.38816 2250 5736
27079 0.71 -0.72 -0.38816 2250 5736
27100 0.71 0.76 0.45743 2250 5736
27101 0.71 0.76 0.45743 2250 5736
27102 0.71 0.73 0.45743 2250 5736
27105 0.61 0.73 0.45743 2250 5736
27106 0.61 0.73 0.45743 2250 5736
27107 0.61 -0.71 -0.38816 2250 5736
27108 0.61 -0.7 -0.38816 2250 5736
27109 0.61 -0.7 -0.38816 2250 5736
27111 0.61 -0.71 -0.38816 2250 5736
27112 0.61 -0.69 -0.38816 2250 5736
27114 0.61 0 2250 5736
27116 -0.58 0.2 -0.38816 2250 5736
27124 -0.7 0.24 -0.38816 2250 5736
27125 -0.7 0.25 -0.38816 2250 5736
27126 -0.7 0.26 -0.38816 2250 5736
27127 -0.7 0.26 -0.38816 2250 5736
27128 -0.7 0.26 -0.38816 2250 5736
27129 -0.7 0.26 -0.38816 2250 5736
27131 -0.7 0.26 -0.38816 2250 5736
27132 -0.7 0.26 -0.38816 2250 5736
27133 -0.7 0.26 -0.38816 2250 5736
27134 -0.7 0.26 -0.38816 2250 5736
27135 -0.7 0.26 -0.38816 2250 5736
27136 -0.7 0.26 -0.38816 2250 5736
27137 -0.7 0.25 -0.38816 2250 5736
27138 -0.7 0.25 -0.38816 2250 5736
27139 -0.7 0.25 -0.38816 2250 5736
27141 -0.7 0.25 -0.38816 2250 5736
27142 -0.7 0.25 -0.38816 2250 5736
27143 -0.7 0.25 -0.38816 2250 5736
27145 -0.7 0.24 -0.38816 2250 5736
27146 -0.7 0.24 -0.38816 2250 5736
27148 -0.7 0.24 -0.38816 2250 5736
27149 -0.7 -0.09 0.45743 2250 5736
27150 -0.7 -0.12 0.45743 2250 5736
27152 -0.64 -0.18 0.45743 2250 5736
27153 -0.64 -0.19 0.45743 2250 5736
27154 -0.64 -0.19 0.45743 2250 5736
27155 -0.64 -0.19 0.45743 2250 5736
27156 -0.64 -0.19 0.45743 2250 5736
27157 -0.64 -0.19 0.45743 2250 5736
27160 -0.64 -0.18 0.45743 2250 5736
27161 -0.64 -0.18 0.45743 2250 5736
27162 -0.64 -0.19 0.45743 2250 5736
27163 -0.64 -0.19 0.45743 2250 5736
27165 -0.64 -0.19 0.45743 2250 5736
27166 -0.64 0.18 -0.38816 2250 5736
27167 -0.64 0.22 -0.38816 2250 5736
27168 -0.64 0.23 -0.38816 2250 5736
27169 -0.64 0.23 -0.38816 2250 5736
27170 -0.64 0.23 -0.38816 2250 5736
27171 -0.64 0.23 -0.38816 2250 5736
27172 -0.64 -0.35 0.45743 2250 5736
27175 -0.7 -0.78 0.45743 2250 5736
27176 -0.7 -0.72 0.45743 2250 5736
27177 -0.7 -0.7 0.45743 2250 5736
27178 -0.7 -0.69 0.45743 2250 5736
27179 -0.7 -0.67 0.45743 2250 5736
27180 -0.7 -0.67 0.45743 2250 5736
27181 -0.7 -0.66 0.45743 2250 5736
27182 -0.7 -0.72 0.45743 2250 5736
27183 -0.7 -0.68 0.45743 2250 5736
27184 -0.7 -0.63 0.45743 2250 5736
27186 -0.7 0.72 -0.38816 2250 5736
27187 -0.7 0.71 -0.38816 2250 5736
27188 -0.7 0.71 -0.38816 2250 5736
27189 -0.7 0.7 -0.38816 2250 5736
27190 -0.7 0.69 -0.38816 2250 5736
27191 -0.7 0.69 -0.38816 2250 5736
27192 -0.7 0.68 -0.38816 2250 5736
27193 -0.7 0.69 -0.38816 2250 5736
27194 -0.7 0.7 -0.38816 2250 5736
27195 -0.7 0.69 -0.38816 2250 5736
27113 0.61 0.68 0.45743 2250 5736


Phi_Angle_CM_Frame for Positive and Negative Target Polarization

Helcode # Negative Target Polarization Positive Target Polarization
1 h>0 PT negative & helcode 1 phi angle cm frame.gif PT positive & helcode 1 phi angle cm frame.gif
2 h<0 PT negative & helcode 2 phi angle cm frame.gif PT positive & helcode 2 phi angle cm frame.gif
3 PT negative & helcode 3 phi angle cm frame.gif PT positive & helcode 3 phi angle cm frame.gif
4 PT negative & helcode 4 phi angle cm frame.gif PT positive & helcode 4 phi angle cm frame.gif


Helicity Difference Plots for Pt>0 and Pt<0

Helicity Difference Negative Target Polarization Positive Target Polarization
h1-h4 Phi angle cm helicity difference for h1 h4 negative PT.gif Phi angle cm helicity difference for h1 h4 positive PT.gif
h3-h2 Phi angle cm helicity difference for h3 h2 negative PT.gif Phi angle cm helicity difference for h3 h2 positive PT.gif



DeltaDoverD_Progress