version 1.1, 2004/04/23 18:13:18
|
version 1.2, 2004/05/04 20:57:17
|
|
|
| |
sig=sig219/1.d+06 !dsig/dtdphicm in microbarns/MeV**2/rad | sig=sig219/1.d+06 !dsig/dtdphicm in microbarns/MeV**2/rad |
| |
|
|
|
C DJG Convert to dsig/dOmega_cm using dt/d(costhetacm) = 2 qcm pcm |
|
|
|
sig = sig*2.*qstar*ppicm |
|
|
******************************************************************************* | ******************************************************************************* |
* GMH: Virtual photon to electron beam flux conversion factor | * GMH: Virtual photon to electron beam flux conversion factor |
gtpr = alpha/2./(pi**2)*vertex.e.E/vertex.Ein*(s-(targ.Mtar_struck/1000.)**2)/2./ | gtpr = alpha/2./(pi**2)*vertex.e.E/vertex.Ein*(s-(targ.Mtar_struck/1000.)**2)/2./ |
|
|
> ( efer+vertex.nu-(vertex.q+pferz*pfer)*vertex.p.E/vertex.p.P*tcos ) | > ( efer+vertex.nu-(vertex.q+pferz*pfer)*vertex.p.E/vertex.p.P*tcos ) |
> - 2*vertex.p.P*pfer | > - 2*vertex.p.P*pfer |
| |
davesig = gtpr*sig*jacobian |
c davesig = gtpr*sig*jacobian |
|
|
|
davesig = gtpr*sig |
sigma_eerho = davesig/1.d3 !ub/GeV-sr --> ub/MeV-sr | sigma_eerho = davesig/1.d3 !ub/GeV-sr --> ub/MeV-sr |
peerho = sigma_eerho | peerho = sigma_eerho |
ntup.sigcm = sig !sig_cm | ntup.sigcm = sig !sig_cm |