#include <SpinFactors4Body_Tensors.h>

Inheritance diagram for SF_DtoAP0_AtoTP1_TtoP2P3:

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

	SF_DtoAP0_AtoTP1_TtoP2P3 (const AssociatedDecayTree &theDecay)

virtual	~SF_DtoAP0_AtoTP1_TtoP2P3 ()

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

MINT::const_counted_ptr< AssociatedDecayTree >	T

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

Constructor & Destructor Documentation

◆ SF_DtoAP0_AtoTP1_TtoP2P3()

SF_DtoAP0_AtoTP1_TtoP2P3::SF_DtoAP0_AtoTP1_TtoP2P3 ( const AssociatedDecayTree & theDecay )

inline

Definition at line 384 of file SpinFactors4Body_Tensors.h.

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

◆ ~SF_DtoAP0_AtoTP1_TtoP2P3()

virtual SF_DtoAP0_AtoTP1_TtoP2P3::~SF_DtoAP0_AtoTP1_TtoP2P3 ( )

inlinevirtual

Definition at line 390 of file SpinFactors4Body_Tensors.h.

390 {}

Member Function Documentation

◆ exampleDecay()

const DecayTree & SF_DtoAP0_AtoTP1_TtoP2P3::exampleDecay ( )

virtual

Implements SpinFactor.

Definition at line 179 of file SpinFactors4Body_Tensors.cpp.

                                                        {
   return getExampleDecay();
 }

◆ getExampleDecay()

const DecayTree & SF_DtoAP0_AtoTP1_TtoP2P3::getExampleDecay ( )

static

Definition at line 172 of file SpinFactors4Body_Tensors.cpp.

                                                           {
   _exampleDecay = new DecayTree(421);
   // example: D->a1(1260) pi, a1 -> f2(1270) pi 
   // (kinematically challenged, but in PDG - broad particles...)
   _exampleDecay->addDgtr(-211, 20213)->addDgtr(211, 225)->addDgtr(211, -211);
   return *_exampleDecay;
 }

◆ getVal()

double SF_DtoAP0_AtoTP1_TtoP2P3::getVal ( IDalitzEvent & evt )

virtual

Implements SpinFactor.

Definition at line 1266 of file SpinFactors4Body_Tensors.cpp.

                                                         {
   bool dbThis=false;
   if(! ( fsPS[0] && fsPS[1] && fsPS[2] && fsPS[3]) ) parseTree(evt.eventPattern());
 
   TLorentzVector pT = p(2, evt) + p(3, evt);
   TLorentzVector qT = p(2, evt) - p(3, evt);
   TLorentzVector pA = pT + p(1, evt);
   TLorentzVector qA = pT - p(1, evt);
   TLorentzVector pD = pA + p(0, evt);
   TLorentzVector qD = pA - p(0, evt);
 
   double MA = mRes(A, evt);
   double MT = mRes(T, evt);
   //double MV = mRes(S);
   
   ZTspin1 tA(qA, pA, MA);
   ZTspin2 tT(qT, pT, MT);
 
   const double units = GeV*GeV*GeV*GeV;
     
   if (_useZemachTensors) {
       ZTspin1 LD(qD,pD,pD.M());
       SpinSumT PT(pT,MT);
       //return (tT.Contract(tA)).Dot(LD);
       return PT.Sandwich(qT,qT,LD,tA)/units;
   }
 
   double returnVal =  (tT.Contract(tA)).Dot(qD) / units;
 
   if(dbThis){
     cout << " SF_DtoAP0_AtoTP1_TtoP2P3::getVal "
      << " returning " << returnVal
      << endl;
   }
   return returnVal;
 }

◆ name()

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

inlinevirtual

Implements SpinFactor.

Definition at line 395 of file SpinFactors4Body_Tensors.h.

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

◆ parseTree()

bool SF_DtoAP0_AtoTP1_TtoP2P3::parseTree ( const DalitzEventPattern & pat )

protectedvirtual

Definition at line 1213 of file SpinFactors4Body_Tensors.cpp.

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

◆ printYourself()

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

virtual

Implements SpinFactor.

Definition at line 1303 of file SpinFactors4Body_Tensors.cpp.

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