THaTrackingDetector.cxx

Go to the documentation of this file.

//*-- Author :    Ole Hansen   7-Sep-00
 
//
// THaTrackingDetector
//
// Abstract base class for a generic tracking detector. 
//
// This is a special THaSpectrometerDetector that is capable of
// finding particle tracks. It does not need to find track vertices
// (target coordinates) - this is usually done by the enclosing apparatus.
//
// Note: The FPP is NOT a "tracking detector" in this scheme because,
// with the present FPP algorithms, it is not used for calculating the 
// fp coordinates. 
//
 
#include "TClonesArray.h"
#include "THaTrackingDetector.h"
#include "THaPIDinfo.h"
#include "THaTrack.h"
#include "THaSpectrometer.h"
#include "TError.h"
 
ClassImp(THaTrackingDetector)
 
//_____________________________________________________________________________
THaTrackingDetector::THaTrackingDetector( const char* name, 
                                          const char* description,
                                          THaApparatus* apparatus )
  : THaSpectrometerDetector(name,description,apparatus)
{
  // Constructor
 
}
 
//_____________________________________________________________________________
THaTrackingDetector::THaTrackingDetector( )
  : THaSpectrometerDetector( )
{
  // Constructor for ROOT I/O only
}
 
//_____________________________________________________________________________
THaTrackingDetector::~THaTrackingDetector() = default;
 
//_____________________________________________________________________________
THaTrack* THaTrackingDetector::AddTrack( TClonesArray& tracks,
                                         Double_t x, Double_t y,
                                         Double_t theta, Double_t phi,
                                         THaTrackID* ID )
{
  // Add a track with the given parameters to the array of tracks
  // Returns the index of the created track in the TClonesArray.
 
  Int_t i = tracks.GetLast()+1;
  THaPIDinfo* pid = nullptr;
  auto* spect = static_cast<THaSpectrometer*>( GetApparatus() );
 
  if( spect ) {
    if( spect->IsPID() ) {
      // Create new PIDinfo and vertex objects if necessary
 
      UInt_t ndet = spect->GetNpidDetectors();
      UInt_t npart = spect->GetNpidParticles();
      TClonesArray& c = *spect->GetTrackPID();
 
      // caution: c[i] creates an empty object at i if slot i was empty,
      // so use At(i) to check.
 
      if( i > c.GetLast() || !c.At(i) ) new(c[i]) THaPIDinfo(ndet, npart);
      pid = static_cast<THaPIDinfo*>( c.At(i) );
    }
  } else if( fDebug>0 ) {
    ::Warning("THaTrackingDetector::AddTrack", "No spectrometer defined for "
              "detector %s. Track has no PID and vertex info.", GetName());
  }
 
  // Create the track
 
  return new( tracks[i] ) THaTrack( x, y, theta, phi, this, ID, pid );
 
}