diff options
Diffstat (limited to 'Donjon/src/DETPLAT.f')
| -rw-r--r-- | Donjon/src/DETPLAT.f | 95 |
1 files changed, 95 insertions, 0 deletions
diff --git a/Donjon/src/DETPLAT.f b/Donjon/src/DETPLAT.f new file mode 100644 index 0000000..cca4764 --- /dev/null +++ b/Donjon/src/DETPLAT.f @@ -0,0 +1,95 @@ +*DECK DETPLAT + SUBROUTINE DETPLAT(DT,RESPON,REF,KC,PDD,REGUL,PDF,PDT,NDDG, + + AP,BP) +* +*----------------------------------------------------------------------- +* +*Purpose: +* Computes platinum detector reading with delay time consideration +* and normalized to a reference value. +* +*Copyright: +* Copyright (C) 2010 Ecole Polytechnique de Montreal. +* +*Author(s): +* E. Varin +* +*Parameters: +* DT +* RESPON +* REF +* KC +* PDD +* REGUL +* PDF +* PDT +* NDDG +* AP +* BP +* +*----------------------------------------------------------------------- +* + INTEGER NDDG,I + REAL RESPON,REF,DT + REAL AEP,AP(NDDG),BP(NDDG),PDT(NDDG),PDF(NDDG+1),PDD(NDDG),PDO + LOGICAL REGUL +* +* FLUX AT DETECTOR SITE +* --------------------- + IF ( KC.EQ.0) THEN +* REFERENCE CALCULATION +* + PDO = RESPON + IF (REGUL) THEN + CALL DETFIL(PDO,RESPON,RESPON,25.,DT,1) + ENDIF + IF (PDO.EQ.0.0) THEN + PDO = 1. + WRITE(6,*)'===> DETECTOR HAS AN INITIAL ZERO VALUE <===' + ENDIF + REF = PDO + ENDIF +* +* NORMALIZATION TO REFERENCE VALUE +* + RESPON = RESPON/REF +* + DO 10 I = 1 , NDDG +* + AEP = DT*PDT(I) + AP(I) = 1.0 +* + IF (AEP.GT.20.0) THEN + AEP = 20.0 + ENDIF + IF (AEP.GT.1.0E-6) THEN + AP(I) = EXP ( -AEP ) + ENDIF + BP(I) = 1.0 - AP(I) +* +10 CONTINUE +* + IF (KC.EQ.0) THEN +* +* INITIALISATION AT REFERENCE CALCULATION +* --------------------------------------- +* + DO 20 I = 1 , NDDG + PDD(I) = PDF(I+1) * RESPON +20 CONTINUE + ENDIF +* +* CALCULATION OF DETECTOR RESPONSE +* -------------------------------- +* + PDO = 0.0 + PDO = PDF(1)*RESPON + DO 30 I = 1 , NDDG + PDD(I) = AP(I)*PDD(I) + BP(I)*PDF(I+1)*RESPON + PDO = PDO + PDD(I) + 30 CONTINUE +* + RESPON = PDO +* + RETURN + END |