TG4ComposedPhysicsList.h

Go to the documentation of this file.

#ifndef TG4_COMPOSED_PHYSICS_LIST_H
#define TG4_COMPOSED_PHYSICS_LIST_H
 
//------------------------------------------------
// The Geant4 Virtual Monte Carlo package
// Copyright (C) 2007 - 2014 Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see geant4_vmc/LICENSE.
// Contact: root-vmc@cern.ch
//-------------------------------------------------

 
#include "TG4ComposedPhysicsMessenger.h"
#include "TG4Verbose.h"
 
#include <G4VUserPhysicsList.hh>
#include <G4HadronicProcessType.hh>
#include <globals.hh>
 
#include <vector>
#include <tuple>

 
class TG4ComposedPhysicsList : public G4VUserPhysicsList, public TG4Verbose
{
 public:
  TG4ComposedPhysicsList();
  virtual ~TG4ComposedPhysicsList();
 
  // methods
  void AddPhysicsList(G4VUserPhysicsList* physicsList);
 
  virtual void ConstructParticle();
  virtual void ConstructProcess();
  virtual void SetCuts();
 
  void SetCutForGamma(G4double cut);
  void SetCutForElectron(G4double cut);
  void SetCutForPositron(G4double cut);
  void SetCutForProton(G4double cut);
 
  void SetProductionCutsTableEnergyRange(G4double min, G4double max);
  void SetGammaToMuonsCrossSectionFactor(G4double value);
  void SetCrossSectionFactor(
         const G4String& particleName, const G4String& processDef,
         G4double value, G4bool isProcessName);
 
  void SetLooperThresholdsLevel(G4int level);
 
  void PrintAllProcesses() const;
  void DumpAllProcesses() const;
 
  virtual G4int VerboseLevel() const;
  virtual void VerboseLevel(G4int level);
 
 private:
  TG4ComposedPhysicsList(const TG4ComposedPhysicsList& right);
  TG4ComposedPhysicsList& operator=(const TG4ComposedPhysicsList& right);
 
  // methods
  void ApplyGammaToMuonsCrossSectionFactor();
  void ApplyCrossSectionFactor(
         const G4String& particleName, const G4String& processDef,
         G4double value, G4bool isProcessName);
  void ApplyCrossSectionFactors();
 
  void SetLooperThresholds();
 
  // static data members
  static const G4double
    fgkDefautLooperThresholdsLevel; 
 
  // data members
  TG4ComposedPhysicsMessenger fMessenger;         
  std::vector<G4VUserPhysicsList*> fPhysicsLists; 

  G4bool fIsProductionCutsTableEnergyRange;

  G4double fProductionCutsTableEnergyMin;

  G4double fProductionCutsTableEnergyMax;

  G4double fGammaToMuonsCrossSectionFactor;

  std::vector<std::tuple<G4String, G4String, G4double, G4bool>> fCrossSectionFactors;

  G4int fLooperThresholdsLevel;
};
 
// inline methods
 
inline void TG4ComposedPhysicsList::SetGammaToMuonsCrossSectionFactor(
  G4double value)
{
  fGammaToMuonsCrossSectionFactor = value;
}
 
inline void TG4ComposedPhysicsList::SetLooperThresholdsLevel(G4int level)
{
  fLooperThresholdsLevel = level;
}
 
#endif // TG4_COMPOSED_MODULAR_PHYSICS_H