G4ProcessHelper.h

Go to the documentation of this file.

#ifndef SimG4Core_CustomPhysics_G4ProcessHelper_H
#define SimG4Core_CustomPhysics_G4ProcessHelper_H

#include "globals.hh"
#include "G4ParticleDefinition.hh"
#include "G4DynamicParticle.hh"
#include "G4Element.hh"
#include "G4Track.hh"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include <vector>
#include <map>

//Typedefs just made to make life easier :-)
typedef std::vector<G4int> ReactionProduct;
typedef std::vector<ReactionProduct> ReactionProductList;
typedef std::map<G4int, ReactionProductList> ReactionMap;

class G4ParticleTable;
class CustomParticleFactory;
class HistoHelper;
class TProfile;
class TH1D;

class G4ProcessHelper {
public:
  G4ProcessHelper(const edm::ParameterSet& p, CustomParticleFactory* ptr);

  ~G4ProcessHelper();

  G4bool ApplicabilityTester(const G4ParticleDefinition& aPart);

  G4double GetInclusiveCrossSection(const G4DynamicParticle* aParticle, const G4Element* anElement);

  //Make sure the element is known (for n/p-decision)
  ReactionProduct GetFinalState(const G4Track& aTrack, G4ParticleDefinition*& aTarget);

  G4ProcessHelper(const G4ProcessHelper&) = delete;
  G4ProcessHelper& operator=(const G4ProcessHelper&) = delete;

private:
  G4double Regge(const double boost);
  G4double Pom(const double boost);

  G4double checkfraction;
  G4int n_22;
  G4int n_23;

  G4ParticleDefinition* theTarget;
  G4ParticleDefinition* theProton;
  G4ParticleDefinition* theNeutron;
  G4ParticleDefinition* theRmesoncloud;
  G4ParticleDefinition* theRbaryoncloud;

  ReactionMap* theReactionMap;

  G4double PhaseSpace(const ReactionProduct& aReaction, const G4DynamicParticle* aDynamicParticle);

  G4double ReactionProductMass(const ReactionProduct& aReaction, const G4DynamicParticle* aDynamicParticle);

  G4bool ReactionGivesBaryon(const ReactionProduct& aReaction);

  G4bool ReactionIsPossible(const ReactionProduct& aReaction, const G4DynamicParticle* aDynamicParticle);

  void ReadAndParse(const G4String& str, std::vector<G4String>& tokens, const G4String& delimiters = " ");

  //Map of applicable particles
  std::map<const G4ParticleDefinition*, G4bool> known_particles;

  //Map for physics parameters, name to value
  std::map<G4String, G4double> parameters;

  //The parameters themselves
  bool resonant;
  double ek_0;
  double gamma;
  double amplitude;
  double suppressionfactor;
  bool reggemodel;
  double hadronlifetime;
  double mixing;

  //Proton-scattering processes
  ReactionMap pReactionMap;

  //Neutron-scattering processes
  ReactionMap nReactionMap;

  CustomParticleFactory* fParticleFactory;
  G4ParticleTable* particleTable;
  HistoHelper* theHistoHelper;
  TProfile* h_xsec_lab;
  TProfile* h_xsec_cms;
  TH1D* h_sqrts;
  TProfile* h_q_p;
  TProfile* h_q_gamma;
};
#endif