ISFParticle.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#include "ISF_Event/ISFParticle.h"
 
#include <cmath> // fabs
 
ISF::ISFParticle::ISFParticle(
                              const Amg::Vector3D& pos,
                              const Amg::Vector3D& mom,
                              double mass,
                              double charge,
                              int pdgCode,
                              int status,
                              double time,
                              const ISFParticle &parent,
                              int id,
                             int barcode,
                              TruthBinding* truth,
                              const HepMcParticleLink* partLink):
  m_position(pos),
  m_momentum(mom),
  m_mass(mass),
  m_charge(charge),
  m_pdgCode(pdgCode),
  m_status(status),
  m_tstamp(time),
  m_history(parent.history()),
  m_barcode(barcode), // TODO drop this soon
  m_uid(id),
  m_truth(truth),
  m_order(ISF::DefaultParticleOrder),
  m_userInfo(nullptr),
  m_partLink(partLink)
{
 
}
 
ISF::ISFParticle::ISFParticle(
                              const HepGeom::Point3D<double>& pos,
                              const HepGeom::Vector3D<double>& mom,
                              double mass,
                              double charge,
                              int pdgCode,
                              int status,
                              double time,
                              const ISFParticle &parent,
                              int id,
                              int barcode,
                              TruthBinding* truth,
                              const HepMcParticleLink* partLink):
  m_position( pos.x(), pos.y(), pos.z()),
  m_momentum( mom.x(), mom.y(), mom.z()),
  m_mass(mass),
  m_charge(charge),
  m_pdgCode(pdgCode),
  m_status(status),
  m_tstamp(time),
  m_history(parent.history()),
  m_barcode(barcode), // TODO drop this soon
  m_uid(id),
  m_truth(truth),
  m_order(ISF::DefaultParticleOrder),
  m_userInfo(nullptr),
  m_partLink(partLink)
{
 
}
 
ISF::ISFParticle::ISFParticle(
                              const Amg::Vector3D& pos,
                              const Amg::Vector3D& mom,
                              double mass,
                              double charge,
                              int pdgCode,
                              int status,
                              double time,
                              const DetRegionSvcIDPair &origin,
                              int id,
                              int barcode,
                              TruthBinding* truth,
                              const HepMcParticleLink* partLink):
  m_position(pos),
  m_momentum(mom),
  m_mass(mass),
  m_charge(charge),
  m_pdgCode(pdgCode),
  m_status(status),
  m_tstamp(time),
  m_history(1, origin),
  m_barcode(barcode), // TODO drop this soon
  m_uid(id),
  m_truth(truth),
  m_order(ISF::DefaultParticleOrder),
  m_userInfo(nullptr),
  m_partLink(partLink)
{
 
}
 
ISF::ISFParticle::ISFParticle(const ISFParticle& isfp):
  m_position(isfp.position()),
  m_momentum(isfp.momentum()),
  m_mass(isfp.mass()),
  m_charge(isfp.charge()),
  m_pdgCode(isfp.pdgCode()),
  m_status(isfp.status()),
  m_tstamp(isfp.timeStamp()),
  m_history(isfp.history()),
  m_barcode(isfp.barcode()), // TODO drop this soon
  m_uid(isfp.id()),
  m_truth(nullptr),
  m_order(ISF::DefaultParticleOrder),
  m_userInfo(nullptr)
{
  const TruthBinding *truth = isfp.getTruthBinding();
  if (truth) {
    m_truth = new TruthBinding(*truth);
  }
  m_partLink = ((isfp.getParticleLink()) ? new HepMcParticleLink(*isfp.getParticleLink()) : nullptr);
}
 
ISF::ISFParticle::ISFParticle(ISFParticle&& isfp):
  m_position(isfp.position()),
  m_momentum(isfp.momentum()),
  m_mass(isfp.mass()),
  m_charge(isfp.charge()),
  m_pdgCode(isfp.pdgCode()),
  m_status(isfp.status()),
  m_tstamp(isfp.timeStamp()),
  m_history(isfp.history()),
  m_barcode(isfp.barcode()), // TODO drop this soon
  m_uid(isfp.id()),
  m_truth(isfp.getTruthBinding()),
  m_order(isfp.getOrder()),
  m_userInfo(isfp.getUserInformation()),
  m_partLink(isfp.getParticleLink())
{
}
 
ISF::ISFParticle::~ISFParticle() {
  delete m_truth;
  delete m_userInfo;
  if (m_partLink) delete m_partLink;
}
 
// cppcheck-suppress operatorEqVarError; m_order not copied
ISF::ISFParticle& ISF::ISFParticle::operator=(const ISF::ISFParticle& rhs)
{
 
  if (this != &rhs) {
    m_position     = rhs.position();
    m_momentum     = rhs.momentum();
    m_mass         = rhs.mass();
    m_charge       = rhs.charge();
    m_pdgCode      = rhs.pdgCode();
    m_status       = rhs.status();
    m_tstamp       = rhs.timeStamp();
    m_history      = rhs.history();
    m_barcode      = rhs.barcode(); // TODO drop this soon
    m_uid           = rhs.id();
 
    delete m_truth;
    m_truth = nullptr;
 
    const TruthBinding *rhsTruth = rhs.getTruthBinding();
    if (rhsTruth) {
      m_truth = new TruthBinding(*rhsTruth);
    }
    m_partLink = ((rhs.getParticleLink()) ? new HepMcParticleLink(*rhs.getParticleLink()): nullptr);
 
    delete m_userInfo;
    m_userInfo = nullptr;
  }
 
  return *this;
}
 
// cppcheck-suppress operatorEqVarError; m_order not copied
ISF::ISFParticle& ISF::ISFParticle::operator=(ISF::ISFParticle&& rhs)
{
  m_position     = rhs.position();
  m_momentum     = rhs.momentum();
  m_mass         = rhs.mass();
  m_charge       = rhs.charge();
  m_pdgCode      = rhs.pdgCode();
  m_status       = rhs.status();
  m_tstamp       = rhs.timeStamp();
  m_history      = rhs.history();
  m_barcode      = rhs.barcode(); // TODO drop this soon
  m_uid           = rhs.id();
  delete m_truth;
  m_truth        = rhs.getTruthBinding();
  delete m_userInfo;
  m_userInfo     = rhs.getUserInformation();
  if (m_partLink) delete m_partLink;
  m_partLink     = rhs.getParticleLink();
 
  return *this;
}
 
bool ISF::ISFParticle::operator==(const ISF::ISFParticle& rhs) const
{
  return isEqual(rhs);
}
 
bool ISF::ISFParticle::isEqual(const ISF::ISFParticle& rhs) const
{
  double epsilon = 1e-6;
 
  bool pass = true;
  pass &= m_position == rhs.position();
  pass &= m_momentum == rhs.momentum();
  pass &= std::fabs(m_mass-rhs.mass()) < epsilon;
  pass &= std::fabs(m_charge-rhs.charge()) < epsilon;
  pass &= m_pdgCode == rhs.pdgCode();
  pass &= m_status == rhs.status();
  pass &= std::fabs(m_tstamp-rhs.timeStamp()) < epsilon;
  pass &= m_history == rhs.history();
  pass &= m_barcode == rhs.barcode(); // TODO drop this soon
  pass &= m_uid == rhs.id();
  {
    const auto rhsTruthPtr = rhs.getTruthBinding();
    if (m_truth && rhsTruthPtr) {
      pass &= *m_truth == *rhsTruthPtr;
    } else {
      pass &= m_truth == rhsTruthPtr; // must be both nullptr to pass
    }
  }
 
  {
    const auto rhsUserInfoPtr = rhs.getUserInformation();
    if (m_userInfo && rhsUserInfoPtr) {
      pass &= *m_userInfo == *rhsUserInfoPtr;
    } else {
      pass &= m_userInfo == rhsUserInfoPtr; // must be both nullptr to pass
    }
  }
 
  {
    const auto rhsPartLinkPtr = rhs.getParticleLink();
    if(m_partLink && rhsPartLinkPtr) {
      pass &= *m_partLink == *rhsPartLinkPtr;
    } else {
      pass &= m_partLink == rhsPartLinkPtr; // must be both nullptr to pass
    }
  }
  return pass;
}
 
bool ISF::ISFParticle::isIdent(const ISF::ISFParticle& rhs) const
{
  bool pass = true;
  pass &= m_position == rhs.position();
  pass &= m_momentum == rhs.momentum();
  pass &= m_mass == rhs.mass();
  pass &= m_charge == rhs.charge();
  pass &= m_pdgCode == rhs.pdgCode();
  pass &= m_status == rhs.status();
  pass &= m_tstamp == rhs.timeStamp();
  pass &= m_history == rhs.history();
  pass &= m_barcode == rhs.barcode(); // TODO drop this soon
  pass &= m_uid == rhs.id();
  pass &= m_truth && rhs.getTruthBinding();
  pass &= m_userInfo == rhs.getUserInformation();
  pass &= m_partLink == rhs.getParticleLink();
  return pass;
}