PndGenfitAdapters.cxx File Reference

#include "PndGenfitAdapters.h"
#include <iostream>
#include "GFTrack.h"
#include "GFAbsTrackRep.h"
#include "GFTrackCand.h"
#include "PndTrack.h"
#include "PndTrackCand.h"
#include "GFDetPlane.h"
#include "GFException.h"
#include "TMatrixT.h"
#include "FairTrackParP.h"
#include "GeaneTrackRep.h"
#include "RKTrackRep.h"
#include <cmath>

Go to the source code of this file.

Functions
PndTrackCand *	GenfitTrackCand2PndTrackCand (const GFTrackCand *cand)
GFTrackCand *	PndTrackCand2GenfitTrackCand (PndTrackCand *cand)
PndTrack *	GenfitTrack2PndTrack (const GFTrack *tr)

Function Documentation

PndTrack* GenfitTrack2PndTrack

(

const GFTrack *

tr

)

Definition at line 43 of file PndGenfitAdapters.cxx.

References GFAbsTrackRep::clone(), count, GenfitTrackCand2PndTrackCand(), GFTrack::getCand(), GFTrack::getCardinalRep(), GFTrack::getChiSqu(), GFAbsTrackRep::getFirstCov(), GFAbsTrackRep::getFirstPlane(), GFAbsTrackRep::getFirstState(), GFAbsTrackRep::getLastCov(), GFAbsTrackRep::getLastPlane(), GFAbsTrackRep::getLastState(), GFAbsTrackRep::getMom(), GFTrack::getNDF(), GFDetPlane::getNormal(), GFDetPlane::getO(), GFDetPlane::getU(), GFDetPlane::getV(), i, PndTrack::SetChi2(), PndTrack::SetFlag(), PndTrack::SetNDF(), and GFException::what().

Referenced by PndLmdKalmanTask::Exec(), PndRecoDafFit::Fit(), and PndRecoKalmanFit::Fit().

                                                 {
  GFAbsTrackRep* clone = tr->getCardinalRep()->clone();
  TMatrixT<double> firstState = clone->getFirstState();
  TMatrixT<double> lastState = clone->getLastState();
  TMatrixT<double> firstCov = clone->getFirstCov();
  TMatrixT<double> lastCov = clone->getLastCov();
  GFDetPlane firstPlane = clone->getFirstPlane();
  GFDetPlane lastPlane = clone->getLastPlane();
  
  GFAbsTrackRep* gtr;
  if (dynamic_cast<GeaneTrackRep*>(clone)!=NULL)
          gtr = dynamic_cast<GeaneTrackRep*>(clone);
  else if (dynamic_cast<RKTrackRep*>(clone)!=NULL)
          gtr = dynamic_cast<RKTrackRep*>(clone);
  else {
    std::cerr << " GenfitGFAbsTrackRep2PndTrack() can currently only handle GeaneTrackRep and RKTrackRep" << std::endl;
    throw;
  }

  //make the FairTrackParP for first and last hit
  double firstCova[15];
  int count=0;;
  for(int i=0; i<5;++i){
     for(int j=i;j<5;++j){
           firstCova[count++]=firstCov[i][j];
     }
  }
  double lastCova[15];
  count=0;;
  for(int i=0; i<5;++i){
    for(int j=i;j<5;++j){
          lastCova[count++]=lastCov[i][j];
    }
  }

  //  calculation of spu = sign[p(DJ x DK)]
  double first_pro(0), last_pro(0), first_spu, last_spu;
  bool exc(false);

  try{
    first_pro = gtr->getMom(firstPlane).Dot(firstPlane.getNormal());
    last_pro = gtr->getMom(lastPlane).Dot(lastPlane.getNormal());
  }
  catch (GFException& e){
    exc=true;
    std::cerr<<"*** PndGenfitAdapters::GenfitTrack2PndTrack" << "\t" << "could not convert GenfitTrack to PndTrack"<<std::endl;
    e.what();
  }

  first_spu = (first_pro>=0) ? 1 : -1;
  last_spu = (last_pro>=0) ? 1 : -1;
    
  FairTrackParP first(firstState[3][0],firstState[4][0],firstState[1][0],firstState[2][0],firstState[0][0],firstCova,firstPlane.getO(),firstPlane.getU(),firstPlane.getV(),first_spu);
  FairTrackParP last(lastState[3][0],lastState[4][0],lastState[1][0],lastState[2][0],lastState[0][0],lastCova,lastPlane.getO(),lastPlane.getU(),lastPlane.getV(),last_spu);
    
  //copy the trackCand
  GFTrackCand genfitCand = tr->getCand();
  PndTrackCand* pndCand = GenfitTrackCand2PndTrackCand(&genfitCand);
  PndTrack* retVal =  new PndTrack(first,last,*pndCand);
  retVal->SetChi2(tr->getChiSqu());
  retVal->SetNDF(tr->getNDF());
  if (tr->getNDF()==0)
    {
      retVal->SetFlag(-15);
    }
  else if (exc) {
    retVal->SetFlag(-1);
  }
  else {
          retVal->SetFlag(1);
  }

  delete pndCand;
  return retVal;
  
}

PndTrackCand* GenfitTrackCand2PndTrackCand

(

const GFTrackCand *

cand

)

Definition at line 19 of file PndGenfitAdapters.cxx.

References PndTrackCand::AddHit(), GFTrackCand::getHit(), GFTrackCand::getMcTrackId(), GFTrackCand::getNHits(), i, nhits, and PndTrackCand::setMcTrackId().

Referenced by GenfitTrack2PndTrack().

                                                                   {
  PndTrackCand* retVal = new PndTrackCand();
  unsigned int nhits = cand->getNHits();
  unsigned detId,hitId;
  double rho;
  for(unsigned int i=0;i<nhits;++i){
    cand->getHit(i,detId,hitId,rho);
    retVal->AddHit(detId,hitId,rho);
  }
  retVal->setMcTrackId(cand->getMcTrackId());
  return retVal;
}

GFTrackCand* PndTrackCand2GenfitTrackCand

(

PndTrackCand *

cand

)

Definition at line 32 of file PndGenfitAdapters.cxx.

References GFTrackCand::addHit(), PndTrackCandHit::GetDetId(), PndTrackCandHit::GetHitId(), PndTrackCand::getMcTrackId(), PndTrackCand::GetNHits(), PndTrackCandHit::GetRho(), PndTrackCand::GetSortedHit(), i, nhits, and GFTrackCand::setMcTrackId().

Referenced by PndLmdKalmanTask::Exec(), PndRecoDafFit::Fit(), and PndRecoKalmanFit::Fit().

                                                             {
  GFTrackCand* retVal = new GFTrackCand();
  unsigned int nhits = cand->GetNHits();
  for(unsigned int i=0;i<nhits;++i){
    PndTrackCandHit candHit = cand->GetSortedHit(i);
    retVal->addHit(candHit.GetDetId(),candHit.GetHitId(),candHit.GetRho(),i);
  }
  retVal->setMcTrackId(cand->getMcTrackId());
  return retVal;
}