TrackDirectionModule.h

Go to the documentation of this file.

//
//
// TrackDirectionModule
//
// Package: DetectorAlignment
//
// Contact: rustem@fnal.gov
//
// Created on: Fri Jul 23 15:26:05 2005
//

#ifndef TRACKDIRECTIONMODULE_H
#define TRACKDIRECTIONMODULE_H

//LOCAL
#include "AlignmentStrip.h"

//MINOS
#include "JobControl/JobCModule.h"
#include "Plex/PlexStripEndId.h"

 //ROOT
#include "TStopwatch.h"

//C++
#include <vector>

using std::vector;

class VldContext;
class MomNavigator;
class NtpAlignmentRecord;
class UgliGeomHandle;
class UgliStripHandle;
class TFile;
class TTree;
class TH1D;
class TH2D;

class TrackDirectionModule : public JobCModule
{
public:
   
   TrackDirectionModule();
   virtual ~TrackDirectionModule(){};
   void BeginJob();              
   JobCResult Ana(const MomNavigator* mom);
   void EndJob();

private:
   
   bool ProcessRecord(const NtpAlignmentRecord* ntprec, 
                      const VldContext& vldc);
   
   bool GetTrackRecoQuality(const NtpAlignmentRecord* ntprec);
   
   void CountNHits(const PlexStripEndId& plexid);

   void ConvertToLocal(AlignmentStrip &astrip, 
                       PlexStripEndId &plexid, 
                       UgliGeomHandle &ugh,
                       UgliStripHandle &ush);

   void PrintJobStatistics();
   
   bool GetDirectionalCosines(Double_t &cosx, Double_t &cosy, Double_t &cosz, 
                              Double_t &length, Double_t &sigma, const VldContext &vld);
   
   TFile *fTrackDirectionFile;

   TStopwatch fTimer;
   TStopwatch fTimerInterval;
   
   int fRun;
   int fSubRun;
   
   vector<AlignmentStrip> fTrackVStrip;
   vector<AlignmentStrip> fTrackUStrip;

   TH1D *fhCosX;
   TH1D *fhCosY;
   TH1D *fhCosZ;
   TH1D *fhTrackLength;
   
   TH2D *f2hCosXvsCosZ;
   TH2D *f2hCosYvsCosZ;
   TH2D *f2hCosYvsCosX;

   //track quality cuts
   const Double_t fMaxTrackChargeCut;
   const Double_t fMinTrackChargeCut;
   const Double_t fMaxSigmaOfTPosCut;
   const Double_t fTrackChargeRatioCut;
   const Double_t fMinCosZCut;
   const Double_t fMinTrackStripChargeCut;
   const Double_t fMaxTrackStripChargeCut;
   
   bool fApplyCuts;  
   
   //tree and tree variables
   TTree *fTrackDirectionTree;
   Double_t fCosX;
   Double_t fCosY;
   Double_t fCosZ;
   Double_t fLength;
   Double_t fSigma;
   Double_t fTime;
   Double_t fAltitude;
   Double_t fAzimuth;
   Double_t fCharge;
   Double_t fNStrip;
   UInt_t fYear, fMonth, fDay, fHour, fMinute, fSec;

   unsigned int fNRecords;
   unsigned int fNFailedRead;
   unsigned int fFailedCut;
   unsigned int fNFailedCutMaxVStripCharge;
   unsigned int fNFailedCutMaxUStripCharge;
   unsigned int fNFailedTrackChargeCut;
   unsigned int fNFailedCosZCut;
   unsigned int fNFailedChargeRatioCut;
   unsigned int fNFailedSigmaTPosCut;
};

#endif //TRACKDIRECTIONMODULE_H

---