EvtPto3PAmpFactory Class Reference

#include <EvtPto3PAmpFactory.hh>

Inheritance diagram for EvtPto3PAmpFactory:

Public Member Functions
	EvtPto3PAmpFactory (const EvtDalitzPlot &dp)
	EvtPto3PAmpFactory (const EvtPto3PAmpFactory &other)
virtual	~EvtPto3PAmpFactory ()
virtual EvtAmpFactory< EvtDalitzPoint > *	clone () const
virtual void	processAmp (EvtComplex c, std::vector< std::string > vv, bool conj)
Public Member Functions inherited from EvtAmpFactory< EvtDalitzPoint >
	EvtAmpFactory ()
virtual	~EvtAmpFactory ()
virtual void	build (const EvtMultiChannelParser &parser, int nItg)
bool	isCPModel () const
double	dm () const
void	setVerbose ()
EvtAmplitudeSum< EvtDalitzPoint > *	getAmp () const
EvtAmplitudeSum< EvtDalitzPoint > *	getAmpConj () const
EvtPdfSum< EvtDalitzPoint > *	getPC () const
const char *	compName (int i) const
EvtComplex	getCoeff (int i) const
double	getTermCoeff (int i) const

Additional Inherited Members
Protected Attributes inherited from EvtAmpFactory< EvtDalitzPoint >
EvtAmplitudeSum< EvtDalitzPoint > *	_amp
EvtAmplitudeSum< EvtDalitzPoint > *	_ampConj
EvtPdfSum< EvtDalitzPoint > *	_pc
std::vector< std::string >	_names
double	_dm
bool	_verbose

Detailed Description

Definition at line 27 of file EvtPto3PAmpFactory.hh.

Constructor & Destructor Documentation

EvtPto3PAmpFactory() [1/2]

EvtPto3PAmpFactory::EvtPto3PAmpFactory ( const EvtDalitzPlot & dp )

inline

Definition at line 30 of file EvtPto3PAmpFactory.hh.

: EvtAmpFactory<EvtDalitzPoint>(), _dp( dp ) {}

Referenced by clone(), and EvtPto3PAmpFactory().

EvtPto3PAmpFactory() [2/2]

EvtPto3PAmpFactory::EvtPto3PAmpFactory ( const EvtPto3PAmpFactory & other )

inline

Definition at line 31 of file EvtPto3PAmpFactory.hh.

      : EvtAmpFactory<EvtDalitzPoint>( other ), _dp( other._dp ) {}

~EvtPto3PAmpFactory()

virtual EvtPto3PAmpFactory::~EvtPto3PAmpFactory ( )

inlinevirtual

Definition at line 33 of file EvtPto3PAmpFactory.hh.

{}

Member Function Documentation

clone()

virtual EvtAmpFactory< EvtDalitzPoint > * EvtPto3PAmpFactory::clone ( ) const

inlinevirtual

Implements EvtAmpFactory< EvtDalitzPoint >.

Definition at line 35 of file EvtPto3PAmpFactory.hh.

                                                       {
    return new EvtPto3PAmpFactory( *this );
  }

processAmp()

void EvtPto3PAmpFactory::processAmp ( EvtComplex c, std::vector< std::string > vv, bool conj )

virtual

Implements EvtAmpFactory< EvtDalitzPoint >.

Definition at line 42 of file EvtPto3PAmpFactory.cc.

                                                                                      {
  if ( _verbose )
  {
 
    printf( "Make amplitude\n" );
    unsigned i;
    for ( i = 0; i < vv.size(); i++ ) printf( "%s\n", vv[i].c_str() );
    printf( "\n" );
  }
 
  EvtAmplitude<EvtDalitzPoint>* amp = 0;
  EvtPdf<EvtDalitzPoint>*       pdf = 0;
  std::string                   name;
 
  int i;
  if ( vv[0] == "PHASESPACE" )
  {
 
    pdf  = new EvtDalitzFlatPdf( _dp );
    amp  = new EvtFlatAmp<EvtDalitzPoint>();
    name = "NR";
  }
  else if ( vv[0] == "RESONANCE" )
  {
 
    EvtPto3PAmp* partAmp = 0;
 
    // RESONANCE stanza
 
    Pair                  pairRes = strToPair( vv[1].c_str() );
    EvtSpinType::spintype spinR;
    double                mR, gR;
    name        = vv[2];
    EvtId resId = EvtPDL::getId( vv[2] );
    if ( _verbose ) printf( "Particles %s form resonance %s\n", vv[1].c_str(), vv[2].c_str() );
 
    // If no valid particle name is given, assume that
    // it is the spin, the mass and the width of the particle.
 
    if ( resId.getId() == -1 )
    {
 
      switch ( atoi( vv[2].c_str() ) )
      {
 
      case 0: {
        spinR = EvtSpinType::SCALAR;
        break;
      }
      case 1: {
        spinR = EvtSpinType::VECTOR;
        break;
      }
      case 2: {
        spinR = EvtSpinType::TENSOR;
        break;
      }
      case 3: {
        spinR = EvtSpinType::SPIN3;
        break;
      }
      case 4: {
        spinR = EvtSpinType::SPIN4;
        break;
      }
      default: {
        assert( 0 );
        break;
      }
      }
 
      mR = strtod( vv[3].c_str(), 0 );
      gR = strtod( vv[4].c_str(), 0 );
      i  = 4;
    }
    else
    {
 
      // For a valid particle get spin, mass and width
 
      spinR = EvtPDL::getSpinType( resId );
      mR    = EvtPDL::getMeanMass( resId );
      gR    = EvtPDL::getWidth( resId );
      i     = 2;
 
      // It's possible to specify mass and width of a particle
      // explicitly
 
      if ( vv[3] != "ANGULAR" )
      {
 
        if ( _verbose )
          printf( "Setting m(%s)=%s g(%s)=%s\n", vv[2].c_str(), vv[3].c_str(), vv[2].c_str(),
                  vv[4].c_str() );
 
        mR = strtod( vv[3].c_str(), 0 );
        gR = strtod( vv[4].c_str(), 0 );
        i  = 4;
      }
    }
 
    // ANGULAR stanza
 
    if ( vv[++i] != "ANGULAR" )
    {
 
      printf( "%s instead of ANGULAR\n", vv[i].c_str() );
      exit( 0 );
    }
    Pair pairAng = strToPair( vv[++i].c_str() );
    if ( _verbose ) printf( "Angle is measured between particles %s\n", vv[i].c_str() );
 
    // TYPE stanza
 
    assert( vv[++i] == "TYPE" );
    std::string type = vv[++i];
    if ( _verbose ) printf( "Propagator type %s\n", vv[i].c_str() );
 
    if ( type == "NBW" )
    {
 
      EvtPropBreitWigner prop( mR, gR );
      partAmp = new EvtPto3PAmp( _dp, pairAng, pairRes, spinR, prop, EvtPto3PAmp::NBW );
    }
    else if ( type == "RBW_ZEMACH" )
    {
 
      EvtPropBreitWignerRel prop( mR, gR );
      partAmp = new EvtPto3PAmp( _dp, pairAng, pairRes, spinR, prop, EvtPto3PAmp::RBW_ZEMACH );
    }
    else if ( type == "RBW_KUEHN" )
    {
 
      EvtPropBreitWignerRel prop( mR, gR );
      partAmp = new EvtPto3PAmp( _dp, pairAng, pairRes, spinR, prop, EvtPto3PAmp::RBW_KUEHN );
    }
    else if ( type == "RBW_CLEO" )
    {
 
      EvtPropBreitWignerRel prop( mR, gR );
      partAmp = new EvtPto3PAmp( _dp, pairAng, pairRes, spinR, prop, EvtPto3PAmp::RBW_CLEO );
    }
    else assert( 0 );
 
    // Optional DVFF, BVFF stanzas
 
    if ( i < vv.size() - 1 )
    {
      if ( vv[i + 1] == "DVFF" )
      {
        i++;
        if ( vv[++i] == "BLATTWEISSKOPF" )
        {
 
          double R = strtod( vv[++i].c_str(), 0 );
          partAmp->set_fd( R );
        }
        else assert( 0 );
      }
    }
 
    if ( i < vv.size() - 1 )
    {
      if ( vv[i + 1] == "BVFF" )
      {
        i++;
        if ( vv[++i] == "BLATTWEISSKOPF" )
        {
 
          if ( _verbose ) printf( "BVFF=%s\n", vv[i].c_str() );
          double R = strtod( vv[++i].c_str(), 0 );
          partAmp->set_fb( R );
        }
        else assert( 0 );
      }
    }
 
    i++;
 
    pdf = new EvtDalitzResPdf( _dp, mR, gR, pairRes );
    amp = partAmp;
  }
 
  assert( amp );
  assert( pdf );
 
  if ( !conj )
  {
 
    _amp->addOwnedTerm( c, amp );
    _pc->addOwnedTerm( abs2( c ), pdf );
  }
  else
  {
 
    _ampConj->addOwnedTerm( c, amp );
    delete pdf;
  }
  _names.push_back( name );
}

---

The documentation for this class was generated from the following files:

• EvtPto3PAmpFactory.hh
• EvtPto3PAmpFactory.cc