#ifndef FSUMCOU_H
#define FSUMCOU_H 1

#include "cfortran.h"
#include "Rtypes.h"

#include "Fdimpar.h"

extern "C" {
//*=== Sumcou ===========================================================*
//*
//*----------------------------------------------------------------------*
//*                                                                      *
//*     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:                        *
//*                                                                      *
//*     (New version of the original STARS comon)                        *
//*                                                                      *
//*                                                                      *
//*     Created on    15 may 1990    by    Alfredo Ferrari & Paola Sala  *
//*                                                   Infn - Milan       *
//*                                                                      *
//*     Last change on  09-jul-05    by    Alfredo Ferrari               *
//*                                                                      *
//*                                                                      *
//*     Variable description:                                            *
//*                                                                      *
//*        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)           *
//*        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 / 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(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)                                  *
//*            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)            *
//*            Edpsco(8,j) = by low energy neutrons (kerma due to low    *
//*                          energy neutrons transport is in effect)     *
//*            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             *
//*                                                                      *
//*----------------------------------------------------------------------*
//*

typedef struct {
   Double_t westar[2][nallwp+7];
   Double_t westop[2][nallwp+7];
   Double_t wedaug[2][nallwp+7];
   Double_t wekill[2][nallwp+7];
   Double_t wedecy[2][nallwp+7];
   Double_t weifis[2][nallwp+7];
   Double_t deccts[2][nallwp+7];
   Double_t wdecct[2][nallwp+7];
   Double_t weprdc[2][nallwp+7];
   Double_t wlwnsc[2][4];
   Double_t wophsc[2][4];
   Double_t edpsco[2][11];
   Double_t woptph;
   Double_t wneulw;
   Double_t weipri;
   Int_t    numdec;
   Int_t    ntstar;
   Int_t    neulow;
   Int_t    numoph;
   Int_t    mumdec;
   Int_t    mtstar;
   Int_t    meulow;
   Int_t    mumoph;
} sumcouCommon;
#define SUMCOU COMMON_BLOCK(SUMCOU,sumcou)
COMMON_BLOCK_DEF(sumcouCommon,SUMCOU);
}

#endif