bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/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 "EvtGenBase/EvtSpinType.hh"
#include "EvtGenBase/EvtId.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