Task class for track finding. More...

#include <PndGemIdealPatternRecoTask.h>

Inheritance diagram for PndGemIdealPatternRecoTask:

Public Member Functions
	PndGemIdealPatternRecoTask ()

virtual	~PndGemIdealPatternRecoTask ()

virtual InitStatus	Init ()

virtual void	Exec (Option_t *opt)

virtual void	SetParContainers ()

Private Member Functions
	ClassDef (PndGemIdealPatternRecoTask, 1)

Private Attributes
TClonesArray *	fGemHitArr

TClonesArray *	fMCTrackArr

TClonesArray *	fGemPointArr

TClonesArray *	fPndTrackCandArr

Detailed Description

Task class for track finding.

Author: A.Bubak arkad.nosp@m.iusz.nosp@m..buba.nosp@m.k@gm.nosp@m.ail.c.nosp@m.om

Date: 27.01.2008 Ideal Pattern reco task

Definition at line 26 of file PndGemIdealPatternRecoTask.h.

Constructor & Destructor Documentation

PndGemIdealPatternRecoTask::PndGemIdealPatternRecoTask ( )

Default constructor

Definition at line 29 of file PndGemIdealPatternRecoTask.cxx.

References fGemHitArr, fGemPointArr, fMCTrackArr, and fPndTrackCandArr.

                                                        {
   fMCTrackArr = NULL;
   fGemHitArr = NULL;
   fGemPointArr = NULL;
   fPndTrackCandArr = NULL;
 }

PndGemIdealPatternRecoTask::~PndGemIdealPatternRecoTask ( )

virtual

Destructor

Definition at line 44 of file PndGemIdealPatternRecoTask.cxx.

References fGemHitArr, fGemPointArr, fMCTrackArr, and fPndTrackCandArr.

                                                         {
   fMCTrackArr->Delete();
   fGemHitArr->Delete();
   fGemPointArr->Delete();
   fPndTrackCandArr->Delete();
   if(fMCTrackArr != NULL) delete fMCTrackArr;
   if(fGemHitArr != NULL) delete fGemHitArr;
   if(fGemPointArr != NULL) delete fGemPointArr;
   if(fPndTrackCandArr != NULL) delete fPndTrackCandArr;
 }

Member Function Documentation

PndGemIdealPatternRecoTask::ClassDef	(	PndGemIdealPatternRecoTask	,
		1
	)

private

void PndGemIdealPatternRecoTask::Exec ( Option_t * opt )

virtual

Task execution

Definition at line 110 of file PndGemIdealPatternRecoTask.cxx.

References PndTrackCand::AddHit(), fGemHitArr, fGemPointArr, fPndTrackCandArr, PndGemHit::GetPosition(), hit, nhits, npoints, and point.

                                                {
   fPndTrackCandArr->Delete();
  
   std::map<int,PndTrackCand*> cands;
 
   int nhits = fGemHitArr->GetEntriesFast();
   int npoints = fGemPointArr->GetEntriesFast();
   for(int ihit=0;ihit<nhits;++ihit){
     PndGemHit* gemHit = (PndGemHit*)fGemHitArr->At(ihit);
     FairHit * hit = gemHit;
     int index = hit->GetRefIndex();
     assert(index<npoints);
     PndGemMCPoint *gemPoint = (PndGemMCPoint*)fGemPointArr->At(index);
     FairMCPoint *point = gemPoint;
     int trackId = point->GetTrackID();
     if(cands[trackId]==NULL) {
       cands[trackId] = new PndTrackCand;
       //      cands[trackId]->setMcTrackId(trackId);
       //PndMCTrack *myMCTrack = (PndMCTrack*)fMCTrackArr->At(trackId); //[R.K. 03/2017] unused variable?
       //int pdg = myMCTrack->GetPdgCode(); //[R.K.02/2017] Unused variable?
       //double charge; //[R.K. 01/2017] unused variable?
       //if(pdg<100000000){ //[R.K. 01/2017] unused variable?
         //charge =  TDatabasePDG::Instance()->GetParticle(pdg)->Charge()/3.; //[R.K. 01/2017] unused variable?
       //} //[R.K. 01/2017] unused variable?
       //else{ //[R.K. 01/2017] unused variable?
         //charge = 0.; //[R.K. 01/2017] unused variable?
       //} //[R.K. 01/2017] unused variable?
 
     }
     cands[trackId]->AddHit(FairRootManager::Instance()->GetBranchId("GEMHit"),ihit,gemHit->GetPosition().Mag());
   }
 
   for(std::map<int,PndTrackCand*>::iterator iter=cands.begin();
       iter!=cands.end();++iter){
     int nentries = fPndTrackCandArr->GetEntriesFast();
     if(iter->second->GetNHits()>0) new ((*fPndTrackCandArr)[nentries]) PndTrackCand(*(iter->second));
     delete iter->second;
   }
 }

InitStatus PndGemIdealPatternRecoTask::Init ( )

virtual

Initialisation at beginning of each event

Definition at line 59 of file PndGemIdealPatternRecoTask.cxx.

References fGemHitArr, fGemPointArr, fMCTrackArr, and fPndTrackCandArr.

                                             {
 
         // Get and check FairRootManager
   FairRootManager* ioman = FairRootManager::Instance();
   if( !ioman ) {
     std::cout << "-E- PndGemIdealPatternRecoTask::Init: "
               << "RootManager not instantised!" << std::endl;
     return kFATAL;
   }
 
   fMCTrackArr  = (TClonesArray*) ioman->GetObject("MCTrack");
   if ( !fMCTrackArr ) {
     std::cout << "-W- PndGemIdealPatternRecoTask::Init: No MCTrack array!" << std::endl;
     return kERROR;
   }
 
   fGemPointArr  = (TClonesArray*) ioman->GetObject("GEMPoint");
   if ( !fGemPointArr ) {
     std::cout << "-W- PndGemIdealPatternRecoTask::Init: No PndGemPoint array!" << std::endl;
     return kERROR;
   }
   fGemHitArr  = (TClonesArray*) ioman->GetObject("GEMHit");
   if ( !fGemHitArr ) {
     std::cout << "-W- PndGemIdealPatternRecoTask::Init: No PndGemHit array!" << std::endl;
     return kERROR;
   }
 
   // Create and register GemTrack array
   fPndTrackCandArr = new TClonesArray("PndTrackCand",100);
   ioman->Register("GemPndTrackCand", "Gem", fPndTrackCandArr  , kTRUE);
   
   std::cout << "PndGemIdealPatternRecoTask::Init() finished" << std::endl;
   return kSUCCESS;  
   
 }

void PndGemIdealPatternRecoTask::SetParContainers ( )

virtual

SetParContainers

Definition at line 99 of file PndGemIdealPatternRecoTask.cxx.

References rtdb.

                                                   {
   FairRunAna* ana = FairRunAna::Instance();
   FairRuntimeDb* rtdb = ana->GetRuntimeDb();
   rtdb->getContainer("FairBaseParSet");
   rtdb->getContainer("PndGeoPassivePar");
 }