[BLAST_ANAWARE] kinematics calculator for deuteron breakup channel in DGen

From: zhangchi (zhangchi@general.lns.mit.edu)
Date: Thu Jun 26 2003 - 15:36:12 EDT


Hi, I added a small function into DGen to calculate center of mass
kinematic variables from lab ones.

The function name is
Kine_Calc_DEEP(double vcms[8], double elab[3], double plab[3],
                double nlab[3], double Ebeam)
elab, plab, nlab are momentun, theta, phi of electron, proton, neutron.
where plab, nlab can be ommitted. Ebeam is beam energy, can be omitted.

output is in vcms and they are from vcms[0] to vcms[7]:
theta_q, phi_q, q vector, Q2, omega, P cms, W (photon+deuteron)^2, p
missing boosted into lab.

Here is how to use it in root:

root [0] .L libDGen_Cint.so
root [1] double vcms[8]
root [2] double elab[3] = { 0.7, 30, 0 }
root [3] double plab[3] = { 0.3, 50, 174 }
root [4] Kine_Calc_DEEP(vcms, elab, plab)
root [5] for(int i=0; i<8; ++i) cout<<vcms[i]<<endl;
55.142
180
0.426532
0.15943
0.15
0.314207
3.92118
0.17474

Please let me know if you hope to add anything to the returned
information. Also please let me know if you find any mistakes in the
codes.

The code itself is in Deuteron_EEP.cc:1089-1159

Chi



This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:29 EST