DEventProcessor_track_hists.h

Go to the documentation of this file.

// $Id$
//
//    File: DEventProcessor_track_hists.h
// Created: Wed Oct 10 13:30:37 EDT 2007
// Creator: davidl (on Darwin fwing-dhcp95.jlab.org 8.10.1 i386)
//

#ifndef _DEventProcessor_track_hists_
#define _DEventProcessor_track_hists_

#include <pthread.h>
#include <map>
using std::map;

#include <TTree.h>
#include <TFile.h>
#include <TH1.h>
#include <TH2.h>

#include <JANA/JFactory.h>
#include <JANA/JEventProcessor.h>
#include <JANA/JEventLoop.h>

#include <PID/DKinematicData.h>
#include <TRACKING/DReferenceTrajectory.h>
#include <TRACKING/DMCTrackHit.h>
#include <TRACKING/DMCThrown.h>
#include <CDC/DCDCTrackHit.h>
#include <FDC/DFDCHit.h>
#include <FDC/DFDCWire.h>
#include <HDGEOMETRY/DLorentzDeflections.h>

#include "dchit.h"
#include "trackpar.h"

class DEventProcessor_track_hists:public JEventProcessor{

   public:
      DEventProcessor_track_hists();
      ~DEventProcessor_track_hists();

      TTree *fdchits, *cdchits, *ttrack;
      dchit cdchit, fdchit;
      dchit *cdchit_ptr, *fdchit_ptr;
      trackpar trk;
      trackpar *trk_ptr;
      
      DCoordinateSystem target;
      const DLorentzDeflections *lorentz_def;//< Correction to FDC cathodes due to Lorentz force

      class hit_info_t{
         public:
            // Inputs
            DReferenceTrajectory *rt;
            const DCoordinateSystem *wire;
            double tdrift;
            
            // Outputs
            double doca;
            double dist;
            double tof;
            double u;
            double u_lorentz;
            int LRfit;
            bool LRis_correct;
            TVector3 pos_doca;
            TVector3 mom_doca;
            TVector3 pos_wire;

            void FindLR(vector<const DMCTrackHit*> &mctrackhits, const DLorentzDeflections *lorentz_def=NULL);
      };
      
      
   private:
      jerror_t init(void); ///< Invoked via DEventProcessor virtual method
      jerror_t brun(JEventLoop *loop, int32_t runnumber);
      jerror_t evnt(JEventLoop *loop, uint64_t eventnumber);   ///< Invoked via DEventProcessor virtual method
      jerror_t erun(void);             ///< Invoked via DEventProcessor virtual method
      jerror_t fini(void);             ///< Invoked via DEventProcessor virtual method
         
      pthread_mutex_t mutex;
      
      int NLRbad, NLRgood, NLRfit_unknown;
      int Nevents;
};

#endif // _DEventProcessor_track_hists_