1 |
* |
* |
2 |
* $Id$ |
* $Id: gpxtr.F,v 3.1.1.1 2002/07/11 16:02:01 cafagna Exp $ |
3 |
|
* |
4 |
|
* $Log: gpxtr.F,v $ |
5 |
|
* Revision 3.1.1.1 2002/07/11 16:02:01 cafagna |
6 |
|
* First GPAMELA release on CVS |
7 |
* |
* |
|
* $Log$ |
|
8 |
* |
* |
9 |
*CMZ : 3.00/00 05/02/2002 12.51.38 by Unknown |
*CMZ : 3.00/00 05/02/2002 12.51.38 by Unknown |
10 |
*-- Author : Marialuigia Ambriola 11/05/2001 |
*-- Author : Marialuigia Ambriola 11/05/2001 |
11 |
SUBROUTINE GPXTR |
SUBROUTINE GPXTR |
12 |
************************************************************************** |
************************************************************************** |
13 |
* * |
* * |
14 |
* Generates and absorbs transition radiation spectrum in TRD * |
* Generates and absorbs transition radiatixon spectrum in TRD * |
15 |
* * |
* * |
16 |
* * |
* * |
17 |
* * |
* * |
36 |
EXTERNAL XTRINTER |
EXTERNAL XTRINTER |
37 |
EXTERNAL XTRINTEG |
EXTERNAL XTRINTEG |
38 |
REAL CAM,RL,G,NF,WR,CR,SPH,CATT,W,E |
REAL CAM,RL,G,NF,WR,CR,SPH,CATT,W,E |
39 |
|
REAL GRNDM0(1) |
40 |
INTEGER INDEX,INDEXAIR,INDEXC,LAY |
INTEGER INDEX,INDEXAIR,INDEXC,LAY |
41 |
c centimeter and GeV: |
c centimeter and GeV: |
42 |
DATA RL1/7.0E-4/,RL2/204.E-4/,N/71/,w1/31.59E-9/,w2/.7E-9/ |
DATA RL1/7.0E-4/,RL2/204.E-4/,N/71/,w1/31.59E-9/,w2/.7E-9/ |
277 |
IF(NPHTR.GT.0) THEN |
IF(NPHTR.GT.0) THEN |
278 |
DO I=1,NPHTR |
DO I=1,NPHTR |
279 |
ccc 10 E=XTRINTER(RNDM(1)*SPH,EINT,ENATT,115) |
ccc 10 E=XTRINTER(RNDM(1)*SPH,EINT,ENATT,115) |
280 |
10 E=XTRINTER(RNDM(1)*SPH,EINT,ENATT,75) |
c ml 26/03/04 |
281 |
|
c 10 E=XTRINTER(RNDM(1)*SPH,EINT,ENATT,75) |
282 |
|
10 CALL GRNDM(GRNDM0,1) |
283 |
|
E=XTRINTER(GRNDM0(1)*SPH,EINT,ENATT,75) |
284 |
|
C ml 26/03/04 |
285 |
C 28/9/2001: TENTATIVO DI FINE TUNING DELLA TR: |
C 28/9/2001: TENTATIVO DI FINE TUNING DELLA TR: |
286 |
******** NEI DATI REALI NON VEDO FOTONI CON ENERGIA MAGGIORE DI 60 keV: |
******** NEI DATI REALI NON VEDO FOTONI CON ENERGIA MAGGIORE DI 60 keV: |
287 |
c fine modifica 28/9/2001 |
c fine modifica 28/9/2001 |