Protect against cos(theta)>1 and/or <1

klendathu2k · genevb · commit 03484404a282 · 2021-08-23T15:21:48.000-04:00
diff --git a/asps/Simulation/geant321/gphys/grayl.F b/asps/Simulation/geant321/gphys/grayl.F
@@ -1,9 +1,13 @@
 *
-* $Id: grayl.F,v 1.1.1.1 2004/01/15 00:12:15 potekhin Exp $
+* $Id: grayl.F,v 1.2 2020/06/05 12:47:55 jwebb Exp $
 *
 * $Log: grayl.F,v $
+* Revision 1.2  2020/06/05 12:47:55  jwebb
+* Protect against cos(theta)>1 and/or <1
+*
 * Revision 1.1.1.1  2004/01/15 00:12:15  potekhin
 *
+*
 * Revision 1.1.1.1  1997/11/03 15:30:43  atlascvs
 * Importing CERNLIB version 08.21.
 *
@@ -82,6 +86,8 @@ SUBROUTINE GRAYL
 C
 C            Find out cos(theta)
       COSTH=1.-0.5*(QNEW*GEKIN1)**2
+      IF (COSTH.LT.-1.) COSTH=-1.
+      IF (COSTH.GT.1.)  COSTH=1.
       IF(RNDM(2).GT.0.5*(1.+COSTH**2)) GO TO 5
       SINTH=SQRT((1.-COSTH)*(1.+COSTH))
 C

Original file line number	Diff line number	Diff line change
`@@ -1,9 +1,13 @@`
`1`	`1`	`*`
`2`		`-* $Id: grayl.F,v 1.1.1.1 2004/01/15 00:12:15 potekhin Exp $`
	`2`	`+* $Id: grayl.F,v 1.2 2020/06/05 12:47:55 jwebb Exp $`
`3`	`3`	`*`
`4`	`4`	`* $Log: grayl.F,v $`
	`5`	`+* Revision 1.2 2020/06/05 12:47:55 jwebb`
	`6`	`+* Protect against cos(theta)>1 and/or <1`
	`7`	`+*`
`5`	`8`	`* Revision 1.1.1.1 2004/01/15 00:12:15 potekhin`
`6`	`9`	`*`
	`10`	`+*`
`7`	`11`	`* Revision 1.1.1.1 1997/11/03 15:30:43 atlascvs`
`8`	`12`	`* Importing CERNLIB version 08.21.`
`9`	`13`	`*`
`@@ -82,6 +86,8 @@ SUBROUTINE GRAYL`
`82`	`86`	`C`
`83`	`87`	`C Find out cos(theta)`
`84`	`88`	`COSTH=1.-0.5(QNEWGEKIN1)**2`
	`89`	`+ IF (COSTH.LT.-1.) COSTH=-1.`
	`90`	`+ IF (COSTH.GT.1.) COSTH=1.`
`85`	`91`	`IF(RNDM(2).GT.0.5(1.+COSTH*2)) GO TO 5`
`86`	`92`	`SINTH=SQRT((1.-COSTH)*(1.+COSTH))`
`87`	`93`	`C`