#include <SpinFactors4Body_Tensors.h>

Inheritance diagram for SF_DtoTP0_TtoVP1_VtoP2P3:

Public Member Functions
virtual void	printYourself (std::ostream &os=std::cout) const

	SF_DtoTP0_TtoVP1_VtoP2P3 (const AssociatedDecayTree &theDecay)

virtual	~SF_DtoTP0_TtoVP1_VtoP2P3 ()

virtual double	getVal (IDalitzEvent &evt)

virtual const DecayTree &	exampleDecay ()

virtual std::string	name () const

Public Member Functions inherited from SpinFactor
virtual double	RealVal (IDalitzEvent &evt)

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

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

double	mRes (const AssociatedDecayTreeItem &adt, IDalitzEvent &evt)

double	mRes (const MINT::const_counted_ptr< AssociatedDecayTree > &adt, IDalitzEvent &evt)

const AssociatedDecayTree &	theDecay () const

const AssociatedDecayTree &	theDecay (const DalitzEventPattern &) const

const AssociatedDecayTree &	theDecay (IDalitzEvent &evt) const

	SpinFactor (const AssociatedDecayTree &decay, int nFinal=3)

	SpinFactor (const SpinFactor &other)

virtual void	printParsing (std::ostream &os=std::cout) const

virtual	~SpinFactor ()

Public Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
virtual	~IReturnRealForEvent ()

Public Member Functions inherited from MINT::IReturnComplexForEvent< IDalitzEvent >
virtual	~IReturnComplexForEvent ()

Public Member Functions inherited from CachedByEvent< std::complex< double > >
	CachedByEvent ()

	CachedByEvent (const CachedByEvent &)

std::complex< double >	getValWithCaching (IDalitzEvent &evt)

std::complex< double >	getValWithCachingPermutation (IDalitzEvent &evt)

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

Public Member Functions inherited from ISpinFactor
virtual	~ISpinFactor ()

Static Public Member Functions
static const DecayTree &	getExampleDecay ()

Protected Member Functions
virtual bool	parseTree (const DalitzEventPattern &pat)

Protected Member Functions inherited from SpinFactor
TLorentzVector	p (int i, IDalitzEvent &evt)

Protected Member Functions inherited from MINT::IReturnRealForEvent< IDalitzEvent >
	IReturnRealForEvent ()

Protected Member Functions inherited from MINT::IReturnComplexForEvent< IDalitzEvent >
	IReturnComplexForEvent ()

Protected Member Functions inherited from CachedByEvent< std::complex< double > >
long int	rememberNumber ()

long int	rememberNumberPermutation (IDalitzEvent &evt)

long int	configNumber ()

std::complex< double >	recalculate (IDalitzEvent &evt)

std::complex< double >	recalculatePermutation (IDalitzEvent &evt)

Protected Member Functions inherited from ISpinFactor
	ISpinFactor ()

Protected Attributes
MINT::const_counted_ptr< AssociatedDecayTree >	T

MINT::const_counted_ptr< AssociatedDecayTree >	V

Protected Attributes inherited from SpinFactor
const AssociatedDecayTree &	_theDecay

int	_nFinal

std::vector< MINT::const_counted_ptr< AssociatedDecayTree > >	fsPS

MINT::NamedParameter< int >	_useZemachTensors

Protected Attributes inherited from CachedByEvent< std::complex< double > >
long int	_rememberNumber

std::vector< long int >	_rememberNumberPermutation

long int	_configNumber

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

Static Protected Attributes
static DecayTree *	_exampleDecay =0

Detailed Description

Definition at line 201 of file SpinFactors4Body_Tensors.h.

Constructor & Destructor Documentation

◆ SF_DtoTP0_TtoVP1_VtoP2P3()

SF_DtoTP0_TtoVP1_VtoP2P3::SF_DtoTP0_TtoVP1_VtoP2P3 ( const AssociatedDecayTree & theDecay )

inline

Definition at line 210 of file SpinFactors4Body_Tensors.h.

      : SpinFactor(theDecay, 4)
      , T(0), V(0)
      {
      }

◆ ~SF_DtoTP0_TtoVP1_VtoP2P3()

virtual SF_DtoTP0_TtoVP1_VtoP2P3::~SF_DtoTP0_TtoVP1_VtoP2P3 ( )

inlinevirtual

Definition at line 216 of file SpinFactors4Body_Tensors.h.

216 {}

Member Function Documentation

◆ exampleDecay()

const DecayTree & SF_DtoTP0_TtoVP1_VtoP2P3::exampleDecay ( )

virtual

Implements SpinFactor.

Definition at line 123 of file SpinFactors4Body_Tensors.cpp.

                                                        {
   return getExampleDecay();
 }

◆ getExampleDecay()

const DecayTree & SF_DtoTP0_TtoVP1_VtoP2P3::getExampleDecay ( )

static

Definition at line 110 of file SpinFactors4Body_Tensors.cpp.

                                                           {
   _exampleDecay = new DecayTree(421);
   _exampleDecay->addDgtr(-211, 215)->addDgtr(211, 113)->addDgtr(-211, 211);
   return *_exampleDecay;
 }

◆ getVal()

double SF_DtoTP0_TtoVP1_VtoP2P3::getVal ( IDalitzEvent & evt )

virtual

Implements SpinFactor.

Definition at line 597 of file SpinFactors4Body_Tensors.cpp.

                                                         {
   bool dbThis=false;
   if(! ( fsPS[0] && fsPS[1] && fsPS[2] && fsPS[3]) ) parseTree(evt.eventPattern());
   
   TLorentzVector pV = p(2, evt) + p(3, evt);
   TLorentzVector qV = p(2, evt) - p(3, evt);
   TLorentzVector pT = pV + p(1, evt);
   TLorentzVector qT = pV - p(1, evt);
   TLorentzVector pD = pT + p(0, evt);
   TLorentzVector qD = pT - p(0, evt);
 
   double MT = mRes(T, evt);
   double MV = mRes(V, evt);
   
   ZTspin2 t2T(qT, pT, MT);
   ZTspin1 tV(qV, pV, MV);
   
   TLorentzVector DT(t2T.Contract(qD));
 
   const double units = GeV*GeV * GeV*GeV * GeV*GeV;
 
   if (_useZemachTensors) {
             ZTspin1 tD(qD,pD,pD.M());
             ZTspin1 t1T(qT, pT, MT);
             SpinSumV P1T(pT,MT);
             SpinSumT P2T(pT,MT);
       
             TLorentzVector tmp = LeviCivita(t1T,pT,P1T.Dot(tV));
             
             return P2T.Sandwich(tD,tD,tmp,t1T)/ units + 1./3. * tD.Dot(tD) / (pT.Dot(pT)) * P2T.Sandwich(pD,pD,tmp,t1T)/ units;
   }  
     
   double returnVal = LeviCivita(pD, qD, DT, tV) / units;
 
   if(dbThis){
     cout << " SF_DtoTP0_TtoVP1_VtoP2P3::getVal "
      << " returning " << returnVal
      << endl;
   }
   return returnVal;
 }

◆ name()

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

inlinevirtual

Implements SpinFactor.

Definition at line 221 of file SpinFactors4Body_Tensors.h.

                                  {
      return "SpinFactor4:SF_DtoTP0_TtoVP1_VtoP2P3("
        + theDecay().oneLiner() + ")";
    }

◆ parseTree()

bool SF_DtoTP0_TtoVP1_VtoP2P3::parseTree ( const DalitzEventPattern & pat )

protectedvirtual

Definition at line 544 of file SpinFactors4Body_Tensors.cpp.

                                                                      {
   //  bool debugThis=false;
   if(fsPS.size() < 4) fsPS.reserve(4);
   if(theDecay(pat).nDgtr() != 2){
     cout << "ERROR in SF_DtoTP0_TtoVP1_VtoP2P3_BASE::parseTree"
      << " expected exactly 2 daughers of D, have "
      << theDecay(pat).nDgtr();
     return false;
   }
   for(int i=0; i< theDecay(pat).nDgtr(); i++){
     const_counted_ptr<AssociatedDecayTree> dgtr= theDecay(pat).getDgtrTreePtr(i);
     if     (dgtr->getVal().SVPAT() == "P" &&   dgtr->isFinalState()) fsPS[0] = dgtr;
     else if(dgtr->getVal().SVPAT() == "T" && ! dgtr->isFinalState()) T = dgtr;
   }
   if(0==T || 0==fsPS[0]){
     cout << "ERROR in SF_DtoTP0_TtoVP1_VtoP2P3_BASE::parseTree"
      << " Didn't find T or P0 " << T.get() << ", " << fsPS[0].get() << endl;
     return false;
   }
   
   if(T->nDgtr() != 2){
     cout << "ERROR in SF_DtoTP0_TtoVP1_VtoP2P3_BASE::parseTree"
      << " T should have 2 daughters, but it says it has "
      << T->nDgtr() << "."
      << endl;
     return false;
   }
   for(int i=0; i< T->nDgtr(); i++){
     const_counted_ptr<AssociatedDecayTree> dgtr= T->getDgtrTreePtr(i);
     if     (dgtr->getVal().SVPAT() == "V" && ! dgtr->isFinalState()) V = dgtr;
     else if(dgtr->getVal().SVPAT() == "P" &&   dgtr->isFinalState()) fsPS[1] = dgtr;
   }
   if(0==V || 0==fsPS[1]){
     cout << "ERROR in SF_DtoTP0_TtoVP1_VtoP2P3_BASE::parseTree"
      << " Didn't find V2 or P1 " << V.get() << ", " << fsPS[1].get() << endl;
     return false;
   }
   
   if(V->nDgtr() != 2){
     cout << "ERROR in SF_DtoTP0_TtoVP1_VtoP2P3_BASE::parseTree"
      << " V should have 2 daughters, but it says it has "
      << V->nDgtr() << "."
      << endl;
     return false;
   }
   fsPS[2] = V->getDgtrTreePtr(0);
   fsPS[3] = V->getDgtrTreePtr(1);
   // normalOrder(fsPS[2], fsPS[3]);
   
   // this->printYourself();
   return true;
 }

◆ printYourself()

void SF_DtoTP0_TtoVP1_VtoP2P3::printYourself ( std::ostream & os = std::cout ) const

virtual

Implements SpinFactor.

Definition at line 639 of file SpinFactors4Body_Tensors.cpp.

                                                              {
   //  bool debugThis = false;
   if(! ( fsPS[0] && fsPS[1] && fsPS[2] && fsPS[3]) ) return;
   os << "spin factor SF_DtoTP0_TtoVP1_VtoP2P3"
      << "\n\t  ZTspin2 tT1(qT, pT, MT);"
      << "\n\t  ZTspin1 tT2(qV, pV, MV);"
      << "\n\t  pV = p(2) + p(3); qV = p(2) - p(3);"
      << "\n\t  pT = pV + p(1); qT = pV - p(1);"
      << "\n\t  pD = pT + p(0); qD = pT - p(0);"
      << "\n\t  return: DT(tT.Contract(qD))/GeV^6"
      << "\n\t    parsed tree " << theDecay().oneLiner()
      << "\n      like this:" << endl;
   this->printParsing(os);
 }