#include <DalitzEvent.h>

Inheritance diagram for DalitzEvent:

Public Member Functions
	DalitzEvent ()

	DalitzEvent (const DalitzEventPattern &pat, const std::vector< TLorentzVector > &mumAndDgtr_p4)

	DalitzEvent (const DalitzEventPattern &pat, double t01_in, double s12_in, double s23_in, double s34_in, double t40_in)
	Constructor for a 4-body decay at a given point in phase space. More...

	DalitzEvent (const DalitzEventPattern &pat, const double s13, const double s23)
	Constructor for a 3-body decay at a given point in phase space. More...

	DalitzEvent (const DalitzEventPattern &pat, const std::vector< TVector3 > &mumAndDgtr_p3)

	DalitzEvent (const DalitzEventPattern &pat)

	DalitzEvent (const DalitzEventPattern &pat, TRandom *rnd)

	DalitzEvent (TNtupleD *ntp)

	DalitzEvent (const IDalitzEvent &ievt)

	DalitzEvent (const IDalitzEvent *ievt)

	DalitzEvent (const DalitzEvent *other)

	DalitzEvent (const DalitzEvent &other)

virtual	~DalitzEvent ()

IDalitzEvent *	clone () const

virtual double	getWeight () const

virtual void	setWeight (double w)

virtual void	setGeneratorPdfRelativeToPhaseSpace (double gpdf)

virtual double	getGeneratorPdfRelativeToPhaseSpace () const

virtual const std::vector< double > &	getVectorOfValues () const

virtual std::vector< double > &	getVectorOfValues ()

virtual const std::vector< double > &	getVectorOfWeights () const

virtual std::vector< double > &	getVectorOfWeights ()

virtual void	setValueInVector (unsigned int i, double value)

virtual void	setWeightInVector (unsigned int i, double weight)

virtual double	getValueFromVector (unsigned int i) const

virtual double	getWeightFromVector (unsigned int i) const

virtual void	P_conjugateYourself ()

virtual void	C_conjugateYourself ()

virtual void	CP_conjugateYourself ()

bool	resetST ()

const Permutation &	currentPermutation () const

void	setPermutationIndex (int i)

int	numPermutations () const

int	permutationIndex () const

bool	fromTree (TTree *tree)

bool	fromNtuple (TNtupleD *ntp)

bool	fromParasTree (TTree *ntp)

bool	fromParasTreeOld (TTree *ntp)

virtual void	setAValue (double a)

virtual double	getAValue () const

virtual const DalitzEventPattern &	eventPattern () const

virtual const TLorentzVector &	p (unsigned int i) const

virtual TLorentzVector &	p (unsigned int i)

virtual double	s (unsigned int i, unsigned int j) const

virtual double	t (unsigned int i, unsigned int j) const

virtual double	sij (const std::vector< int > &indices) const

virtual double	sij (const MINT::PolymorphVector< int > &indices) const

virtual void	setMothers3Momentum (const TVector3 &mp3)

virtual double	phaseSpace () const

virtual bool	retrieveValue (int i, std::complex< double > &value, long int configNumber=1)

virtual void	setValue (int i, const std::complex< double > &value, long int configNumber=1)

virtual bool	retrieveValue (int i, double value, long int configNumber=1)

virtual void	setValue (int i, double value, long int configNumber=1)

double	sijMin (const std::vector< int > &indices) const

double	sijMin (const MINT::PolymorphVector< int > &indices) const

double	sijMax (const std::vector< int > &indices) const

double	sijMax (const MINT::PolymorphVector< int > &indices) const

int	nDgtr () const

double	m (unsigned int i) const

double	m2 (unsigned int i) const

TLorentzVector	p_dgtr_sum () const

bool	kinematicallyAllowed (double epsilon=1.e-9) const

void	print (std::ostream &os=std::cout) const

void	generateThisToPhaseSpace (TRandom *rnd=0)

void	generateThisToPhaseSpace (double p_of_D, TRandom *rnd=0)

void	generateThisToPhaseSpace (const TVector3 &pVec_of_D, TRandom *rnd=0)

std::string	makeNtupleVarnames () const

bool	fillNtupleVarArray (std::vector< Double_t > &array) const

unsigned int	ntupleVarArraySize () const

Calculate4BodyProps	makeCalculate4BodyProps () const

Public Member Functions inherited from IDalitzEvent
virtual	~IDalitzEvent ()

Static Public Member Functions
static long int	assignUniqueRememberNumber ()

static long int	eventCounter ()

Protected Member Functions
const TLorentzVector &	p_intern (unsigned int i) const

TLorentzVector &	p_intern (unsigned int i)

double	s_intern (unsigned int i, unsigned int j) const

double	t_intern (unsigned int i, unsigned int j) const

double	evalsij_intern (const std::vector< int > &intern_indices) const

double	sijMap_intern (const std::vector< int > &intern_indices) const

double	sijMap (const std::vector< int > &indices) const

bool	shoutAndKill ()

double	BDet () const

double	Delta4 () const

double	phaseSpace3 (double epsilon=1.e-9) const

double	phaseSpace4 () const

bool	setMomenta (const std::vector< TLorentzVector > &mumAndDgtr_p4)

bool	setMomenta (const std::vector< TVector3 > &mumAndDgtr_p3)

void	setP (unsigned int i, const TLorentzVector &p4)

Protected Member Functions inherited from IDalitzEvent
	IDalitzEvent ()

Static Protected Member Functions
static std::string	prtToNtpName (const std::string &s_in)

static std::string	ntpToPrtName (const std::string &s_in)

static int	singleParticleNtpArraySize ()

static bool	parseNtpEntryName (const std::string &entry, std::string &part1, std::string &part2)

Protected Attributes
DalitzEventPattern	_pat

std::vector< TLorentzVector >	_p

double	_rememberPhaseSpace

RememberAnythingFast< std::complex< double > >	_rememberComplexFast

RememberAnythingFast< double >	_rememberDoubleFast

double	_aValue

double	_weight

double	_generatorPdfRelativeToPhaseSpace

std::vector< double >	_vectorOfValues

std::vector< double >	_vectorOfWeights

std::vector< std::vector< double > >	_s

std::vector< std::vector< double > >	_t

std::map< std::vector< int >, double >	_sijMap

int	_permutationIndex

Permutator	_perm

Static Protected Attributes
static const char	ntpNameChars [] = { '#', '~', '\0' }

static const char	prtNameChars [] = { '+', '-', '\0' }

static long int	_eventCounter =0

static long int	_rememberVectorCounter =0

Detailed Description

Definition at line 27 of file DalitzEvent.h.

Constructor & Destructor Documentation

◆ DalitzEvent() [1/12]

DalitzEvent::DalitzEvent ( )

Definition at line 33 of file DalitzEvent.cpp.

   : _rememberPhaseSpace(-9999.)
   , _aValue(-9999.)
   , _weight(1)
   , _generatorPdfRelativeToPhaseSpace(1)
   , _vectorOfValues()
   , _vectorOfWeights()
   , _permutationIndex(0)
   , _perm()
 {
   _eventCounter++;
 }

◆ DalitzEvent() [2/12]

DalitzEvent::DalitzEvent	(	const DalitzEventPattern &	pat,
		const std::vector< TLorentzVector > &	mumAndDgtr_p4
	)

Definition at line 45 of file DalitzEvent.cpp.

   : _pat(pat)
   , _rememberPhaseSpace(-9999.)
   , _aValue(-9999.)
   , _weight(1)
   , _generatorPdfRelativeToPhaseSpace(1)
   , _vectorOfValues()
   , _vectorOfWeights()
   , _s(pat.size())
   , _t(pat.size())
   , _permutationIndex(0)
   , _perm(pat)
 {
   setMomenta(mumAndDgtrs_p4);
   if(_pat.size() < 3 || _pat.size() != _p.size()) shoutAndKill();
   resetST();
   _eventCounter++;
 }

◆ DalitzEvent() [3/12]

DalitzEvent::DalitzEvent	(	const DalitzEventPattern &	pat,
		double	t01_in,
		double	s12_in,
		double	s23_in,
		double	s34_in,
		double	t40_in
	)

Constructor for a 4-body decay at a given point in phase space.

Definition at line 66 of file DalitzEvent.cpp.

   : _pat(pat)
   , _rememberPhaseSpace(-9999.)
   , _aValue(-9999.)
   , _weight(1)
   , _generatorPdfRelativeToPhaseSpace(1)
   , _vectorOfValues()
   , _vectorOfWeights()
   , _s(pat.size())
   , _t(pat.size())
   , _permutationIndex(0)
   , _perm(pat)
 {
   bool dbThis=false;
 
   if(_pat.size() != 5){
     cout << "DalitzEvent::DalitzEvent ERROR "
      << " event pattern " << pat
      << " doesn't fit 4-body constructor"
      << endl;
     throw "shit happens";
   }
   Calculate4BodyProps c4bp(t01_in, s12_in, s23_in, s34_in, t40_in
                , _pat[0].mass()
                , _pat[1].mass(), _pat[2].mass()
                , _pat[3].mass(), _pat[4].mass()
                );
   if(dbThis){
     cout << " c4bp phase space " << c4bp.phaseSpaceFactor() << endl;
   }
   std::vector<TLorentzVector> p4List(5);
   if(c4bp.phaseSpaceFactor() > 0){ 
     for(int i=0; i<=4; i++){
       p4List[i] = c4bp.pVec(i);
       if(dbThis)cout << "p4List[" << i << "] " << p4List[i] << endl;
     }
   }else{
     for(int i=0; i<=4; i++){
       TLorentzVector p4(-9999, -9999, -9999, -9999);
       p4List[i] = p4;
     }
   }
   setMomenta(p4List);
   if(_pat.size() < 3 || _pat.size() != _p.size()) shoutAndKill();
   resetST();
 
   if(dbThis){
     cout << "constructed event from this pattern: " << _pat << "\n and "
      << t01_in << ", " << s12_in << ", " << s23_in << ", " << s34_in << ", " << t40_in
      << endl;
     
     cout << " and got out: "
      << t_intern(0,1) << ", " << s_intern(1,2) << ", " << s_intern(2,3)
      << ", " << s_intern(3,4) << ", " << t_intern(4,0)
      << endl;
   }
   _eventCounter++;
 }

◆ DalitzEvent() [4/12]

DalitzEvent::DalitzEvent	(	const DalitzEventPattern &	pat,
		const double	s13,
		const double	s23
	)

Constructor for a 3-body decay at a given point in phase space.

Definition at line 129 of file DalitzEvent.cpp.

   : _pat(pat)
   , _rememberPhaseSpace(-9999.)
   , _aValue(-9999.)
   , _weight(1)
   , _generatorPdfRelativeToPhaseSpace(1)
   , _vectorOfValues()
   , _vectorOfWeights()
   , _s(pat.size())
   , _t(pat.size())
   , _permutationIndex(0)
   , _perm(pat)
 {
   bool dbThis=false;
 
   if(_pat.size() != 4){
     cout << "DalitzEvent::DalitzEvent ERROR "
      << " event pattern " << pat
      << " doesn't fit 4-body constructor"
      << endl;
     throw "shit happens";
   }
 
   double m0 = pat[0].mass() ;
   double m1 = pat[1].mass() ;
   double m2 = pat[2].mass() ;
   double m3 = pat[3].mass() ;
 
   // Calculate energies of pi+, pi- from known parameters
   double E1 = ( m0 * m0 + m1 * m1 - s23 ) / (2 * m0) ;
   double E2 = ( m0 * m0 + m2 * m2 - s13 ) / (2 * m0) ;
 
   // Get E3 from conservation of energy
   double E3 = m0 - E1 - E2 ;
 
   // Reject events where kinematics are impossible
   double pz1Sq = E1 * E1 - m1 * m1 ;
   if(pz1Sq < 0){
     cout << "Invalid phase space point for the given pattern (s13, s23) = ("
      << s13 << ", " << s23 << ")" << endl ;
     throw out_of_range("") ;
   }
 
   // Get pz1 from mass/energy relation (co-ordinate system defined so that px1 = py1 = 0)
   double pz1 = sqrt(pz1Sq) ;
 
   // Calculate pz3 from conservation of momentum 
   double pz3 = ( E2 * E2 - m2 * m2 - E3 * E3 + m3 * m3 - pz1 * pz1 ) / ( 2*pz1 ) ;
 
   // Use pz3 to get py3 (co-ordinate system designed so that px3 = 0)
   double py3Sq =  E3 * E3 - m3 * m3 - pz3 * pz3 ;
 
   // Reject events where kinematics are impossible
   if(py3Sq < 0){
     cout << "Invalid phase space point for the given pattern (s13, s23) = ("
      << s13 << ", " << s23 << ")" << endl ;
     throw out_of_range("") ;
   }
 
   double py3 = sqrt(py3Sq) ;
   // Conservation of momentum => py2 = -py3 and pz2 = -(pz1 + pz3)
   double py2 = -1 * py3 ;
   double pz2 = -1 * (pz1 + pz3) ;
   std::vector<TLorentzVector> p4List(4, TLorentzVector(0, 0, 0, m0));
   p4List[1] = TLorentzVector(0., 0., pz1, E1) ;
   p4List[2] = TLorentzVector(0., py2, pz2, E2) ;
   p4List[3] = TLorentzVector(0., py3, pz3, E3) ;
 
   setMomenta(p4List);
   if(_pat.size() < 3 || _pat.size() != _p.size()) shoutAndKill();
   resetST();
 
   if(!kinematicallyAllowed()){
     cout << "Event not kinematically allowed! Was given (s13, s23) = ("
      << s13 << ", " << s23 << "), got (" << s_intern(1, 3) << ", "
      << s_intern(2, 3) << ")" << endl ;
   throw out_of_range("") ;
   }
   _eventCounter++;
 }

◆ DalitzEvent() [5/12]

DalitzEvent::DalitzEvent	(	const DalitzEventPattern &	pat,
		const std::vector< TVector3 > &	mumAndDgtr_p3
	)

Definition at line 212 of file DalitzEvent.cpp.

   : _pat(pat)
   , _rememberPhaseSpace(-9999.)
   , _aValue(-9999.)
   , _weight(1)
   , _generatorPdfRelativeToPhaseSpace(1)
   , _vectorOfValues()
   , _vectorOfWeights()
   , _s(pat.size())
   , _t(pat.size())
   , _permutationIndex(0)
   , _perm(pat)
 {
   setMomenta(mumAndDgtrs_p3);
   if(_pat.size() < 3 || _pat.size() != mumAndDgtrs_p3.size()) shoutAndKill();
   resetST();
   _eventCounter++;
 }

◆ DalitzEvent() [6/12]

DalitzEvent::DalitzEvent ( const DalitzEventPattern & pat )

Definition at line 232 of file DalitzEvent.cpp.

   : _pat(pat)
   , _rememberPhaseSpace(-9999.)
   , _aValue(-9999.)
   , _weight(1)
   , _generatorPdfRelativeToPhaseSpace(1)
   , _vectorOfValues()
   , _vectorOfWeights()
   , _s(pat.size())
   , _t(pat.size())
   , _permutationIndex(0)
   , _perm(pat)
 {
   if (_pat.size() < 3) shoutAndKill();
   _p.resize(_pat.size());
   resetST();
   _eventCounter++;
 } // generate or add momemta later.

◆ DalitzEvent() [7/12]

DalitzEvent::DalitzEvent	(	const DalitzEventPattern &	pat,
		TRandom *	rnd
	)

Definition at line 251 of file DalitzEvent.cpp.

   : _pat(pat)
   , _rememberPhaseSpace(-9999.)
   , _aValue(-9999.)
   , _weight(1)
   , _generatorPdfRelativeToPhaseSpace(1)
   , _vectorOfValues()
   , _vectorOfWeights()
   , _s(pat.size())
   , _t(pat.size())
   , _permutationIndex(0)
   , _perm(pat)
 {
   if (_pat.size() < 3) shoutAndKill();
   resetST();
   generateThisToPhaseSpace(rnd);
   _eventCounter++;
 }

◆ DalitzEvent() [8/12]

DalitzEvent::DalitzEvent ( TNtupleD * ntp )

Definition at line 365 of file DalitzEvent.cpp.

   : _rememberPhaseSpace(-9999.)
   , _aValue(-9999.)
   , _weight(1)
   , _generatorPdfRelativeToPhaseSpace(1)
   , _vectorOfValues()
   , _vectorOfWeights()
   , _permutationIndex(0)
 {
   if(! fromTree(ntp)){
     cout << "ERROR in DalitzEvent constructor from ntuple"
      << " something went wrong!"
      << endl;
   }
   _eventCounter++;
 }

◆ DalitzEvent() [9/12]

DalitzEvent::DalitzEvent ( const IDalitzEvent & ievt )

Definition at line 293 of file DalitzEvent.cpp.

   : _pat(other.eventPattern())
   , _rememberPhaseSpace(other.phaseSpace())
   , _rememberComplexFast()
   , _rememberDoubleFast()
   , _aValue(other.getAValue())
   , _weight(other.getWeight())
   , _generatorPdfRelativeToPhaseSpace(other.getGeneratorPdfRelativeToPhaseSpace())
   , _vectorOfValues(other.getVectorOfValues())
   , _vectorOfWeights(other.getVectorOfWeights())
   , _permutationIndex(0)
   , _perm(other.eventPattern())
 {
   int np = other.eventPattern().size();
   _p.resize(np);
   for(int i=0; i < np; i++) _p[i] = other.p(i);
   resetST();
   _eventCounter++;
 }

◆ DalitzEvent() [10/12]

DalitzEvent::DalitzEvent ( const IDalitzEvent * ievt )

Definition at line 270 of file DalitzEvent.cpp.

   : _pat(other->eventPattern())
     //  , _p(other.)
   , _rememberPhaseSpace(other->phaseSpace())
   , _rememberComplexFast()
   , _rememberDoubleFast()
   , _aValue(other->getAValue())
   , _weight(other->getWeight())
   , _generatorPdfRelativeToPhaseSpace(other->getGeneratorPdfRelativeToPhaseSpace())
   , _vectorOfValues(other->getVectorOfValues())
   , _vectorOfWeights(other->getVectorOfWeights())
     //  , _s()
     //, _t(other._t)
   , _permutationIndex(0)
   , _perm(other->eventPattern())
 {
   int np = other->eventPattern().size();
   _p.resize(np);
   for(int i=0; i < np; i++) _p[i] = other->p(i);
   resetST();
   _eventCounter++;
 }

◆ DalitzEvent() [11/12]

DalitzEvent::DalitzEvent ( const DalitzEvent * other )

Definition at line 315 of file DalitzEvent.cpp.

   : IDalitzEvent() // the compiler wants it.
   , _pat(other->_pat)
   , _p(other->_p)
   , _rememberPhaseSpace(other->_rememberPhaseSpace)
   , _rememberComplexFast(other->_rememberComplexFast)
   , _rememberDoubleFast(other->_rememberDoubleFast)
   , _aValue(other->_aValue)
   , _weight(other->_weight)
   , _generatorPdfRelativeToPhaseSpace(other->_generatorPdfRelativeToPhaseSpace)
   , _vectorOfValues(other->_vectorOfValues)
   , _vectorOfWeights(other->_vectorOfWeights)
   , _s(other->_s)
   , _t(other->_t)
   , _sijMap(other->_sijMap)
   , _permutationIndex(other->_permutationIndex)
   , _perm(other->_perm)
 {
   bool dbThis=false;
   if(dbThis) cout << " copy constructor of DalitzEvent called."
           << " weight before " << other->getWeight()
           << " and after " << this->getWeight()
           << endl;
   _eventCounter++;
 }

◆ DalitzEvent() [12/12]

DalitzEvent::DalitzEvent ( const DalitzEvent & other )

Definition at line 340 of file DalitzEvent.cpp.

   : IDalitzEvent() // the compiler wants it.
   , _pat(other._pat)
   , _p(other._p)
   , _rememberPhaseSpace(other._rememberPhaseSpace)
   , _rememberComplexFast(other._rememberComplexFast)
   , _rememberDoubleFast(other._rememberDoubleFast)
   , _aValue(other._aValue)
   , _weight(other._weight)
   , _generatorPdfRelativeToPhaseSpace(other._generatorPdfRelativeToPhaseSpace)
   , _vectorOfValues(other._vectorOfValues)
   , _vectorOfWeights(other._vectorOfWeights)
   , _s(other._s)
   , _t(other._t)
   , _sijMap(other._sijMap)
   , _permutationIndex(other._permutationIndex)
   , _perm(other._perm)
 {
   bool dbThis=false;
   if(dbThis) cout << " copy constructor of DalitzEvent called."
           << " weight before " << other.getWeight()
           << " and after " << this->getWeight()
           << endl;
   _eventCounter++;
 }

◆ ~DalitzEvent()

DalitzEvent::~DalitzEvent ( )

virtual

Definition at line 382 of file DalitzEvent.cpp.

                          {
   _eventCounter--;
   //  cout << "DalitzEvent destroyed\n" << endl;
 }

Member Function Documentation

◆ assignUniqueRememberNumber()

long int DalitzEvent::assignUniqueRememberNumber ( )

static

Definition at line 29 of file DalitzEvent.cpp.

                                                 {
   return _rememberVectorCounter++;
 }

◆ BDet()

double DalitzEvent::BDet ( ) const

protected

Definition at line 883 of file DalitzEvent.cpp.

                              {
   /*
   cout << " for sij " << t(0,1) << ", " << s(1,2) << ", " << s(2,3)
        << ", " << s(3,4) << ", " << t(4,0) << endl;
   cout << "\n " << m2(0) << ", " << m2(1) << ", " << m2(2)
        << ", " << m2(3) << m2(4) << endl;
   */
 
   if(nDgtr() < 4){
     cout << "ERROR in DalitzEvent::BDet(): nDgtr() = "
      << nDgtr() << " < 4. How did this happen?"
      << endl;
     return +9999;
   }
 
   TMatrixDSym B(6);
   B(0,0)=0;
   B(1,0)=B(0,1)=m2(2); B(1,1) = 0;
   B(2,0)=B(0,2)=s(2,3) ; B(2,1)=B(1,2) = m2(3); B(2,2)=0;
   B(3,0)=B(0,3)=t(0,1) ; B(3,1)=B(1,3) = s(3,4) ; B(3,2)=B(2,3)=m2(4); B(3,3)=0;
   B(4,0)=B(0,4)=m2(1); B(4,1)=B(1,4) = s(1,2) ; B(4,2)=B(2,4)=t(4,0) ; B(4,3)=B(3,4)=m2(0); B(4,4)=0;
   B(5,0)=B(0,5)=B(5,1)=B(1,5)=B(5,2)=B(2,5)=B(5,3)=B(3,5)=B(5,4)=B(4,5)=1; 
   B(5,5)=0;
 
   /*
   double Bdet = B.Determinant();
   cout << "Bdet = " << Bdet << endl;
   */
 
   return B.Determinant();
 }

◆ C_conjugateYourself()

void DalitzEvent::C_conjugateYourself ( )

virtual

Definition at line 808 of file DalitzEvent.cpp.

                                      {
   _pat = _pat.makeCPConjugate();
 }

◆ clone()

IDalitzEvent * DalitzEvent::clone ( ) const

virtual

Implements IDalitzEvent.

Definition at line 387 of file DalitzEvent.cpp.

                                       {
   return new DalitzEvent(*this);
 }

◆ CP_conjugateYourself()

void DalitzEvent::CP_conjugateYourself ( )

virtual

Definition at line 811 of file DalitzEvent.cpp.

                                       {
   P_conjugateYourself();
   C_conjugateYourself();
 }

◆ currentPermutation()

const Permutation& DalitzEvent::currentPermutation ( ) const

inline

Definition at line 148 of file DalitzEvent.h.

                                                {
     return _perm[_permutationIndex];
   }

◆ Delta4()

double DalitzEvent::Delta4 ( ) const

protected

Definition at line 914 of file DalitzEvent.cpp.

                                 {
   return -BDet()/16;
 }

◆ evalsij_intern()

double DalitzEvent::evalsij_intern ( const std::vector< int > & intern_indices ) const

protected

Definition at line 699 of file DalitzEvent.cpp.

                                         {
   TLorentzVector psum(0.0, 0.0, 0.0, 0.0);
   for(unsigned int i=0; i< intern_indices.size(); i++){
     //cout << " p(" << i << ") " << p(indices[i]);
     psum += p_intern(intern_indices[i]);
   }
   return psum.Mag2();
 }

◆ eventCounter()

static long int DalitzEvent::eventCounter ( )

inlinestatic

Definition at line 125 of file DalitzEvent.h.

125 { return _eventCounter;}

DalitzEvent::_eventCounter

static long int _eventCounter

Definition: DalitzEvent.h:37

◆ eventPattern()

virtual const DalitzEventPattern& DalitzEvent::eventPattern ( ) const

inlinevirtual

Implements IDalitzEvent.

Definition at line 173 of file DalitzEvent.h.

173 {return _pat;}

DalitzEvent::_pat

DalitzEventPattern _pat

Definition: DalitzEvent.h:42

◆ fillNtupleVarArray()

bool DalitzEvent::fillNtupleVarArray ( std::vector< Double_t > & array ) const

Definition at line 1014 of file DalitzEvent.cpp.

                                                                    {
 
   if(eventPattern().empty()){
     return false;
   }
   int counter=0;
   int maxCounter = (array.size() - 2) - (singleParticleNtpArraySize()-1);
   for(unsigned int i=0; i< _p.size() && counter < maxCounter; i++){
     array.at(counter++) = _pat[i];
     array.at(counter++) = _p[i].E();
     array.at(counter++) = _p[i].X();
     array.at(counter++) = _p[i].Y();
     array.at(counter++) = _p[i].Z();
   }
   array.at(counter++) = getWeight();
   array.at(counter++) = getGeneratorPdfRelativeToPhaseSpace();
   return true;
 }

◆ fromNtuple()

bool DalitzEvent::fromNtuple ( TNtupleD * ntp )

Definition at line 1061 of file DalitzEvent.cpp.

                                          {
   // assumes ntuple is set to the correct
   // entry (using ntp->GetEntry(int i))
   bool dbThis=false;
   if(dbThis) cout << "DalitzEvent::fromNtuple(" << ntp->GetName() << ") called" << endl;
   int arraySize = ntp->GetNvar();
   //  int arraySize = ntp->GetNbranches();
   if(dbThis) cout << " arraySize " << arraySize << endl;
   bool newFormat = false;
   if(0 == (arraySize-2)%singleParticleNtpArraySize() ){
     if(dbThis) cout << "New Format" << endl;
     newFormat = true; // includes weight & generatorPdf
   }else if( 0 == (arraySize)%singleParticleNtpArraySize()){
     newFormat = false; // excludes weight & generatorPdf
   }else{
     cout << "ERROR in DalitzEvent::fromNtuple length of array"
      << " in ntuple = " << arraySize
      << "\n   >  not a multiple of # entries per particle = "
      << singleParticleNtpArraySize()
      << endl;
     return false;
   }
   const Double_t* array = ntp->GetArgs();
 
   if(dbThis) cout << " array ptr " << array << endl;
   int counter = 0;
   int maxCounter = arraySize - (singleParticleNtpArraySize()-1);
   if(newFormat) maxCounter -= 2;
   int nParticles;
   if(newFormat){
     nParticles = (arraySize-2)/singleParticleNtpArraySize();
   }else{
     nParticles = arraySize/singleParticleNtpArraySize();
   }
   if(dbThis){
     cout << "nParticles = " << nParticles
      << " arraySize = " << arraySize
      << endl;
   }
 
   _pat.resize(nParticles);
   _p.resize(nParticles);
   for(int i=0; i< nParticles && counter < maxCounter; i++){
     if(dbThis) cout << "filling " << i << " th particle: " << endl;
     _pat[i] = nearestInt(array[counter++]);
     _p[i].SetE(array[counter++]);
     _p[i].SetX(array[counter++]);
     _p[i].SetY(array[counter++]);
     _p[i].SetZ(array[counter++]);
     
     // check units (should be MeV, but will recognise if its GeV and 
     // translate it to MeV)
 
     double units(1);
     double dMeV = abs(_pat[i].mass() - _p[i].M()*MeV);
     double dGeV = abs(_pat[i].mass() - _p[i].M()*GeV);
     if(dMeV < dGeV) units=MeV;
     else units=GeV;
     _p[i] *= units;
     
     if(dbThis){
       cout << "\t" << _pat[i] << " " << _p[i] 
        << ", m= " << _p[i].M() << endl;
     }
   }
   if(newFormat){
     setWeight(array[counter++]);
     setGeneratorPdfRelativeToPhaseSpace(array[counter++]);
   }
   if(dbThis){
     cout << " weight, pdf "
      << getWeight() << ", "
      << getGeneratorPdfRelativeToPhaseSpace() << endl;
   }
   resetST();
   _perm.setPattern(_pat);
   if(dbThis){
     cout << "this is me now " << endl;
     this->print();
     cout << "done it all, returning" << endl;
   }
   return true;
 }

◆ fromParasTree()

bool DalitzEvent::fromParasTree ( TTree * ntp )

Definition at line 1271 of file DalitzEvent.cpp.

                                          {
   bool dbThis=false;
   bool forcePDGMassForFinalState=true;
 
   if(dbThis) cout << "DalitzEvent::fromTREE(" << ntp << ") called" << endl;
   
   int arraySize = ntp->GetNbranches();
   if(dbThis) cout << " arraySize " << arraySize << endl;
   const TObjArray* branchArray = ntp->GetListOfBranches();
 
   map<string, map<std::string, double> > prtMap;
 
   double weight=1;
   double genPdf=1;
   for(int i=0; i< arraySize; i++){
     TObjArray* leafArray=
       ((TBranch*) (*branchArray)[i])->GetListOfLeaves();
     double leafVal = ((TLeaf*) (*leafArray)[0])->GetValue();
     std::string leafName = ((TLeaf*) (*leafArray)[0])->GetName();
     
     std::string parName, element;
     parseNtpEntryName(leafName, parName, element);
     if(dbThis){
       cout << " leaf[ " << i << "] " 
        << leafName << "; " << leafVal << endl;
     }
 
     if("weight" == leafName) weight=leafVal;
     else if("genPdf" == leafName) genPdf = leafVal;
     else{
       std::string parName, element;
       if(!parseNtpEntryName(leafName, parName, element)){
     // leafName doesn't contain an underscore.
     continue ;
       }
       (prtMap[parName])[element] = leafVal;
       if(dbThis){
     cout << "\t leaf[ " << i << "] name parsed: " 
          << parName << ", " << element
          << endl;
       }
     }    
   }
   int nParticles = prtMap.size();
   if(dbThis) cout << "nParticles = " << nParticles << endl;
   _pat.resize(nParticles);
   _p.resize(nParticles);
   
   int counter=0;
   for(map<string, map<string, double> >::iterator it = 
     prtMap.begin();
       it != prtMap.end();
       it++, counter++){
     _pat[counter] = (it->second)["pdg"];
     if(dbThis){
       cout << "just filled _pat[" << counter << "] with: "
        << (it->second)["pdg"] << "; see: " << _pat[counter]
        << endl;
     }
     _p[counter].SetE((it->second)["E"]);
     _p[counter].SetPx((it->second)["Px"]);
     _p[counter].SetPy((it->second)["Py"]);
     _p[counter].SetPz((it->second)["Pz"]);
 
     // check units (Paras usually GeV, else usually MeV)
     double units;
     double dMeV = abs(_pat[counter].mass() - _p[counter].M()*MeV);
     double dGeV = abs(_pat[counter].mass() - _p[counter].M()*GeV);
     if(dMeV < dGeV) units=MeV;
     else units=GeV;
     _p[counter] *= units;
     
     if(forcePDGMassForFinalState && counter > 0){
       double m  = _pat[counter].mass();
       double p3 = _p[counter].Rho();
       _p[counter].SetE(sqrt(m*m + p3*p3));
     }
   }
   setWeight(weight);
   setGeneratorPdfRelativeToPhaseSpace(genPdf);
   if(dbThis){
     cout << " weight, pdf "
      << getWeight() << ", "
      << getGeneratorPdfRelativeToPhaseSpace() << endl;
     for(int i=0; i < nParticles; i++){
       cout << _pat[i] << " ) " << _p[i] << " m = " << _p[i].M() << endl;
     }
   }
   resetST();
   _perm.setPattern(_pat);
   if(dbThis) cout << "done it all, returning" << endl;
   return true;
   }

◆ fromParasTreeOld()

bool DalitzEvent::fromParasTreeOld ( TTree * ntp )

Definition at line 1144 of file DalitzEvent.cpp.

                                             {
   bool dbThis=false;
   bool forcePDGMassForFinalState=true;
 
   if(dbThis) cout << "DalitzEvent::fromTREE(" << ntp << ") OLD called" << endl;
   
   int arraySize = ntp->GetNbranches();
   if(dbThis) cout << " arraySize " << arraySize << endl;
   bool newFormat = false;
   if(0 == (arraySize-2)%singleParticleNtpArraySize() ){
     newFormat = true; // includes weight & generatorPdf
     if(dbThis) cout << " new Format" << endl;
   }else if( 0 == (arraySize)%singleParticleNtpArraySize()){
     newFormat = false; // excludes weight & generatorPdf
     if(dbThis) cout << " old Format" << endl;
   }else{
     cout << "ERROR in DalitzEvent::fromTree length of array"
      << " in ntuple = " << arraySize
      << "\n   >  not a multiple of # entries per particle = "
      << singleParticleNtpArraySize()
      << endl;
     return false;
   }
   const TObjArray* branchArray = ntp->GetListOfBranches();
 
   if(dbThis) cout << " array ptr " << branchArray << endl;
   int counter = 0;
   int maxCounter = arraySize - (singleParticleNtpArraySize()-1);
   if(newFormat) maxCounter -= 2;
   int nParticles;
   if(newFormat){
     nParticles = (arraySize-2)/singleParticleNtpArraySize();
   }else{
     nParticles = arraySize/singleParticleNtpArraySize();
   }
   if(dbThis){
     cout << "nParticles = " << nParticles
      << " arraySize = " << arraySize
      << endl;
   }
 
   _pat.resize(nParticles);
   _p.resize(nParticles);
   for(int i=0; i< nParticles && counter < maxCounter; i++){
     if(dbThis) cout << "filling " << i << " th particle: " << endl;
     TObjArray* leafArray;
     double leafVal = -9999;
     std::string leafName;
     leafArray=
       ((TBranch*) (*branchArray)[counter++])->GetListOfLeaves();
     leafVal = ((TLeaf*) (*leafArray)[0])->GetValue();
     leafName = ((TLeaf*) (*leafArray)[0])->GetName();
     if(dbThis) cout << " leaf[ " << counter-1 << "] " 
             << leafVal << endl;
     _p[i].SetE(leafVal*GeV);
     leafArray = 
       ((TBranch*) (*branchArray)[counter++])->GetListOfLeaves();
     leafVal = ((TLeaf*) (*leafArray)[0])->GetValue();
     leafName = ((TLeaf*) (*leafArray)[0])->GetName();
     if(dbThis){
       cout << " leaf[ " << counter-1 << "] " << leafName 
        << ": " << leafVal << endl;
     }
     _p[i].SetX(leafVal*GeV);
     leafArray = 
       ((TBranch*) (*branchArray)[counter++])->GetListOfLeaves();
     leafVal = ((TLeaf*) (*leafArray)[0])->GetValue();
     if(dbThis) cout << " leaf[ " << counter-1 << "] " << leafVal << endl;
     _p[i].SetY(leafVal*GeV);
     leafArray = 
       ((TBranch*) (*branchArray)[counter++])->GetListOfLeaves();
     leafVal = ((TLeaf*) (*leafArray)[0])->GetValue();
     if(dbThis) cout << " leaf[ " << counter-1 << "] " << leafVal << endl;
     _p[i].SetZ(leafVal*GeV);
     leafArray = 
       ((TBranch*) (*branchArray)[counter++])->GetListOfLeaves();
     leafVal = ((TLeaf*) (*leafArray)[0])->GetValue();
     if(dbThis) cout << " leaf[ " << counter-1 << "] " << leafVal << endl;
     _pat[i] = nearestInt(leafVal);
 
     if(forcePDGMassForFinalState && i > 0){
       double m=_pat[i].mass();
       _p[i].SetE(sqrt(m*m + _p[i].Rho()*_p[i].Rho()));
     }
     if(dbThis) cout << " got pattern " << (int) _pat[i] << endl;
 
     if(dbThis){
       cout << "\t" << _pat[i] << " " << _p[i] 
        << ", m= " << _p[i].M() << endl;
     }
   }
   if(newFormat){
     TObjArray* leafArray;
     double leafVal = -9999;
     leafArray = 
       ((TBranch*) (*branchArray)[counter++])->GetListOfLeaves();
     leafVal = ((TLeaf*) (*leafArray)[0])->GetValue();
     if(dbThis) cout << " leaf[ " << counter-1 << "] " << leafVal << endl;
     setWeight(leafVal);
 
     leafArray = 
       ((TBranch*) (*branchArray)[counter++])->GetListOfLeaves();
     leafVal = ((TLeaf*) (*leafArray)[0])->GetValue();
     setGeneratorPdfRelativeToPhaseSpace(leafVal);
   }
   if(dbThis){
     cout << " weight, pdf "
      << getWeight() << ", "
      << getGeneratorPdfRelativeToPhaseSpace() << endl;
   }
   resetST();
   _perm.setPattern(_pat);
   if(dbThis) cout << "done it all, returning" << endl;
   return true;
 }

◆ fromTree()

bool DalitzEvent::fromTree ( TTree * tree )

Definition at line 1047 of file DalitzEvent.cpp.

                                      {
   bool dbThis=false;
   // figure out format: TNtuple or what?
   if(0 == tree) return 0;
   std::string treeName(tree->ClassName());
   if(dbThis) cout << "ClassName \"" << treeName << "\"" << endl;
   if(treeName == (std::string) "TNtupleD"){
     if(dbThis)cout << "DalitzEvent::fromTree: tree is TNtupleD" << endl;
     return fromNtuple((TNtupleD*)tree);
   }else{
     if(dbThis) cout << "DalitzEvent::fromTree: tree is ParasTree" << endl;
     return fromParasTree(tree);
   }
 }

◆ generateThisToPhaseSpace() [1/3]

void DalitzEvent::generateThisToPhaseSpace ( TRandom * rnd = 0 )

Definition at line 546 of file DalitzEvent.cpp.

                                                       {
   generateThisToPhaseSpace(0.0, rnd);
 }

◆ generateThisToPhaseSpace() [2/3]

void DalitzEvent::generateThisToPhaseSpace	(	double	p_of_D,
		TRandom *	rnd = `0`
	)

Definition at line 549 of file DalitzEvent.cpp.

                                                                 {
   TVector3 p3(0,0,p);
   generateThisToPhaseSpace(p3, rnd);
 }

◆ generateThisToPhaseSpace() [3/3]

void DalitzEvent::generateThisToPhaseSpace	(	const TVector3 &	pVec_of_D,
		TRandom *	rnd = `0`
	)

Definition at line 553 of file DalitzEvent.cpp.

                                       {
   bool dbthis=false;
   if(dbthis) cout << "got called" << endl;
   TRandom* remember_gRandom = gRandom;
   bool needToReset_gRandom=false;
   if(0 != rnd){
     needToReset_gRandom=true;
     gRandom = rnd; // horrible fix because
     // TGenPhaseSpace can only use gRandom
   }
   if(dbthis) cout << "got so far" << endl;
   TGenPhaseSpace gps; // uses gRandom, so no worries
   // about re-starting with the same seed (won't happen), 
   // but worries about what kind of
   // random number generator gRandom is.
 
   TLorentzVector mumsP4;
   mumsP4.SetXYZM(mumsp3.X()
          , mumsp3.Y()
          , mumsp3.Z()
          , _pat[0].mass());
 
   Double_t* mdgt = new Double_t[nDgtr()];
   if(dbthis) cout << "booked " << nDgtr() << " daughters " << endl;
   for(int i=1; i<= nDgtr(); i++){
     mdgt[i-1] = this->m(i);
   }
   
   gps.SetDecay(mumsP4, nDgtr() , mdgt);
   double maxWeight = 1.0;//gps.GetWtMax();
   double weight=0;
   if(dbthis) cout << " maxWeight " << maxWeight << endl;
   //normal:
   do{
     weight=gps.Generate();
   }while(weight < gRandom->Rndm() * maxWeight);
   // for debug:
   //  weight=gps.Generate();
 
   if(needToReset_gRandom){
     gRandom = remember_gRandom;
   }
   if(weight > maxWeight){
     cout << "ERROR in DalitzEvent::generateThisToPhaseSpace"
      << " weight = " << weight << " > " << maxWeight << " = maxWeight"
      << "\n    > this is a programming error that needs to be fixed."
      << endl;
     throw "...so I got maxWeight wrong.";
   }
   if(dbthis) cout << " this weight " << weight << endl;
   if(dbthis) cout << " max weight " << gps.GetWtMax() << endl;
 
   _p.resize(nDgtr()+1);
   if(dbthis) cout << "generated event" << endl;
   setP(0, mumsP4);
   if(dbthis)cout << "set p0" << endl;
   for(unsigned int i=1; i< _pat.size(); i++){
     if(dbthis)cout << "setting P (" << i << ") " 
            << (*(gps.GetDecay(i-1)))
            << "..." << endl;
     setP(i, (*(gps.GetDecay(i-1))));
     if(dbthis)cout << " and look what I've set: " << p(i) << endl;
   }
   delete[] mdgt;
   _weight = 1.0;
   _generatorPdfRelativeToPhaseSpace = 1.0;
   return;
 }

◆ getAValue()

virtual double DalitzEvent::getAValue ( ) const

inlinevirtual

Implements IDalitzEvent.

Definition at line 171 of file DalitzEvent.h.

171 {return _aValue;}

DalitzEvent::_aValue

double _aValue

Definition: DalitzEvent.h:50

◆ getGeneratorPdfRelativeToPhaseSpace()

double DalitzEvent::getGeneratorPdfRelativeToPhaseSpace ( ) const

virtual

Implements IDalitzEvent.

Definition at line 396 of file DalitzEvent.cpp.

                                                              {
   return _generatorPdfRelativeToPhaseSpace;
 }

◆ getValueFromVector()

double DalitzEvent::getValueFromVector ( unsigned int i ) const

virtual

Implements IDalitzEvent.

Definition at line 425 of file DalitzEvent.cpp.

                                                           {
     return _vectorOfValues[i];
 }

◆ getVectorOfValues() [1/2]

const std::vector< double > & DalitzEvent::getVectorOfValues ( ) const

virtual

Implements IDalitzEvent.

Definition at line 403 of file DalitzEvent.cpp.

403 {

404 return _vectorOfValues;}

DalitzEvent::_vectorOfValues

std::vector< double > _vectorOfValues

Definition: DalitzEvent.h:54

◆ getVectorOfValues() [2/2]

std::vector< double > & DalitzEvent::getVectorOfValues ( )

virtual

Implements IDalitzEvent.

Definition at line 406 of file DalitzEvent.cpp.

406 {

407 return _vectorOfValues;}

DalitzEvent::_vectorOfValues

std::vector< double > _vectorOfValues

Definition: DalitzEvent.h:54

◆ getVectorOfWeights() [1/2]

const std::vector< double > & DalitzEvent::getVectorOfWeights ( ) const

virtual

Implements IDalitzEvent.

Definition at line 409 of file DalitzEvent.cpp.

409 {

410 return _vectorOfWeights;}

DalitzEvent::_vectorOfWeights

std::vector< double > _vectorOfWeights

Definition: DalitzEvent.h:55

◆ getVectorOfWeights() [2/2]

std::vector< double > & DalitzEvent::getVectorOfWeights ( )

virtual

Implements IDalitzEvent.

Definition at line 412 of file DalitzEvent.cpp.

412 {

413 return _vectorOfWeights;}

DalitzEvent::_vectorOfWeights

std::vector< double > _vectorOfWeights

Definition: DalitzEvent.h:55

◆ getWeight()

double DalitzEvent::getWeight ( ) const

virtual

Implements IDalitzEvent.

Definition at line 390 of file DalitzEvent.cpp.

                                    {
   return _weight;
 }

◆ getWeightFromVector()

double DalitzEvent::getWeightFromVector ( unsigned int i ) const

virtual

Implements IDalitzEvent.

Definition at line 429 of file DalitzEvent.cpp.

                                                            {
     return _vectorOfWeights[i];
 }

◆ kinematicallyAllowed()

bool DalitzEvent::kinematicallyAllowed ( double epsilon = 1.e-9 ) const

Definition at line 466 of file DalitzEvent.cpp.

                                                           {
   bool dbThis=false;
 
   TLorentzVector pdiff = (p(0) - p_dgtr_sum());
   TLorentzVector psum  = (p(0) + p_dgtr_sum());
 
   double diffsum = 0;
   double sumsum  = 0;
   for(unsigned int i=0; i<4; i++){
     diffsum += pdiff[i]*pdiff[i];
     sumsum  += psum[i] *psum[i];
   }
   double avsum = sumsum/2.0;
 
   if(diffsum > epsilon*epsilon * avsum){
     if(dbThis){
       cout << "DalitzEvent::kinematicallyAllowed? "
        << " failing diffsum: " << diffsum
        << " max allowed: " << epsilon*epsilon*avsum << endl;
     } 
     return false;
   }
 
   for(unsigned int i = 0; i<_pat.size(); i++){
     double mreco = p(i).M();
     if(mreco < 0){
       if(dbThis) cout << "failing mreco " << mreco << " < 0" << endl;
       return false;
     }
     double avmass = (mreco + m(i))/2.0;
     if( fabs(mreco - m(i)) > epsilon*avmass ){
       if(dbThis){
     cout << " failing mreco - m(i) fabs(" 
          << mreco << " - " << m(i) << ") = " 
          << fabs(mreco - m(i)) << " > " << epsilon*avmass
          << endl;
       }
       return false;
     }
   }
   return true;
 }

◆ m()

double DalitzEvent::m ( unsigned int i ) const

Definition at line 433 of file DalitzEvent.cpp.

                                          {
   // easy access to the nominal masses
   if(i >= _pat.size()){
     cout << "ERROR in DalitzEvent::m(unsigned int i)"
      << " index i= " << i << " out of range!"
      << " allowed range: [0," << _pat.size()-1 << "]"
      << " returning -9999"
      << endl;
     return -9999;
   }
   return _pat[i].mass();
 }

◆ m2()

double DalitzEvent::m2 ( unsigned int i ) const

Definition at line 445 of file DalitzEvent.cpp.

                                           {
   if(i >= _pat.size()){
     cout << "ERROR in DalitzEvent::m2(unsigned int i)"
      << " index i= " << i << " out of range!"
      << " allowed range: [0," << _pat.size()-1 << "]"
      << " returning -9999"
      << endl;
     return -9999;
   }  
   return m(i)*m(i);
 }

◆ makeCalculate4BodyProps()

Calculate4BodyProps DalitzEvent::makeCalculate4BodyProps ( ) const

Definition at line 751 of file DalitzEvent.cpp.

                                                               {
   if(nDgtr() != 4){
     cout << "ERROR in DalitzEvent::makeCalculate4BodyProps()"
      << "\n   > You can't expect a sensible Calculate4BodyProps"
      << " Object for a " << nDgtr() << " body decay!"
      << endl;
     return Calculate4BodyProps(0,0,0,0,0, 0,0,0,0,0);
   }
   return Calculate4BodyProps(t(0,1), s(1,2), s(2,3), s(3,4), t(4,0)
                  , m(0), m(1), m(2), m(3), m(4)
                  );
 }

◆ makeNtupleVarnames()

std::string DalitzEvent::makeNtupleVarnames ( ) const

Definition at line 987 of file DalitzEvent.cpp.

                                               {
   std::string s="";
   if(eventPattern().empty()){
     return s;
   }
 
   for(unsigned int i=0; i< _pat.size(); i++){
     std::string name = _pat[i].name();
     name = "_" + anythingToString((int) i) + "_" + name;
     name = prtToNtpName(name);
 
     s += name + "_pdg:";
     s += name + "_E:";
     s += name + "_Px:";
     s += name + "_Py:";
     s += name + "_Pz";
     //    if(i != _pat.size() -1) s+= ":";
     s+= ":";
   }
   s += "weight:";
   s += "genPdf";
   return s;
 }

◆ nDgtr()

int DalitzEvent::nDgtr ( ) const

inline

Definition at line 222 of file DalitzEvent.h.

                    {
     return _pat.size()-1;
   }

◆ ntpToPrtName()

std::string DalitzEvent::ntpToPrtName ( const std::string & s_in )

staticprotected

Definition at line 979 of file DalitzEvent.cpp.

                                                         {
   std::string s(s_in);
   for(int i=0; ntpNameChars[i] != '\0'; i++){
     replace(s.begin(), s.end(), ntpNameChars[i], prtNameChars[i]);
   }
   return s;
 }

◆ ntupleVarArraySize()

unsigned int DalitzEvent::ntupleVarArraySize ( ) const

Definition at line 1032 of file DalitzEvent.cpp.

                                                   {
   return _p.size() * singleParticleNtpArraySize() + 2;
 }

◆ numPermutations()

int DalitzEvent::numPermutations ( ) const

virtual

Implements IDalitzEvent.

Definition at line 843 of file DalitzEvent.cpp.

                                       {
   return _perm.size();
 }

◆ p() [1/2]

const TLorentzVector & DalitzEvent::p ( unsigned int i ) const

virtual

Implements IDalitzEvent.

Definition at line 673 of file DalitzEvent.cpp.

                                                         {
   bool dbThis=false;
   if(dbThis){
     cout << "p(" << i << ") called. Current permutation: "
      << currentPermutation()[i]
      << endl;
   }
   return p_intern(currentPermutation()[i]);
 }

◆ p() [2/2]

TLorentzVector & DalitzEvent::p ( unsigned int i )

virtual

Definition at line 682 of file DalitzEvent.cpp.

                                             {
   bool dbThis=false;
   if(dbThis){
     cout << "p(" << i << ") called. Current permutation: "
      << currentPermutation()[i]
      << endl;
   }
   return p_intern(currentPermutation()[i]);
 }

◆ P_conjugateYourself()

void DalitzEvent::P_conjugateYourself ( )

virtual

Definition at line 795 of file DalitzEvent.cpp.

                                      {
   // p-conjugates in mums restframe, but keeps old D momentum.
 
   resetST(); // shouldn't be necessary, but to be save.
   if(_p.empty()) return;
   TVector3 mums3Momentum(_p[0].X(), _p[0].Y(), _p[0].Z());
   for(unsigned int i=0; i < _p.size(); i++){
     _p[i].SetX( - _p[i].X() );
     _p[i].SetY( - _p[i].Y() );
     _p[i].SetZ( - _p[i].Z() );
   }
   setMothers3Momentum(mums3Momentum);
 }

◆ p_dgtr_sum()

TLorentzVector DalitzEvent::p_dgtr_sum ( ) const

Definition at line 457 of file DalitzEvent.cpp.

                                             {
   TLorentzVector pd(0.0, 0.0, 0.0, 0.0);
 
   for(unsigned int i=1; i< _p.size(); i++){
     pd += _p[i];
   }
   return pd;
 }

◆ p_intern() [1/2]

const TLorentzVector & DalitzEvent::p_intern ( unsigned int i ) const

protected

Definition at line 623 of file DalitzEvent.cpp.

                                                                {
   if(i >= _p.size()){
     cout << "in TLorentzVector& DalitzEvent::p(unsigned int i):"
      << " index out of range. DalitzEventPattern:"
      << "\n  > \"" << eventPattern() << "\""
      << "\n  > num mother + daughters = " << _p.size()
      << endl;
     throw "index out or range in  TLorentzVector& DalitzEvent::p(unsigned int i)";
   }
   return _p[i];
 }

◆ p_intern() [2/2]

TLorentzVector & DalitzEvent::p_intern ( unsigned int i )

protected

Definition at line 635 of file DalitzEvent.cpp.

                                                    {
   if(i >= _p.size()){
     cout << "in TLorentzVector& DalitzEvent::p(unsigned int i):"
      << " index out of range. DalitzEventPattern:"
      << "\n  > \"" << eventPattern() << "\""
      << "\n  > num mother + daughters = " << _p.size()
      << endl;
     throw "index out or range in  TLorentzVector& DalitzEvent::p(unsigned int i)";
   }
   return _p[i];
 }

◆ parseNtpEntryName()

bool DalitzEvent::parseNtpEntryName	(	const std::string &	entry,
		std::string &	part1,
		std::string &	part2
	)

staticprotected

Definition at line 1260 of file DalitzEvent.cpp.

                                        {
   int posOf_ = entry.find_last_of('_');
   if(std::string::npos == posOf_)
     return false ;
   part1 = entry.substr(0, posOf_);
   part2 = entry.substr(posOf_ +1);
   return true;
 }

◆ permutationIndex()

int DalitzEvent::permutationIndex ( ) const

inlinevirtual

Implements IDalitzEvent.

Definition at line 154 of file DalitzEvent.h.

                               {
     return _permutationIndex;
   }

◆ phaseSpace()

double DalitzEvent::phaseSpace ( ) const

virtual

Implements IDalitzEvent.

Definition at line 510 of file DalitzEvent.cpp.

                                     {
   if(nDgtr() <= 2){
     _rememberPhaseSpace = 1;
     return 1;
   }
   if(_p.size() != _pat.size()){
     cout << "ERROR in DalitzEvent::phaseSpace: inconsistency!"
      << " pattern array size != momentum array size"
      << " " << _pat.size() << " != " << _p.size()
      << ". returning -9999"
      << endl;
     return -9999;
   }
 
   if(_rememberPhaseSpace > -9998){
     //    cout << " remembering... " << _rememberPhaseSpace << endl;
     return _rememberPhaseSpace;
   }
 
   if(nDgtr() == 3){
     _rememberPhaseSpace = phaseSpace3(); // one or zero
     return _rememberPhaseSpace;
   }else if(nDgtr() == 4){
     _rememberPhaseSpace = phaseSpace4();
     return _rememberPhaseSpace;
   }else{
     cout << "ERROR in DalitzEvent::phaseSpace inconsistency!"
      << "\n   > Sorry, don't know how to calculate phase-space"
      << " factor for nDgtr = " << nDgtr()
      << "\n   >  - can only do up to nDgtr <= 4"
      << ". returning -9999"
      << endl;
     return -9999;
   }
 }

◆ phaseSpace3()

double DalitzEvent::phaseSpace3 ( double epsilon = 1.e-9 ) const

protected

Definition at line 952 of file DalitzEvent.cpp.

                                                    {
   bool dbThis=false;
   // return either one or zero.
   // one if it's kinematically allowed, zero otherwise.
   if(kinematicallyAllowed(epsilon)){
     if(dbThis) {
       cout << "DalitzEvent::phaseSpace3(" << epsilon << ")"
        << " returning 1" << endl;
     }
     //return pi*pi/(4.0*m(0)*m(0));
     return 1;
   }else{
     if(dbThis) {
       cout << "DalitzEvent::phaseSpace3(" << epsilon << ")"
        << " returning 0" << endl;
     }
     return 0;
   }
 }

◆ phaseSpace4()

double DalitzEvent::phaseSpace4 ( ) const

protected

Definition at line 917 of file DalitzEvent.cpp.

                                      { 
   bool dbThis=false;
   if(dbThis)cout << "\n\n called phaseSpace4() " << endl;
   if(nDgtr() < 4){
     cout << "ERROR in DalitzEvent::phaseSpace4: nDgtr() = "
      << nDgtr() << " < 4. How did this happen?"
      << endl;
     return -9999;
   }
   if(dbThis){
     Calculate4BodyProps c4bp(t_intern(0,1), s_intern(1,2), s_intern(2,3)
                  , s_intern(3,4), t_intern(4,0)
                  , m(0), m(1), m(2), m(3), m(4)
                  );
     
     cout << "4-body probs gets this phase space: "
      << c4bp.phaseSpaceFactor() << endl;
     cout << " ... for t01 <<... = "
      << t_intern(0,1) << ", " << s_intern(1,2) << ", " << s_intern(2,3)
      << ", " << s_intern(3,4) << ", " << t_intern(4,0)
      << endl;
   }
 
   // Note that this calculation assumes that this is a "physical" event.
   // It doesn't check if it is outside the physically allowed
   // phase space area, for which there are more requirements than
   // only Delta4 < 0.
   double invPhsp2 = -Delta4();
   if(invPhsp2 <= 0) return 0;
   double phsp = pi*pi/(32.0*m(0)*m(0)) *  1./sqrt(invPhsp2);
   if(dbThis) cout << "  ... I get: phsp = " << phsp << endl;
   if(phsp < 0) return 0;
   else return phsp;
 }

◆ print()

void DalitzEvent::print ( std::ostream & os = std::cout ) const

virtual

Implements IDalitzEvent.

Definition at line 764 of file DalitzEvent.cpp.

                                            {
   os << "Dalitz event for pattern "
      << _pat;
   for(unsigned int i=0; i<_p.size(); i++){
     cout << "\n\t p(" << i << ") " << p(i);
   }
   os << "\n\t m(i) ";
   for(unsigned int i=0; i<_p.size(); i++){
     cout << m(i) << " ; ";
   }
   if(nDgtr() == 4){
     os << "\n\t sij :"
        << t(0,1) << ", " << s(1,2) << ", " << s(2,3) 
        << ", " << s(3,4) << ", " << t(4, 0);
   }
   os << "\n\t ps = " << this->phaseSpace();
   os << "\t weight = " << this->getWeight();
   os << endl;
 }

◆ prtToNtpName()

std::string DalitzEvent::prtToNtpName ( const std::string & s_in )

staticprotected

Definition at line 972 of file DalitzEvent.cpp.

                                                         {
   std::string s(s_in);
   for(int i=0; prtNameChars[i] != '\0'; i++){
     replace(s.begin(), s.end(), prtNameChars[i], ntpNameChars[i]);
   }
   return s;
 }

◆ resetST()

bool DalitzEvent::resetST ( )

Definition at line 816 of file DalitzEvent.cpp.

                          { 
   // resets (possibly previously calculated) values for sij, tij
 
   _s.resize(_pat.size());
   for(unsigned int i=0; i< _s.size(); i++){
     _s[i].clear();
     _s[i].resize(_pat.size());
     for(unsigned int j=0; j< _s[i].size(); j++){
       _s[i][j] = -1;
     }
   }
 
   _t.resize(_pat.size());
   for(unsigned int i=0; i< _t.size(); i++){
     _t[i].clear();
     _t[i].resize(_pat.size());
     for(unsigned int j=0; j< _t[i].size(); j++){
       _t[i][j] = -1;
     }
   }
   return true;
 }

◆ retrieveValue() [1/2]

virtual bool DalitzEvent::retrieveValue	(	int	i,
		std::complex< double > &	value,
		long int	configNumber = `1`
	)

inlinevirtual

Implements IDalitzEvent.

Definition at line 187 of file DalitzEvent.h.

                                            {
     /*std::cout << "I remember these numbers: " << std::endl;
     for(unsigned int j=0; j < _rememberComplexFast.size(); j++){
       std::complex<double> tempVal;
       _rememberComplexFast.get(j, tempVal, configNumber);
       std::cout << j << ") " << tempVal << std::endl;
     }
     std::cout << "and I am returning this: " << 
       _rememberComplexFast.get(i, value, configNumber);
     std::cout << " " << i << ") " << value << std::endl;
     */
     return _rememberComplexFast.get(i, value, configNumber);
   }

◆ retrieveValue() [2/2]

virtual bool DalitzEvent::retrieveValue	(	int	i,
		double	value,
		long int	configNumber = `1`
	)

inlinevirtual

Implements IDalitzEvent.

Definition at line 206 of file DalitzEvent.h.

                                                                           {
     return _rememberDoubleFast.get(i, value, configNumber);
   }

◆ s()

double DalitzEvent::s	(	unsigned int	i,
		unsigned int	j
	)		const

virtual

Implements IDalitzEvent.

Definition at line 691 of file DalitzEvent.cpp.

                                                          {
   return s_intern(currentPermutation()[i], currentPermutation()[j]);
 }

◆ s_intern()

double DalitzEvent::s_intern	(	unsigned int	i,
		unsigned int	j
	)		const

protected

Definition at line 647 of file DalitzEvent.cpp.

                                                                 {
   if(_s[i][j] < 0) _s[i][j] = (p_intern(i) + p_intern(j)).Mag2();
   return _s[i][j];
 }

◆ setAValue()

virtual void DalitzEvent::setAValue ( double a )

inlinevirtual

Implements IDalitzEvent.

Definition at line 170 of file DalitzEvent.h.

170 {_aValue=a;}

DalitzEvent::_aValue

double _aValue

Definition: DalitzEvent.h:50

◆ setGeneratorPdfRelativeToPhaseSpace()

void DalitzEvent::setGeneratorPdfRelativeToPhaseSpace ( double gpdf )

virtual

Implements IDalitzEvent.

Definition at line 399 of file DalitzEvent.cpp.

                                                                 {
   _generatorPdfRelativeToPhaseSpace = gpdf;
 }

◆ setMomenta() [1/2]

bool DalitzEvent::setMomenta ( const std::vector< TLorentzVector > & mumAndDgtr_p4 )

protected

Definition at line 847 of file DalitzEvent.cpp.

                                                                    {
   _rememberPhaseSpace = -9999;
   if(_pat.size() != p4List.size()){
     return false;
   }
   _p.resize(p4List.size());
   for(unsigned int i=0; i<p4List.size(); i++){
     _p[i] = p4List[i];
   }
   return true;
 }

◆ setMomenta() [2/2]

bool DalitzEvent::setMomenta ( const std::vector< TVector3 > & mumAndDgtr_p3 )

protected

Definition at line 859 of file DalitzEvent.cpp.

                                                              {
   _rememberPhaseSpace = -9999;
 
   _p.resize(p3List.size());
   if(_pat.size() != p3List.size()){
     return false;
   }
   for(unsigned int i=0; i< _pat.size(); i++){
     const TVector3& p3 = p3List[i];
     double m = _pat[i].mass();
     _p[i].SetXYZM(p3.X(), p3.Y(), p3.Z(), m);
   }
   return true;
 }

◆ setMothers3Momentum()

void DalitzEvent::setMothers3Momentum ( const TVector3 & mp3 )

virtual

Implements IDalitzEvent.

Definition at line 737 of file DalitzEvent.cpp.

                                                         {
   TLorentzVector newMum4;
   newMum4.SetXYZM(mp3.X(), mp3.Y(), mp3.Z(), _p[0].M());
   
   TVector3 boostToOldMumsRest(- _p[0].BoostVector());
   TVector3 boostToNewMumsLab(newMum4.BoostVector());
   for(unsigned int i=0; i < _p.size(); i++){
     _p[i].Boost(boostToOldMumsRest); // to D0 rest frame
     _p[i].Boost(boostToNewMumsLab); //  to lab frame with new D0 momentum   
   }
   
 }

◆ setP()

void DalitzEvent::setP	(	unsigned int	i,
		const TLorentzVector &	p4
	)

protected

Definition at line 874 of file DalitzEvent.cpp.

                                                               {
   _rememberPhaseSpace = -9999;
   if(i >= _p.size()){
     _p.resize(i+1);
     //    cout << "had to resize" << endl;
   }
   _p[i] = p4;
 }

◆ setPermutationIndex()

void DalitzEvent::setPermutationIndex ( int i )

virtual

Implements IDalitzEvent.

Definition at line 839 of file DalitzEvent.cpp.

                                           {
   if(i >= numPermutations()) return;
   _permutationIndex=i;
 }

◆ setValue() [1/2]

virtual void DalitzEvent::setValue	(	int	i,
		const std::complex< double > &	value,
		long int	configNumber = `1`
	)

inlinevirtual

Implements IDalitzEvent.

Definition at line 201 of file DalitzEvent.h.

                                       {
     //std::cout << "setting: " << i << ") " << value << std::endl;
     _rememberComplexFast.set(i,value, configNumber);
   }

◆ setValue() [2/2]

virtual void DalitzEvent::setValue	(	int	i,
		double	value,
		long int	configNumber = `1`
	)

inlinevirtual

Implements IDalitzEvent.

Definition at line 209 of file DalitzEvent.h.

                                                                      {
     _rememberDoubleFast.set(i,value, configNumber);
   }

◆ setValueInVector()

void DalitzEvent::setValueInVector	(	unsigned int	i,
		double	value
	)

virtual

Implements IDalitzEvent.

Definition at line 415 of file DalitzEvent.cpp.

                                                               {
     if(_vectorOfValues.size() <= i) _vectorOfValues.resize(i+1);
     _vectorOfValues[i]=value;
 }

◆ setWeight()

void DalitzEvent::setWeight ( double w )

virtual

Implements IDalitzEvent.

Definition at line 393 of file DalitzEvent.cpp.

                                    {
   _weight = w;
 }

◆ setWeightInVector()

void DalitzEvent::setWeightInVector	(	unsigned int	i,
		double	weight
	)

virtual

Implements IDalitzEvent.

Definition at line 420 of file DalitzEvent.cpp.

                                                                 {
     if(_vectorOfWeights.size() <= i) _vectorOfWeights.resize(i+1);
     _vectorOfWeights[i]=weight;
 }

◆ shoutAndKill()

bool DalitzEvent::shoutAndKill ( )

protected

Definition at line 787 of file DalitzEvent.cpp.

                               {
   cout << "FATAL ERROR in DalitzEvent:"
        << "\n  > pattern: \"" << _pat << "\""
        << "\n  > 4-vector list size: " << _p.size()
        << endl;
   throw "probably insonsistent pattern in DalitzEvent";
 }

◆ sij() [1/2]

double DalitzEvent::sij ( const std::vector< int > & indices ) const

virtual

Definition at line 709 of file DalitzEvent.cpp.

                                                           {
   bool dbThis=false;
 
   if(1 == indices.size()){
     return p(indices[0]).Mag2();
   }else if(2 == indices.size()){
     return s(indices[0], indices[1]);
   }else{
     if(dbThis){
       TLorentzVector psum(0.0, 0.0, 0.0, 0.0);
       for(unsigned int i=0; i< indices.size(); i++){
     //cout << " p(" << i << ") " << p(indices[i]);
     psum += p(indices[i]);
       }
       cout << "DalitzEvent sij: compare "  << psum.Mag2()
        << " = " << sijMap(indices) << endl;
     }
     return sijMap(indices);
   }
 }

◆ sij() [2/2]

virtual double DalitzEvent::sij ( const MINT::PolymorphVector< int > & indices ) const

inlinevirtual

Implements IDalitzEvent.

Definition at line 179 of file DalitzEvent.h.

                                                                    {
     return sij(indices.theVector());
   }

◆ sijMap()

double DalitzEvent::sijMap ( const std::vector< int > & indices ) const

protected

Definition at line 665 of file DalitzEvent.cpp.

                                                              {
   vector<int> intern_indices(indices);
   for(unsigned int i=0; i < intern_indices.size(); i++){
     intern_indices[i] = currentPermutation()[ indices[i] ];
   }
   return sijMap_intern(intern_indices);
 }

◆ sijMap_intern()

double DalitzEvent::sijMap_intern ( const std::vector< int > & intern_indices ) const

protected

Definition at line 655 of file DalitzEvent.cpp.

                                                                            {
   bool successFlag=false;
   double returnVal = MINT::keyFinder(intern_indices, _sijMap
                      , (double) 0, successFlag);
   if(! successFlag){
     returnVal = evalsij_intern(intern_indices);
     _sijMap[intern_indices] = returnVal;
   }
   return returnVal;
 }

◆ sijMax() [1/2]

double DalitzEvent::sijMax ( const std::vector< int > & indices ) const

Definition at line 733 of file DalitzEvent.cpp.

                                                              {
   return _pat.sijMax(indices);
 }

◆ sijMax() [2/2]

double DalitzEvent::sijMax ( const MINT::PolymorphVector< int > & indices ) const

inlinevirtual

Implements IDalitzEvent.

Definition at line 219 of file DalitzEvent.h.

219 {

220 return sijMax(indices.theVector());}

MINT::PolymorphVector::theVector

std::vector< T > & theVector()

Definition: PolymorphVector.h:22

DalitzEvent::sijMax

double sijMax(const std::vector< int > &indices) const

Definition: DalitzEvent.cpp:733

◆ sijMin() [1/2]

double DalitzEvent::sijMin ( const std::vector< int > & indices ) const

Definition at line 730 of file DalitzEvent.cpp.

                                                              {
   return _pat.sijMin(indices);
 }

◆ sijMin() [2/2]

double DalitzEvent::sijMin ( const MINT::PolymorphVector< int > & indices ) const

inlinevirtual

Implements IDalitzEvent.

Definition at line 215 of file DalitzEvent.h.

215 {

216 return sijMin(indices.theVector());}

DalitzEvent::sijMin

double sijMin(const std::vector< int > &indices) const

Definition: DalitzEvent.cpp:730

MINT::PolymorphVector::theVector

std::vector< T > & theVector()

Definition: PolymorphVector.h:22

◆ singleParticleNtpArraySize()

int DalitzEvent::singleParticleNtpArraySize ( )

staticprotected

Definition at line 1011 of file DalitzEvent.cpp.

                                            {
   return 5;
 }

◆ t()

double DalitzEvent::t	(	unsigned int	i,
		unsigned int	j
	)		const

virtual

Implements IDalitzEvent.

Definition at line 694 of file DalitzEvent.cpp.

                                                          {
   return t_intern(currentPermutation()[i], currentPermutation()[j]);
 }

◆ t_intern()

double DalitzEvent::t_intern	(	unsigned int	i,
		unsigned int	j
	)		const

protected

Definition at line 651 of file DalitzEvent.cpp.

                                                                 {
   if(_t[i][j] < 0) _t[i][j] = (p_intern(i) - p_intern(j)).Mag2();
   return _t[i][j];
 }

Member Data Documentation

◆ _aValue

double DalitzEvent::_aValue

protected

Definition at line 50 of file DalitzEvent.h.

◆ _eventCounter

long int DalitzEvent::_eventCounter =0

staticprotected

Definition at line 37 of file DalitzEvent.h.

◆ _generatorPdfRelativeToPhaseSpace

double DalitzEvent::_generatorPdfRelativeToPhaseSpace

protected

Definition at line 52 of file DalitzEvent.h.

◆ _p

std::vector<TLorentzVector> DalitzEvent::_p

protected

Definition at line 43 of file DalitzEvent.h.

◆ _pat

DalitzEventPattern DalitzEvent::_pat

protected

Definition at line 42 of file DalitzEvent.h.

◆ _perm

Permutator DalitzEvent::_perm

protected

Definition at line 63 of file DalitzEvent.h.

◆ _permutationIndex

int DalitzEvent::_permutationIndex

protected

Definition at line 62 of file DalitzEvent.h.

◆ _rememberComplexFast

RememberAnythingFast<std::complex<double> > DalitzEvent::_rememberComplexFast

protected

Definition at line 47 of file DalitzEvent.h.

◆ _rememberDoubleFast

RememberAnythingFast<double > DalitzEvent::_rememberDoubleFast

protected

Definition at line 48 of file DalitzEvent.h.

◆ _rememberPhaseSpace

double DalitzEvent::_rememberPhaseSpace

mutableprotected

Definition at line 45 of file DalitzEvent.h.

◆ _rememberVectorCounter

long int DalitzEvent::_rememberVectorCounter =0

staticprotected

Definition at line 39 of file DalitzEvent.h.

◆ _s

std::vector< std::vector<double> > DalitzEvent::_s

mutableprotected

Definition at line 57 of file DalitzEvent.h.

◆ _sijMap

std::map<std::vector<int>, double> DalitzEvent::_sijMap

mutableprotected

Definition at line 59 of file DalitzEvent.h.

◆ _t

std::vector< std::vector<double> > DalitzEvent::_t

mutableprotected

Definition at line 58 of file DalitzEvent.h.

◆ _vectorOfValues

std::vector<double> DalitzEvent::_vectorOfValues

protected

Definition at line 54 of file DalitzEvent.h.

◆ _vectorOfWeights

std::vector<double> DalitzEvent::_vectorOfWeights

protected

Definition at line 55 of file DalitzEvent.h.

◆ _weight

double DalitzEvent::_weight

protected

Definition at line 51 of file DalitzEvent.h.

◆ ntpNameChars

const char DalitzEvent::ntpNameChars = { '#', '~', '\0' }

staticprotected

Definition at line 30 of file DalitzEvent.h.

◆ prtNameChars

const char DalitzEvent::prtNameChars = { '+', '-', '\0' }

staticprotected

Definition at line 31 of file DalitzEvent.h.

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

Mint/DalitzEvent.h
src/Mojito/DalitzEvents/DalitzEvent.cpp

Public Member Functions

Static Public Member Functions

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DalitzEvent() [1/12]

◆ DalitzEvent() [2/12]

◆ DalitzEvent() [3/12]

◆ DalitzEvent() [4/12]

◆ DalitzEvent() [5/12]

◆ DalitzEvent() [6/12]

◆ DalitzEvent() [7/12]

◆ DalitzEvent() [8/12]

◆ DalitzEvent() [9/12]

◆ DalitzEvent() [10/12]

◆ DalitzEvent() [11/12]

◆ DalitzEvent() [12/12]

◆ ~DalitzEvent()

Member Function Documentation

◆ assignUniqueRememberNumber()

◆ BDet()

◆ C_conjugateYourself()

◆ clone()

◆ CP_conjugateYourself()

◆ currentPermutation()

◆ Delta4()

◆ evalsij_intern()

◆ eventCounter()

◆ eventPattern()

◆ fillNtupleVarArray()

◆ fromNtuple()

◆ fromParasTree()

◆ fromParasTreeOld()

◆ fromTree()

◆ generateThisToPhaseSpace() [1/3]

◆ generateThisToPhaseSpace() [2/3]

◆ generateThisToPhaseSpace() [3/3]

◆ getAValue()

◆ getGeneratorPdfRelativeToPhaseSpace()

◆ getValueFromVector()

◆ getVectorOfValues() [1/2]

◆ getVectorOfValues() [2/2]

◆ getVectorOfWeights() [1/2]

◆ getVectorOfWeights() [2/2]

◆ getWeight()

◆ getWeightFromVector()

◆ kinematicallyAllowed()

◆ m()

◆ m2()

◆ makeCalculate4BodyProps()

◆ makeNtupleVarnames()

◆ nDgtr()

◆ ntpToPrtName()

◆ ntupleVarArraySize()

◆ numPermutations()

◆ p() [1/2]

◆ p() [2/2]

◆ P_conjugateYourself()

◆ p_dgtr_sum()

◆ p_intern() [1/2]

◆ p_intern() [2/2]

◆ parseNtpEntryName()

◆ permutationIndex()

◆ phaseSpace()

◆ phaseSpace3()

◆ phaseSpace4()

◆ print()

◆ prtToNtpName()

◆ resetST()

◆ retrieveValue() [1/2]

◆ retrieveValue() [2/2]

◆ s()

◆ s_intern()

◆ setAValue()

◆ setGeneratorPdfRelativeToPhaseSpace()

◆ setMomenta() [1/2]

◆ setMomenta() [2/2]