PTSim.cxx

Go to the documentation of this file.

//
// PTSim
//
// Package: PTSim.
//
// S. Kasahara 05/04
//
// Purpose: Define data types and provide general purpose utility routines 
//          to the persistency package.
//

#include <TString.h>

#include "ParticleTransportSim/PTSim.h"
#include "MessageService/MsgService.h"

CVSID("$Id: PTSim.cxx,v 1.7 2007/04/07 04:05:23 schubert Exp $");

//_____________________________________________________________________________
const char* PTSim::AsString(EStdHepSelectMode mode) {
  //  Purpose:  Convert enumerated stdhep select mode to a string.
  //
  //  Return:  mode string. If unknown, returns "Unknown".
  //

  switch ( mode ) {

    case kUnknown:
      return "Unknown";
    case kMomentum:
      return "Momentum";
    case kHit:
      return "Hit";
    default:
      MSG("PTSim",Msg::kWarning) 
      << "PTSim::AsString called with unknown EStdHepSelectMode" 
      << (int)mode << endl;
      return "Unknown";

  }//end of switch

}

//_____________________________________________________________________________
int PTSim::GetStdHepSelectMode(const char* mode) {
  // Purpose: Convert text string cut to an enumerated mode.
  //
  // Return: returns -1 if no match, else EStdHepSelectMode
  //

  TString tmpstr(mode);
  tmpstr.ToLower();
  if ( strcmp(tmpstr.Data(),"momentum") == 0 ) return PTSim::kMomentum;
  else if ( strcmp(tmpstr.Data(),"hit") == 0 ) return PTSim::kHit;

  return -1;

}

//_____________________________________________________________________________
const char* PTSim::AsString(ECut cut) {
  //  Purpose:  Convert enumerated mc cut to a cut string.
  //
  //  Argument: cut: enumerated ECut described in PTSim.h.
  //
  //  Return:  cut string. If unknown, returns "Unknown".
  //
  //  Contact:   S. Kasahara
  //

  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";

    default:
      MSG("PTSim",Msg::kWarning) 
      << "PTSim::AsString called with unknown cut" 
      << (int)cut << endl;
      return "Unknown";
  }//end of switch

}

//_____________________________________________________________________________
int PTSim::GetCut(const char* cut) {
  // Purpose: Convert text string cut to an enumerated code.
  // 
  // Argument: cut string
  //
  // Return: returns -1 if no match
  //
  // Contact: S. Kasahara
  //

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

  return -1;

}

//_____________________________________________________________________________
const char* PTSim::AsString(EProcess process) {
  //  Purpose:  Convert enumerated mc process to a process string.
  //
  //  Argument: process: enumerated EProcess described in PTSim.h.
  //
  //  Return:  process string. If unknown, returns "Unknown".
  //
  //  Contact:   S. Kasahara
  //

  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";

    default:
      MSG("PTSim",Msg::kWarning) 
      << "PTSim::AsString called with unknown process" << (int)process << endl;
      return "Unknown";
  }//end of switch

}

//_____________________________________________________________________________
int PTSim::GetProcess(const char* process) {
  // Purpose: Convert text string process to an enumerated code.
  // 
  // Argument: process string
  //
  // Return: returns -1 if no match
  //
  // Contact: S. Kasahara
  //

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

  return -1;

}

---