gl3Track.h

//:>------------------------------------------------------------------
//: FILE:       gl3Tracks.h
//: HISTORY:
//:              6dec1999 version 1.00
//:              2feb2000 add sector to add methods
//:             27jul2000 add methods to drop hits
//:<------------------------------------------------------------------
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#ifndef TRG_VERSION
#define TRG_VERSION 0x32
#endif

#include "daqFormats.h"
#include "FtfTrack.h"
//#include "l3List.h"
#include "gl3Hit.h"

#ifndef GL3TRACK
#define GL3TRACK


class gl3Track: public FtfBaseTrack {
private:
   gl3Track* getNextTrack ( )    { return (gl3Track *)nextTrack ; } ;
public:
   void*     nextTrack ;
   int       sector ;
   float     dca;
   
   inline virtual   void nextHit ()
       { currentHit = ((gl3Hit *)currentHit)->nextHit ; } ;

   int addTrack ( gl3Track* ) ;
   void      dropHits ( int rest, int rowMin, int rowMax ) ;
   gl3Track* merge ( FtfContainer *trackArea ) ;
   void      Print ( int level ) ;

   void      setDca(Ftf3DHit vertex);
   
   //
   //#######################################################################
   float getRealEta ( ) {
      float theta = atan2(1.,(double)tanl);
      float rEta  = -1. * log (tan(theta/2.)) ;
      return rEta ;
   }
   //
   //#######################################################################
   void print ( ) {
      float gl3ToDeg = 180./acos(-1.);
      printf ( "pt %f tanl %f psi %f r0 %f z0 %f phi0 %f nHits %d\n", 
                pt, tanl, psi, r0, z0, phi0*gl3ToDeg, nHits ) ;
   }
   //
   //#######################################################################
   void set ( short sectorIn, local_track* trk ) {
       //    id          = sectorIn * 10000 + abs(trk->id) ;
       nHits       = trk->nHits ;
       nDedx       = trk->ndedx ;
       chi2[0]     = float(trk->xy_chisq)/10. ;
       chi2[1]     = float(trk->sz_chisq)/10. ;
       dedx        = trk->dedx ; 
       pt          = fabs(trk->pt) ;
       psi         = trk->psi ;
       tanl        = trk->tanl ;
       eta         = getRealEta();
       z0          = trk->z0 ;
       length      = trk->trackLength ;
       innerMostRow= trk->innerMostRow ;
       outerMostRow= trk->outerMostRow ;
       r0          = trk->r0 ;
       q           = (short )(trk->pt/fabs(trk->pt)) ;
       phi0        = trk->phi0 ;
       dpt         = float(trk->dpt)/32768. * pt ; 
       dpsi        = DecompressOver1(trk->dpsi,psi);
       dtanl       = DecompressOver1(trk->dtanl,tanl)/64.;
       dz0         = float(trk->dz0)/1024. ;
       nextTrack   = 0 ;
       firstHit    = 0 ;
       lastHit     = 0 ;

       dca = -1;
       //
       // Check errors are not zero
       //
       if ( dpt   == 0 ) dpt   = 1.e-5 * pt ;
       if ( dpsi  == 0 ) dpsi  = 1.e-5 ;
       if ( dtanl == 0 ) dtanl = 1.e-5 ;
       if ( dz0   == 0 ) dz0   = 1.e-3 ; 
       
   }
   //
   //#######################################################################
   void set ( FtfTrack* trk ) {
       id          = abs(trk->id) ;
       printf ( "set %d\n", id ) ;
       nHits       = trk->nHits;
       chi2[0]     = trk->chi2[0] ;
       chi2[1]     = trk->chi2[1] ; 
       dedx        = trk->dedx ; 
       pt          = fabs(trk->pt) ;
       psi         = trk->psi ;
       tanl        = trk->tanl ;
       z0          = trk->z0 ;
       length      = trk->length ;
       r0          = trk->r0 ;
       phi0        = trk->phi0 ;
       q           = trk->q ;
       dtanl = dpsi = dpt = 0. ;
              
       dca = -1;
   }
};

#endif