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