PndHypIdealTrackingTask Class Reference

#include <PndHypIdealTrackingTask.h>

Inheritance diagram for PndHypIdealTrackingTask:

Public Member Functions
	PndHypIdealTrackingTask ()
	~PndHypIdealTrackingTask ()
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
Double_t	GetTrackDip (PndMCTrack *myTrack)
Double_t	GetTrackCurvature (PndMCTrack *myTrack)

Private Member Functions
void	Register ()
void	Reset ()
void	ProduceHits ()
	ClassDef (PndHypIdealTrackingTask, 1)

Private Attributes
TString	fBranchName
TClonesArray *	fHitArray
TClonesArray *	fDigiStripArray
TClonesArray *	fMctruthArray
TClonesArray *	fClusterStripArray
TClonesArray *	fPointArray
TClonesArray *	fTrackOutputArray

Detailed Description

Definition at line 28 of file PndHypIdealTrackingTask.h.

Constructor & Destructor Documentation

PndHypIdealTrackingTask::PndHypIdealTrackingTask

(

)

Default constructor

Definition at line 33 of file PndHypIdealTrackingTask.cxx.

References fBranchName.

                                                 :
FairTask("Digitization task for PANDA PndHyp") 
{
  //fBranchName         = "HYPStripClusterHit";
        fBranchName     = "";
}

PndHypIdealTrackingTask::~PndHypIdealTrackingTask

(

)

Destructor

Definition at line 44 of file PndHypIdealTrackingTask.cxx.

{
}

Member Function Documentation

PndHypIdealTrackingTask::ClassDef ( PndHypIdealTrackingTask  , 1    )

private

void PndHypIdealTrackingTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

detector number hardcoded TODO: convention

Definition at line 105 of file PndHypIdealTrackingTask.cxx.

References Double_t, fClusterStripArray, fDigiStripArray, fHitArray, fMctruthArray, fPointArray, fTrackOutputArray, fVerbose, PndHypHit::GetBotIndex(), GetTrackCurvature(), GetTrackDip(), hit, mom, point, and pos.

{
        if ( ! fTrackOutputArray ) Fatal("Exec", "No fTrackOutputArray");
        fTrackOutputArray->Delete();
  
  // CREATE MONTECARLO MAP
  std::map<Int_t, PndMCTrack*> mcmap;
  Int_t nmc=fMctruthArray->GetEntriesFast();
  for(Int_t imc=0;imc<nmc;++imc){
    PndMCTrack* mc=(PndMCTrack*)fMctruthArray->At(imc);
    mcmap[imc]=mc;
    //if(mc->GetPdgCode()>100)cout<<imc<<" "<<mc->GetPdgCode()<<endl;
    
  }
  
  // BUILD TRACK Candidates
  std::map<Int_t, GFTrackCand*> trackMap;
  
  Int_t nPndHypHits=fHitArray->GetEntriesFast();
  cout<<nPndHypHits<<endl;
  
  
  for(Int_t iHit=0;iHit<nPndHypHits;++iHit)
  {
    PndHypHit* hit=(PndHypHit*)fHitArray->At(iHit);
    Int_t clustid=hit->GetRefIndex();
    Int_t digiid,pointid,MCid;
    
    
    if(hit->GetBotIndex()>0)
    {
      digiid=((PndHypCluster*)fClusterStripArray->At(clustid))->GetDigiIndex(0);
      pointid=((PndHypDigi*)fDigiStripArray->At(digiid))->GetIndex();
    }
    
    //cout<<iHit<<" "<<pointid<<" "<<fPointArray->GetEntriesFast()<<endl;
    
    
    //if(pointid >fPointArray->GetEntriesFast())continue;
    
    
    PndHypPoint* point = (PndHypPoint*)fPointArray->At(pointid);
    MCid = point->GetTrackID();
    
    //if(point)cout<<point->GetVolumeID()<<endl;
    
    // cut on secondaries (deltas) etc
    if(MCid<0)continue;
    
    cout<<MCid<<endl;
    TVector3 pos;
    hit->Position(pos);
    
    // build tracks here:
    
    // check if track already candidated:
    if(trackMap[MCid]==0)
    {Int_t size = fTrackOutputArray->GetEntriesFast();
      //caution this should be an intrinsic assignment, so trackMap and
      //fTrackOutputArray contain the same pointers.
      trackMap[MCid]=new ((*fTrackOutputArray)[size]) GFTrackCand();
      
      Int_t pdgcode;
      if( mcmap[MCid]==0)continue;
      
      pdgcode = mcmap[MCid]->GetPdgCode();
      cout<<pdgcode<<endl;
      
      TParticlePDG pdgpart(pdgcode);
      Int_t mcCharge = (Int_t)pdgpart.Charge();
      
      TVector3 vtx = mcmap[MCid]->GetStartVertex();
      TVector3 mom = mcmap[MCid]->GetMomentum();
      
      
      Double_t invp = 1 / mom.Mag();
      Double_t dxdz = invp*mom.x()*mom.z();
      Double_t dydz = invp*mom.y()*mom.z();
      Double_t sigx = 0.01, sigy = 0.01, sigdxdz = 0.01, sigdydz = 0.01, siginvp = 0.01;
      
      sigx*=vtx.x();
      sigy*=vtx.y();
      sigdxdz*=dxdz;
      sigdydz*=dydz;
      siginvp*=invp;
      
      if(fVerbose>1) std::cout<<"PndHypIdealTrackingTask::Exec(), Particle Info: "
                                <<pdgpart.GetName()<<" ("<<pdgcode<<") Q="<<mcCharge<<std::endl;
      
      trackMap[MCid]->setCurv(GetTrackCurvature(mcmap[MCid]));
      trackMap[MCid]->setDip(GetTrackDip(mcmap[MCid]));
      trackMap[MCid]->setInverted(false);
      
      
      
    }
    
    // add the hit index
    //       trackMap[MCid]->addHit(,hit->GetDetectorID(),iHit);
    //       PndHypRecoHit arecohit(hit);
    
    trackMap[MCid]->addHit(3,iHit);
    cout<<iHit<<endl;
    
    
    
  }//end for hit
  
  
  std::cout<<fTrackOutputArray->GetEntriesFast()
  <<" Tracks created"<<std::endl;
  
  
  
  
}

Double_t PndHypIdealTrackingTask::GetTrackCurvature

(

PndMCTrack *

myTrack

)

Definition at line 223 of file PndHypIdealTrackingTask.cxx.

References PndMCTrack::GetMomentum(), p, and CAMath::Sqrt().

Referenced by Exec().

{
  TVector3 p = myTrack->GetMomentum();
  return (2/TMath::Sqrt(p.Px()*p.Px() + p.Py()*p.Py()));
}

Double_t PndHypIdealTrackingTask::GetTrackDip

(

PndMCTrack *

myTrack

)

Definition at line 229 of file PndHypIdealTrackingTask.cxx.

References PndMCTrack::GetMomentum(), p, and CAMath::Sqrt().

Referenced by Exec().

{
  TVector3 p= myTrack->GetMomentum();
  return (p.Mag()/TMath::Sqrt(p.Px()*p.Px() + p.Py()*p.Py()));
}

InitStatus PndHypIdealTrackingTask::Init ( )

virtual

Definition at line 49 of file PndHypIdealTrackingTask.cxx.

References fClusterStripArray, fDigiStripArray, fHitArray, fMctruthArray, fPointArray, and fTrackOutputArray.

{
  // Get RootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) {
    std::cout << "-E- PndHypIdealTrackingTask::Init: "<< "RootManager not instantiated!" << std::endl;
  return kFATAL;  }
  
  // Get input array
  fHitArray = (TClonesArray*) ioman->GetObject("HYPHit");
  if ( ! fHitArray ) {
    std::cout << "-W- PndHypIdealTrackingTask::Init: "<< "No fHitarray!" << std::endl;
  return kERROR;  }
  
  fClusterStripArray = (TClonesArray*) ioman->GetObject("HYPStripClusterCand");
  if ( ! fClusterStripArray ) {
    std::cout << "-W- PndMvdIdealTrackingTask::Init: "<< "No fClusterStripArray!" << std::endl;
  return kERROR;  }
  
  fDigiStripArray = (TClonesArray*) ioman->GetObject("HypStripDigis");
  if ( ! fDigiStripArray ) {
    std::cout << "-W- PndMvdIdealTrackingTask::Init: "<< "No fDigiStripArray array!" << std::endl;
  return kERROR;  }
  
  fPointArray = (TClonesArray*) ioman->GetObject("HypPoint");
  if ( ! fPointArray )
  {
    std::cout << "-W- PndHypStripHitProducer::Init: "
    << "No MVDPoint array!" << std::endl;
    return kERROR;
  }
  
  // Get MCTruth collection
  fMctruthArray=(TClonesArray*) ioman->GetObject("MCTrack");
  if(fMctruthArray==0)  {
    std::cout << "-W- PndHypIdealTrackingTask::Init: "<< "No McTruth array!" << std::endl;
  return kERROR;  }
  
  // Create and register output array
  fTrackOutputArray = new TClonesArray("GFTrackCand");
  ioman->Register("HypTrackCand", "PndHyp ideal tracklets", fTrackOutputArray, kTRUE);
  
  
  return kSUCCESS;
}

void PndHypIdealTrackingTask::ProduceHits ( )

private

void PndHypIdealTrackingTask::Register ( )

private

void PndHypIdealTrackingTask::Reset ( )

private

void PndHypIdealTrackingTask::SetParContainers ( )

virtual

Virtual method Init

Definition at line 95 of file PndHypIdealTrackingTask.cxx.

References rtdb.

{
  // Get Base Container
  FairRunAna* ana = FairRunAna::Instance();
  FairRuntimeDb* rtdb=ana->GetRuntimeDb();
  
}

Member Data Documentation

TString PndHypIdealTrackingTask::fBranchName

private

Definition at line 53 of file PndHypIdealTrackingTask.h.

Referenced by PndHypIdealTrackingTask().

TClonesArray* PndHypIdealTrackingTask::fClusterStripArray

private

Definition at line 56 of file PndHypIdealTrackingTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypIdealTrackingTask::fDigiStripArray

private

Definition at line 55 of file PndHypIdealTrackingTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypIdealTrackingTask::fHitArray

private

Input array of Points

Definition at line 55 of file PndHypIdealTrackingTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypIdealTrackingTask::fMctruthArray

private

Definition at line 56 of file PndHypIdealTrackingTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypIdealTrackingTask::fPointArray

private

Definition at line 57 of file PndHypIdealTrackingTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndHypIdealTrackingTask::fTrackOutputArray

private

Output array of Hits

Definition at line 60 of file PndHypIdealTrackingTask.h.

Referenced by Exec(), and Init().

---

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

• PndHypIdealTrackingTask.h
• PndHypIdealTrackingTask.cxx