PREAN/TFndTrack.h

// @(#)fROOT/PREAN:$Name:  $:$Id: TFndTrack.h,v 1.9 2007/10/02 09:00:08 Diego_Faso Exp $
// Author: Diego Faso <mailto:faso@to.infn.it>, 2007/01/08

#ifndef FROOT_TFndTrack
#define FROOT_TFndTrack

//                             TFndTrack                                //
//                                                          //
// Vertices are needed by tracks, so they are defined here: //
//  see TFndVertex                                          //
//                                                          //

#include "FROOT.h" // very important (procedure definitions)
#include "FIN_PHYS.h"

#include <Riostream.h>
#include <TObject.h>
#include <TString.h>
#include "TDatime.h"
#include "TTree.h"
#include "TFile.h"

class TFndVertex: public TObject{

  // --- data-members
 protected:
  Int_t fPID; // particle ID (geant convention) [see FIN_PHYS namespace]
  Float_t fPosition[3]; // MRS position [x,y,z]
  Float_t fDistFromPrimary; // distance from K-(or K+) vertex
  Float_t fZtollerance; // z-tollerance

 public:
  TFndVertex();
  ~TFndVertex();

  void Reset();
  // ---
  friend class TFndPrean;
  friend class TFndPreanHistos;
  friend class TFndEvd;
  ClassDef(TFndVertex,1) // FINUDA vertices definition

};


class TFndTrack: public TObject{
  
 public:
  enum Trk_Charge_e {
    E_Trk_Neut = 0,
    E_Trk_Neg = -1,
    E_Trk_Pos =  1,
  };

  enum Trk_Vect_e { // momentum vector definition
    E_Trk_X = 0,
    E_Trk_Y = 1,
    E_Trk_Z = 2,
    E_Trk_T = 3,
    E_Trk_cosX = 4,
    E_Trk_cosY = 5,
    E_Trk_cosZ = 6,
    E_Trk_Vect_End = 7,
  };

  enum Trk_Helix_Geo_e {
    E_Hel_Xc = 0, // X of the center of the helix 
    E_Hel_Yc = 1, // Y of the center of the helix
    E_Hel_R  = 2, // Radius of curvature
    E_Hel_Geo_End  = 3, // Used for arrays
  };

  enum Trk_Helix_Mom_e {
    E_Hel_Mom = 0, // Momentum (Gev/c)
    E_Hel_Lam = 1, // Lambda angle
    E_Hel_Phi = 2, // Phi angle
    E_Hel_Mom_End  = 3, // Used for arrays
  };

  enum Trk_FitError_e {
    E_Trk_Fit_Success = 0,       // Spline fit successfull
    E_Trk_Fit_FewPoints = -2,    // Less than 4 points for a long track
    E_Trk_Fit_BadMom = 1,        // Bad momentum reconstruction IERMOM=1
    E_Trk_Fit_UndefCharge = 10,  // Charge not defined. ICHARG=0
    E_Trk_Fit_UndefMom = 11,     // Momentum not defined. PMAG=0
    //  (>0) Numerical error. Value of IERMAT
  };
  
  enum Trk_ChiSquared_e {
    E_Trk_ChiSq_DEV = 0, // Spline (DEVNEW**2) or numerical fitting algorithm
    E_Trk_ChiSq_RES = 1, // Spline fit residuals in x,y plane (RESDYX)
    E_Trk_ChiSq_LOSTstb   = 2, // Chi2 for left/right ambiguity from LOST
    E_Trk_ChiSq_STEREOstb = 3, // Chi2 from stereo straw fitting (FMIN)
    E_Trk_ChiSq_End = 4, // Used for arrays
  };

  enum Trk_BackTrack_e {
    E_Trk_Bcktrk_NormAng          = 0, // Track incidence angle relative to the normal to the target
    E_Trk_Bcktrk_SurfApproach     = 1, // length from the tgt surface to point of max approach to vertex
    E_Trk_Bcktrk_MaxApproach      = 2, // Distance of max. approach of the track to the vertex
    E_Trk_Bcktrk_SurfPlanePath    = 3, // Crossing path from the target surface to the vertex plane
    E_Trk_Bcktrk_ExtrDist         = 4, // Distance in the plane from vertex to extrapolated point
    E_Trk_Bcktrk_dPtgt_spline     = 5, // DeltaP = P(target surface)-P(spline)
    E_Trk_Bcktrk_dPtgt_spline_cor = 6, // DeltaP = P(corrected for crossing)-P(target surface)
    E_Trk_Bcktrk_pid              = 7, // Particle type used for backtracking (geant convention)
    E_Trk_Bcktrk_End              = 8, // Used for arrays
  };

  // --- data-members
 protected:

  Int_t fVertexID; // ID-number of associated vertex (-1 if K_stop)

  // General features
  Float_t fCharge;
  Bool_t fChargeAmbig; // charge opposite to the expected one after fit

  Bool_t fFwd; // true for fwd tracks

  Float_t fStart[E_Trk_Vect_End]; // start point (see TFndTrack::TrkVect_e)
  Float_t fStop[E_Trk_Vect_End];  // stop point  (see TFndTrack::TrkVect_e)

  Float_t fTOF; // time of flight (used mainly for neutral tracks)
  Float_t fEnergy; // Energy (evaluated on TOFONE for neutral tracks)
  // ---
  
  // by pattern recognition
  Short_t fNumOfHits; // number of hits used for reconstruction
  UInt_t fPath; // sequence of hit detectors
  // --- description of "path" [hex code] (zero in case of missing hit):
  //      digit 0: TOFINO slab (1-12) [0x1  - 0xC]
  //      digit 1: ISIM module (1- 8) [0x1  - 0x8]
  //      digit 2: OSIM module (1-10) [0x11 - 0xA]
  //      digit 3: INN LMD dch (1- 8) [0x1  - 0x8]
  //      digit 4: OUT LMD dch (1- 8) [0x1  - 0x8]
  //      digits 5-6: TOFONE slab (1-72) [0x1  - 0x48]
  //      digit 7: free
  

  Float_t fPR_Helix[E_Hel_Geo_End]; // pattern-recognized helix
  Float_t fPR_MomVert[E_Hel_Mom_End]; // pattern-recognized momentum (used also for neutral tracks)

  // ---

  // by fitting process and extrapolation
  Int_t fFitErr; // (see TFndTrack::Trk_FitError_e)
  Float_t fLength; // track length (cm)
  Float_t fChiSq[E_Trk_ChiSq_End];

  Float_t fFit_MomVert[E_Hel_Mom_End]; // fitted momentum
  Float_t fFit_MomVertErr[E_Hel_Mom_End]; // errors on fitted momentum

  Float_t fDeDx[2]; // dE/dx [ISIM,OSIM]
  Float_t fCorr_DeDx[4]; //  dE/dx [ISIM,OSIM,LMD1,LMD2] corrected for attack-angle

  Bool_t fExtrapolated; // true if back-extrapoletion has been performed
  Float_t fBackTrack[E_Trk_Bcktrk_End]; // note: E_Trk_Bcktrk_pdg must be converted from "int"
        
  Int_t fType; // type of tracks (number of hits per layer)
  Float_t fDistVert; // distance from vertex

  // director cosines (track helix/detector)
  Float_t fIsimCos;
  Float_t fOsimCos;
  Float_t fLmd1Cos;
  Float_t fLmd2Cos;

  Float_t fSpecMom; // momentum extimated inside the spectrometer

  void EvalTrack();

 public:
  TFndTrack();
  ~TFndTrack();
  
  void Reset(); // reset all data-members
  void SetPath(UInt_t tofi,UInt_t isim,UInt_t osim,UInt_t lmd_inn,UInt_t lmd_out,UInt_t tofo);
  UInt_t GetPath(UInt_t &tofi,UInt_t &isim,UInt_t &osim,UInt_t &lmd_inn,UInt_t &lmd_out,UInt_t &tofo);

  // --- following methods can be
  //     used in order to get the
  //     single-detector-channel
  UInt_t GetPath_Tofino();
  UInt_t GetPath_Isim();
  UInt_t GetPath_Osim();
  UInt_t GetPath_Lmd_Inn();
  UInt_t GetPath_Lmd_Out();
  UInt_t GetPath_Tofone();

  void PrintData(Int_t depth = 1);

  friend class TFndPrean;
  friend class TFndPreanHistos;
  friend class TFndEvd;
  ClassDef(TFndTrack,4) // FINUDA Track general structure

    // streamer versions:
    //  1: valid up to fidarc-602 (long tracks only)
    //  2: used for test purposes (while adding short tracks support)
    //  3: used for test purposes (while adding short tracks support)
    //  4: valid from fidarc-603 (short tracks are used) / corrected dE/dx added
};

#endif // FROOT_TFndTrack

---