PndIdealTrackFinder Class Reference

Ideal track finder for all types of tracking detectors The PndIdealTrackFinder combines all hits in given branches into PndTrackCand and calculates the PndTrack based on MC information. More...

#include <PndIdealTrackFinder.h>

Inheritance diagram for PndIdealTrackFinder:

Public Member Functions
	PndIdealTrackFinder ()
virtual	~PndIdealTrackFinder ()
virtual InitStatus	Init ()
virtual void	AddBranchName (TString name)
	Search for tracks only in given branches. If no BranchName is given all tracking detectors are taken. More...
virtual void	AddBranchName (TString name, TString mcName)
	Search for tracks only in given branches. If no BranchName is given all tracking detectors are taken. More...
virtual void	SetOutputBranchName (TString name)
virtual void	Exec (Option_t *opt)
void	SetMomentumSmearing (Double_t sigmax=-1., Double_t sigmay=-1., Double_t sigmaz=-1.)
void	SetRelativeMomentumSmearing (Double_t dpop=-1.)
void	SetVertexSmearing (Double_t sigmax=-1., Double_t sigmay=-1., Double_t sigmaz=-1.)
void	SetTrackingEfficiency (Double_t eff=1.)
void	SetTrackSelector (TString selector)
void	SetPersistence (Bool_t persistence)
void	SetPersistency (Bool_t val=kTRUE)
Bool_t	GetPersistency ()

Protected Member Functions
virtual void	CreateTrackCands ()
virtual void	CreateTracks ()
virtual void	FilterTrackCands ()
virtual FairMCPoint *	GetFairMCPoint (TString hitBranch, FairMultiLinkedData_Interface *links, FairMultiLinkedData &array)
virtual void	SmearVector (TVector3 &vec, const TVector3 &sigma)
	ClassDef (PndIdealTrackFinder, 2)

Protected Attributes
TString	fOutBranchName
TClonesArray *	fTrackCand
TClonesArray *	fTrack
TClonesArray *	fMCTrack
PndTrackFunctor *	fTrackSelector
std::map< TString, TClonesArray * >	fBranchMap
std::vector< TString >	fBranchNames
std::map< FairLink, PndTrackCand >	fTrackCandMap
std::map< FairLink, FairMCPoint >	fFirstPointMap
std::map< FairLink, FairMCPoint >	fLastPointMap
std::map< TString, TString >	fPointBranchMap
TDatabasePDG *	fPdg
	! Particle DB More...
Int_t	fHitCount
TVector3	fMomSigma
	Momentum smearing sigma [GeV]. More...
Double_t	fDPoP
	Relative momentum Smearing. More...
Bool_t	fRelative
	flag More...
TVector3	fVtxSigma
	Vertex smearing sigma [cm]. More...
Double_t	fEfficiency
	Tracking efficiency - if (0 <= e < 1), some tracks will be discarded. More...

Detailed Description

Ideal track finder for all types of tracking detectors The PndIdealTrackFinder combines all hits in given branches into PndTrackCand and calculates the PndTrack based on MC information.

PndIdealTrackFinder.h

Author

Tobias Stockmanns t.sto.nosp@m.ckma.nosp@m.nns@f.nosp@m.z-ju.nosp@m.elich.nosp@m..de

Date

Apr 12, 2010

Definition at line 25 of file PndIdealTrackFinder.h.

Constructor & Destructor Documentation

PndIdealTrackFinder::PndIdealTrackFinder

(

)

Definition at line 22 of file PndIdealTrackFinder.cxx.

References fPointBranchMap, and PndPersistencyTask::SetPersistency().

                                         :
  fOutBranchName("IdealTrack"), fTrackCand(0), fTrack(0), fMCTrack(0), fTrackSelector(0), fPdg(0), fHitCount(0), fMomSigma(0,0,0), fDPoP(0.), fRelative (kFALSE), fVtxSigma(0,0,0), fEfficiency(1.)
{
        fPointBranchMap["MVDHitsPixel"] = "MVDPoint";
        fPointBranchMap["MVDHitsStrip"] = "MVDPoint";
        fPointBranchMap["STTHit"] = "STTPoint";
        fPointBranchMap["GEMHit"] = "GEMPoint";
        fPointBranchMap["FTSHit"] = "FTSPoint";

        fPointBranchMap["SciTHit"] = "SciTPoint";
        fPointBranchMap["MdtHit"] = "MdtPoint";

        fPointBranchMap["DircHit"] = "";         // no FairLinks Provided!
        fPointBranchMap["FTofHit"] = "";         // no FairLinks Provided!
        fPointBranchMap["RichHit"] = "";         // no FairLinks Provided!

        SetPersistency(kTRUE);
}

PndIdealTrackFinder::~PndIdealTrackFinder ( )

virtual

Definition at line 41 of file PndIdealTrackFinder.cxx.

References fTrackSelector.

                                          {
        if (fTrackSelector != 0)
                delete fTrackSelector;
}

Member Function Documentation

virtual void PndIdealTrackFinder::AddBranchName ( TString  name )

inlinevirtual

Search for tracks only in given branches. If no BranchName is given all tracking detectors are taken.

Definition at line 33 of file PndIdealTrackFinder.h.

References fBranchNames.

Referenced by idealbarreltracking(), idealforwardtracking(), Init(), PndMasterRecoTask::PndMasterRecoTask(), reco_complete(), reco_complete_newSTT(), recoqa_ca(), recoqa_complete(), and trackingQANewLinks().

virtual void PndIdealTrackFinder::AddBranchName ( TString  name, TString  mcName  )

inlinevirtual

Search for tracks only in given branches. If no BranchName is given all tracking detectors are taken.

Definition at line 34 of file PndIdealTrackFinder.h.

References fBranchNames, fPointBranchMap, and name.

PndIdealTrackFinder::ClassDef ( PndIdealTrackFinder  , 2    )

protected

void PndIdealTrackFinder::CreateTrackCands ( )

protectedvirtual

Definition at line 103 of file PndIdealTrackFinder.cxx.

References fBranchMap, fFirstPointMap, fHitCount, fLastPointMap, fTrackCandMap, GetFairMCPoint(), i, point, and TString.

Referenced by PndCloneGenerator::Exec(), and Exec().

{
        fHitCount = 0;
        for (std::map<TString, TClonesArray*>::iterator iter = fBranchMap.begin(); iter != fBranchMap.end(); iter++){
                //std::cout << "AddLinks from Branch: " << iter->first << std::endl;
                for (int i = 0; i < iter->second->GetEntriesFast(); i++){
                        
                        FairMultiLinkedData array;
                        FairMultiLinkedData_Interface* links = (FairMultiLinkedData_Interface*)iter->second->At(i);
                        TString hitBranch = iter->first;

                        FairMCPoint *point = GetFairMCPoint(hitBranch, links, array);
                        if (point == 0) {
                                continue;
                        }
                        FairMCPoint firstpoint = *point;
                        FairMCPoint lastpoint = *point;

                        double tof = point->GetTime();
                        delete(point);

                        for (int ipnt = 1; ipnt < array.GetNLinks(); ipnt++){
                          point = (FairMCPoint *) FairRootManager::Instance()->GetCloneOfLinkData(array.GetLink(ipnt));
                          tof += point->GetTime();
                          //  std::cout << ipnt << " " << tof << std::endl;
                          if( point->GetTime() < firstpoint.GetTime()) firstpoint = *point;
                          if( point->GetTime() > lastpoint.GetTime()) lastpoint = *point;
                          delete(point);
                        }
                        tof /= array.GetNLinks();
                        // std::cout << i << " " << tof << std::endl;
                        // .............................................
                        FairMultiLinkedData mctracks = links->GetLinksWithType(FairRootManager::Instance()->GetBranchId("MCTrack"));
                        for (int trackIndex = 0; trackIndex < mctracks.GetNLinks(); trackIndex++){
                                if (!fTrackCandMap.count(mctracks.GetLink(trackIndex))){
                                        fTrackCandMap[mctracks.GetLink(trackIndex)] = PndTrackCand();
                                        fTrackCandMap[mctracks.GetLink(trackIndex)].SetInsertHistory(kTRUE);
                                        fFirstPointMap[mctracks.GetLink(trackIndex)] = firstpoint;
                                        // fFirstPointMap[mctracks.GetLink(trackIndex)].SetInsertHistory(kTRUE);
                                        fLastPointMap[mctracks.GetLink(trackIndex)] = lastpoint;
                                        // fLastPointMap[mctracks.GetLink(trackIndex)].SetInsertHistory(kTRUE);
                                }
                                else {
                                  FairMCPoint tmpfirstpoint = fFirstPointMap[mctracks.GetLink(trackIndex)];
                                  if(firstpoint.GetTime() < tmpfirstpoint.GetTime()) fFirstPointMap[mctracks.GetLink(trackIndex)] = firstpoint;
                                  FairMCPoint tmplastpoint = fLastPointMap[mctracks.GetLink(trackIndex)];
                                  if(lastpoint.GetTime() > tmplastpoint.GetTime()) fLastPointMap[mctracks.GetLink(trackIndex)] = lastpoint;
                                }
                                FairLink link(-1, FairRootManager::Instance()->GetEntryNr(),FairRootManager::Instance()->GetBranchId(iter->first),i);
                                //std::cout << "CreateTrackCands " << mctracks.GetLink(trackIndex) << " : " << link << std::endl;
                                fTrackCandMap[mctracks.GetLink(trackIndex)].SetInsertHistory(kTRUE);
                                // fTrackCandMap[mctracks.GetLink(trackIndex)].AddHit(link, fHitCount++);                       //todo Rho is not properly calculated!
                                fTrackCandMap[mctracks.GetLink(trackIndex)].AddHit(link, tof);
                        }
                }
        }
}

void PndIdealTrackFinder::CreateTracks ( )

protectedvirtual

Definition at line 175 of file PndIdealTrackFinder.cxx.

References fBranchMap, fDPoP, fEfficiency, fFirstPointMap, fLastPointMap, fMomSigma, fPdg, fRelative, fTrack, fTrackCand, fTrackCandMap, fVtxSigma, PndTrackCandHit::GetDetId(), PndTrackCandHit::GetHitId(), PndTrackCand::GetNHits(), PndMCTrack::GetPdgCode(), PndTrackCand::GetSortedHit(), PndGemMCPoint::PositionOut(), PndTrackCand::setMcTrackId(), and SmearVector().

Referenced by PndCloneGenerator::Exec(), and Exec().

{
        for (std::map<FairLink, PndTrackCand>::iterator iter = fTrackCandMap.begin(); iter != fTrackCandMap.end(); iter++){
                PndTrackCand* myTrackCand = new((*fTrackCand)[fTrackCand->GetEntriesFast()]) PndTrackCand(iter->second);
                myTrackCand->setMcTrackId(iter->first.GetIndex());
                myTrackCand->AddLink(iter->first);
                myTrackCand->SetTimeStamp(FairRootManager::Instance()->GetEventTime());
                //std::cout << myTrackCand->GetLinksWithType(FairRootManager::Instance()->GetBranchId("MCTrack")) << " : " << std::endl;
                //myTrackCand->Print();
                //std::cout << *myTrackCand << std::endl;

                // ....... track
                PndMCTrack *mc = (PndMCTrack *) FairRootManager::Instance()->GetCloneOfLinkData(iter->first);
                int charge = 0;
                if (mc->GetPdgCode()<100000000) charge = (Int_t)TMath::Sign(1.0, ((TParticlePDG*) fPdg->GetParticle(mc->GetPdgCode()))->Charge());
                else charge = 1;


                if(0 < fEfficiency && fEfficiency < 1){
                  if(gRandom->Rndm() > fEfficiency) continue;
                }

                // first
                FairMCPoint firstpoint = fFirstPointMap[iter->first];
                TVector3 firstpos(0, 0, 0), firstmom(0, 0, 0);
                if(myTrackCand->GetSortedHit(0).GetDetId() == FairRootManager::Instance()->GetBranchId("GEMHit")) {
                  TClonesArray *gemhitarray = fBranchMap["GEMHit"];
                  Int_t hitid = myTrackCand->GetSortedHit(0).GetHitId();
                  PndGemHit *gemhit = (PndGemHit*) gemhitarray->At(hitid);
                  FairMultiLinkedData gemhitlink = gemhit->GetLinksWithType(FairRootManager::Instance()->GetBranchId("GEMPoint"));
                  PndGemMCPoint *gempoint = (PndGemMCPoint *) FairRootManager::Instance()->GetCloneOfLinkData(gemhitlink.GetLink(0));

                  TVector3 posin(0, 0, 0), posout(0, 0, 0);
                  gempoint->Position(posin);
                  gempoint->PositionOut(posout);
                  firstpos = 0.5 * (posin + posout);
                }
                else firstpoint.Position(firstpos);

                SmearVector(firstpos, fVtxSigma);
                firstpoint.Momentum(firstmom);

                if (fRelative) fMomSigma.SetXYZ(fDPoP*firstmom.Mag(),fDPoP*firstmom.Mag(),fDPoP*firstmom.Mag());
                SmearVector(firstmom, fMomSigma);
                FairTrackParP firstPar(firstpos, firstmom,
                                                          fVtxSigma, fMomSigma,
                                                          charge, firstpos,
                                                          TVector3(1.,0.,0.), TVector3(0.,1.,0.));
                // last
                FairMCPoint lastpoint = fLastPointMap[iter->first];
                TVector3 lastpos(0, 0, 0), lastmom(0, 0, 0);

                if(myTrackCand->GetSortedHit(myTrackCand->GetNHits() - 1).GetDetId() == FairRootManager::Instance()->GetBranchId("GEMHit")) {
                  TClonesArray *gemhitarray = fBranchMap["GEMHit"];
                  Int_t hitid = myTrackCand->GetSortedHit(myTrackCand->GetNHits() - 1).GetHitId();
                  PndGemHit *gemhit = (PndGemHit*) gemhitarray->At(hitid);
                  FairMultiLinkedData gemhitlink = gemhit->GetLinksWithType(FairRootManager::Instance()->GetBranchId("GEMPoint"));
                  PndGemMCPoint *gempoint = (PndGemMCPoint *) FairRootManager::Instance()->GetCloneOfLinkData(gemhitlink.GetLink(0));

                  TVector3 posin(0, 0, 0), posout(0, 0, 0);
                  gempoint->Position(posin);
                  gempoint->PositionOut(posout);
                  lastpos = 0.5 * (posin + posout);
                }
                else lastpoint.Position(lastpos);

                SmearVector(lastpos, fVtxSigma);
                lastpoint.Momentum(lastmom);
                SmearVector(lastmom, fMomSigma);
                FairTrackParP lastPar(lastpos, lastmom,
                                                         fVtxSigma, fMomSigma,
                                                         charge, lastpos,
                                                         TVector3(1.,0.,0.), TVector3(0.,1.,0.));



                new((*fTrack)[fTrack->GetEntriesFast()]) PndTrack(firstPar, lastPar, *myTrackCand, 0,0,1,mc->GetPdgCode(), -1,-1);
        }
}

void PndIdealTrackFinder::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Reimplemented in PndCloneGenerator.

Definition at line 89 of file PndIdealTrackFinder.cxx.

References CreateTrackCands(), CreateTracks(), FilterTrackCands(), fTrack, fTrackCand, and fTrackCandMap.

{
        fTrackCand->Delete();
        fTrack->Delete();
        fTrackCandMap.clear();

//      std::cout << "Event #" << FairRootManager::Instance()->GetEntryNr() << std::endl;
        CreateTrackCands();
        FilterTrackCands();

//      std::cout << "PndIdealTrackFinder:Found Tracks:" << std::endl;
        CreateTracks();
}

void PndIdealTrackFinder::FilterTrackCands ( )

protectedvirtual

Definition at line 161 of file PndIdealTrackFinder.cxx.

References fTrackCandMap, and fTrackSelector.

Referenced by Exec().

{
        if (fTrackSelector == 0){
                return;
        }
        for (std::map<FairLink, PndTrackCand>::iterator iter = fTrackCandMap.begin(); iter != fTrackCandMap.end();){
                if (!(*fTrackSelector)(iter->second.GetPointerToLinks(), true)){
                        fTrackCandMap.erase(iter++);
                } else {
                        ++iter;
                }
        }
}

FairMCPoint * PndIdealTrackFinder::GetFairMCPoint ( TString  hitBranch, FairMultiLinkedData_Interface *  links, FairMultiLinkedData &  array  )

protectedvirtual

Definition at line 255 of file PndIdealTrackFinder.cxx.

References fPointBranchMap.

Referenced by CreateTrackCands(), and PndCloneGenerator::SetFirstLastHit().

{
        // get the mc point(s) from each reco hit ......
        FairMultiLinkedData mcpoints = links->GetLinksWithType(FairRootManager::Instance()->GetBranchId(fPointBranchMap[hitBranch]));
        // std::cout << "hit " << i << " connected to points " << mvdpoints.GetNLinks() << " " << sttpoints.GetNLinks() << " " << gempoints.GetNLinks() << std::endl;

        // There seems to be a bug with ghost hits from the GEM stations. If more than one
        // MC point is associated to a hit, there is a good chance for false assignments
        // leading to wrong tracks. For the moment, skip hits with more than 1 GEM point.

        if (hitBranch == "GEMHit" && mcpoints.GetNLinks() > 1) return 0;
//      if ((*iter).first == "MVDHitsStrip" && mvdpoints.GetNLinks() > 1) return 0;

        array = mcpoints;


        if (array.GetNLinks() == 0){

                return 0;
        }
        return (FairMCPoint *) FairRootManager::Instance()->GetCloneOfLinkData(array.GetLink(0));
}

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

Referenced by PndLmdPixelHitProducerFast::GetPersistance(), PndMdtDigitization::Init(), PndMdtHitProducerIdeal::Init(), PndMdtClusterTask::Init(), PndFtsHitProducerRealFast::Init(), PndSttHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndRichHitProducer::Init(), PndSttHelixHitProducer::Init(), PndDiscTaskPID::Init(), Init(), PndSttMvdGemTracking::Init(), PndMdtTrkProducer::Init(), PndFtsHitProducerRealFull::Init(), PndLmdPixelClusterTask::Init(), PndSttHitProducerRealFull::Init(), PndLmdStripClusterTask::Init(), PndEmcApdHitProducer::Init(), PndMissingPzCleanerTask::Init(), PndEmcMakeRecoHit::Init(), PndEmcMakeClusterOnline::Init(), PndTrackSmearTask::Init(), PndEmcFWEndcapTimebasedWaveforms::Init(), PndSttHitProducerIdeal::Init(), PndEmcFWEndcapDigi::Init(), PndFtsHitProducerIdeal::Init(), PndEmcMakeCluster::Init(), PndMdtPointsToWaveform::Init(), PndDiscTaskDigitization::Init(), PndEmcMakeDigi::Init(), PndSdsTimeWalkCorrTask::Init(), PndLmdPixelHitProducerFast::Init(), PndDrcHitFinder::Init(), PndRichHitFinder::Init(), PndEmcMakeCorr::Init(), PndFtofHitProducerIdeal::Init(), PndEmcHitsToWaveform::Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), PndEmcHitProducer::Init(), PndRecoMultiKalmanTask2::Init(), PndDrcHitProducerReal::Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndDrcTimeDigiTask::Init(), PndRecoKalmanTask2::Init(), PndEmcExpClusterSplitter::Init(), PndSdsNoiseProducer::Init(), PndFtsHoughTrackerTask::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsHybridHitProducer::Init(), PndSdsIdealRecoTask::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndSdsStripHitProducer::Init(), PndGemDigitize::Init(), PndGemFindHits::Init(), PndSdsPixelClusterTask::Init(), PndSdsStripClusterTask::Init(), PndMvdGemTrackFinderOnHits::Init(), PndBarrelTrackFinder::Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

{ return fPersistency; }

InitStatus PndIdealTrackFinder::Init ( )

virtual

Virtual method Init

Definition at line 47 of file PndIdealTrackFinder.cxx.

References AddBranchName(), fBranchMap, fBranchNames, fMCTrack, fOutBranchName, fPdg, fPointBranchMap, fTrack, fTrackCand, fTrackSelector, PndPersistencyTask::GetPersistency(), and i.

{

        FairRootManager* ioman = FairRootManager::Instance();
        if (!ioman) {
                std::cout << "-E- PndMCTestHitCompare::Init: "
                                << "RootManager not instantiated!" << std::endl;
                return kFATAL;
        }

        if (fBranchNames.size() == 0){
            // Use hits of all tracking subsystems if nothing is given
                AddBranchName("MVDHitsPixel");
                AddBranchName("MVDHitsStrip");
                AddBranchName("STTHit");
        AddBranchName("GEMHit");
        AddBranchName("FTSHit");
        }

        for (size_t i = 0; i < fBranchNames.size(); i++){
                if (ioman->GetObject(fBranchNames[i]) != 0){
                        fBranchMap[fBranchNames[i]] = (TClonesArray*)ioman->GetObject(fBranchNames[i]);
                        ioman->GetObject(fPointBranchMap[fBranchNames[i]]);  // initialise the used FairMcPoint Branches
                }
        }
        fMCTrack = (TClonesArray*)ioman->GetObject("MCTrack");


        fTrackCand = new TClonesArray("PndTrackCand");
        ioman->Register(fOutBranchName + "Cand", "MC", fTrackCand, GetPersistency());
        fTrack = new TClonesArray("PndTrack");
        ioman->Register(fOutBranchName, "MC", fTrack, GetPersistency());

        if (fTrackSelector == 0){
                std::cout << "-W- PndIdealTrackFinder::Init() no fTrackSelector set! All possible tracks will be taken!" << std::endl;
        }
 
        fPdg = new TDatabasePDG();

        return kSUCCESS;
}

void PndIdealTrackFinder::SetMomentumSmearing ( Double_t  sigmax = -1., Double_t  sigmay = -1., Double_t  sigmaz = -1.  )

inline

Definition at line 44 of file PndIdealTrackFinder.h.

References fabs(), fMomSigma, and fRelative.

                                                                                                      { 
            fMomSigma.SetXYZ(fabs(sigmax),fabs(sigmay),fabs(sigmaz)); 
            fRelative=kFALSE; 
          }; // in GeV

virtual void PndIdealTrackFinder::SetOutputBranchName ( TString  name )

inlinevirtual

Definition at line 36 of file PndIdealTrackFinder.h.

References fOutBranchName, and name.

Referenced by idealbarreltracking(), idealforwardtracking(), PndMasterRecoIdealTask::PndMasterRecoIdealTask(), PndMasterRecoTask::PndMasterRecoTask(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), reco_rich(), and recoideal_complete().

{ fOutBranchName = name; };

void PndIdealTrackFinder::SetPersistence ( Bool_t  persistence )

inline

Definition at line 63 of file PndIdealTrackFinder.h.

References PndPersistencyTask::SetPersistency().

Referenced by idealbarreltracking(), idealforwardtracking(), PndMasterRecoIdealTask::PndMasterRecoIdealTask(), reco_complete(), reco_complete_newSTT(), reco_rich(), and recoideal_complete().

{ SetPersistency(persistence); }

void PndPersistencyTask::SetPersistency ( Bool_t  val = kTRUE )

inlineinherited

Definition at line 31 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency, and val.

Referenced by barrelTrackFinder(), digi_complete(), digi_complete_newSTT(), digiOnly_complete(), PndBarrelTrackFinder::PndBarrelTrackFinder(), PndCATracking::PndCATracking(), PndDrcHitFinder::PndDrcHitFinder(), PndEmc2DLocMaxFinder::PndEmc2DLocMaxFinder(), PndEmcExpClusterSplitter::PndEmcExpClusterSplitter(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcFWEndcapDigi::PndEmcFWEndcapDigi(), PndEmcFWEndcapTimebasedWaveforms::PndEmcFWEndcapTimebasedWaveforms(), PndEmcHitProducer::PndEmcHitProducer(), PndEmcHitsToWaveform::PndEmcHitsToWaveform(), PndEmcMakeBump::PndEmcMakeBump(), PndEmcMakeCluster::PndEmcMakeCluster(), PndEmcMakeClusterOnline::PndEmcMakeClusterOnline(), PndEmcMakeDigi::PndEmcMakeDigi(), PndEmcMakeRecoHit::PndEmcMakeRecoHit(), PndEmcMultiWaveformToCalibratedDigi::PndEmcMultiWaveformToCalibratedDigi(), PndEmcPhiBumpSplitter::PndEmcPhiBumpSplitter(), PndEmcTmpWaveformToDigi::PndEmcTmpWaveformToDigi(), PndEmcWaveformToCalibratedDigi::PndEmcWaveformToCalibratedDigi(), PndEmcWaveformToDigi::PndEmcWaveformToDigi(), PndFtofHitProducerIdeal::PndFtofHitProducerIdeal(), PndFtsCATracking::PndFtsCATracking(), PndFtsHitProducerIdeal::PndFtsHitProducerIdeal(), PndFtsHitProducerRealFast::PndFtsHitProducerRealFast(), PndFtsHitProducerRealFull::PndFtsHitProducerRealFull(), PndFtsHoughTrackerTask::PndFtsHoughTrackerTask(), PndGemDigitize::PndGemDigitize(), PndGemFindHits::PndGemFindHits(), PndIdealTrackFinder(), PndLmdPixelClusterTask::PndLmdPixelClusterTask(), PndLmdPixelHitProducerFast::PndLmdPixelHitProducerFast(), PndMdtClusterTask::PndMdtClusterTask(), PndMdtDigitization::PndMdtDigitization(), PndMdtHitProducerIdeal::PndMdtHitProducerIdeal(), PndMdtPointsToWaveform::PndMdtPointsToWaveform(), PndMdtTrkProducer::PndMdtTrkProducer(), PndMissingPzCleanerTask::PndMissingPzCleanerTask(), PndMvdGemTrackFinderOnHits::PndMvdGemTrackFinderOnHits(), PndMvdHitProducerIdeal::PndMvdHitProducerIdeal(), PndMvdPixelClusterTask::PndMvdPixelClusterTask(), PndMvdTimeWalkCorrTask::PndMvdTimeWalkCorrTask(), PndMvdToPix4ClusterTask::PndMvdToPix4ClusterTask(), PndRecoKalmanTask::PndRecoKalmanTask(), PndRecoKalmanTask2::PndRecoKalmanTask2(), PndRecoMultiKalmanTask::PndRecoMultiKalmanTask(), PndRecoMultiKalmanTask2::PndRecoMultiKalmanTask2(), PndRichHitFinder::PndRichHitFinder(), PndRichHitProducer::PndRichHitProducer(), PndSciTDigiTask::PndSciTDigiTask(), PndSciTHitProducerIdeal::PndSciTHitProducerIdeal(), PndSdsHitProducerIdeal::PndSdsHitProducerIdeal(), PndSdsHybridHitProducer::PndSdsHybridHitProducer(), PndSdsIdealClusterTask::PndSdsIdealClusterTask(), PndSdsIdealRecoTask::PndSdsIdealRecoTask(), PndSdsNoiseProducer::PndSdsNoiseProducer(), PndSdsPixelClusterTask::PndSdsPixelClusterTask(), PndSdsStripClusterTask::PndSdsStripClusterTask(), PndSdsStripHitProducer::PndSdsStripHitProducer(), PndSdsTimeWalkCorrTask::PndSdsTimeWalkCorrTask(), PndSttFindTracks::PndSttFindTracks(), PndSttHelixHitProducer::PndSttHelixHitProducer(), PndSttHitProducerIdeal::PndSttHitProducerIdeal(), PndSttHitProducerRealFast::PndSttHitProducerRealFast(), PndSttHitProducerRealFull::PndSttHitProducerRealFull(), PndSttMatchTracks::PndSttMatchTracks(), PndSttMvdGemTracking::PndSttMvdGemTracking(), PndTrackSmearTask::PndTrackSmearTask(), PndTrkTracking2::PndTrkTracking2(), reco(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), recoideal_complete(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndLmdPixelHitProducerFast::SetPersistance(), PndSdsHitProducerIdeal::SetPersistance(), PndSttMvdGemTracking::SetPersistenc(), PndMdtClusterTask::SetPersistence(), PndSttHelixHitProducer::SetPersistence(), PndMissingPzCleanerTask::SetPersistence(), PndFtsHitProducerRealFast::SetPersistence(), PndFtsHitProducerRealFull::SetPersistence(), PndSttHitProducerIdeal::SetPersistence(), PndSttHitProducerRealFull::SetPersistence(), PndSttHitProducerRealFast::SetPersistence(), PndFtsHitProducerIdeal::SetPersistence(), PndTrackSmearTask::SetPersistence(), PndSciTHitProducerIdeal::SetPersistence(), SetPersistence(), PndSttMatchTracks::SetPersistence(), PndSttFindTracks::SetPersistence(), PndFtsHoughTrackerTask::SetPersistence(), PndTrkTracking2::SetPersistence(), PndEmcMakeRecoHit::SetStorageOfData(), PndEmcMakeClusterOnline::SetStorageOfData(), PndEmcFWEndcapDigi::SetStorageOfData(), PndEmcFWEndcapTimebasedWaveforms::SetStorageOfData(), PndEmcMakeDigi::SetStorageOfData(), PndMdtPointsToWaveform::SetStorageOfData(), PndEmc2DLocMaxFinder::SetStorageOfData(), PndEmcMakeCluster::SetStorageOfData(), PndEmcHitsToWaveform::SetStorageOfData(), PndEmcMakeBump::SetStorageOfData(), PndEmcTmpWaveformToDigi::SetStorageOfData(), PndEmcWaveformToDigi::SetStorageOfData(), PndEmcWaveformToCalibratedDigi::SetStorageOfData(), PndEmcMultiWaveformToCalibratedDigi::SetStorageOfData(), PndEmcExpClusterSplitter::SetStorageOfData(), PndEmcPhiBumpSplitter::SetStorageOfData(), standard_tracking(), and PndEmcFullDigiTask::StoreDigi().

{ fPersistency = val; }

void PndIdealTrackFinder::SetRelativeMomentumSmearing ( Double_t  dpop = -1. )

inline

Definition at line 48 of file PndIdealTrackFinder.h.

References fabs(), fDPoP, and fRelative.

Referenced by PndMasterRecoIdealTask::PndMasterRecoIdealTask(), PndMasterRecoTask::PndMasterRecoTask(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), reco_rich(), and recoideal_complete().

                                                                { 
            fDPoP=fabs(dpop); 
            fRelative=kTRUE;
          }; // in GeV

void PndIdealTrackFinder::SetTrackingEfficiency ( Double_t  eff = 1. )

inline

Definition at line 55 of file PndIdealTrackFinder.h.

References fEfficiency.

Referenced by idealbarreltracking(), idealforwardtracking(), PndMasterRecoIdealTask::PndMasterRecoIdealTask(), PndMasterRecoTask::PndMasterRecoTask(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), reco_rich(), and recoideal_complete().

                                                        { 
            fEfficiency=eff; 
          };

void PndIdealTrackFinder::SetTrackSelector ( TString  selector )

inline

< The track selector decides if a track is taken or discarded. You can set them via a string. Valid strings can be found in PndTrackFunctor.h

Definition at line 59 of file PndIdealTrackFinder.h.

References fTrackSelector, and PndTrackFunctor::make_PndTrackFunctor().

Referenced by idealbarreltracking(), idealforwardtracking(), PndMasterRecoIdealTask::PndMasterRecoIdealTask(), PndMasterRecoTask::PndMasterRecoTask(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), reco_rich(), and recoideal_complete().

                                                 {                      
                  fTrackSelector = PndTrackFunctor::make_PndTrackFunctor(selector.Data());
          }

void PndIdealTrackFinder::SetVertexSmearing ( Double_t  sigmax = -1., Double_t  sigmay = -1., Double_t  sigmaz = -1.  )

inline

Definition at line 52 of file PndIdealTrackFinder.h.

References fabs(), and fVtxSigma.

Referenced by PndMasterRecoIdealTask::PndMasterRecoIdealTask(), PndMasterRecoTask::PndMasterRecoTask(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), reco_rich(), and recoideal_complete().

                                                                                                    { 
            fVtxSigma.SetXYZ(fabs(sigmax),fabs(sigmay),fabs(sigmaz)); 
          }; // in cm

void PndIdealTrackFinder::SmearVector ( TVector3 &  vec, const TVector3 &  sigma  )

protectedvirtual

Definition at line 278 of file PndIdealTrackFinder.cxx.

References Double_t.

Referenced by CreateTracks().

{
  // gaussian smearing
  Double_t rannn=0.;
  rannn = gRandom->Gaus(vec.X(),sigma.X());
  vec.SetX(rannn);
  
  rannn = gRandom->Gaus(vec.Y(),sigma.Y());
  vec.SetY(rannn);
  
  rannn = gRandom->Gaus(vec.Z(),sigma.Z());
  vec.SetZ(rannn);
  
  return;
}

Member Data Documentation

std::map<TString, TClonesArray*> PndIdealTrackFinder::fBranchMap

protected

Definition at line 80 of file PndIdealTrackFinder.h.

Referenced by CreateTrackCands(), CreateTracks(), Init(), and PndCloneGenerator::SetFirstLastHit().

std::vector<TString> PndIdealTrackFinder::fBranchNames

protected

Definition at line 81 of file PndIdealTrackFinder.h.

Referenced by AddBranchName(), and Init().

Double_t PndIdealTrackFinder::fDPoP

protected

Relative momentum Smearing.

Definition at line 93 of file PndIdealTrackFinder.h.

Referenced by CreateTracks(), and SetRelativeMomentumSmearing().

Double_t PndIdealTrackFinder::fEfficiency

protected

Tracking efficiency - if (0 <= e < 1), some tracks will be discarded.

Definition at line 96 of file PndIdealTrackFinder.h.

Referenced by CreateTracks(), and SetTrackingEfficiency().

std::map<FairLink, FairMCPoint > PndIdealTrackFinder::fFirstPointMap

protected

Definition at line 83 of file PndIdealTrackFinder.h.

Referenced by CreateTrackCands(), CreateTracks(), and PndCloneGenerator::SetFirstLastHit().

Int_t PndIdealTrackFinder::fHitCount

protected

Definition at line 89 of file PndIdealTrackFinder.h.

Referenced by CreateTrackCands().

std::map<FairLink, FairMCPoint > PndIdealTrackFinder::fLastPointMap

protected

Definition at line 84 of file PndIdealTrackFinder.h.

Referenced by CreateTrackCands(), CreateTracks(), and PndCloneGenerator::SetFirstLastHit().

TClonesArray* PndIdealTrackFinder::fMCTrack

protected

Definition at line 78 of file PndIdealTrackFinder.h.

Referenced by Init().

TVector3 PndIdealTrackFinder::fMomSigma

protected

Momentum smearing sigma [GeV].

Definition at line 92 of file PndIdealTrackFinder.h.

Referenced by CreateTracks(), and SetMomentumSmearing().

TString PndIdealTrackFinder::fOutBranchName

protected

Definition at line 75 of file PndIdealTrackFinder.h.

Referenced by Init(), and SetOutputBranchName().

TDatabasePDG* PndIdealTrackFinder::fPdg

protected

! Particle DB

Definition at line 87 of file PndIdealTrackFinder.h.

Referenced by CreateTracks(), and Init().

std::map<TString, TString> PndIdealTrackFinder::fPointBranchMap

protected

Definition at line 85 of file PndIdealTrackFinder.h.

Referenced by AddBranchName(), GetFairMCPoint(), Init(), and PndIdealTrackFinder().

Bool_t PndIdealTrackFinder::fRelative

protected

flag

Definition at line 94 of file PndIdealTrackFinder.h.

Referenced by CreateTracks(), SetMomentumSmearing(), and SetRelativeMomentumSmearing().

TClonesArray* PndIdealTrackFinder::fTrack

protected

Definition at line 77 of file PndIdealTrackFinder.h.

Referenced by CreateTracks(), PndCloneGenerator::Exec(), Exec(), and Init().

TClonesArray* PndIdealTrackFinder::fTrackCand

protected

Definition at line 76 of file PndIdealTrackFinder.h.

Referenced by CreateTracks(), PndCloneGenerator::Exec(), Exec(), and Init().

std::map<FairLink, PndTrackCand> PndIdealTrackFinder::fTrackCandMap

protected

Definition at line 82 of file PndIdealTrackFinder.h.

Referenced by PndCloneGenerator::CreateTrackCandClones(), CreateTrackCands(), CreateTracks(), PndCloneGenerator::Exec(), Exec(), and FilterTrackCands().

PndTrackFunctor* PndIdealTrackFinder::fTrackSelector

protected

Definition at line 79 of file PndIdealTrackFinder.h.

Referenced by FilterTrackCands(), Init(), SetTrackSelector(), and ~PndIdealTrackFinder().

TVector3 PndIdealTrackFinder::fVtxSigma

protected

Vertex smearing sigma [cm].

Definition at line 95 of file PndIdealTrackFinder.h.

Referenced by CreateTracks(), and SetVertexSmearing().

---

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

• PndIdealTrackFinder.h
• PndIdealTrackFinder.cxx