Utils.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "SimTracker/TrackHistory/interface/CMSProcessTypes.h"
#include "SimTracker/TrackHistory/interface/Utils.h"

G4toCMSLegacyProcTypeMap::G4toCMSLegacyProcTypeMap() {
  // Geant4Process -> CmsProcess
  m_map[0] = CMS::Primary;                // "Primary"         -> "Primary"
  m_map[91] = CMS::Unknown;               // "Transportation"  -> "Unknown"
  m_map[92] = CMS::Unknown;               // "CoupleTrans"     -> "Unknown" => DOUBLE-CHECK THIS
  m_map[1] = CMS::Unknown;                // "CoulombScat"     -> "Unknown" => DOUBLE-CHECK THIS
  m_map[2] = CMS::EIoni;                  // "Ionisation"      -> "EIoni" => WHAT ABOUT MuIoni?
  m_map[3] = CMS::EBrem;                  // "Brems"           -> "EBrem" => WHAT ABOUT MuBrem?
  m_map[4] = CMS::MuPairProd;             // "PairProdCharged" -> "MuPairProd" => DOUBLE-CHECK THIS
  m_map[5] = CMS::Annihilation;           // "Annih"           -> "Annihilation"
  m_map[6] = CMS::Annihilation;           // "AnnihToMuMu"     -> "Annihilation"
  m_map[7] = CMS::Annihilation;           // "AnnihToHad"      -> "Annihilation"
  m_map[8] = CMS::Hadronic;               // "NuclearStopp"    -> "Hadronic" => DOUBLE-CHECK THIS
  m_map[10] = CMS::Unknown;               // "Msc"             -> "Unknown" => DOUBLE-CHECK THIS
  m_map[11] = CMS::Unknown;               // "Rayleigh"        -> "Unknown" => DOUBLE-CHECK THIS
  m_map[12] = CMS::Photon;                // "PhotoElectric"   -> "Photon"
  m_map[13] = CMS::Compton;               // "Compton"         -> "Compton"
  m_map[14] = CMS::Conversions;           // "Conv"            -> "Conversions"
  m_map[15] = CMS::Conversions;           // "ConvToMuMu"      -> "Conversions"
  m_map[21] = CMS::Unknown;               // "Cerenkov"        -> "Unknown" => DOUBLE-CHECK THIS
  m_map[22] = CMS::Unknown;               // "Scintillation"   -> "Unknown" => DOUBLE-CHECK THIS
  m_map[23] = CMS::SynchrotronRadiation;  // "SynchRad"        -> "SynchrotronRadiation"
  m_map[24] = CMS::SynchrotronRadiation;  // "TransRad"        -> "SynchrotronRadiation"
                                          // => DOUBLE-CHECK THIS
  m_map[31] = CMS::Unknown;               // "OpAbsorp"        -> "Unknown" => NOT USED IN CMS
  m_map[32] = CMS::Unknown;               // "OpBoundary"      -> "Unknown" => NOT USED IN CMS
  m_map[33] = CMS::Unknown;               // "OpRayleigh"      -> "Unknown" => NOT USED IN CMS
  m_map[34] = CMS::Unknown;               // "OpWLS"           -> "Unknown" => NOT USED IN CMS
  m_map[35] = CMS::Unknown;               // "OpMieHG"         -> "Unknown" => NOT USED IN CMS
  m_map[51] = CMS::Unknown;               // "DNAElastic"      -> "Unknown" => NOT USED IN CMS
  m_map[52] = CMS::Unknown;               // "DNAExcit"        -> "Unknown" => NOT USED IN CMS
  m_map[53] = CMS::Unknown;               // "DNAIonisation"   -> "Unknown" => NOT USED IN CMS
  m_map[54] = CMS::Unknown;               // "DNAVibExcit"     -> "Unknown" => NOT USED IN CMS
  m_map[55] = CMS::Unknown;               // "DNAAttachment"   -> "Unknown" => NOT USED IN CMS
  m_map[56] = CMS::Unknown;               // "DNAChargeDec"    -> "Unknown" => NOT USED IN CMS
  m_map[57] = CMS::Unknown;               // "DNAChargeInc"    -> "Unknown" => NOT USED IN CMS
  m_map[111] = CMS::Hadronic;             // "HadElastic"      -> "Hadronic"
  m_map[121] = CMS::Hadronic;             // "HadInelastic"    -> "Hadronic"
  m_map[131] = CMS::Hadronic;             // "HadCapture"      -> "Hadronic"
  m_map[141] = CMS::Hadronic;             // "HadFission"      -> "Hadronic"
  m_map[151] = CMS::Hadronic;             // "HadAtRest"       -> "Hadronic"
  m_map[161] = CMS::Hadronic;             // "HadCEX"          -> "Hadronic"
  m_map[201] = CMS::Decay;                // "Decay"           -> "Decay"
  m_map[202] = CMS::Decay;                // "DecayWSpin"      -> "Decay"
  m_map[203] = CMS::Decay;                // "DecayPiWSpin"    -> "Decay"
  m_map[210] = CMS::Decay;                // "DecayRadio"      -> "Decay"
  m_map[211] = CMS::Decay;                // "DecayUnKnown"    -> "Decay"
  m_map[231] = CMS::Decay;                // "DecayExt"        -> "Decay"
  m_map[301] = CMS::Unknown;              // "GFlash"          -> "Unknown"
  m_map[401] = CMS::Unknown;              // "StepLimiter"     -> "Unknown"
  m_map[402] = CMS::Unknown;              // "UsrSpecCuts"     -> "Unknown"
  m_map[403] = CMS::Unknown;              // "NeutronKiller"   -> "Unknown"
}

const unsigned int G4toCMSLegacyProcTypeMap::processId(unsigned int g4ProcessId) const {
  MapType::const_iterator it = m_map.find(g4ProcessId);

  if (it == m_map.end()) {
    edm::LogError("UnknownProcessType") << "Encountered an unknown process type: " << g4ProcessId
                                        << ". Mapping it to 'Undefined'.";
    return CMS::Undefined;
  } else
    return it->second;
}