G4HadParticles.cc

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.          *
// ********************************************************************
//
//
// Geant4 class G4HadParticles
//
// Author V.Ivanchenko 09.05.2020
//
 
#include "G4HadParticles.hh"
 
// p, n, pi+, pi- 
const std::vector<G4int> G4HadParticles::sLightHadrons = {
  2212, 2112, 211, -211
};
 
// Lambda, Sigma+, Sigma-, Xi0, Xi-, Omega-
// (note that Sigma0 has not been included because it decays very quickly)
const std::vector<G4int> G4HadParticles::sHyperons = {
  3122, 3222, 3112, 3322, 3312, 3334
};
 
// anti_Lambda, anti_Sigma+, anti_Sigma-, anti_Xi0, anti_Xi-, anti_Omega- 
// (note that anti_Sigma0 has not been included because it decays very quickly)
const std::vector<G4int> G4HadParticles::sAntiHyperons = {
  -3122, -3222, -3112, -3322, -3312, -3334
};
 
// K+, K-, KS, KL
const std::vector<G4int> G4HadParticles::sKaons = {
  321, -321, 310, 130
};
 
// Note: etac, JPsi, SigmaC++, SigmaC+, SigmaC0, Upsilon,
//       SigmaB+, SigmaB0, SigmaB- are not included because
// they decay very quickly (therefore, their hadronic
// interactions can be neglected, as for pi0 and Sigma0).
const std::vector<G4int> G4HadParticles::sBCHadrons = {
  // D+, D0, D-, D0bar, Ds+, Ds-
  411, 421, -411, -421, 431, -431,
  // B+, B0, B-, B0bar, Bs0, Bs0bar, Bc+, Bc-,
  521, 511, -521, -511, 531, -531, 541, -541,
  // LambdaC+, XiC+, XiC0, OmegaC0
  4122, 4232, 4132, 4332,  
  // LambdaB, XiB0, XiB-, OmegaB-
  5122, 5232, 5132, 5332,
  // corresponding anti_baryons
  -4122, -4232, -4132, -4332,  
  -5122, -5232, -5132, -5332
};
 
// d, t, He3, alpha
const std::vector<G4int> G4HadParticles::sLightIons = {
  1000010020, 1000010030, 1000020030, 1000020040 
};
 
// pbar, nbar, light anti-ions
const std::vector<G4int> G4HadParticles::sLightAntiIons = {
  -2212, -2112, -1000010020, -1000010030, -1000020030, -1000020040
};
 
// hyper_t, hyper_H4, hyper_He4, hyder_He5, 2-hyper-2n, 2-hyper_H4 
const std::vector<G4int> G4HadParticles::sHyperNuclei = {
  1010010030, 1010010040, 1010020040, 1010020050, 1020000040, 1020010040
};
 
// anti-hyper-nuclei
const std::vector<G4int> G4HadParticles::sHyperAntiNuclei = {
  -1010010030, -1010010040, -1010020040, -1010020050, -1020000040, -1020010040 
};
 
//
// charged particles for EM physics
//
const std::vector<G4int> G4HadParticles::sHeavyChargedPart = {
  // Sigma+, Sigma-, Xi-, Omega-, anti_hyperons 
  3222, 3112, 3312, 3334, -3222, -3112, -3312, -3334,
  // light anti_ions 
  -1000010020, -1000010030, -1000020030, -1000020040,
  // tau+-
  15, -15
};
 
const std::vector<G4int> G4HadParticles::sBCChargedHadrons = {
  // D+, D-, Ds+, Ds-
  411, -411, 431, -431,
  // B+, B-, Bc+, Bc-,
  521, -521, 541, -541,
  // LambdaC+, SigmaC++, SigmaC+, XiC+
  4122, 4222, 4212, 4232, 
  // SigmaB+, SigmaB-, XiB-, OmegaB-
  5222, 5112, 5132, 5332,
  // anti_baryons
  -4122, -4222, -4212, -4232, -5222, -5112, -5132, -5332 
};
 
// hyper_t
const std::vector<G4int> G4HadParticles::sChargedHyperNuclei = {
  1010010030,  1010010040,  1010020040,  1010020050,  1020010040,
 -1010010030, -1010010040, -1010020040, -1010020050, -1020010040
};
 
const std::vector<G4int>& G4HadParticles::GetLightHadrons()
{
  return sLightHadrons;
}
 
const std::vector<G4int>& G4HadParticles::GetHyperons()
{
  return sHyperons;
}
 
const std::vector<G4int>& G4HadParticles::GetAntiHyperons()
{
  return sAntiHyperons;
}
 
const std::vector<G4int>& G4HadParticles::GetKaons()
{
  return sKaons;
}
 
const std::vector<G4int>& G4HadParticles::GetBCHadrons()
{
  return sBCHadrons;
}
 
const std::vector<G4int>& G4HadParticles::GetLightIons()
{
  return sLightIons;
}
 
const std::vector<G4int>& G4HadParticles::GetLightAntiIons()
{
  return sLightAntiIons;
}
 
const std::vector<G4int>& G4HadParticles::GetHyperNuclei()
{
  return sHyperNuclei;
}
 
const std::vector<G4int>& G4HadParticles::GetHyperAntiNuclei()
{
  return sHyperAntiNuclei;
}
 
const std::vector<G4int>& G4HadParticles::GetHeavyChargedParticles()
{
  return sHeavyChargedPart;
}
 
const std::vector<G4int>& G4HadParticles::GetBCChargedHadrons()
{
  return sBCChargedHadrons;
}
 
const std::vector<G4int>& G4HadParticles::GetChargedHyperNuclei()
{
  return sChargedHyperNuclei;
}