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