EvtPDL.hh

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
// This software is part of the EvtGen package developed jointly
// for the BaBar and CLEO collaborations.  If you use all or part
// of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
// Copyright (C) 1998 Caltech, UCSB
//
// Module: EvtGen/EvtPDL.hh
//
// Description:Class to keep track of particle properties.
//
// Modification history:
//
// DJL/RYD September 25, 1996 Module created
//
//------------------------------------------------------------------------
 
#ifndef EVTPDL_HH
#define EVTPDL_HH
 
#include "EvtId.hh"
#include "EvtPartProp.hh"
#include "EvtSpinType.hh"
#include "EvtStringHash.hh"
#include <map>
#include <vector>
 
const int SPIN_NAME_LENGTH = 100;
 
class EvtPDL {
 
public:
  EvtPDL();
 
  ~EvtPDL();
 
  void read( const char* fname );
  void readPDT( const std::string fname );
 
  static double getMeanMass( EvtId i ) { return partlist()[i.getId()].getMass(); }
  static double getMass( EvtId i ) { return partlist()[i.getId()].rollMass(); }
  static double getRandMass( EvtId i, EvtId* parId, int nDaug, EvtId* dauId, EvtId* othDaugId,
                             double maxMass, double* dauMasses ) {
    return partlist()[i.getId()].getRandMass( parId, nDaug, dauId, othDaugId, maxMass,
                                              dauMasses );
  }
  static double getMassProb( EvtId i, double mass, double massPar, int nDaug,
                             double* massDau ) {
    return partlist()[i.getId()].getMassProb( mass, massPar, nDaug, massDau );
  }
 
  static double getMaxMass( EvtId i ) { return partlist()[i.getId()].getMassMax(); }
  static double getMinMass( EvtId i ) { return partlist()[i.getId()].getMassMin(); }
  // the number we got from PDT
  static double getMaxRange( EvtId i ) { return partlist()[i.getId()].getMaxRange(); }
  static double getWidth( EvtId i ) { return partlist()[i.getId()].getWidth(); }
  static double getctau( EvtId i ) { return partlist()[i.getId()].getctau(); }
  static int    getStdHep( EvtId id ) { return partlist()[id.getId()].getStdHep(); }
  static int    getLundKC( EvtId id ) { return partlist()[id.getId()].getLundKC(); }
  static EvtId  evtIdFromStdHep( int stdhep );
  static EvtId  chargeConj( EvtId id );
  static int    chg3( EvtId i ) { return partlist()[i.getId()].getChg3(); }
  static EvtSpinType::spintype getSpinType( EvtId i ) {
    return partlist()[i.getId()].getSpinType();
  }
  static EvtId       getId( const std::string& name );
  static std::string name( EvtId i ) { return partlist()[i.getAlias()].getName(); }
  static void        alias( EvtId num, const std::string& newname );
  static void        aliasChgConj( EvtId a, EvtId abar );
  static int         entries() { return partlist().size(); }
  static void reSetMass( EvtId i, double mass ) { partlist()[i.getId()].reSetMass( mass ); }
  static void reSetWidth( EvtId i, double width ) {
    partlist()[i.getId()].reSetWidth( width );
  }
  static void reSetMassMin( EvtId i, double mass ) {
    partlist()[i.getId()].reSetMassMin( mass );
  }
  static void reSetMassMax( EvtId i, double mass ) {
    partlist()[i.getId()].reSetMassMax( mass );
  }
  static void reSetBlatt( EvtId i, double blatt ) {
    partlist()[i.getId()].reSetBlatt( blatt );
  }
  static void includeBirthFactor( EvtId i, bool yesno ) {
    partlist()[i.getId()].includeBirthFactor( yesno );
  }
  static void includeDecayFactor( EvtId i, bool yesno ) {
    partlist()[i.getId()].includeDecayFactor( yesno );
  }
  static void addFactorPn( EvtId i, double factor ) {
    partlist()[i.getId()].addFactorPn( factor );
  }
  static void changeLS( EvtId i, std::string& newLS ) {
    partlist()[i.getId()].newLineShape( newLS );
  }
  static void setPWForDecay( EvtId i, int spin, EvtId d1, EvtId d2 ) {
    partlist()[i.getId()].setPWForDecay( spin, d1, d2 );
  }
  static void setPWForBirthL( EvtId i, int spin, EvtId par, EvtId othD );
  static void fixLSForSP8( EvtId i ) { partlist()[i.getId()].fixLSForSP8(); }
 
private:
  void setUpConstsPdt();
 
  static int _firstAlias;
  static int _nentries;
 
  static std::vector<EvtPartProp>& partlist() {
    static std::vector<EvtPartProp> s_partlist;
    return s_partlist;
  }
 
  static std::map<std::string, int> _particleNameLookup;
 
}; // EvtPDL.h
 
#endif