[u/mrichter/AliRoot.git] / TFluka / Fsumcou.h

#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
Commit	Line	Data
380c2ff0	1	#ifndef FSUMCOU_H
380c2ff0	2	#define FSUMCOU_H 1
b9d0a01d	3
	4	#include "cfortran.h"
	5	#include "Rtypes.h"
	6
	7	#include "Fdimpar.h"
	8
	9	extern "C" {
380c2ff0	10	//=== Sumcou ===========================================================
b9d0a01d	11	//*
	12	//----------------------------------------------------------------------
	13	//* *
380c2ff0	14	//* All additions wrt the original one are: *
	15	//* Copyright (C) 2005-2005 by Alfredo Ferrari & Paola Sala *
	16	//* All Rights Reserved. *
	17	//* *
	18	//* *
	19	//* SUMmary COUnters common for FLUKA9x/200x: *
b9d0a01d	20	//* *
380c2ff0	21	//* (New version of the original STARS comon) *
b9d0a01d	22	//* *
b9d0a01d	23	//* *
380c2ff0	24	//* Created on 15 may 1990 by Alfredo Ferrari & Paola Sala *
380c2ff0	25	//* Infn - Milan *
b9d0a01d	26	//* *
380c2ff0	27	//* Last change on 09-jul-05 by Alfredo Ferrari *
b9d0a01d	28	//* *
b9d0a01d	29	//* *
380c2ff0	30	//* Variable description: *
b9d0a01d	31	//* *
380c2ff0	32	//* Ntstar = total number of stars generated (modulo 1000000000) *
380c2ff0	33	//* Neulow = total number of low energy neutron interactions *
b9d0a01d	34	//* (modulo 1000000000) *
380c2ff0	35	//* Numdec = total number of decays (modulo 1000000000) *
	36	//* Numoph = total number of optical photons generated *
	37	//* (modulo 1000000000) *
	38	//* Mtstar = total number of stars generated / 100000000 *
	39	//* Meulow = total number of low energy neutron interactions *
b9d0a01d	40	//* / 1000000000 *
380c2ff0	41	//* Mumdec = total number of decays / 1000000000 *
	42	//* Mumoph = total number of optical photons generated /1000000000*
	43	//* Wneulw(k,j) = total weight of the low energy neutron interactions *
	44	//* Woptph(k,j) = total weight of generated optical photons *
	45	//* Westar(k,j) = weight of the stars generated by *
b9d0a01d	46	//* different particle types *
380c2ff0	47	//* Westop(k,j) = weight of particles of different types stopped *
	48	//* Wedaug(k,j) = weight of the decay products *
	49	//* Wekill(k,j) = weight of the time-killed particles *
	50	//* Wedecy(k,j) = weight of the particles decayed *
	51	//* Weifis(k,j) = weight of high energy fissions generated by particles*
	52	//* of different kind *
	53	//* Deccts(k,j) = decay c tau scoring *
	54	//* Wdecct(k,j) = weight of the particles decayed for which a c tau *
	55	//* scoring is done *
	56	//* Weprdc(k,j) = weight of produced particles of different kind *
	57	//* Wlwnsc(k,j) = weight of the low energy neutrons interaction secon- *
	58	//* daries *
	59	//* Wophsc(k,j) = weight of optical photons production/interaction *
	60	//* secondaries *
	61	//* Weipri = total weight of the primaries handled *
	62	//* Edpsco(i,j) = energy deposition scoring (j = 1: prompt particles, *
	63	//* j = 2: decay particles) *
	64	//* Edpsco(1,j) = by ionisation *
	65	//* Edpsco(2,j) = by pi-zeros and/or EM cascade *
	66	//* Edpsco(3,j) = by nuclear excitation (or nuclear recoil *
b9d0a01d	67	//* and heavies if the evaporation module is *
b9d0a01d	68	//* activated) *
380c2ff0	69	//* Edpsco(4,j) = by stopped particles *
	70	//* Edpsco(5,j) = energy leaving the system *
	71	//* Edpsco(6,j) = energy carried by discarded particles *
	72	//* Edpsco(7,j) = by residual excitation energy (only if the *
b9d0a01d	73	//* evaporation module is activated) *
380c2ff0	74	//* Edpsco(8,j) = by low energy neutrons (kerma due to low *
b9d0a01d	75	//* energy neutrons transport is in effect) *
380c2ff0	76	//* Edpsco(9,j) = energy carried by time killed particles *
	77	//* Edpsco(10,j)= energy wasted for nuclear binding energy *
	78	//* effects for (low energy neutrons excluded) *
	79	//* Edpsco(11,j)= energy wasted for nuclear binding energy *
b9d0a01d	80	//* effects for low energy neutrons *
380c2ff0	81	//* *
b9d0a01d	82	//----------------------------------------------------------------------
	83	//*
	84
	85	typedef struct {
6177a6a6	86	Double_t westar[2][nallwp+7];
	87	Double_t westop[2][nallwp+7];
	88	Double_t wedaug[2][nallwp+7];
	89	Double_t wekill[2][nallwp+7];
	90	Double_t wedecy[2][nallwp+7];
	91	Double_t weifis[2][nallwp+7];
	92	Double_t deccts[2][nallwp+7];
	93	Double_t wdecct[2][nallwp+7];
	94	Double_t weprdc[2][nallwp+7];
	95	Double_t wlwnsc[2][4];
	96	Double_t wophsc[2][4];
	97	Double_t edpsco[2][11];
cd4c194d	98	Double_t woptph;
b9d0a01d	99	Double_t wneulw;
380c2ff0	100	Double_t weipri;
b9d0a01d	101	Int_t numdec;
	102	Int_t ntstar;
	103	Int_t neulow;
cd4c194d	104	Int_t numoph;
b9d0a01d	105	Int_t mumdec;
	106	Int_t mtstar;
	107	Int_t meulow;
cd4c194d	108	Int_t mumoph;
6177a6a6	109	} sumcouCommon;
	110	#define SUMCOU COMMON_BLOCK(SUMCOU,sumcou)
	111	COMMON_BLOCK_DEF(sumcouCommon,SUMCOU);
b9d0a01d	112	}
	113
	114	#endif