#include <SpinFactors4Body_Tensors.h>

Inheritance diagram for SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3:

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

	SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3 (const AssociatedDecayTree &theDecay)

virtual	~SF_DtoPseudoTP0_PseudoTtoTP1_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 >	PT

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

Constructor & Destructor Documentation

◆ SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3()

SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3::SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3 ( const AssociatedDecayTree & theDecay )

inline

Definition at line 353 of file SpinFactors4Body_Tensors.h.

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

◆ ~SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3()

virtual SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3::~SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3 ( )

inlinevirtual

Definition at line 359 of file SpinFactors4Body_Tensors.h.

359 {}

Member Function Documentation

◆ exampleDecay()

const DecayTree & SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3::exampleDecay ( )

virtual

Implements SpinFactor.

Definition at line 150 of file SpinFactors4Body_Tensors.cpp.

                                                                    {
   return getExampleDecay();
 }

◆ getExampleDecay()

const DecayTree & SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3::getExampleDecay ( )

static

Definition at line 145 of file SpinFactors4Body_Tensors.cpp.

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

◆ getVal()

double SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3::getVal ( IDalitzEvent & evt )

virtual

Implements SpinFactor.

Definition at line 1152 of file SpinFactors4Body_Tensors.cpp.

                                                                     {
   bool dbThis=false;
   if(! ( fsPS[0] && fsPS[1] && fsPS[2] && fsPS[3]) ) parseTree(evt.eventPattern());
 
   double invGeV = 1;// 1.0/GeV;
   
   TLorentzVector pT = p(2, evt) + p(3, evt);
   TLorentzVector qT = p(2, evt) - p(3, evt);
   TLorentzVector pPT = pT + p(1, evt);
   TLorentzVector qPT = pT - p(1, evt);
   TLorentzVector pD = pPT + p(0, evt);
   TLorentzVector qD = pPT - p(0, evt);
 
   double MPT = mRes(PT, evt);
   double MT  = mRes(T , evt);
   //double MV = mRes(S);
   
   ZTspin2 tPT(qD*invGeV, pPT*invGeV, MPT*invGeV);
   ZTspin2 tT(qT*invGeV, pT*invGeV, MT*invGeV);
 
   const double units = GeV*GeV*GeV*GeV;
     
   if (_useZemachTensors) {
             SpinSumT PT1(pPT,MPT);
             ZTspin1 LT(qT,pT,MT);
       
             return (PT1.Sandwich(qD,qD,LT,LT) + 1./3. * LT.Dot(LT)/(pT.Dot(pT)) * PT1.Sandwich(qD,qD,pT,pT))/units;
   }
 
   double returnVal =  tPT.Contract_2(tT) / units;
 
   if(dbThis){
     cout << " SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3::getVal "
      << " returning " << returnVal
      << " = " << tPT.Contract_2(tT) << " / " << units
      << "\n MPT = " << MPT << ", MT = " << MT
      << endl;
   }
   return returnVal;
 }

◆ name()

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

inlinevirtual

Implements SpinFactor.

Definition at line 364 of file SpinFactors4Body_Tensors.h.

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

◆ parseTree()

bool SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3::parseTree ( const DalitzEventPattern & pat )

protectedvirtual

Definition at line 1099 of file SpinFactors4Body_Tensors.cpp.

                                                                                  {
   //  bool debugThis=false;
   if(fsPS.size() < 4) fsPS.reserve(4);
   if(theDecay(pat).nDgtr() != 2){
     cout << "ERROR in SF_DtoPseudoTP0_PseudoTtoTP1_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() == "PT" && ! dgtr->isFinalState()) PT = dgtr;
   }
   if(0==PT || 0==fsPS[0]){
     cout << "ERROR in SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3_BASE::parseTree"
      << " Didn't find T or P0 " << PT.get() << ", " << fsPS[0].get() << endl;
     return false;
   }
   
   if(PT->nDgtr() != 2){
     cout << "ERROR in SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3_BASE::parseTree"
      << " T should have 2 daughters, but it says it has "
      << PT->nDgtr() << "."
      << endl;
     return false;
   }
   for(int i=0; i< PT->nDgtr(); i++){
     const_counted_ptr<AssociatedDecayTree> dgtr= PT->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_DtoPseudoTP0_PseudoTtoTP1_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_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3_BASE::parseTree"
      << " V 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_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3::printYourself ( std::ostream & os = std::cout ) const

virtual

Implements SpinFactor.

Definition at line 1193 of file SpinFactors4Body_Tensors.cpp.

                                                                          {
   //  bool debugThis = false;
   if(! ( fsPS[0] && fsPS[1] && fsPS[2] && fsPS[3]) ) return;
   os << "spin factor SF_DtoPseudoTP0_PseudoTtoTP1_TtoP2P3"
      << "\n\t  ZTspin2 tPT(qD, pPT, MPT);"
      << "\n\t  ZTspin2 tT(qT, pT, MT);"
      << "\n\t  pT = p(2) + p(3); qT = p(2) - p(3);"
      << "\n\t  pPT = pT + p(1); qPT = pT - p(1);"
      << "\n\t  pD = pPT + p(0); qD = pPT - p(0);"
      << "\n\t  tPT.Contract_2(tT)/ GeV^8"
      << "\n\t    parsed tree " << theDecay().oneLiner()
      << "\n      like this:" << endl;
   this->printParsing(os);
 }