#include <SpinFactors4Body_Tensors.h>

Inheritance diagram for SF_DtoTS_TtoP0P1_StoP2P3:

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

	SF_DtoTS_TtoP0P1_StoP2P3 (const AssociatedDecayTree &theDecay)

virtual	~SF_DtoTS_TtoP0P1_StoP2P3 ()

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 >	S

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 258 of file SpinFactors4Body_Tensors.h.

Constructor & Destructor Documentation

◆ SF_DtoTS_TtoP0P1_StoP2P3()

SF_DtoTS_TtoP0P1_StoP2P3::SF_DtoTS_TtoP0P1_StoP2P3 ( const AssociatedDecayTree & theDecay )

inline

Definition at line 267 of file SpinFactors4Body_Tensors.h.

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

◆ ~SF_DtoTS_TtoP0P1_StoP2P3()

virtual SF_DtoTS_TtoP0P1_StoP2P3::~SF_DtoTS_TtoP0P1_StoP2P3 ( )

inlinevirtual

Definition at line 273 of file SpinFactors4Body_Tensors.h.

273 {}

Member Function Documentation

◆ exampleDecay()

const DecayTree & SF_DtoTS_TtoP0P1_StoP2P3::exampleDecay ( )

virtual

Implements SpinFactor.

Definition at line 140 of file SpinFactors4Body_Tensors.cpp.

                                                        {
   return getExampleDecay();
 }

◆ getExampleDecay()

const DecayTree & SF_DtoTS_TtoP0P1_StoP2P3::getExampleDecay ( )

static

Definition at line 129 of file SpinFactors4Body_Tensors.cpp.

                                                           {
   if(0==_exampleDecay){
     _exampleDecay = new DecayTree(421);
     // remark: addDgtr always returns a pointer to the 
     // last daughter that was added, thus allowing these
     // chains:
     _exampleDecay->addDgtr( 225)->addDgtr( 221, -211);
     _exampleDecay->addDgtr( 9010221)->addDgtr(-321,  321);
   }
   return *_exampleDecay;
 }

◆ getVal()

double SF_DtoTS_TtoP0P1_StoP2P3::getVal ( IDalitzEvent & evt )

virtual

Implements SpinFactor.

Definition at line 821 of file SpinFactors4Body_Tensors.cpp.

                                                         {
   bool dbThis=false;
   if(! ( fsPS[0] && fsPS[1] && fsPS[2] && fsPS[3]) ) parseTree(evt.eventPattern());
   
   TLorentzVector pT = p(0, evt) + p(1, evt);
   TLorentzVector qT = p(0, evt) - p(1, evt);
   TLorentzVector pS = p(2, evt) + p(3, evt);
   TLorentzVector qD = pT - pS;
   TLorentzVector pD = pT + pS;
   
   double MT = mRes(T, evt);
   
   ZTspin2 tT(qT, pT, MT);
   
   const double units = GeV*GeV * GeV*GeV;
 
   if (_useZemachTensors) {
             ZTspin2 tD(qD,pD,pD.M());
             SpinSumT P2T(pT,MT);        
             SpinSumV P1T(pT,MT);        
             ZTspin1 t1T(qT, pT, MT);
 
             SymmLorentzMatrix qq(qD);
             //return tT.Contract_2(qq);
             if(dbThis){
                 cout << "" << endl;
                 cout << P2T.Sandwich(qD,qD,qT,qT)/units << endl;
                 cout << tT.Contract_2(qq)/units << endl;
                 cout << tT.Contract(qD).Dot(qD)/units << endl;
                 cout << (qD.Dot(t1T)*qD.Dot(t1T) + 1./3. * t1T.Dot(t1T) * P1T.Sandwich(qD,qD))/units << endl;
                 cout << "" << endl;
             }
             //return P2T.Sandwich(qD,qD,qT,qT)/units;
             return tT.Contract_2(tD)/units;
   }  
     
   double returnVal = (tT.Contract(qD)).Dot(qD) / units;
 
   if(dbThis){
       
       
     cout << " SF_DtoTS_TtoP0P1_StoP2P3::getVal "
      << " returning " << returnVal
      << endl;
   }
   return returnVal;
   
 }

◆ name()

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

inlinevirtual

Implements SpinFactor.

Definition at line 278 of file SpinFactors4Body_Tensors.h.

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

◆ parseTree()

bool SF_DtoTS_TtoP0P1_StoP2P3::parseTree ( const DalitzEventPattern & pat )

protectedvirtual

Definition at line 771 of file SpinFactors4Body_Tensors.cpp.

                                                                      {
   //  bool debugThis=false;
   if(fsPS.size() < 4) fsPS.reserve(4);
   if(theDecay(pat).nDgtr() != 2){
     cout << "ERROR in SF_DtoTS_TtoP0P1_StoP2P3::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() == "T" && ! dgtr->isFinalState()) T = dgtr;
     else if(dgtr->getVal().SVPAT() == "S" && ! dgtr->isFinalState()) S = dgtr;
   }
 
   if(0==S || 0==T){
     cout << "ERROR in SF_DtoTS_TtoP0P1_StoP2P3::parseTree"
      << " Didn't find T or S " << T.get() << ", " << S.get() << endl;
     return false;
   }
   if(T->nDgtr() != 2){
     cout << "ERROR in SF_DtoTS_TtoP0P1_StoP2P3::parseTree"
      << " T should have 2 daughters, but it says it has "
      << T->nDgtr() << "."
      << endl;
     return false;
   }
   fsPS[0] = T->getDgtrTreePtr(0);
   fsPS[1] = T->getDgtrTreePtr(1);
   // normalOrder(fsPS[0], fsPS[1]);
 
   if(S->nDgtr() != 2){
     cout << "ERROR in SF_DtoTS_TtoP0P1_StoP2P3::parseTree"
      << " S should have 2 daughters, but it says it has "
      << S->nDgtr() << "."
      << endl;
     return false;
   }
   fsPS[2] = S->getDgtrTreePtr(0);
   fsPS[3] = S->getDgtrTreePtr(1);
   // normalOrder(fsPS[2], fsPS[3]);
 
   // this->printYourself();
   return true;
 }

◆ printYourself()

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

virtual

Implements SpinFactor.

Definition at line 870 of file SpinFactors4Body_Tensors.cpp.

                                                              {
   //  bool debugThis = false;
   if(! ( fsPS[0] && fsPS[1] && fsPS[2] && fsPS[3]) ) return;
   os << "spin factor SF_DtoTS_TtoP0P1_StoP2P3"
      << "\n\t  ZTspin2 tT(qT, pT, MT);"
      << "\n\t  return: (tT.Contract(qD)).Dot(qD) / GeV^4"
      << "\n\t with pT = p(0) + p(1), qT = p(0) - p(1), pS = p(2) + p(3), qD = pT - qV;"
      << "\n\t    parsed tree " << theDecay().oneLiner()
      << "\n      like this:" << endl;
   this->printParsing(os);
 }