PndBarrelTrackFinder Class Reference

#include <PndBarrelTrackFinder.h>

Inheritance diagram for PndBarrelTrackFinder:

Public Member Functions
	PndBarrelTrackFinder ()
	PndBarrelTrackFinder (Int_t iVerbose)
	PndBarrelTrackFinder (const char *name, Int_t iVerbose)
virtual	~PndBarrelTrackFinder ()
virtual void	Exec (Option_t *opt)
void	UseMvdSttGem (const Bool_t useMvd, const Bool_t useStt, const Bool_t useGem)
void	DrawDetails (Bool_t tb=kTRUE)
void	DrawHistos (Bool_t tb=kTRUE)
void	SetPersistency (Bool_t val=kTRUE)
Bool_t	GetPersistency ()

Private Member Functions
Bool_t	MatchHitToTrack (FairHit *thisHit, Int_t detId, Int_t hitNo, Int_t trackNo)
Bool_t	MatchSkewedSttHitTT (FairHit *thisHit, Int_t detId, Int_t hitNo, Int_t trackNo)
Bool_t	MatchParallelSttHitTT (FairHit *thisHit, Int_t detId, Int_t hitNo, Int_t trackNo)
Bool_t	MatchHitWithZInfoTT (FairHit *thisHit, Int_t detId, Int_t hitNo, Int_t trackNo)
Bool_t	ExtractMeanRPhiFromTrack (Int_t trackNo)
Bool_t	ExtractMeanZ_PFromTrack (Int_t trackNo)
Bool_t	MatchHitToHit (FairHit *thisHit, Int_t detId, Int_t hitNo, Int_t prevHNo)
void	RemoveHitFromPreviousHits (Int_t hitNo)
Bool_t	HitBelongsToTrack (Int_t detId, Int_t hitNo, Int_t trackNo)
void	AddHitToPreviousHits (FairHit *thisHit, Int_t detId, Int_t hitNo)
void	DrawTracks ()
void	PrintTracks ()
Int_t	CleanTracks ()
Int_t	WriteTracks ()
void	RemoveShortTracks ()
void	RemoveCloneTracks ()
Bool_t	FindCircPar (Double_t *c1, Double_t *c2, Double_t *c3, Int_t cno, Double_t *cl)
Double_t	FindCircDist (Double_t *c1, Double_t *c2)
Int_t	FindInterestingRegions (Double_t *circ, Double_t *tube, Double_t *reg)
Double_t	CalcPhi (Double_t x, Double_t y)
Double_t	CalcZ_P (Double_t *circ, Double_t hx, Double_t hy, Double_t hz)
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
void	Reset ()
virtual void	Finish ()
	ClassDef (PndBarrelTrackFinder, 1)

Private Attributes
Bool_t	fDrawDetails
Bool_t	fDrawHistos
TEllipse *	fEllipse [1000]
TCanvas *	fCanvas
TLatex *	fStageDesc
TH1F *	fClonesXDiff
TH1F *	fClonesYDiff
TH1F *	fClonesRDiff
TH1F *	fClonesZ_PDiff
TH1F *	fClonesXDiffRel
TH1F *	fClonesYDiffRel
TH1F *	fClonesRDiffRel
TH1F *	fClonesZ_PDiffRel
std::vector< TrackBasis >	fTracksVector
std::vector< Int_t >	fHitDetId
std::vector< Int_t >	fHitDetNo
std::vector< FairHit * >	fHitVector
std::vector< Int_t >	fHitVectDI
std::vector< Int_t >	fHitVectHN
Bool_t	fIncludeDet [4]
TClonesArray *	fHitArray [4]
TString	fDetName [4]
Int_t	fDetType [4]
TString	fHitArrayName [4]
TClonesArray *	fBarrelTrackArray
	Output array of PndTrackCands. More...
TClonesArray *	fBarrelTrackCandArray
Int_t	fTNofEvents
	event counter More...
Int_t	fTNofTracks
Double_t	fMaximalDist
Double_t	fReasonableDist
Double_t	fMaximalZ_PD
Double_t	fReasonableZ_PD
Double_t	fMaximalRadDiff
Double_t	fMaximalPhiDiff
TClonesArray *	fTubeArray
PndGeoSttPar *	fSttParameters
Int_t	fGeoType

Detailed Description

Definition at line 65 of file PndBarrelTrackFinder.h.

Constructor & Destructor Documentation

PndBarrelTrackFinder::PndBarrelTrackFinder

(

)

Default constructor

Definition at line 67 of file PndBarrelTrackFinder.cxx.

References fBarrelTrackArray, fBarrelTrackCandArray, fDetName, fDetType, fDrawDetails, fDrawHistos, fHitArray, fHitArrayName, fHitDetId, fHitDetNo, fHitVectDI, fHitVectHN, fHitVector, fIncludeDet, fMaximalDist, fMaximalPhiDiff, fMaximalRadDiff, fMaximalZ_PD, fReasonableDist, fReasonableZ_PD, fSttParameters, fTNofEvents, fTNofTracks, fTracksVector, fTubeArray, and Reset().

                                           : PndPersistencyTask("Barrel Track Finder", 1), fGeoType(0) {
  for ( Int_t idet = 0 ; idet < 4 ; idet++ ) {
    fIncludeDet  [idet] = kFALSE;
    fHitArray    [idet] = NULL;
    fDetName     [idet] = "";
    fDetType     [idet] = 0;
    fHitArrayName[idet] = "";
  }
  fBarrelTrackArray          = NULL;
  fBarrelTrackCandArray      = NULL;

  fTNofEvents    = 0;
  fTNofTracks    = 0;
  fTracksVector.clear();

  fHitDetId.clear();
  fHitDetNo.clear();
  fHitVector.clear();
  fHitVectDI.clear();
  fHitVectHN.clear();

  fMaximalDist    = 0.01;
  fReasonableDist = 0.2;

  fMaximalZ_PD    = 2.;
  fReasonableZ_PD = 5.;

  fMaximalRadDiff = 0.1;
  fMaximalPhiDiff = 0.1;

  fTubeArray = NULL;
  fSttParameters = NULL;

  fDrawDetails = kFALSE;
  fDrawHistos  = kFALSE;
 
  Reset();
}

PndBarrelTrackFinder::PndBarrelTrackFinder

(

Int_t

iVerbose

)

Standard constructor

Definition at line 110 of file PndBarrelTrackFinder.cxx.

References fBarrelTrackArray, fBarrelTrackCandArray, fDetName, fDetType, fHitArray, fHitArrayName, fHitDetId, fHitDetNo, fHitVectDI, fHitVectHN, fHitVector, fIncludeDet, fMaximalDist, fMaximalPhiDiff, fMaximalRadDiff, fMaximalZ_PD, fReasonableDist, fReasonableZ_PD, fSttParameters, fTNofEvents, fTNofTracks, fTracksVector, fTubeArray, and Reset().

  : PndPersistencyTask("Barrel Track Finder", iVerbose) {
  for ( Int_t idet = 0 ; idet < 4 ; idet++ ) {
    fIncludeDet[idet] = kFALSE;
    fHitArray  [idet] = NULL;
    fDetName     [idet] = "";
    fDetType     [idet] = 0;
    fHitArrayName[idet] = "";
  }
  fBarrelTrackArray          = NULL;
  fBarrelTrackCandArray      = NULL;

  fTNofEvents    = 0;
  fTNofTracks    = 0;
  fTracksVector.clear();

  fHitDetId.clear();
  fHitDetNo.clear();
  fHitVector.clear();
  fHitVectDI.clear();
  fHitVectHN.clear();

  fMaximalDist    = 0.01;
  fReasonableDist = 0.2;

  fMaximalZ_PD    = 2.;
  fReasonableZ_PD = 5.;

  fMaximalRadDiff = 0.1;
  fMaximalPhiDiff = 0.1;

  fTubeArray = NULL;
  fSttParameters = NULL;
 
  Reset();
}

PndBarrelTrackFinder::PndBarrelTrackFinder	(	const char *	name,
		Int_t	iVerbose
	)

Constructor with name

Definition at line 151 of file PndBarrelTrackFinder.cxx.

References fBarrelTrackArray, fBarrelTrackCandArray, fDetName, fDetType, fHitArray, fHitArrayName, fHitDetId, fHitDetNo, fHitVectDI, fHitVectHN, fHitVector, fIncludeDet, fMaximalDist, fMaximalPhiDiff, fMaximalRadDiff, fMaximalZ_PD, fReasonableDist, fReasonableZ_PD, fSttParameters, fTNofEvents, fTNofTracks, fTracksVector, fTubeArray, Reset(), and PndPersistencyTask::SetPersistency().

  : PndPersistencyTask(name, iVerbose) {
  for ( Int_t idet = 0 ; idet < 4 ; idet++ ) {
    fIncludeDet[idet] = kFALSE;
    fHitArray  [idet] = NULL;
    fDetName     [idet] = "";
    fDetType     [idet] = 0;
    fHitArrayName[idet] = "";
  }
  fBarrelTrackArray          = NULL;
  fBarrelTrackCandArray      = NULL;

  fTNofEvents    = 0;
  fTNofTracks    = 0;
  fTracksVector.clear();

  fHitDetId.clear();
  fHitDetNo.clear();
  fHitVector.clear();
  fHitVectDI.clear();
  fHitVectHN.clear();

  fMaximalDist    = 0.01;
  fReasonableDist = 0.2;
  fMaximalZ_PD    = 2.;
  fReasonableZ_PD = 5.;

  fMaximalRadDiff = 0.1;
  fMaximalPhiDiff = 0.1;

  fTubeArray = NULL;
  fSttParameters = NULL;
 
  SetPersistency(kTRUE);
  Reset();
}

PndBarrelTrackFinder::~PndBarrelTrackFinder ( )

virtual

Destructor

Definition at line 192 of file PndBarrelTrackFinder.cxx.

References fBarrelTrackArray, fBarrelTrackCandArray, and Reset().

                                            { 
  fBarrelTrackArray    ->Delete();
  fBarrelTrackCandArray->Delete();
  Reset();
}

Member Function Documentation

void PndBarrelTrackFinder::AddHitToPreviousHits ( FairHit *  thisHit, Int_t  detId, Int_t  hitNo  )

private

Definition at line 1567 of file PndBarrelTrackFinder.cxx.

References fHitVectDI, fHitVectHN, fHitVector, and fVerbose.

Referenced by Exec(), and RemoveShortTracks().

                                                                                               {
  if ( fVerbose > 3 )
    cout << "!!! ADDING HIT " << detId << "." << hitNo << endl;
  fHitVector.push_back(thisHit);
  fHitVectDI.push_back(detId);
  fHitVectHN.push_back(hitNo);
}

Double_t PndBarrelTrackFinder::CalcPhi ( Double_t  x, Double_t  y  )

private

Definition at line 1687 of file PndBarrelTrackFinder.cxx.

References Double_t, p, Pi, r, and CAMath::Sqrt().

Referenced by CalcZ_P(), ExtractMeanRPhiFromTrack(), ExtractMeanZ_PFromTrack(), MatchHitToHit(), MatchHitWithZInfoTT(), MatchParallelSttHitTT(), MatchSkewedSttHitTT(), PrintTracks(), RemoveCloneTracks(), and WriteTracks().

                                                             {
  Double_t r = TMath::Sqrt(x*x+y*y);
  Double_t p = TMath::ACos(x/r);
  if ( y < 0 )
    p = 2.*TMath::Pi() - p;
  return p;
}

Double_t PndBarrelTrackFinder::CalcZ_P ( Double_t *  circ, Double_t  hx, Double_t  hy, Double_t  hz  )

private

Definition at line 1697 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), CalcPhi(), CAMath::Cos(), Double_t, Pi, CAMath::Sin(), tmpx, and tmpy.

Referenced by ExtractMeanZ_PFromTrack(), MatchHitToHit(), MatchHitWithZInfoTT(), and MatchParallelSttHitTT().

                                                                                            {
  Double_t tmpx = hx-circ[0];
  Double_t tmpy = hy-circ[1];
  Double_t rotA = TMath::Pi()-CalcPhi(circ[0],circ[1]);
  Double_t tmpx2 = tmpx*TMath::Cos(rotA)-tmpy*TMath::Sin(rotA);
  Double_t tmpy2 = tmpx*TMath::Sin(rotA)+tmpy*TMath::Cos(rotA);
  Double_t phiT = CalcPhi(tmpx2,tmpy2);
  if ( phiT > TMath::Pi() ) phiT = -(2.*TMath::Pi()-phiT);
  return hz/TMath::Abs(phiT);
  //Double_t phiT = CalcPhi(hx-circ[0],hy-circ[1]);
  //  return phiT/hz;
}

PndBarrelTrackFinder::ClassDef ( PndBarrelTrackFinder  , 1    )

private

Int_t PndBarrelTrackFinder::CleanTracks ( )

private

Definition at line 1372 of file PndBarrelTrackFinder.cxx.

                                         {
  return -1;
}

void PndBarrelTrackFinder::DrawDetails ( Bool_t  tb = kTRUE )

inline

Definition at line 91 of file PndBarrelTrackFinder.h.

References fDrawDetails.

{fDrawDetails=tb;}

void PndBarrelTrackFinder::DrawHistos ( Bool_t  tb = kTRUE )

inline

Definition at line 92 of file PndBarrelTrackFinder.h.

References fDrawHistos.

{fDrawHistos =tb;}

void PndBarrelTrackFinder::DrawTracks ( )

private

Definition at line 1327 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), fEllipse, fGeoType, fTracksVector, fTubeArray, PndSttTubeMap::GetTube(), PndSttTube::GetWireDirection(), and PndSttTubeMap::Instance().

Referenced by Exec().

                                        {
  for ( Int_t iell = 0 ; iell < 1000 ; iell++ ) {
    fEllipse[iell]->SetX1(100);
    fEllipse[iell]->SetY1(100);
    fEllipse[iell]->SetR1(0);
    fEllipse[iell]->SetR2(0);
    fEllipse[iell]->Paint();
  }

  for ( size_t itr = 0 ; itr < fTracksVector.size() ; itr++ ) {
    if ( TMath::Abs(fTracksVector[itr].meanR) > 0.1 ) {
      for ( size_t ihit = 0 ; ihit < fTracksVector[itr].trackHitD.size() ; ihit++ ) {
        //Bool_t skewedSttHit = kFALSE; //[R.K. 01/2017] unused variable?
        if ( fTracksVector[itr].trackHitD[ihit] == 2 ) {
          Int_t iTube = ((PndSttHit*) fTracksVector[itr].trackHits[ihit])->GetTubeID(); 
            PndSttTube *sttTube = nullptr;
            if (fGeoType == 1)
                sttTube = (PndSttTube*) fTubeArray->At(iTube);
            else if (fGeoType == 2){
                sttTube = PndSttTubeMap::Instance()->GetTube(iTube);
            }
          if ( sttTube->GetWireDirection().Z() < 1. ) {
            //skewedSttHit = kTRUE; //[R.K. 01/2017] unused variable?
          }
        }
      }
      
      //Double_t thisPhi = CalcPhi(fTracksVector[itr].meanX,
                                 //fTracksVector[itr].meanY); //[R.K. 01/2017] unused variable?
      //Double_t thisRad =         fTracksVector[itr].meanR; //[R.K. 01/2017] unused variable?
      
      fEllipse[itr]->SetX1(fTracksVector[itr].meanX);
      fEllipse[itr]->SetY1(fTracksVector[itr].meanY);
      fEllipse[itr]->SetR1(fTracksVector[itr].meanR);
      fEllipse[itr]->SetR2(fTracksVector[itr].meanR);
      fEllipse[itr]->SetFillStyle(0);
      fEllipse[itr]->SetLineWidth(2);
      fEllipse[itr]->SetLineColor(51+2*itr);
      fEllipse[itr]->Paint(); 
    }
  }
}

void PndBarrelTrackFinder::Exec ( Option_t *  opt )

virtual

Execution

Definition at line 209 of file PndBarrelTrackFinder.cxx.

References AddHitToPreviousHits(), At, Bool_t, Double_t, DrawTracks(), fBarrelTrackArray, fBarrelTrackCandArray, fCanvas, fDetName, fEllipse, fGeoType, fHitArray, fHitDetId, fHitDetNo, fHitVectDI, fHitVectHN, fHitVector, fIncludeDet, fStageDesc, fTNofEvents, fTNofTracks, fTracksVector, fTubeArray, fVerbose, PndSttTubeMap::GetTube(), PndSttTube::GetWireDirection(), HitBelongsToTrack(), PndSttTubeMap::Instance(), MatchHitToHit(), MatchHitWithZInfoTT(), MatchParallelSttHitTT(), MatchSkewedSttHitTT(), PrintTracks(), RemoveCloneTracks(), RemoveHitFromPreviousHits(), RemoveShortTracks(), Reset(), and WriteTracks().

                                         {
  Reset();
  fBarrelTrackArray    ->Delete();
  fBarrelTrackCandArray->Delete();

  if ( fVerbose > 0 ) 
    cout << "=============== EVENT " << fTNofEvents << " =================" << endl;

  //Int_t nofCreatedTracks = 0; //[R.K. 01/2017] unused variable?

  Int_t nofHits = 0;
  Int_t firstDH[4];
  
  // putting hits into one common array
  for ( Int_t idet = 0 ; idet < 4 ; idet++ ) {
    firstDH[idet] = nofHits;
    if ( fIncludeDet[idet] == kTRUE ) {
      nofHits += fHitArray[idet]->GetEntriesFast();
      //      cout << fDetName[idet].Data() << " hits from " << firstDH[idet] << " to " << nofHits << "   ( " << fHitArray[idet]->GetEntriesFast() << " hits )" << endl;
      for ( Int_t ihit = firstDH[idet] ; ihit < nofHits ; ihit++ ) { // hit mixing could go here
        fHitDetId.push_back(idet);
        fHitDetNo.push_back(ihit-firstDH[idet]);
      }
    }
  }
  std::vector<Bool_t> hitUsed(fHitDetId.size(),kFALSE);

  FairHit* detHit;
  // DRAW

  if ( fCanvas ) {
    fCanvas->Clear();
    Double_t minX = -60.;
    Double_t maxX = +60.;
    Double_t minY = -60.;
    Double_t maxY = +60.;
    
    fCanvas->Range(minX,minY,maxX,maxY);
    fCanvas->SetFillColor(0);
    fCanvas->SetBorderSize(0);
    fCanvas->SetBorderMode(0);
    fCanvas->SetFrameFillColor(0);
    fCanvas->SetGrid(0,0);
    fCanvas->SetLogy(0);
    
    // middle cross
    TLine* horLft = new TLine(-5., 0.,-2., 0.);  horLft->SetLineWidth(2);  horLft->Draw();
    TLine* horRgt = new TLine( 2., 0., 5., 0.);  horRgt->SetLineWidth(2);  horRgt->Draw();
    TLine* verTop = new TLine( 0., 2., 0., 5.);  verTop->SetLineWidth(2);  verTop->Draw();
    TLine* verBot = new TLine( 0.,-5., 0.,-2.);  verBot->SetLineWidth(2);  verBot->Draw();
    for ( Int_t xline = -100 ; xline <= 100 ; xline+=10 ) {
      TLine* verGrid = new TLine(xline,-100,xline,100);
      verGrid->SetLineStyle(2);
      verGrid->Draw();
    }
    for ( Int_t yline = -100 ; yline <= 100 ; yline+=10 ) {
      TLine* horGrid = new TLine(-100,yline,100,yline);
      horGrid->SetLineStyle(2);
      horGrid->Draw();
    }
    
    for ( Int_t iell = 0 ; iell < 1000 ; iell++ ) {
      fEllipse[iell]->SetX1(100);
      fEllipse[iell]->SetY1(100);
      fEllipse[iell]->SetR1(0);
      fEllipse[iell]->SetR2(0);
      fEllipse[iell]->Draw();
      fStageDesc->Draw();
    }
    
    for ( Int_t ihit = 0 ; ihit < nofHits ; ihit++ ) {
      
      if ( fHitDetId[ihit] == 0 || fHitDetId[ihit] == 1 ) {
        detHit = (PndSdsHit*)fHitArray[fHitDetId[ihit]]->At(fHitDetNo[ihit]);
        TMarker* marker = new TMarker(detHit->GetX(),detHit->GetY(),20);
        marker->SetMarkerColor(3);
        marker->Draw();
      }
      else                    if ( fHitDetId[ihit] == 2 ) {
        detHit = (PndSttHit*)fHitArray[fHitDetId[ihit]]->At(fHitDetNo[ihit]);
        Int_t iTube = ((PndSttHit*)detHit)->GetTubeID(); 
        PndSttTube *sttTube = nullptr;
        if (fGeoType == 1)
        sttTube = (PndSttTube*) fTubeArray->At(iTube);
    else if (fGeoType == 2){
        sttTube = PndSttTubeMap::Instance()->GetTube(iTube);
    }

        if ( sttTube->GetWireDirection().Z() < 1. )
          continue; // don't draw skewed hits (how else to put them?)

        TEllipse* ellipse = new TEllipse(detHit->GetX(),detHit->GetY(),((PndSttHit*)detHit)->GetIsochrone());
        ellipse->SetFillStyle(0);
        ellipse->SetLineWidth(2);
        ellipse->SetLineColor(2);
        ellipse->Draw(); 
      }
      else                    if ( fHitDetId[ihit] == 3 ) {
        detHit = (PndGemHit*)fHitArray[fHitDetId[ihit]]->At(fHitDetNo[ihit]);
        TMarker* marker = new TMarker(detHit->GetX(),detHit->GetY(),20);
        marker->SetMarkerColor(6);
        marker->Draw();
      }
    }
    fStageDesc->SetText(-50.,50.,"HITS");
    fStageDesc->Paint();
    fCanvas->Modified();
    fCanvas->Update();
    gSystem->Sleep(1000);
  }

  Int_t hitN = 0;
  // looping over hits as they are in the common array

  for ( Int_t ihit = 0 ; ihit < nofHits ; ihit++ ) {

    hitN = (hitN+gRandom->Integer(nofHits))%nofHits;
    while ( 1==1 ) {
      if ( hitUsed[hitN] == kFALSE ) break;
      hitN = (hitN+1)%nofHits;
    }
    hitUsed[hitN] = kTRUE;
    //    cout << " for hit " << ihit << " would like to take hit " << hitN << " >> " << fHitDetId[hitN] << "." << fHitDetNo[hitN] << endl;

    //    hitN = (firstDH[2]+ihit)%nofHits; // hit mixing could go here

    if ( fHitDetId[hitN] == 0 || fHitDetId[hitN] == 1 )
      detHit = (PndSdsHit*)fHitArray[fHitDetId[hitN]]->At(fHitDetNo[hitN]);
    else                    if ( fHitDetId[hitN] == 2 )
      detHit = (PndSttHit*)fHitArray[fHitDetId[hitN]]->At(fHitDetNo[hitN]);
    else                    if ( fHitDetId[hitN] == 3 )
      detHit = (PndGemHit*)fHitArray[fHitDetId[hitN]]->At(fHitDetNo[hitN]);
    
    if ( fVerbose > 3 ) {
      cout << "--------------------------------------------------------" << endl;
      cout << ihit << " (" << hitN << ") > " << fHitDetId[hitN] << "." << fHitDetNo[hitN] << " in " << fDetName[fHitDetId[hitN]].Data() << " at " 
           << detHit->GetX() << " " << detHit->GetY() << " " << detHit->GetZ() 
           << "    error " << detHit->GetDx() << " " << detHit->GetDy() << " " << detHit->GetDz() << endl;
    }
    if ( fVerbose > 3 )
      PrintTracks();

    Bool_t       sttHit = kFALSE;
    Bool_t skewedSttHit = kFALSE;
    if ( fHitDetId[hitN] == 2 ) {
      sttHit = kTRUE;
      Int_t iTube = ((PndSttHit*)detHit)->GetTubeID(); 
      PndSttTube *sttTube = nullptr;
      if (fGeoType == 1)
          sttTube = (PndSttTube*) fTubeArray->At(iTube);
      else if (fGeoType == 2){
          sttTube = PndSttTubeMap::Instance()->GetTube(iTube);
      }
      if (sttTube == nullptr)
          std::cout << "-E- PndBarrelTrackFinder::Exec no sttTube " << iTube << std::endl;
      if ( sttTube->GetWireDirection().Z() < 1. ) {
        skewedSttHit = kTRUE;
        if ( fVerbose > 3 ) 
          cout << "SKEWED HIT" << endl;
      }
    }

    Bool_t hitBelongsToKnownTrack = kFALSE;
    for ( size_t itr = 0 ; itr < fTracksVector.size() ; itr++ ) {
      if   ( !sttHit )      { if ( MatchHitWithZInfoTT  (detHit,fHitDetId[hitN],fHitDetNo[hitN],itr) ) hitBelongsToKnownTrack = kTRUE; }
      else 
        if ( skewedSttHit ) { if ( MatchSkewedSttHitTT  (detHit,fHitDetId[hitN],fHitDetNo[hitN],itr) ) hitBelongsToKnownTrack = kTRUE; }
        else                { if ( MatchParallelSttHitTT(detHit,fHitDetId[hitN],fHitDetNo[hitN],itr) ) hitBelongsToKnownTrack = kTRUE; }
    }
    if ( hitBelongsToKnownTrack ) continue;

    if ( skewedSttHit ) {
      AddHitToPreviousHits     (detHit,fHitDetId[hitN],fHitDetNo[hitN]);
      continue;
    }
    
    Bool_t hitMatchedToHit = kFALSE;
    if ( fVerbose > 3 )
      cout << "TRYING TO MATCH TO PREVIOUS " << fHitVector.size() << " HITS" << endl;
    for ( Int_t ihitP = fHitVector.size()-1 ; ihitP >= 0 ; ihitP-- ) {
      if ( MatchHitToHit(detHit,fHitDetId[hitN],fHitDetNo[hitN],ihitP) ) hitMatchedToHit = kTRUE;
    }
    if ( hitMatchedToHit ) continue;
    
    AddHitToPreviousHits     (detHit,fHitDetId[hitN],fHitDetNo[hitN]);
  }
  if ( fCanvas ) {
    DrawTracks();
    fStageDesc->SetText(-50.,50.,"TRACKS FOUND");
    fStageDesc->Paint();
    fCanvas->Modified();
    fCanvas->Update();
    gSystem->Sleep(1000);
  }
//   PrintTracks();
//  cout << "REMOVING TRACKS!" << endl;
  RemoveShortTracks();
  //  PrintTracks();
  if ( fVerbose > 0 )
    PrintTracks();
  if ( fCanvas ) {
    DrawTracks();
    fStageDesc->SetText(-50.,50.,"SHORT TRACKS REMOVED");
    fStageDesc->Paint();
    fCanvas->Modified();
    fCanvas->Update();
    gSystem->Sleep(1000);
  }

  // deal with left previous hits
  for ( Int_t iuh = fHitVector.size()-1 ; iuh >= 0 ; iuh-- ) {
    if ( fVerbose > 3 )
      cout << "PREVIOUS HIT " << iuh << " out of " << fHitVector.size() << endl;
    Bool_t       sttHit = kFALSE;
    Bool_t skewedSttHit = kFALSE;
    if ( fHitVectDI[iuh] == 2 ) {
      sttHit = kTRUE;
      Int_t iTube = ((PndSttHit*)fHitVector[iuh])->GetTubeID(); 
      PndSttTube *sttTube = nullptr;
      if (fGeoType == 1)
          sttTube = (PndSttTube*) fTubeArray->At(iTube);
      else if (fGeoType == 2){
          sttTube = PndSttTubeMap::Instance()->GetTube(iTube);
      }
      if ( sttTube->GetWireDirection().Z() < 1. ) {
        skewedSttHit = kTRUE;
      }
    }

    Bool_t hitBelongsToKnownTrack = kFALSE;
    for ( size_t itr = 0 ; itr < fTracksVector.size() ; itr++ ) {
      if ( HitBelongsToTrack(fHitVectDI[iuh],fHitVectHN[iuh],itr) ) continue;
      if   ( !sttHit )      { if ( MatchHitWithZInfoTT  (fHitVector[iuh],fHitVectDI[iuh],fHitVectHN[iuh],itr) ) hitBelongsToKnownTrack = kTRUE; }
      else 
        if ( skewedSttHit ) { if ( MatchSkewedSttHitTT  (fHitVector[iuh],fHitVectDI[iuh],fHitVectHN[iuh],itr) ) hitBelongsToKnownTrack = kTRUE; }
        else                { if ( MatchParallelSttHitTT(fHitVector[iuh],fHitVectDI[iuh],fHitVectHN[iuh],itr) ) hitBelongsToKnownTrack = kTRUE; }
    }
    if ( hitBelongsToKnownTrack ) {
      RemoveHitFromPreviousHits(iuh);
      continue;
    }

  }

  // cout << "==============================================================" << endl;
  // cout << "==============================================================" << endl;
  // cout << "there are still " << fHitVector.size() << " unused hits: " << flush;
  // for ( Int_t iuh = 0 ; iuh < fHitVector.size() ; iuh++ ) 
  //   cout << fHitVectDI[iuh] << "." << fHitVectHN[iuh] << " " << flush;
  // cout << endl;
  // cout << "==============================================================" << endl;
  //  cout << "==============================================================" << endl;
  //  cout << "FOUND " << fTracksVector.size() << " tracks:" << endl;
  if ( fVerbose > 0 )
    PrintTracks();
  if ( fCanvas ) {
    DrawTracks();
    fStageDesc->SetText(-50.,50.,"TRACKS REORGANIZED");
    fStageDesc->Paint();
    fCanvas->Modified();
    fCanvas->Update();
    gSystem->Sleep(1000);
  }
  //  cout << "==============================================================" << endl;
  //  cout << "REMOVE CLONE TRACKS" << endl;
  RemoveCloneTracks();
  if ( fVerbose > 0 )
    PrintTracks();
  if ( fCanvas ) {
    DrawTracks();
    fStageDesc->SetText(-50.,50.,"CLONE TRACKS REMOVED");
    fStageDesc->Paint();
    fCanvas->Modified();
    fCanvas->Update();
    gSystem->Sleep(5000);
  }

  WriteTracks();
  
  fTNofEvents++;
  fTNofTracks += fTracksVector.size();

  fTracksVector.clear();
  fHitDetId.clear();
  fHitDetNo.clear();
  fHitVector.clear();
  fHitVectDI.clear();
  fHitVectHN.clear();

  return;
}

Bool_t PndBarrelTrackFinder::ExtractMeanRPhiFromTrack ( Int_t  trackNo )

private

Definition at line 1053 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), CalcPhi(), CAMath::Cos(), Double_t, ExtractMeanZ_PFromTrack(), FindCircDist(), fReasonableDist, fTracksVector, fVerbose, n, CAMath::Sin(), x, y, and z_p.

Referenced by MatchHitWithZInfoTT(), and MatchParallelSttHitTT().

                                                                        {
  if ( fVerbose > 3 )
    cout << "extracting mean rphi from track " << trackNo << " with " << fTracksVector[trackNo].trackHits.size() << " hits." << endl;
  if ( fTracksVector[trackNo].trackHits.size() < 4 ) return kFALSE;
  if ( TMath::Abs(fTracksVector[trackNo].meanR) > 0.1 ) return kTRUE;
  
  Double_t meanPhi = 0.;
  Double_t meanRad = 0.;
  Double_t goodPNR = 0.;
  for ( size_t ipar = 0 ; ipar < fTracksVector[trackNo].trackPars.size() ; ipar++ ) {
    if ( fTracksVector[trackNo].trackPars[ipar].n < 0.9*fTracksVector[trackNo].trackHitD.size() ) continue;
    meanPhi += CalcPhi(fTracksVector[trackNo].trackPars[ipar].x,fTracksVector[trackNo].trackPars[ipar].y);
    meanRad += fTracksVector[trackNo].trackPars[ipar].r;
    goodPNR += 1.;
  }
  if ( goodPNR < 0.5 ) return kFALSE;
  meanRad = meanRad/goodPNR;
  meanPhi = meanPhi/goodPNR;
  if ( fVerbose > 3 )
    cout << "mean rad phi is " << meanRad << " " << meanPhi << endl;

  Double_t sH1[3] = {meanRad*TMath::Cos(meanPhi),
                     meanRad*TMath::Sin(meanPhi),
                     meanRad};
  Double_t sH2[3] = {0.,0.,0.};

  Int_t nofCloseHits = 0;
  for ( size_t ihit = 0 ; ihit < fTracksVector[trackNo].trackHits.size() ; ihit++ ) {
    sH2[0] = fTracksVector[trackNo].trackHits[ihit]->GetX();
    sH2[1] = fTracksVector[trackNo].trackHits[ihit]->GetY();
    sH2[2] = 0.;
    if ( fTracksVector[trackNo].trackHitD[ihit] == 2 ) {
      sH2[2] = ((PndSttHit*)fTracksVector[trackNo].trackHits[ihit])->GetIsochrone();
    }

    if ( FindCircDist(sH1,sH2) < fReasonableDist ) {
      nofCloseHits++;
    }
  }
  if ( fVerbose > 3 )
    cout << "there are " << nofCloseHits << " close hits." << endl;
  if ( nofCloseHits < 0.7*fTracksVector[trackNo].trackHitD.size() ) return kFALSE;

  fTracksVector[trackNo].meanX = sH1[0];
  fTracksVector[trackNo].meanY = sH1[1];
  fTracksVector[trackNo].meanR = sH1[2];

  for ( size_t ipar = 0 ; ipar < fTracksVector[trackNo].trackPars.size() ; ipar++ ) {
    if ( fTracksVector[trackNo].trackPars[ipar].n < 0.9*fTracksVector[trackNo].trackHitD.size() ) continue;
    if ( fTracksVector[trackNo].trackPars[ipar].z_p > 1.e6 ) continue;
    fTracksVector[trackNo].trackPars[ipar].n = 400;
  }

  ExtractMeanZ_PFromTrack(trackNo);

  return kTRUE;
}

Bool_t PndBarrelTrackFinder::ExtractMeanZ_PFromTrack ( Int_t  trackNo )

private

Definition at line 1113 of file PndBarrelTrackFinder.cxx.

References a, CAMath::Abs(), CalcPhi(), CalcZ_P(), CAMath::Cos(), Double_t, fGeoType, FindInterestingRegions(), fReasonableZ_PD, fTracksVector, fTubeArray, fVerbose, PndSttTube::GetHalfLength(), PndSttTube::GetPosition(), PndSttTubeMap::GetTube(), PndSttTube::GetWireDirection(), PndSttTubeMap::Instance(), and CAMath::Sqrt().

Referenced by ExtractMeanRPhiFromTrack(), MatchHitWithZInfoTT(), and MatchSkewedSttHitTT().

                                                                       {
  if ( fVerbose > 3 )
    cout << "extracting mean Z_P from track " << trackNo << " with " << fTracksVector[trackNo].trackHits.size() << " hits." << endl;
  if ( fTracksVector[trackNo].meanZ_P < 1.e6 ) return kTRUE;

  // count hits with z information
  Int_t nofZHits = 0;  
  Double_t circPar[4] = {fTracksVector[trackNo].meanX,
                         fTracksVector[trackNo].meanY,
                         fTracksVector[trackNo].meanR,
                         CalcPhi(fTracksVector[trackNo].meanX,fTracksVector[trackNo].meanY)};

  std::vector<Double_t> trackHitsZ_P;

  for ( size_t ihit = 0 ; ihit < fTracksVector[trackNo].trackHits.size() ; ihit++ ) {
    if ( fTracksVector[trackNo].trackHitD[ihit] != 2 ) {
      nofZHits++;
      trackHitsZ_P.push_back(CalcZ_P(circPar,
                                     fTracksVector[trackNo].trackHits[ihit]->GetX(),
                                     fTracksVector[trackNo].trackHits[ihit]->GetY(),
                                     fTracksVector[trackNo].trackHits[ihit]->GetZ()));
    }
    else {
      Int_t tubeId = ((PndSttHit*)fTracksVector[trackNo].trackHits[ihit])->GetTubeID();
      PndSttTube *sttTube = nullptr;
      if (fGeoType == 1)
          sttTube = (PndSttTube*) fTubeArray->At(tubeId);
      else if (fGeoType == 2){
          sttTube = PndSttTubeMap::Instance()->GetTube(tubeId);
      }
      
      if ( sttTube->GetWireDirection().Z() < 1. ) {
        nofZHits++;

        Double_t ra = TMath::ATan(TMath::Sqrt(sttTube->GetWireDirection().X()*sttTube->GetWireDirection().X()+
                                              sttTube->GetWireDirection().Y()*sttTube->GetWireDirection().Y())/
                                  sttTube->GetWireDirection().Z());
        Double_t a  = 0.5*((PndSttHit*)fTracksVector[trackNo].trackHits[ihit])->GetIsochrone()*(1.+1./TMath::Cos(ra));
        
        Double_t tubePar[8];
        tubePar[0] = sttTube->GetPosition().X();
        tubePar[1] = sttTube->GetPosition().Y();
        tubePar[2] = sttTube->GetPosition().Z();
        tubePar[3] = sttTube->GetHalfLength();
        tubePar[4] = sttTube->GetWireDirection().X();
        tubePar[5] = sttTube->GetWireDirection().Y();
        tubePar[6] = sttTube->GetWireDirection().Z();
        tubePar[7] = a;

        Double_t intReg[4];
        
        Int_t nofReg = FindInterestingRegions(circPar,tubePar,intReg);
        for ( Int_t ireg = 0 ; ireg < nofReg ; ireg++ ) {
          trackHitsZ_P.push_back(intReg[ireg]);
        }
      }
    }
  }
  if ( nofZHits < 4 ) return kFALSE;

  if ( fVerbose > 3 )
    cout << "trying to etract mean z_p. There are " << nofZHits << " hits with Z information, parameters z_p: " << flush;

  Int_t nofClosePars =  0;
  Int_t highestNPars =  2;
  Int_t highestNPPar = -1;
  Double_t meanZ_P     = 0.;
  Double_t sumOfDist   = 0.;
  Double_t smallestSOD = 1000.;
  Double_t bestMeanZ_P = 0.;
  if ( fVerbose > 3 )
    cout << " got " << trackHitsZ_P.size() << " z_p parameters: " << flush;
  for ( size_t ipar = 0 ; ipar < trackHitsZ_P.size() ; ipar++ ) {

    meanZ_P      = 0.;
    sumOfDist    = 0.;
    nofClosePars = 0;
    
    for ( size_t ipar2 = 0 ; ipar2 < trackHitsZ_P.size() ; ipar2++ ) {
      if ( TMath::Abs(trackHitsZ_P[ipar]-trackHitsZ_P[ipar2]) > fReasonableZ_PD ) continue;
      meanZ_P += trackHitsZ_P[ipar2];
      sumOfDist += TMath::Abs(trackHitsZ_P[ipar]-trackHitsZ_P[ipar2]);
      nofClosePars++;
    }
    
    if ( nofClosePars >= highestNPars ) {
      if ( nofClosePars > highestNPars ) smallestSOD = 1000.;
      if ( sumOfDist < smallestSOD ) {
        highestNPars = nofClosePars;
        highestNPPar = ipar;
        smallestSOD  = sumOfDist;
        bestMeanZ_P  = meanZ_P/((Double_t)(nofClosePars));
      }
    }
    
    if ( fVerbose > 3 )
      cout << trackHitsZ_P[ipar] << "(" << ipar << "/" << nofClosePars << "/" << sumOfDist << ") " << endl;
  }
  if ( fVerbose > 3 ) {
    cout << "best parameter number " << highestNPPar << " with " << highestNPars << " close parameters, sumOfDist = " << smallestSOD << " and meanZ_P = " << bestMeanZ_P << endl;
  }
  
  if ( highestNPars < 3 )
    return kFALSE;

  fTracksVector[trackNo].meanZ_P = bestMeanZ_P;

  return kTRUE;
}

Double_t PndBarrelTrackFinder::FindCircDist ( Double_t *  c1, Double_t *  c2  )

private

Definition at line 1625 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), Double_t, and CAMath::Sqrt().

Referenced by ExtractMeanRPhiFromTrack(), MatchHitWithZInfoTT(), and MatchParallelSttHitTT().

                                                                       {
  Double_t circDist = TMath::Sqrt((c1[0]-c2[0])*(c1[0]-c2[0])+
                                  (c1[1]-c2[1])*(c1[1]-c2[1]));
  if ( circDist >=            c1[2]+c2[2]  ) { return circDist-c1[2]-c2[2];}
  if ( circDist <= TMath::Abs(c1[2]-c2[2]) ) { return TMath::Abs(c1[2]-c2[2])-circDist;}
  if ( c1[2] > c2[2] ) { return c2[2]-TMath::Abs(circDist-c1[2]);}
  return c1[2]-TMath::Abs(circDist-c2[2]);
}

Bool_t PndBarrelTrackFinder::FindCircPar ( Double_t *  c1, Double_t *  c2, Double_t *  c3, Int_t  cno, Double_t *  cl  )

private

Definition at line 1580 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), b2, Double_t, and CAMath::Sqrt().

Referenced by MatchHitToHit(), MatchHitWithZInfoTT(), and MatchParallelSttHitTT().

                                                                                                          {
  if ( cno < 0 || cno > 7 ) return kFALSE;
//   cout << "c1 : " << c1[0] << "," << c1[1] << " - " << c1[2] << endl;
//   cout << "c2 : " << c2[0] << "," << c2[1] << " - " << c2[2] << endl;
//   cout << "c3 : " << c3[0] << "," << c3[1] << " - " << c3[2] << endl;
  Double_t s1 = -1+2*(cno%2); cno/=2;
  Double_t s2 = -1+2*(cno%2); cno/=2;
  Double_t s3 = -1+2*(cno%2);
  
  Double_t a2 = 2.*(c1[0]-c2[0]);
  Double_t b2 = 2.*(c1[1]-c2[1]);
  Double_t e2 = 2.*(s1*c1[2]-s2*c2[2]);
  Double_t d2 = (c1[0]*c1[0]+c1[1]*c1[1]-c1[2]*c1[2]) - (c2[0]*c2[0]+c2[1]*c2[1]-c2[2]*c2[2]);

  Double_t a3 = 2.*(c1[0]-c3[0]);
  Double_t b3 = 2.*(c1[1]-c3[1]);
  Double_t e3 = 2.*(s1*c1[2]-s3*c3[2]);
  Double_t d3 = (c1[0]*c1[0]+c1[1]*c1[1]-c1[2]*c1[2]) - (c3[0]*c3[0]+c3[1]*c3[1]-c3[2]*c3[2]);
  //  cout << "a2 " << a2 << " a3 " << a3 << " b2 " << b2 << " b3 " << b3 << " /// a2*b3-a3*b2 = " << a2*b3-a3*b2 << endl;
  if ( TMath::Abs(a2*b3-a3*b2) < 0.001 ) return kFALSE;

  Double_t ax = (d2*b3-d3*b2)/(a2*b3-a3*b2);
  Double_t bx = (e2*b3-e3*b2)/(a2*b3-a3*b2);
  Double_t ay = (a2*d3-a3*d2)/(a2*b3-a3*b2);
  Double_t by = (a2*e3-a3*e2)/(a2*b3-a3*b2);

  Double_t sa = bx*bx+by*by-1.;
  Double_t sb = 2.*bx*(ax-c1[0])+2.*by*(ay-c1[1])+2.*s1*c1[2];
  Double_t sc = (ax-c1[0])*(ax-c1[0])+(ay-c1[1])*(ay-c1[1])-c1[2]*c1[2];
  if ( sa == 0 ) return kFALSE;
  Double_t dt = sb*sb-4.*sa*sc;
  if ( dt < 0 ) return kFALSE;

  cl[2] = (-sb-TMath::Sqrt(dt))/(2.*sa);
  cl[0] = ax+bx*cl[2];
  cl[1] = ay+by*cl[2];
  cl[2] = TMath::Abs(cl[2]);
  if ( TMath::IsNaN(cl[2]) ) {
    return kFALSE;
  }
  return kTRUE;
}

Int_t PndBarrelTrackFinder::FindInterestingRegions ( Double_t *  circ, Double_t *  tube, Double_t *  reg  )

private

Definition at line 1636 of file PndBarrelTrackFinder.cxx.

References a, CAMath::Abs(), b, c, Double_t, Pi, CAMath::Sqrt(), x, y, z, and z_p.

Referenced by ExtractMeanZ_PFromTrack(), and MatchSkewedSttHitTT().

                                                                                                {
  Double_t ax = tube[4]/tube[6];
  Double_t ay = tube[5]/tube[6];
  Double_t a  = ax*ax+ay*ay;
  Double_t bx = (tube[0]-ax*tube[2]-circ[0]);
  Double_t by = (tube[1]-ay*tube[2]-circ[1]);
  Double_t b  = 2.*ax*bx+2.*ay*by;
  Double_t c  = bx*bx+by*by;

  Double_t delta[2];
  delta[0] = b*b-4.*a*(c-(circ[2]+tube[7])*(circ[2]+tube[7]));
  delta[1] = b*b-4.*a*(c-(circ[2]-tube[7])*(circ[2]-tube[7]));

  Int_t nofReg = 0;

  for ( Int_t idel = 0 ; idel < 2 ; idel++ ) {
    if ( delta[idel] < 0 ) continue;
    Double_t z[2];
    z[0] = (-b-TMath::Sqrt(delta[idel]))/2./a;
    z[1] = (-b+TMath::Sqrt(delta[idel]))/2./a;
    //    cout << " tube crosses circle in " << z[0] << " or " << z[1] << endl;
    
    for ( Int_t iz = 0 ; iz < 2 ; iz++ ) {
      Double_t z_p = 2.e6;
      if ( TMath::Abs(z[iz]-tube[2]) > tube[3] ) continue;
      //      cout << "checking z = " << z[iz] << endl;
      Double_t x  = tube[0] + ax * ( z[iz] - tube[2] );
      Double_t y  = tube[1] + ay * ( z[iz] - tube[2] );
      Double_t xd = x - circ[0];
      Double_t yd = y - circ[1];
      Double_t rd = TMath::Sqrt(xd*xd+yd*yd);
      Double_t pointPhi = TMath::ACos(xd/rd);
      if ( yd < 0 )
        pointPhi = 2.*TMath::Pi() - pointPhi;
      pointPhi += TMath::Pi();
      if ( pointPhi > 2.*TMath::Pi() ) pointPhi -= 2.*TMath::Pi();
      Double_t phiDiffr = (pointPhi-circ[3]);
      if ( phiDiffr >  TMath::Pi() ) phiDiffr -= 2.*TMath::Pi();
      if ( phiDiffr < -TMath::Pi() ) phiDiffr += 2.*TMath::Pi();
      z_p = z[iz]/TMath::Abs(phiDiffr);
      //      cout << "the z_p = " << z_p << " cause dz = " << z[iz] << " and dphi = " << phiDiffr << endl;
      //      cout << "got " << z_p << " while from other calc. it is: " << CalcZ_P(circ,x,y,z[iz]) << endl;
      reg[nofReg] = z_p;
      nofReg++;
    }
  }
  return nofReg;
}

void PndBarrelTrackFinder::Finish ( )

privatevirtual

Finish at the end of each event

Definition at line 1852 of file PndBarrelTrackFinder.cxx.

References Double_t, fBarrelTrackArray, fBarrelTrackCandArray, fClonesRDiff, fClonesRDiffRel, fClonesXDiff, fClonesXDiffRel, fClonesYDiff, fClonesYDiffRel, fClonesZ_PDiff, fClonesZ_PDiffRel, fDrawHistos, fTNofEvents, and fTNofTracks.

                                  {
  fBarrelTrackArray    ->Clear();
  fBarrelTrackCandArray->Clear();

  if ( fDrawHistos ) {
    new TCanvas("xdiff","xdiff",10,10,1000,800);
    fClonesXDiff     ->Draw();
    new TCanvas("ydiff","ydiff",10,10,1000,800);
    fClonesYDiff     ->Draw();
    new TCanvas("rdiff","rdiff",10,10,1000,800);
    fClonesRDiff     ->Draw();
    new TCanvas("zdiff","`diff",10,10,1000,800);
    fClonesZ_PDiff   ->Draw();
    new TCanvas("xdiffR","xdiffR",10,10,1000,800);
    fClonesXDiffRel  ->Draw();
    new TCanvas("ydiffR","ydiffR",10,10,1000,800);
    fClonesYDiffRel  ->Draw();
    new TCanvas("rdiffR","rdiffR",10,10,1000,800);
    fClonesRDiffRel  ->Draw();
    new TCanvas("zdiffR","zdiffR",10,10,1000,800);
    fClonesZ_PDiffRel->Draw();
  }

  cout << "-------------------- " << fName.Data() << " : Summary -----------------------" << endl;
  cout << " Events:        " << setw(10) << fTNofEvents << endl;
  cout << " Tracks:     " << setw(10) << fTNofTracks << "    ( " << (Double_t)fTNofTracks/((Double_t)fTNofEvents) << " per event )" << endl;
  cout << "--------------------------------------------------------------------------------" << endl; 
}

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(), PndIdealTrackFinder::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(), Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

{ return fPersistency; }

Bool_t PndBarrelTrackFinder::HitBelongsToTrack ( Int_t  detId, Int_t  hitNo, Int_t  trackNo  )

private

Definition at line 510 of file PndBarrelTrackFinder.cxx.

References fTracksVector.

Referenced by Exec().

                                                                                      {
  for ( size_t ihitP = 0 ; ihitP < fTracksVector[trackNo].trackHitD.size() ; ihitP++ ) {
    if ( fTracksVector[trackNo].trackHitD[ihitP] == detId &&
         fTracksVector[trackNo].trackHitN[ihitP] == hitNo ) 
      return kTRUE;
  }
  return kFALSE;
}

InitStatus PndBarrelTrackFinder::Init ( )

privatevirtual

Intialisation

Definition at line 1730 of file PndBarrelTrackFinder.cxx.

References fBarrelTrackArray, fBarrelTrackCandArray, fCanvas, fClonesRDiff, fClonesRDiffRel, fClonesXDiff, fClonesXDiffRel, fClonesYDiff, fClonesYDiffRel, fClonesZ_PDiff, fClonesZ_PDiffRel, fDetName, fDetType, fDrawDetails, fEllipse, fGeoType, fHitArray, fHitArrayName, PndSttMapCreator::FillTubeArray(), fIncludeDet, fStageDesc, fSttParameters, fTubeArray, PndGeoSttPar::GetGeometryType(), PndPersistencyTask::GetPersistency(), kGemHit, kMVDHitsPixel, kMVDHitsStrip, kSttHit, and rtdb.

                                      {
  //  fTrackSegments.clear();
  
  // Get and check FairRootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if( !ioman ) {
    cout << "-E- "<< GetName() <<"::Init: "
         << "RootManager not instantised!" << endl;
    return kERROR;
  }
  
  // Get the pointer to the singleton FairRunAna object
  FairRunAna* ana = FairRunAna::Instance();
  if(NULL == ana) {
    cout << "-E- "<< GetName() <<"::Init :"
         <<" no FairRunAna object!" << endl;
    return kERROR;
  }
  // Get the pointer to run-time data base
  FairRuntimeDb* rtdb = ana->GetRuntimeDb();
  if(NULL == rtdb) {
    cout << "-E- "<< GetName() <<"::Init :"
         <<" no runtime database!" << endl;
    return kERROR;
  }
  
  fDetName[0] = "MVD Pixel";
  fDetName[1] = "MVD Strip";
  fDetName[2] = "STT";
  fDetName[3] = "GEM";

  fHitArrayName[0] = "MVDHitsPixel";
  fHitArrayName[1] = "MVDHitsStrip";
  fHitArrayName[2] = "STTHit";
  fHitArrayName[3] = "GEMHit";

  fDetType[0] = kMVDHitsPixel;// FIXME: Use FairRootManager::Instance()->GetBranchId()
  fDetType[1] = kMVDHitsStrip;// FIXME: Use FairRootManager::Instance()->GetBranchId()
  fDetType[2] = kSttHit;
  fDetType[3] = kGemHit;


  for ( Int_t idet = 0 ; idet < 4 ; idet++ ) {
    if ( fIncludeDet[idet] == kTRUE ) {
      fHitArray[idet] = (TClonesArray*) ioman->GetObject(fHitArrayName[idet].Data());
      
      if ( !fHitArray[idet] ) {
        cout << "-E- " << GetName() << "::Init: No " << fHitArrayName[idet].Data() << " array!" << endl;
        return kERROR;
      }
      std::cout << "-I- " << GetName() << ": " << fHitArrayName[idet].Data() << " array found" << std::endl;
    }
  }
  
  // Create and register PndTrack and PndTrackCand arrays
  fBarrelTrackArray = new TClonesArray("PndTrack",100);
  ioman->Register("BarrelTrack", "Barrel Track", fBarrelTrackArray, GetPersistency());

  fBarrelTrackCandArray = new TClonesArray("PndTrackCand",100);
  ioman->Register("BarrelTrackCand", "Barrel TrackCand", fBarrelTrackCandArray, GetPersistency());

  std::cout << "-I- " << GetName() << ": Initialization successfull" << std::endl;
  std::cout << "-I- " << GetName() << ": Looking for tracks in " << flush;
  for ( Int_t idet = 0 ; idet < 4 ; idet++ ) {
    if ( fIncludeDet[idet] == kFALSE ) continue;
    cout << fDetName[idet].Data() << ", ";
  }
  cout << "\b\b. " << endl;
  cout << "================================================================================" << endl;

  fGeoType = fSttParameters->GetGeometryType();
  if ( fGeoType == 1){
      PndSttMapCreator *mapper = new PndSttMapCreator(fSttParameters);
      fTubeArray = mapper->FillTubeArray();
  }
//   cout << "****************************************************" << endl;
//   cout << "**  " << fTubeArray->GetEntriesFast() << "  TUBES  ***********************************" << endl;
//   cout << "****************************************************" << endl;
 
  // DRAW
  fCanvas = NULL;
  if ( fDrawDetails ) 
    fCanvas = new TCanvas("AnimatedCanvas","AnimatedCanvas",10,10,800,800);

  if ( fCanvas ) {
    for ( Int_t iell = 0 ; iell < 1000 ; iell++ ) {
      fEllipse[iell] = new TEllipse(100,100,0);
    }
    fStageDesc = new TLatex(0,0,"Init");
  }

  fClonesXDiff      = new TH1F("fClonesXDiff"     ,"fClonesXDiff"     ,20000,-100.,100.);
  fClonesYDiff      = new TH1F("fClonesYDiff"     ,"fClonesYDiff"     ,20000,-100.,100.);
  fClonesRDiff      = new TH1F("fClonesRDiff"     ,"fClonesRDiff"     ,20000,-100.,100.);
  fClonesZ_PDiff    = new TH1F("fClonesZ_PDiff"   ,"fClonesZ_PDiff"   ,20000,-100.,100.);
  fClonesXDiffRel   = new TH1F("fClonesXDiffRel"  ,"fClonesXDiffRel"  ,20000,-100.,100.);
  fClonesYDiffRel   = new TH1F("fClonesYDiffRel"  ,"fClonesYDiffRel"  ,20000,-100.,100.);
  fClonesRDiffRel   = new TH1F("fClonesRDiffRel"  ,"fClonesRDiffRel"  ,20000,-100.,100.);
  fClonesZ_PDiffRel = new TH1F("fClonesZ_PDiffRel","fClonesZ_PDiffRel",20000,-100.,100.);

  return kSUCCESS;
}

Bool_t PndBarrelTrackFinder::MatchHitToHit ( FairHit *  thisHit, Int_t  detId, Int_t  hitNo, Int_t  prevHNo  )

private

Definition at line 913 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), Bool_t, CalcPhi(), CalcZ_P(), Double_t, fGeoType, fHitVectDI, fHitVectHN, fHitVector, FindCircPar(), fReasonableZ_PD, fTracksVector, fTubeArray, fVerbose, PndSttTubeMap::GetTube(), PndSttTube::GetWireDirection(), PndSttTubeMap::Instance(), TrackBasis::meanR, TrackBasis::meanX, TrackBasis::meanY, TrackBasis::meanZ_P, TrackParameter::n, Pi, TrackParameter::r, RemoveHitFromPreviousHits(), CAMath::Sqrt(), TrackBasis::trackHitD, TrackBasis::trackHitN, TrackBasis::trackHits, TrackBasis::trackPars, TrackParameter::x, TrackParameter::y, and TrackParameter::z_p.

Referenced by Exec().

                                                                                                       {
  Double_t sH1[3] = {thisHit->GetX(),thisHit->GetY(),0.};
  if ( detId == 2 ) 
    sH1[2] = ((PndSttHit*)thisHit)->GetIsochrone();

  Double_t sH2[3] = {fHitVector[prevHNo]->GetX(),fHitVector[prevHNo]->GetY(),0.};
  if ( fHitVectDI[prevHNo] == 2 ) 
    sH2[2] = ((PndSttHit*)fHitVector[prevHNo])->GetIsochrone();

  Double_t sH3[3] = {0.,0.,0.}; // set to origin now, should be another hit (double for loop over hits)
  Double_t circPar[3];

  //      cout << fHitVector[prevHNo]->GetZ() << " " << flush;
  if ( fHitVectDI[prevHNo] == 2 ) {
    Int_t prevTube = ((PndSttHit*)fHitVector[prevHNo])->GetTubeID(); 
    PndSttTube *sttTube = nullptr;
    if (fGeoType == 1)
        sttTube = (PndSttTube*) fTubeArray->At(prevTube);
    else if (fGeoType == 2){
        sttTube = PndSttTubeMap::Instance()->GetTube(prevTube);
    }
    if ( sttTube->GetWireDirection().Z() < 1. )
      return kFALSE;
    sH2[2] = ((PndSttHit*)fHitVector[prevHNo])->GetIsochrone();
  }
  else {
    if ( detId == fHitVectDI[prevHNo] ) 
      if ( TMath::Abs(thisHit->GetZ()-fHitVector[prevHNo]->GetZ()) < 3. ) 
        return kFALSE;
    if ( detId != 2 ) {
      Double_t thisRad = TMath::Sqrt(sH1[0]*sH1[0]+sH1[1]*sH1[1]);
      Double_t prevRad = TMath::Sqrt(sH2[0]*sH2[0]+sH2[1]*sH2[1]);
      if ( fVerbose > 3 ) {
        cout << "matching r z " << thisRad << " " << thisHit->GetZ() << " (" << sH1[0] << " , " << sH1[1] << ") " << endl;
        cout << "    with r z " << prevRad << " " << fHitVector[prevHNo]->GetZ() << " (" << sH2[0] << " , " << sH2[1] << ") " << endl;
      }
    }
  }  


  Bool_t hitsMatched = kFALSE;

  Double_t thisHitPhi = CalcPhi(sH1[0],sH1[1]);
  Double_t prevHitPhi = CalcPhi(sH2[0],sH2[1]);
  
  Double_t phiDiffr = TMath::Abs(thisHitPhi-prevHitPhi);
  if ( phiDiffr > TMath::Pi()/10. && phiDiffr < TMath::Pi()*19./10. ) return kFALSE;
  
  Int_t stepCno = 1; 
  if ( sH1[2] < 1.e-5 ) stepCno  = 2;
  if ( sH2[2] < 1.e-5 ) stepCno += 2; 
//   cout << "stepCno = " << stepCno << " cause " << sH1[2] << " and " << sH2[2] << endl;
  
  Int_t fTN = -1;
  for ( Int_t icirc = 0 ; icirc < 4 ; icirc += stepCno ) {
    //    cout << icirc << endl;
    if ( FindCircPar(sH1,sH2,sH3,icirc,circPar) ) {
      //      cout << " ---> " << circPar[2] << endl;
      if ( circPar[2] < 10. || circPar[2] > 100000. ) continue;
      
      TrackParameter tp0;
      tp0.x = circPar[0];
      tp0.y = circPar[1];
      tp0.r = circPar[2];
      tp0.z_p = 2.e6;
      // temporary recognition of non-STT hit
      // should take only parameters created using hit with VALID z information
      if ( sH1[2] < 1.e-5 && sH2[2] < 1.e-5 ) {
        if ( TMath::Abs(CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ()) - CalcZ_P(circPar,sH2[0],sH2[1],fHitVector[prevHNo]->GetZ())) > 3.*fReasonableZ_PD ) {
          continue;
        }
        tp0.z_p = (CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ()) + CalcZ_P(circPar,sH2[0],sH2[1],fHitVector[prevHNo]->GetZ()))/2.;
      }
      else {
        if      ( sH1[2] < 1.e-5 ) tp0.z_p = CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ());
        else if ( sH2[2] < 1.e-5 ) tp0.z_p = CalcZ_P(circPar,sH2[0],sH2[1],fHitVector[prevHNo]->GetZ());
      }

      if ( fVerbose > 3 )
        if ( fHitVectDI[prevHNo] != 2 && detId != 2 ) {
          cout << " combining hit " << fHitVectDI[prevHNo] << "." << fHitVectHN[prevHNo] 
               << " with hit "      << detId << "." << hitNo << " --> z_p(1) = "
               << CalcZ_P(circPar,sH2[0],sH2[1],fHitVector[prevHNo]->GetZ()) << " --> z_p(2) = "
               << CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ()) << endl;
        }
      //          cout << "and new z_p = " << tp0.z_p << endl;

      tp0.n = 2;
      
      if ( fTN == -1 ) {
        TrackBasis tb0;
        tb0.trackPars.push_back(tp0);
        tb0.trackHits.push_back(fHitVector[prevHNo]);
        tb0.trackHitD.push_back(fHitVectDI[prevHNo]);
        tb0.trackHitN.push_back(fHitVectHN[prevHNo]);
        tb0.trackHits.push_back(thisHit);
        tb0.trackHitD.push_back(detId);
        tb0.trackHitN.push_back(hitNo);
        tb0.meanX   = 0.;
        tb0.meanY   = 0.;
        tb0.meanR   = 0.;
        tb0.meanZ_P = 2.e6;
        
        fTN = fTracksVector.size();
        fTracksVector.push_back(tb0);

        hitsMatched = kTRUE;
        //          cout << "creating track " << fTracksVector.size() << " from hits " 
        //               << fHitVectDI[prevHNo] << "." << fHitVectHN[prevHNo] << " and " 
        //               << fHitDetId[hitN] << "." << fHitDetNo[hitN] << endl;
        if ( fVerbose > 3 ) {
          cout << "!!! created track " << fTracksVector.size()-1 << " with " << fHitVectDI[prevHNo] << "." << fHitVectHN[prevHNo] << endl;
          //        if ( ihit == 54 ) {
          cout << "   >> with " << CalcPhi(tp0.x,tp0.y) << " (" << tp0.x << "," << tp0.y << ") " << tp0.r << " " << tp0.z_p << endl; //"(" << CalcPhi(tp0.x-sH1[0],tp0.y-sH1[1]) << " on " << detHit->GetZ() << ")" << endl;
          //          Double_t tmpx = sH1[0]-tp0.x;
          //          Double_t tmpy = sH1[1]-tp0.y;
          //          cout << " should be: " << sH1[0] << " " << sH1[1] << " ---> " << tmpx << " " << tmpy << endl;
          //          Double_t rotA = TMath::Pi()-CalcPhi(tp0.x,tp0.y);
          //          cout << "rotate by " << rotA*TMath::RadToDeg() << " deg" << endl;
          //          Double_t tmpx2 = tmpx*TMath::Cos(rotA)-tmpy*TMath::Sin(rotA);
          //          Double_t tmpy2 = tmpx*TMath::Sin(rotA)+tmpy*TMath::Cos(rotA);
          //          cout << " and then : " << tmpx   << " " << tmpy   << " ---> " << tmpx2 << " " << tmpy2 << endl;
          //          cout << " and finally: " << CalcPhi(tmpx2,tmpy2) << endl;
        }
      }
      else {
        fTracksVector[fTN].trackPars.push_back(tp0);
        if ( fVerbose > 3 ) {
          cout << "   >> with " << CalcPhi(tp0.x,tp0.y) << " (" << tp0.x << "," << tp0.y << ") " << tp0.r << " " << tp0.z_p << "(" << CalcPhi(tp0.x-sH1[0],tp0.y-sH1[1]) << " on " << thisHit->GetZ() << ")" << endl;
        }
      }
    }
  }

  if ( hitsMatched ) RemoveHitFromPreviousHits(prevHNo);
  return hitsMatched;
}

Bool_t PndBarrelTrackFinder::MatchHitToTrack ( FairHit *  thisHit, Int_t  detId, Int_t  hitNo, Int_t  trackNo  )

private

Definition at line 503 of file PndBarrelTrackFinder.cxx.

                                                                                                       {
  if ( thisHit ) cout << detId << "." << hitNo << " to " << trackNo << endl;
  return kFALSE;
}

Bool_t PndBarrelTrackFinder::MatchHitWithZInfoTT ( FairHit *  thisHit, Int_t  detId, Int_t  hitNo, Int_t  trackNo  )

private

Definition at line 748 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), Bool_t, CalcPhi(), CalcZ_P(), Double_t, ExtractMeanRPhiFromTrack(), ExtractMeanZ_PFromTrack(), FindCircDist(), FindCircPar(), fMaximalDist, fMaximalZ_PD, fReasonableDist, fReasonableZ_PD, fTracksVector, fVerbose, n, TrackParameter::n, Pi, TrackParameter::r, TrackParameter::x, x, TrackParameter::y, y, z_p, and TrackParameter::z_p.

Referenced by Exec().

                                                                                                               {
  Bool_t hitBelongsToTrack = kFALSE;

  Double_t sH1[3] = {thisHit->GetX(),thisHit->GetY(),0.};
  Double_t sH2[3]; // in the loop over track previous hits
  Double_t sH3[3] = {0.,0.,0.}; // set to origin now, should be another hit (double for loop over hits)
  Double_t sH4[3]; // another loop over hits
  Double_t circPar[3];
 
  Double_t thisHitPhi   = CalcPhi(sH1[0],sH1[1]);
  Double_t thisTrackPhi = CalcPhi(fTracksVector[trackNo].trackHits[0]->GetX(),fTracksVector[trackNo].trackHits[0]->GetY());
  Double_t phiHTDiffr   = TMath::Abs(thisHitPhi-thisTrackPhi);
  if ( phiHTDiffr > TMath::Pi()/4. && phiHTDiffr < TMath::Pi()*7./4. ) return kFALSE; // checking if hit and track match approximately, should make it radius dependent
  
  if ( TMath::Abs(fTracksVector[trackNo].meanR) > 1. ) {
    circPar[0] = fTracksVector[trackNo].meanX;
    circPar[1] = fTracksVector[trackNo].meanY;
    circPar[2] = fTracksVector[trackNo].meanR;

    if ( fVerbose > 3 )
      cout << "hit " << detId << "." << hitNo << " distance to track " << trackNo << " is " << FindCircDist(circPar,sH1) << " with errors ( "
           << ((PndGemHit*)thisHit)->GetDr() << " , " << ((PndGemHit*)thisHit)->GetDp() << " ) " << endl;

    if ( FindCircDist(circPar,sH1) < fReasonableDist )
      hitBelongsToTrack = kTRUE;
    if ( !hitBelongsToTrack ) return kFALSE;

    if ( fTracksVector[trackNo].meanZ_P > 1.e6 ) {
      TrackParameter tp0;
      tp0.x = circPar[0];
      tp0.y = circPar[1];
      tp0.r = circPar[2];
      tp0.z_p = CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ());
      tp0.n   = 400;
      fTracksVector[trackNo].trackPars.push_back(tp0);

      // putting thisHit into track
      fTracksVector[trackNo].trackHits.push_back(thisHit);
      fTracksVector[trackNo].trackHitD.push_back(detId);
      fTracksVector[trackNo].trackHitN.push_back(hitNo);
      ExtractMeanZ_PFromTrack(trackNo);
    }
    else {
      //      cout << "should compare if new found z_p " << CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ()) << " is close enough to mean: " << fTracksVector[trackNo].meanZ_P << endl; 

      if ( TMath::Abs(CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ())-fTracksVector[trackNo].meanZ_P) < fReasonableZ_PD ) {
        // putting thisHit into track
        fTracksVector[trackNo].trackHits.push_back(thisHit);
        fTracksVector[trackNo].trackHitD.push_back(detId);
        fTracksVector[trackNo].trackHitN.push_back(hitNo);
      }
      else
        return kFALSE;
    }
    if ( fVerbose > 3 )
      cout << " This hit belongs to track " << trackNo << " with " << CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ()) << endl;
  }
  else {
    // track has no mean information yet
    for ( size_t ipar = 0 ; ipar < fTracksVector[trackNo].trackPars.size() ; ipar++ ) {
      if ( fTracksVector[trackNo].trackPars[ipar].n < 0.6*fTracksVector[trackNo].trackHits.size() ) continue;
      
      circPar[0] = fTracksVector[trackNo].trackPars[ipar].x;
      circPar[1] = fTracksVector[trackNo].trackPars[ipar].y;
      circPar[2] = fTracksVector[trackNo].trackPars[ipar].r;
      
      if ( fVerbose > 4 ) {
        cout << " ..... to parameter (" << trackNo << "." << ipar << ") " << CalcPhi(circPar[0],circPar[1]) << " " << circPar[2] << ": distance is " << FindCircDist(sH1,circPar) << flush;
        cout << "   new z_p is " << CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ()) << " vs " << fTracksVector[trackNo].trackPars[ipar].z_p << " old" << flush;
        cout << endl;
      }
      
      if ( ( FindCircDist(sH1,circPar) < fMaximalDist &&
             TMath::Abs(CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ())-fTracksVector[trackNo].trackPars[ipar].z_p) < 3.*fReasonableZ_PD ) ||
           ( FindCircDist(sH1,circPar) < 3.*fReasonableDist &&
             TMath::Abs(CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ())-fTracksVector[trackNo].trackPars[ipar].z_p) < fMaximalZ_PD ) ) { 
        hitBelongsToTrack = kTRUE;
        
        fTracksVector[trackNo].trackPars[ipar].n += 1; // because thisHit is close to parameter ipar
        
        if ( fVerbose > 4 ) 
          cout << " par: " << CalcPhi(circPar[0],circPar[1]) << " | " << circPar[2] << " of track " << trackNo 
               << " matches to this hit with dist " << FindCircDist(sH1,circPar) << endl;
      }     
    }
    if ( !hitBelongsToTrack ) return kFALSE;
    
    // calculating track parameters with previous hits belonging to the track
    for ( size_t ihitP = 0 ; ihitP < fTracksVector[trackNo].trackHitD.size() ; ihitP++ ) {
      sH2[0] = fTracksVector[trackNo].trackHits[ihitP]->GetX();
      sH2[1] = fTracksVector[trackNo].trackHits[ihitP]->GetY();
      sH2[2] = 0.;
      if ( fTracksVector[trackNo].trackHitD[ihitP] == 2 )
        sH2[2] = ((PndSttHit*)fTracksVector[trackNo].trackHits[ihitP])->GetIsochrone();
      
      Int_t stepCno = 2;
      if ( sH2[2] < 1.e-5 ) stepCno += 2; 
      
      for ( Int_t icirc = 0 ; icirc < 4 ; icirc += stepCno ) {
        if ( FindCircPar(sH1,sH2,sH3,icirc,circPar) ) {
          if ( circPar[2] < 10. || circPar[2] > 100000. ) continue; // for now remove too egzotic tracks
          
          Double_t newParZ_P = CalcZ_P(circPar,sH1[0],sH1[1],thisHit->GetZ());
          
          Int_t nofCloseHits = 2;
          for ( size_t ihitF = 0 ; ihitF < fTracksVector[trackNo].trackHitD.size() ; ihitF++ ) {
            if ( ihitF == ihitP ) continue;
            sH4[0] = fTracksVector[trackNo].trackHits[ihitF]->GetX();
            sH4[1] = fTracksVector[trackNo].trackHits[ihitF]->GetY();
            sH4[2] = 0.;
            if ( fTracksVector[trackNo].trackHitD[ihitF] == 2 )
              sH4[2] = ((PndSttHit*)fTracksVector[trackNo].trackHits[ihitF])->GetIsochrone();
            
//          if ( FindCircDist(sH4,circPar) < fMaximalDist ||     // if the hit is very very close to the new parameter in xy plane
//               ( FindCircDist(sH4,circPar) < fReasonableDist &&   // if the hit is close to the new parameter in xy plane and pz space
//                 TMath::Abs(CalcZ_P(circPar,sH4[0],sH4[1],fTracksVector[trackNo].trackHits[ihitF]->GetZ())-newParZ_P) < fMaximalZ_PD ) ) { 
            if ( ( FindCircDist(sH4,circPar) < fMaximalDist &&
                   TMath::Abs(CalcZ_P(circPar,sH4[0],sH4[1],fTracksVector[trackNo].trackHits[ihitF]->GetZ())-newParZ_P) < 3.*fReasonableZ_PD ) ||
                 ( FindCircDist(sH4,circPar) < 3.*fReasonableDist &&
                   TMath::Abs(CalcZ_P(circPar,sH4[0],sH4[1],fTracksVector[trackNo].trackHits[ihitF]->GetZ())-newParZ_P) < fMaximalZ_PD ) ) { 
              nofCloseHits += 1;
            }
          }
          
          if ( fVerbose > 5 ) 
            cout << "there are " << nofCloseHits << " with limit set to " << 0.6*fTracksVector[trackNo].trackHits.size() << " (0.6*" << fTracksVector[trackNo].trackHits.size() << ")" << endl;
          if ( nofCloseHits < 0.6*fTracksVector[trackNo].trackHits.size() ) continue; // the track param has not enough close hits among previous track hits
          
          TrackParameter tp0;
          tp0.x   = circPar[0];
          tp0.y   = circPar[1];
          tp0.r   = circPar[2];
          tp0.z_p = newParZ_P;
          tp0.n   = nofCloseHits;
          
          fTracksVector[trackNo].trackPars.push_back(tp0);
          
        }
      }
    }

    // putting thisHit into track
    fTracksVector[trackNo].trackHits.push_back(thisHit);
    fTracksVector[trackNo].trackHitD.push_back(detId);
    fTracksVector[trackNo].trackHitN.push_back(hitNo);

    if ( fVerbose > 5 ) {
      cout << "HIT DOES " << (hitBelongsToTrack?"":"NOT ") << "BELONG TO TRACK " << trackNo << " WITH "
           << fTracksVector[trackNo].trackHits.size() << " HITS, " << fTracksVector[trackNo].trackPars.size() << " PARS: phi,rad,z_p,n: " << endl;
      for ( size_t ipar = 0 ; ipar < fTracksVector[trackNo].trackPars.size() ; ipar++ ) 
        cout << CalcPhi(fTracksVector[trackNo].trackPars[ipar].x,fTracksVector[trackNo].trackPars[ipar].y) << " " 
             << fTracksVector[trackNo].trackPars[ipar].r << " "
             << fTracksVector[trackNo].trackPars[ipar].z_p << " "
             << fTracksVector[trackNo].trackPars[ipar].n << "  |  " << flush;
      cout << endl;
    }

    ExtractMeanRPhiFromTrack(trackNo);
  }  

  return kTRUE;
}

Bool_t PndBarrelTrackFinder::MatchParallelSttHitTT ( FairHit *  thisHit, Int_t  detId, Int_t  hitNo, Int_t  trackNo  )

private

Definition at line 615 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), Bool_t, CalcPhi(), CalcZ_P(), Double_t, ExtractMeanRPhiFromTrack(), FindCircDist(), FindCircPar(), fMaximalDist, fMaximalZ_PD, fReasonableDist, fTracksVector, fVerbose, n, TrackParameter::n, Pi, TrackParameter::r, TrackParameter::x, x, TrackParameter::y, y, and TrackParameter::z_p.

Referenced by Exec().

                                                                                                                {
  Bool_t hitBelongsToTrack = kFALSE;

  Double_t sH1[3] = {thisHit->GetX(),thisHit->GetY(),0.};
  sH1[2] = ((PndSttHit*)thisHit)->GetIsochrone();
  Double_t sH2[3]; // in the loop over track previous hits
  Double_t sH3[3] = {0.,0.,0.}; // set to origin now, should be another hit (double for loop over hits)
  Double_t sH4[3]; // another loop over hits
  Double_t circPar[3];
 
  Double_t thisHitPhi   = CalcPhi(sH1[0],sH1[1]);
  Double_t thisTrackPhi = CalcPhi(fTracksVector[trackNo].trackHits[0]->GetX(),fTracksVector[trackNo].trackHits[0]->GetY());
  Double_t phiHTDiffr   = TMath::Abs(thisHitPhi-thisTrackPhi);
  if ( phiHTDiffr > TMath::Pi()/4. && phiHTDiffr < TMath::Pi()*7./4. ) return kFALSE; // checking if hit and track match approximately, should make it radius dependent
  
  if ( TMath::Abs(fTracksVector[trackNo].meanR) > 1. ) {
    circPar[0] = fTracksVector[trackNo].meanX;
    circPar[1] = fTracksVector[trackNo].meanY;
    circPar[2] = fTracksVector[trackNo].meanR;

    if ( FindCircDist(circPar,sH1) < fReasonableDist )
      hitBelongsToTrack = kTRUE;
    if ( !hitBelongsToTrack ) return kFALSE;
    
    // putting thisHit into track
    fTracksVector[trackNo].trackHits.push_back(thisHit);
    fTracksVector[trackNo].trackHitD.push_back(detId);
    fTracksVector[trackNo].trackHitN.push_back(hitNo);

    if ( fVerbose > 3 )
      cout << "this stt hit belongs to track " << trackNo << " (cause dist = " << FindCircDist(circPar,sH1) << ")" << endl;
  }
  else {
    // track has no mean information yet
    for ( size_t ipar = 0 ; ipar < fTracksVector[trackNo].trackPars.size() ; ipar++ ) {
      if ( fTracksVector[trackNo].trackPars[ipar].n < 0.9*fTracksVector[trackNo].trackHits.size() ) continue;
      
      circPar[0] = fTracksVector[trackNo].trackPars[ipar].x;
      circPar[1] = fTracksVector[trackNo].trackPars[ipar].y;
      circPar[2] = fTracksVector[trackNo].trackPars[ipar].r;
      
      if ( fVerbose > 4 ) {
        cout << " ..... to parameter (" << trackNo << "." << ipar << ") " << CalcPhi(circPar[0],circPar[1]) << " " << circPar[2] << ": distance is " << FindCircDist(sH1,circPar) << flush;
        cout << endl;
      }
    
      if ( FindCircDist(sH1,circPar) < fMaximalDist ) {
        hitBelongsToTrack = kTRUE;
      
        fTracksVector[trackNo].trackPars[ipar].n += 1; // because thisHit is close to parameter ipar
      
        if ( fVerbose > 4 ) 
          cout << " par: " << CalcPhi(circPar[0],circPar[1]) << " | " << circPar[2] << " of track " << trackNo 
               << " matches to this hit with dist " << FindCircDist(sH1,circPar) << endl;
      }     
    }
    if ( !hitBelongsToTrack ) return kFALSE;

    // calculating track parameters with previous hits belonging to the track
    for ( size_t ihitP = 0 ; ihitP < fTracksVector[trackNo].trackHitD.size() ; ihitP++ ) {
      sH2[0] = fTracksVector[trackNo].trackHits[ihitP]->GetX();
      sH2[1] = fTracksVector[trackNo].trackHits[ihitP]->GetY();
      sH2[2] = 0.;
      if ( fTracksVector[trackNo].trackHitD[ihitP] == 2 )
        sH2[2] = ((PndSttHit*)fTracksVector[trackNo].trackHits[ihitP])->GetIsochrone();
      
      Int_t stepCno = 1; 
      if ( sH2[2] < 1.e-5 ) stepCno += 2; 
      
      for ( Int_t icirc = 0 ; icirc < 4 ; icirc += stepCno ) {
        if ( FindCircPar(sH1,sH2,sH3,icirc,circPar) ) {
          if ( circPar[2] < 10. || circPar[2] > 100000. ) continue; // for now remove too egzotic tracks
          
          Double_t newParZ_P = 2.e6;
          if ( sH2[2] < 1.e-5 ) newParZ_P = CalcZ_P(circPar,sH2[0],sH2[1],fTracksVector[trackNo].trackHits[ihitP]->GetZ());
          
          Int_t nofCloseHits = 2;
          for ( size_t ihitF = 0 ; ihitF < fTracksVector[trackNo].trackHitD.size() ; ihitF++ ) {
            if ( ihitF == ihitP ) continue;
            sH4[0] = fTracksVector[trackNo].trackHits[ihitF]->GetX();
            sH4[1] = fTracksVector[trackNo].trackHits[ihitF]->GetY();
            sH4[2] = 0.;
            if ( fTracksVector[trackNo].trackHitD[ihitF] == 2 )
              sH4[2] = ((PndSttHit*)fTracksVector[trackNo].trackHits[ihitF])->GetIsochrone();
            
            if ( FindCircDist(sH4,circPar) < fMaximalDist ||     // if the hit is very very close to the new parameter in xy plane
                 ( newParZ_P < 1.e6 && FindCircDist(sH4,circPar) < fReasonableDist &&   // if the hit is close to the new parameter in xy plane and pz space
                   TMath::Abs(CalcZ_P(circPar,sH4[0],sH4[1],fTracksVector[trackNo].trackHits[ihitF]->GetZ())-newParZ_P) < fMaximalZ_PD ) ) { 
              nofCloseHits += 1;
            }
          }
          
          if ( fVerbose > 5 ) 
            cout << "there are " << nofCloseHits << " with limit set to " << 0.6*fTracksVector[trackNo].trackHits.size() << " (0.6*" << fTracksVector[trackNo].trackHits.size() << ")" << endl;
          if ( nofCloseHits < 0.6*fTracksVector[trackNo].trackHits.size() ) continue; // the track param has not enough close hits among previous track hits
          
          TrackParameter tp0;
          tp0.x   = circPar[0];
          tp0.y   = circPar[1];
          tp0.r   = circPar[2];
          tp0.z_p = newParZ_P;
          tp0.n   = nofCloseHits;
          
          fTracksVector[trackNo].trackPars.push_back(tp0);
          
        }
      }
    }

    // putting thisHit into track
    fTracksVector[trackNo].trackHits.push_back(thisHit);
    fTracksVector[trackNo].trackHitD.push_back(detId);
    fTracksVector[trackNo].trackHitN.push_back(hitNo);

    if ( fVerbose > 5 ) {
      cout << "HIT DOES " << (hitBelongsToTrack?"":"NOT ") << "BELONG TO TRACK " << trackNo << " WITH "
           << fTracksVector[trackNo].trackHits.size() << " HITS, " << fTracksVector[trackNo].trackPars.size() << " PARS: phi,rad,z_p,n: " << endl;
      for ( size_t ipar = 0 ; ipar < fTracksVector[trackNo].trackPars.size() ; ipar++ ) 
        cout << CalcPhi(fTracksVector[trackNo].trackPars[ipar].x,fTracksVector[trackNo].trackPars[ipar].y) << " " 
             << fTracksVector[trackNo].trackPars[ipar].r << " "
             << fTracksVector[trackNo].trackPars[ipar].z_p << " "
             << fTracksVector[trackNo].trackPars[ipar].n << "  |  " << flush;
      cout << endl;
    }

    ExtractMeanRPhiFromTrack(trackNo);
  }

  return kTRUE;
}

Bool_t PndBarrelTrackFinder::MatchSkewedSttHitTT ( FairHit *  thisHit, Int_t  detId, Int_t  hitNo, Int_t  trackNo  )

private

Definition at line 521 of file PndBarrelTrackFinder.cxx.

References a, CAMath::Abs(), Bool_t, CalcPhi(), CAMath::Cos(), Double_t, ExtractMeanZ_PFromTrack(), fGeoType, FindInterestingRegions(), fReasonableZ_PD, fTracksVector, fTubeArray, fVerbose, PndSttTube::GetHalfLength(), PndSttTube::GetPosition(), PndSttTubeMap::GetTube(), PndSttTube::GetWireDirection(), PndSttTubeMap::Instance(), TrackParameter::n, TrackParameter::r, CAMath::Sqrt(), TrackParameter::x, TrackParameter::y, and TrackParameter::z_p.

Referenced by Exec().

                                                                                                               {
  Bool_t printInfo = kFALSE;//kTRUE;
//   if ( fVerbose > 4 ) 
//     cout << "skewed hit" << endl;
  if ( TMath::Abs(fTracksVector[trackNo].meanR) < 1. ) return kFALSE;

  Int_t iTube = ((PndSttHit*)thisHit)->GetTubeID(); 
  PndSttTube *sttTube = nullptr;
  if (fGeoType == 1)
      sttTube = (PndSttTube*) fTubeArray->At(iTube);
  else if (fGeoType == 2){
      sttTube = PndSttTubeMap::Instance()->GetTube(iTube);
  }

  Double_t ra = TMath::ATan(TMath::Sqrt(sttTube->GetWireDirection().X()*sttTube->GetWireDirection().X()+
                                        sttTube->GetWireDirection().Y()*sttTube->GetWireDirection().Y())/
                            sttTube->GetWireDirection().Z());
  Double_t a  = 0.5*((PndSttHit*)thisHit)->GetIsochrone()*(1.+1./TMath::Cos(ra));
      
  Double_t tubePar[8];
  tubePar[0] = sttTube->GetPosition().X();
  tubePar[1] = sttTube->GetPosition().Y();
  tubePar[2] = sttTube->GetPosition().Z();
  tubePar[3] = sttTube->GetHalfLength();
  tubePar[4] = sttTube->GetWireDirection().X();
  tubePar[5] = sttTube->GetWireDirection().Y();
  tubePar[6] = sttTube->GetWireDirection().Z();
  tubePar[7] = a;
  
  if ( fVerbose > 4 || printInfo )
    cout << "matching to track " << trackNo << " with " << fTracksVector[trackNo].trackHits.size() << " hits" << endl;
  
  Bool_t trackFits = kFALSE;
  
  Double_t meanPhi = CalcPhi(fTracksVector[trackNo].meanX,fTracksVector[trackNo].meanY);
  Double_t meanCirc[4] = {fTracksVector[trackNo].meanX,
                          fTracksVector[trackNo].meanY,
                          fTracksVector[trackNo].meanR,
                          meanPhi};
  Double_t intReg[4];
  
  Int_t nofReg = FindInterestingRegions(meanCirc,tubePar,intReg);
  if ( nofReg==0 ) return kFALSE;
  
  if ( fVerbose > 4 || printInfo ) {
    cout << " --> connects to "   << trackNo << " with (" 
         << meanPhi << " " 
         << meanCirc[2] << " "
         << fTracksVector[trackNo].meanZ_P << ") at z_p = " << flush;
      for ( Int_t ireg = 0 ; ireg < nofReg ; ireg++ ) 
        cout << intReg[ireg] << "   " << flush;
      cout << endl;
  }

  // there is no mean z_p, attach all we can
  if ( fTracksVector[trackNo].meanZ_P > 1.e6 ) {
    for ( Int_t ireg = 0 ; ireg < nofReg ; ireg++ ) {
      TrackParameter tp0;
      tp0.x = meanCirc[0];
      tp0.y = meanCirc[1];
      tp0.r = meanCirc[2];
      tp0.z_p = intReg[ireg];
      tp0.n   = 400/nofReg;//fTracksVector[trackNo].trackHits.size();
      fTracksVector[trackNo].trackPars.push_back(tp0);
      if ( fVerbose > 4 || printInfo )
        cout << "creating another trackPar to track " << trackNo << " with z_p -> " << intReg[ireg] << endl;
    }
    fTracksVector[trackNo].trackHits.push_back(thisHit);
    fTracksVector[trackNo].trackHitD.push_back(detId);
    fTracksVector[trackNo].trackHitN.push_back(hitNo);

    ExtractMeanZ_PFromTrack(trackNo);
  }
  // there is already mean z_p !!!, check if any is good
  else {
    //    cout << "should compare if any of new found z_p: " << flush;
    for ( Int_t ireg = 0 ; ireg < nofReg ; ireg++ ) {
      //cout << intReg[ireg] << " " << flush;
      if ( TMath::Abs(intReg[ireg]-fTracksVector[trackNo].meanZ_P) < fReasonableZ_PD )
        trackFits = kTRUE;
    }
    //    cout << " is close enough to mean: " << fTracksVector[trackNo].meanZ_P << endl; 
    if ( !trackFits ) return kFALSE;

    fTracksVector[trackNo].trackHits.push_back(thisHit);
    fTracksVector[trackNo].trackHitD.push_back(detId);
    fTracksVector[trackNo].trackHitN.push_back(hitNo);
  }

  return kTRUE;
}

void PndBarrelTrackFinder::PrintTracks ( )

private

Definition at line 1225 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), Bool_t, CalcPhi(), Double_t, fDetName, fGeoType, fIncludeDet, fTracksVector, fTubeArray, PndSttTubeMap::GetTube(), PndSttTube::GetWireDirection(), PndSttTubeMap::Instance(), Pi, CAMath::Sqrt(), x, and y.

Referenced by Exec().

                                         {
  if ( fTracksVector.size() > 0 ) {
    cout << "          --------- " << flush;
    for ( Int_t idet = 0 ; idet < 4 ; idet++ ) {
      if ( fIncludeDet[idet] == kTRUE ) {
        cout << "\033[" << 91+idet << "m" 
             << fDetName[idet].Data() 
             << "\033[0m " << flush;
        if ( idet == 2 ) 
          cout << "\033[33m(skewed)\033[0m " << flush;
        
        //      cout << "\033[" << (skewedSttHit?33:91+fTracksVector[itr].trackHitD[ihit]) << "m" 
        //           << fTracksVector[itr].trackHitN[ihit] << " " 
        //           << "\033[0m" << flush;
      }
    }
    cout << "-------------" << endl;
  }

  for ( size_t itr = 0 ; itr < fTracksVector.size() ; itr++ ) {
    cout << "          |   track " << setw(2) << itr << ": " << flush;
    
    if ( TMath::Abs(fTracksVector[itr].meanR) > 0.1 ) {
      for ( size_t ihit = 0 ; ihit < fTracksVector[itr].trackHitD.size() ; ihit++ ) {
        Bool_t skewedSttHit = kFALSE;
        if ( fTracksVector[itr].trackHitD[ihit] == 2 ) {
          Int_t iTube = ((PndSttHit*) fTracksVector[itr].trackHits[ihit])->GetTubeID(); 
            PndSttTube *sttTube = nullptr;
            if (fGeoType == 1)
                sttTube = (PndSttTube*) fTubeArray->At(iTube);
            else if (fGeoType == 2){
                sttTube = PndSttTubeMap::Instance()->GetTube(iTube);
            }
          if ( sttTube->GetWireDirection().Z() < 1. ) {
            skewedSttHit = kTRUE;
          }
        }

        cout << "\033[" << (skewedSttHit?33:91+fTracksVector[itr].trackHitD[ihit]) << "m" 
             << fTracksVector[itr].trackHitN[ihit] << " " 
             << "\033[0m" << flush;
      }
      cout << endl;
      
      Double_t thisPhi = CalcPhi(fTracksVector[itr].meanX,
                                 fTracksVector[itr].meanY);
      Double_t thisRad =         fTracksVector[itr].meanR;
      Double_t thisZ_P =         fTracksVector[itr].meanZ_P;
      if ( TMath::IsNaN(thisZ_P) ) thisZ_P = 2.e6;
      
      Double_t calcPt = thisRad/165.974;
      Double_t calcPz = thisZ_P*0.0058997;
      Double_t calcP  = TMath::Sqrt(calcPt*calcPt+calcPz*calcPz);
      Double_t calcThe = TMath::ACos(calcPz/calcP);
      
      Double_t calcPhi = thisPhi-TMath::Pi()/2.;
      //Double_t calcChg = -1.; //[R.K. 01/2017] unused variable?
      if ( calcPhi < 0. ) calcPhi += TMath::Pi()*2.;
      Double_t firpPhi = CalcPhi(fTracksVector[itr].trackHits[0]->GetX(),fTracksVector[itr].trackHits[0]->GetY());
      if ( TMath::Abs(calcPhi-firpPhi) > TMath::Pi()/2. && TMath::Abs(calcPhi-firpPhi) < 3.*TMath::Pi()/2. ) {
        calcPhi = thisPhi+TMath::Pi()/2.;
        if ( calcPhi > 2.*TMath::Pi() ) calcPhi -= TMath::Pi()*2.;
        //calcChg = 1.; //[R.K. 01/2017] unused variable?
      }
      if ( calcPhi > TMath::Pi() )
        calcPhi -= 2.*TMath::Pi();
      
      TVector3 trackMomentum;
      trackMomentum.SetMagThetaPhi(calcP,calcThe,calcPhi);

      cout << "          | " 
           << "\033[96m" 
           << setw(7) << setprecision(6) << thisPhi << " " 
           << setw(7) << setprecision(6) << thisRad << " " 
           << setw(7) << setprecision(6) << thisZ_P << " " 
           << setw(7) << setprecision(6) << calcPt  << " " 
           << setw(7) << setprecision(6) << calcPz  << " " 
           << setw(7) << setprecision(6) << calcP   << " " 
           << setw(7) << setprecision(6) << calcPhi << " " 
           << setw(7) << setprecision(6) << calcThe << "\033[0m" << endl;
    }

    else {
      for ( size_t ihit = 0 ; ihit < fTracksVector[itr].trackHitD.size() ; ihit++ )
        cout << setw(1) << fTracksVector[itr].trackHitD[ihit] << "." << fTracksVector[itr].trackHitN[ihit] << " " << flush;
      cout << endl;

      for ( size_t ipar = 0 ; ipar < fTracksVector[itr].trackPars.size() ; ipar++ ) {
        cout << "          |           /" << setw(2) << ipar << " " 
             << setw(6) << setprecision(6) << CalcPhi(fTracksVector[itr].trackPars[ipar].x,fTracksVector[itr].trackPars[ipar].y) << " " 
             << setw(6) << setprecision(6) << fTracksVector[itr].trackPars[ipar].r << " " 
             << setw(6) << setprecision(6) << fTracksVector[itr].trackPars[ipar].z_p << " "
             << setw(6) << setprecision(6) << fTracksVector[itr].trackPars[ipar].n << endl;
      }
    }
    
    cout << "          ------------------------------------------------------------" << endl;
  }
}

InitStatus PndBarrelTrackFinder::ReInit ( )

privatevirtual

Reinitialisation

Definition at line 1837 of file PndBarrelTrackFinder.cxx.

                                        {

  return kSUCCESS;

}

void PndBarrelTrackFinder::RemoveCloneTracks ( )

private

Definition at line 1395 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), Bool_t, CalcPhi(), Double_t, fClonesRDiff, fClonesRDiffRel, fClonesXDiff, fClonesXDiffRel, fClonesYDiff, fClonesYDiffRel, fClonesZ_PDiff, fClonesZ_PDiffRel, fTracksVector, and fVerbose.

Referenced by Exec().

                                              {
  Bool_t printInfo = kFALSE;
  for ( Int_t itr1 = fTracksVector.size()-1 ; itr1 >= 0 ; itr1-- ) {
    for ( Int_t itr2 = fTracksVector.size()-1 ; itr2 > itr1 ; itr2-- ) {
      fClonesXDiff     ->Fill( fTracksVector[itr1].meanX  -fTracksVector[itr2].meanX);
      fClonesYDiff     ->Fill( fTracksVector[itr1].meanY  -fTracksVector[itr2].meanY);
      fClonesRDiff     ->Fill( fTracksVector[itr1].meanR  -fTracksVector[itr2].meanR);
      fClonesZ_PDiff   ->Fill( fTracksVector[itr1].meanZ_P-fTracksVector[itr2].meanZ_P);
      fClonesXDiffRel  ->Fill((fTracksVector[itr1].meanX  -fTracksVector[itr2].meanX  )/fTracksVector[itr2].meanR);
      fClonesYDiffRel  ->Fill((fTracksVector[itr1].meanY  -fTracksVector[itr2].meanY  )/fTracksVector[itr2].meanR);
      fClonesRDiffRel  ->Fill((fTracksVector[itr1].meanR  -fTracksVector[itr2].meanR  )/fTracksVector[itr2].meanR);
      fClonesZ_PDiffRel->Fill((fTracksVector[itr1].meanZ_P-fTracksVector[itr2].meanZ_P)/fTracksVector[itr2].meanZ_P);
    }
  }

  for ( Int_t itr1 = fTracksVector.size()-1 ; itr1 >= 0 ; itr1-- ) {
    Double_t tr1Phi  = CalcPhi(fTracksVector[itr1].meanX,fTracksVector[itr1].meanY);
    Int_t nofTr1Hits = fTracksVector[itr1].trackHitD.size();
    for ( Int_t itr2 = fTracksVector.size()-1 ; itr2 > itr1 ; itr2-- ) {
      Double_t tr2Phi = CalcPhi(fTracksVector[itr2].meanX,fTracksVector[itr2].meanY);
      
//       if (           (mcMag-reMag) <  0.3*mcMag &&
//                   (mcMag-reMag) > -0.3*mcMag &&
//         TMath::Abs(mcPhi-rePhi) <  5. &&
//         TMath::Abs(mcThe-reThe) <  5. ) {

//       if ( TMath::Abs(fTracksVector[itr1].meanR  -fTracksVector[itr2].meanR  ) < 0.1*fTracksVector[itr1].meanR &&
//         TMath::Abs(fTracksVector[itr1].meanZ_P-fTracksVector[itr2].meanZ_P) < fReasonableZ_PD+0.05*TMath::Abs(fTracksVector[itr1].meanZ_P) && 
//         TMath::Abs(tr1Phi-tr2Phi) < 0.1*TMath::RadToDeg() ) {
      if ( TMath::Abs(fTracksVector[itr1].meanX  -fTracksVector[itr2].meanX  ) < 0.2*TMath::Abs(fTracksVector[itr1].meanR) &&
           TMath::Abs(fTracksVector[itr1].meanY  -fTracksVector[itr2].meanY  ) < 0.2*TMath::Abs(fTracksVector[itr1].meanR) &&
           TMath::Abs(fTracksVector[itr1].meanR  -fTracksVector[itr2].meanR  ) < 0.3*fTracksVector[itr1].meanR &&
           TMath::Abs(fTracksVector[itr1].meanZ_P-fTracksVector[itr2].meanZ_P) < 0.4*TMath::Abs(fTracksVector[itr1].meanZ_P) && 
           TMath::Abs(tr1Phi-tr2Phi) < 0.2*TMath::RadToDeg() ) {
        if ( fVerbose > 0 || printInfo ) {
          cout << "seems that tracks " << itr1 << " (" << fTracksVector[itr1].meanX << "," << fTracksVector[itr1].meanY << "," << fTracksVector[itr1].meanR << "," << fTracksVector[itr1].meanZ_P << " )" << endl;
          cout << "              and " << itr2 << " (" << fTracksVector[itr2].meanX << "," << fTracksVector[itr2].meanY << "," << fTracksVector[itr2].meanR << "," << fTracksVector[itr2].meanZ_P << " ) ARE THE SAME!" << endl;
        }

        Int_t nofTr2Hits = fTracksVector[itr2].trackHitD.size();
        for ( Int_t ith2 = 0 ; ith2 < nofTr2Hits ; ith2++ ) {
          Bool_t hitExists = kFALSE;
          for ( Int_t ith1 = 0 ; ith1 < nofTr1Hits ; ith1++ ) {
            if ( fTracksVector[itr1].trackHitD[ith1] == fTracksVector[itr2].trackHitD[ith2] &&
                 fTracksVector[itr1].trackHitN[ith1] == fTracksVector[itr2].trackHitN[ith2] ) {
              hitExists = kTRUE;
              break;
            }
          }
          if ( hitExists ) continue;
          fTracksVector[itr1].trackHits.push_back(fTracksVector[itr2].trackHits[ith2]);
          fTracksVector[itr1].trackHitD.push_back(fTracksVector[itr2].trackHitD[ith2]);
          fTracksVector[itr1].trackHitN.push_back(fTracksVector[itr2].trackHitN[ith2]);
          //      cout << endl;
        }
        for ( size_t intr = itr2 ; intr < fTracksVector.size()-1 ; intr++ ) {
          fTracksVector[intr] = fTracksVector[intr+1];
        }
        fTracksVector.pop_back();
      }
    }
  }
}

void PndBarrelTrackFinder::RemoveHitFromPreviousHits ( Int_t  hitNo )

private

Definition at line 1553 of file PndBarrelTrackFinder.cxx.

References fHitVectDI, fHitVectHN, and fHitVector.

Referenced by Exec(), and MatchHitToHit().

                                                                {
  //    cout << endl << "There are " << fTracksVector.size() << " segments in the event" << endl;
  for ( size_t irh = hitNo ; irh < fHitVector.size()-1 ; irh++ ) {
    fHitVector[irh] = fHitVector[irh+1];
    fHitVectDI[irh] = fHitVectDI[irh+1];
    fHitVectHN[irh] = fHitVectHN[irh+1];
  }
  fHitVector.pop_back();
  fHitVectDI.pop_back();
  fHitVectHN.pop_back();
}

void PndBarrelTrackFinder::RemoveShortTracks ( )

private

Definition at line 1378 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), AddHitToPreviousHits(), and fTracksVector.

Referenced by Exec().

                                              {
  for ( Int_t itr = fTracksVector.size()-1 ; itr >= 0 ; itr-- ) {
    if ( TMath::Abs(fTracksVector[itr].meanZ_P) < 1.e6 && TMath::Abs(fTracksVector[itr].meanR) > 1. ) continue;
    for ( size_t ihit = 0 ; ihit < fTracksVector[itr].trackHits.size() ; ihit++ ) {
      AddHitToPreviousHits     (fTracksVector[itr].trackHits[ihit],
                                fTracksVector[itr].trackHitD[ihit],
                                fTracksVector[itr].trackHitN[ihit]);
    }
    for ( size_t intr = itr ; intr < fTracksVector.size()-1 ; intr++ ) {
      fTracksVector[intr] = fTracksVector[intr+1];
    }
    fTracksVector.pop_back();
  }
}

void PndBarrelTrackFinder::Reset ( )

private

Reset eventwise counters

Definition at line 1847 of file PndBarrelTrackFinder.cxx.

Referenced by Exec(), PndBarrelTrackFinder(), and ~PndBarrelTrackFinder().

                                 {
}

void PndBarrelTrackFinder::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 1712 of file PndBarrelTrackFinder.cxx.

References fSttParameters, and run.

                                            {
  
  // Get run and runtime database
  FairRunAna* run = FairRunAna::Instance();
  if ( ! run ) Fatal("SetParContainers", "No analysis run");

  FairRuntimeDb* db = run->GetRuntimeDb();
  if ( ! db ) Fatal("SetParContainers", "No runtime database");

  fSttParameters = (PndGeoSttPar*) db->getContainer("PndGeoSttPar");

  cout << "-I- PndBarrelTrackFinder::SetParContainers()." << endl;
}

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(), 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::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(), PndIdealTrackFinder::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 PndBarrelTrackFinder::UseMvdSttGem	(	const Bool_t	useMvd,
		const Bool_t	useStt,
		const Bool_t	useGem
	)

Definition at line 200 of file PndBarrelTrackFinder.cxx.

References fIncludeDet.

Referenced by barrelTrackFinder(), digi_BARREL(), reco_complete_newSTT(), tracks_BARREL_1000(), tracks_BARREL_1000sep_local(), and tracks_BARREL_1000sep_proof().

                                                                                                     {
  fIncludeDet[0] = useMvd;
  fIncludeDet[1] = useMvd;
  fIncludeDet[2] = useStt;
  fIncludeDet[3] = useGem;
};

Int_t PndBarrelTrackFinder::WriteTracks ( )

private

Definition at line 1461 of file PndBarrelTrackFinder.cxx.

References CAMath::Abs(), PndTrackCand::AddHit(), CalcPhi(), Double_t, fHitArrayName, fTracksVector, Pi, PndTrackCand::Sort(), and CAMath::Sqrt().

Referenced by Exec().

                                         {
  Int_t nofCreatedTracks = 0;

  for ( size_t itr = 0 ; itr < fTracksVector.size() ; itr++ ) {
    Double_t thisPhi = CalcPhi(fTracksVector[itr].meanX,
                               fTracksVector[itr].meanY);
    Double_t thisRad =         fTracksVector[itr].meanR;
    Double_t thisZ_P =         fTracksVector[itr].meanZ_P;
    if ( TMath::IsNaN(thisZ_P) ) thisZ_P = 2.e6;

    //       for ( Int_t ihit = 0 ; ihit < fTracksVector[itr].trackHitD.size() ; ihit++ )
    //  cout << fTracksVector[itr].trackHitD[ihit] << "." << fTracksVector[itr].trackHitN[ihit] << " " << flush;
    //       cout << endl;
    Double_t calcPt = thisRad/165.974;
    Double_t calcPz = thisZ_P*0.0058997;
    Double_t calcP  = TMath::Sqrt(calcPt*calcPt+calcPz*calcPz);
    Double_t calcThe = TMath::ACos(calcPz/calcP);
    
    Double_t calcPhi = thisPhi-TMath::Pi()/2.;
    Int_t calcChg = -1;
    if ( calcPhi < 0. ) calcPhi += TMath::Pi()*2.;
    Double_t firpPhi = CalcPhi(fTracksVector[itr].trackHits[0]->GetX(),fTracksVector[itr].trackHits[0]->GetY());
    if ( TMath::Abs(calcPhi-firpPhi) > TMath::Pi()/2. && TMath::Abs(calcPhi-firpPhi) < 3.*TMath::Pi()/2. ) {
      calcPhi = thisPhi+TMath::Pi()/2.;
      if ( calcPhi > 2.*TMath::Pi() ) calcPhi -= TMath::Pi()*2.;
      calcChg = 1;
    }
    if ( calcPhi > TMath::Pi() )
      calcPhi -= 2.*TMath::Pi();

    //    cout << "F " << thisPhi << " " << thisRad << " " << thisZ_P << "    " << calcPt << "  " << calcPhi << " " << calcPz << " " << calcP << " " << calcPhi << " " << calcThe << endl;//flush;
    
    //    cout << "                      MTP: " << calcP << " " << calcPhi << " " << calcThe << " and charge = " << calcChg << flush;
    
    TVector3 trackPosition(0.,0.,0.);
    TVector3 trackMomentum;
    trackMomentum.SetMagThetaPhi(calcP,calcThe,calcPhi);
    
    PndTrackCand* trackCand = new((*fBarrelTrackCandArray)[nofCreatedTracks]) PndTrackCand();

    std::vector<Int_t> trackHitsPerDet(4,0);
    for ( size_t ihit = 0 ; ihit < fTracksVector[itr].trackHitD.size() ; ihit++ ) {
      trackHitsPerDet[fTracksVector[itr].trackHitD[ihit]] += 1;
      Double_t tempPos = TMath::Sqrt(fTracksVector[itr].trackHits[ihit]->GetX()*fTracksVector[itr].trackHits[ihit]->GetX()+
                                     fTracksVector[itr].trackHits[ihit]->GetY()*fTracksVector[itr].trackHits[ihit]->GetY()+
                                     fTracksVector[itr].trackHits[ihit]->GetZ()*fTracksVector[itr].trackHits[ihit]->GetZ());
//       cout << ">>> " << fDetName[fTracksVector[itr].trackHitD[ihit]] << " hit will be written as "
//         << fDetType[fTracksVector[itr].trackHitD[ihit]] << "." << fTracksVector[itr].trackHitN[ihit] << endl;

      trackCand->AddHit(FairRootManager::Instance()->GetBranchId(fHitArrayName[fTracksVector[itr].trackHitD[ihit]].Data()),
                        fTracksVector[itr].trackHitN[ihit],
                        tempPos);
    }
    trackCand->Sort();

//     cout << "track hits in :" << flush;
//     for ( Int_t itch = 0 ; itch < trackCand->GetNHits() ; itch++ ) {
//       cout << trackCand->GetSortedHit(itch).GetDetId() << " " << flush;
//     }
//     cout << endl;

    FairTrackParP           firstPar(trackPosition,trackMomentum,
                                     TVector3(0.5, 0.5, 0.5),
                                     0.1*trackMomentum,
                                     calcChg,
                                     trackPosition,
                                     TVector3(1.,0.,0.),
                                     TVector3(0.,1.,0.));                                        
    TVector3 tempVect1(trackHitsPerDet[0],trackHitsPerDet[1],trackHitsPerDet[2]);
    TVector3 tempVect2(trackHitsPerDet[3],-1.,-1.); 
    
    //    FairTrackParP*            lastPar = new FairTrackParP(trackPosition,trackMomentum,
    FairTrackParP           lastPar (tempVect1,tempVect2,
                                     TVector3(0.5, 0.5, 0.5),
                                     0.1*trackMomentum,
                                     calcChg,
                                     trackPosition,
                                     TVector3(1.,0.,0.),
                                     TVector3(0.,1.,0.));                                        
    
    new((*fBarrelTrackArray)[nofCreatedTracks]) PndTrack(firstPar,lastPar,*trackCand);

    nofCreatedTracks++;

    //    cout << endl;
  }

  return nofCreatedTracks;
}

Member Data Documentation

TClonesArray* PndBarrelTrackFinder::fBarrelTrackArray

private

Output array of PndTrackCands.

Output array of global tracks

Definition at line 126 of file PndBarrelTrackFinder.h.

Referenced by Exec(), Finish(), Init(), PndBarrelTrackFinder(), and ~PndBarrelTrackFinder().

TClonesArray* PndBarrelTrackFinder::fBarrelTrackCandArray

private

Definition at line 127 of file PndBarrelTrackFinder.h.

Referenced by Exec(), Finish(), Init(), PndBarrelTrackFinder(), and ~PndBarrelTrackFinder().

TCanvas* PndBarrelTrackFinder::fCanvas

private

Definition at line 99 of file PndBarrelTrackFinder.h.

Referenced by Exec(), and Init().

TH1F* PndBarrelTrackFinder::fClonesRDiff

private

Definition at line 103 of file PndBarrelTrackFinder.h.

Referenced by Finish(), Init(), and RemoveCloneTracks().

TH1F* PndBarrelTrackFinder::fClonesRDiffRel

private

Definition at line 107 of file PndBarrelTrackFinder.h.

Referenced by Finish(), Init(), and RemoveCloneTracks().

TH1F* PndBarrelTrackFinder::fClonesXDiff

private

Definition at line 101 of file PndBarrelTrackFinder.h.

Referenced by Finish(), Init(), and RemoveCloneTracks().

TH1F* PndBarrelTrackFinder::fClonesXDiffRel

private

Definition at line 105 of file PndBarrelTrackFinder.h.

Referenced by Finish(), Init(), and RemoveCloneTracks().

TH1F* PndBarrelTrackFinder::fClonesYDiff

private

Definition at line 102 of file PndBarrelTrackFinder.h.

Referenced by Finish(), Init(), and RemoveCloneTracks().

TH1F* PndBarrelTrackFinder::fClonesYDiffRel

private

Definition at line 106 of file PndBarrelTrackFinder.h.

Referenced by Finish(), Init(), and RemoveCloneTracks().

TH1F* PndBarrelTrackFinder::fClonesZ_PDiff

private

Definition at line 104 of file PndBarrelTrackFinder.h.

Referenced by Finish(), Init(), and RemoveCloneTracks().

TH1F* PndBarrelTrackFinder::fClonesZ_PDiffRel

private

Definition at line 108 of file PndBarrelTrackFinder.h.

Referenced by Finish(), Init(), and RemoveCloneTracks().

TString PndBarrelTrackFinder::fDetName[4]

private

Definition at line 121 of file PndBarrelTrackFinder.h.

Referenced by Exec(), Init(), PndBarrelTrackFinder(), and PrintTracks().

Int_t PndBarrelTrackFinder::fDetType[4]

private

Definition at line 122 of file PndBarrelTrackFinder.h.

Referenced by Init(), and PndBarrelTrackFinder().

Bool_t PndBarrelTrackFinder::fDrawDetails

private

Definition at line 96 of file PndBarrelTrackFinder.h.

Referenced by DrawDetails(), Init(), and PndBarrelTrackFinder().

Bool_t PndBarrelTrackFinder::fDrawHistos

private

Definition at line 97 of file PndBarrelTrackFinder.h.

Referenced by DrawHistos(), Finish(), and PndBarrelTrackFinder().

TEllipse* PndBarrelTrackFinder::fEllipse[1000]

private

Definition at line 98 of file PndBarrelTrackFinder.h.

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

Int_t PndBarrelTrackFinder::fGeoType

private

Definition at line 176 of file PndBarrelTrackFinder.h.

Referenced by DrawTracks(), Exec(), ExtractMeanZ_PFromTrack(), Init(), MatchHitToHit(), MatchSkewedSttHitTT(), and PrintTracks().

TClonesArray* PndBarrelTrackFinder::fHitArray[4]

private

Definition at line 120 of file PndBarrelTrackFinder.h.

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

TString PndBarrelTrackFinder::fHitArrayName[4]

private

Definition at line 123 of file PndBarrelTrackFinder.h.

Referenced by Init(), PndBarrelTrackFinder(), and WriteTracks().

std::vector<Int_t> PndBarrelTrackFinder::fHitDetId

private

Definition at line 112 of file PndBarrelTrackFinder.h.

Referenced by Exec(), and PndBarrelTrackFinder().

std::vector<Int_t> PndBarrelTrackFinder::fHitDetNo

private

Definition at line 113 of file PndBarrelTrackFinder.h.

Referenced by Exec(), and PndBarrelTrackFinder().

std::vector<Int_t> PndBarrelTrackFinder::fHitVectDI

private

Definition at line 116 of file PndBarrelTrackFinder.h.

Referenced by AddHitToPreviousHits(), Exec(), MatchHitToHit(), PndBarrelTrackFinder(), and RemoveHitFromPreviousHits().

std::vector<Int_t> PndBarrelTrackFinder::fHitVectHN

private

Definition at line 117 of file PndBarrelTrackFinder.h.

Referenced by AddHitToPreviousHits(), Exec(), MatchHitToHit(), PndBarrelTrackFinder(), and RemoveHitFromPreviousHits().

std::vector<FairHit*> PndBarrelTrackFinder::fHitVector

private

Definition at line 115 of file PndBarrelTrackFinder.h.

Referenced by AddHitToPreviousHits(), Exec(), MatchHitToHit(), PndBarrelTrackFinder(), and RemoveHitFromPreviousHits().

Bool_t PndBarrelTrackFinder::fIncludeDet[4]

private

Definition at line 119 of file PndBarrelTrackFinder.h.

Referenced by Exec(), Init(), PndBarrelTrackFinder(), PrintTracks(), and UseMvdSttGem().

Double_t PndBarrelTrackFinder::fMaximalDist

private

Definition at line 133 of file PndBarrelTrackFinder.h.

Referenced by MatchHitWithZInfoTT(), MatchParallelSttHitTT(), and PndBarrelTrackFinder().

Double_t PndBarrelTrackFinder::fMaximalPhiDiff

private

Definition at line 138 of file PndBarrelTrackFinder.h.

Referenced by PndBarrelTrackFinder().

Double_t PndBarrelTrackFinder::fMaximalRadDiff

private

Definition at line 137 of file PndBarrelTrackFinder.h.

Referenced by PndBarrelTrackFinder().

Double_t PndBarrelTrackFinder::fMaximalZ_PD

private

Definition at line 135 of file PndBarrelTrackFinder.h.

Referenced by MatchHitWithZInfoTT(), MatchParallelSttHitTT(), and PndBarrelTrackFinder().

Double_t PndBarrelTrackFinder::fReasonableDist

private

Definition at line 134 of file PndBarrelTrackFinder.h.

Referenced by ExtractMeanRPhiFromTrack(), MatchHitWithZInfoTT(), MatchParallelSttHitTT(), and PndBarrelTrackFinder().

Double_t PndBarrelTrackFinder::fReasonableZ_PD

private

Definition at line 136 of file PndBarrelTrackFinder.h.

Referenced by ExtractMeanZ_PFromTrack(), MatchHitToHit(), MatchHitWithZInfoTT(), MatchSkewedSttHitTT(), and PndBarrelTrackFinder().

TLatex* PndBarrelTrackFinder::fStageDesc

private

Definition at line 100 of file PndBarrelTrackFinder.h.

Referenced by Exec(), and Init().

PndGeoSttPar* PndBarrelTrackFinder::fSttParameters

private

Definition at line 175 of file PndBarrelTrackFinder.h.

Referenced by Init(), PndBarrelTrackFinder(), and SetParContainers().

Int_t PndBarrelTrackFinder::fTNofEvents

private

event counter

Event counter

Definition at line 130 of file PndBarrelTrackFinder.h.

Referenced by Exec(), Finish(), and PndBarrelTrackFinder().

Int_t PndBarrelTrackFinder::fTNofTracks

private

Definition at line 131 of file PndBarrelTrackFinder.h.

Referenced by Exec(), Finish(), and PndBarrelTrackFinder().

std::vector<TrackBasis> PndBarrelTrackFinder::fTracksVector

private

Definition at line 110 of file PndBarrelTrackFinder.h.

Referenced by DrawTracks(), Exec(), ExtractMeanRPhiFromTrack(), ExtractMeanZ_PFromTrack(), HitBelongsToTrack(), MatchHitToHit(), MatchHitWithZInfoTT(), MatchParallelSttHitTT(), MatchSkewedSttHitTT(), PndBarrelTrackFinder(), PrintTracks(), RemoveCloneTracks(), RemoveShortTracks(), and WriteTracks().

TClonesArray* PndBarrelTrackFinder::fTubeArray

private

from parameters array of PndSttTube

Definition at line 174 of file PndBarrelTrackFinder.h.

Referenced by DrawTracks(), Exec(), ExtractMeanZ_PFromTrack(), Init(), MatchHitToHit(), MatchSkewedSttHitTT(), PndBarrelTrackFinder(), and PrintTracks().

---

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

• PndBarrelTrackFinder.h
• PndBarrelTrackFinder.cxx