DecayTreeFitter::DecayChain Class Reference

#include <DecayChain.h>

Public Member Functions
	DecayChain ()
	DecayChain (RhoCandidate *bc, const Configuration &config)
	DecayChain (RhoCandidate *bc, const RhoVector3Err &pv, const Configuration &config)
	DecayChain (RhoCandidate *bc, const RhoLorentzVectorErr &lv, const Configuration &config)
	DecayChain (RhoCandidate *bc, const RhoLorentzVectorErr &lv, const RhoVector3Err &pv, const Configuration &config)
virtual	~DecayChain ()
int	dim () const
void	initConstraintList ()
ErrCode	init (FitParams *par)
ErrCode	filter (FitParams *par, bool firstpass=true)
double	chiSquare (const FitParams *par) const
ParticleBase *	mother ()
const ParticleBase *	cand ()
const ParticleBase *	mother () const
const ParticleBase *	locate (RhoCandidate *bc) const
int	index (RhoCandidate *bc) const
int	posIndex (RhoCandidate *bc) const
int	momIndex (RhoCandidate *bc) const
int	lenIndex (RhoCandidate *bc) const
void	setOwner (bool b)
void	printConstraints (std::ostream &os=std::cout) const
void	setMassConstraint (RhoCandidate *bc, bool add=true)
void	setMassConstraint (RhoCandidate *bc, double mass)
ChiSquare	chiSquare (RhoCandidate *bc, const FitParams *fitpars) const

Private Types
typedef std::map< RhoCandidate *, const ParticleBase * >	ParticleMap

Private Attributes
int	m_dim
ParticleBase *	m_mother
const ParticleBase *	m_cand
ParticleBase::constraintlist	m_constraintlist
std::vector< Constraint * >	m_mergedconstraintlist
MergedConstraint	m_mergedconstraint
ParticleMap	m_particleMap
bool	m_isOwner

Detailed Description

Definition at line 26 of file DecayChain.h.

Member Typedef Documentation

typedef std::map<RhoCandidate*,const ParticleBase*> DecayTreeFitter::DecayChain::ParticleMap

private

Definition at line 72 of file DecayChain.h.

Constructor & Destructor Documentation

DecayTreeFitter::DecayChain::DecayChain ( )

inline

Definition at line 29 of file DecayChain.h.

: m_mother(0) {}

DecayTreeFitter::DecayChain::DecayChain	(	RhoCandidate *	bc,
		const Configuration &	config
	)

Definition at line 30 of file DecayChain.cxx.

References locate(), m_cand, m_dim, m_mother, and DecayTreeFitter::ParticleBase::updateIndex().

: m_dim(0),m_mother(0),m_isOwner(true)
{
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate* bc="<<bc<<")"<<std::endl;
  m_mother = (ParticleBase*) new InteractionPoint(bc,config) ;
  //m_mother = ParticleBase::createParticle(bc,0,config) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*) - updateIndex() now"<<std::endl;
  m_mother->updateIndex(m_dim) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*) - locate() now"<<std::endl;
  m_cand   = locate(bc) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*) - done"<<std::endl;
}

DecayTreeFitter::DecayChain::DecayChain	(	RhoCandidate *	bc,
		const RhoVector3Err &	pv,
		const Configuration &	config
	)

Definition at line 43 of file DecayChain.cxx.

References locate(), m_cand, m_dim, m_mother, and DecayTreeFitter::ParticleBase::updateIndex().

: m_dim(0),m_mother(0),m_isOwner(true)
{
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate* bc="<<bc<<",RhoVector3Err v("<<pv.x()<<","<<pv.y()<<","<<pv.z()<<"))"<<std::endl;
  m_mother = (ParticleBase*) new InteractionPoint(pv,bc,config) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - updateIndex() now"<<std::endl;
  m_mother->updateIndex(m_dim) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - locate() now"<<std::endl;
  m_cand   = locate(bc) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - done"<<std::endl;
}

DecayTreeFitter::DecayChain::DecayChain	(	RhoCandidate *	bc,
		const RhoLorentzVectorErr &	lv,
		const Configuration &	config
	)

Definition at line 55 of file DecayChain.cxx.

References locate(), m_cand, m_dim, m_mother, and DecayTreeFitter::ParticleBase::updateIndex().

: m_dim(0),m_mother(0),m_isOwner(true)
{
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate* bc="<<bc<<",RhoLorentzVectorErr p("<<lv.X()<<","<<lv.Y()<<","<<lv.Z()<<","<<lv.E()<<"))"<<std::endl;
  m_mother = (ParticleBase*) new InteractionPoint(lv,bc,config) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - updateIndex() now"<<std::endl;
  m_mother->updateIndex(m_dim) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - locate() now"<<std::endl;
  m_cand   = locate(bc) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - done"<<std::endl;
}

DecayTreeFitter::DecayChain::DecayChain	(	RhoCandidate *	bc,
		const RhoLorentzVectorErr &	lv,
		const RhoVector3Err &	pv,
		const Configuration &	config
	)

Definition at line 67 of file DecayChain.cxx.

References locate(), m_cand, m_dim, m_mother, and DecayTreeFitter::ParticleBase::updateIndex().

: m_dim(0),m_mother(0),m_isOwner(true)
{
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate* bc="<<bc<<",RhoLorentzVectorErr p("<<lv.X()<<","<<lv.Y()<<","<<lv.Z()<<","<<lv.E()<<")"<<",RhoVector3Err v("<<pv.x()<<","<<pv.y()<<","<<pv.z()<<"))"<<std::endl;
  m_mother = (ParticleBase*) new InteractionPoint(lv,pv,bc,config) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - updateIndex() now"<<std::endl;
  m_mother->updateIndex(m_dim) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - locate() now"<<std::endl;
  m_cand   = locate(bc) ;
  //std::cout<<"DecayTreeFitter::DecayChain::DecayChain(RhoCandidate*,RhoVector3Err) - done"<<std::endl;
}

DecayTreeFitter::DecayChain::~DecayChain ( )

virtual

Definition at line 79 of file DecayChain.cxx.

{
  if(m_mother && m_isOwner) delete m_mother ;
}

Member Function Documentation

const ParticleBase* DecayTreeFitter::DecayChain::cand ( )

inline

Definition at line 46 of file DecayChain.h.

References m_cand.

{ return m_cand ; }

double DecayTreeFitter::DecayChain::chiSquare

(

const FitParams *

par

)

const

Definition at line 195 of file DecayChain.cxx.

{
  return m_mother->chiSquareD(par) ;
}

ChiSquare DecayTreeFitter::DecayChain::chiSquare	(	RhoCandidate *	bc,
		const FitParams *	fitpars
	)		const

Definition at line 298 of file DecayChain.cxx.

References DecayTreeFitter::ParticleBase::chiSquare(), DecayTreeFitter::ParticleBase::mother(), and DecayTreeFitter::ParticleBase::posIndex().

                                                                                      {
  const ParticleBase* base = locate(bc) ;
  ChiSquare chisq ;
  if( base ) {
    chisq = base->chiSquare(fitpars) ;
    // uhh: need to subtract the lifetime dof if this particle has a
    // mother with a different vertex
    if( base->mother() && base->posIndex()>=0 &&
       base->posIndex() != base->mother()->posIndex() )
      chisq += ChiSquare(0,1) ;
  }
  return chisq ;
}

int DecayTreeFitter::DecayChain::dim ( ) const

inline

Definition at line 38 of file DecayChain.h.

References m_dim.

{ return m_dim ; }

ErrCode DecayTreeFitter::DecayChain::filter	(	FitParams *	par,
		bool	firstpass = true
	)

Definition at line 145 of file DecayChain.cxx.

References DecayTreeFitter::ErrCode::failure(), DecayTreeFitter::ErrCode::Print(), DecayTreeFitter::FitParams::resetCov(), status, DecayTreeFitter::FitParams::testCov(), and vtxverbose.

{
  if(vtxverbose>=4){std::cout<<"DecayChain::filter(): - ok, let's do this! Btw. we're "<<( firstpass ? "" : "not ")<<"on first pass."<<std::endl;}
  ErrCode status ;
  if(m_constraintlist.empty()){ 
    if(vtxverbose>=4){std::cout<<"DecayChain::filter(): - init constrains"<<std::endl;}
    initConstraintList() ;
  }
  
  par->resetCov(1000) ; //FIXME: this scales the cov with a factor 1000. WHY???
  //par->resetCov(100) ; //FIXME: this scales the cov with a factor 1000. WHY???
  if(firstpass || !par->testCov()){
    status |= m_mother->initCov(par) ;
    if(vtxverbose>=3){std::cout<<"DecayChain::filter(): - initcov done"<<std::endl;}
  }

  if( vtxverbose>=4 || (vtxverbose>=2&&firstpass) ) {
    std::cout << "DecayTreeFitter::DecayChain::filter, after initialization: "  << std::endl ;
    m_mother->print(par) ;
  }

  FitParams* reference = new FitParams(*par) ;

  if(m_mergedconstraintlist.empty() ) {
    for( ParticleBase::constraintlist::const_iterator it = m_constraintlist.begin() ;
        it != m_constraintlist.end(); ++it) {
      status |= it->filter(par,reference) ;
      if( vtxverbose>=2 && status.failure() ) {
        std::cout << "DecayChain::filter(): status is failure after parsing constraint: " ;
        it->print() ;
      }
    }
  } else {
    for( std::vector<Constraint*>::const_iterator it = m_mergedconstraintlist.begin() ;
        it != m_mergedconstraintlist.end(); ++it) {
      status |= (*it)->filter(par,reference) ;
      if( vtxverbose>=2 && status.failure() ) {
        std::cout << "DecayChain::filter(): status is failure after parsing constraint: " ;
        (*it)->print() ;
      }
    }
  }


  if(vtxverbose>=3) {std::cout << "DecayChain::filter: status = "; status.Print(std::cout); std::cout<< std::endl ;}
  if(reference) delete reference;
  return status ;
}

int DecayTreeFitter::DecayChain::index

(

RhoCandidate *

bc

)

const

Definition at line 272 of file DecayChain.cxx.

References DecayTreeFitter::ParticleBase::index().

                                                       {
  int rc = -1 ;
  const ParticleBase* base = locate(bc) ;
  if(base) rc = base->index() ;
  return rc ;
}

ErrCode DecayTreeFitter::DecayChain::init

(

FitParams *

par

)

FIXME Here we have some bad behavior: We're using automatic copy constructor, implicit copied objects etc. How to do cleanly? Pointers directly with delete/new?

Definition at line 114 of file DecayChain.cxx.

References DecayTreeFitter::FitParams::dim(), DecayTreeFitter::ErrCode::Print(), DecayTreeFitter::FitParams::reset(), DecayTreeFitter::FitParams::resetCov(), DecayTreeFitter::FitParams::resetPar(), status, and vtxverbose.

{
  if(vtxverbose>5){std::cout<<"DecayChain::init(): - start"<<std::endl;}
  ErrCode status ;
  if( m_dim==0 ){
  if(vtxverbose>5){std::cout<<"DecayChain::init(): - update index"<<std::endl;}
     m_mother->updateIndex(m_dim) ;
  }
  //if( par.dim() != m_dim ) par = FitParams(m_dim) ;
  if( par->dim() != m_dim ){
    if(vtxverbose>5){std::cout<<"DecayChain::init(): - reset 1"<<std::endl;}
     par->reset(m_dim);
  }  else {
    // set everything to 0
    if(vtxverbose>5){std::cout<<"DecayChain::init(): - reset 2"<<std::endl;}
    par->resetPar() ;
    par->resetCov() ;
  }
  // let the mother do it
  if(vtxverbose>5){std::cout<<"DecayChain::init(): - initpar"<<std::endl;}
  status |= m_mother->initPar1(par) ;
  if(vtxverbose>5){std::cout<<"DecayChain::init(): - initcov"<<std::endl;}
  status |= m_mother->initCov(par) ;

  if(vtxverbose>=2){ std::cout << "status after init: "; status.Print(std::cout); std::cout<< std::endl ;}

  return status ;
}

void DecayTreeFitter::DecayChain::initConstraintList

(

)

Definition at line 85 of file DecayChain.cxx.

References vtxverbose.

{
  if( m_dim==0 ) m_mother->updateIndex(m_dim) ;
  m_constraintlist.clear() ;
  m_mother->addToConstraintList(m_constraintlist,0) ;

  // the order of the constraints is a rather delicate thing
  std::sort(m_constraintlist.begin(),m_constraintlist.end()) ;

  // merge all non-lineair constraints
  m_mergedconstraintlist.clear() ;
  if(false) {
    m_mergedconstraint = MergedConstraint() ;
    for( ParticleBase::constraintlist::iterator it =  m_constraintlist.begin() ;
        it != m_constraintlist.end(); ++it) {
      if( it->isLineair() ) m_mergedconstraintlist.push_back(&(*it)) ;
      else  m_mergedconstraint.push_back(&(*it)) ;
    }
    if( m_mergedconstraint.dim()>0 )
      m_mergedconstraintlist.push_back(&m_mergedconstraint) ;
  }
  if(vtxverbose>=2) {
    std::cout << "DecayTreeFitter::DecayChain::initConstraintList(): " << m_constraintlist.size() << " "
    << m_mergedconstraintlist.size()<< std::endl ;
    printConstraints() ;
  }
}

int DecayTreeFitter::DecayChain::lenIndex

(

RhoCandidate *

bc

)

const

Definition at line 292 of file DecayChain.cxx.

References DecayTreeFitter::ParticleBase::lenIndex().

                                                          {
  const ParticleBase* base = locate(bc) ;
  return base ? base->lenIndex() : 0 ;
}

const ParticleBase * DecayTreeFitter::DecayChain::locate

(

RhoCandidate *

bc

)

const

Definition at line 201 of file DecayChain.cxx.

References DecayTreeFitter::ParticleBase::locate(), and DecayTreeFitter::ParticleBase::particle().

Referenced by DecayChain().

                                                        {
  const ParticleBase* rc(0) ;
  // return m_mother->locate(bc) ;
  ParticleMap::const_iterator it = m_particleMap.find(bc) ;
  if( it == m_particleMap.end() ) {
    rc = m_mother->locate(bc) ;
    // used to add every candidate here, but something goes wrong
    // somewhere. now only cache pointers that we internally reference.
    if(rc) m_particleMap[rc->particle()] = rc ;
  } else {
    rc = it->second ;
  }
  return rc ;
}

int DecayTreeFitter::DecayChain::momIndex

(

RhoCandidate *

bc

)

const

Definition at line 286 of file DecayChain.cxx.

References DecayTreeFitter::ParticleBase::momIndex().

                                                          {
  const ParticleBase* base = locate(bc) ;
  return base ? base->momIndex() : -1 ;
}

ParticleBase* DecayTreeFitter::DecayChain::mother ( )

inline

Definition at line 45 of file DecayChain.h.

References m_mother.

{ return m_mother ; }

const ParticleBase* DecayTreeFitter::DecayChain::mother ( ) const

inline

Definition at line 47 of file DecayChain.h.

References m_mother.

{ return m_mother ; }

int DecayTreeFitter::DecayChain::posIndex

(

RhoCandidate *

bc

)

const

Definition at line 280 of file DecayChain.cxx.

References DecayTreeFitter::ParticleBase::posIndex().

                                                          {
  const ParticleBase* base = locate(bc) ;
  return base ? base->posIndex() : -1 ;
}

void DecayTreeFitter::DecayChain::printConstraints

(

std::ostream &

os = std::cout

)

const

Definition at line 21 of file DecayChain.cxx.

{
  os << "Constraints in decay tree: " << m_constraintlist.size() << std::endl ;
  for( ParticleBase::constraintlist::const_iterator
      it = m_constraintlist.begin() ;
      it != m_constraintlist.end(); ++it)
    it->print(os) ;
}

void DecayTreeFitter::DecayChain::setMassConstraint	(	RhoCandidate *	bc,
		bool	add = true
	)

Definition at line 217 of file DecayChain.cxx.

References DecayTreeFitter::ParticleBase::setMassConstraint().

{
  ParticleBase* part = const_cast<ParticleBase*>(locate(bc)) ;
  if(part && part->setMassConstraint(add)) {
    m_dim = 0 ;
    m_constraintlist.clear() ;
  }
}

void DecayTreeFitter::DecayChain::setMassConstraint	(	RhoCandidate *	bc,
		double	mass
	)

Definition at line 227 of file DecayChain.cxx.

References DecayTreeFitter::ParticleBase::setMassConstraint().

{
  ParticleBase* part = const_cast<ParticleBase*>(locate(bc)) ;
  if(part) {
    part->setMassConstraint(mass) ;
    m_dim = 0 ;
    m_constraintlist.clear() ;
  }
}

void DecayTreeFitter::DecayChain::setOwner ( bool  b )

inline

Definition at line 55 of file DecayChain.h.

References b, and m_isOwner.

{ m_isOwner=b ;}

Member Data Documentation

const ParticleBase* DecayTreeFitter::DecayChain::m_cand

private

Definition at line 68 of file DecayChain.h.

Referenced by cand(), and DecayChain().

ParticleBase::constraintlist DecayTreeFitter::DecayChain::m_constraintlist

private

Definition at line 69 of file DecayChain.h.

int DecayTreeFitter::DecayChain::m_dim

private

Definition at line 66 of file DecayChain.h.

Referenced by DecayChain(), and dim().

bool DecayTreeFitter::DecayChain::m_isOwner

private

Definition at line 74 of file DecayChain.h.

Referenced by setOwner().

MergedConstraint DecayTreeFitter::DecayChain::m_mergedconstraint

private

Definition at line 71 of file DecayChain.h.

std::vector<Constraint*> DecayTreeFitter::DecayChain::m_mergedconstraintlist

private

Definition at line 70 of file DecayChain.h.

ParticleBase* DecayTreeFitter::DecayChain::m_mother

private

Definition at line 67 of file DecayChain.h.

Referenced by DecayChain(), and mother().

ParticleMap DecayTreeFitter::DecayChain::m_particleMap

mutableprivate

Definition at line 73 of file DecayChain.h.

---

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

• DecayChain.h
• DecayChain.cxx