Difference between revisions of "Geant3 (How to)"

From Hall A Wiki
Jump to: navigation, search
(Created page with "<b>Back to main Moller</b> =Geant3 Simulation= The simulation package can be found here (Where should I put it?) The package needs to run on ifarm1101 ==Simple kumac with...")
 
Line 11: Line 11:
 
*>paw (workstation 3)
 
*>paw (workstation 3)
 
*>exe azz
 
*>exe azz
 +
 
*This produces an output called azz.dat  
 
*This produces an output called azz.dat  
 
===azz.dat===   
 
===azz.dat===   
Line 23: Line 24:
  
 
  -1.19  1.75  0.51  1.73  6.74    8.43      0.33    100000.  59737.  49799.  0.7530    0.448      89.621                  3.00
 
  -1.19  1.75  0.51  1.73  6.74    8.43      0.33    100000.  59737.  49799.  0.7530    0.448      89.621                  3.00
 +
 +
==How to get Levchuk correction==
 +
The simulation needs to be run three times.  One for scattering off of each shell M (free e-), L and K.
 +
*>paw (workstation 3)
 +
*>call sim_levch.f(N1,N2,N3,K1,K2,K3)
 +
Where N1,N2,N3 are the number of simulated events and K1,K2,K3 and the number of coincidence from the simulation for the M(free e-), L, and K shells respectively.

Revision as of 11:11, 20 November 2018

Back to main Moller

Geant3 Simulation

The simulation package can be found here (Where should I put it?)

The package needs to run on ifarm1101

Simple kumac with visualization

  • >./comg_inter.exe
  • >exe gmag_test

How to get analyzing power

  • >paw (workstation 3)
  • >exe azz
  • This produces an output called azz.dat

azz.dat

Q1    Q2    Q3    Q4   Dipole theta(RMS)  phi(RMS)  #sim  singles   COIN      Azz    Azz err  singles(mean left arm)  Holding Field 


-1.19  1.75  0.51  1.73   6.74    9.53       0.33    100000.  64193.  59208.  0.7460    0.411      90.210                  3.00


-1.19  1.75  0.51  1.73   6.74    9.37       0.33    100000.  63526.  58107.  0.7472    0.415      90.112                  3.00


-1.19  1.75  0.51  1.73   6.74    8.43       0.33    100000.  59737.  49799.  0.7530    0.448      89.621                  3.00

How to get Levchuk correction

The simulation needs to be run three times. One for scattering off of each shell M (free e-), L and K.

  • >paw (workstation 3)
  • >call sim_levch.f(N1,N2,N3,K1,K2,K3)

Where N1,N2,N3 are the number of simulated events and K1,K2,K3 and the number of coincidence from the simulation for the M(free e-), L, and K shells respectively.