EvtAbsLineShape.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/EvtAbsLineShape.hh
//
// Description: Class to keep the particle properties for
//              one particle
//
// Modification history:
//
//    Lange     March 10, 2001         Module created
//
//------------------------------------------------------------------------
 
#ifndef EVTABSLINESHAPE_HH
#define EVTABSLINESHAPE_HH
 
#include "EvtId.hh"
#include "EvtSpinType.hh"
#include <vector>
 
class EvtId;
 
class EvtAbsLineShape {
 
public:
  EvtAbsLineShape();
  EvtAbsLineShape( double mass, double width, double maxRange, EvtSpinType::spintype sp );
  virtual ~EvtAbsLineShape();
  EvtAbsLineShape& operator=( const EvtAbsLineShape& x );
  EvtAbsLineShape( const EvtAbsLineShape& x );
 
  double                   getMass() { return _mass; }
  double                   getMassMin() { return _massMin; }
  double                   getMassMax() { return _massMax; }
  double                   getMaxRange() { return _maxRange; }
  double                   getWidth() { return _width; }
  EvtSpinType::spintype    getSpinType() { return _spin; }
  virtual double           rollMass();
  virtual EvtAbsLineShape* clone();
 
  void reSetMass( double mass ) { _mass = mass; }
  void reSetWidth( double width ) {
    _width = width;
    // <--- added by L.L. Wang to fix a bug
    double maxdelta = 15.0 * width;
    if ( _maxRange > 0.00001 )
    {
      _massMax = _mass + maxdelta;
      _massMin = _mass - _maxRange;
    }
    else
    {
      _massMax = _mass + maxdelta;
      _massMin = _mass - 15.0 * width;
    }
    if ( _massMin < 0. ) _massMin = 0.;
    _massMax = _mass + maxdelta;
    // ---> //
  }
  void         reSetMassMin( double mass ) { _massMin = mass; }
  void         reSetMassMax( double mass ) { _massMax = mass; }
  virtual void reSetBlatt( double blatt ){};
  void         includeBirthFactor( bool yesno ) { _includeBirthFact = yesno; }
  void         addFactorPn( double factor = 0. ) { _addFactorPn = factor; }
  void         includeDecayFactor( bool yesno ) { _includeDecayFact = yesno; }
  void         setPWForDecay( int spin, EvtId d1, EvtId d2 ) {
    _userSetPW.push_back( spin );
    _userSetPWD1.push_back( d1 );
    _userSetPWD2.push_back( d2 );
  }
  void setPWForBirthL( int spin, EvtId par, EvtId othD ) {
    _userSetBirthPW.push_back( spin );
    _userSetBirthOthD.push_back( othD );
    _userSetBirthPar.push_back( par );
  }
 
  virtual double getRandMass( EvtId* parId, int nDaug, EvtId* dauId, EvtId* othDaugId,
                              double maxMass, double* dauMasses );
  virtual double getMassProb( double mass, double massPar, int nDaug, double* massDau );
 
  void fixForSP8() { _applyFixForSP8 = true; }
 
protected:
  bool   _includeDecayFact;
  bool   _includeBirthFact;
  double _addFactorPn;
  double _mass;
  double _massMin;
  double _massMax;
  double _width;
  double _maxRange;
 
  // allow for special cases where the default method of picking the
  // lowest allowed partial wave for a decay is not the right answer.
  // string is "<spin> <daughter1> <daughter2>"
  // new 9/12/2003 Lange
  std::vector<EvtId> _userSetPWD1, _userSetPWD2;
  std::vector<int>   _userSetPW;
 
  // also do it for birth factors
  std::vector<EvtId> _userSetBirthPar, _userSetBirthOthD;
  std::vector<int>   _userSetBirthPW;
 
  EvtSpinType::spintype _spin;
 
  bool _applyFixForSP8;
};
 
#endif