EvtFlatLineShape.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 "EvtFlatLineShape.hh"
#include "EvtPDL.hh"
#include "EvtPatches.hh"
#include "EvtPropBreitWigner.hh"
#include "EvtPropBreitWignerRel.hh"
#include "EvtRandom.hh"
#include "EvtSpinType.hh"
#include "EvtTwoBodyVertex.hh"
 
EvtFlatLineShape::EvtFlatLineShape() {}
 
EvtFlatLineShape::~EvtFlatLineShape() {}
 
EvtFlatLineShape::EvtFlatLineShape( 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.;
}
 
EvtFlatLineShape::EvtFlatLineShape( const EvtFlatLineShape& x ) {
  _mass     = x._mass;
  _width    = x._width;
  _spin     = x._spin;
  _massMax  = x._massMax;
  _massMin  = x._massMin;
  _maxRange = x._maxRange;
}
 
EvtFlatLineShape& EvtFlatLineShape::operator=( const EvtFlatLineShape& x ) {
  _mass     = x._mass;
  _massMax  = x._massMax;
  _massMin  = x._massMin;
  _width    = x._width;
  _maxRange = x._maxRange;
  _spin     = x._spin;
  return *this;
}
 
EvtAbsLineShape* EvtFlatLineShape::clone() { return new EvtFlatLineShape( *this ); }
 
double EvtFlatLineShape::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 EvtFlatLineShape::getRandMass( EvtId* parId, int nDaug, EvtId* dauId, EvtId* othDaugId,
                                      double massMax, double* dauMasses ) {
 
  return EvtRandom::Flat( _massMin, _massMax );
}