#include <PndLmdSimplePixelClusterFinder.h>

Inheritance diagram for PndLmdSimplePixelClusterFinder:

Public Member Functions
	PndLmdSimplePixelClusterFinder (Int_t verbose=0)

std::vector< std::vector< Int_t > >	GetClusters (std::vector< PndSdsDigiPixel > hits)

virtual void	SetMaxCols (Int_t col)

virtual void	SetMaxRows (Int_t row)

virtual void	SetRadius (Double_t rad)

void	Print ()

void	Print (std::vector< PndSdsDigiPixel > hits)

void	PrintResult (std::vector< std::vector< Int_t > > clusters)

void	SetVerbose (Int_t level)

	ClassDef (PndSdsPixelClusterFinder, 2)

Public Attributes
std::vector< PndSdsDigiPixel >	fHits

Int_t	fVerbose

PndSdsChargeConversion *	fChargeConverter

Protected Attributes
Int_t	fcols

Int_t	frows

Double_t	fradius

Detailed Description

PndLmdSimplePixelClusterFinder.h Simple cluster finder for pixel detectors.

It takes the first hit in the DigiArray and moves it into a new cluster. Then it looks through the remaining hits in the DigiArray if there are hits which are in the area around the first hit and moves them into the cluster. The radius is given by the first parameter. Then it takes the new hits in the cluster and looks if there are other hits which are in the range of these hits and moves them into the cluster, too. This process is repeated until no hits around the cluster are found The the first hit of the remaining DigiHits is taken to create a new hit and the cluster search process is started again.

Definition at line 18 of file PndLmdSimplePixelClusterFinder.h.

Constructor & Destructor Documentation

PndLmdSimplePixelClusterFinder::PndLmdSimplePixelClusterFinder ( Int_t verbose = 0 )

Definition at line 8 of file PndLmdSimplePixelClusterFinder.cxx.

References PndSdsPixelClusterFinder::fChargeConverter, PndSdsSimplePixelClusterFinder::fcols, PndSdsSimplePixelClusterFinder::fradius, PndSdsSimplePixelClusterFinder::frows, PndSdsPixelClusterFinder::fVerbose, PndSdsPixelDigiPar::GetChargeConvMethod(), PndSdsTotDigiPar::GetChargingTime(), PndSdsTotDigiPar::GetClockFrequency(), PndSdsPixelDigiPar::GetClustRadius(), PndSdsTotDigiPar::GetConstCurrent(), PndSdsPixelDigiPar::GetFECols(), PndSdsPixelDigiPar::GetFERows(), PndSdsPixelDigiPar::GetNoise(), PndSdsPixelDigiPar::GetThreshold(), rtdb, and verbose.

     : PndSdsSimplePixelClusterFinder() {
   fVerbose = verbose;
   FairRun* ana = FairRun::Instance();
   FairRuntimeDb* rtdb = ana->GetRuntimeDb();
   PndSdsPixelDigiPar* fDigiPar =
       (PndSdsPixelDigiPar*)(rtdb->getContainer("LMDPixelDigiPar"));
   PndSdsTotDigiPar* fTotDigiPar =
       (PndSdsTotDigiPar*)(rtdb->getContainer("LMDPixelTotDigiPar"));
   if (fDigiPar->GetChargeConvMethod() == 0) {
     if (fVerbose > 0)
       std::cout
           << "Info in <PndLmdSimplePixelClusterFinder>: ideal charge conversion"
           << std::endl;
     fChargeConverter = new PndSdsIdealChargeConversion(fDigiPar->GetNoise());
   } else if (fDigiPar->GetChargeConvMethod() == 1) {
     if (fVerbose > 0)
       std::cout
           << "Info in <PndLmdSimplePixelClusterFinder>: TOT charge conversion"
           << std::endl;
     fChargeConverter = new PndSdsTotChargeConversion(
         fTotDigiPar->GetChargingTime(), fTotDigiPar->GetConstCurrent(),
         fDigiPar->GetThreshold(), fTotDigiPar->GetClockFrequency(), fVerbose);
   } else {
     std::cout << "Error in <PndLmdSimplePixelClusterFinder>: charge conversion "
                  "method not defined, use ideal converter..."
               << std::endl;
     fChargeConverter = new PndSdsIdealChargeConversion(fDigiPar->GetNoise());
   }
   fcols = fDigiPar->GetFECols();
   frows = fDigiPar->GetFERows();
   fradius = fDigiPar->GetClustRadius();
 };

Member Function Documentation

PndSdsPixelClusterFinder::ClassDef	(	PndSdsPixelClusterFinder	,
		2
	)

inherited

std::vector< std::vector< Int_t > > PndSdsSimplePixelClusterFinder::GetClusters ( std::vector< PndSdsDigiPixel > hits )

virtualinherited

Main method which searches for the clusters. It returns a matrix of ints where a column corresponds to a cluster and an integer to a hit in the DigiArray

Implements PndSdsPixelClusterFinder.

Definition at line 6 of file PndSdsSimplePixelClusterFinder.cxx.

References PndSdsChargeConversion::DigiValueToCharge(), PndSdsPixelClusterFinder::fChargeConverter, PndSdsPixelClusterFinder::fHits, PndSdsPixelClusterFinder::fVerbose, hits, i, PndSdsSimplePixelClusterFinder::IsInRange(), PndSdsSimplePixelClusterFinder::MoveHit(), PndSdsPixelClusterFinder::PrintResult(), and push_back().

 {
         fHits = hits;
 //      if (fVerbose > 1) Print();
         std::vector<Int_t> posHits;
         for (UInt_t i = 0; i < fHits.size(); i++) posHits.push_back(i);
         std::vector< std::vector< Int_t> > result;
         Int_t sizeTempHits = posHits.size();
         //Int_t actHit = 0;
         while (sizeTempHits != 0){
                 std::vector<Int_t> tempInt;
                 if (fChargeConverter->DigiValueToCharge(fHits[posHits[0]]) == 0){
                         MoveHit(&posHits,0);
                         sizeTempHits = posHits.size();
                 }
                 else {
                         tempInt.push_back(MoveHit(&posHits,0));
                         result.push_back(tempInt);
                         Int_t sizeResultI = (result.end()-1)->size(); // size of last vector in array
                         for (Int_t i = 0; i < sizeResultI; i++){
                                 sizeTempHits = posHits.size();
                                 for (Int_t j = 0; j < sizeTempHits; j++){
                                         if (fChargeConverter->DigiValueToCharge(fHits[posHits[j]]) == 0){
                                                                 MoveHit(&posHits,j);
                                                                 if (fVerbose > 1)
                                                                         std::cout << "Charge too low!" << std::endl;
                                                                 j--;
                                         }
                                         else{
                                                 if (fHits[(*(result.end()-1))[i]].GetSensorID() == fHits[posHits[j]].GetSensorID()){
                                                         if (IsInRange(fHits[(*(result.end()-1))[i]], fHits[posHits[j]])) {
                                                                         (result.end()-1)->push_back(MoveHit(&posHits,j));
                                                                         j--;
                                                                         if (fVerbose > 1)
                                                                                 std::cout << "Hit added to cluster: " << result.size()-1 << std::endl;
                                                         }
                                                 }
 
                                         }
                                         sizeTempHits = posHits.size();
                                 }
                                 sizeResultI = (result.end()-1)->size();
                         }
                 }
         }
         if (fVerbose > 1) PrintResult(result);
         return result;
 }

void PndSdsPixelClusterFinder::Print ( )

inherited

Definition at line 20 of file PndSdsPixelClusterFinder.cxx.

References PndSdsPixelClusterFinder::fHits.

Referenced by PndSdsIdealPixelClusterFinder::GetClusters(), and PndSdsPixelClusterFinder::PrintResult().

 {
         Print(fHits);
 }

void PndSdsPixelClusterFinder::Print ( std::vector< PndSdsDigiPixel > hits )

inherited

Definition at line 12 of file PndSdsPixelClusterFinder.cxx.

References i.

 {
         std::cout << "--------------" << std::endl;
         for (UInt_t i = 0; i < hits.size(); i++){
                 std::cout << hits[i] << std::endl;
         }
 }

void PndSdsPixelClusterFinder::PrintResult ( std::vector< std::vector< Int_t > > clusters )

inherited

Definition at line 25 of file PndSdsPixelClusterFinder.cxx.

References PndSdsPixelClusterFinder::fHits, i, and PndSdsPixelClusterFinder::Print().

Referenced by PndSdsIdealPixelClusterFinder::GetClusters(), and PndSdsSimplePixelClusterFinder::GetClusters().

 {
         for (UInt_t i = 0; i < clusters.size(); i++){
                 std::cout << "PndSdsPixelClusterFinder::PrintResult(): Cluster " << i << ": " << std::endl;
                 for (UInt_t j = 0; j < clusters[i].size();j++)
                         fHits[clusters[i][j]].Print();
         }
 }