HadronicProcessHelper.h

Go to the documentation of this file.

#ifndef SimG4Core_CustomPhysics_HadronicProcessHelper_H
#define SimG4Core_CustomPhysics_HadronicProcessHelper_H
 
#include "globals.hh"
#include "G4ParticleDefinition.hh"
#include "G4DynamicParticle.hh"
#include "G4Element.hh"
#include "G4Track.hh"
 
#include <vector>
#include <map>
#include <string>
 
class G4ParticleTable;
 
class HadronicProcessHelper {
public:
  //Typedefs just made to make life easier :-)
  typedef std::vector<G4int> ReactionProduct;
  typedef std::vector<ReactionProduct> ReactionProductList;
  typedef std::map<G4int, ReactionProductList> ReactionMap;
 
  HadronicProcessHelper(const std::string& fileName);
 
  HadronicProcessHelper(const HadronicProcessHelper&);
 
  HadronicProcessHelper& operator=(const HadronicProcessHelper&);
 
  HadronicProcessHelper* instance() { return this; };
 
  G4bool applicabilityTester(const G4ParticleDefinition& particle);
 
  G4double inclusiveCrossSection(const G4DynamicParticle* particle, const G4Element* element);
 
  //Make sure the element is known (for n/p-decision)
  ReactionProduct finalState(const G4Track& track, G4ParticleDefinition*& target) {
    return finalState(track.GetDynamicParticle(), track.GetMaterial(), target);
  }
  ReactionProduct finalState(const G4DynamicParticle* particle,
                             const G4Material* material,
                             G4ParticleDefinition*& target);
 
private:
  G4ParticleDefinition* m_proton;
  G4ParticleDefinition* m_neutron;
 
  //  ReactionMap* m_reactionMap;
 
  G4double m_phaseSpace(const ReactionProduct& aReaction,
                        const G4DynamicParticle* aDynamicParticle,
                        G4ParticleDefinition* target);
 
  G4double m_reactionProductMass(const ReactionProduct& aReaction,
                                 const G4DynamicParticle* aDynamicParticle,
                                 G4ParticleDefinition* target);
 
  G4bool m_reactionIsPossible(const ReactionProduct& aReaction,
                              const G4DynamicParticle* aDynamicParticle,
                              G4ParticleDefinition* target);
 
  void m_readAndParse(const G4String& str, std::vector<G4String>& tokens, const G4String& delimiters = " ");
 
  //Map of applicable particles
  std::map<const G4ParticleDefinition*, G4bool> m_knownParticles;
 
  //Proton-scattering processes
  ReactionMap m_protonReactionMap;
 
  //Neutron-scattering processes
  ReactionMap m_neutronReactionMap;
 
  G4ParticleTable* m_particleTable;
 
  G4double m_checkFraction;
  G4int m_n22;
  G4int m_n23;
};
 
#endif