#include <BW_BW.h>

Inheritance diagram for BW_BW:

Public Member Functions
virtual double	prSq () const

virtual double	prSqForGofM () const

virtual double	pABSq ()

virtual double	prSqMax () const

	BW_BW (const AssociatedDecayTree &decay, const std::string &lineshapePrefix="", MINT::MinuitParameterSet *mps=0)

	BW_BW (const BW_BW &other)

virtual	~BW_BW ()

virtual std::complex< double >	getVal (IDalitzEvent &evt)

virtual void	print (IDalitzEvent &evt, std::ostream &out=std::cout)

virtual void	print (std::ostream &out=std::cout) const

virtual DalitzCoordinate	getDalitzCoordinate (double nSigma=3) const

virtual MINT::counted_ptr< IGenFct >	generatingFunction () const

virtual void	setGenerationLimits (double mi, double ma)

virtual std::string	name () const

virtual int	twoLPlusOne () const

virtual std::complex< double >	EvtGenValue (IDalitzEvent &evt)

Public Member Functions inherited from ILineshape
virtual	~ILineshape ()

Public Member Functions inherited from MINT::FitParDependent
virtual unsigned int	size () const

virtual const FitParRef &	operator[] (unsigned int i) const

virtual FitParRef &	operator[] (unsigned int i)

virtual bool	changedSinceLastCall () const

virtual void	rememberFitParValues ()

virtual bool	registerFitParDependence (const IFitParDependent &fpd)

bool	registerFitParDependence (const FitParRef &fpr)

void	removeAllFitParDependencies ()

	FitParDependent (IFitParRegister *daddy=0)

	FitParDependent (const FitParDependent &other, IFitParRegister *newDaddy=0)

void	listFitParDependencies (std::ostream &os=std::cout) const

Public Member Functions inherited from MINT::PolymorphVector< FitParRef >
	PolymorphVector ()

	PolymorphVector (unsigned int N)

	PolymorphVector (unsigned int N, const FitParRef &c)

	PolymorphVector (const PolymorphVector &other)

	PolymorphVector (const typename std::vector< FitParRef > &other)

virtual	~PolymorphVector ()

std::vector< FitParRef > &	theVector ()

const std::vector< FitParRef > &	theVector () const

FitParRef &	operator[] (unsigned int i)

const FitParRef &	operator[] (unsigned int i) const

FitParRef &	at (unsigned int i)

const FitParRef &	at (unsigned int i) const

std::vector< FitParRef >::iterator	begin ()

std::vector< FitParRef >::const_iterator	begin () const

std::vector< FitParRef >::iterator	end ()

std::vector< FitParRef >::const_iterator	end () const

std::vector< FitParRef >::iterator	find (const FitParRef &c)

std::vector< FitParRef >::const_iterator	find (const FitParRef &c) const

FitParRef &	front ()

const FitParRef &	front () const

FitParRef &	back ()

const FitParRef &	back () const

unsigned int	size () const

bool	empty () const

void	push_back (const FitParRef &c)

void	pop_back ()

void	erase (typename std::vector< FitParRef >::iterator pos)

void	erase (typename std::vector< FitParRef >::iterator first, typename std::vector< FitParRef >::iterator last)

PolymorphVector< FitParRef > &	operator= (const PolymorphVector< FitParRef > &other)

void	clear ()

void	resize (unsigned int N)

void	resize (unsigned int N, const FitParRef &c)

	operator const typename std::vector< FitParRef > & () const

	operator typename std::vector< FitParRef > & ()

bool	operator== (const MINT::PolymorphVector< FitParRef > &v2) const

bool	operator!= (const MINT::PolymorphVector< FitParRef > &v2) const

bool	operator< (const MINT::PolymorphVector< FitParRef > &v2) const

bool	operator> (const MINT::PolymorphVector< FitParRef > &v2) const

Protected Member Functions
const MINT::MinuitParameterSet *	getMinuitParameterSet () const

MINT::MinuitParameterSet *	getMinuitParameterSet ()

const std::string &	prefix () const

bool	substitutePDGForReco () const

virtual double	twoBody_dgtPsq_in_MumsFrame (double mumsMass, double mA, double mB) const

virtual double	twoBody_dgtPsq_in_MumsPDGFrame () const

virtual double	twoBody_recodgtPsq_in_MumsPDGFrame () const

virtual double	twoBody_dgtPsq_in_MumsRecoFrame ()

bool	setEventPtr (IDalitzEvent &evt) const

IDalitzEvent *	getEvent () const

bool	setAllFitParameters ()

ResonancePropertiesList *	resonancePropertiesList () const

const ResonanceProperties *	resonanceProperties () const

virtual const ParticleProperties *	mumsProperties () const

virtual ResonancePropertiesFitRef &	mumsFittableProperties () const

virtual int	mumsPID () const

virtual double	mumsMass () const

virtual double	mumsWidth () const

virtual double	mumsRadius () const

virtual double	globalRadius () const

virtual double	Radius () const

double	GetAlpha () const

virtual int	lowestPossibleTwoLPlusOne () const

virtual int	maxDaughterSpinSum () const

virtual int	minDaughterSpinSum () const

virtual int	minDaughterSpinSum2 () const

virtual int	minDaughterSpinSum3 () const

virtual int	minDaughterSpinSum4 () const

virtual int	maxDaughterPairSpinSum (int i, int j) const

virtual int	minDaughterPairSpinSum (int i, int j) const

virtual std::string	mumsSpin () const

virtual int	mumsSpinValue () const

virtual int	mumsParity () const

virtual bool	startOfDecayChain () const

virtual bool	parityConservingL (int L) const

virtual double	mumsRecoMass2 () const

virtual double	mumsRecoMass () const

virtual MultiQuarkContent	mumsQuarkContent () const

bool	nonResonant () const

virtual int	numDaughters () const

virtual TLorentzVector	daughterP4 (int i) const

virtual double	daughterPDGMass (const int &i) const

virtual double	daughterWidth (int i) const

virtual double	daughterRecoMass2 (int i) const

virtual double	daughterRecoMass (int i) const

virtual std::string	daughterSpin (int i) const

virtual int	daughterSpinValue (int i) const

virtual MultiQuarkContent	daughterQuarkContent (int i) const

virtual int	daughterP (int i) const

virtual int	dgtrsInternalParity () const

virtual MultiQuarkContent	dgtrsQuarkContent () const

virtual bool	isWeakDecay () const

virtual double	Fr_PDG_BL ()

virtual double	Fr_BELLE (double prSquared)

virtual double	Fr_BELLE_Max ()

virtual double	Fr ()

virtual double	FrForGofM ()

virtual double	FrMax ()

virtual double	GofM ()

virtual std::complex< double >	BreitWigner ()

virtual void	resetInternals ()

virtual void	resetPDG ()

Protected Attributes
MINT::MinuitParameterSet *	_mps

std::string	_prefix

MINT::NamedParameter< int >	_normBF

bool	_useGlobalRadius

const AssociatedDecayTree &	_theDecay

double	_gen_s_mi

double	_gen_s_ma

ResonancePropertiesList *	_RPL

ResonancePropertiesFitRef *	_fittableResonancePropertiesPtr

MINT::FitParRef *	_fittableGlobalRadiusPtr

Protected Attributes inherited from MINT::PolymorphVector< FitParRef >
std::vector< FitParRef >	_vec

Private Member Functions
double	peakShift () const

double	peakPosition () const

void	makeGeneratingFunction () const

Private Attributes
IDalitzEvent *	_eventPtr

double	_prSq

double	_prSqForGofM

double	_pABSq

double	_mumsRecoMass2

double	_mumsRecoMass

double	_Fr_BELLE

double	_Fr_PDG_BL

double	_GofM

int	_mumsPID

bool	_mumsPID_set

int	_mumsParity

int	_dgtrsInternalParity

int	_twoLPlusOne

std::vector< double >	_daughterRecoMass2

std::vector< double >	_daughterPDGMass

std::vector< double >	_daughterWidth

bool	_substitutePDGForReco

MINT::counted_ptr< IGenFct >	_genFct

Detailed Description

Definition at line 30 of file BW_BW.h.

Constructor & Destructor Documentation

◆ BW_BW() [1/2]

BW_BW::BW_BW	(	const AssociatedDecayTree &	decay,
		const std::string &	lineshapePrefix = `""`,
		MINT::MinuitParameterSet *	mps = `0`
	)

Definition at line 26 of file BW_BW.cpp.

   : _eventPtr(0)
   , _prSq(-9999.0)
   , _prSqForGofM(-9999.0)
   , _pABSq(-9999.0)
     //  , _mumsPDGMass(-9999.0)
     //  , _mumsPDGWidth(-9999.0) 
   , _mumsRecoMass2(-9999.0)
   , _mumsRecoMass(-9999.0)
   , _Fr_BELLE(-9999.0)
   , _Fr_PDG_BL(-9999.0)
   , _GofM(-9999.0)
   , _mumsPID(0)
   , _mumsPID_set(false)
     //  , _mumsPDGRadius(-9999.0)
   , _mumsParity(0)
   , _dgtrsInternalParity(0)
   , _twoLPlusOne(-9999)
   , _daughterRecoMass2(2, -9999.0)
   , _daughterPDGMass(2, -9999.0)
   , _daughterWidth(2, -9999.0)
   , _substitutePDGForReco(false)
   , _genFct(0)
   , _mps(mps)
   , _prefix(lineshapePrefix)
   , _normBF("NormBF", 2)
   , _useGlobalRadius(true)
   , _theDecay(decay)
   , _gen_s_mi(0)
   , _gen_s_ma(0)
   , _RPL(0)
   , _fittableResonancePropertiesPtr(0)
   , _fittableGlobalRadiusPtr(0)
 {
   //  const fitSetup* fs = fitSetup::getMe();
   resetPDG();
   if(0 == _mps) _mps = MinuitParameterSet::getDefaultSet();
 
   if(_theDecay.nDgtr() > 2 && ! (startOfDecayChain())){
     std::cout << "WARNING: BW_BW can only properly handle"
           << "\n    > two body decays."
           << "\n    > This: " << _theDecay.oneLiner()
           << "\n    > is a " << _theDecay.nDgtr() << " body decay."
           << "\n    > Please improve me."
           << std::endl;
     std::cout << " startOfDecayChain()? " << startOfDecayChain()
           << "  PID " << mumsPID()
           << std::endl;
   }
   /*    
     _RPL = ResonancePropertiesList::getMe();
     if(_RPL->get(mumsPID())==0){
         ResonanceProperties* rp= new ResonanceProperties(mumsProperties()->name());
         _RPL->AddToList(rp);
     }  
   */
   setAllFitParameters();
 }

◆ BW_BW() [2/2]

BW_BW::BW_BW ( const BW_BW & other )

Definition at line 84 of file BW_BW.cpp.

   : ILineshape()
   , FitParDependent(other)
   , _prSq(other._prSq)
   , _prSqForGofM(other._prSqForGofM)
   , _pABSq(other._pABSq)
     //  , _mumsPDGMass(other._mumsPDGMass)
     //  , _mumsPDGWidth(other._mumsPDGWidth) 
   , _mumsRecoMass2(other._mumsRecoMass2)
   , _mumsRecoMass(other._mumsRecoMass)
   , _Fr_BELLE(other._Fr_BELLE)
   , _Fr_PDG_BL(other._Fr_PDG_BL)
   , _GofM(other._GofM)
   , _mumsPID(other._mumsPID)
   , _mumsPID_set(other._mumsPID_set)
     //  , _mumsPDGRadius(other._mumsPDGRadius)
   , _mumsParity(other._mumsParity)
   , _dgtrsInternalParity(other._dgtrsInternalParity)
   , _twoLPlusOne(other._twoLPlusOne)
   , _daughterRecoMass2(other._daughterRecoMass2)
   , _daughterPDGMass(other._daughterPDGMass)
   , _daughterWidth(other._daughterWidth)
   , _substitutePDGForReco(other._substitutePDGForReco)
   , _genFct(other._genFct)
   , _mps(other._mps)
   , _prefix(other._prefix)
   , _normBF(other._normBF)
   , _useGlobalRadius(other._useGlobalRadius)
   , _theDecay(other._theDecay)
   , _gen_s_mi(other._gen_s_mi)
   , _gen_s_ma(other._gen_s_ma)
   , _RPL(other._RPL)
   , _fittableResonancePropertiesPtr(0)
   , _fittableGlobalRadiusPtr(0)
 {
   
   /*  if(_RPL->get(mumsPID())==0){
         ResonanceProperties* rp= new ResonanceProperties(mumsProperties()->name());
         _RPL->AddToList(rp);
     }  
   */
   setAllFitParameters();
 }

◆ ~BW_BW()

BW_BW::~BW_BW ( )

virtual

Definition at line 128 of file BW_BW.cpp.

              {
   delete _fittableResonancePropertiesPtr;
   delete  _fittableGlobalRadiusPtr;
 }

Member Function Documentation

◆ BreitWigner()

std::complex< double > BW_BW::BreitWigner ( )

protectedvirtual

Reimplemented in Histo_BW, Bugg_BW, FocusFlatte, Flatte, SBW, and GLass.

Definition at line 1093 of file BW_BW.cpp.

                                       {
   double mass = mumsMass();
   double width = mumsWidth();
 
   double gamma = sqrt(mass*mass*(mass*mass+width*width));
   double k = mass*width*gamma/sqrt(mass*mass+gamma);
 
   std::complex<double> BW(mass*mass - mumsRecoMass2(), mass * GofM());
   double den = (mass*mass - mumsRecoMass2())*(mass*mass - mumsRecoMass2()) + mass * GofM() * mass * GofM();
 
   //std::complex<double> invBW(mass*mass - mumsRecoMass2(), - mass * GofM());
   //return 1.*GeV*GeV/invBW;
   return sqrt(k)*BW/den;
   //  return std::complex<Double_t>(ReBreitWigner(), ImBreitWigner());
 }

◆ daughterP()

int BW_BW::daughterP ( int i ) const

protectedvirtual

Definition at line 255 of file BW_BW.cpp.

                                {
   if( i >= numDaughters() || i < 0){
     cout << " ERROR in BW_BW::daughterSpinValue:"
      << " You requested the spin of dgtr number " << i
      << ". There are " << numDaughters()
      << " daughters." << endl;
     return 0;
   }
   const ParticleProperties* pp = _theDecay.getDgtrVal(i).props();
   if(0==pp){
     cout << "ERROR in BW_BW::daughterP(" << i <<")"
      << " can't find properties for first element"
      << " in this decay tree\n" << _theDecay << endl;
     throw "invalid decay tree in BW_BW::daughterP()";
   }
   return pp->P();
 }

◆ daughterP4()

TLorentzVector BW_BW::daughterP4 ( int i ) const

protectedvirtual

Definition at line 566 of file BW_BW.cpp.

                                            {
   // for diagnostics only.
   if( i >= _theDecay.nDgtr() || i < 0){
     cout << " ERROR in BW_BW::daughterP4:"
      << " You requested the 4-momentum of dgtr number " << i
      << ". There are " << _theDecay.nDgtr() 
      << " daughters." << endl;
     exit(1);
   }
   const_counted_ptr<AssociatedDecayTree> dgtr = _theDecay.getDgtrTreePtr(i);
   std::vector<int> asi = dgtr->getVal().asi();
   if(asi.size() < 2){
     return getEvent()->p(asi[0]);
   }else{
     return TLorentzVector();
   }
 }

◆ daughterPDGMass()

double BW_BW::daughterPDGMass ( const int & i ) const

protectedvirtual

Definition at line 648 of file BW_BW.cpp.

                                                  {
   if( i >= _theDecay.nDgtr() || i < 0 ){
     cout << " ERROR in BW_BW::daughterPDGMass:"
      << " You requested the mass of dgtr number " << i
      << ". There are " << _theDecay.nDgtr()
      << " daughters." << endl;
     return -9999;
   }
 
   if( _daughterPDGMass[i] < 0 ){
     _daughterPDGMass[i] = _theDecay.getDgtrVal(i).mass();
   }
 
   if( _theDecay.getDgtrVal(i).isNonResonant() ){
     _daughterPDGMass[i] = daughterRecoMass(i);
   }
 
   return _daughterPDGMass[i];
 }

◆ daughterQuarkContent()

MultiQuarkContent BW_BW::daughterQuarkContent ( int i ) const

protectedvirtual

Definition at line 692 of file BW_BW.cpp.

                                                         {
   if( i >= _theDecay.nDgtr() || i < 0){
     cout << " ERROR in BW_BW::daughterQuarkContent"
      << " You requested the spin of dgtr number " << i
      << ". There are " << _theDecay.nDgtr() 
      << " daughters." 
      << "\n\t I'll crash now. Sorry."
      << endl;
     throw "can't handle this";
   }
 
   return netQuarkContent(*(_theDecay.getDgtrTreePtr(i)));
   //return _theDecay.getDgtrVal(i).props()->netQuarkContent();
 }

◆ daughterRecoMass()

double BW_BW::daughterRecoMass ( int i ) const

protectedvirtual

Definition at line 617 of file BW_BW.cpp.

                                          {
   if(substitutePDGForReco()) return daughterPDGMass(i);
 
   double m2 = daughterRecoMass2(i);
   if(m2 < 0){
     cout << " ERROR negative reconstructed mass2"
      << " in decay:\n " << _theDecay << endl;
     return - sqrt(-m2);
   }
   //  cout << " for index " << i << endl;
   //  cout << " m2 = " << m2 << " mass = " << sqrt(m2) << endl;
   return sqrt(m2);
 }

◆ daughterRecoMass2()

double BW_BW::daughterRecoMass2 ( int i ) const

protectedvirtual

Definition at line 583 of file BW_BW.cpp.

                                           {
   if(substitutePDGForReco()) return daughterPDGMass(i)*daughterPDGMass(i);
 
   if( i >= _theDecay.nDgtr() || i < 0){
     cout << " ERROR in BW_BW::daughterRecoMass2:"
      << " You requested the mass of dgtr number " << i
      << ". There are " << _theDecay.nDgtr() 
      << " daughters." << endl;
     return -9999;
   }
   if(_daughterRecoMass2[i] < 0){
     
     const_counted_ptr<AssociatedDecayTree> dgtr = _theDecay.getDgtrTreePtr(i);
     if(dgtr->isFinalState()){
       double m = daughterPDGMass(i);
       _daughterRecoMass2[i] = m*m;
       //_daughterRecoMass2[i] = daughterP4(i).M2();
     }else{
       std::vector<int> asi = dgtr->getVal().asi();
       if(asi.size() < 2){
     cout << "ERROR in BW_BW::daughterRecoMass2() "
          << dgtr->getVal()._pdg_id
          << ", " << dgtr->getVal().name()
          << " is associated to " << asi.size() << " particles"
          << " and is not stable? "
          << _theDecay
          << endl;
     
       } 
       _daughterRecoMass2[i] = getEvent()->sij(asi);
     }
   }
   return _daughterRecoMass2[i];
 }

◆ daughterSpin()

std::string BW_BW::daughterSpin ( int i ) const

protectedvirtual

Definition at line 682 of file BW_BW.cpp.

                                         {
   if( i >= _theDecay.nDgtr() || i < 0){
     cout << " ERROR in BW_BW::daughterSpinValue:"
      << " You requested the spin of dgtr number " << i
      << ". There are " << _theDecay.nDgtr() 
      << " daughters." << endl;
     return "-9999";
   }
   return _theDecay.getDgtrVal(i).J();
 }

◆ daughterSpinValue()

int BW_BW::daughterSpinValue ( int i ) const

protectedvirtual

Definition at line 707 of file BW_BW.cpp.

                                        {
   return atoi(daughterSpin(i).c_str());
 }

◆ daughterWidth()

double BW_BW::daughterWidth ( int i ) const

protectedvirtual

Definition at line 668 of file BW_BW.cpp.

                                       {
   if( i >= _theDecay.nDgtr() || i < 0){
     cout << " ERROR in BW_BW::daughterWidth:"
      << " You requested the mass of dgtr number " << i
      << ". There are " << _theDecay.nDgtr() 
      << " daughters." << endl;
     return -9999;
   }
   if(_daughterWidth[i] < 0){
     _daughterWidth[i] =  _theDecay.getDgtrVal(i).width();
   }
   return _daughterWidth[i];
 }

◆ dgtrsInternalParity()

int BW_BW::dgtrsInternalParity ( ) const

protectedvirtual

Definition at line 273 of file BW_BW.cpp.

                                     {
   if( 0 == _dgtrsInternalParity ){
     _dgtrsInternalParity = 1;
     for(int i=0; i < numDaughters(); i++){
       _dgtrsInternalParity *= daughterP(i);
     }
   }
   return _dgtrsInternalParity;
 
 }

◆ dgtrsQuarkContent()

MultiQuarkContent BW_BW::dgtrsQuarkContent ( ) const

protectedvirtual

Definition at line 225 of file BW_BW.cpp.

                                                 {
   MultiQuarkContent sum;
   for(int i=0; i < numDaughters(); i++){
     sum += daughterQuarkContent(i);
   }
   return sum;
 }

◆ EvtGenValue()

std::complex< double > BW_BW::EvtGenValue ( IDalitzEvent & evt )

virtual

Definition at line 1304 of file BW_BW.cpp.

                                                       {
   // For debugging only.  
   // This is the EvtGen version of the full 3-body amplitude
   // (cut'n'paste with minimal changes).  
   // BW_BW only returns only the Breit Wigner plus form factor for one
   // resonance. So don't compare with BW_BW, but with Amplitude.  This
   // function is part of BW_BW for convencience, BW_BW provides all
   // the "services" needed to compute this.
 
   setEventPtr(evt);
   double pi180inv = pi/108.0;//1.0/EvtConst::radToDegrees;
 
   complex<double> ampl(0,0);
   double magn =1;
   double theta = 0;
   
   int spin = (twoLPlusOne() -1)/2;
 
   // get parent
   if(0 == _theDecay.hasParent()) {
     return -9999;
   }  
   TLorentzVector p4_p = getEvent()->p(0)  * (1./GeV);
   TLorentzVector p4_d1 = getEvent()->p(_theDecay.getDgtrTreePtr(0)->getVal().asi()[0])   * (1./GeV);
   TLorentzVector p4_d2 = getEvent()->p(_theDecay.getDgtrTreePtr(1)->getVal().asi()[0])   * (1./GeV);
   TLorentzVector p4_d3 = p4_p-p4_d1-p4_d2;
 
   //get cos of the angle between the daughters from their 4-momenta
   //and the 4-momentum of the parent
 
   //in general, EvtDecayAngle(parent, part1+part2, part1) gives the angle
   //the missing particle (not listed in the arguments) makes
   //with part2 in the rest frame of both
   //listed particles (12)
  
   //angle 3 makes with 2 in rest frame of 12 (CS3)  
   //double cos_phi_0 = EvtDecayAngle(p4_p, p4_d1+p4_d2, p4_d1);
   //angle 3 makes with 1 in 12 is, of course, -cos_phi_0
 
   //first compute several quantities...follow CLEO preprint 00-23
 
   bool invmass_angdenom = false;
   double bwm = mumsMass()  / GeV;
   double gamma = mumsWidth()  / GeV;
 
   double mAB=(p4_d1+p4_d2).M();
   double mBC=(p4_d2+p4_d3).M();
   double mAC=(p4_d1+p4_d3).M();
   double mA=p4_d1.M(); 
   double mB=p4_d2.M(); 
   double mD=p4_p.M();
   double mC=p4_d3.M();
   
   double mR=bwm;
   double mdenom = invmass_angdenom ? mAB : mR;
   double gammaR = gamma;
   double pAB=sqrt( (((mAB*mAB-mA*mA-mB*mB)*(mAB*mAB-mA*mA-mB*mB)/4.0) -
             mA*mA*mB*mB)/(mAB*mAB));
   double pR=sqrt( (((mR*mR-mA*mA-mB*mB)*(mR*mR-mA*mA-mB*mB)/4.0) -
            mA*mA*mB*mB)/(mR*mR));
 
   double pD= (((mD*mD-mR*mR-mC*mC)*(mD*mD-mR*mR-mC*mC)/4.0) -
            mR*mR*mC*mC)/(mD*mD);
   if ( pD>0 ) { pD=sqrt(pD); } else {pD=0;}
   double pDAB=sqrt( (((mD*mD-mAB*mAB-mC*mC)*(mD*mD-mAB*mAB-mC*mC)/4.0) -
            mAB*mAB*mC*mC)/(mD*mD));
 
 
   //  report(INFO,"EvtGen")  << mAB<<" "<< mBC<<" "<< mAC<<" "<< mA<<" "<< mB<<" "<< mC<<" "
   //        << mD<<" "<< mR<<" "<< gammaR<<" "<< pAB<<" "<< pR<<" "<< pD<<" "<<pDAB<<endl;
 
   double fR=1;
   double fD=1;
   int power=0;
   switch (spin) {
   case 0:
     fR=1.0;
     fD=1.0;
     power=1;
     //report(INFO,"EvtGen") << "fR="<<fR<<" fD="<<fD<<endl;
     break;
   case 1:
     fR=sqrt(1.0+1.5*1.5*pR*pR)/sqrt(1.0+1.5*1.5*pAB*pAB);
     fD=sqrt(1.0+5.0*5.0*pD*pD)/sqrt(1.0+5.0*5.0*pDAB*pDAB);
     //report(INFO,"EvtGen") << "fR="<<fR<<" fD="<<fD<<endl;
     power=3;
     break;
   case 2:
     fR = sqrt( (9+3*pow((1.5*pR),2)+pow((1.5*pR),4))/(9+3*pow((1.5*pAB),2)+pow((1.5*pAB),4)) );
     fD = sqrt( (9+3*pow((5.0*pD),2)+pow((5.0*pD),4))/(9+3*pow((5.0*pDAB),2)+pow((5.0*pDAB),4)) );
     power=5;
     //report(INFO,"EvtGen") << "fR="<<fR<<" fD="<<fD<<std::endl;
     break;
 
     
 
   default:
     ;//report(INFO,"EvtGen") << "Incorrect spin in EvtResonance22.cc\n";
   }
 
   cout << "EvtGen: fR = " << fR << ", fD = " << fD << endl;
 
   //  report(INFO,"EvtGen") << "fR = " << fR << " fD = " << fD << endl;
   
   double gammaAB= gammaR*pow(pAB/pR,power)*(mR/mAB)*fR*fR;
   cout << "EvtGen: gamma AB = " << gammaAB << endl;
   //report(INFO,"EvtGen") << gammaAB<<endl;
   double sf;
   switch (spin) {
   case 0:
     ampl=magn*complex<double>(cos(theta*pi180inv),sin(theta*pi180inv))*
           fR*fD/(mR*mR-mAB*mAB-complex<double>(0.0,mR*gammaAB));
     sf=1;
     break;
   case 1:
     ampl=magn*complex<double>(cos(theta*pi180inv),sin(theta*pi180inv))*
       (fR*fD*(mAC*mAC-mBC*mBC+((mD*mD-mC*mC)*(mB*mB-mA*mA)/(mdenom*mdenom)))/
        (mR*mR-mAB*mAB-complex<double>(0.0,mR*gammaAB)));
     
     sf = mAC*mAC-mBC*mBC+((mD*mD-mC*mC)*(mB*mB-mA*mA)/(mdenom*mdenom));
     cout << "EvtGen: sf spin 1 : " <<  sf << endl;
     break;
   case 2:
 //     ampl=magn*complex<double>(cos(theta*pi180inv),sin(theta*pi180inv))*
 //  fR*fD/(mR*mR-mAB*mAB-complex<double>(0.0,mR*gammaAB))*
 //        (pow((mBC*mBC-mAC*mAC+(mD*mD-mC*mC)*(mA*mA-mB*mB)/(mAB*mAB)),2)-
 //  (1.0/3.0)*(mAB*mAB-2*mD*mD-2*mC*mC+pow((mD*mD- mC*mC)/mAB, 2))*
 //  (mAB*mAB-2*mA*mA-2*mB*mB+pow((mA*mA-mB*mB)/mAB,2)));
      ampl=magn*complex<double>(cos(theta*pi180inv),sin(theta*pi180inv))*
     fR*fD/(mR*mR-mAB*mAB-complex<double>(0.0,mR*gammaAB))*
         (pow((mBC*mBC-mAC*mAC+(mD*mD-mC*mC)*(mA*mA-mB*mB)/(mdenom*mdenom)),2)-
     (1.0/3.0)*(mAB*mAB-2*mD*mD-2*mC*mC+pow((mD*mD- mC*mC)/mdenom, 2))*
     (mAB*mAB-2*mA*mA-2*mB*mB+pow((mA*mA-mB*mB)/mdenom,2)));
      
      sf =  pow((mBC*mBC-mAC*mAC+(mD*mD-mC*mC)*(mA*mA-mB*mB)/(mdenom*mdenom)),2)-
        (1.0/3.0)*(mAB*mAB-2*mD*mD-2*mC*mC+pow((mD*mD- mC*mC)/mdenom, 2))*
        (mAB*mAB-2*mA*mA-2*mB*mB+pow((mA*mA-mB*mB)/mdenom,2));
      cout << "EvtGen: sf spin 2 " << sf << endl;
    break;
 
   default:
     ;//report(INFO,"EvtGen") << "Incorrect spin in EvtResonance22.cc\n";
   }
 
   //report(INFO,"EvtGen") <<"The amplitude is "<<ampl<<endl;
   return ampl;
 
 }

◆ Fr()

double BW_BW::Fr ( )

protectedvirtual

Definition at line 710 of file BW_BW.cpp.

                 {
   return Fr_BELLE(prSq());
 }

◆ Fr_BELLE()

double BW_BW::Fr_BELLE ( double prSquared )

protectedvirtual

Definition at line 720 of file BW_BW.cpp.

                                       {
   // (ratio of) Blatt-Weisskopf penetration factors
   // as in BELLE note hep-ex/0406067, pg 5.
   bool dbThis=false;
   //if(startOfDecayChain()) dbThis=true;
 
   if(dbThis){
     cout << " mums name " << mumsProperties()->name() << endl;
     cout << " decay chain " << _theDecay.oneLiner() << endl;
     cout << " twoLPlusOne: " << twoLPlusOne() << endl;
   }
 
   double Fr_BELLE;  
     
   //if(_Fr_BELLE < 0){
     if(twoLPlusOne() == 1){                        // spin = 0
       Fr_BELLE = 1;
       if(dbThis){
        cout << " _Fr_BELLE = " << Fr_BELLE << endl;
     }
       return Fr_BELLE;
     }
     double R = Radius(); // or mumsRadius()  ???
     double R_pr_sq  = R*R * prSquared;
     double R_pAB_sq = R*R * pABSq();
     if(twoLPlusOne() == 3){                        // spin = 1
       if(dbThis) cout << "in Fr() " << sqrt(prSqForGofM())
               << ", " << sqrt(prSquared) << endl;
       
       if(dbThis)cout << "BW Fr R = " << R << "prSquared = " << prSqForGofM()
              << " pABSq() " << prSquared
              << endl;
       Fr_BELLE = sqrt( (1. + R_pr_sq)/(1. + R_pAB_sq) );
     }else if(twoLPlusOne() == 5){                   // spin == 2
       double R_pr_4  = R_pr_sq*R_pr_sq;
       double R_pAB_4 = R_pAB_sq*R_pAB_sq;
       
       Fr_BELLE = sqrt( (  9. + 3.*R_pr_sq  + R_pr_4)
           /( 9. + 3.*R_pAB_sq + R_pAB_4) 
           );
     } else{
       std::cout << "BW_BW for decay \n" << _theDecay 
         << "\n ERROR! Can't deal with L == " 
         << twoLPlusOne()
         << ". Please improve me!"
         << std::endl;
       Fr_BELLE = -9999;
     }
   //}
     
     
   if(dbThis){
     cout << "returning _Fr_BELLE = " << Fr_BELLE << endl;
   }
   return Fr_BELLE;
 }

◆ Fr_BELLE_Max()

double BW_BW::Fr_BELLE_Max ( )

protectedvirtual

Definition at line 776 of file BW_BW.cpp.

                           {
   // Max value of Blatt-Weisskopf penetration factors
   // as in BELLE note hep-ex/0406067, pg 5.
   bool dbThis=false;
   //if(startOfDecayChain()) dbThis=true;
 
   if(dbThis){
     cout << " mums name " << mumsProperties()->name() << endl;
     cout << " decay chain " << _theDecay.oneLiner() << endl;
     cout << " twoLPlusOne: " << twoLPlusOne() << endl;
   }
 
   double FrBelleMax=-9999;
 
   if(twoLPlusOne() == 1){                        // spin = 0
     FrBelleMax = 1;
     if(dbThis){
       cout << " FrBelleMax = " << FrBelleMax << endl;
     }
     return FrBelleMax;
   }
   double R = Radius(); // or mumsRadius()  ???
   //  double R_pr_sq  = R*R * prSq();
   double R_pr_sq  = R*R * prSq();
   double R_pAB_sq = 0;//R*R * pABSq();
   if(twoLPlusOne() == 3){                        // spin = 1
     if(dbThis) cout << "in FrMax() " << sqrt(prSq())
             << ", " << sqrt(pABSq()) << endl;
     
     if(dbThis)cout << "BW FrMax R = " << R << "prSq() = " << prSq() 
            << " pABSq() " << pABSq()
            << endl;
     FrBelleMax = sqrt( (1. + R_pr_sq)/(1. + R_pAB_sq) );
   }else if(twoLPlusOne() == 5){                   // spin == 2
     double R_pr_4  = R_pr_sq*R_pr_sq;
     double R_pAB_4 = 0;//R_pAB_sq*R_pAB_sq;
     
     FrBelleMax = sqrt( (  9. + 3.*R_pr_sq  + R_pr_4)
                /( 9. + 3.*R_pAB_sq + R_pAB_4) 
                );
   } else{
     std::cout << "BW_BW for decay \n" << _theDecay 
           << "\n ERROR! Can't deal with L == " 
           << twoLPlusOne()
           << ". Please improve me!"
           << std::endl;
     FrBelleMax = -9999;
   }
   if(dbThis){
     cout << "returning FrBelleMax = " << FrBelleMax << endl;
     cout << " compare: Fr: " << Fr() << endl;
   }
   return FrBelleMax;
 }

◆ Fr_PDG_BL()

double BW_BW::Fr_PDG_BL ( )

protectedvirtual

Definition at line 830 of file BW_BW.cpp.

                        {
   // (ratio of) Blatt-Weisskopf penetration factors
   // as David Asner's PDG review, Table 1 (page 3),
   // column 2 (BELLE, above is col 3 in the same table)
 
   bool dbThis = false;
   //if(startOfDecayChain()) dbThis=true;
 
   if(dbThis){
     cout << " BW_BW::Fr mums name " << mumsProperties()->name() << endl;
     cout << " decay chain " << _theDecay.oneLiner() << endl;
     cout << " twoLPlusOne: " << twoLPlusOne() << endl;
   }
 
   if(_Fr_PDG_BL < 0){
     if(twoLPlusOne() == 1){                        // spin = 0
       _Fr_PDG_BL = 1;
       if(dbThis){
     cout << " _Fr_PDG_BL = " << _Fr_PDG_BL << endl;
       }
       return _Fr_PDG_BL;
     }
     double R = Radius(); // or mumsRadius()  ???
     //    double R_pr_sq  = R*R * prSq();
     double R_pAB_sq = R*R * pABSq();
 
     double z = R_pAB_sq;
 
     double FrSq = -9999;
 
     if(twoLPlusOne() == 3){                        // spin = 1
       if(dbThis) cout << "in Fr() " << sqrt(prSq())
               << ", " << sqrt(pABSq()) << endl;
       
       if(dbThis)cout << "BW Fr R = " << R << "prSq() = " << prSq() 
              << " pABSq() " << pABSq()
              << endl;
       FrSq = 2.0*z/(1.0 + z);
       if(FrSq >=0) _Fr_PDG_BL = sqrt(FrSq);
       else _Fr_PDG_BL = 0;
     }else if(twoLPlusOne() == 5){                   // spin == 2
       FrSq = 13.0*z*z / ( (z-3.0)*(z-3.0) + 9*z );
       if(FrSq >=0) _Fr_PDG_BL = sqrt(FrSq);
       else _Fr_PDG_BL = 0;     
     }else{
       std::cout << "BW_BW for decay \n" << _theDecay 
         << "\n ERROR! Can't deal with L == " 
         << twoLPlusOne()
         << ". Please improve me!"
         << std::endl;
       _Fr_PDG_BL = -9999;
     }
   }
   if(dbThis){
     cout << "returning _Fr_PDG_BL = " << _Fr_PDG_BL << endl;
   }
   return _Fr_PDG_BL;
 }

◆ FrForGofM()

double BW_BW::FrForGofM ( )

protectedvirtual

Definition at line 713 of file BW_BW.cpp.

                        { 
   return Fr_BELLE(prSqForGofM());
 }

◆ FrMax()

double BW_BW::FrMax ( )

protectedvirtual

Definition at line 716 of file BW_BW.cpp.

                    {
   return Fr_BELLE_Max();
 }

◆ generatingFunction()

counted_ptr< IGenFct > BW_BW::generatingFunction ( ) const

virtual

Implements ILineshape.

Reimplemented in CrystalBallRhoOmega.

Definition at line 1277 of file BW_BW.cpp.

                                                     {
   if(! _genFct){
     makeGeneratingFunction();
   }
   return _genFct;
 }

◆ GetAlpha()

double BW_BW::GetAlpha ( ) const

inlineprotected

Definition at line 116 of file BW_BW.h.

                          {
       return _RPL->get(mumsPID())->alpha();
   }

◆ getDalitzCoordinate()

DalitzCoordinate BW_BW::getDalitzCoordinate ( double nSigma = 3 ) const

virtual

Implements ILineshape.

Reimplemented in CrystalBallRhoOmega.

Definition at line 537 of file BW_BW.cpp.

                                                               {
   DalitzCoordinate coord(_theDecay.getVal().asi());
   double meanM = mumsMass();
   double width = mumsWidth();
   
   double minMass = meanM - nSigma*width/2.0;
   double maxMass = meanM + nSigma*width/2.0;
 
   coord.setMinMax(minMass*minMass, maxMass*maxMass);
   return coord;
 }

◆ getEvent()

IDalitzEvent * BW_BW::getEvent ( ) const

protected

Definition at line 147 of file BW_BW.cpp.

                                    {
   return _eventPtr;
 }

◆ getMinuitParameterSet() [1/2]

const MINT::MinuitParameterSet * BW_BW::getMinuitParameterSet ( ) const

protected

Definition at line 133 of file BW_BW.cpp.

                                                                {
   return _mps;
 }

◆ getMinuitParameterSet() [2/2]

MINT::MinuitParameterSet * BW_BW::getMinuitParameterSet ( )

protected

Definition at line 136 of file BW_BW.cpp.

                                                     {
   if(0 == _mps) _mps = MinuitParameterSet::getDefaultSet();
   return _mps;
 }

◆ getVal()

std::complex< double > BW_BW::getVal ( IDalitzEvent & evt )

virtual

Implements ILineshape.

Reimplemented in RhoOmegaGS, CrystalBallRhoOmega, Lass, GounarisSakurai, CrystalBarrelFOCUS, Model_independent, NonRes, and Rho0Omega.

Definition at line 1138 of file BW_BW.cpp.

                                                  {
   bool dbThis= false; // bigDbg;
 
   setEventPtr(evt);
   resetInternals();
 
   double formFactor= 1.;
   if( _normBF == 1 ) formFactor = Fr();
   else if( _normBF == 0 ) formFactor = Fr_PDG_BL();  
   else if(_normBF == 2 ) formFactor = Fr_BELLE(0.);  
   else if(_normBF == 3 ) formFactor =  pow(pABSq(),GetAlpha());
  
   if( nonResonant() ){
         return formFactor;
   }
     
   if(startOfDecayChain()){
     // in principle there is no need to distinguish the start
     // of the decay chain from the rest - it could just get
     // a Breit Wigner (with a constant value of the width of 
     // the D is zero, as usual). However, 
     // this is to comply with the usual convention: Only the
     // form factor, not the BW-propagator.
     if (compareToOldRooFit){
       return 1;
     }
     if(_theDecay.nDgtr() > 2){
       // all calculations of Fr etc are meaningless in this case
       // assume Fr=1;
       if(dbThis){
           cout << "BW_BW::getVal() for " << name() << endl;
           cout << "instead of Fr() for 2l+1 = " << twoLPlusOne() << endl;
           cout << " I'll return 1 " << endl;
       }
       return 1;
     }
     return formFactor;
   }
   
   if(dbThis){
     if(0 != getEvent()){
       cout << "BW pz " << getEvent()->p(1).Z() << endl;
     }
     cout << " BW_BW for " 
      << _theDecay.oneLiner() << endl; // dbg
     cout << "\n    >  nominalMass " << mumsMass()
      << "\n    > nominalWidth " << mumsWidth()
      << "\n    > GoM() " << GofM()
      << "\n    > Fr() " << Fr()
       << "\n    > FrForGofM " <<FrForGofM()
       << "\n    > Fr_BELLE() " <<Fr_BELLE(0.)
      << "\n    > BW   " << BreitWigner()
      << "\n    > FR*BW = " << Fr() * BreitWigner()
      << "\n    > recoMass " << mumsRecoMass()
      << "\n    > EvtGenValue " << EvtGenValue(*(getEvent()))
      << endl;
   }
     
   const std::complex<double> returnVal = formFactor*BreitWigner();
 
   if(dbThis) cout << " value = " << returnVal 
        << "|A|^2 | " << returnVal.real()*returnVal.real() 
            + returnVal.imag()*returnVal.imag()
        << endl; //dbg
 
   return returnVal;
 }

◆ globalRadius()

double BW_BW::globalRadius ( ) const

protectedvirtual

Definition at line 470 of file BW_BW.cpp.

                                 {
   if(0 == _fittableGlobalRadiusPtr){
     cout << "ERROR in BW_BW::Radius()"
      << " _fittableGlobalRadiusPtr is 0"
      << ", although it should have been set during construction."
      << " Looking at  this decay tree\n" << _theDecay << endl;
     throw "invalid decay tree in BW_BW::Radius()";
   }
   return *_fittableGlobalRadiusPtr;
 }

◆ GofM()

double BW_BW::GofM ( )

protectedvirtual

Reimplemented in FermiPS_BW, Histo_BW, FocusFlatte, Flatte, and SBW.

Definition at line 1035 of file BW_BW.cpp.

                   {
   //  return mumsWidth(); // debug only!
   bool dbThis=false;
   if(_GofM < 0){
     double densq = prSqForGofM();  
     double numsq = pABSq();
     if(numsq < 0){
       cout << "ERROR in BW_BW::GofM() for decay\n" 
        << _theDecay
        << "\n ERROR: pABSq() = " << pABSq() << endl;
       _GofM = -9999;
       return _GofM;
     }
     double epsilon = 1 * eV*eV;
     if(0.0 == densq ) densq = epsilon;
     
     double pratio = 1;
     if(densq > 0 && numsq >=0){
       pratio = sqrt(numsq/densq);
     }
     if(dbThis){
       cout << "pratio = sqrt( " << numsq << " / " << densq << " )" << endl;
     }
     double pratio_to_2Jplus1 = 1;
     if(dbThis) cout  << " pratio " 
                    << pratio 
                //  << "^  2*mumsSpinValue() + 1 = "; //dbg
                    << "^ " << twoLPlusOne() << " = "; //dbg
     //    for(int i=0; i < 2*mumsSpinValue() + 1; i++){
     for(int i=0; i < twoLPlusOne(); i++){
       pratio_to_2Jplus1 *= pratio;
     }
     if(dbThis) cout << pratio_to_2Jplus1 << endl;//dbg
     
     densq = mumsRecoMass2();
     if(densq <= 0){
       cout << "ERROR in BW_BW::GofM() for decay\n"
        << _theDecay
        << ": ERROR: mumsRecoMass() = " 
        << mumsRecoMass() << endl;
       _GofM = -9999;
       return _GofM;
     }
     double mRatio = mumsMass()/sqrt(densq);
     // this factor varies for different final states - check!!
     
     double thisFr = FrForGofM();
     if(dbThis) cout << "GofM() = " << mumsWidth()  << " * "  // dbg
                   << pratio_to_2Jplus1 << " * "
                   << mRatio << " * " 
                   << thisFr << " * " << thisFr
                   << endl;
     _GofM = mumsWidth() * pratio_to_2Jplus1 * mRatio * thisFr*thisFr;
     if(dbThis) cout << " so, GofM/G0 = " << _GofM/mumsWidth() << endl;
   }
   return _GofM;
 }

◆ isWeakDecay()

bool BW_BW::isWeakDecay ( ) const

protectedvirtual

Definition at line 211 of file BW_BW.cpp.

                              {
   return ! mumsQuarkContent().compatible(dgtrsQuarkContent());
 }

◆ lowestPossibleTwoLPlusOne()

int BW_BW::lowestPossibleTwoLPlusOne ( ) const

protectedvirtual

Definition at line 163 of file BW_BW.cpp.

                                          {
   bool dbThis=false;
   int l_1 = abs(mumsSpinValue() - maxDaughterSpinSum());
   int l_2 = abs(mumsSpinValue() - minDaughterSpinSum());
 
   int l = (l_1 < l_2 ? l_1 : l_2);
   int maxL = mumsSpinValue() + maxDaughterSpinSum();
 
   if(dbThis){
     cout << "BW_BW::lowestPossibleTwoLPlusOne() "
      << _theDecay.oneLiner()
      << " before considering parity: " << 2*l+1
      << " for " << mumsSpinValue() 
      << (mumsParity() > 0 ? "+" : "-")
      << " -> " 
      << daughterSpinValue(0)
      << (daughterP(0) > 0 ? "+" : "-")
      << ", "
      << daughterSpinValue(1)
      << (daughterP(1) > 0 ? "+" : "-")
      << " l = " << l
      << " and maxL " << maxL
      << " is weak decay? " << isWeakDecay()
      << endl;
   }
 
   if(! isWeakDecay()){
     while((! parityConservingL(l)) && l < maxL) l++;
   }
   if(dbThis){
     if(isWeakDecay()){
       cout << " is weak decay: mum " << mumsQuarkContent()
        << " daughters: " << dgtrsQuarkContent() << endl;
       cout << "is weak decay - would it have made a difference?"
        << " now l = " << l;
       int tmpL=l;
       while((! parityConservingL(tmpL)) && tmpL < maxL) tmpL++;
       cout << ", would have been: " << tmpL << endl;
     }
     cout << "\tBW_BW::lowestPossibleTwoLPlusOne() "
      << " returning " << 2*l+1
      << " (l=" << l <<")"
      << endl;
   }
 
   return 2*l + 1;
 }

◆ makeGeneratingFunction()

void BW_BW::makeGeneratingFunction ( ) const

private

Definition at line 1245 of file BW_BW.cpp.

                                         {
     if(! _genFct){
         
         std::vector< const_counted_ptr<DDTree<AssociatedDecayTreeItem> > > sister = _theDecay.getSistersTrees() ; 
         
         if(nonResonant()){
             counted_ptr<IGenFct> gptr(new FlatFct(getDalitzCoordinate()));
             _genFct = gptr;
         }
         else if(sister.size()>0){
             for (unsigned int i = 0; i< sister.size(); i++) {
                 if(sister[i]->getVal().props()->isNonResonant()){
                     counted_ptr<IGenFct> gptr(new FlatFct(getDalitzCoordinate()));
                     _genFct = gptr;
                     return;
                 }
             }
         }
         if (_theDecay.nDgtr()==2 && (_theDecay.getDgtrVal(0).isNonResonant() || _theDecay.getDgtrVal(1).isNonResonant()) ){
             counted_ptr<IGenFct> gptr(new BWFct(getDalitzCoordinate()
                                                 , mumsMass()
                                                 , mumsWidth() * 1.02 ));
             _genFct = gptr;
         }else{
             counted_ptr<IGenFct> gptr(new BWFct(getDalitzCoordinate()
                                                 , peakPosition()
                                                 , mumsWidth() * 2. ));
             _genFct = gptr;
         }
     }
 }

◆ maxDaughterPairSpinSum()

int BW_BW::maxDaughterPairSpinSum	(	int	i,
		int	j
	)		const

protectedvirtual

Definition at line 371 of file BW_BW.cpp.

                                                   {
   return daughterSpinValue(i) + daughterSpinValue(j);
 }

◆ maxDaughterSpinSum()

int BW_BW::maxDaughterSpinSum ( ) const

protectedvirtual

Definition at line 284 of file BW_BW.cpp.

                                   {
   int sum=0;
   for(int i=0; i < numDaughters(); i++){
     sum += daughterSpinValue(i);
   }
   return sum;
 }

◆ minDaughterPairSpinSum()

int BW_BW::minDaughterPairSpinSum	(	int	i,
		int	j
	)		const

protectedvirtual

Definition at line 375 of file BW_BW.cpp.

                                                    {
   return abs(daughterSpinValue(i) - daughterSpinValue(j));
 }

◆ minDaughterSpinSum()

int BW_BW::minDaughterSpinSum ( ) const

protectedvirtual

Definition at line 292 of file BW_BW.cpp.

                                    {
   if(2 == numDaughters()){
     return minDaughterSpinSum2();
   }else if(3 == numDaughters()){
     return minDaughterSpinSum3();
   }if(4 == numDaughters()){
     return minDaughterSpinSum4();
   }else{
     cout << "BW_BW::minDaughterSpinSum() WARNING! "
      << " can't handle more than 4 daughters - have "
      << numDaughters() << ". Please improve me." << endl;
     return 0;
   }
 }

◆ minDaughterSpinSum2()

int BW_BW::minDaughterSpinSum2 ( ) const

protectedvirtual

Definition at line 340 of file BW_BW.cpp.

                                     {
   return abs(daughterSpinValue(1) - daughterSpinValue(0));
 }

◆ minDaughterSpinSum3()

int BW_BW::minDaughterSpinSum3 ( ) const

protectedvirtual

Definition at line 343 of file BW_BW.cpp.

                                     {
   int minSpinPair = minDaughterPairSpinSum(0,1);
   int maxSpinPair = maxDaughterPairSpinSum(0,1);
   int otherSpin = daughterSpinValue(2);
   int min = maxSpinPair + otherSpin;
   for(int s = minSpinPair; s <= maxSpinPair; s++){
     int thisMinSpin = abs(s - otherSpin);
     if(thisMinSpin < min) min=thisMinSpin;
   }
   return min;
 }

◆ minDaughterSpinSum4()

int BW_BW::minDaughterSpinSum4 ( ) const

protectedvirtual

Definition at line 354 of file BW_BW.cpp.

                                     {
   int minSpinPair1 = minDaughterPairSpinSum(0,1);
   int maxSpinPair1 = maxDaughterPairSpinSum(0,1);
   
   int minSpinPair2 = minDaughterPairSpinSum(2,3);
   int maxSpinPair2 = maxDaughterPairSpinSum(2,3);
   
   int min = maxSpinPair1 + maxSpinPair2;
   for(int s1 = minSpinPair1; s1 <= maxSpinPair1; s1++){
     for(int s2 = minSpinPair2; s2 <= maxSpinPair2; s2++){
       int thisMinSpin = abs(s1 - s2);
       if(thisMinSpin < min) min=thisMinSpin;
     }
   }
   return min;
 }

◆ mumsFittableProperties()

ResonancePropertiesFitRef & BW_BW::mumsFittableProperties ( ) const

protectedvirtual

Definition at line 414 of file BW_BW.cpp.

                                                               {
   if(0 == _fittableResonancePropertiesPtr){
     cout << "something went wrong in BW_BW::mumsFittableProperties() "
      << " _fittableResonancePropertiesPtr is 0 although it should be"
      << " set at construction."
      << " Looking at  this decay tree\n" << _theDecay << endl;
 
     throw "error in BW_BW::mumsFittableProperties()";
   }
   return *_fittableResonancePropertiesPtr;
 }

◆ mumsMass()

double BW_BW::mumsMass ( ) const

protectedvirtual

Definition at line 453 of file BW_BW.cpp.

                             {
   bool dbThis=false;
   return mumsFittableProperties().mass();
   // return _RPL->get(mumsPID())->mass();
   (void)dbThis;
 }

◆ mumsParity()

int BW_BW::mumsParity ( ) const

protectedvirtual

Definition at line 238 of file BW_BW.cpp.

                            {
   if( 0 == _mumsParity ){
     const ParticleProperties* pp = _theDecay.getVal().props();
     if(0==pp){
       cout << "ERROR in BW_BW::mumsParity()"
        << " can't find properties for first element"
        << " in this decay tree\n" << _theDecay << endl;
       throw "invalid decay tree in BW_BW::mumsParity()";
     }
     _mumsParity = pp->P();
   }
   return _mumsParity;
 }

◆ mumsPID()

int BW_BW::mumsPID ( ) const

protectedvirtual

Definition at line 517 of file BW_BW.cpp.

                         {
   if(! _mumsPID_set){
     _mumsPID = mumsProperties()->pdg_id();
     _mumsPID_set = true;
   }
   return _mumsPID;
 }

◆ mumsProperties()

const ParticleProperties * BW_BW::mumsProperties ( ) const

protectedvirtual

Reimplemented in BW_BW_DifferentMother.

Definition at line 380 of file BW_BW.cpp.

                                                      {
   const ParticleProperties* pp = _theDecay.getVal().props();
   if(0==pp){
     cout << "ERROR in BW_BW::mumsProperties()"
      << " can't find properties for first element"
      << " in this decay tree\n" << _theDecay << endl;
     throw "invalid decay tree in BW_BW::mumsProperties()";
   }
   return pp;
 }

◆ mumsQuarkContent()

MultiQuarkContent BW_BW::mumsQuarkContent ( ) const

protectedvirtual

Definition at line 215 of file BW_BW.cpp.

                                                {
   return netQuarkContent(_theDecay); // handles non-res
 
   //return mumsProperties()->netQuarkContent();
 }

◆ mumsRadius()

double BW_BW::mumsRadius ( ) const

protectedvirtual

Definition at line 465 of file BW_BW.cpp.

                               {
   return mumsFittableProperties().radius();
   //    return _RPL->get(mumsPID())->radius();
 }

◆ mumsRecoMass()

double BW_BW::mumsRecoMass ( ) const

protectedvirtual

Definition at line 631 of file BW_BW.cpp.

                                 {
   if(substitutePDGForReco()){
     return mumsMass();
   }
   if(_mumsRecoMass < 0){
     double m2 = mumsRecoMass2();
     if(m2 < 0){
       cout << " ERROR negative reconstructed mass2"
        << " in decay:\n " << _theDecay << endl;
       _mumsRecoMass = - sqrt(-m2);
     }else{
       _mumsRecoMass = sqrt(m2);
     }
   }
   return _mumsRecoMass;
 }

◆ mumsRecoMass2()

double BW_BW::mumsRecoMass2 ( ) const

protectedvirtual

Definition at line 548 of file BW_BW.cpp.

                                  {
   if(substitutePDGForReco()) return mumsMass()*mumsMass();
 
   if(_mumsRecoMass2 < 0){
     std::vector<int> asi = _theDecay.getVal().asi();
     if(asi.size() < 2){
       cout << "ERROR in BW_BW::mumsRecoMass2() "
        << _theDecay.getVal()._pdg_id
        << ", " << _theDecay.getVal().name()
        << "  decays into " << asi.size() << " particles?\n" 
        << _theDecay
        << endl;
     }
     _mumsRecoMass2 = getEvent()->sij(asi);
   }
   return _mumsRecoMass2;
 }

◆ mumsSpin()

std::string BW_BW::mumsSpin ( ) const

protectedvirtual

Definition at line 509 of file BW_BW.cpp.

                                {
   return mumsProperties()->J();
 }

◆ mumsSpinValue()

int BW_BW::mumsSpinValue ( ) const

protectedvirtual

Definition at line 513 of file BW_BW.cpp.

                               {
   return atoi(mumsSpin().c_str());
 }

◆ mumsWidth()

double BW_BW::mumsWidth ( ) const

protectedvirtual

Definition at line 460 of file BW_BW.cpp.

                              {
   return mumsFittableProperties().width();
   //    return _RPL->get(mumsPID())->width();
 }

◆ name()

virtual std::string BW_BW::name ( ) const

inlinevirtual

Implements ILineshape.

Reimplemented in FocusFlatte, Flatte, CrystalBallRhoOmega, Histo_BW, Lass, GounarisSakurai, GLass, CrystalBarrelFOCUS, FermiPS_BW, NonRes, Rho0Omega, Bugg_BW, RhoOmegaGS, Model_independent, and SBW.

Definition at line 198 of file BW_BW.h.

                                 {
     return "BW_BW("+_theDecay.oneLiner() +")";
   }

◆ nonResonant()

bool BW_BW::nonResonant ( ) const

protected

Definition at line 221 of file BW_BW.cpp.

                              {
   return mumsProperties()->isNonResonant();
 }

◆ numDaughters()

int BW_BW::numDaughters ( ) const

protectedvirtual

Definition at line 252 of file BW_BW.cpp.

                              {
   return _theDecay.nDgtr();
 }

◆ pABSq()

double BW_BW::pABSq ( )

virtual

Definition at line 1014 of file BW_BW.cpp.

                    {
   bool dbThis=false;
   if(_pABSq < 0 || dbThis){
     _pABSq = twoBody_dgtPsq_in_MumsRecoFrame();
     if(_pABSq < 0 || dbThis){ //dbg
       cout << "In " << _theDecay.oneLiner()
        << " BW_BW::pABSq()"
        << " suspicious... pABSq < 0. This is for:"
        << "\n   mumM = " << mumsRecoMass() 
        << "\n, m1 " << daughterRecoMass(0)
        << ", p1 " << daughterP4(0) << ", " << daughterP4(0).M()
        << "\n , m2 " << daughterRecoMass(1)
        << ", p2 " << daughterP4(1) << ", " << daughterP4(1).M()
        << " returnVal " << _pABSq
        << endl;
     }
     if(_pABSq < 0) _pABSq = 0;
   }
   return _pABSq;
 }

◆ parityConservingL()

bool BW_BW::parityConservingL ( int L ) const

protectedvirtual

Definition at line 233 of file BW_BW.cpp.

                                         {
   int PL = (L%2 == 0 ? 1 : -1);
   return ((PL * dgtrsInternalParity()) == mumsParity());
 }

◆ peakPosition()

double BW_BW::peakPosition ( ) const

private

Definition at line 1237 of file BW_BW.cpp.

                                 {
   // from: J.D. Jackson, Il Nuovo Cimento, Siere X, Vol 34, pag 1644-1666
   // 16 Dicembre 1964
   // equation 10. (refers to (4)).
   //
   return mumsMass() + peakShift();
 }

◆ peakShift()

double BW_BW::peakShift ( ) const

private

Definition at line 1207 of file BW_BW.cpp.

                              {
   // from: J.D. Jackson, Il Nuovo Cimento, Siere X, Vol 34, pag 1644-1666
   // 16 Dicembre 1964
   // equation 10. (refers to (4)).
   //
   double m1 = daughterPDGMass(0);
   double m2 = daughterPDGMass(1);
 
   double m1s = m1*m1;
   double m2s = m2*m2;
 
   double om = mumsMass();
   double om2 = om*om;
   double om4 = om2*om2;
 
   double G  = mumsWidth();
 
   double mpm = m1s + m2s;
   double mmm = m1s - m2s;
   double mmm2 = mmm*mmm;
 
   double num = om4 - mmm2;
   double den = om4 - 2*mpm*om2 + mmm2;
 
   double f1 = ((double) twoLPlusOne())/8.0;
   double f2 = G*G/om;
   double f3 = num/den;
   
   return f1*f2*f3;
 }

◆ prefix()

const std::string& BW_BW::prefix ( ) const

inlineprotected

Definition at line 72 of file BW_BW.h.

72 {return _prefix;}

BW_BW::_prefix

std::string _prefix

Definition: BW_BW.h:64

◆ print() [1/2]

void BW_BW::print	(	IDalitzEvent &	evt,
		std::ostream &	out = `std::cout`
	)

virtual

Reimplemented in CrystalBallRhoOmega, CrystalBarrelFOCUS, NonRes, and Rho0Omega.

Definition at line 1287 of file BW_BW.cpp.

                                                     {
   setEventPtr(evt);
   resetInternals();
   out << name()
       << "\n\t> co-ordinate: " << getDalitzCoordinate()
       << "\n\t> This is the decay I'm looking at:"
       << "\n" << _theDecay
       << "\n\t> These are a few values: "
       << " startOfDecayChain? " << startOfDecayChain()
       << ", BreitWigner = " 
       << BreitWigner()
       << ", Blatt-Weisskopf penetration factor: "
       << Fr()
       << ", total BW_BW: " 
       << getVal(evt);
 }

◆ print() [2/2]

void BW_BW::print ( std::ostream & out = std::cout ) const

virtual

Implements ILineshape.

Reimplemented in CrystalBallRhoOmega, CrystalBarrelFOCUS, NonRes, and Rho0Omega.

Definition at line 1284 of file BW_BW.cpp.

                                       {
   out << name();
 }

◆ prSq()

double BW_BW::prSq ( ) const

virtual

Definition at line 969 of file BW_BW.cpp.

                         {
   bool dbThis=false;
   if(_prSq <  0){
     _prSq = twoBody_dgtPsq_in_MumsPDGFrame();
     if(dbThis && _prSq < 0){ //dbg
       cout << " BW_BW::prSq()"
        << " suspicious... pr < 0. This is for:"
        << "\n   mumM PDG " << mumsMass()
        << ", m1_nominal " << daughterPDGMass(0)
        << ", m2_nominal " << daughterPDGMass(1)
        << ", m1_reco " << daughterRecoMass(0)
        << ", m2_reco " << daughterRecoMass(1)
        << " _prSq " << _prSq
        << " will return 0 "
        << endl;
     }
   }
   if(_prSq < 0) _prSq = 0;
   return _prSq;
 }

◆ prSqForGofM()

double BW_BW::prSqForGofM ( ) const

virtual

Definition at line 989 of file BW_BW.cpp.

                                {
   bool dbThis=false;
   if(DoAsLaurenDid) return prSq();
 
   if(_prSqForGofM <  0){
     _prSqForGofM = fabs(twoBody_recodgtPsq_in_MumsPDGFrame()); // forces > 0.
   }
   if(dbThis || _prSqForGofM < 0){ //dbg
     cout << " compare : " << _prSqForGofM
      << ", " << fabs(twoBody_recodgtPsq_in_MumsPDGFrame())
      << endl;
     cout << " BW_BW::prSqForGofM()"
      << " suspicious... pr < 0. This is for:"
      << "\n   mumM PDG " << mumsMass()
      << ", m1_nominal " << daughterPDGMass(0)
      << ", m2_nominal " << daughterPDGMass(1)
      << ", m1_reco " << daughterRecoMass(0)
      << ", m2_reco " << daughterRecoMass(1)
      << " _prSq " << _prSqForGofM
      << " will return 0 "
      << endl;
   }
   if(_prSqForGofM < 0) _prSqForGofM = 0; // can not happen as I'm using fabs above...
   return _prSqForGofM;
 }

◆ prSqMax()

double BW_BW::prSqMax ( ) const

virtual

Definition at line 944 of file BW_BW.cpp.

                            {
 
 
   std::vector<const AssociatedDecayTreeItem*> dgtrOne = 
     _theDecay.getDgtrTreePtr(0)->finalState();
   double m1_min=0;
   for(unsigned int i=0; i < dgtrOne.size(); i++) m1_min += dgtrOne[i]->mass();
 
   std::vector<const AssociatedDecayTreeItem*> dgtrTwo = 
     _theDecay.getDgtrTreePtr(1)->finalState();
   double m2_min=0;
   for(unsigned int i=0; i < dgtrTwo.size(); i++) m2_min += dgtrTwo[i]->mass();
 
   /*
   double m1 = daughterPDGMass(0) - 0.5*daughterWidth(0);
   if(m1 < m1_min) m1 = m1_min;
 
   double m2 = daughterPDGMass(1) - 0.5*daughterWidth(1);
   if(m2 < m2_min) m2 = m2_min;
   */
 
   return twoBody_dgtPsq_in_MumsFrame(mumsMass(), m1_min, m2_min);
 
   //  return mumsPDGMass()*mumsPDGMass()/4.0;
 }

◆ Radius()

double BW_BW::Radius ( ) const

protectedvirtual

Definition at line 481 of file BW_BW.cpp.

                           {
     if(_useGlobalRadius) return globalRadius();
     else return mumsRadius();
 }

◆ resetInternals()

void BW_BW::resetInternals ( )

protectedvirtual

Definition at line 1114 of file BW_BW.cpp.

                           {
   _prSq =  _prSqForGofM = _pABSq = _mumsRecoMass2 =  _mumsRecoMass = _Fr_BELLE = _Fr_PDG_BL =  _GofM = -1;
 
   _daughterRecoMass2.resize(_theDecay.nDgtr());
   for(int i=0; i< _theDecay.nDgtr(); i++){
     _daughterRecoMass2[i] = -1;
   }
 }

◆ resetPDG()

void BW_BW::resetPDG ( )

protectedvirtual

Definition at line 1123 of file BW_BW.cpp.

                     {
   //  _mumsPDGMass = _mumsPDGWidth = -1;
   _mumsPID_set = false;
   //  _mumsPDGRadius  = -9999.0;
 
   _daughterPDGMass.resize(_theDecay.nDgtr());
   for(int i=0; i< _theDecay.nDgtr(); i++){
     _daughterPDGMass[i] = -1;
   }
   _daughterWidth.resize(_theDecay.nDgtr());
   for(int i=0; i< _theDecay.nDgtr(); i++){
     _daughterWidth[i] = -1;
   }
 }

◆ resonanceProperties()

const ResonanceProperties * BW_BW::resonanceProperties ( ) const

protected

Definition at line 403 of file BW_BW.cpp.

                                                            {
   const ResonanceProperties* rp = resonancePropertiesList()->AddToListIfMissing(mumsPID());
   if(0 == rp){
     cout << "ERROR in BW_BW::ResonanceProperties()"
      << " can't find properties for first element"
      << " in this decay tree\n" << _theDecay << endl;
     throw "invalid decay tree in BW_BW::mumsFittableProperties()";
   }
   return rp;
 }

◆ resonancePropertiesList()

ResonancePropertiesList * BW_BW::resonancePropertiesList ( ) const

protected

Definition at line 391 of file BW_BW.cpp.

                                                               {
   //cout << "in BW_BW::resonancePropertiesList(): prefix() = " << prefix() << endl;
   if(0 == _RPL) _RPL = ResonancePropertiesList::getMe(prefix(), _mps);
   if(0 == _RPL){
     cout << "ERROR in BW_BW::resonanceProperties()"
      << " can't find properties for first element"
      << " in this decay tree\n" << _theDecay << endl;
     throw "invalid decay tree in BW_BW::resonanceProperties()";
   }
   return _RPL;
 }

◆ setAllFitParameters()

bool BW_BW::setAllFitParameters ( )

protected

Definition at line 426 of file BW_BW.cpp.

                                {
   bool dbThis = false;
   if(dbThis)cout << "BW_BW::setAllFitParameters() called" << endl;
   bool s=true;
   if(0 == resonancePropertiesList()){
     cout << "big problem in BW_BW::setAllFitParameters"
      << ", resonancePropertiesList is zero" << endl;
     throw "BW_BW::setAllFitParameters can't find resonancePropertiesList()";
   }
   if(dbThis)cout << "getting _fittableGlobalRadiusPtr" << endl;
   const FitParameter& rad(resonancePropertiesList()->radius());
   if(dbThis)cout << "called radius, no crash, now pringint it:" << endl;
   if(dbThis)cout << "rad is" << rad << endl;
   if(dbThis)cout << " radius is: " << resonancePropertiesList()->radius() << endl;
   _fittableGlobalRadiusPtr = new FitParRef(resonancePropertiesList()->radius(), this);
   if(dbThis)cout << "hurray, got _fittableGlobalRadiusPtr" << endl;
   s &= (0 != _fittableGlobalRadiusPtr);
   _fittableResonancePropertiesPtr = new ResonancePropertiesFitRef(*resonanceProperties(), this);
   s &= (0 != _fittableResonancePropertiesPtr);
   NamedParameter<int> useGlobalRadius(("UseGlobalRadiusFor_"+ resonanceProperties()->nameFromPid(mumsPID()) ).c_str(),1,NamedParameterBase::QUIET);
   if(useGlobalRadius == 0) _useGlobalRadius = 0;  
   if(dbThis && _useGlobalRadius)cout << "I'll use the global radius for this decay tree\n" << _theDecay << endl;
   if(dbThis && !_useGlobalRadius)cout << "I'll use the individual radius for this decay tree\n" << _theDecay << endl;
   return s;
 }

◆ setEventPtr()

bool BW_BW::setEventPtr ( IDalitzEvent & evt ) const

protected

Definition at line 141 of file BW_BW.cpp.

                                               {
   //_oldPointers.push(getEvent());
   _eventPtr = &(evt);
   return true;
 }

◆ setGenerationLimits()

void BW_BW::setGenerationLimits	(	double	mi,
		double	ma
	)

virtual

Definition at line 1109 of file BW_BW.cpp.

                                                    {
   _gen_s_mi = mi;
   _gen_s_ma = ma;
 }

◆ startOfDecayChain()

bool BW_BW::startOfDecayChain ( ) const

protectedvirtual

Definition at line 532 of file BW_BW.cpp.

                                    {
   return ! (_theDecay.hasParent());
 }

◆ substitutePDGForReco()

bool BW_BW::substitutePDGForReco ( ) const

inlineprotected

Definition at line 73 of file BW_BW.h.

                                    {
     return _substitutePDGForReco;
   }

◆ twoBody_dgtPsq_in_MumsFrame()

double BW_BW::twoBody_dgtPsq_in_MumsFrame	(	double	mumsMass,
		double	mA,
		double	mB
	)		const

protectedvirtual

Definition at line 917 of file BW_BW.cpp.

                                        {
   bool dbThis=false;
 
   double Msq = mumsMass*mumsMass;
   if(Msq <=0) return -9999;
   
   double mpsq = (mA + mB)*(mA + mB);
   double mmsq = (mA - mB)*(mA - mB);
   double num = (Msq - mpsq)*(Msq - mmsq);
   
   double returnVal = num/(4*Msq);
 
   if(dbThis){
     double mAB = mumsMass;
     double EvtGenVal = (((mAB*mAB-mA*mA-mB*mB)*(mAB*mAB-mA*mA-mB*mB)/4.0) -
             mA*mA*mB*mB)/(mAB*mAB);
     cout << " BW_BW::twoBody_dgtPsq_in_MumsFrame"
      << " compare: cleo: " << EvtGenVal
      << " MINT: " << returnVal
      << endl;
   }
   return returnVal;
   
 }

◆ twoBody_dgtPsq_in_MumsPDGFrame()

double BW_BW::twoBody_dgtPsq_in_MumsPDGFrame ( ) const

protectedvirtual

Definition at line 889 of file BW_BW.cpp.

                                                   {
   if(! PDGWithReco){
     return twoBody_dgtPsq_in_MumsFrame(mumsMass()
                        , daughterPDGMass(0)
                        , daughterPDGMass(1)
                        );
   }else{
     return twoBody_dgtPsq_in_MumsFrame(mumsMass()
                        , daughterRecoMass(0)
                        , daughterRecoMass(1)
                        );
   }
 }

◆ twoBody_dgtPsq_in_MumsRecoFrame()

double BW_BW::twoBody_dgtPsq_in_MumsRecoFrame ( )

protectedvirtual

Definition at line 910 of file BW_BW.cpp.

                                              {
   return twoBody_dgtPsq_in_MumsFrame(mumsRecoMass()
                      , daughterRecoMass(0)
                      , daughterRecoMass(1)
                      );
 }

◆ twoBody_recodgtPsq_in_MumsPDGFrame()

double BW_BW::twoBody_recodgtPsq_in_MumsPDGFrame ( ) const

protectedvirtual

Definition at line 903 of file BW_BW.cpp.

                                                       {
   return twoBody_dgtPsq_in_MumsFrame(mumsMass()
                      , daughterRecoMass(0)
                      , daughterRecoMass(1)
                      );
 }

◆ twoLPlusOne()

int BW_BW::twoLPlusOne ( ) const

virtual

Definition at line 151 of file BW_BW.cpp.

                             {
   if( _twoLPlusOne < 0 ){
     int externallySetL = _theDecay.getVal().L();
     if( externallySetL >= 0){
       _twoLPlusOne = 2* externallySetL + 1;
     }else{
       _twoLPlusOne = lowestPossibleTwoLPlusOne();
     }
   }
   return _twoLPlusOne;
 }

Member Data Documentation

◆ _daughterPDGMass

std::vector<double> BW_BW::_daughterPDGMass

mutableprivate

Definition at line 47 of file BW_BW.h.

◆ _daughterRecoMass2

std::vector<double> BW_BW::_daughterRecoMass2

mutableprivate

Definition at line 47 of file BW_BW.h.

◆ _daughterWidth

std::vector<double> BW_BW::_daughterWidth

mutableprivate

Definition at line 47 of file BW_BW.h.

◆ _dgtrsInternalParity

int BW_BW::_dgtrsInternalParity

mutableprivate

Definition at line 42 of file BW_BW.h.

◆ _eventPtr

IDalitzEvent* BW_BW::_eventPtr

mutableprivate

Definition at line 32 of file BW_BW.h.

◆ _fittableGlobalRadiusPtr

MINT::FitParRef* BW_BW::_fittableGlobalRadiusPtr

protected

Definition at line 99 of file BW_BW.h.

◆ _fittableResonancePropertiesPtr

ResonancePropertiesFitRef* BW_BW::_fittableResonancePropertiesPtr

protected

Definition at line 98 of file BW_BW.h.

◆ _Fr_BELLE

double BW_BW::_Fr_BELLE

mutableprivate

Definition at line 35 of file BW_BW.h.

◆ _Fr_PDG_BL

double BW_BW::_Fr_PDG_BL

mutableprivate

Definition at line 35 of file BW_BW.h.

◆ _gen_s_ma

double BW_BW::_gen_s_ma

protected

Definition at line 70 of file BW_BW.h.

◆ _gen_s_mi

double BW_BW::_gen_s_mi

protected

Definition at line 70 of file BW_BW.h.

◆ _genFct

MINT::counted_ptr<IGenFct> BW_BW::_genFct

mutableprivate

Definition at line 51 of file BW_BW.h.

◆ _GofM

double BW_BW::_GofM

mutableprivate

Definition at line 35 of file BW_BW.h.

◆ _mps

MINT::MinuitParameterSet* BW_BW::_mps

mutableprotected

Definition at line 59 of file BW_BW.h.

◆ _mumsParity

int BW_BW::_mumsParity

mutableprivate

Definition at line 42 of file BW_BW.h.

◆ _mumsPID

int BW_BW::_mumsPID

mutableprivate

Definition at line 37 of file BW_BW.h.

◆ _mumsPID_set

bool BW_BW::_mumsPID_set

mutableprivate

Definition at line 38 of file BW_BW.h.

◆ _mumsRecoMass

double BW_BW::_mumsRecoMass

mutableprivate

Definition at line 35 of file BW_BW.h.

◆ _mumsRecoMass2

double BW_BW::_mumsRecoMass2

mutableprivate

Definition at line 35 of file BW_BW.h.

◆ _normBF

MINT::NamedParameter<int> BW_BW::_normBF

protected

Definition at line 66 of file BW_BW.h.

◆ _pABSq

double BW_BW::_pABSq

mutableprivate

Definition at line 35 of file BW_BW.h.

◆ _prefix

std::string BW_BW::_prefix

protected

Definition at line 64 of file BW_BW.h.

◆ _prSq

double BW_BW::_prSq

mutableprivate

Definition at line 35 of file BW_BW.h.

◆ _prSqForGofM

double BW_BW::_prSqForGofM

mutableprivate

Definition at line 35 of file BW_BW.h.

◆ _RPL

ResonancePropertiesList* BW_BW::_RPL

mutableprotected

Definition at line 97 of file BW_BW.h.

◆ _substitutePDGForReco

bool BW_BW::_substitutePDGForReco

mutableprivate

Definition at line 49 of file BW_BW.h.

◆ _theDecay

const AssociatedDecayTree& BW_BW::_theDecay

protected

Definition at line 69 of file BW_BW.h.

◆ _twoLPlusOne

int BW_BW::_twoLPlusOne

mutableprivate

Definition at line 44 of file BW_BW.h.

◆ _useGlobalRadius

bool BW_BW::_useGlobalRadius

protected

Definition at line 67 of file BW_BW.h.

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

Mint/BW_BW.h
src/Mojito/Lineshapes/BW_BW.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ BW_BW() [1/2]

◆ BW_BW() [2/2]

◆ ~BW_BW()

Member Function Documentation

◆ BreitWigner()

◆ daughterP()

◆ daughterP4()

◆ daughterPDGMass()

◆ daughterQuarkContent()

◆ daughterRecoMass()

◆ daughterRecoMass2()

◆ daughterSpin()

◆ daughterSpinValue()

◆ daughterWidth()

◆ dgtrsInternalParity()

◆ dgtrsQuarkContent()

◆ EvtGenValue()

◆ Fr()

◆ Fr_BELLE()

◆ Fr_BELLE_Max()

◆ Fr_PDG_BL()

◆ FrForGofM()

◆ FrMax()

◆ generatingFunction()

◆ GetAlpha()

◆ getDalitzCoordinate()

◆ getEvent()

◆ getMinuitParameterSet() [1/2]

◆ getMinuitParameterSet() [2/2]

◆ getVal()

◆ globalRadius()

◆ GofM()

◆ isWeakDecay()

◆ lowestPossibleTwoLPlusOne()

◆ makeGeneratingFunction()

◆ maxDaughterPairSpinSum()

◆ maxDaughterSpinSum()

◆ minDaughterPairSpinSum()

◆ minDaughterSpinSum()

◆ minDaughterSpinSum2()

◆ minDaughterSpinSum3()

◆ minDaughterSpinSum4()

◆ mumsFittableProperties()

◆ mumsMass()

◆ mumsParity()

◆ mumsPID()

◆ mumsProperties()

◆ mumsQuarkContent()

◆ mumsRadius()

◆ mumsRecoMass()

◆ mumsRecoMass2()

◆ mumsSpin()

◆ mumsSpinValue()

◆ mumsWidth()

◆ name()

◆ nonResonant()

◆ numDaughters()

◆ pABSq()

◆ parityConservingL()

◆ peakPosition()

◆ peakShift()

◆ prefix()

◆ print() [1/2]

◆ print() [2/2]

◆ prSq()

◆ prSqForGofM()

◆ prSqMax()

◆ Radius()

◆ resetInternals()

◆ resetPDG()

◆ resonanceProperties()

◆ resonancePropertiesList()

◆ setAllFitParameters()