THaVDCChamber.h

Go to the documentation of this file.

#ifndef Podd_THaVDCChamber_h_
#define Podd_THaVDCChamber_h_
 
//                                                                           //
// THaVDCChamber                                                             //
//                                                                           //
 
#include "THaSubDetector.h"
#include "TClonesArray.h"
#include "THaVDCPlane.h"
#include "THaVDCPoint.h"
#include <cassert>
 
class THaEvData;
 
// Coordinates of a point in this chamber
class PointCoords_t {
public:
  Double_t x;
  Double_t y;
  Double_t theta;  // TRANSPORT x' = dx/dz
  Double_t phi;    // TRANSPORT y' = dy/dz
};
 
class THaVDCChamber : public THaSubDetector {
 
public:
 
  explicit THaVDCChamber( const char* name="", const char* description="",
                          THaDetectorBase* parent = nullptr );
  virtual ~THaVDCChamber();
 
  virtual void    Clear( Option_t* opt="" );    // Reset event-by-event data
  virtual Int_t   Decode( const THaEvData& evData );
  virtual Int_t   CoarseTrack();          // Find clusters & estimate track
  virtual Int_t   FineTrack();            // More precisely calculate track
  virtual EStatus Init( const TDatime& date );
  virtual void    SetDebug( Int_t level );
 
  PointCoords_t   CalcDetCoords( const THaVDCCluster* u,
                                 const THaVDCCluster* v ) const;
 
  // Get and Set Functions
  THaVDCPlane*    GetUPlane()      const { return fU; }
  THaVDCPlane*    GetVPlane()      const { return fV; }
  Int_t           GetNPoints()     const { return fPoints->GetLast()+1; }
  TClonesArray*   GetPoints()      const { return fPoints; }
  Double_t        GetSpacing()     const { return fSpacing;}
  THaVDCPoint*    GetPoint( Int_t i ) const
    { assert( i>=0 && i<GetNPoints() );
      return static_cast<THaVDCPoint*>( fPoints->UncheckedAt(i) ); }
  Double_t        GetZ()           const { return fU->GetZ(); }
 
protected:
 
  // Subdetectors
  THaVDCPlane*  fU;           // The U plane
  THaVDCPlane*  fV;           // The V plane
 
  // Event data
  TClonesArray* fPoints;      // Pairs of possibly matching U and V clusters
 
  // Geometry
  Double_t fSpacing;          // Space between U & V planes (m)
  Double_t fSin_u;            // Trig functions for the U plane axis angle
  Double_t fCos_u;            //
  Double_t fSin_v;            // Trig functions for the V plane axis angle
  Double_t fCos_v;            //
  Double_t fInv_sin_vu;       // 1/Sine of the difference between the
                              // V axis angle and the U axis angle
 
  void  ApplyTimeCorrection(); // Correct hit drift times by common offset
  void  FindClusters();       // Find clusters in U and V planes
  void  FitTracks();          // Fit local tracks for each cluster
  Int_t MatchUVClusters();    // Match clusters in U with clusters in V
  Int_t CalcPointCoords() const;
 
  Double_t UVtoX( Double_t u, Double_t v ) const;
  Double_t UVtoY( Double_t u, Double_t v ) const;
 
  virtual Int_t DefineVariables( EMode mode = kDefine );
 
  ClassDef(THaVDCChamber,0)   // VDC chamber (pair of a U and a V plane)
};
 
//_____________________________________________________________________________
inline Double_t THaVDCChamber::UVtoX( Double_t u, Double_t v ) const
{
  return (u*fSin_v - v*fSin_u) * fInv_sin_vu;
}
 
//_____________________________________________________________________________
inline Double_t THaVDCChamber::UVtoY( Double_t u, Double_t v ) const
{
  return (v*fCos_u - u*fCos_v) * fInv_sin_vu;
}
 
 
#endif