G4INCLEventInfo.hh

Go to the documentation of this file.

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
// INCL++ intra-nuclear cascade model
// Pekka Kaitaniemi, CEA and Helsinki Institute of Physics
// Davide Mancusi, CEA
// Alain Boudard, CEA
// Sylvie Leray, CEA
// Joseph Cugnon, University of Liege
//
#define INCLXX_IN_GEANT4_MODE 1

#include "globals.hh"

#ifndef G4INCLEVENTINFO_HH
#define G4INCLEVENTINFO_HH 1

#include "G4INCLParticleType.hh"
#ifdef INCL_ROOT_USE
#include <Rtypes.h>
#endif
#include <string>
#include <vector>
#include <algorithm>

namespace G4INCL {
#ifndef INCL_ROOT_USE
    typedef G4int Int_t;
    typedef short Short_t;
    typedef G4float Float_t;
    typedef G4double Double_t;
    typedef G4bool Bool_t;
#endif

    struct EventInfo {
      EventInfo() :
        projectileType(UnknownParticle),
        At(0), Zt(0), Ap(0), Zp(0),
        Ep(0.),
        impactParameter(0.0), nCollisions(0), stoppingTime(0.0),
        EBalance(0.0), pLongBalance(0.0), pTransBalance(0.0),
        nCascadeParticles(0), nRemnants(0), nParticles(0),
        transparent(true),
        forcedCompoundNucleus(false),
        nucleonAbsorption(false), pionAbsorption(false), nDecays(0),
        nBlockedCollisions(0), nBlockedDecays(0),
        effectiveImpactParameter(0.0),
        deltasInside(false),
        forcedDeltasInside(false),
        forcedDeltasOutside(false),
        clusterDecay(false),
        firstCollisionTime(0.),
        firstCollisionXSec(0.),
        nReflectionAvatars(0),
        nCollisionAvatars(0),
        nDecayAvatars(0),
        nUnmergedSpectators(0)
      {
        std::fill_n(ARem, maxSizeRemnants, 0);
        std::fill_n(ZRem, maxSizeRemnants, 0);
        std::fill_n(EStarRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(JRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(EKinRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(pxRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(pyRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(pzRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(thetaRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(phiRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(jxRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(jyRem, maxSizeRemnants, ((Float_t)0.));
        std::fill_n(jzRem, maxSizeRemnants, ((Float_t)0.));

        std::fill_n(A, maxSizeParticles, 0);
        std::fill_n(Z, maxSizeParticles, 0);
        std::fill_n(emissionTime, maxSizeParticles, ((Float_t)0.));
        std::fill_n(EKin, maxSizeParticles, ((Float_t)0.));
        std::fill_n(px, maxSizeParticles, ((Float_t)0.));
        std::fill_n(py, maxSizeParticles, ((Float_t)0.));
        std::fill_n(pz, maxSizeParticles, ((Float_t)0.));
        std::fill_n(theta, maxSizeParticles, ((Float_t)0.));
        std::fill_n(phi, maxSizeParticles, ((Float_t)0.));
        std::fill_n(origin, maxSizeParticles, 0);
      };

      static Int_t eventNumber;
      ParticleType projectileType;

      Short_t At;
      Short_t Zt;

      Short_t Ap;
      Short_t Zp;
      Float_t Ep;

      Float_t impactParameter;
      Int_t nCollisions;
      Float_t stoppingTime;

      Float_t EBalance;
      Float_t pLongBalance;
      Float_t pTransBalance;

      Short_t nCascadeParticles;
      Int_t nRemnants;
      Int_t nParticles;

      Bool_t transparent;
      Bool_t forcedCompoundNucleus;
      Bool_t nucleonAbsorption;
      Bool_t pionAbsorption;
      Int_t nDecays;
      Int_t nBlockedCollisions;
      Int_t nBlockedDecays;
      Float_t effectiveImpactParameter;

      Bool_t deltasInside;
      Bool_t forcedDeltasInside;
      Bool_t forcedDeltasOutside;
      Bool_t clusterDecay;

      Float_t firstCollisionTime;
      Float_t firstCollisionXSec;

      Int_t nReflectionAvatars;
      Int_t nCollisionAvatars;
      Int_t nDecayAvatars;

      Int_t nUnmergedSpectators;

      static const Short_t maxSizeRemnants = 10;
      Short_t ARem[maxSizeRemnants];
      Short_t ZRem[maxSizeRemnants];
      Float_t EStarRem[maxSizeRemnants];
      Float_t JRem[maxSizeRemnants];
      Float_t EKinRem[maxSizeRemnants];
      Float_t pxRem[maxSizeRemnants];
      Float_t pyRem[maxSizeRemnants];
      Float_t pzRem[maxSizeRemnants];
      Float_t thetaRem[maxSizeRemnants];
      Float_t phiRem[maxSizeRemnants];
      Float_t jxRem[maxSizeRemnants];
      Float_t jyRem[maxSizeRemnants];
      Float_t jzRem[maxSizeRemnants];

      static const Short_t maxSizeParticles = 1000;
      Short_t A[maxSizeParticles];
      Short_t Z[maxSizeParticles];
      Float_t emissionTime[maxSizeParticles];
      Float_t EKin[maxSizeParticles];
      Float_t px[maxSizeParticles];
      Float_t py[maxSizeParticles];
      Float_t pz[maxSizeParticles];
      Float_t theta[maxSizeParticles];
      Float_t phi[maxSizeParticles];
      Short_t origin[maxSizeParticles];
      std::vector<std::string> history;

#ifdef INCL_INVERSE_KINEMATICS

      Float_t EKinPrime[maxSizeParticles];
      Float_t pzPrime[maxSizeParticles];
      Float_t thetaPrime[maxSizeParticles];
#endif // INCL_INVERSE_KINEMATICS

      void reset() {
        Ap = 0;
        Zp = 0;
        At = 0;
        Zt = 0;
        impactParameter = 0.0;
        effectiveImpactParameter = 0.0;
        stoppingTime = 0.0;
        EBalance = 0.0;
        pLongBalance = 0.0;
        pTransBalance = 0.0;
        nCollisions = 0;
        nBlockedCollisions = 0;
        nDecays = 0;
        nBlockedDecays= 0;
        nDecays = 0;
        nCascadeParticles = 0;
        nRemnants = 0;
        nParticles = 0;
        transparent = true;
        forcedCompoundNucleus = false;
    nucleonAbsorption = false;
    pionAbsorption = false;
        forcedDeltasInside = false;
        forcedDeltasOutside = false;
        deltasInside = false;
        clusterDecay = false;
        nUnmergedSpectators = 0;
      }

#ifdef INCL_INVERSE_KINEMATICS
      void fillInverseKinematics(const Double_t gamma);
#endif // INCL_INVERSE_KINEMATICS
    };
}

#endif /* G4INCLEVENTINFO_HH */