COSMOS v7.655  COSMOSv7655
(AirShowerMC)
cfinTracking.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine cfintracking
 

Function/Subroutine Documentation

◆ cfintracking()

subroutine cfintracking ( )

Definition at line 2 of file cfinTracking.f.

References false, and true.

Referenced by ceventloop(), and chookenevent().

2  implicit none
3 #include "Ztrack.h"
4 #include "Ztrackp.h"
5 #include "Ztrackv.h"
6 #include "Zobs.h"
7 #include "Zobsp.h"
8 #include "Zobsv.h"
9 
10  integer i, ib, ie, is
11  logical dvlpd
12 
13  if( (trace .gt. 0 .and. trace .lt. 60) .or. trace .gt. 100) then
14  close(tracedev) ! close default trace file for 1 event.
15  elseif(trace .gt. 60 .and. trace .lt. 100) then
16  call cputcerenkove ! cerenkov output trailer
17  endif
18 ! for each event
19 ! get size weighted age
20  dvlpd = .false.
21  if(observeas) then
22  if(upgoing) then ! upgoing primary
23  ib = noofassites
24  ie = 1
25  is = -1
26  else
27  ib = 1
28  ie = noofassites
29  is = 1
30  endif
31  do i = ib, ie, is
32  if(asobssites(i)%esize .gt. 0.) then
33  dvlpd = .true.
34  asobssites(i)%age = asobssites(i)%age/asobssites(i)%esize
35  elseif(dvlpd) then
36  asobssites(i)%age = 5.0 ! died out.
37  else
38  asobssites(i)%age = 0.0 ! not yet developed
39  endif
40  enddo
41  endif
nodes i
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon true
Definition: cblkElemag.h:7
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon ! knockon is considered Obsolete *PhotoProd false
Definition: cblkElemag.h:7
block data cblkIncident data *Za1ry is
Definition: cblkIncident.h:5
Here is the caller graph for this function: