LMU/musrsim
LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Added position offset for positron counters.

Browse Source
This commit is contained in:
prokscha 2008-05-20 15:45:34 +00:00
parent b466f067ff
commit 93e78a2062
2 changed files with 12 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
geant3/src/lemsr/gukine.F
View File

@ -124,11 +124,12 @@ c
pz0 = p0 * sintheta * sinphi pz0 = p0 * sintheta * sinphi
c c
c quick and dirty, introduce long. polarization in -z c quick and dirty, introduce long. polarization in -z
c rotate spin in z direction by 180 degree around y axis c rotate spin in z direction by 180 degree around y axis; have to keep
c right-handed cartesian system. I could also rotate around the x-axis
c c
c pz0 = -p0 * costheta c pz0 = -p0 * costheta
c px0 = -p0 * sintheta * cosphi c px0 = -p0 * sintheta * cosphi
c py0 = p0 * sintheta * sinphi c py0 = p0 * sintheta * sinphi
c c
c quick and dirty, introduce long. polarization in +z c quick and dirty, introduce long. polarization in +z
c c

17
geant3/src/lemsr/ugeom.F
View File

@ -755,8 +755,9 @@ c--------------------------
c c
c the LEFT detector c the LEFT detector
c c
call gspos('SCIS',1,'MARS',0.,0.,0.,0,'ONLY') off = 0.
call gspos('SCOS',1,'MARS',0.,0.,0.,0,'ONLY') call gspos('SCIS',1,'MARS',0.,0.,off,0,'ONLY')
call gspos('SCOS',1,'MARS',0.,0.,off,0,'ONLY')
c c
c-------------------------- c--------------------------
c c
@ -774,8 +775,8 @@ c GSROTM(97...) defines rotation matrix 97
c c
call gsrotm(97,theta1,phi1,theta2,phi2,theta3,phi3) call gsrotm(97,theta1,phi1,theta2,phi2,theta3,phi3)
c c
call gspos('SCIS',2,'MARS',0.,0.,0.,97,'ONLY') call gspos('SCIS',2,'MARS',0.,0., off,97,'ONLY')
call gspos('SCOS',2,'MARS',0.,0.,0.,97,'ONLY') call gspos('SCOS',2,'MARS',0.,0., off, 97,'ONLY')
c c
c-------------------------------- c--------------------------------
c c
@ -815,8 +816,8 @@ c GSROTM(98...) defines rotation matrix 98
c c
call gsrotm(98,theta1,phi1,theta2,phi2,theta3,phi3) call gsrotm(98,theta1,phi1,theta2,phi2,theta3,phi3)
c c
call gspos('SCIS',3,'MARS',0.,0.,0.,98,'ONLY') call gspos('SCIS',3,'MARS',0.,0., off, 98,'ONLY')
call gspos('SCOS',3,'MARS',0.,0.,0.,98,'ONLY') call gspos('SCOS',3,'MARS',0.,0., off, 98,'ONLY')
c c
c-------------------------- c--------------------------
c c
@ -834,8 +835,8 @@ c GSROTM(99...) defines rotation matrix 99
c c
call gsrotm(99,theta1,phi1,theta2,phi2,theta3,phi3) call gsrotm(99,theta1,phi1,theta2,phi2,theta3,phi3)
c c
call gspos('SCIS',4,'MARS',0.,0.,0.,99,'ONLY') call gspos('SCIS',4,'MARS',0.,0., off, 99,'ONLY')
call gspos('SCOS',4,'MARS',0.,0.,0.,99,'ONLY') call gspos('SCOS',4,'MARS',0.,0., off, 99,'ONLY')
c c
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c c