LoadMinosPDG.cxx

Go to the documentation of this file.

// $Id: LoadMinosPDG.cxx,v 1.2 2009/05/12 22:36:19 rhatcher Exp $
//
// Load PDG database and insert MINOS data
//
// rhatcher@fnal.gov

#include "MessageService/MsgService.h"
CVSID("$Id: LoadMinosPDG.cxx,v 1.2 2009/05/12 22:36:19 rhatcher Exp $");

//#include "TString.h"
#include "TSystem.h"
#include "TEnv.h"
#include "TDatabasePDG.h"
#include "TParticle.h"


void LoadMinosPDG(bool onlyRootVersion, bool fake)
{

  static bool first = true;
  if (!first) return;
  first = false;

  if (fake) {
    MSG("DataUtil", Msg::kInfo) 
      << "LoadMinosPDG fake: NOT initializing PDG database " 
      << ", assume user will do so via alternative method (e.g. GENIE)"
      << endl;
    return;
  }

  MSG("DataUtil", Msg::kInfo) 
    << "LoadMinosPDG: initializing PDG database " 
    << (onlyRootVersion ? "(root version only)" : "(minos version)" )
    << endl;

  TDatabasePDG *pdg = TDatabasePDG::Instance();

  // read default $ROOTSYS/etc/pdg_table.txt 
  // (unless overridden by Root.DatabasePDG)
  // Only do it once, since it crashes the program the second time. NJT Sept 03

  const char* minos_search_path =
    ".:$SRT_PRIVATE_CONTEXT/etc:$SRT_PUBLIC_CONTEXT/etc";
  const char* minos_file = "minos_pdg_table.txt";  

  char *minos_pdg = (onlyRootVersion) ? 0 :
    gSystem->Which(minos_search_path,minos_file,kReadPermission);
  if ( minos_pdg ) {
    MSG("DataUtil", Msg::kInfo) 
      << "LoadMinosPDG: reading " << minos_pdg << endl;

    pdg->ReadPDGTable(minos_pdg);
    // gSystem->Which returns a string the user must delete
    delete [] minos_pdg;
  }
  else {
    if (!onlyRootVersion)
      MSG("DataUtil", Msg::kInfo)
        << "No \"" << minos_file << "\" file found " << endl
        << "search path was: " << minos_search_path << endl;
    
    const char* root_search_path =
      ".:$ROOTSYS/etc:/etc/root";
    const char* root_file = "pdg_table.txt";

    char *root_pdg =
      gSystem->Which(root_search_path,root_file,kReadPermission);
    if ( root_pdg ) {

      MSG("DataUtil", Msg::kInfo) 
        << "LoadMinosPDG: reading " << root_pdg << endl;

      pdg->ReadPDGTable(root_pdg);
      // gSystem->Which returns a string the user must delete
      delete [] root_pdg;
    }
    else {
      MSG("DataUtil", Msg::kInfo)
        << "No \"" << root_file << "\" file found " << endl
        << "search path was: " << root_search_path << endl;
    }
  }
}

---