#include <SpinFactor3.h>

Inheritance diagram for SpinFactor3:

Public Member Functions
	SpinFactor3 (const AssociatedDecayTree &theDecay)

virtual double	getVal (IDalitzEvent &evt)

virtual const DecayTree &	exampleDecay ()

virtual std::string	name () const

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

virtual	~SpinFactor3 ()

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
bool	parseTree (const DalitzEventPattern &pat)

double	nonResVal () const

double	spinZeroVal () const

double	GSSpinFactor (IDalitzEvent &evt)

double	spinOneVal (IDalitzEvent &evt)

double	spinTwoVal (IDalitzEvent &evt)

double	spinOneFromMasses (IDalitzEvent &evt)

double	spinTwoFromMasses (IDalitzEvent &evt)

double	spinOneFromZemach (IDalitzEvent &evt)

bool	setSpin ()

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

int	_spin

bool	_nonResonant

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 10 of file SpinFactor3.h.

Constructor & Destructor Documentation

◆ SpinFactor3()

SpinFactor3::SpinFactor3 ( const AssociatedDecayTree & theDecay )

inline

Definition at line 32 of file SpinFactor3.h.

     : SpinFactor(theDecay, 3)
     , R(0)
     , _spin(-1)
     , _nonResonant(false)
     {
       // we'll do the parsing now when we see the 1st event.
       //if( ! parseTree()){
       //throw "error in parseTree in constructor of SpinFactor3";
       //}
     }

◆ ~SpinFactor3()

virtual SpinFactor3::~SpinFactor3 ( )

inlinevirtual

Definition at line 51 of file SpinFactor3.h.

51 {}

Member Function Documentation

◆ exampleDecay()

const DecayTree & SpinFactor3::exampleDecay ( )

virtual

Implements SpinFactor.

Definition at line 31 of file SpinFactor3.cpp.

                                           {
   return getExampleDecay();
 }

◆ getExampleDecay()

const DecayTree & SpinFactor3::getExampleDecay ( )

static

Definition at line 22 of file SpinFactor3.cpp.

                                              {
   if(0==_exampleDecay){
     _exampleDecay = new DecayTree(421);
     _exampleDecay->addDgtr(310, 113)->addDgtr(211, -211);// Ks rho - why not
     _exampleDecay->addDgtr(-211, 200321)->addDgtr(321, -10323);
   }
   return *_exampleDecay;
 }

◆ getVal()

double SpinFactor3::getVal ( IDalitzEvent & evt )

virtual

Implements SpinFactor.

Reimplemented in GSSpinFactor.

Definition at line 106 of file SpinFactor3.cpp.

                                            {
   if(_nonResonant){
     return nonResVal();
   }
   
   if(! ( fsPS[0] && fsPS[1] && fsPS[2])) parseTree(evt.eventPattern());
   
   if(0 == R){
     cout << "ERROR in SpinFactor3::getVal(): 0 == R"
      << endl;
     throw "such things shouldn't happen";
   }
 
   if(_spin == 0){
     return spinZeroVal();
   }else if(_spin == 1){
     return spinOneVal(evt);
   }else if(_spin == 2){
     return spinTwoVal(evt);
   }else{    
     cout << "SpinFactor3::getVal() Don't know how to"
      << " handle resonance with spin = " << R->getVal().J()
      << "\n    > looking at decay\n" << theDecay()
      << "\n    > Sorry, will crash now!" << endl;
     throw "sorry";
   }
   return -9999;
   
 }

◆ GSSpinFactor()

double SpinFactor3::GSSpinFactor ( IDalitzEvent & evt )

protected

Definition at line 136 of file SpinFactor3.cpp.

                                                  {
 
   double m2AC = (p(0, evt) + p(1, evt)).M2();
   double m2BC = (p(0, evt) + p(2, evt)).M2();
 
   //     cout << "GS spin factor " << m2AC - m2BC << endl;
   return (m2AC - m2BC)/(GeV*GeV);
 }

◆ name()

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

inlinevirtual

Implements SpinFactor.

Reimplemented in GSSpinFactor.

Definition at line 47 of file SpinFactor3.h.

                                 {
     return "SpinFactor3(" + theDecay().oneLiner() + ")";
   }

◆ nonResVal()

double SpinFactor3::nonResVal ( ) const

inlineprotected

Definition at line 18 of file SpinFactor3.h.

18 {return 1;}

◆ parseTree()

bool SpinFactor3::parseTree ( const DalitzEventPattern & pat )

protected

Definition at line 58 of file SpinFactor3.cpp.

                                                         {
   if(fsPS.size() < 3) fsPS.reserve(3);
 
   if(theDecay(pat).nDgtr() == (int)(theDecay(pat).finalState().size())){
     _nonResonant = true;
     return true;
   }
 
   for(int i=0; i< theDecay(pat).nDgtr(); i++){
     const_counted_ptr<AssociatedDecayTree> dgtr= theDecay(pat).getDgtrTreePtr(i);
     if   (! dgtr->isFinalState())      R = dgtr;
     else                         fsPS[0] = dgtr;
   }
   if(0==R || 0==fsPS[0]){
     cout << "SpinFactor3::parseTree"
      << " Didn't find R or P1 " << R << ", " << fsPS[0] << endl;
     return false;
   }
   if(R->nDgtr() != 2){
     cout << "ERROR in SpinFactor3::parseTree"
      << " Resonance should have 2 daughters, but it says it has "
      << R->nDgtr() << "."
      << endl;
     return false;
   }
   fsPS[1] = R->getDgtrTreePtr(0);
   fsPS[2] = R->getDgtrTreePtr(1);
 
   //normalOrder(fsPS[1], fsPS[2]);
 
 //  printYourself();
 
   setSpin();
   return true;
 
 }

◆ printYourself()

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

virtual

Implements SpinFactor.

Definition at line 95 of file SpinFactor3.cpp.

                                                   {
   os << "INFO from SpinFactor3::printYourself():\n"
      << "    > parsed the following tree "
      << theDecay().oneLiner() << " like this:\n"
      << "    > R =  " << R->getVal().name()
      << ", fsPS[0]=C= " << fsPS[0]->getVal().name()
      << ", fsPS[1]=A= " << fsPS[1]->getVal().name()
      << ", fsPS[2]=B= " << fsPS[2]->getVal().name()
      << endl;
 }

◆ setSpin()

bool SpinFactor3::setSpin ( )

protected

Definition at line 35 of file SpinFactor3.cpp.

                          {
   if(0 == R){
     cout << " SpinFactor3::setSpin can't set spin, R==0"
      << " probably need to parse tree, first."
      << endl;
     return false;
   }
   if(R->getVal().J() == "0"){
     _spin = 0;
   }else if(R->getVal().J() == "1"){
     _spin = 1;
   }else if(R->getVal().J() == "2"){
     _spin = 2;
   }else{    
     cout << "SpinFactor3::setSpin() Don't know how to"
      << " handle resonance with spin = " << R->getVal().J()
      << "\n    > looking at decay\n" << theDecay()
      << "\n    > Sorry, will crash now!" << endl;
     throw "sorry";
   }
   return true;
 }

◆ spinOneFromMasses()

double SpinFactor3::spinOneFromMasses ( IDalitzEvent & evt )

protected

Definition at line 209 of file SpinFactor3.cpp.

                                                       {
   // parsed as:
   // D -> V P0; V->P1, P2;
   // fsPS[0] && fsPS[1] && fsPS[2]
   
   double MV = mRes(R, evt);
   double mA = fsPS[1]->getVal().mass();
   double mB = fsPS[2]->getVal().mass();
 
   double mC = fsPS[0]->getVal().mass();
 
   double mD = (p(0, evt)+p(1, evt)+p(2, evt)).M();
 
   double m2AC = (p(0, evt) + p(1, evt)).M2();
   double m2BC = (p(0, evt) + p(2, evt)).M2();
 
   return (m2AC - m2BC + (mD*mD - mC*mC)*(mB*mB-mA*mA)/(MV*MV))/(GeV*GeV);
 
 }

◆ spinOneFromZemach()

double SpinFactor3::spinOneFromZemach ( IDalitzEvent & evt )

protected

Definition at line 181 of file SpinFactor3.cpp.

                                                       {
   // all wrong - just experimenting - don't use!!
 
   TLorentzVector pV = p(1, evt) + p(2, evt);
   TLorentzVector qV = p(1, evt) - p(2, evt);
   TLorentzVector pD = pV + p(0, evt);
   TLorentzVector qD = pV - p(0, evt);
   //TLorentzVector qD = pD + p(0);
 
   double mr = mRes(R, evt);
   double mD = pD.M();
 
   ZTspin1 tV(qV, pV, mr);
   ZTspin1 tD(qD, pD, mD);
   //  ZTspin1 tP(pD + p(0), pD-p(0), mD);
   ZTspin1 tP(p(0, evt), pD, mD);
   //double norm = 1.;//tV.M() * tD.M();
   //return tV.Contract(pD + p(0)); << this works
   return tV.Contract(tP );
   //  return tV.Contract(tD)/norm;
 
   //ZTspin1 tVinD(tV, pD, mD);
 
   //double norm = 1.;// tVinD.M() * tD.M();
   //return tD.Contract(tVinD)/norm;
 
 }

◆ spinOneVal()

double SpinFactor3::spinOneVal ( IDalitzEvent & evt )

protected

Definition at line 145 of file SpinFactor3.cpp.

                                                {
   // parsed as:
   // D -> V P0; V->P1, P2;
   // fsPS[0] && fsPS[1] && fsPS[2]
 
   bool dbThis=false;
 
   TLorentzVector pV = p(1, evt) + p(2, evt);
   TLorentzVector pD = pV + p(0, evt);
   double mr = mRes(R, evt);
   SpinSumV spin_sum(pV, mr);
 
   TLorentzVector lhs = pD + p(0, evt);
   TLorentzVector rhs = p(1, evt) - p(2, evt);
 
   if(dbThis){
     cout << " spinOneVal for " << theDecay().oneLiner()
      << "\n   > compare: Sandwich: "  
      << -spin_sum.Sandwich(lhs, rhs)
      << "\n   > from masses             "  
      << spinOneFromMasses(evt) 
      << "\n   > ratio sw/m         "  
      << -spin_sum.Sandwich(lhs, rhs)/spinOneFromMasses(evt) 
      << "\n   > from Zemach " << spinOneFromZemach(evt)
      << "\n   > ratios sw/Zemach "
      << spin_sum.Sandwich(lhs, rhs)/spinOneFromZemach(evt)
      << endl;
   }
 
   // note '-' sign. This makes my decay-tree
   // sorting algorithm and this spin factor
   // comparible with BaBar/BELLE spin
   // factor conventions.
   return -spin_sum.Sandwich(lhs, rhs)/(GeV*GeV);
 }

◆ spinTwoFromMasses()

double SpinFactor3::spinTwoFromMasses ( IDalitzEvent & evt )

protected

Definition at line 229 of file SpinFactor3.cpp.

                                                       {
   // parsed as:
   // D -> V P0; V->P1, P2;
   // fsPS[0] && fsPS[1] && fsPS[2]
   
   double MV = mRes(R, evt);
   double mA = fsPS[1]->getVal().mass();
   double mB = fsPS[2]->getVal().mass();
 
   double mC = fsPS[0]->getVal().mass();
 
   double mD = (p(0, evt)+p(1, evt)+p(2, evt)).M();
 
   double m2AB = (p(1, evt) + p(2, evt)).M2();
   double m2AC = (p(0, evt) + p(1, evt)).M2();
   double m2BC = (p(0, evt) + p(2, evt)).M2();
   
   Double_t term1 = m2BC-m2AC+(mD*mD-mC*mC)*(mA*mA-mB*mB)/(MV*MV);
   Double_t term2 = m2AB-2.0*mD*mD-2.0*mC*mC+pow(mD*mD-mC*mC,2)/(MV*MV);
   Double_t term3 = m2AB-2.0*mA*mA-2.0*mB*mB+pow(mA*mA-mB*mB,2)/(MV*MV);
   
   return (pow(term1,2)-(1.0/3.0)*term2*term3)/(GeV*GeV*GeV*GeV);
 }

◆ spinTwoVal()

double SpinFactor3::spinTwoVal ( IDalitzEvent & evt )

protected

Definition at line 253 of file SpinFactor3.cpp.

                                                {
   // parsed as:
   // D -> T P0; T->P1, P2;
   bool dbThis=false;
 
   TLorentzVector pT = p(1, evt) + p(2, evt);
   TLorentzVector pD = pT + p(0, evt);
   double mr = mRes(R, evt);
   SpinSumT spin_sum(pT, mr);
 
   TLorentzVector lhs = pD + p(0, evt);
   TLorentzVector rhs = p(1, evt) - p(2, evt);
 
   double returnVal = spin_sum.Sandwich(lhs, lhs, rhs, rhs)/(GeV*GeV*GeV*GeV);
   if(dbThis){
     cout << "spin-2 spin factor got called for "
      << theDecay().oneLiner()
      << "\n   >  returning " << returnVal
      << "\n   >  compare " << spinTwoFromMasses(evt)
      << "\n   >  ratio " << returnVal/spinTwoFromMasses(evt)
      << endl;
   }
 
   return returnVal;
 }

◆ spinZeroVal()

double SpinFactor3::spinZeroVal ( ) const

inlineprotected

Definition at line 19 of file SpinFactor3.h.

19 { return 1;}

Member Data Documentation

◆ _exampleDecay

DecayTree * SpinFactor3::_exampleDecay =0

staticprotected

Definition at line 12 of file SpinFactor3.h.

◆ _nonResonant

bool SpinFactor3::_nonResonant

protected

Definition at line 15 of file SpinFactor3.h.

◆ _spin

int SpinFactor3::_spin

protected

Definition at line 14 of file SpinFactor3.h.

◆ R

MINT::const_counted_ptr<AssociatedDecayTree> SpinFactor3::R

protected

Definition at line 13 of file SpinFactor3.h.

The documentation for this class was generated from the following files:

Mint/SpinFactor3.h
src/Mojito/SpinFactors/SpinFactor3.cpp

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ SpinFactor3()

◆ ~SpinFactor3()

Member Function Documentation

◆ exampleDecay()

◆ getExampleDecay()

◆ getVal()

◆ GSSpinFactor()

◆ name()

◆ nonResVal()

◆ parseTree()

◆ printYourself()

◆ setSpin()

◆ spinOneFromMasses()

◆ spinOneFromZemach()

◆ spinOneVal()

◆ spinTwoFromMasses()

◆ spinTwoVal()

◆ spinZeroVal()

Member Data Documentation

◆ _exampleDecay

◆ _nonResonant

◆ _spin

◆ R