EvtManyDeltaFuncLineShape.cc

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: EvtLineShape.cc
//
// Description: Store particle properties for one particle.
//
// Modification history:
//
//    Lange      March 10, 2001        Module created
//    Dvoretskii June  03, 2002        Reimplemented rollMass()
//
//------------------------------------------------------------------------
#include "EvtPatches.hh"
 
#include "EvtBlattWeisskopf.hh"
#include "EvtManyDeltaFuncLineShape.hh"
#include "EvtPDL.hh"
#include "EvtPatches.hh"
#include "EvtPropBreitWigner.hh"
#include "EvtPropBreitWignerRel.hh"
#include "EvtRandom.hh"
#include "EvtSpinType.hh"
#include "EvtTwoBodyVertex.hh"
 
EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape() {}
 
EvtManyDeltaFuncLineShape::~EvtManyDeltaFuncLineShape() {}
 
EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape( double mass, double width,
                                                      double                maxRange,
                                                      EvtSpinType::spintype sp ) {
 
  _mass     = mass;
  _width    = width;
  _spin     = sp;
  _maxRange = maxRange;
 
  double maxdelta = width;
 
  _massMax = mass + maxdelta;
  _massMin = mass - maxdelta;
 
  if ( _massMin < 0. ) _massMin = 0.;
}
 
EvtManyDeltaFuncLineShape::EvtManyDeltaFuncLineShape( const EvtManyDeltaFuncLineShape& x ) {
  _mass     = x._mass;
  _width    = x._width;
  _spin     = x._spin;
  _massMax  = x._massMax;
  _massMin  = x._massMin;
  _maxRange = x._maxRange;
}
 
EvtManyDeltaFuncLineShape&
EvtManyDeltaFuncLineShape::operator=( const EvtManyDeltaFuncLineShape& x ) {
  _mass     = x._mass;
  _massMax  = x._massMax;
  _massMin  = x._massMin;
  _width    = x._width;
  _maxRange = x._maxRange;
  _spin     = x._spin;
  return *this;
}
 
EvtAbsLineShape* EvtManyDeltaFuncLineShape::clone() {
 
  return new EvtManyDeltaFuncLineShape( *this );
}
 
double EvtManyDeltaFuncLineShape::getMassProb( double mass, double massPar, int nDaug,
                                               double* massDau ) {
 
  double dTotMass = 0.;
 
  int i;
  for ( i = 0; i < nDaug; i++ ) { dTotMass += massDau[i]; }
  if ( ( mass < dTotMass ) ) return 0.;
 
  if ( massPar > 0.0000000001 )
  {
    if ( mass > massPar ) return 0.;
  }
 
  return 1.;
}
 
double EvtManyDeltaFuncLineShape::getRandMass( EvtId* parId, int nDaug, EvtId* dauId,
                                               EvtId* othDaugId, double maxMass,
                                               double* dauMasses ) {
 
  int nDelta = int( ( _massMax - _massMin ) / _width );
  nDelta++;
  double rand  = EvtRandom::Flat( 0., float( nDelta ) );
  int    randI = int( rand );
  return _massMin + randI * _width;
}