DTrackHitSelectorALT2.h

Go to the documentation of this file.

// $Id$
//
//    File: DTrackHitSelectorALT2.h
// Created: Fri Feb  6 08:22:58 EST 2009
// Creator: davidl (on Darwin harriet.jlab.org 9.6.0 i386)
//

#ifndef _DTrackHitSelectorALT2_
#define _DTrackHitSelectorALT2_
#include <TMath.h>
#include <TTree.h>
#include <JANA/jerror.h>
#include <DANA/DApplication.h>


#include <TRACKING/DTrackHitSelector.h>

class DTrackHitSelectorALT2:public DTrackHitSelector{
   public:
   DTrackHitSelectorALT2(jana::JEventLoop *loop, int32_t runnumber);
      virtual ~DTrackHitSelectorALT2();
      
      void GetCDCHits(fit_type_t fit_type, const DReferenceTrajectory *rt, const vector<const DCDCTrackHit*> &cdchits_in, vector<const DCDCTrackHit*> &cdchits_out,int N=20) const;
      void GetFDCHits(fit_type_t fit_type, const DReferenceTrajectory *rt, const vector<const DFDCPseudo*> &fdchits_in, vector<const DFDCPseudo*> &fdchits_out, int N=20) const;
      void GetCDCHits(double Bz,double q,const vector<DTrackFitter::Extrapolation_t> &extrapolations, const vector<const DCDCTrackHit*> &cdchits_in, vector<const DCDCTrackHit*> &cdchits_out, int N) const;
      void GetFDCHits(double Bz,double q,
            const vector<DTrackFitter::Extrapolation_t> &extrapolations, const vector<const DFDCPseudo*> &fdchits_in, vector<const DFDCPseudo*> &fdchits_out,int N) const;

   private:
      const DMagneticFieldMap *bfield;

      int HS_DEBUG_LEVEL;
      bool MAKE_DEBUG_TREES;
      double MIN_HIT_PROB_CDC;
      double MIN_HIT_PROB_FDC;
      double MIN_FDC_SIGMA_ANODE_CANDIDATE;
      double MIN_FDC_SIGMA_CATHODE_CANDIDATE;
      double MIN_FDC_SIGMA_ANODE_WIREBASED;
      double MIN_FDC_SIGMA_CATHODE_WIREBASED;
      double MAX_DOCA;
      
      TTree *cdchitsel;
      TTree *fdchitsel;

      typedef struct{
         int fit_type;
         float p;
         float theta;
         float mass;
         float sigma;
         float x;
         float y;
         float z;
         float s;
         float itheta02;
         float itheta02s;
         float itheta02s2;
         float dist;
         float doca;
         float resi;
         float chisq;
         float prob;
         float sig_phi;
         float sig_lambda;
         float sig_pt;
      }cdchitdbg_t;
      mutable cdchitdbg_t cdchitdbg;
      
      typedef struct{
         int fit_type;
         float p;
         float theta;
         float mass;
         float sigma_anode;
         float sigma_cathode;
         float x;
         float y;
         float z;
         float s;
         float itheta02;
         float itheta02s;
         float itheta02s2;
         float dist;
         float doca;
         float resi;
         float u;
         float u_cathodes;
         float resic;
         float chisq;
         float prob;
         float sig_phi;
         float sig_lambda;
         float sig_pt;
      }fdchitdbg_t;
      mutable fdchitdbg_t fdchitdbg;
};

#endif // _DTrackHitSelectorALT2_