| 1 |
* |
| 2 |
* $Id: gpgig.F,v 1.1 2005/12/20 12:21:05 cafagna Exp $ |
| 3 |
* |
| 4 |
* $Log: gpgig.F,v $ |
| 5 |
* Revision 1.1 2005/12/20 12:21:05 cafagna |
| 6 |
* gpnd directory added along with ND files |
| 7 |
* |
| 8 |
* Revision 3.1.1.1 2002/07/11 16:02:00 cafagna |
| 9 |
* First GPAMELA release on CVS |
| 10 |
* |
| 11 |
* |
| 12 |
*CMZ : 3.00/00 11/05/2001 13.44.09 by Marialuigia Ambriola |
| 13 |
*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola |
| 14 |
*CMZ : 2.00/00 03/03/2000 15.39.05 by Francesco Cafagna |
| 15 |
*CMZ : 1.01/00 30/04/96 18.28.33 by Francesco Cafagna |
| 16 |
*CMZ : 1.00/03 29/04/96 16.56.08 by Francesco Cafagna |
| 17 |
*-- Author : Francesco Cafagna 29/04/96 |
| 18 |
SUBROUTINE GPGIG |
| 19 |
************************************************************************ |
| 20 |
* * |
| 21 |
* User action in GUSTEP * |
| 22 |
* * |
| 23 |
* Variables definition: * |
| 24 |
* IN: * |
| 25 |
* IFLAG = 1, just entered GUSTEP * |
| 26 |
* IFLAG = 2, leaving GUSTEP * |
| 27 |
* * |
| 28 |
* Called by: GUSTEP * |
| 29 |
* Author: Francesco Cafagna, 29/04/96 16.56.08 * |
| 30 |
* * |
| 31 |
************************************************************************ |
| 32 |
#include "gctmed.inc" |
| 33 |
#include "gcking.inc" |
| 34 |
#include "gctrak.inc" |
| 35 |
#include "gckine.inc" |
| 36 |
#include "gpgig.inc" |
| 37 |
#include "gpgneut.inc" |
| 38 |
* |
| 39 |
LOGICAL W,GPLOOK |
| 40 |
INTEGER IFLAG |
| 41 |
INTEGER IFL,IN,MECNAM(MAXMEC) |
| 42 |
c$$$ IN = INWVOL |
| 43 |
c$$$ IF(IN.NE.0) RETURN |
| 44 |
|
| 45 |
IF ((GETOT.GT.0.009097).AND.(GETOT.LT.0.028608) |
| 46 |
+ .AND.IPART.EQ.1) THEN |
| 47 |
* |
| 48 |
* Call the Gigantic resonance routine. |
| 49 |
* |
| 50 |
CALL GPGRES(STEP,GETOT,VECT(4)*VECT(7), |
| 51 |
+ VECT(5)*VECT(7),VECT(6)*VECT(7)) |
| 52 |
* |
| 53 |
* First of all let's stop the gamma tracking if neutron(s) is(are) generated |
| 54 |
* |
| 55 |
IF(Number_N.gt.0) THEN |
| 56 |
ISTOP = 1 |
| 57 |
DO I=1,Number_n |
| 58 |
* |
| 59 |
* Increment the produced secondaries number |
| 60 |
* |
| 61 |
NGKINE = NGKINE + 1 |
| 62 |
* |
| 63 |
* Store the mechnism number |
| 64 |
* |
| 65 |
KASE = IGMEC |
| 66 |
* |
| 67 |
* Store the secondaries position in the GCKING common |
| 68 |
* |
| 69 |
GPOS(1,NGKINE) = VECT(1) |
| 70 |
GPOS(2,NGKINE) = VECT(2) |
| 71 |
GPOS(2,NGKINE) = VECT(3) |
| 72 |
* |
| 73 |
* Now the other parameters: px, py, pz, E, Particle ID, TOF |
| 74 |
* |
| 75 |
GKIN(1,NGKINE) = SQRT(NEUT_EN(I)**2-AMGIG**2)* |
| 76 |
+ SIN(TETA_N(I))*COS(FI_N(I)) |
| 77 |
GKIN(2,NGKINE) = SQRT(NEUT_EN(I)**2-AMGIG**2)* |
| 78 |
+ SIN(TETA_N(I))*SIN(FI_N(I)) |
| 79 |
GKIN(3,NGKINE) = SQRT(NEUT_EN(I)**2-AMGIG**2)* |
| 80 |
+ COS(TETA_N(I)) |
| 81 |
GKIN(4,NGKINE) = NEUT_EN(I) |
| 82 |
GKIN(5,NGKINE) = IPGIG |
| 83 |
TOFD(NGKINE) = TOFG |
| 84 |
ENDDO |
| 85 |
ENDIF |
| 86 |
ENDIF |
| 87 |
999 CONTINUE |
| 88 |
|
| 89 |
RETURN |
| 90 |
END |
| 91 |
|
Parent Directory
| 
Revision Log