#include <SpinFactors4Body_Tensors.h>

Inheritance diagram for SF_DtoT1P0_T1toT2P1_T2toP2P3:

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

	SF_DtoT1P0_T1toT2P1_T2toP2P3 (const AssociatedDecayTree &theDecay)

virtual	~SF_DtoT1P0_T1toT2P1_T2toP2P3 ()

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

MINT::const_counted_ptr< AssociatedDecayTree >	T2

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

Constructor & Destructor Documentation

◆ SF_DtoT1P0_T1toT2P1_T2toP2P3()

SF_DtoT1P0_T1toT2P1_T2toP2P3::SF_DtoT1P0_T1toT2P1_T2toP2P3 ( const AssociatedDecayTree & theDecay )

inline

Definition at line 238 of file SpinFactors4Body_Tensors.h.

      : SpinFactor(theDecay, 4)
      , T1(0), T2(0)
      {
      }

◆ ~SF_DtoT1P0_T1toT2P1_T2toP2P3()

virtual SF_DtoT1P0_T1toT2P1_T2toP2P3::~SF_DtoT1P0_T1toT2P1_T2toP2P3 ( )

inlinevirtual

Definition at line 244 of file SpinFactors4Body_Tensors.h.

244 {}

Member Function Documentation

◆ exampleDecay()

const DecayTree & SF_DtoT1P0_T1toT2P1_T2toP2P3::exampleDecay ( )

virtual

Implements SpinFactor.

Definition at line 115 of file SpinFactors4Body_Tensors.cpp.

                                                            {
   return getExampleDecay();
 }

◆ getExampleDecay()

const DecayTree & SF_DtoT1P0_T1toT2P1_T2toP2P3::getExampleDecay ( )

static

Definition at line 118 of file SpinFactors4Body_Tensors.cpp.

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

◆ getVal()

double SF_DtoT1P0_T1toT2P1_T2toP2P3::getVal ( IDalitzEvent & evt )

virtual

Implements SpinFactor.

Definition at line 707 of file SpinFactors4Body_Tensors.cpp.

                                                             {
   bool dbThis=false;
   if(! ( fsPS[0] && fsPS[1] && fsPS[2] && fsPS[3]) ) parseTree(evt.eventPattern());
   
   TLorentzVector pT2 = p(2, evt) + p(3, evt);
   TLorentzVector qT2 = p(2, evt) - p(3, evt);
   TLorentzVector pT1 = pT2 + p(1, evt);
   TLorentzVector qT1 = pT2 - p(1, evt);
   TLorentzVector pD = pT1 + p(0, evt);
   TLorentzVector qD = pT1 - p(0, evt);
 
   double MT1 = mRes(T1 ,evt);
   double MT2 = mRes(T2, evt);
   
   ZTspin2 tT1(qT1, pT1, MT1);
   ZTspin2 tT2(qT2, pT2, MT2);
   
   TLorentzVector DT(tT1.Contract(qT1));
 
   const double units = GeV*GeV * GeV*GeV * GeV*GeV;
     
   if (_useZemachTensors) {
       ZTspin1 LT1(qT1, pT1, MT1);
       ZTspin1 LT2(qT2, pT2, MT2);
       SpinSumT PT1(pT1,MT1);
       
       TLorentzVector tmp1 = LeviCivita(LT1,pT1,LT2);
       TLorentzVector tmp2 = LeviCivita(LT1,pT1,pT2);
       tmp2 *= LT2.Dot(LT2)*1./(3.*pT2.Dot(pT2));
             
       return (PT1.Sandwich(qD,qD,tmp1,LT2)+PT1.Sandwich(qD,qD,tmp2,pT2))/units;
   }  
 
   //returns always 0 since tT2 is a symmetric matrix ! 
   //should not be used !  
   double returnVal = LeviCivita(pD, DT, tT2) / units;   
 
   if(dbThis){
     cout << " SF_DtoT1P0_T1toT2P1_T2toP2P3::getVal "
      << " returning " << returnVal
      << endl;
   }
   return returnVal;
 }

◆ name()

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

inlinevirtual

Implements SpinFactor.

Definition at line 249 of file SpinFactors4Body_Tensors.h.

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

◆ parseTree()

bool SF_DtoT1P0_T1toT2P1_T2toP2P3::parseTree ( const DalitzEventPattern & pat )

protectedvirtual

Definition at line 654 of file SpinFactors4Body_Tensors.cpp.

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

◆ printYourself()

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

virtual

Implements SpinFactor.

Definition at line 752 of file SpinFactors4Body_Tensors.cpp.

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