InternalParticle.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

// ******************************************************
// DecayTreeFitter Package
// We thank the original author Wouter Hulsbergen
// (BaBar, LHCb) for providing the sources.
// http://arxiv.org/abs/physics/0503191v1 (2005)
// Adaptation & Development for PANDA: Ralf Kliemt (2015)
// Adjustment for fully neutral reactions: Jean Noel (2023)
// **********************************************************************************
#ifndef INTERNALPARTICLE_H
#define INTERNALPARTICLE_H 1

#include "ParticleBase.h"
#include <vector>
#include "Rtypes.h"

extern bool fitvertex;
extern bool addfreeZ;

namespace DecayTreeFitter {

class InternalParticle : public ParticleBase {
 public:
  InternalParticle(RhoCandidate *bc, const ParticleBase *mother, const Configuration &config);

  virtual int dim() const { return fitvertex ? (mother() ? 8 : 7) : (addfreeZ ? 5 : 4); }

  virtual ErrCode initPar1(FitParams *);
  virtual ErrCode initPar2(FitParams *);
  virtual int type() const { return kInternalParticle; }

  // parameter definition
  virtual int posIndex() const { return fitvertex ? index() : (addfreeZ ? index() : -1); }
  virtual int lenIndex() const { return fitvertex ? (mother() ? index() + 3 : -1) : -1; }
  virtual int momIndex() const { return fitvertex ? (mother() ? index() + 4 : index() + 3) : (addfreeZ ? index() + 1: index()); }
  virtual bool hasEnergy() const { return true; }
  virtual bool hasPosition() const { return fitvertex ? true : (addfreeZ ? true : false); }         
  virtual std::string parname(int index) const;

  // constraints
  ErrCode projectKineConstraint(const FitParams *, Projection &) const;
  ErrCode projectLifeTimeConstraint(const FitParams *, Projection &) const;
  ErrCode projectConversionConstraint(const FitParams *, Projection &p) const;
  ErrCode projectMassConstraintTwoBody(const FitParams *fitparams, Projection &p) const;
  virtual ErrCode projectConstraint(Constraint::Type type, const FitParams *fitparams, Projection &p) const;

  // some of that other stuff
  virtual void addToConstraintList(constraintlist &alist, int depth) const;

  // bool swapMotherDaughter(FitParams* fitparams, const ParticleBase* newmother) ;

 protected:
  ErrCode initMom(FitParams *fitparams) const;

 private:
  bool m_lifetimeconstraint;
  bool m_isconversion;
  ClassDef(InternalParticle, 1)
};

} // namespace DecayTreeFitter

#endif