Skip to content
Snippets Groups Projects
Commit 91f31e40 authored by Sylvester Joosten's avatar Sylvester Joosten
Browse files

minor update, added debug info

parent 3e99a8d6
Branches
No related tags found
No related merge requests found
...@@ -5,8 +5,8 @@ ...@@ -5,8 +5,8 @@
################################################################################ ################################################################################
cmake_minimum_required (VERSION 3.5) cmake_minimum_required (VERSION 3.5)
project (simc) project (simc)
set (SIMC_MAJOR_VERSION 1) set (SIMC_MAJOR_VERSION 0)
set (SIMC_MINOR_VERSION 0) set (SIMC_MINOR_VERSION 1)
set (SIMC_PATCH_VERSION 0) set (SIMC_PATCH_VERSION 0)
set (SIMC_VERSION set (SIMC_VERSION
${SIMC_MAJOR_VERSION}.${SIMC_MINOR_VERSION}.${SIMC_PATCH_VERSION}) ${SIMC_MAJOR_VERSION}.${SIMC_MINOR_VERSION}.${SIMC_PATCH_VERSION})
......
...@@ -308,6 +308,11 @@ C modified 5/15/06 for poinct ...@@ -308,6 +308,11 @@ C modified 5/15/06 for poinct
vertex%p%delta = 100.*(vertex%p%P-spec%p%P)/spec%p%P vertex%p%delta = 100.*(vertex%p%P-spec%p%P)/spec%p%P
vertex%e%E = sqrt(vertex%e%P*vertex%e%P + 0.511*0.511) vertex%e%E = sqrt(vertex%e%P*vertex%e%P + 0.511*0.511)
vertex%p%E = sqrt(vertex%p%P*vertex%p%P + 0.511*0.511) vertex%p%E = sqrt(vertex%p%P*vertex%p%P + 0.511*0.511)
if(debug(5)) then
write(*,*) ' E and Delta: '
write(*,*) ' HMS: ',vertex%e%E,vertex%e%delta
write(*,*) ' SHMS: ',vertex%p%E,vertex%p%delta
endif !debug
endif endif
! Calculate the electron and proton PHYSICS angles from the spectrometer angles. ! Calculate the electron and proton PHYSICS angles from the spectrometer angles.
...@@ -318,6 +323,11 @@ C modified 5/15/06 for poinct ...@@ -318,6 +323,11 @@ C modified 5/15/06 for poinct
& vertex%e%xptar,vertex%e%yptar,vertex%e%theta,vertex%e%phi) & vertex%e%xptar,vertex%e%yptar,vertex%e%theta,vertex%e%phi)
call physics_angles(spec%p%theta,spec%p%phi, call physics_angles(spec%p%theta,spec%p%phi,
& vertex%p%xptar,vertex%p%yptar,vertex%p%theta,vertex%p%phi) & vertex%p%xptar,vertex%p%yptar,vertex%p%theta,vertex%p%phi)
if(debug(5)) then
write(*,*) ' PHYSICS ANGLES: '
write(*,*) ' HMS(deg): ',vertex%e%theta/pi*180,vertex%e%phi/pi*180
write(*,*) ' SHMS(deg): ',vertex%p%theta/pi*180,vertex%p%phi/pi*180
endif !debug
! Generate Fermi Momentum and Em for A(e,e'pi) and A(e,e'K). ! Generate Fermi Momentum and Em for A(e,e'pi) and A(e,e'K).
pfer=0.0 pfer=0.0
......
...@@ -13,6 +13,7 @@ c dydzp : yptar for proton ...@@ -13,6 +13,7 @@ c dydzp : yptar for proton
c p_mom : proton momentum ( MeV) c p_mom : proton momentum ( MeV)
c c
implicit none implicit none
include 'simulate.inc'
c c
real*8 e_mom real*8 e_mom
real*8 p_mom real*8 p_mom
...@@ -43,13 +44,21 @@ c ...@@ -43,13 +44,21 @@ c
read(51,*,end=999,err=999) dxdz,dydz,e_mom,dxdzp,dydzp,p_mom,w read(51,*,end=999,err=999) dxdz,dydz,e_mom,dxdzp,dydzp,p_mom,w
count = count + 1 count = count + 1
c c
p_mom = p_mom *1. p_mom = p_mom*1000.
e_mom = e_mom *1. e_mom = e_mom*1000.
dxdz = dxdz/1000. dxdz = dxdz/1.
dydz = dydz/1000. dydz = dydz/1.
dxdzp = dxdzp/1000.
dydzp = dydzp/1000. dxdzp = dxdzp/1.
dydzp = dydzp/1.
weight = w weight = w
if(debug(5)) then
write(*,*) ' '
write(*,*) ' NEW EVENT: ',count
write(*,*) ' HMS: ',dxdz,dydz,e_mom
write(*,*) ' SHMS: ',dxdzp,dydzp,p_mom
endif !debug
c c
return return
c c
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment