SimG4Core/CustomPhysics/interface/CustomProcessHelper.h - master - githubmirror/cmssw

SimG4Core/CustomPhysics/interface/CustomProcessHelper.h
91 строка3 KB
Vladimir
rename G4 -> CMS
06 июл 2025, 14:31
06 июл 2025, 14:316987118
#ifndef SimG4Core_CustomPhysics_CustomProcessHelper_H
#define SimG4Core_CustomPhysics_CustomProcessHelper_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 CustomProcessHelper {
public:
  CustomProcessHelper(const edm::ParameterSet& p, CustomParticleFactory* ptr);

  ~CustomProcessHelper() = default;

  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);

  CustomProcessHelper(const CustomProcessHelper&) = delete;
  CustomProcessHelper& operator=(const CustomProcessHelper&) = 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;
};

#endif