UtilIstHEP.cxx

Go to the documentation of this file.

//
// UtilIstHEP
//
// Package: Util
//
// S. Kasahara 08/07
//
// Purpose: Utility methods and enumerated types for purposes of 
//          interpreting IstHEP Codes used in the MINOS Monte Carlo.
//
#include <cassert>

#include <TString.h>
#include <TMCProcess.h>

#include "Util/UtilIstHEP.h"
#include "MessageService/MsgService.h"

CVSID("$Id: UtilIstHEP.cxx,v 1.7 2009/02/28 21:46:17 gmieg Exp $");

// Definition of methods 
// *********************

//_____________________________________________________________________________
const char* UtilIstHEP::AsString(ENeugenIstHEP neugenisthep) {
  //  Purpose:  Convert enumerated Neugen IstHEP code to a string.
  //  Return:  isthep string. If unknown, returns "Unknown".

  switch ( neugenisthep ) {

    case kNInitial:
      return "Initial";
    case kNFinal:
      return "Final";
    case kNIntermediate:
      return "Intermediate";
    case kNDecayed:
      return "Decayed";
    case kNTarget:
      return "Target";
    case kNDIS:
      return "DIS";
    case kNRes:
      return "Res";
    case kNRescatter:
      return "Rescatter";
    default:
      MSG("Util",Msg::kWarning) 
      << "UtilIstHEP::AsString called with unknown ENeugenIstHEP " 
      << (Int_t)neugenisthep << endl;
      return "Unknown";

  }//end of switch

}

//_____________________________________________________________________________
UtilIstHEP::EProdMethod UtilIstHEP::GetProdMethod(TMCProcess process) {
  // Purpose: Interpret a TMCProcess as an EProdMethod

  switch (process) {

  case kPPrimary:
    return kMPrimary;
  case kPMultipleScattering:
    return kMMultipleScattering;
  case kPEnergyLoss:         
    return kMEnergyLoss;
  case kPMagneticFieldL:
    return kMMagneticFieldL;
  case kPDecay:
    return kMDecay;
  case kPPair:
    return kMPair;
  case kPCompton:
    return kMCompton;
  case kPPhotoelectric:
    return kMPhotoelectric;
  case kPBrem:
    return kMBrem;
  case kPDeltaRay:
    return kMDeltaRay;
  case kPAnnihilation:
    return kMAnnihilation;
  case kPHadronic:
    return kMHadronic;
  case kPEvaporation:
    return kMEvaporation;
  case kPNuclearFission:
    return kMNuclearFission;
  case kPNuclearAbsorption:
    return kMNuclearAbsorption;
  case kPPbarAnnihilation:
    return kMPbarAnnihilation;  
  case kPNCapture:
    return kMNCapture;
  case kPHElastic:
    return kMHElastic;
  case kPHInhelastic:
    return kMHInhelastic;
  case kPMuonNuclear:
    return kMMuonNuclear;
  case kPTOFlimit:
    return kMTOFlimit;
  case kPPhotoFission:
    return kMPhotoFission;
  case kPRayleigh:
    return kMRayleigh;
  case kPNull:
    return kMNull;
  case kPStop:
    return kMStop;
  case kPLightAbsorption:
    return kMLightAbsorption;
  case kPLightDetection:
    return kMLightDetection;
  case kPLightScattering:
    return kMLightScattering;
  case kStepMax:
    return kMStepMax;
  case kPCerenkov:
    return kMCerenkov;
  case kPFeedBackPhoton:
    return kMFeedBackPhoton;
  case kPLightReflection:
    return kMLightReflection;
  case kPLightRefraction:
    return kMLightRefraction;
  case kPSynchrotron:
    return kMSynchrotron;
  case kPTransportation:
    return kMTransportation;
  case kPNoProcess:
    return kMNoProcess;
#if ROOT_VERSION_CODE >= ROOT_VERSION(5,17,06)
  // These only appeared in TMCProcess beginning with ROOT v5.17/06
  case kPAnnihilationRest:
    return kMAnnihilationRest;
  case kPAnnihilationFlight:
    return kMAnnihilationFlight;
  case kPNbarAnnihilation:
    return kMNbarAnnihilation;
  case kPHIElastic:
    return kMHIElastic;
  case kPHCElastic:
    return kMHCElastic;
  case kPPhotonInhelastic:
    return kMPhotonInhelastic;
  case kPElectronNuclear:
    return kMElectronNuclear;
  case kPPositronNuclear:
    return kMPositronNuclear;
  case kPScintillation:
    return kMScintillation;
#endif
  default:
    MSG("Util",Msg::kError) << "UtilIstHEP::GetProdMethod for TMCProcess " 
                            << process << "/" << TMCProcessName[process] 
                            << " failed!\nUpdate UtilIstHEP to include " 
                            << "missing process. Abort." << endl;
    abort();
  }
    
}

//_____________________________________________________________________________
TMCProcess UtilIstHEP::GetTMCProcess(EProdMethod prodmethod) {
  // Purpose: Interpret a EProdMethod as a TMCProcess

  switch (prodmethod) {

  case kMPrimary:
    return kPPrimary;
  case kMMultipleScattering:
    return kPMultipleScattering;
  case kMEnergyLoss:         
    return kPEnergyLoss;
  case kMMagneticFieldL:
    return kPMagneticFieldL;
  case kMDecay:
    return kPDecay;
  case kMPair:
    return kPPair;
  case kMCompton:
    return kPCompton;
  case kMPhotoelectric:
    return kPPhotoelectric;
  case kMBrem:
    return kPBrem;
  case kMDeltaRay:
    return kPDeltaRay;
  case kMAnnihilation:
    return kPAnnihilation;
  case kMHadronic:
    return kPHadronic;
  case kMEvaporation:
    return kPEvaporation;
  case kMNuclearFission:
    return kPNuclearFission;
  case kMNuclearAbsorption:
    return kPNuclearAbsorption;
  case kMPbarAnnihilation:
    return kPPbarAnnihilation;  
  case kMNCapture:
    return kPNCapture;
  case kMHElastic:
    return kPHElastic;
  case kMHInhelastic:
    return kPHInhelastic;
  case kMMuonNuclear:
    return kPMuonNuclear;
  case kMTOFlimit:
    return kPTOFlimit;
  case kMPhotoFission:
    return kPPhotoFission;
  case kMRayleigh:
    return kPRayleigh;
  case kMNull:
    return kPNull;
  case kMStop:
    return kPStop;
  case kMLightAbsorption:
    return kPLightAbsorption;
  case kMLightDetection:
    return kPLightDetection;
  case kMLightScattering:
    return kPLightScattering;
  case kMStepMax:
    return kStepMax;
  case kMCerenkov:
    return kPCerenkov;
  case kMFeedBackPhoton:
    return kPFeedBackPhoton;
  case kMLightReflection:
    return kPLightReflection;
  case kMLightRefraction:
    return kPLightRefraction;
  case kMSynchrotron:
    return kPSynchrotron;
  case kMTransportation:
    return kPTransportation;
  case kMNoProcess:
    return kPNoProcess;
#if ROOT_VERSION_CODE >= ROOT_VERSION(5,17,06)
  // These only appeared in TMCProcess beginning with ROOT v5.17/06
  case kMAnnihilationRest:
    return kPAnnihilationRest;
  case kMAnnihilationFlight:
    return kPAnnihilationFlight;
  case kMNbarAnnihilation:
    return kPNbarAnnihilation;
  case kMHIElastic:
    return kPHIElastic;
  case kMHCElastic:
    return kPHCElastic;
  case kMPhotonInhelastic:
    return kPPhotonInhelastic;
  case kMElectronNuclear:
    return kPElectronNuclear;
  case kMPositronNuclear:
    return kPPositronNuclear;
  case kMScintillation:
    return kPScintillation;
#endif
  default:
    MSG("Util",Msg::kError) << "UtilIstHEP::GetTMCProcess for EProdMethod " 
                            << prodmethod
                            << " failed!\nUpdate UtilIstHEP to include " 
                            << "missing production method. Abort." << endl;
    abort();
  }
    
}

  
//_____________________________________________________________________________
Int_t UtilIstHEP::GetNeugenIstHEP(const char* isthepname) {
  // Purpose: Convert text string isthepname to an enumerated code.
  // 
  // Argument: isthepname string
  //
  // Return: returns -1 if no match
  //

  TString tmpstr(isthepname);
  tmpstr.ToLower();
  if ( strcmp(tmpstr.Data(),"initial") == 0 ) return UtilIstHEP::kNInitial;
  else if ( strcmp(tmpstr.Data(),"final") == 0 ) return UtilIstHEP::kNFinal;
  else if ( strcmp(tmpstr.Data(),"intermediate") == 0 ) 
                                           return UtilIstHEP::kNIntermediate;
  else if ( strcmp(tmpstr.Data(),"decayed") == 0 ) 
                                           return UtilIstHEP::kNDecayed;
  else if ( strcmp(tmpstr.Data(),"target") == 0 ) return UtilIstHEP::kNTarget;
  else if ( strcmp(tmpstr.Data(),"dis") == 0 ) return UtilIstHEP::kNDIS;
  else if ( strcmp(tmpstr.Data(),"res") == 0 ) return UtilIstHEP::kNRes;
  else if ( strcmp(tmpstr.Data(),"rescatter") == 0 ) 
                                               return UtilIstHEP::kNRescatter;

  return -1;

}

//_____________________________________________________________________________
const char* UtilIstHEP::AsString(EProdMethod prodmethod) {
  // Purpose: Convert production method status code to a string

  TMCProcess processcode = GetTMCProcess(prodmethod);
  return TMCProcessName[processcode];
  
}

//_____________________________________________________________________________
const char* UtilIstHEP::AsString(Int_t isthep) {
  //  Purpose:  Convert isthep status code to a string.
  //
  //  Argument: isthep code.
  //
  //  Return:  isthep string. If unknown, returns "Unknown".
  //
  //  Notes: If isthep <  200 => Neugen IstHEP code
  //         If isthep >= kProdMethodOffset 
  //         && isthep <  kProdMethodOffset + kNProdMethod  
  //                          => EProdMethod production method code
  //         If isthep == 999 => "Info" 
  //         If isthep == kProdMethodOffset + kMDecay + 1000
  //                      => "Decay Parent" 
  //         Otherwise return "Unknown"


  if ( isthep < 200 ) {
    ENeugenIstHEP neugenisthep = (ENeugenIstHEP)isthep;
    return AsString(neugenisthep);
  }
  else if ( isthep >= kProdMethodOffset 
         && isthep <  kProdMethodOffset + kNProdMethod ) {
    EProdMethod prodmethod = (EProdMethod)(isthep-kProdMethodOffset);
    return AsString(prodmethod);
  }
  else if ( isthep == 999 ) {
    return "Info";
  }
  else if ( isthep == (kProdMethodOffset + kMDecay + 1000) ) {
    return "Decay Parent";
  }
  else {
    //MSG("Util",Msg::kWarning) 
    //  << "UtilIstHEP::AsString called with unknown isthep code " 
    //  << isthep << endl;
    return "Unknown";
  }

}

---