PdtEntry.cxx

Go to the documentation of this file.

//
// PdtEntry.cc - data class for a particle
//
// Copyright (C)  1993  The Board of Trustees of The Leland Stanford
//
// History:
//  Migration for BESIII MDC
//                      Junior University.  All Rights Reserved.
//
// $Id: PdtEntry.cxx,v 1.2 2009/12/23 02:59:56 zhangy Exp $
//
// A PdtEntry holds the Particle Data Table data for a single particle.
//
// Modified:
//    Luca Lista     04 oct 96    lookup by different types, not by integer
//
// See Also
//    Pdt, DecayMode
 
extern "C" {
#include <float.h>
#include <string.h>
}
 
// This is for systems that define the BSD string functions as macros:
#ifdef index
#  undef index
#endif
#ifdef rindex
#  undef rindex
#endif
 
#include "MdcRecoUtil/DecayMode.h"
#include "MdcRecoUtil/Pdt.h"
#include "MdcRecoUtil/PdtEntry.h"
#include <iomanip>
#include <iostream>
#include <vector>
using std::endl;
using std::ios;
using std::ostream;
using std::setw;
using std::vector;
 
#define HBARC ( 197.327 * 1.e-3 * 1.e-13 ) // GeV*cm
 
PdtEntry::PdtEntry( const char* name, PdtLund::LundType code, float spin, float charge,
                    float mass, float width, float massCut )
    : _name( new char[strlen( name ) + 1] )
    , _mass( mass )
    , _width( width )
    , _lifetime( ( width > 0.0 ) ? ( HBARC / width ) : FLT_MAX )
    , _spin( spin )
    , _charge( charge )
    , _widthCut( ( massCut > 0.0 ) ? massCut : 2.0 * width )
    , _sumBR( 0.0 )
    , _decayList( new vector<DecayMode*> )
    , _lundId( code )
    , _pdgId( Pdt::pdgId( code ) )
    , _geantId( Pdt::geantId( code ) )
    , _pidId( Pdt::pidId( code ) )
    , _pidNeutId( Pdt::pidNeutId( code ) )
    , _conjugate( 0 ) {
  strcpy( _name, name );
}
 
PdtEntry::PdtEntry( const char* name, PdtGeant::GeantType code, float spin, float charge,
                    float mass, float width, float massCut )
    : _name( new char[strlen( name ) + 1] )
    , _mass( mass )
    , _width( width )
    , _lifetime( ( width > 0.0 ) ? ( HBARC / width ) : FLT_MAX )
    , _spin( spin )
    , _charge( charge )
    , _widthCut( ( massCut > 0.0 ) ? massCut : 2.0 * width )
    , _sumBR( 0.0 )
    , _decayList( new vector<DecayMode*> )
    , _lundId( Pdt::lundId( code ) )
    , _pdgId( Pdt::pdgId( code ) )
    , _geantId( code )
    , _pidId( Pdt::pidId( code ) )
    , _pidNeutId( Pdt::pidNeutId( code ) )
    , _conjugate( 0 ) {
  strcpy( _name, name );
}
 
PdtEntry::PdtEntry( const char* name, PdtPdg::PdgType code, float spin, float charge,
                    float mass, float width, float massCut )
    : _name( new char[strlen( name ) + 1] )
    , _mass( mass )
    , _width( width )
    , _lifetime( ( width > 0.0 ) ? ( HBARC / width ) : FLT_MAX )
    , _spin( spin )
    , _charge( charge )
    , _widthCut( ( massCut > 0.0 ) ? massCut : 2.0 * width )
    , _sumBR( 0.0 )
    , _decayList( new vector<DecayMode*> )
    , _lundId( Pdt::lundId( code ) )
    , _pdgId( code )
    , _geantId( Pdt::geantId( code ) )
    , _pidId( Pdt::pidId( code ) )
    , _pidNeutId( Pdt::pidNeutId( code ) )
    , _conjugate( 0 ) {
  strcpy( _name, name );
}
 
void PdtEntry::addDecay( float bf, vector<PdtEntry*>* kids ) {
  _decayList->push_back( new DecayMode( bf, kids ) );
  _sumBR += bf;
}
 
#ifndef HP1022
void PdtEntry::printOn( ostream& str ) const {
  // pkg coordinator:  we need to see if form() really is part of
  // iostream.  In any case, i've recoded things below simply to
  // avoid it.  Perhaps this is good enough.
  //
 
  //   str.form("%6d %-8s %8g ", (int) _lundId, _name, _mass);
  //   Not sure if the is exactly the same as using form(), but on
  //   sun, it's claimed that you get %g behavior from ios as default.
  //   I'm ignoring the left-field justification for now.
  //
  str << setw( 6 ) << (int)_lundId << " " << setw( 8 ) << _name << " " << setw( 8 ) << _mass
      << " ";
 
  if ( _width > 0.0 ) str << " " << _width;
  else str << "Stable      ";
  str << setw( 3 ) << _spin << "   " << setw( 2 ) << _charge << "   ";
  if ( _width > 0.0 ) str << " " << _lifetime << endl;
  else str << " Stable\n";
}
#else
void PdtEntry::printOn( ostream& str ) const {
  ios::sync_with_stdio();
  printf( "%6d %-8s %8g ", (int)_lundId, _name, _mass );
  if ( _width > 0.0 )
  {
    printf( "%g ", _width );
    str << " ";
  }
  else str << "Stable      ";
 
  printf( "%3g   %2g   ", _spin, _charge );
  if ( _width > 0.0 )
  {
    printf( " %g\n", _lifetime );
    str << " ";
  }
  else str << " Stable\n";
}
#endif
 
void PdtEntry::printBFOn( ostream& str ) const {
  int l = _decayList->size();
  for ( int i = 0; i < l; i++ ) ( *_decayList )[i]->printOn( str );
}
 
PdtEntry::~PdtEntry() {
  vector<DecayMode*>::iterator iter = _decayList->begin();
  while ( iter != _decayList->end() )
  {
    delete *iter;
    ++iter;
  }
  _decayList->clear();
  delete _decayList;
  delete[] _name;
}
 
const PdtEntry* PdtEntry::conjugate() const {
  if ( !_conjugate ) { const_cast<PdtEntry*>( this )->_conjugate = Pdt::conjugate( this ); }
  return _conjugate;
}
 
bool PdtEntry::operator==( const PdtEntry& theOther ) const {
  return ( theOther._lundId == _lundId && theOther._geantId == _geantId );
}
 
bool PdtEntry::operator<( const PdtEntry& theOther ) const {
  return ( _lundId < theOther._lundId ||
           ( _lundId == theOther._lundId && _geantId < theOther._geantId ) );
}