EvtRelBreitWignerBarrierFact Class Reference

#include <EvtRelBreitWignerBarrierFact.hh>

Inheritance diagram for EvtRelBreitWignerBarrierFact:

Public Member Functions
	EvtRelBreitWignerBarrierFact ()
	EvtRelBreitWignerBarrierFact (double mass, double width, double maxRange, EvtSpinType::spintype sp)
	~EvtRelBreitWignerBarrierFact ()
EvtRelBreitWignerBarrierFact &	operator= (const EvtRelBreitWignerBarrierFact &x)
	EvtRelBreitWignerBarrierFact (const EvtRelBreitWignerBarrierFact &x)
EvtAbsLineShape *	clone ()
double	getMassProb (double mass, double massPar, int nDaug, double *massDau)
double	getRandMass (EvtId *parId, int nDaug, EvtId *dauId, EvtId *othDaugId, double maxMass, double *dauMasses)
virtual void	reSetBlatt (double blatt)
Public Member Functions inherited from EvtAbsLineShape
	EvtAbsLineShape ()
	EvtAbsLineShape (double mass, double width, double maxRange, EvtSpinType::spintype sp)
virtual	~EvtAbsLineShape ()
EvtAbsLineShape &	operator= (const EvtAbsLineShape &x)
	EvtAbsLineShape (const EvtAbsLineShape &x)
double	getMass ()
double	getMassMin ()
double	getMassMax ()
double	getMaxRange ()
double	getWidth ()
EvtSpinType::spintype	getSpinType ()
virtual double	rollMass ()
virtual EvtAbsLineShape *	clone ()
void	reSetMass (double mass)
void	reSetWidth (double width)
void	reSetMassMin (double mass)
void	reSetMassMax (double mass)
virtual void	reSetBlatt (double blatt)
void	includeBirthFactor (bool yesno)
void	addFactorPn (double factor=0.)
void	includeDecayFactor (bool yesno)
void	setPWForDecay (int spin, EvtId d1, EvtId d2)
void	setPWForBirthL (int spin, EvtId par, EvtId othD)
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 ()

Protected Attributes
double	_blatt
bool	_errorCond
Protected Attributes inherited from EvtAbsLineShape
bool	_includeDecayFact
bool	_includeBirthFact
double	_addFactorPn
double	_mass
double	_massMin
double	_massMax
double	_width
double	_maxRange
std::vector< EvtId >	_userSetPWD1
std::vector< EvtId >	_userSetPWD2
std::vector< int >	_userSetPW
std::vector< EvtId >	_userSetBirthPar
std::vector< EvtId >	_userSetBirthOthD
std::vector< int >	_userSetBirthPW
EvtSpinType::spintype	_spin
bool	_applyFixForSP8

Detailed Description

Definition at line 27 of file EvtRelBreitWignerBarrierFact.hh.

Constructor & Destructor Documentation

EvtRelBreitWignerBarrierFact() [1/3]

EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact

(

)

Definition at line 36 of file EvtRelBreitWignerBarrierFact.cc.

                                                           {
 
}

Referenced by clone().

EvtRelBreitWignerBarrierFact() [2/3]

EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact	(	double	mass,
		double	width,
		double	maxRange,
		EvtSpinType::spintype	sp
	)

Definition at line 43 of file EvtRelBreitWignerBarrierFact.cc.

                                                                                                                             :
  EvtAbsLineShape(mass,width,maxRange,sp)
{ // double mDaug1, double mDaug2, int l) {
 
  _includeDecayFact=true;
  _includeBirthFact=true;
  _mass=mass;
  _width=width;
  _spin=sp;
  _blatt=3.0;
  _maxRange=maxRange;
  _errorCond=false;
 
  double maxdelta = 15.0*width;
 
  if ( maxRange > 0.00001 ) {
    _massMax=mass+maxdelta;
    _massMin=mass-maxRange;
  }
  else{
    _massMax=mass+maxdelta;
    _massMin=mass-15.0*width;
  }
 
  _massMax=mass+maxdelta;
  if ( _massMin< 0. ) _massMin=0.;
}

~EvtRelBreitWignerBarrierFact()

EvtRelBreitWignerBarrierFact::~EvtRelBreitWignerBarrierFact

(

)

Definition at line 40 of file EvtRelBreitWignerBarrierFact.cc.

                                                            {
}

EvtRelBreitWignerBarrierFact() [3/3]

EvtRelBreitWignerBarrierFact::EvtRelBreitWignerBarrierFact

(

const EvtRelBreitWignerBarrierFact &

x

)

Definition at line 71 of file EvtRelBreitWignerBarrierFact.cc.

                                                                                                :
  EvtAbsLineShape(x)
{
  _massMax=x._massMax;
  _massMin=x._massMin;
  _blatt=x._blatt;
  _maxRange=x._maxRange;
  _includeDecayFact=x._includeDecayFact;
  _includeBirthFact=x._includeBirthFact;
  _errorCond=x._errorCond;
 
}

Member Function Documentation

clone()

EvtAbsLineShape * EvtRelBreitWignerBarrierFact::clone ( )

virtual

Reimplemented from EvtAbsLineShape.

Definition at line 99 of file EvtRelBreitWignerBarrierFact.cc.

                                                     {
 
  return new EvtRelBreitWignerBarrierFact(*this);
}

getMassProb()

double EvtRelBreitWignerBarrierFact::getMassProb ( double  mass, double  massPar, int  nDaug, double *  massDau  )

virtual

Reimplemented from EvtAbsLineShape.

Definition at line 105 of file EvtRelBreitWignerBarrierFact.cc.

                                                                                                       {
 
  _errorCond=false;
  //return EvtAbsLineShape::getMassProb(mass,massPar,nDaug,massDau);
  if (nDaug!=2) return EvtAbsLineShape::getMassProb(mass,massPar,nDaug,massDau);
 
  double dTotMass=0.;
 
  int i;
  for (i=0; i<nDaug; i++) {
    dTotMass+=massDau[i];
  }
  //report(INFO,"EvtGen") << mass << " " << massPar << " " << dTotMass << " "<< endl;
  //    if ( (mass-dTotMass)<0.0001 ) return 0.;
  //report(INFO,"EvtGen") << mass << " " << dTotMass << endl;
  if ( (mass<dTotMass) ) return 0.;
 
  if ( _width< 0.0001) return 1.;
 
  if ( massPar>0.0000000001 ) {
    if ( mass > massPar) return 0.;
  }
 
  if ( _errorCond ) return 0.;
 
  // we did all the work in getRandMass
  return 1.;
}

getRandMass()

double EvtRelBreitWignerBarrierFact::getRandMass ( EvtId *  parId, int  nDaug, EvtId *  dauId, EvtId *  othDaugId, double  maxMass, double *  dauMasses  )

virtual

Reimplemented from EvtAbsLineShape.

Definition at line 134 of file EvtRelBreitWignerBarrierFact.cc.

                                                                                                                                          {
  if ( nDaug!=2) {
    if(fabs(_addFactorPn) >0.00000001) EvtAbsLineShape::addFactorPn(_addFactorPn);
    double themass = EvtAbsLineShape::getRandMass(parId,nDaug,dauId,othDaugId,maxMass,dauMasses);
    return themass;
  }
 
  if ( _width< 0.00001) return _mass;
 
  //first figure out L - take the lowest allowed.
 
  EvtSpinType::spintype spinD1=EvtPDL::getSpinType(dauId[0]);
  EvtSpinType::spintype spinD2=EvtPDL::getSpinType(dauId[1]);
 
  int t1=EvtSpinType::getSpin2(spinD1);
  int t2=EvtSpinType::getSpin2(spinD2);
  int t3=EvtSpinType::getSpin2(_spin);
 
  int Lmin=-10;
 
 
  // the user has overridden the partial wave to use.
  for ( int vC=0; vC<_userSetPW.size(); vC++) {
    if ( dauId[0]==_userSetPWD1[vC] &&  dauId[1]==_userSetPWD2[vC] ) Lmin=2*_userSetPW[vC]; 
    if ( dauId[0]==_userSetPWD2[vC] &&  dauId[1]==_userSetPWD1[vC] ) Lmin=2*_userSetPW[vC]; 
  }
  
  // allow for special cases.
  if (Lmin<-1 ) {
    
    //There are some things I don't know how to deal with
    if ( t3>4) return EvtAbsLineShape::getRandMass(parId,nDaug,dauId,othDaugId,maxMass,dauMasses);
    if ( t1>4) return EvtAbsLineShape::getRandMass(parId,nDaug,dauId,othDaugId,maxMass,dauMasses);
    if ( t2>4) return EvtAbsLineShape::getRandMass(parId,nDaug,dauId,othDaugId,maxMass,dauMasses);
    
    //figure the min and max allowwed "spins" for the daughters state
    Lmin=std::max(t3-t2-t1,std::max(t2-t3-t1,t1-t3-t2));
    if (Lmin<0) Lmin=0;
    assert(Lmin==0||Lmin==2||Lmin==4);
  }
 
  //double massD1=EvtPDL::getMeanMass(dauId[0]);
  //double massD2=EvtPDL::getMeanMass(dauId[1]);
  double massD1=dauMasses[0];
  double massD2=dauMasses[1];
 
  // I'm not sure how to define the vertex factor here - so retreat to nonRel code.
  if ( (massD1+massD2)> _mass ) return  EvtAbsLineShape::getRandMass(parId,nDaug,dauId,othDaugId,maxMass,dauMasses);
 
  //parent vertex factor not yet implemented
  double massOthD=-10.;
  double massParent=-10.;
  int birthl=-10;
  if ( othDaugId) {
    EvtSpinType::spintype spinOth=EvtPDL::getSpinType(*othDaugId);
    EvtSpinType::spintype spinPar=EvtPDL::getSpinType(*parId);
    
    int tt1=EvtSpinType::getSpin2(spinOth);
    int tt2=EvtSpinType::getSpin2(spinPar);
    int tt3=EvtSpinType::getSpin2(_spin);
    
    
    //figure the min and max allowwed "spins" for the daughters state
    if ( (tt1<=4) && ( tt2<=4) ) {
      birthl=std::max(tt3-tt2-tt1,std::max(tt2-tt3-tt1,tt1-tt3-tt2));
      if (birthl<0) birthl=0;
    
      massOthD=EvtPDL::getMeanMass(*othDaugId);
      massParent=EvtPDL::getMeanMass(*parId);
    
    }
 
 
     // allow user to override
     for (size_t vC=0; vC<_userSetBirthPW.size(); vC++) {
       if ( *othDaugId==_userSetBirthOthD[vC] && *parId==_userSetBirthPar[vC]){
    birthl=2*_userSetBirthPW[vC];
       } 
     }
 
  }
  double massM=_massMax;
  if ( (maxMass > -0.5) && (maxMass < massM) ) massM=maxMass;
 
  //special case... if the parent mass is _fixed_ we can do a little better
  //and only for a two body decay as that seems to be where we have problems
 
  // Define relativistic propagator amplitude
 
  EvtTwoBodyVertex vd(massD1,massD2,_mass,Lmin/2);
  vd.set_f(_blatt);
  EvtPropBreitWignerRel bw(_mass,_width);
  EvtMassAmp amp(bw,vd);
  if ( _includeDecayFact) {
    amp.addDeathFact();
    amp.addDeathFactFF();
  }
 
  if(fabs(_addFactorPn) >0.00000001){
    //    std::cout<<"EvtRelBreitWignerBarrierFact "<< _addFactorPn<<std::endl;
    amp.addFactorPn( _addFactorPn);
  }
  if ( massParent>-1.) {
    if ( _includeBirthFact ) {
 
      EvtTwoBodyVertex vb(_mass,massOthD,massParent,birthl/2);
      if ( _applyFixForSP8 ) vb.set_f(_blatt);
      amp.setBirthVtx(vb);
      amp.addBirthFact();
      amp.addBirthFactFF();
    }
  }
 
 
 
  EvtAmpPdf<EvtPoint1D> pdf(amp);
 
  // Estimate maximum and create predicate for accept reject
 
 
  double tempMaxLoc=_mass;
  if ( maxMass>-0.5 && maxMass<_mass) tempMaxLoc=maxMass;
  double tempMax=_massMax;
  if ( maxMass>-0.5 && maxMass<_massMax) tempMax=maxMass;
  double tempMinMass=_massMin;
  if ( massD1+massD2 > _massMin) tempMinMass=massD1+massD2;
 
  //redo sanity check - is there a solution to our problem.
  //if not return an error condition that is caught by the
  //mass prob calculation above.
  if ( tempMinMass > tempMax ) {
    _errorCond=true;
    return tempMinMass;
  }
  
  if ( tempMaxLoc < tempMinMass) tempMaxLoc=tempMinMass;
 
  double safetyFactor=1.2;
  if ( _applyFixForSP8 ) safetyFactor=1.4;
 
  EvtPdfMax<EvtPoint1D> max(safetyFactor*pdf.evaluate(EvtPoint1D(tempMinMass,tempMax,tempMaxLoc)));
 
  EvtPdfPred<EvtPoint1D> pred(pdf);
  pred.setMax(max);
 
  EvtIntervalFlatPdf flat(tempMinMass,tempMax);
  EvtPdfGen<EvtPoint1D> gen(flat);
  EvtPredGen<EvtPdfGen<EvtPoint1D>,EvtPdfPred<EvtPoint1D> > predgen(gen,pred);
 
  EvtPoint1D point = predgen();
  return point.value();
 
}

operator=()

EvtRelBreitWignerBarrierFact & EvtRelBreitWignerBarrierFact::operator=

(

const EvtRelBreitWignerBarrierFact &

x

)

Definition at line 84 of file EvtRelBreitWignerBarrierFact.cc.

                                                                                                           {
  _mass=x._mass;
  _width=x._width;
  _spin=x._spin;
  _massMax=x._massMax;
  _massMin=x._massMin;
  _blatt=x._blatt;
  _maxRange=x._maxRange;
  _includeDecayFact=x._includeDecayFact;
  _includeBirthFact=x._includeBirthFact;
  _errorCond=x._errorCond;
 
  return *this;
}

reSetBlatt()

virtual void EvtRelBreitWignerBarrierFact::reSetBlatt ( double  blatt )

inlinevirtual

Reimplemented from EvtAbsLineShape.

Definition at line 46 of file EvtRelBreitWignerBarrierFact.hh.

{ _blatt = blatt; }

Member Data Documentation

_blatt

double EvtRelBreitWignerBarrierFact::_blatt

protected

Definition at line 50 of file EvtRelBreitWignerBarrierFact.hh.

Referenced by EvtRelBreitWignerBarrierFact(), getRandMass(), operator=(), and reSetBlatt().

_errorCond

bool EvtRelBreitWignerBarrierFact::_errorCond

protected

Definition at line 51 of file EvtRelBreitWignerBarrierFact.hh.

Referenced by EvtRelBreitWignerBarrierFact(), getMassProb(), getRandMass(), and operator=().

---

The documentation for this class was generated from the following files:

• EvtRelBreitWignerBarrierFact.hh
• EvtRelBreitWignerBarrierFact.cc