#include "CbmVMCSettings.h"

#include <FairLogger.h>

#include <TVirtualMC.h>

void CbmVMCSettings::Init(TVirtualMC* vmc) {

  LOG(info) << ": Configuring global VMC settings";
  assert(vmc);

  // TODO: These settings were taken over from SetCuts.C. Their meanings
  // will have to be looked up and documented.

  
  // The processes with the following names are defined for TGeant3 and TGeant4
  // such that they can be used similarily for both transport engines
  // PAIR, COMP, PHOT, PFIS, DRAY, ANNI, BREM, HADR, MUNU,
  // DCAY, LOSS, MULS

  // Also the following processes are defined for both transport engines. To unknown
  // reasons they are only applied for Geant3  
  // RAYL, CKOV

  // The following process is not set for any of the two transport engines
  // For Geant3 the default value is used
  // SYNC
  
  // This ppocess is only defined for TGeant4
  // Somehow it is also available in the Geant3 Fortran code but it is not forwarded
  // from the TGeant3 C++ Wraper
  // LABS

  // The energy loss straggling can be only set for Geant3. For Geant4 there is a code fragment
  // which produces a warning informing the user that this property is only
  // supported via /mcPhysics/emModel commands of Geant4
  // STRA
   
  // Processes
  vmc->SetProcess("PAIR", fProcessPairProduction);         // pair production
  vmc->SetProcess("COMP", fProcessComptonScattering);      // Compton scattering
  vmc->SetProcess("PHOT", fProcessPhotoEffect);            // photo electric effect
  vmc->SetProcess("PFIS", fProcessPhotoFission);           // photofission
  vmc->SetProcess("DRAY", fProcessDeltaRay);               // delta-ray
  vmc->SetProcess("ANNI", fProcessAnnihilation);           // annihilation
  vmc->SetProcess("BREM", fProcessBremsstrahlung);         // bremsstrahlung
  vmc->SetProcess("HADR", fProcessHadronicInteraction);    // hadronic process
  vmc->SetProcess("MUNU", fProcessMuonNuclearInteraction); // muon nuclear interaction
  vmc->SetProcess("DCAY", fProcessDecay);                  // decay
  vmc->SetProcess("LOSS", fProcessEnergyLossModel);        // energy loss
  vmc->SetProcess("MULS", fProcessMultipleScattering);     // multiple scattering

  // The cuts with the following names are defined for TGeant3 and TGeant4
  // such that they can be used similarily for both transport engines
  // The values given will be taken as they are for Geant3
  // For Geant4 the values are used to calculate the proper values in the format
  // needed by Geant
  // CUTGAM, CUTELE, CUTNEU ,CUTHAD, CUTMUO, BCUTE, BCUTM
  // DCUTE, DCUTM, PPCUTM, TOFMAX

  vmc->SetCut("CUTGAM", fEnergyCutGammas);                   // gammas (GeV)
  vmc->SetCut("CUTELE", fEnergyCutElectrons);                // electrons (GeV)
  vmc->SetCut("CUTNEU", fEnergyCutNeutralHadrons);           // neutral hadrons (GeV)
  vmc->SetCut("CUTHAD", fEnergyCutChargedHadrons);           // charged hadrons (GeV)
  vmc->SetCut("CUTMUO", fEnergyCutMuons);                    // muons (GeV)
  vmc->SetCut("BCUTE",  fEnergyCutElectronBremsstrahlung);   // electron bremsstrahlung (GeV)
  vmc->SetCut("BCUTM",  fEnergyCutMuonHadronBremsstrahlung); // muon/hadron bremsstrahlung(GeV)
  vmc->SetCut("DCUTE",  fEnergyCutElectronDeltaRay);         // delta-rays by electrons (GeV)
  vmc->SetCut("DCUTM",  fEnergyCutMuonDeltaRay);             // delta-rays by muons (GeV)
  vmc->SetCut("PPCUTM", fEnergyCutMuonPairProduction); // direct pair production by muons (GeV)
  vmc->SetCut("TOFMAX", fTimeCutTof);                        // time of flight cut in seconds
}

ClassImp(CbmVMCSettings);