UtilMCFlag.cxx

Go to the documentation of this file.

//
// UtilMCFlag
//
// Package: Util
//
// S. Kasahara 01/2008
//
// Purpose: Utility methods and enumerated types for purposes of 
//          interpreting MC flags used in configuring the MINOS VMC
//          Monte Carlo.
//

#include <TString.h>
#include "Util/UtilMCFlag.h"
#include "MessageService/MsgService.h"

CVSID("$Id: UtilMCFlag.cxx,v 1.2 2008/06/12 19:52:10 rhatcher Exp $");

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

//_____________________________________________________________________________
const char* UtilMCFlag::AsString(ECut cut) {
  //  Purpose:  Convert enumerated mc cut to a cut string.
  //  If unknown, return "UnknownCut".

  switch ( cut ) {

  case kCUTGAM:
    return "CUTGAM";
  case kCUTELE:
    return "CUTELE";
  case kCUTNEU:
    return "CUTNEU";
  case kCUTHAD:
    return "CUTHAD";
  case kCUTMUO:
    return "CUTMUO";
  case kBCUTE:
    return "BCUTE";
  case kBCUTM:
    return "BCUTM";
  case kDCUTE:
    return "DCUTE";
  case kDCUTM:
    return "DCUTM";
  case kPPCUTM:
    return "PPCUTM";
  case kTOFMAX:
    return "TOFMAX";
  case kUnknownCut:
    return "UnknownCut";
  default:
    MSG("UtilMCFlag",Msg::kWarning) 
      << "UtilMCFlag::AsString called with unknown cut " << (int)cut << endl;
    return "UnknownCut";
  }//end of switch

}

//_____________________________________________________________________________
UtilMCFlag::ECut UtilMCFlag::GetCut(const char* cut) {
  // Purpose: Convert text string cut to an enumerated code.
  // Return kUnknownCut if unknown

  TString tmpstr(cut);
  tmpstr.ToLower();
  if ( strcmp(tmpstr.Data(),"cutgam") == 0 ) return UtilMCFlag::kCUTGAM;
  else if ( strcmp(tmpstr.Data(),"cutele") == 0 ) return UtilMCFlag::kCUTELE;
  else if ( strcmp(tmpstr.Data(),"cutneu") == 0 ) return UtilMCFlag::kCUTNEU;
  else if ( strcmp(tmpstr.Data(),"cuthad") == 0 ) return UtilMCFlag::kCUTHAD;
  else if ( strcmp(tmpstr.Data(),"cutmuo") == 0 ) return UtilMCFlag::kCUTMUO;
  else if ( strcmp(tmpstr.Data(),"bcute") == 0 ) return UtilMCFlag::kBCUTE;
  else if ( strcmp(tmpstr.Data(),"bcutm") == 0 ) return UtilMCFlag::kBCUTM;
  else if ( strcmp(tmpstr.Data(),"dcute") == 0 ) return UtilMCFlag::kDCUTE;
  else if ( strcmp(tmpstr.Data(),"dcutm") == 0 ) return UtilMCFlag::kDCUTM;
  else if ( strcmp(tmpstr.Data(),"ppcutm") == 0 ) return UtilMCFlag::kPPCUTM;
  else if ( strcmp(tmpstr.Data(),"tofmax") == 0 ) return UtilMCFlag::kTOFMAX;
  return kUnknownCut;

}

//_____________________________________________________________________________
const char* UtilMCFlag::AsString(EProcess process) {
  //  Purpose:  Convert enumerated mc process to a process string.
  //  If unknown process, return "UnknownProcess".

  switch ( process ) {

  case kANNI:
    return "ANNI";
  case kBREM:
    return "BREM";
  case kCOMP:
    return "COMP";
  case kDCAY:
    return "DCAY";
  case kDRAY:
    return "DRAY";
  case kHADR:
    return "HADR";
  case kLOSS:
    return "LOSS";
  case kMULS:
    return "MULS";
  case kMUNU:
    return "MUNU";
  case kPAIR:
    return "PAIR";
  case kPFIS:
    return "PFIS";
  case kPHOT:
    return "PHOT";
  case kUnknownProcess:
    return "UnknownProcess";
  default:
    MSG("UtilMCFlag",Msg::kWarning) 
      << "UtilMCFlag::AsString called with unknown process "  << (int)process 
      << endl;
    return "UnknownProcess";
  } //end of switch

}

//_____________________________________________________________________________
UtilMCFlag::EProcess UtilMCFlag::GetProcess(const char* process) {
  // Purpose: Convert text string process to an enumerated code.
  // If no match, return kUnknownProcess.

  TString tmpstr(process);
  tmpstr.ToLower();
  if ( strcmp(tmpstr.Data(),"anni") == 0 ) return UtilMCFlag::kANNI;
  else if ( strcmp(tmpstr.Data(),"brem") == 0 ) return UtilMCFlag::kBREM;
  else if ( strcmp(tmpstr.Data(),"comp") == 0 ) return UtilMCFlag::kCOMP;
  else if ( strcmp(tmpstr.Data(),"dcay") == 0 ) return UtilMCFlag::kDCAY;
  else if ( strcmp(tmpstr.Data(),"dray") == 0 ) return UtilMCFlag::kDRAY;
  else if ( strcmp(tmpstr.Data(),"hadr") == 0 ) return UtilMCFlag::kHADR;
  else if ( strcmp(tmpstr.Data(),"loss") == 0 ) return UtilMCFlag::kLOSS;
  else if ( strcmp(tmpstr.Data(),"muls") == 0 ) return UtilMCFlag::kMULS;
  else if ( strcmp(tmpstr.Data(),"munu") == 0 ) return UtilMCFlag::kMUNU;
  else if ( strcmp(tmpstr.Data(),"pair") == 0 ) return UtilMCFlag::kPAIR;
  else if ( strcmp(tmpstr.Data(),"pfis") == 0 ) return UtilMCFlag::kPFIS;
  else if ( strcmp(tmpstr.Data(),"phot") == 0 ) return UtilMCFlag::kPHOT;

  return kUnknownProcess;

}

//_____________________________________________________________________________
const char* UtilMCFlag::AsString(ETracking tracking) {
  //  Purpose:  Convert enumerated tracking parameter to a tracking string.
  //  If unknown tracking parameter, return "UnknownTracking".

  switch ( tracking ) {

  case kISVOL:
    return "ISVOL";
  case kIFIELD:
    return "IFIELD";
  case kFIELDM:
    return "FIELDM";
  case kTMAXFD:
    return "TMAXFD";
  case kSTEMAX:
    return "STEMAX";
  case kDEEMAX:
    return "DEEMAX";
  case kEPSIL:
    return "EPSIL";
  case kSTMIN:
    return "STMIN";
  case kUnknownTracking:
    return "UnknownTracking";
  default:
    MSG("UtilMCFlag",Msg::kWarning) 
      << "UtilMCFlag::AsString called with unknown tracking parameter "  
      << (int)tracking << endl;
    return "UnknownTracking";
  } //end of switch

}

//_____________________________________________________________________________
UtilMCFlag::ETracking UtilMCFlag::GetTracking(const char* tracking) {
  // Purpose: Convert text string tracking parameter to an enumerated code.
  // If no match, return kUnknownTracking.

  TString tmpstr(tracking);
  tmpstr.ToLower();
  if ( strcmp(tmpstr.Data(),"isvol") == 0 ) return kISVOL;
  else if ( strcmp(tmpstr.Data(),"ifield") == 0 ) return kIFIELD;
  else if ( strcmp(tmpstr.Data(),"fieldm" ) == 0 ) return kFIELDM;
  else if ( strcmp(tmpstr.Data(),"tmaxfd" ) == 0 ) return kTMAXFD;
  else if ( strcmp(tmpstr.Data(),"stemax" ) == 0 ) return kSTEMAX;
  else if ( strcmp(tmpstr.Data(),"deemax" ) == 0 ) return kDEEMAX;
  else if ( strcmp(tmpstr.Data(),"epsil" ) == 0 ) return kEPSIL;
  else if ( strcmp(tmpstr.Data(),"stmin" ) == 0 ) return kSTMIN;
  return kUnknownTracking;

}

---