GGSCRMCPhysicsList Class Reference

Physics list based on CRMC for high-energy hadron-hadron interactions. More...

#include <GGSCRMCPhysicsList.h>

Inheritance diagram for GGSCRMCPhysicsList:

Collaboration diagram for GGSCRMCPhysicsList:

Public Member Functions
	GGSCRMCPhysicsList ()
	Constructor.
void	ConstructProcess ()
	Constructs the physical processes. More...

Detailed Description

Physics list based on CRMC for high-energy hadron-hadron interactions.

This physics list inherits from FTFP_BERT, and smoothly interpolates the hadronic model from low energies where FTFP_BERT is used to high energies where a CRMC model is used. The user can define the upper and lower limits of the transition region and the CMC model (default values: 1 teV, 10 TeV and DPMJET3).

Definition at line 21 of file GGSCRMCPhysicsList.h.

Member Function Documentation

void GGSCRMCPhysicsList::ConstructProcess

(

)

Constructs the physical processes.

Calls FTFP_Bert::ConstructProcess and then creates the CRMC model.

Definition at line 30 of file GGSCRMCPhysicsList.cpp.

                                          {
  const std::string routineName("GGSCRMCPhysicsList::ConstructProcess");

  if (G4Threading::IsMultithreadedApplication()) {
    G4Exception("GGSCRMCPhysicsList::ConstructProcess", "CRMC is not thread-safe", G4ExceptionSeverity::FatalException,
                "GGSCRMC physics list can't be used in a MT run; please launch GGSPenny with the -st option.");
  }

  FTFP_BERT::ConstructProcess();

  // Set CRMC underlying model
  //   0 - EPOS LHC   1 - EPOS 1.99   12 - DPMJET3
  //   see CRMC documentation for the complete list of
  //   available models
  int crmcModelID = -1;

  // Post-LHC
  if (_modelName == "EPOS-LHC") {
    crmcModelID = 0;
  } else if (_modelName == "QGSJETII-04") {
    crmcModelID = 7;
  } else if (_modelName == "SIBYLL-2.3") {
    crmcModelID = 6;
  }
  // Pre-LHC
  else if (_modelName == "DPMJET") {
    crmcModelID = 12;
  } else if (_modelName == "EPOS-1.99") {
    crmcModelID = 1;
  } else if (_modelName == "QGSJET-01") {
    crmcModelID = 2;
  } else if (_modelName == "QGSJETII-03") {
    crmcModelID = 11;
  }

  else {
    G4Exception("GGSCRMCPhysicsList::ConstructProcess", "0", G4ExceptionSeverity::FatalException,
                std::string("Unknown model: ").append(_modelName).c_str());
  }

  GGSCOUT(INFO) << "CRMC hadronic interaction engine. Model: " << _modelName << ", thresholds: ("
                << _modelLowerThreshold / CLHEP::TeV << "," << _modelUpperThreshold / CLHEP::TeV << ") TeV." << GGSENDL;
  G4HadronPhysicsCRMC *crmcModel = new G4HadronPhysicsCRMC(crmcModelID);
  crmcModel->SetEnergyThreshold(_modelLowerThreshold, _modelUpperThreshold);
  crmcModel->ConstructProcess();
}

---

The documentation for this class was generated from the following files:

• GGSCRMCPhysicsList.h
• GGSCRMCPhysicsList.cpp