/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Generator/BesEvtGen/BesEvtGen-00-03-58/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

---