*********************************************************************** c called by guinit c c adapted to EMC by K.Lapidus c SUBROUTINE EMCINIT(PARAMFILE) *********************************************************************** IMPLICIT NONE #include "geantdef.h" #include "emctuple.inc" CHARACTER*200 PARAMFILE INTEGER GETMEDNUM CHARACTER*15 EMCMED1 DATA EMCMED1 /'LEAD_GLASS$'/ c initialize PROFAST CALL PROFASTINIT(PARAMFILE, POINT1, POINT2) c PMT spectral efficiency EMCPHOTLEN(1) = 1.907 EMCPHOTEFF(1) = 0.03 EMCPHOTLEN(2) = 2.066 EMCPHOTEFF(2) = 0.07 EMCPHOTLEN(3) = 2.254 EMCPHOTEFF(3) = 0.12 EMCPHOTLEN(4) = 2.48 EMCPHOTEFF(4) = 0.20 EMCPHOTLEN(5) = 2.755 EMCPHOTEFF(5) = 0.24 EMCPHOTLEN(6) = 3.100 EMCPHOTEFF(6) = 0.25 EMCPHOTLEN(7) = 3.542 EMCPHOTEFF(7) = 0.24 EMCPHOTLEN(8) = 3.82 EMCPHOTEFF(8) = 0.14 EMCPHOTLEN(9) = 4.133 EMCPHOTEFF(9) = 0.05 NUMED_LEADGLASS = GETMEDNUM(EMCMED1) ! get medium number of lead glass material IF(NUMED_LEADGLASS.LE.0) THEN WRITE(6,'(/'' *** Unknown medium in EMCINIT: '',A15)') EMCMED1 ENDIF RETURN END