CascadeEvents.C File Reference

to link an ouput from cascade and discret gamma-rays. More...

#include <stdio.h>
#include <time.h>
#include "Riostream.h"
#include <cstdlib>
#include <iostream>
#include <fstream>
#include <string>
#include <istream>
#include "TRandom.h"
#include "TTimeStamp.h"
#include "TMatrixD.h"

Go to the source code of this file.

Macros
#define	AAS   23500.0
#define	AC   720.0
#define	APa   33500.0
#define	AS   18560.0
#define	AV   15677.0
#define	MALFA   3727409.0
#define	MN   939573.0
#define	MP   938280.0

Functions
void	CascadeEvents (const Char_t *finput, Int_t nbcas=100, const Char_t *nameOUT="toGEANT101.event")
void	cascataAGS (const Int_t, Int_t, FILE, const Double_t, const BaseGEM)
	This function calculates and writes in a file the gamma cascades. More...
void	cascataAGS (const Int_t r, Int_t nbcas=1, FILE *fileOUT, const Double_t eccemitter, const BaseGEM random[])
Int_t	eccresiduoa (const Double_t &, Int_t &, Int_t &, Double_t &, const Double_t *)
	This function calculates the excitation energy after the emission of an alpha particle. More...
Int_t	eccresiduon (const Double_t &, Int_t &, Int_t &, Double_t &, const Double_t *)
	This function calculate the excitation energy after the emission of a neutron. More...
Int_t	eccresiduop (const Double_t &, Int_t &, Int_t &, Double_t &, const Double_t *)
	This function calculate the excitation energy after the emission of a proton. More...
void	matriceRES (FILE, Int_t *, Double_t **, Double_t *, Int_t, Int_t, Int_t *)
	This function loads the residue nuclei matrix from CASCADE output. More...
void	matriceRES (FILE *fileIN, Int_t massa[29], Double_t pop[29][6], Double_t popsum[29], Int_t &numpari, Int_t &numdispari, Int_t NmenoZ[11])
void	parabola (const Int_t *, const Double_t *, const Int_t &, Double_t &, Double_t &, Double_t &)
	This function calculates the coefficients of the parabola for the calculation of the nucleus mass from tabulated masses of nuclei with the same A. More...
void	parabola (const Int_t x[], const Double_t y[], const Int_t &Nnuc, Double_t &a, Double_t &b, Double_t &c)
Double_t	parabolame (Int_t &, Int_t &, Int_t &)
	This function calculates the mass of a nucleus using the 3 coefficients calculated with void parabola(...) More...
void	posizioneAin_me ()
	A bit of output for debug. More...
void	trovaAZ (const Int_t &, const Int_t &, Int_t &, Int_t &)
	This function calculates A and Z of the nucleus from its position in the residue nuclei matrix. More...
Double_t	trovame (Int_t &, Int_t &)
	This function return nucleus mass for a given A and Z. More...
Double_t	wei1 (Int_t &, Int_t &)
	This function calculates the mass of a nucleus using a semiempirical formula. More...

Variables
Int_t	A [3000]
Int_t	AGS [200]
Int_t	amin
Int_t	massa [29]
Double_t	me [3000]
Int_t	NmenoZ [11]
Int_t	numA
Int_t	numnuclei
Int_t	posizA [3000]
Int_t	rettA
Int_t	rettZ
Int_t	Z [3000]
Int_t	zmin

Detailed Description

to link an ouput from cascade and discret gamma-rays.

The cascade program produces an ouput given the light particle distributions, the entrance (excitation energy, spin) and the cross section for the residual nuclei. This macro reads this input to generate residual nuclei and their associated particles. The GammaWare is used then to add discret gamma-rays.

    To use it in interactive mode you have to: 
    -# Generate a cascade output file with the modificated cascade program (ask Franceso for more details) 
    -# Start a root session: \n
       <em> root </em>
    -# Load the macro : \n
       <em> .L CascadeEvents.C </em>
    -# Start the run : \n
       <em> CascadeEvents("cascade.out", number_of_events, "Agata.input") </em> \n
       where cascade.out is the cascade output file, number_of_events is the number
       of events to simulate and Agata.input is the CascadeEvents output 
       file.

    To use it in batch mode you have to:
    -# Write a batch script (you can find an example: ExCascadeEvents.C)
    -# Start the run: \n
       <em> root -b -q ExCascadeEvents.C > output.txt& </em>

Author

Francesco Recchia (recchia(AT)pd.infn.it)

Definition in file CascadeEvents.C.