-#ifndef FSTARS_H
-#define FSTARS_H 1
+#ifndef FSUMCOU_H
+#define FSUMCOU_H 1
#include "cfortran.h"
#include "Rtypes.h"
#include "Fdimpar.h"
extern "C" {
-//*$ create stars.add
-//*copy stars
-//*
-//*=== stars ============================================================*
+//*=== Sumcou ===========================================================*
//*
//*----------------------------------------------------------------------*
//* *
-//* stars common for fluka: *
+//* All additions wrt the original one are: *
+//* Copyright (C) 2005-2005 by Alfredo Ferrari & Paola Sala *
+//* All Rights Reserved. *
+//* *
+//* *
+//* SUMmary COUnters common for FLUKA9x/200x: *
//* *
-//* version for fluka90/99/... of the original one of fluka86 *
+//* (New version of the original STARS comon) *
//* *
//* *
-//* created on 15 may 1990 by alfredo ferrari & paola sala *
-//* infn - milan *
+//* Created on 15 may 1990 by Alfredo Ferrari & Paola Sala *
+//* Infn - Milan *
//* *
-//* last change on 07-oct-99 by alfredo ferrari *
+//* Last change on 09-jul-05 by Alfredo Ferrari *
//* *
//* *
-//* description of the common variable(s): *
+//* Variable description: *
//* *
-//* ntstar = total number of stars generated (modulo 1000000000) *
-//* neulow = total number of low energy neutron interactions *
+//* Ntstar = total number of stars generated (modulo 1000000000) *
+//* Neulow = total number of low energy neutron interactions *
//* (modulo 1000000000) *
-//* numdec = total number of decays (modulo 1000000000) *
-//* mtstar = total number of stars generated / 100000000 *
-//* meulow = total number of low energy neutron interactions *
+//* Numdec = total number of decays (modulo 1000000000) *
+//* Numoph = total number of optical photons generated *
+//* (modulo 1000000000) *
+//* Mtstar = total number of stars generated / 100000000 *
+//* Meulow = total number of low energy neutron interactions *
//* / 1000000000 *
-//* mumdec = total number of decays (modulo 1000000000) *
-//* wstars = total weight of the stars generated *
-//* wneulw = total weight of the low energy neutron interactions *
-//* westar = weight of the stars generated by *
+//* Mumdec = total number of decays / 1000000000 *
+//* Mumoph = total number of optical photons generated /1000000000*
+//* Wneulw(k,j) = total weight of the low energy neutron interactions *
+//* Woptph(k,j) = total weight of generated optical photons *
+//* Westar(k,j) = weight of the stars generated by *
//* different particle types *
-//* westop = total weight of the particles stopped *
-//* wstop = weight of particles of different types stopped *
-//* weisec = total weight of the secondaries created *
-//* weifis = total weight of high energy fissions *
-//* weipri = total weight of the primaries handled *
-//* edi = deposited energy *
-//* edi(1) = by ionisation *
-//* edi(2) = by pi-zeros and/or em cascade *
-//* edi(3) = by nuclear excitation (or nuclear recoil *
+//* Westop(k,j) = weight of particles of different types stopped *
+//* Wedaug(k,j) = weight of the decay products *
+//* Wekill(k,j) = weight of the time-killed particles *
+//* Wedecy(k,j) = weight of the particles decayed *
+//* Weifis(k,j) = weight of high energy fissions generated by particles*
+//* of different kind *
+//* Deccts(k,j) = decay c tau scoring *
+//* Wdecct(k,j) = weight of the particles decayed for which a c tau *
+//* scoring is done *
+//* Weprdc(k,j) = weight of produced particles of different kind *
+//* Wlwnsc(k,j) = weight of the low energy neutrons interaction secon- *
+//* daries *
+//* Wophsc(k,j) = weight of optical photons production/interaction *
+//* secondaries *
+//* Weipri = total weight of the primaries handled *
+//* Edpsco(i,j) = energy deposition scoring (j = 1: prompt particles, *
+//* j = 2: decay particles) *
+//* Edpsco(1,j) = by ionisation *
+//* Edpsco(2,j) = by pi-zeros and/or EM cascade *
+//* Edpsco(3,j) = by nuclear excitation (or nuclear recoil *
//* and heavies if the evaporation module is *
//* activated) *
-//* edi(4) = by stopped particles *
-//* edi(5) = energy leaving the system *
-//* edi(6) = energy carried by discarded particles *
-//* edi(7) = by residual excitation energy (only if the *
+//* Edpsco(4,j) = by stopped particles *
+//* Edpsco(5,j) = energy leaving the system *
+//* Edpsco(6,j) = energy carried by discarded particles *
+//* Edpsco(7,j) = by residual excitation energy (only if the *
//* evaporation module is activated) *
-//* edi(8) = by low energy neutrons (kerma due to low *
+//* Edpsco(8,j) = by low energy neutrons (kerma due to low *
//* energy neutrons transport is in effect) *
-//* edi(9) = energy carried by time killed particles *
-//* edi(10)= energy wasted for nuclear binding energy *
-//* effects for e > 50 mev *
-//* edi(11)= energy wasted for nuclear binding energy *
+//* Edpsco(9,j) = energy carried by time killed particles *
+//* Edpsco(10,j)= energy wasted for nuclear binding energy *
+//* effects for (low energy neutrons excluded) *
+//* Edpsco(11,j)= energy wasted for nuclear binding energy *
//* effects for low energy neutrons *
-//* wdec = weight of the particles decayed *
-//* wdau = weight of the decay-products *
-//* wtkll = weight of the time-killed particles *
-//* wtdec = total weight of the particles decayed *
-//* wtdau = total weight of the decay products *
-//* wttkll = total weight of the time killed particles *
-//* wlwnsc = weight of the low energy neutrons interaction secon- *
-//* daries *
-//* wdecct = weight of the particles decayed for which a c tau *
-//* scoring is done *
-//* deccts = decay c tau scoring *
+//* *
//*----------------------------------------------------------------------*
//*
Double_t edpsco[2][11];
Double_t woptph;
Double_t wneulw;
- Double_t weipri;
+ Double_t weipri;
Int_t numdec;
Int_t ntstar;
Int_t neulow;