ParticlePropertySvc.h

Go to the documentation of this file.

#ifndef PARTICLEPROPERTYSVC_PARTICLEPROPERTYSVC_H
#define PARTICLEPROPERTYSVC_PARTICLEPROPERTYSVC_H
// ============================================================================
// Include Files
// ============================================================================
// STD&STL
// ============================================================================
#include <map>
#include <set>
// ============================================================================
// GaudiKernel
// ============================================================================
#include "GaudiKernel/Service.h"
#include "GaudiKernel/IParticlePropertySvc.h"

class IFileAccess;

namespace Gaudi {
// ============================================================================
class ParticlePropertySvc: public extends<Service,
                                          IParticlePropertySvc>
{
public:
  // typedefs for container type
  typedef ParticleProperty* mapped_type;
  typedef std::map< int,         mapped_type > MapID       ;
  typedef std::map< std::string, mapped_type > MapName     ;
  typedef std::map< int,         mapped_type > MapStdHepID ;
  typedef std::map< int,         mapped_type > MapPythiaID ;
  using IParticlePropertySvc::VectPP;
  using IParticlePropertySvc::const_iterator;
  using IParticlePropertySvc::iterator;
  // Inherited Service overrides:

  StatusCode initialize() override;
  StatusCode finalize() override;
  StatusCode push_back
  ( const std::string& particle ,
    int                geantId  ,
    int                jetsetId ,
    double             charge   ,
    double             mass     ,
    double             tlife    ,
    const std::string& evtName  ,
    int                pythiaId ,
    double             maxWidth ) override ;
  StatusCode push_back( ParticleProperty* pp ) override;
  const_iterator begin() const override { return m_vectpp.begin() ; }
  const_iterator end  () const override { return m_vectpp.end() ; }
  int size() const override { return m_vectpp.size() ; };
  ParticleProperty* find( int geantId ) override
  { return m_idmap[ geantId ] ; }
  ParticleProperty* find( const std::string& name ) override
  { return m_namemap[ name ] ; }
  ParticleProperty* findByStdHepID( int stdHepId ) override
  { return m_stdhepidmap[ stdHepId ] ; }
  ParticleProperty* findByPythiaID( int pythiaID ) override
  { return m_pythiaidmap[ pythiaID ]; }
  StatusCode erase( int geantId ) override
  { return erase ( find ( geantId ) ) ; }
  StatusCode erase( const std::string& name ) override
  { return erase ( find ( name ) ) ; }
  StatusCode eraseByStdHepID( int stdHepId ) override
  { return erase( findByStdHepID ( stdHepId ) ) ; }
  ParticlePropertySvc ( const std::string& name , ISvcLocator* svc );
  ~ParticlePropertySvc() override = default;
protected:
  const ParticleProperty* anti
  ( const ParticleProperty* pp ) const ;
  StatusCode setAntiParticles() ;
  StatusCode rebuild() ;
  StatusCode erase ( const ParticleProperty* pp ) ;
  StatusCode parse ( const std::string& file ) ;
  // Text file parsing routine.
  StatusCode parse();
  // treat additional particles
  StatusCode addParticles() ;
  // compare 2 entries
  bool diff
  ( const ParticleProperty* o              ,
    const ParticleProperty* n              ,
    const MSG::Level l = MSG::DEBUG ) const ;
private:
  typedef std::vector<std::string> Files ;
  typedef std::vector<std::string> Particles ;

  std::string m_filename = "ParticleTable.txt"  ; 
  Files m_other ; 
  // properties to be redefined explicitly
  Particles m_particles ;

  VectPP      m_vectpp;          
  MapID       m_idmap;           
  MapName     m_namemap;         
  MapStdHepID m_stdhepidmap;     
  MapPythiaID m_pythiaidmap;     

  // local storage of ALL properties
  std::set<std::unique_ptr<ParticleProperty>>   m_owned    ;
  std::set<std::string>         m_replaced ;

  SmartIF<IFileAccess> m_fileAccess;
};
} // namespace Gaudi
// =============================================================================
// The END
// =============================================================================
#endif
// =============================================================================