MakeAlignmentModule.h

Go to the documentation of this file.

//
//
// MakeAlignmentModule
//
// Package: DetectorAlignment
//
// Alignment engine that processes ntuples files
//
// Contact: rustem@fnal.gov
//
// Created on: Fri Sep 24 15:26:05 2004
//

#ifndef MAKEALIGNMENTMODULE_H
#define MAKEALIGNMENTMODULE_H

//LOCAL
#include "AlignmentStrip.h"
#include "AlignmentAlgorithm.h"
#include "AlignmentHistograms.h"
#include "AlignmentRunSummary.h"

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

 //ROOT
#include "TStopwatch.h"

//C++
#include <iostream>
#include <string>
#include <map>
#include <vector>

using std::map;
using std::vector;

class MomNavigator;
class UgliGeomHandle;
class UgliStripHandle;

class AlignmentHistograms;
class NtpAlignmentRecord;
class NtpAlignTrackStrip;
class NtpAlignCandStrip;

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

private:
   
   bool ReadRecord(const NtpAlignmentRecord* ntprec);

   void ProcessRecord(const NtpAlignmentRecord* ntprec);

   bool RecalculateResiduals();
   
   bool GetTrackQuality(const NtpAlignmentRecord* ntprec);
   
   double FitTrackLessOne(vector<AlignmentStrip>::const_iterator begin, 
                          vector<AlignmentStrip>::const_iterator end, 
                          vector<AlignmentStrip>::const_iterator skip);
   
   void ConvertToLocal(AlignmentStrip &astrip, 
                       PlexStripEndId &plexid, 
                       UgliGeomHandle &ugh,
                       UgliStripHandle &ush);

   void PrintJobStatistics();

   void PrintAlignmentStrip(const AlignmentStrip& astrip);

   //three major alignment objects
   AlignmentAlgorithm* fAlgorithm;
   AlignmentHistograms* fHistogram;
   AlignmentRunSummary fRunSummary;
   
   const VldTimeStamp fNearCoilOn; //Time when ND coil was first turned on
   
   TFile *fRootFile;

   int fIteration; //Alignment iteration number
   VldContext fValidityContext; //VldContext of current record

   int fCurrentRun;
   int fCurrentSubRun;
   int fNProcessedSubruns;
   const int fNSubRunsPerSubset;

   vector<AlignmentStrip> fTrackVStrip;
   vector<AlignmentStrip> fTrackUStrip;
   vector<AlignmentStrip> fCandVStrip;
   vector<AlignmentStrip> fCandUStrip;

   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;

   const Int_t fSpecialNDRunNumber;
   
   bool fRecalculateResidual;
   bool fApplyCuts;

   const unsigned int fStatisticalSubsetSize;   

   unsigned int fNFailedFit;
   unsigned int fNRecords;
   unsigned int fNPassedRecords;
   unsigned int fNSubsetRecords;
   unsigned int fNFailedRead;
   unsigned int fNMissedVPlane;
   unsigned int fNMissedUPlane;
   unsigned int fTotalVPlaneHits;
   unsigned int fTotalUPlaneHits;
   unsigned int fFailedCut;
   unsigned int fNHitsOutsideStrip;
   unsigned int fNFailedCutMaxVStripCharge;
   unsigned int fNFailedCutMaxUStripCharge;
   
   Double_t fMissedVPlaneMeanCharge;
   Double_t fMissedUPlaneMeanCharge;
   Double_t fMissedVPlaneChargeSigma;
   Double_t fMissedUPlaneChargeSigma;
   
   unsigned int fNFailedTrackChargeCut;
   unsigned int fNFailedStripChargeCut;
   unsigned int fNFailedCosZCut;
   unsigned int fNFailedChargeRatioCut;
   unsigned int fNFailedSigmaTPosCut;

   TStopwatch fTimer;
   TStopwatch fTimerInterval;
};


#endif

---