GGSCRMCPhysicsList.cpp

/*
 * GGSCRMCPhysicsList.cpp
 *
 *  Created on: 15 Feb 2019
 *      Author: Nicola Mori
 */

#include "montecarlo/physicslists/GGSCRMCPhysicsList.h"
#include "montecarlo/pluginmanagers/GGSMCPluginMacros.h"
#include "utils/GGSSmartLog.h"

#include "G4HadronPhysicsCRMC.hh"

PhysicsListPlugin(GGSCRMCPhysicsList)

    GGSCRMCPhysicsList::GGSCRMCPhysicsList()
    : FTFP_BERT(), _modelName("DPMJET"), _modelUpperThreshold(1. * CLHEP::TeV), _modelLowerThreshold(10. * CLHEP::TeV),
      _messenger(this, "/GGS/physicsList/crmc/") {
  _messenger.DeclareProperty("modelName", _modelName)
      .SetCandidates("EPOS-LHC QGSJETII-04 SIBYLL-2.3 DPMJET EPOS-1.99 QGSJET-01 QGSJETII-03")
      .SetGuidance("The CRMC model to be used for high-energy hadron-hadron interactions");
  _messenger.DeclarePropertyWithUnit(
      "modelLowerThreshold", "TeV", _modelLowerThreshold,
      "The lower energy threshold of the smooth transition region from FTFP_BERT to the CRMC model");
  _messenger.DeclarePropertyWithUnit(
      "modelUpperThreshold", "TeV", _modelUpperThreshold,
      "The upper energy threshold of the smooth transition region from FTFP_BERT to the CRMC model");
}

void GGSCRMCPhysicsList::ConstructProcess() {
  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();
}