ciniSPrim.f

Go to the documentation of this file.

!       *******************************************************
!       *  ciniSPrim: initialize csampPrimary
!       *******************************************************
!   Initializes a table for sampling the primary energy
!         and species.  For data to be prepared for the table, see
!         sample.d in Data/Primary directory.
!   NOTE:  Init for angle sampling is not done here. It is 
!          managed in the Tracking directory.
!
        subroutine cinisprim(fn)
!          fn:  character string.  input. file name where primary spectrum
!                                         data is given.  
!                                         See sample.d in Data/Primary
!
        implicit none
#include  "Zptcl.h"
#include  "Zprimary.h"
#include  "Zprimaryc.h"
#include  "Zprimaryv.h"
        character*(*) fn
!
!     
        call csetprimtbl      ! make PrimaryIdTbl
        call cseteunittbl     ! make ErgUnitTbl
        call cinisprim0(prim, fn)
        end
!       *******************************************************
!       *  ciniSPrim0: initialize csampPrimary for a given rec
!       *******************************************************
!   Initializes a table for sampling the primary energy
!         and species.  For data to be prepared for the table, see
!         sample.d in Data/Primary directory.
!
        subroutine cinisprim0(prm, fn)
!          fn:  character string.  input. file name where primary spectrum
!                                         data is given.  
!                                         See sample.d in Data/Primary
!
        implicit none
#include  "Zmanagerp.h"
#include  "Zptcl.h"
#include  "Zprimary.h"
#include  "Zprimaryc.h"

        character*(*) fn
        type(primaries):: prm
        integer icon
!
!     
        call copenf(tempdev,  fn, icon)
        call cskipcomment(tempdev, icon)
        call crdprimdata(prm )
        close(tempdev)
        call cprocprimdt(prm )
        end
!
!     ******************************************     
      subroutine csetprimtbl
         implicit none

#include  "Zcode.h"
#include  "Zptcl.h"
#include  "Zprimary.h"
#include  "Zprimaryc.h"
!
        primaryidtbl(1)%symb = 'gamma'
        primaryidtbl(1)%code = kphoton
        primaryidtbl(1)%subcode = 0
        primaryidtbl(1)%charge = 0

        primaryidtbl(2)%symb = 'photon'
        primaryidtbl(2)%code = kphoton
        primaryidtbl(2)%subcode = 0
        primaryidtbl(2)%charge = 0

        primaryidtbl(3)%symb = 'e'
        primaryidtbl(3)%code = kelec
        primaryidtbl(3)%subcode = regptcl
        primaryidtbl(3)%charge = -1

        primaryidtbl(4)%symb = 'e-'
        primaryidtbl(4)%code = kelec
        primaryidtbl(4)%subcode = regptcl
        primaryidtbl(4)%charge = -1

        primaryidtbl(5)%symb = 'electron'
        primaryidtbl(5)%code = kelec
        primaryidtbl(5)%subcode = regptcl
        primaryidtbl(5)%charge = -1

        primaryidtbl(6)%symb = 'e+'
        primaryidtbl(6)%code = kelec
        primaryidtbl(6)%subcode = antip
        primaryidtbl(6)%charge = 1

        primaryidtbl(7)%symb = 'positron'
        primaryidtbl(7)%code = kelec
        primaryidtbl(7)%subcode = antip
        primaryidtbl(7)%charge = 1

        primaryidtbl(8)%symb = 'mu-'
        primaryidtbl(8)%code = kmuon
        primaryidtbl(8)%subcode = regptcl
        primaryidtbl(8)%charge = -1

        primaryidtbl(9)%symb = 'mu+'
        primaryidtbl(9)%code = kmuon
        primaryidtbl(9)%subcode = antip
        primaryidtbl(9)%charge = 1

        primaryidtbl(10)%symb = 'pi+'
        primaryidtbl(10)%code = kpion
        primaryidtbl(10)%subcode = regptcl
        primaryidtbl(10)%charge = 1

        primaryidtbl(11)%symb = 'pi-'
        primaryidtbl(11)%code = kpion
        primaryidtbl(11)%subcode = antip
        primaryidtbl(11)%charge = -1

        primaryidtbl(12)%symb = 'pi0'
        primaryidtbl(12)%code = kpion
        primaryidtbl(12)%subcode = 0
        primaryidtbl(12)%charge = 0

        primaryidtbl(13)%symb = 'k+'
        primaryidtbl(13)%code = kkaon
        primaryidtbl(13)%subcode = regptcl
        primaryidtbl(13)%charge = 1

        primaryidtbl(14)%symb = 'k-'
        primaryidtbl(14)%code = kkaon
        primaryidtbl(14)%subcode = antip
        primaryidtbl(14)%charge = -1

        primaryidtbl(15)%symb = 'k0l'
        primaryidtbl(15)%code = kkaon
        primaryidtbl(15)%subcode = k0l
        primaryidtbl(15)%charge = 0

        primaryidtbl(16)%symb = 'k0s'
        primaryidtbl(16)%code = kkaon
        primaryidtbl(16)%subcode = k0s
        primaryidtbl(16)%charge = 0

        primaryidtbl(17)%symb = 'proton'
        primaryidtbl(17)%code = knuc
        primaryidtbl(17)%subcode = regptcl
        primaryidtbl(17)%charge = 1

        primaryidtbl(18)%symb = 'p'
        primaryidtbl(18)%code = knuc
        primaryidtbl(18)%subcode = regptcl
        primaryidtbl(18)%charge = 1

        primaryidtbl(19)%symb = 'n'
        primaryidtbl(19)%code = knuc
        primaryidtbl(19)%subcode = kneutron
        primaryidtbl(19)%charge = 0

        primaryidtbl(20)%symb = 'neutron'
        primaryidtbl(20)%code = knuc
        primaryidtbl(20)%subcode = kneutron
        primaryidtbl(20)%charge = 0

        primaryidtbl(21)%symb = 'rho'
        primaryidtbl(21)%code = krho
        primaryidtbl(21)%subcode = 0
        primaryidtbl(21)%charge = 0

        primaryidtbl(22)%symb = 'omega'
        primaryidtbl(22)%code = komega
        primaryidtbl(22)%subcode = 0
        primaryidtbl(22)%charge = 0

        primaryidtbl(23)%symb = 'phi'
        primaryidtbl(23)%code = kphi
        primaryidtbl(23)%subcode = 0
        primaryidtbl(23)%charge = 0

        primaryidtbl(24)%symb = 'd+'
        primaryidtbl(24)%code = kdmes
        primaryidtbl(24)%subcode = regptcl
        primaryidtbl(24)%charge = 1

        primaryidtbl(25)%symb = 'd-'
        primaryidtbl(25)%code = kdmes
        primaryidtbl(25)%subcode = antip
        primaryidtbl(25)%charge = -1

        primaryidtbl(26)%symb = 'd0'
        primaryidtbl(26)%code = kdmes
        primaryidtbl(26)%subcode = regptcl
        primaryidtbl(26)%charge = 0

        primaryidtbl(27)%symb = 'dd~'
        primaryidtbl(27)%code = kddb
        primaryidtbl(27)%subcode = 0
        primaryidtbl(27)%charge = 0

        primaryidtbl(28)%symb = 'nn~'
        primaryidtbl(28)%code = knnb
        primaryidtbl(28)%subcode = 0
        primaryidtbl(28)%charge = 0

        primaryidtbl(29)%symb = 'neu_e'
        primaryidtbl(29)%code = kneue
        primaryidtbl(29)%subcode = regptcl
        primaryidtbl(29)%charge = 0

        primaryidtbl(30)%symb = 'neu_mu'
        primaryidtbl(30)%code = kneumu
        primaryidtbl(30)%subcode = regptcl
        primaryidtbl(30)%charge = 0

        primaryidtbl(31)%symb = 'deuteron'
        primaryidtbl(31)%code = kgnuc
        primaryidtbl(31)%subcode =2
        primaryidtbl(31)%charge = 1

        primaryidtbl(32)%symb = 'alfa'
        primaryidtbl(32)%code = kalfa
        primaryidtbl(32)%subcode = regptcl
        primaryidtbl(32)%charge = 1

        primaryidtbl(33)%symb = 'he'
        primaryidtbl(33)%code = kalfa
        primaryidtbl(33)%subcode = regptcl
        primaryidtbl(33)%charge = 1

        primaryidtbl(34)%symb = 'alpha'
        primaryidtbl(34)%code = kalfa
        primaryidtbl(34)%subcode = regptcl
        primaryidtbl(34)%charge = 1          ! charge is internally set

        primaryidtbl(35)%symb = 'l'
        primaryidtbl(35)%code = klibe
        primaryidtbl(35)%subcode = regptcl
        primaryidtbl(35)%charge = 1

        primaryidtbl(36)%symb = 'libeb'
        primaryidtbl(36)%code = klibe
        primaryidtbl(36)%subcode = regptcl
        primaryidtbl(36)%charge = 1

        primaryidtbl(37)%symb = 'm'
        primaryidtbl(37)%code = kcno
        primaryidtbl(37)%subcode = regptcl
        primaryidtbl(37)%charge = 1

        primaryidtbl(38)%symb = 'cno'
        primaryidtbl(38)%code = kcno
        primaryidtbl(38)%subcode = regptcl
        primaryidtbl(38)%charge = 1

        primaryidtbl(39)%symb = 'h'
        primaryidtbl(39)%code = khvy
        primaryidtbl(39)%subcode = regptcl
        primaryidtbl(39)%charge = 1

        primaryidtbl(40)%symb = 'namgsi'
        primaryidtbl(40)%code = khvy
        primaryidtbl(40)%subcode = regptcl
        primaryidtbl(40)%charge = 1

        primaryidtbl(41)%symb = 'vh'
        primaryidtbl(41)%code = kvhvy
        primaryidtbl(41)%subcode = regptcl
        primaryidtbl(41)%charge = 1

        primaryidtbl(42)%symb = 'sclar'
        primaryidtbl(42)%code = kvhvy
        primaryidtbl(42)%subcode = regptcl
        primaryidtbl(42)%charge = 1

        primaryidtbl(43)%symb = 'fe'
        primaryidtbl(43)%code = kiron
        primaryidtbl(43)%subcode = regptcl
        primaryidtbl(43)%charge = 1

        primaryidtbl(44)%symb = 'iron'
        primaryidtbl(44)%code = kiron
        primaryidtbl(44)%subcode = regptcl
        primaryidtbl(44)%charge = 1

        primaryidtbl(45)%symb = 'iso'
        primaryidtbl(45)%code = kgnuc
        primaryidtbl(45)%subcode = 3    ! to be determined later
        primaryidtbl(45)%charge = 2   ! tbd

        primaryidtbl(46)%symb = 'sigma0'
        primaryidtbl(46)%code = ksigma
        primaryidtbl(46)%subcode = -1
        primaryidtbl(46)%charge = 0

        primaryidtbl(47)%symb = 'sigma+'
        primaryidtbl(47)%code = ksigma
        primaryidtbl(47)%subcode = -1
        primaryidtbl(47)%charge = 1

        primaryidtbl(48)%symb = 'sigma-'
        primaryidtbl(48)%code = ksigma
        primaryidtbl(48)%subcode = -1
        primaryidtbl(48)%charge = -1

        primaryidtbl(49)%symb = 'gzai0'
        primaryidtbl(49)%code = kgzai
        primaryidtbl(49)%subcode = -1
        primaryidtbl(49)%charge = 0

        primaryidtbl(50)%symb = 'gzai+'
        primaryidtbl(50)%code = kgzai
        primaryidtbl(50)%subcode = -1
        primaryidtbl(50)%charge = 1

        primaryidtbl(51)%symb = 'bomega'
        primaryidtbl(51)%code = kbomega
        primaryidtbl(51)%subcode = -1
        primaryidtbl(51)%charge = -1

        primaryidtbl(52)%symb = 'lambda'
        primaryidtbl(52)%code = klambda
        primaryidtbl(52)%subcode = -1
        primaryidtbl(52)%charge = 0

        primaryidtbl(53)%symb = 'lambdac'
        primaryidtbl(53)%code = klambdac
        primaryidtbl(53)%subcode = -1
        primaryidtbl(53)%charge = 1

        primaryidtbl(54)%symb = 'pbar'
        primaryidtbl(54)%code = knuc
        primaryidtbl(54)%subcode = antip
        primaryidtbl(54)%charge = -1

        
        primaryidtbl(55)%symb = 'nbar'
        primaryidtbl(55)%code = knuc
        primaryidtbl(55)%subcode = antip
        primaryidtbl(55)%charge = 0

        primaryidtbl(56)%symb = 'light'
        primaryidtbl(56)%code = klight
        primaryidtbl(56)%subcode = 1
        primaryidtbl(56)%charge = 0

        primaryidtbl(57)%symb = 'Edepo'
        primaryidtbl(57)%code =  kedepo
        primaryidtbl(57)%subcode = 0
        primaryidtbl(57)%charge = 0

        primaryidtbl(58)%symb = 'chgPath'
        primaryidtbl(58)%code =  kchgpath
        primaryidtbl(58)%subcode = 0
        primaryidtbl(58)%charge = 1

        primaryidtbl(59)%symb = 'tau+'
        primaryidtbl(59)%code =  ktau
        primaryidtbl(59)%subcode = antip
        primaryidtbl(59)%charge = 1

        primaryidtbl(60)%symb = 'tau-'
        primaryidtbl(60)%code =  ktau
        primaryidtbl(60)%subcode =  regptcl
        primaryidtbl(60)%charge = -1

        primaryidtbl(61)%symb = 'neu_tau'
        primaryidtbl(61)%code =  kneutau
        primaryidtbl(61)%subcode =  regptcl
        primaryidtbl(61)%charge = 0

        primaryidtbl(62)%symb = 'neu_e~'
        primaryidtbl(62)%code = kneue
        primaryidtbl(62)%subcode = antip
        primaryidtbl(62)%charge = 0

        primaryidtbl(63)%symb = 'neu_mu~'
        primaryidtbl(63)%code = kneumu
        primaryidtbl(63)%subcode = antip
        primaryidtbl(63)%charge = 0

        primaryidtbl(64)%symb = 'neu_tau~'
        primaryidtbl(64)%code =  kneutau
        primaryidtbl(64)%subcode =  antip
        primaryidtbl(64)%charge = 0

        primaryidtbl(65)%symb = 'pdg'
        primaryidtbl(65)%code = 0       ! later fixed
        primaryidtbl(65)%subcode = 0    ! //
        primaryidtbl(65)%charge = 0     !

        primaryidtbl(66)%symb = 'p~'
        primaryidtbl(66)%code = knuc
        primaryidtbl(66)%subcode = antip
        primaryidtbl(66)%charge = -1

        primaryidtbl(67)%symb = 'n~'
        primaryidtbl(67)%code = knuc
        primaryidtbl(67)%subcode = antip
        primaryidtbl(67)%charge = 0
        
        
       end
!***************************************************
       subroutine cseteunittbl
         implicit none

#include  "Zptcl.h"
#include  "Zprimary.h"
#include  "Zprimaryc.h"

         ergunittbl(1)%symb = 'ev'      
         ergunittbl(1)%togev = 1.d-9
        
         ergunittbl(2)%symb = 'kev'
         ergunittbl(2)%togev = 1.d-6
        
         ergunittbl(3)%symb = 'mev'
         ergunittbl(3)%togev = 1.d-3
        
         ergunittbl(4)%symb = 'gev'
         ergunittbl(4)%togev = 1.d0
        
         ergunittbl(5)%symb = 'tev'
         ergunittbl(5)%togev = 1.d3
        
         ergunittbl(6)%symb = 'pev'
         ergunittbl(6)%togev = 1.d6
        
         ergunittbl(7)%symb = 'eev'
         ergunittbl(7)%togev = 1.d9
       end