PartPropSvc Namespace Reference

Functions
	PartPropSvc (const std::string &name, ISvcLocator *svc)
	~PartPropSvc ()
StatusCode	initialize ()
StatusCode	finalize ()

---

Function Documentation

StatusCode PartPropSvc::finalize

(

)

Definition at line 137 of file PartPropSvc.cpp.

                                 {

  MsgStream log( msgSvc(), name() );
  StatusCode status = Service::finalize();

  if ( status.isSuccess() )
    log << MSG::INFO << "Service finalised successfully" << endmsg;

  return status;
}

StatusCode PartPropSvc::initialize

(

)

Definition at line 44 of file PartPropSvc.cpp.

                                   {

  MsgStream log( msgSvc(), name() );
  std::vector<std::string>::const_iterator itr;

  StatusCode status = Service::initialize();

  log << MSG::INFO << "PDT file(s): " << endmsg;
  for (itr=m_pdtFiles.begin(); itr!=m_pdtFiles.end(); ++itr) {
    log << MSG::INFO << "    " << *itr << endmsg;
  }
  log << MSG::INFO << "Type:     " << m_inputType << endmsg;

  if ( status.isFailure() ) {
    log << MSG::ERROR << "Could not initialize main svc" << endmsg;
    return StatusCode::FAILURE;
  }

  bool (*pF)  (std::istream &,
               HepPDT::TableBuilder &);


  // Determine type of input
  if (m_inputType == "PDG") {
    pF = &HepPDT::addPDGParticles;
  } else if (m_inputType == "Pythia") {
    pF = &HepPDT::addPythiaParticles;
  } else if (m_inputType == "EvtGen") {
    pF = &HepPDT::addEvtGenParticles;
  } else if (m_inputType == "Herwig") {
    pF = &HepPDT::addHerwigParticles;
  } else if (m_inputType == "IsaJet") {
    pF = &HepPDT::addIsajetParticles;
  } else if (m_inputType == "QQ") {
    pF = &HepPDT::addQQParticles;
  } else {
    log << MSG::ERROR << "Unknown Particle Data file type: \""
        << m_inputType << "\"" << endmsg;
    return StatusCode::FAILURE;
  }

  // Make sure we have at least one file
  if (m_pdtFiles.size() == 0) {
    log << MSG::ERROR << "Need at least 1 PDT file" << endmsg;
    log << MSG::ERROR << "Check value of property \"InputFile\"" << endmsg;
    return StatusCode::FAILURE;
  }

  m_pdt = new HepPDT::ParticleDataTable;

  {
    // Construct table builder
    HepPDT::TableBuilder  tb( *m_pdt );

    // read the input
    int good(0);
    for (itr=m_pdtFiles.begin(); itr!=m_pdtFiles.end(); ++itr) {

      std::string rfile = System::PathResolver::find_file(*itr,"DATAPATH");
      if (rfile == "") {
        log << MSG::ERROR << "Could not find PDT file: \"" << *itr
            << "\" in $DATAPATH" << endmsg;
        continue;
      }

      std::ifstream pdfile( rfile.c_str() );
      if (!pdfile) {
        log << MSG::ERROR << "Could not open PDT file: \"" << rfile
            << "\"" << endmsg;
        continue;
      }

      if ( ! pF(pdfile,tb) ) {
        log << MSG::ERROR << "Error reading PDT file: \"" << rfile
            << "\"" << endmsg;
        return StatusCode::FAILURE;
      }
      ++good;
    }
    if (0 == good) {
      log << MSG::ERROR << "Unable to access any PDT file" <<endmsg;
      return StatusCode::FAILURE;
    }

  }   // the tb destructor fills datacol



  return status;
}

PartPropSvc::PartPropSvc	(	const std::string &	name,
		ISvcLocator *	svc
	)

Definition at line 21 of file PartPropSvc.cpp.

  : base_class( name, svc ), m_pdt(0) {

  declareProperty( "InputType", m_inputType="PDG");
  declareProperty( "InputFile", m_pdtFiles);

  if (m_pdtFiles.empty() ) {
    m_pdtFiles.push_back("PDGTABLE.MeV");
  }

}

PartPropSvc::~PartPropSvc

(

)

Definition at line 35 of file PartPropSvc.cpp.

                          {
  if (m_pdt != 0) {
    delete m_pdt;
    m_pdt = 0;
  }
}