PndGemMatchHits Class Reference

#include <PndGemMatchHits.h>

Inheritance diagram for PndGemMatchHits:

Public Member Functions
	PndGemMatchHits ()
	PndGemMatchHits (Int_t iVerbose)
	PndGemMatchHits (const char *name, Int_t iVerbose)
virtual	~PndGemMatchHits ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()

Private Member Functions
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
void	Reset ()
	ClassDef (PndGemMatchHits, 1)

Private Attributes
PndGemDigiPar *	fDigiPar
TClonesArray *	fPoints
TClonesArray *	fHits
Int_t	fNHits
Int_t	fNMatchedHits
Int_t	fNFakeHits
Int_t	fNMultiHits

Detailed Description

Definition at line 32 of file PndGemMatchHits.h.

Constructor & Destructor Documentation

PndGemMatchHits::PndGemMatchHits

(

)

Default constructor

Definition at line 46 of file PndGemMatchHits.cxx.

References fDigiPar, fHits, fPoints, and Reset().

                                 : FairTask("GEM MatchHits", 1) {
  fDigiPar     = NULL;
  fPoints      = NULL;
  fHits        = NULL;
  Reset();
}

PndGemMatchHits::PndGemMatchHits

(

Int_t

iVerbose

)

Standard constructor

Definition at line 57 of file PndGemMatchHits.cxx.

References fDigiPar, fHits, fPoints, and Reset().

  : FairTask("GEM MatchHitsr", iVerbose) { 
  fDigiPar     = NULL;
  fPoints      = NULL;
  fHits        = NULL;
  Reset();
}

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

Constructor with name

Definition at line 69 of file PndGemMatchHits.cxx.

References fDigiPar, fHits, fPoints, and Reset().

  : FairTask(name, iVerbose) { 
  fDigiPar     = NULL;
  fPoints      = NULL;
  fHits        = NULL;
  Reset();
}

PndGemMatchHits::~PndGemMatchHits ( )

virtual

Destructor

Definition at line 81 of file PndGemMatchHits.cxx.

References fDigiPar, and Reset().

                                  { 
  if ( fDigiPar)   delete fDigiPar;
  Reset();
}

Member Function Documentation

PndGemMatchHits::ClassDef ( PndGemMatchHits  , 1    )

private

void PndGemMatchHits::Exec ( Option_t *  opt )

virtual

Execution

Definition at line 88 of file PndGemMatchHits.cxx.

References CAMath::Abs(), Bool_t, Double_t, fHits, fNFakeHits, fNHits, fNMatchedHits, fNMultiHits, fPoints, fVerbose, PndGemHit::GetDp(), PndGemHit::GetDr(), PndGemMCPoint::GetSensorId(), nHits, Pi, PndGemHit::SetBotIndex(), CAMath::Sqrt(), and CAMath::Tan().

                                    {
  Int_t nofPoints = fPoints->GetEntriesFast();
  Int_t nofHits = fHits->GetEntriesFast();

  if ( fVerbose )
    cout << "PndGemMatchHits::Exec() with " << nofPoints << " points and " << nofHits << " hits." << endl;

  Int_t nHits = 0;
  Int_t nMatchedHits = 0;
  Int_t nFakeHits = 0;
  Int_t nMultiHits = 0;

  vector<Double_t> pointZ;
  vector<Double_t> pointR;
  vector<Double_t> pointP;
  for ( Int_t iPoint = 0 ; iPoint < nofPoints ; iPoint++ ) {
    PndGemMCPoint* currentPndGemMCPoint = (PndGemMCPoint*)fPoints->At(iPoint);
    
    Double_t pointX = currentPndGemMCPoint->GetX();
    Double_t pointY = currentPndGemMCPoint->GetY();

    if ( fVerbose > 1 )
      cout << " .... " << pointX << "  " << pointY << " " << currentPndGemMCPoint->GetZ() << endl;

    Double_t phiAValue = TMath::ATan(pointX/pointY);
    if ( pointY < 0 ) phiAValue += TMath::Pi();
    else if ( pointX < 0 ) phiAValue +=  2.*TMath::Pi();

    pointZ.push_back(currentPndGemMCPoint->GetZ());
    pointR.push_back(TMath::Sqrt(pointX*pointX+pointY*pointY));
    pointP.push_back(phiAValue);
    
    if ( fVerbose > 1 )
      cout << "point " << iPoint << ", sensor Id = " << currentPndGemMCPoint->GetSensorId() << " (" << pointZ[pointZ.size()-1] << "," << pointR[pointR.size()-1] << "," << pointP[pointP.size()-1] << ")" << endl;
  }

  for ( Int_t iHit = 0 ; iHit < nofHits ; iHit++ ) {
    PndGemHit* currentPndGemHit = (PndGemHit*)fHits->At(iHit);

    Double_t hitX = currentPndGemHit->GetX();
    Double_t hitY = currentPndGemHit->GetY();
    Double_t hitZ = currentPndGemHit->GetZ();

    if ( fVerbose > 1 )
      cout << " .... " << hitX << "  " << hitY << " " << currentPndGemHit->GetZ() << endl;

    Double_t hitP = TMath::ATan(hitX/hitY);
    if ( hitY < 0 ) hitP += TMath::Pi();
    else if ( hitX < 0 ) hitP +=  2.*TMath::Pi();
    Double_t hitR = TMath::Sqrt(hitX*hitX+hitY*hitY);
    
    if ( fVerbose > 1 ) {
      cout << "-------------------------------------------" << endl;
      cout << "hit " << iHit << " (" << hitZ << "," << hitR << "," << hitP << ")" << endl;
    }

    Int_t matchPoint = -1;
    Double_t closestDistance = 1000.;
    Bool_t multiHit = kFALSE;
    for ( size_t iPoint = 0 ; iPoint < pointZ.size() ; iPoint++ ) {
      //     if ( fVerbose > 1 )
      if ( TMath::Abs(pointZ[iPoint]-hitZ) > currentPndGemHit->GetDz() ) { /*cout << "FAILED Z" << endl;*/ continue; }
      if ( TMath::Abs(pointR[iPoint]-hitR) > currentPndGemHit->GetDr()*TMath::Sqrt(3.) ) { /*cout << "FAILED R" << endl;*/ continue; }
    
      if ( hitP < 1.0 && pointP[iPoint] > 5.5 ) hitP += 2.*TMath::Pi();
      if ( hitP > 5.5 && pointP[iPoint] < 1.0 ) hitP -= 2.*TMath::Pi();

          if ( TMath::Abs(TMath::Tan(pointP[iPoint]-hitP)) > currentPndGemHit->GetDp()*TMath::Sqrt(3.)/hitR ) { /*cout << "FAILED PHI" << endl;*/ continue; }
      if ( fVerbose > 1 ) {
        cout << "matched with " << " (" << pointZ[iPoint] << "," << pointR[iPoint] << "," << pointP[iPoint] << ") " << endl;
        cout << "DP = " << currentPndGemHit->GetDp() << " after transf. = " << currentPndGemHit->GetDp()*TMath::Sqrt(3.)/hitR << " while p_p = " << pointP[iPoint] << " " << " h_p = " << hitP << endl;
        cout << "PASSED WITH POINT " << iPoint << endl;
      }
      Double_t distance = TMath::Sqrt((pointR[iPoint]-hitR)*(pointR[iPoint]-hitR)/currentPndGemHit->GetDr()/currentPndGemHit->GetDr()+
                                      (pointP[iPoint]-hitP)*(pointP[iPoint]-hitP)/currentPndGemHit->GetDp()/currentPndGemHit->GetDp());
      if ( matchPoint != -1 ) multiHit = kTRUE;
      if ( distance > closestDistance ) continue;
      closestDistance = distance;
      matchPoint = iPoint;   
    }
    currentPndGemHit->SetRefIndex(matchPoint);
    if ( matchPoint != -1 ) {
      PndGemMCPoint* matchp = (PndGemMCPoint*)fPoints->At(matchPoint);
      currentPndGemHit->SetBotIndex(matchp->GetTrackID());
    }
    else {
      currentPndGemHit->SetBotIndex(-1);
    }

    nHits++;
    if ( matchPoint != -1 ) nMatchedHits++;
    else nFakeHits++;
    if ( multiHit ) nMultiHits++;
  }

  fNHits        += nHits;
  fNMatchedHits += nMatchedHits;
  fNFakeHits    += nFakeHits;
  fNMultiHits   += nMultiHits;
  
  if ( fVerbose ) {
    cout << "************PndGemMatchHits**************" << endl;
    cout << " Number of all hits " << nHits << endl;
    cout << " Number of matched hits " << nMatchedHits << " -> " << 100.*(Double_t)nMatchedHits/(Double_t)nHits << endl;
    cout << " Number of fake hits " << nFakeHits << " -> " << 100.*(Double_t)nFakeHits/(Double_t)nHits << endl;
    cout << " Number of multi hits " << nMultiHits << " -> " << 100.*(Double_t)nMultiHits/(Double_t)nHits << endl;
    cout << "*****************************************" << endl;
  }

}

void PndGemMatchHits::Finish ( )

virtual

Finish

Definition at line 254 of file PndGemMatchHits.cxx.

References Double_t, fNFakeHits, fNHits, fNMatchedHits, and fNMultiHits.

                             {
  cout << "************PndGemMatchHits summary**************" << endl;
  cout << " Number of all hits " << fNHits << endl;
  cout << " Number of matched hits " << fNMatchedHits << " -> " << 100.*(Double_t)fNMatchedHits/(Double_t)fNHits << endl;
  cout << " Number of fake hits " << fNFakeHits << " -> " << 100.*(Double_t)fNFakeHits/(Double_t)fNHits << endl;
  cout << " Number of multi hits " << fNMultiHits << " -> " << 100.*(Double_t)fNMultiHits/(Double_t)fNHits << endl;
  cout << "*************************************************" << endl;
}

InitStatus PndGemMatchHits::Init ( )

privatevirtual

Intialisation

Definition at line 221 of file PndGemMatchHits.cxx.

References fHits, and fPoints.

                                 {

  // Get input array 
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) Fatal("Init", "No FairRootManager");
  fPoints = (TClonesArray*) ioman->GetObject("GEMPoint");

  fHits = (TClonesArray*) ioman->GetObject("GEMHit");

  return kSUCCESS;

}

InitStatus PndGemMatchHits::ReInit ( )

privatevirtual

Reinitialisation

Definition at line 238 of file PndGemMatchHits.cxx.

                                   {

  return kSUCCESS;

}

void PndGemMatchHits::Reset ( )

private

Reset eventwise counters

Definition at line 248 of file PndGemMatchHits.cxx.

References fNFakeHits, fNHits, fNMatchedHits, and fNMultiHits.

Referenced by PndGemMatchHits(), and ~PndGemMatchHits().

                            {
  fNHits = fNMatchedHits = fNFakeHits = fNMultiHits = 0;
}

void PndGemMatchHits::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 203 of file PndGemMatchHits.cxx.

References fDigiPar, 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");

  // Get GEM digitisation parameter container
  fDigiPar = (PndGemDigiPar*)(db->getContainer("PndGemDetectors"));

}

Member Data Documentation

PndGemDigiPar* PndGemMatchHits::fDigiPar

private

Definition at line 56 of file PndGemMatchHits.h.

Referenced by PndGemMatchHits(), SetParContainers(), and ~PndGemMatchHits().

TClonesArray* PndGemMatchHits::fHits

private

Input array of PndGemPoint

Definition at line 58 of file PndGemMatchHits.h.

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

Int_t PndGemMatchHits::fNFakeHits

private

Definition at line 62 of file PndGemMatchHits.h.

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

Int_t PndGemMatchHits::fNHits

private

Input array of PndGemHit

Definition at line 60 of file PndGemMatchHits.h.

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

Int_t PndGemMatchHits::fNMatchedHits

private

Definition at line 61 of file PndGemMatchHits.h.

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

Int_t PndGemMatchHits::fNMultiHits

private

Definition at line 63 of file PndGemMatchHits.h.

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

TClonesArray* PndGemMatchHits::fPoints

private

Definition at line 57 of file PndGemMatchHits.h.

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

---

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

• PndGemMatchHits.h
• PndGemMatchHits.cxx