DeltaDoverD Progress

From New IAC Wiki
Revision as of 13:01, 21 September 2008 by Oborn (talk | contribs) (→‎SIDIS Analysis)
Jump to navigation Jump to search

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

make plots for [math]1 GeV^2 \lt  Q^2\lt  2 GeV^2[/math].


Detector Construction

   4 chambers are built.   
TGEM:
Foils have been mounted on the TGEM comparison chamber.  Both charge collectors are mounted on the TGEM test detectors.  One TGEM test detector has the Thick PC board GEM foils which need much higher HV than the regular GEM foils.  The second TGEM test chamber has 3 GEM foils from CERN.

Need a min of 32 1 Meg Ohm resistors to complete the output termination connectors. Need 32 more termination connectors made from 16 wires.

Need to grind down 16, M3 bolts for mounting the GEM foils and TGEM PCboards.

Need 24 washers for GEM foils. Check mounting of the TGEM boards. Look up spacing and HV for the TGEM boards[1], Media:01352098.pdf .

Need to etch 2 cathodes for the TGEm boards.


Qweak:

  a.) Need to do final outer footprint machining so there is no interference with the Electron profile of the other octant.
  b.)  Need to machining back of the chamber for the Charge collector
 c.) Need to machine thick frames for the cathode and maybe GEM foils.
 d.) Apply electrical insulation to HV distribution boards
 e.) Need to mount GEM foils on the Qweak chambers.

SIS3610 I/O software

Objectives: a.) The first step will be to read 16 of the I/O input channel into a CODA data file.

b.) Display the 16 input channels on a GUI. Unfortunately, only 2 of the 16 will be used to read in the GEM output. The GEM output will transfer 128 hit/no hit signals to a single I/O channel in a serial fashion. The data from one I/O channel needs to be decoded according to the data structure described in Figure 8 and 9 of the VFAT manual.

c.) The final task will be to write a multiple trigger function so the I/O can be triggered by several different interrupt trigger signals and label those trigger signals.


Tasks:

Inject a signal into the I/O board input connector and use a Read function from the ROC to determine if the signal is high or low.

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 lower cross section might be binning


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 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

Detector work

TGEM assembled and ready for testing.

Need to assemble GEM comparison detector.

Made 2 thick frames for Qweak Cathode.

Get Fe source from TSO on loan for many months.
Machine hole punch for Qweak charge collector holes
Drill holes in Qweak chamber for Charge collector mounting
 Measure sag of Qweak foils and cathode.  Try using a string stretched across the frame


SIS3610 I/O software

Tasks:

Inject a signal into the I/O board input connector and use a Read function from the ROC to determine if the signal is high or low.


The SIS module latches input when a VME read is initiated

The command below sets a low constant output level on the SIS output which is then conencted to one of the SIS input line directly.

-> s3610WriteOutput(0,12288)
value = 0 = 0x0

I initiated a read function and saw the following bits set

-> s3610ReadInput(0)        
value = 18464 = 0x4820

Now I zero the output and read call the read function

-> s3610WriteOutput(0,0)    
value = 0 = 0x0
-> s3610ReadInput(0)    
value = 2080 = 0x820


I am clearly turning bits on and off but there is some randomness to other channels. Perhaps terminationg the other channels in 100 Ohms will solve this problem.

Check


created subroutine in SIS3610 library to generate a single pulse according to an given bit pattern passed as a decimal number

void TDpulse(int id, unsigned int val) {

 if((id<0) || (s3610p[id] == NULL)) {
   logMsg("s3610WriteOutput: ERROR : SIS3610 id %d not initialized \n",id,0,0,0,0,0);
   return;
 }
 s3610p[id]->d_out = val;
 s3610p[id]->d_out = 0;
 return;

}

-> TDpulse(0,12288) value = 0 = 0x0


Now we tried to use

-> TDpulse(0,12288) value = 0 = 0x0 -> s3610ReadInput(0) value = 2080 = 0x820

2080 d = 100011000 b => channel 4,5 and 9 were high

The input did not change

I think we need to Latch it.

18464 d = 100100000100000 b

Read manual to see how to latch the input.

Go Back