TrackHandle_TruthTrack::Imp Class Reference

Collaboration diagram for TrackHandle_TruthTrack::Imp:

Public Member Functions
	Imp (TrackHandle_TruthTrack *tc, const SimBarCode &sbc, const SimHitList &shl, HepMC::ConstGenParticlePtr p)
void	ensureInitAscObjs ()
void	ensureInitTrkTracks ()
void	createTrack (Trk::TrackStates *trackStateOnSurfaces)

Static Public Member Functions
static Trk::Perigee *	createTrkPerigeeFromProdVertex (HepMC::ConstGenParticlePtr p, const double &charge)
static Trk::TrackParameters *	createTrkParamFromDecayVertex (HepMC::ConstGenParticlePtr p, const double &charge)
static Trk::TrackStateOnSurface *	createTSOS (Trk::TrackParameters *pars)
static void	addPars (Trk::TrackStates *dv, Trk::TrackParameters *pars)

Public Attributes
TrackHandle_TruthTrack *	theclass
SimBarCode	simBarCode
SimHitList	simHitList
HepMC::ConstGenParticlePtr	genParticle
bool	ascObjVis
std::vector< AscObj_TruthPoint * > *	ascObjs
const Trk::Track *	trkTrack

Detailed Description

Definition at line 28 of file TrackHandle_TruthTrack.cxx.

Constructor & Destructor Documentation

Imp()

TrackHandle_TruthTrack::Imp::Imp ( TrackHandle_TruthTrack * tc, const SimBarCode & sbc, const SimHitList & shl, HepMC::ConstGenParticlePtr p )

inline

Definition at line 30 of file TrackHandle_TruthTrack.cxx.

    : theclass(tc),
      simBarCode(sbc),
      simHitList(shl),
      genParticle(p),
      ascObjVis(false),
      ascObjs(nullptr),
      trkTrack(nullptr) {}

Member Function Documentation

addPars()

void TrackHandle_TruthTrack::Imp::addPars ( Trk::TrackStates * dv, Trk::TrackParameters * pars )

inlinestatic

Definition at line 92 of file TrackHandle_TruthTrack.cxx.

  {
    if (!pars)
      return;
    Trk::TrackStateOnSurface * tsos = createTSOS(pars);
    if (tsos)
      dv->push_back(tsos);
  }

createTrack()

void TrackHandle_TruthTrack::Imp::createTrack ( Trk::TrackStates * trackStateOnSurfaces )

inline

Definition at line 101 of file TrackHandle_TruthTrack.cxx.

  {
    if (!trackStateOnSurfaces) {
      VP1Msg::messageDebug("TrackHandle_TruthTrack WARNING: Could not create track due to null TSOS vector");
      return;
    }
    if (trackStateOnSurfaces->empty()) {
      VP1Msg::messageDebug("TrackHandle_TruthTrack WARNING: Could not create track due to empty TSOS vector");
      delete trackStateOnSurfaces;
      return;
    }
    if (trkTrack) {
      VP1Msg::messageDebug("TrackHandle_TruthTrack ERROR: Already create trkTrack previously!");
      delete trackStateOnSurfaces;
      return;
    }
 
    Trk::TrackInfo ti(Trk::TrackInfo::Unknown,theclass->extrapolationParticleHypothesis());
    std::unique_ptr<Trk::TrackStates> sink(trackStateOnSurfaces);
    trkTrack = new Trk::Track(ti,
                              std::move(sink),
                              nullptr /*fitquality*/);
 
  }

createTrkParamFromDecayVertex()

Trk::TrackParameters * TrackHandle_TruthTrack::Imp::createTrkParamFromDecayVertex ( HepMC::ConstGenParticlePtr p, const double & charge )

inlinestatic

Definition at line 65 of file TrackHandle_TruthTrack.cxx.

  {
    if (!p)
      return nullptr;//Fixme: message!
    HepMC::ConstGenVertexPtr v = p->end_vertex();
    if (!v)
      return nullptr;//Fixme: message!
    Amg::Vector3D mom(p->momentum().px(),p->momentum().py(),p->momentum().pz());
//     double absmom(mom.mag());
//     if (absmom<=0)
//       return 0;//Fixme: message!
    Amg::Vector3D pos(v->position().x(),v->position().y(),v->position().z());
 
    Amg::Translation3D amgtranslation(pos.x(),pos.y(),pos.z());
    Amg::Transform3D amgTransf(amgtranslation * Amg::RotationMatrix3D::Identity());
 
    return new Trk::AtaPlane(pos,mom,charge, *(new Trk::PlaneSurface(amgTransf)));
  }

createTrkPerigeeFromProdVertex()

Trk::Perigee * TrackHandle_TruthTrack::Imp::createTrkPerigeeFromProdVertex ( HepMC::ConstGenParticlePtr p, const double & charge )

inlinestatic

Definition at line 50 of file TrackHandle_TruthTrack.cxx.

  {
    if (!p)
      return nullptr;//Fixme: message!
    HepMC::ConstGenVertexPtr v = p->production_vertex();
    if (!v)
      return nullptr;//Fixme: message!
    Amg::Vector3D mom(p->momentum().px(),p->momentum().py(),p->momentum().pz());
    double absmom(mom.mag());
    if (absmom<=0)
      return nullptr;//Fixme: message!
    Amg::Vector3D pos(v->position().x(),v->position().y(),v->position().z());
    return new Trk::Perigee(0.,0.,mom.phi(), mom.theta(), charge/absmom, pos);
   }

createTSOS()

Trk::TrackStateOnSurface * TrackHandle_TruthTrack::Imp::createTSOS ( Trk::TrackParameters * pars )

inlinestatic

Definition at line 84 of file TrackHandle_TruthTrack.cxx.

  {
    return pars ? new Trk::TrackStateOnSurface(
                    nullptr,
                    std::unique_ptr<Trk::TrackParameters>(pars),
                    nullptr)
                : nullptr;
  }

ensureInitAscObjs()

void TrackHandle_TruthTrack::Imp::ensureInitAscObjs

(

)

Definition at line 307 of file TrackHandle_TruthTrack.cxx.

{
  if (ascObjs)
    return;
  ascObjs = new std::vector<AscObj_TruthPoint*>;
  HepMC::ConstGenVertexPtr vprod{nullptr};
  HepMC::ConstGenVertexPtr vend{nullptr};
  if (genParticle) {
   vprod=genParticle->production_vertex();
   vend=genParticle->end_vertex();
  }
  ascObjs->reserve((vprod?1:0)+(vend?1:simHitList.size()));
  if (vprod)
    ascObjs->push_back(new AscObj_TruthPoint(theclass,vprod,genParticle));
 
  if (vend) {
    ascObjs->push_back(new AscObj_TruthPoint(theclass,vend,genParticle));
  } else {
    SimHitList::const_iterator it, itE(simHitList.end());
    for ( it = simHitList.begin(); it != itE; ++it )
      ascObjs->push_back(new AscObj_TruthPoint(theclass,it->second));
  }
  std::vector<AscObj_TruthPoint*>::iterator it, itE(ascObjs->end());
  for (it=ascObjs->begin();it!=itE;++it)
    theclass->registerAssocObject(*it);
  for (it=ascObjs->begin();it!=itE;++it)
    (*it)->setVisible(ascObjVis);
}

ensureInitTrkTracks()

void TrackHandle_TruthTrack::Imp::ensureInitTrkTracks

(

)

Definition at line 157 of file TrackHandle_TruthTrack.cxx.

{
  if (trkTrack)
    return;
 
  //The GenParticle part is used if it is available with a production
  //vertex. The sim. hits are used if present and the genparticle does
  //not have an end vertex:
  bool useGenParticle = genParticle && genParticle->production_vertex();
  bool decayedGenParticle = useGenParticle && genParticle->end_vertex();
  bool useSimHits = !decayedGenParticle && !simHitList.empty();
 
  if (!useGenParticle&&!useSimHits) {
    VP1Msg::message("TrackHandle_TruthTrack ERROR: Track has neither a genparticle or sim. hits!!");
    return;
  }
 
  if (!theclass->hasCharge()) {
    VP1Msg::message("TrackHandle_TruthTrack ERROR: Could not determine particle charge (pdg="
            +QString::number(theclass->pdgCode())+").");//Fixme: I guess we could show non-extrapolated version?
    return;
  }
  const double charge = theclass->charge();
 
  Trk::TrackStates* trackStateOnSurfaces = new Trk::TrackStates;
 
  if (useGenParticle) {
    addPars(trackStateOnSurfaces,createTrkPerigeeFromProdVertex(genParticle,charge));
 
    if (decayedGenParticle) {
      addPars(trackStateOnSurfaces,createTrkParamFromDecayVertex(genParticle,charge));
      createTrack(trackStateOnSurfaces);
      return;
    }
  }
 
  if (useSimHits) {
    //Add parameters from simhits (yes, if !useGenParticle, we get no perigee).
    SimHitList::const_iterator it, itE(simHitList.end());
    for ( it = simHitList.begin(); it != itE; ++it ) {
      //Fixme: momentum() < 0 (i.e. not present);
      //Fixme: Possibly add points for both posStart() and posEnd() (and use energy loss information to get different momenta?)
      addPars(trackStateOnSurfaces,it->second->createTrackParameters());
    }
 
  }
 
  createTrack(trackStateOnSurfaces);
}

Member Data Documentation

ascObjs

std::vector<AscObj_TruthPoint*>* TrackHandle_TruthTrack::Imp::ascObjs

Definition at line 45 of file TrackHandle_TruthTrack.cxx.

ascObjVis

bool TrackHandle_TruthTrack::Imp::ascObjVis

Definition at line 44 of file TrackHandle_TruthTrack.cxx.

genParticle

HepMC::ConstGenParticlePtr TrackHandle_TruthTrack::Imp::genParticle

Definition at line 42 of file TrackHandle_TruthTrack.cxx.

simBarCode

SimBarCode TrackHandle_TruthTrack::Imp::simBarCode

Definition at line 40 of file TrackHandle_TruthTrack.cxx.

simHitList

SimHitList TrackHandle_TruthTrack::Imp::simHitList

Definition at line 41 of file TrackHandle_TruthTrack.cxx.

theclass

TrackHandle_TruthTrack* TrackHandle_TruthTrack::Imp::theclass

Definition at line 39 of file TrackHandle_TruthTrack.cxx.

trkTrack

const Trk::Track* TrackHandle_TruthTrack::Imp::trkTrack

Definition at line 47 of file TrackHandle_TruthTrack.cxx.

---

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

• TrackHandle_TruthTrack.cxx