AlignmentAlgorithm.h

Go to the documentation of this file.

//
//
// AlignmentAlgorithm
//
// Package: DetectorAlignment
//
// An alignment engine's class 
//
// Contact: rustem@fnal.gov
//
// Created: Wed Jul 27 21:25:05 2005
//

#ifndef ALIGNMENTALGORITHM_H
#define ALIGNMENTALGORITHM_H

//Local package headers
#include "SubsetStat.h"
#include "AlignmentStripHistograms.h"
#include "StatisticalFluctuations.h"

//MINOS headers
#include "Plex/PlexStripEndId.h"

//C++ headers
#include <map>
#include <cstring>

using std::map;

class TDirectory;
class TH1D;
class TProfile;
class TFile;
class NtpAlignmentRecord;
class VldContext;

class AlignmentAlgorithm
{
   
public:
   
   AlignmentAlgorithm(TFile *rootfile, const VldContext vldc);

   virtual ~AlignmentAlgorithm();
   
   void Fill(const AlignmentStrip& astrip);

   void FillCandStrip(const AlignmentStrip& astrip);

   void ResetSubset(const VldContext &vldc, const int ntracks) {fFluctuations.Reset(vldc, ntracks);};

   void FitStripAttenuationHistograms();

   void RunAlgorithm();

private:

   TFile *fRootFile;
   TDirectory *fdNonTrackStripChargeDir;
   TDirectory *fdMdlResDir;
   TDirectory *fdMdlRotResDir;
   TDirectory *fdPlaneRotResDir;

   TDirectory *fdResidualDir;
   TDirectory *fdChargeDir;
   TDirectory *fdAttenDir;
   TDirectory *fdOccupancyDir;
   TDirectory *fdRotResDir;

   TProfile* GetMdlRotationalResidualHistogram(int plane, int module, double low, double high);
   TProfile* GetPlaneRotationalResidualHistogram(int plane, double low, double high);

   void AnalyzeStripData();

   void AnalyzeModuleResiduals();

   void FitPlaneRotRes();
   
   StatisticalFluctuations fFluctuations;

   //Module residuals
   map<PlexScintMdlId, TH1D *> fmModuleResidual; 

   //Object that and initilizes and holds pointers to strip histograms
   map<PlexStripEndId, AlignmentStripHistograms> fmStripHistograms; 
   
   //Rotational module residual TProfile histograms
   map<PlexScintMdlId, TProfile *> fpMdlRotationalResidual; 

   //Rotational plane residual TProfile histograms
   map<PlexPlaneId, TProfile *> fpPlaneRotationalResidual; 

   //Candidate strip charge histograms
   map<PlexStripEndId, TH1D *> fhCandStripCharge;

   bool fAlignmentAlgorithmCompleted;

   const Int_t fNStripChargeBins;
   const Int_t fNMdlResBins;
   const Double_t fStripChargeLimit;
   const Double_t fMdlResLimit;
   const Double_t fMdlRotResBinWidth;
   const Double_t fMinResHistEntriesForFit;
   const Double_t fMinAttenHistEntriesForFit;
   const Double_t fMinRotResHistEntriesForFit;
   const Double_t fMinStripRotResHistEntriesForFit;
   const Double_t fPlaneRotResBinWidth;
   
   //initial values for strip attenuation exponential fit
   const Double_t fStripAttenPar0;
   const Double_t fStripAttenPar1;
};

#endif //ALIGNMENTHISTOGRAMS_H

---