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-3.06"), _modelThreshold(300. * CLHEP::GeV), _messenger(this,
        "/GGS/physicsList/crmc/") {
  _messenger.DeclareProperty("modelName", _modelName).SetCandidates(
      "EPOS-LHC QGSJETII-04 SIBYLL-2.3 DPMJET-3.06 EPOS-1.99 QGSJET-01 QGSJETII-03").SetGuidance(
      "The CRMC model to be used for high-energy hadron-hadron interactions");
  _messenger.DeclareProperty("modelThreshold", _modelThreshold).SetUnitCategory("Energy").SetDefaultUnit("GeV").SetGuidance(
      "The lower energy threshold of the CRMC model");
}

void GGSCRMCPhysicsList::ConstructProcess() {
  const std::string routineName("GGSCRMCPhysicsList::ConstructProcess");
  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-3.06") {
    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());
  }

  COUT(INFO) << "CRMC hadronic interaction engine. Model: " << _modelName << ", threshold: "
      << _modelThreshold / CLHEP::GeV << " GeV." << ENDL;
  G4HadronPhysicsCRMC* crmcModel = new G4HadronPhysicsCRMC(crmcModelID);
  crmcModel->SetEnergyThreshold(_modelThreshold);
  crmcModel->ConstructProcess();

}