Meeting No.5: Nitrogen Pressure Curve
I still trying to explain why the ration between the entrance- and exit-window peak changes
when we increase the gas pressure inside the cell.
1.) 
Nitrogen run vs. the production run
2.) 
Nitrogen run vs. the empty-cell run
3He-reference cell run vs. production run:
Old results (from last week):
3.) 
4.) 
======> 50% correction factor
New results (with React.z instead of Tg.y):
Took only those events, where beam current is not zero. Since beam current is not measured for each event, I used the epics
informations, which is put into raw data file every few seconds (approx. 5sec). This correction is equivalent to the
accumulated charge correction but allows me to analyze smaller portions of run.
5.) 
6.) 
======> approx. 10% correction factor
How cell pressure profile depends on the temperature profile?
According to RTD readouts, there is a temperature profile inside the target cell. How this
affects the pressure profile inside the target. (Try to simulate this using Comsol?)
-----
RTD1 (upstream) : 78.3301
RTD2 : 46.6895
RTD3 : 40.7422
RTD4 : 41.7188
RTD5 (downstream) : 52.4121
RTD6 (pump cell up) : 224.707
RTD7 (pump cell down) : 250.781
HacAi_hpchan7: Temp at Ref. Cell : -100
HacAi_hpchan8: Temp at Ref. Cell : -19.3945
HacAi_hpchan9: Temp at Ref. Cell : -20.3125
Reference Cell Gas Pressure (psi) : -277.5
Reference Cell System Vacuum (Torr) : -282
Cooling Jets Flow Rate : 0
Target position : -9.37
Oven Air Flow : 230.6
-----
Mceep simulation of 3He reference cell run:
Mceep configuration file:
---------------------------------------------------------------------------
100000 # tries
4,4,4,4,4,4 for default ranges
938.27,1,2.2 m_eject,z_eject,em_bound
2425.49,0.,0.,2277.01,14.5,0.,100,-100.0,0. kinematics
4.5,-4.5,4.5,-4.5 momentum acceptances
'R','R',60.0, 130.0, 60.0, 130.0 nominal solid angles
126.75,1.,1. luminosity,time,spec_fac
45.,2.2,2.2 for singles only
3.016,2.0,0.00121,4,3 targ: a,z,dens,targ_mod,eloss_mod
-0.2,0.2 targ: cell start/end
1.109,1.100 drift to aperture - nom. sld. ang.
0.0,0.,0.,0.,0. beam: pol, vert, disp, df, tof_win
0.2,0.2,0.,0.,0. beam: FWHM in cm,cm,mr,mr,%
0.,0.,0.,0.,0. beam: offset in cm,cm,mr,mr,%
'N',0.0,0.0 beam: raster shape, X size, Y size
'E',T,2,-90.,0.,0.,0. ELECTRON ARM
'HRS','E',T,T,T,T,T aperture tests and tgt-->fp
'HRI','E' inverse map: fp-->tgt
'P',F,0,-90.,0.,0.,0. HADRON ARM
0 # global cuts
0 # specific cuts
6 # plots
'P1D',1,356,362,200,1.,0.,0,6,'elastic_scattered_electron_momentum.top'
'P1D',1,-100,100,200,1.,0.,0,81,'elastic_scattered_electron_xfp.top'
'P1D',1,-30,30,200,1.,0.,0,66,'elastic_scattered_electron_delta.top'
'P1D',1,-1000,1000,200,1.,0.,0,24,'elastic_scattered_missing_mass.top'
'NTU',1,27,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,45,46,47,62,63,64,66,81,82,83,84,'elastic_tree.ntu'
'P2D',0,-100,100,-20,20,25,25,1.,1.,0.,0.,0,81,83,'elastic_xyFP.top'
Comments: 3He(e,e'p) with JLAB Hall A HRS
---------------------------------------------------------------------------
The way I executed the simulation:
---------------------------------------------------------------------------
miham@gaia:~/JLab/e05102/Mceep/analiza_helium$ ./do_elastic
Enter input file prefix (.inp) >
elastic
Elast (E), Bnd State (B), Contin (C), Acceptance (A)?
C
Radiate? (0=No, 1=Yes, 2=Yes & multiphoton corr.)
2
NOTE: Cross sections will be additionally
differential in the photon energy.
Cutoff energy in MeV (RET for 1 MeV default) >
S(p,E) X Sigma_CC1 (unpolarized) --------------- 100
S(p,E) X Sigma_CC2 (unpolarized) --------------- 101
S(p,E) X Sigma_CC1 (polarized - Van Orden) ----- 200
Gross/Dmitrasinovic Deuterium IA --------------- 300
S. Jeschonnek Deuterium PWBA ------------------- 350
p(e,e'p)pi0 ------------------------------------ 500
p(e,e'pi+)n ------------------------------------ 600
p(e,e'K+)Y ------------------------------------- 700
d(e,e'p)n Arenhoevel unpolarized responses ----- 804
d(e,e'p)n Laget unpolarized responses ----- 814
Enter Option>
100
(e,e'p) or (e,e'n)? (P/N)>
p
Unit Weighting: ---------------------------------- 1
Enter Option or for more options>
Deuteron: Krautschneider ------------------------- 10
Saha parametric fit to data ------------ 11
from EPC code -------------------------- 12
Bernheim data -------------------------- 13
Van Orden parameterization ------------- 14
Paris ---------------------------------- 15
Bonn ----------------------------------- 16
Argonne V18 ---------------------------- 17
CD Bonn -------------------------------- 18
Gross ---------------------------------- 19
Enter Option or for more options>
3He: Meier-Hajduk et al. ------------------- 20
Argonne + VII ------------------------- 21
Urbana + VII -------------------------- 22
Jans/Marchand data -------------------- 23
Salme --------------------------------- 24
Enter Option or for more options>
20
Compute singles and accidentals? (Y/N)>
y
Warning: Duty Factor = 0, Setting to 100%
---------------------------------------------------------------------------
Simulation results (do not agree well with the measurements):
7.) 
8.) 
9.) 
10.) 
11.) 
Measured results:
12.) 
13.) 
14.) 
15.) 
16.) 
Last Modified: 09/28/09
Nitrogen run vs. the empty-cell run 
4.) 
======> 50% correction factor
6.) 
======> approx. 10% correction factor
8.) 
9.) 
10.) 
11.) 
13.) 
14.) 
15.) 
16.) 