PndMvdRiemannTrackFinderTaskEff Class Reference

#include <PndMvdRiemannTrackFinderTaskEff.h>

Inheritance diagram for PndMvdRiemannTrackFinderTaskEff:

Public Member Functions
	PndMvdRiemannTrackFinderTaskEff ()
virtual	~PndMvdRiemannTrackFinderTaskEff ()
	PndMvdRiemannTrackFinderTaskEff (const PndMvdRiemannTrackFinderTaskEff &)=delete
PndMvdRiemannTrackFinderTaskEff &	operator= (const PndMvdRiemannTrackFinderTaskEff &)=delete
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
virtual void	Exec (Option_t *opt)
virtual void	FinishEvent ()
void	SetMaxSZChi2 (double val)
void	SetMaxSZDist (double val)
void	SetMinPointDist (double val)
void	SetMaxDist (double val)
void	SetVerbose (Int_t verbose)

Public Attributes
TH2F *	eff0H
TH2F *	effH
TH2F *	GhH

Private Member Functions
void	Register ()
void	Reset ()
void	ProduceHits ()
void	ComparingFandR (std::vector< PndTrackCand * > RecoT)
bool	CheckRecoTrack (PndTrackCand *cand, PndMCTrack *myTrack)
void	AddGhostTrack (int trackF)

Private Attributes
TString	fHitBranch
TString	fHitBranch2
TString	fMCTrackBranch
TString	fIdealTrackBranch
TString	fFTrackBranch
int	fEventNr
double	fMaxSZChi2
double	fMaxSZDist
double	fMinPointDist
double	fMaxDist
TClonesArray *	fHitArray
TClonesArray *	fHitArray2
TClonesArray *	fTrackCandArray
TClonesArray *	fIdealTrackCandArray
TClonesArray *	fMCTracksArray

Detailed Description

Definition at line 12 of file PndMvdRiemannTrackFinderTaskEff.h.

Constructor & Destructor Documentation

PndMvdRiemannTrackFinderTaskEff::PndMvdRiemannTrackFinderTaskEff

(

)

Definition at line 32 of file PndMvdRiemannTrackFinderTaskEff.cxx.

                                                                 : FairTask("MVD Riemann Track Finder"),
    eff0H(NULL),
    effH(NULL),
    GhH(NULL),
          fHitBranch("MVDHitsPixel"),
          fHitBranch2("MVDHitsStrip"),
    fMCTrackBranch("MCTrack"),
    fIdealTrackBranch("MVDIdealTrackCand"),
    fFTrackBranch("MVDRiemannTrackCand"),
    fEventNr(0),
    fMaxSZChi2(1),
    fMaxSZDist(10),
    fMinPointDist(1),
    fMaxDist(1),
          fHitArray(NULL),
          fHitArray2(NULL),
          fTrackCandArray(NULL),
          fIdealTrackCandArray(NULL),
          fMCTracksArray(NULL)
{
}

PndMvdRiemannTrackFinderTaskEff::~PndMvdRiemannTrackFinderTaskEff ( )

virtual

Definition at line 54 of file PndMvdRiemannTrackFinderTaskEff.cxx.

{
}

PndMvdRiemannTrackFinderTaskEff::PndMvdRiemannTrackFinderTaskEff ( const PndMvdRiemannTrackFinderTaskEff &  )

delete

Member Function Documentation

void PndMvdRiemannTrackFinderTaskEff::AddGhostTrack ( int  trackF )

private

Definition at line 231 of file PndMvdRiemannTrackFinderTaskEff.cxx.

References PndRiemannTrack::addHit(), CAMath::Cos(), PndRiemannTrack::dip(), Double_t, fHitArray, fHitArray2, fHitBranch, fHitBranch2, fTrackCandArray, GhH, hit, i, P, Pi, point, R, PndRiemannTrack::r(), PndRiemannTrack::refit(), PndRiemannHit::setDXYZ(), PndRiemannHit::setXYZ(), sign(), PndRiemannTrack::szFit(), CAMath::Tan(), and track.

Referenced by ComparingFandR().

{
        PndRiemannTrack track;
        unsigned int detID, hitID;
        bool sign=true;
        for(UInt_t i=0;i<((PndTrackCand*)fTrackCandArray->At(trackF))->GetNHits();i++){
                        detID=((PndTrackCand*)fTrackCandArray->At(trackF))->GetSortedHit(i).GetDetId();
                        hitID=((PndTrackCand*)fTrackCandArray->At(trackF))->GetSortedHit(i).GetHitId();
                        PndSdsHit *point;
                        if ((int)detID == FairRootManager::Instance()->GetBranchId(fHitBranch))
                                 point = (PndSdsHit*)fHitArray->At(hitID);
                        else if ((int)detID == FairRootManager::Instance()->GetBranchId(fHitBranch2))
                                 point = (PndSdsHit*)fHitArray2->At(hitID);
                        else point = 0;
                        if (point!=0){
                                PndRiemannHit hit;
                                hit.setXYZ(point->GetX(),point->GetY(),point->GetZ());
                                hit.setDXYZ(point->GetDx(),point->GetDy(),point->GetDz());
                                track.addHit(hit);
                        }
                        if (i==1){
                                if (point->GetZ()>0){
                                        sign=true;
                                }
                                else sign=false;
                        }
        }
        track.refit();
        track.szFit();
        Double_t R=track.r();
        Double_t Pt=R*(2*3*1E8)/(1E9*100);
        Double_t dip=track.dip();
        Double_t Theta;
         if (sign) Theta=(TMath::ATan(TMath::Power(TMath::Tan(TMath::ACos(dip)),-1)))*180/TMath::Pi();
     else Theta=(TMath::Pi()-TMath::ATan(TMath::Power(TMath::Tan(TMath::ACos(dip)),-1)))*180/TMath::Pi();
        Double_t P=Pt/TMath::Cos(Theta);
        GhH->Fill(P,Theta);

}

bool PndMvdRiemannTrackFinderTaskEff::CheckRecoTrack ( PndTrackCand *  cand, PndMCTrack *  myTrack  )

private

Definition at line 147 of file PndMvdRiemannTrackFinderTaskEff.cxx.

References a, count, fHitArray, fHitArray2, fHitBranch, fHitBranch2, fMinPointDist, PndTrackCandHit::GetDetId(), PndTrackCandHit::GetHitId(), PndTrackCand::GetNHits(), PndSdsHit::GetPosition(), PndTrackCand::GetSortedHit(), and i.

Referenced by Exec().

{
        if (/*myTrack->GetPdgCode()==211 &&*/ (cand->GetNHits()>2)/* && (myTrack->GetMotherID()==-1)*/){
                int count=0;
                unsigned int detIDi, hitIDi;
                unsigned int detIDj, hitIDj;
                for(unsigned int i=0;i<cand->GetNHits();i++){
                        detIDi=cand->GetSortedHit(i).GetDetId();
                        hitIDi=cand->GetSortedHit(i).GetHitId();
                        PndSdsHit *pointI;
                        if ((int)detIDi == FairRootManager::Instance()->GetBranchId(fHitBranch))
                                 pointI = (PndSdsHit*)fHitArray->At(hitIDi);
                        else if ((int)detIDi == FairRootManager::Instance()->GetBranchId(fHitBranch2))
                                 pointI = (PndSdsHit*)fHitArray2->At(hitIDi);
                        else pointI = 0;

                                for(unsigned int j=0;j<cand->GetNHits();j++){
                                        detIDj=cand->GetSortedHit(j).GetDetId();
                                        hitIDj=cand->GetSortedHit(j).GetHitId();
                                        PndSdsHit *pointJ;
                                        if ((int)detIDj == FairRootManager::Instance()->GetBranchId(fHitBranch))
                                                 pointJ = (PndSdsHit*)fHitArray->At(hitIDj);
                                        else if ((int)detIDj == FairRootManager::Instance()->GetBranchId(fHitBranch2))
                                                 pointJ = (PndSdsHit*)fHitArray2->At(hitIDj);
                                        else pointJ = 0;

                                        if ((pointI!=0) && (pointJ!=0) && (i!=j)){
                                                TVector3 a=pointI->GetPosition() - pointJ->GetPosition();
                                                if (a.Mag()<fMinPointDist){
                                                        count++;
                                                }
                                        }
                                }
                }
                if ((cand->GetNHits()-count/2)<3){
                        std::cout <<"Less then 3 base points in RecoTrack"<<"   "<<cand->GetNHits()<<"   "<<count <<std::endl;
                        return false;
                }
                else return true;
        }
        else
                return false;
}

void PndMvdRiemannTrackFinderTaskEff::ComparingFandR ( std::vector< PndTrackCand * >  RecoT )

private

Definition at line 190 of file PndMvdRiemannTrackFinderTaskEff.cxx.

References AddGhostTrack(), effH, fMCTracksArray, fTrackCandArray, PndMCTrack::GetMomentum(), P, and Pi.

Referenced by Exec().

{
          unsigned int detidRC, hitidRC;
          unsigned int detidF, hitidF;
                for(int trackF=0;trackF<fTrackCandArray->GetEntriesFast();trackF++){
                        for(unsigned int trackRC=0;trackRC<RecoT.size();trackRC++){
                                int cSame=0;
                          for(unsigned int iF=0;iF<((PndTrackCand*)fTrackCandArray->At(trackF))->GetNHits();iF++){
                                  detidF=((PndTrackCand*)fTrackCandArray->At(trackF))->GetSortedHit(iF).GetDetId();
                                  hitidF=((PndTrackCand*)fTrackCandArray->At(trackF))->GetSortedHit(iF).GetHitId();
                                        for(unsigned int iRC=0;iRC<RecoT[trackRC]->GetNHits();iRC++){
                                                detidRC=RecoT[trackRC]->GetSortedHit(iRC).GetDetId();
                                                hitidRC=RecoT[trackRC]->GetSortedHit(iRC).GetHitId();
                                                if (detidRC==detidF && hitidRC==hitidF){
                                                        cSame++;
                                                }
                                        }
                                }

                              PndMCTrack* myTrack = (PndMCTrack*)fMCTracksArray->At(RecoT[trackRC]->getMcTrackId());
                                  if (  (cSame>=4) && ((int)(((PndTrackCand*)fTrackCandArray->At(trackF))->GetNHits())==cSame)){
                                          std::cout<<"TRUE FOUND"<<std::endl;
                                             TVector3 Pvec=myTrack->GetMomentum();
                                                 double Theta=180.0*Pvec.Theta()/TMath::Pi();
                                                 double P=Pvec.Mag();
                                                 effH->Fill(P,Theta);
                                                 RecoT.erase(RecoT.begin()+trackRC);
                                             break;
                                  }
                                  else{
                                          if ((trackRC)==RecoT.size()-1){
                                        std::cout<<"        FALSE FOUND"<<std::endl;
                                        AddGhostTrack(trackF);
                                          }
                                  }

                         }
                }
}

void PndMvdRiemannTrackFinderTaskEff::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 122 of file PndMvdRiemannTrackFinderTaskEff.cxx.

References PndTrackCand::AddHit(), CheckRecoTrack(), ComparingFandR(), eff0H, fEventNr, fIdealTrackCandArray, fMCTracksArray, PndTrackCand::getMcTrackId(), PndMCTrack::GetMomentum(), PndMCTrack::GetMotherID(), i, P, and Pi.

{

  std::vector<PndTrackCand*>  RecoT;

          for(int i=0;i<fIdealTrackCandArray->GetEntriesFast();i++){
                 PndTrackCand *cand=(PndTrackCand*)fIdealTrackCandArray->At(i);
                 PndMCTrack* myTrack = (PndMCTrack*)fMCTracksArray->At(cand->getMcTrackId());
                 if (CheckRecoTrack(cand,myTrack)){
                         if (myTrack->GetMotherID()==-1) cand->AddHit(0,0,0);
                         RecoT.push_back(cand);
                         TVector3 Pvec=myTrack->GetMomentum();
                         double Theta=180.0*Pvec.Theta()/TMath::Pi();
                         double P=Pvec.Mag();
                         eff0H->Fill(P,Theta);
             }
          }
          ComparingFandR(RecoT);
          std::cout<<"                             Event n: "<<fEventNr++<<std::endl;

}

void PndMvdRiemannTrackFinderTaskEff::FinishEvent ( )

virtual

Definition at line 144 of file PndMvdRiemannTrackFinderTaskEff.cxx.

{
}

InitStatus PndMvdRiemannTrackFinderTaskEff::Init ( )

virtual

Definition at line 69 of file PndMvdRiemannTrackFinderTaskEff.cxx.

References eff0H, effH, fFTrackBranch, fHitArray, fHitArray2, fHitBranch, fHitBranch2, fIdealTrackBranch, fIdealTrackCandArray, fMCTrackBranch, fMCTracksArray, fTrackCandArray, and GhH.

{

  FairRootManager* ioman = FairRootManager::Instance();

  if ( ! ioman )
    {
      std::cout << "-E- PndMvdRiemannTrackFinderTask::Init: "
     << "RootManager not instantiated!" << std::endl;
      return kFATAL;
    }

  // Get input array
  fHitArray = (TClonesArray*) ioman->GetObject(fHitBranch);
  if ( !fHitArray){
    std::cout << "-W- PndMvdRiemannTrackFinderTask::Init: " << "No hitArray!" << std::endl;
    return kERROR;
  }

  fHitArray2 = (TClonesArray*) ioman->GetObject(fHitBranch2);
   if ( !fHitArray2){
     std::cout << "-W- PndMvdRiemannTrackFinderTask::Init: " << "No hitArray2!" << std::endl;
     return kERROR;
   }

   fMCTracksArray = (TClonesArray*) ioman->GetObject(fMCTrackBranch);
   if ( !fMCTracksArray){
        std::cout << "-W- PndMvdRiemannTrackFinderTask::Init: " << "No MCtracks!" << std::endl;
        return kERROR;
   }

  fIdealTrackCandArray = (TClonesArray*) ioman->GetObject(fIdealTrackBranch);
  if ( !fIdealTrackCandArray){
     std::cout << "-W- PndMvdRiemannTrackFinderTask::Init: " << "No IdealTrackCands!" << std::endl;
     return kERROR;
  }
  fTrackCandArray = (TClonesArray*) ioman->GetObject(fFTrackBranch);
  if ( !fTrackCandArray){
     std::cout << "-W- PndMvdRiemannTrackFinderTask::Init: " << "No FTrackCands!" << std::endl;
     return kERROR;
  }

    int Nbin=10;
    eff0H = new TH2F("eff0H","eff0H",Nbin,0.1,1,Nbin,15,150);
    effH = new TH2F("effH","effH",Nbin,0.1,1,Nbin,15,150);

    GhH = new TH2F("GhH","GhH",Nbin,0.1,1,Nbin,15,150);

  std::cout << "-I- PndMvdRiemannTrackFinderTask: Initialisation successfull" << std::endl;
  return kSUCCESS;
}

PndMvdRiemannTrackFinderTaskEff& PndMvdRiemannTrackFinderTaskEff::operator= ( const PndMvdRiemannTrackFinderTaskEff &  )

delete

void PndMvdRiemannTrackFinderTaskEff::ProduceHits ( )

private

void PndMvdRiemannTrackFinderTaskEff::Register ( )

private

InitStatus PndMvdRiemannTrackFinderTaskEff::ReInit ( )

virtual

Definition at line 62 of file PndMvdRiemannTrackFinderTaskEff.cxx.

{
  InitStatus stat=kSUCCESS;
  return stat;
}

void PndMvdRiemannTrackFinderTaskEff::Reset ( )

private

void PndMvdRiemannTrackFinderTaskEff::SetMaxDist ( double  val )

inline

Definition at line 32 of file PndMvdRiemannTrackFinderTaskEff.h.

References fMaxDist, and val.

{fMaxDist = val;}

void PndMvdRiemannTrackFinderTaskEff::SetMaxSZChi2 ( double  val )

inline

Definition at line 29 of file PndMvdRiemannTrackFinderTaskEff.h.

References fMaxSZChi2, and val.

{fMaxSZChi2 = val;}

void PndMvdRiemannTrackFinderTaskEff::SetMaxSZDist ( double  val )

inline

Definition at line 30 of file PndMvdRiemannTrackFinderTaskEff.h.

References fMaxSZDist, and val.

{fMaxSZDist = val;}

void PndMvdRiemannTrackFinderTaskEff::SetMinPointDist ( double  val )

inline

Definition at line 31 of file PndMvdRiemannTrackFinderTaskEff.h.

References fMinPointDist, and val.

{fMinPointDist = val;}

void PndMvdRiemannTrackFinderTaskEff::SetParContainers ( )

virtual

Virtual method Init

Definition at line 58 of file PndMvdRiemannTrackFinderTaskEff.cxx.

{
}

void PndMvdRiemannTrackFinderTaskEff::SetVerbose ( Int_t  verbose )

inline

Definition at line 35 of file PndMvdRiemannTrackFinderTaskEff.h.

References fVerbose, and verbose.

{ fVerbose = verbose;};

Member Data Documentation

TH2F* PndMvdRiemannTrackFinderTaskEff::eff0H

Definition at line 35 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by Exec(), and Init().

TH2F* PndMvdRiemannTrackFinderTaskEff::effH

Definition at line 38 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by ComparingFandR(), and Init().

int PndMvdRiemannTrackFinderTaskEff::fEventNr

private

Definition at line 49 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by Exec().

TString PndMvdRiemannTrackFinderTaskEff::fFTrackBranch

private

Definition at line 47 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by Init().

TClonesArray* PndMvdRiemannTrackFinderTaskEff::fHitArray

private

Definition at line 56 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by AddGhostTrack(), CheckRecoTrack(), and Init().

TClonesArray* PndMvdRiemannTrackFinderTaskEff::fHitArray2

private

Definition at line 57 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by AddGhostTrack(), CheckRecoTrack(), and Init().

TString PndMvdRiemannTrackFinderTaskEff::fHitBranch

private

Definition at line 42 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by AddGhostTrack(), CheckRecoTrack(), and Init().

TString PndMvdRiemannTrackFinderTaskEff::fHitBranch2

private

Definition at line 43 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by AddGhostTrack(), CheckRecoTrack(), and Init().

TString PndMvdRiemannTrackFinderTaskEff::fIdealTrackBranch

private

Definition at line 46 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by Init().

TClonesArray* PndMvdRiemannTrackFinderTaskEff::fIdealTrackCandArray

private

Definition at line 59 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by Exec(), and Init().

double PndMvdRiemannTrackFinderTaskEff::fMaxDist

private

Definition at line 54 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by SetMaxDist().

double PndMvdRiemannTrackFinderTaskEff::fMaxSZChi2

private

Definition at line 51 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by SetMaxSZChi2().

double PndMvdRiemannTrackFinderTaskEff::fMaxSZDist

private

Definition at line 52 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by SetMaxSZDist().

TString PndMvdRiemannTrackFinderTaskEff::fMCTrackBranch

private

Definition at line 44 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by Init().

TClonesArray* PndMvdRiemannTrackFinderTaskEff::fMCTracksArray

private

Definition at line 60 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by ComparingFandR(), Exec(), and Init().

double PndMvdRiemannTrackFinderTaskEff::fMinPointDist

private

Definition at line 53 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by CheckRecoTrack(), and SetMinPointDist().

TClonesArray* PndMvdRiemannTrackFinderTaskEff::fTrackCandArray

private

Definition at line 58 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by AddGhostTrack(), ComparingFandR(), and Init().

TH2F* PndMvdRiemannTrackFinderTaskEff::GhH

Definition at line 39 of file PndMvdRiemannTrackFinderTaskEff.h.

Referenced by AddGhostTrack(), and Init().

---

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

• PndMvdRiemannTrackFinderTaskEff.h
• PndMvdRiemannTrackFinderTaskEff.cxx