AlignmentHistograms.h

Go to the documentation of this file.

//
//
// AlignmentHistograms
//
// Package: DetectorAlignment
//
// An alignment engine's histograming class 
//
// Contact: rustem@fnal.gov
//
// Created: Tuesday Sep 21 20:40:05 2004
//

#ifndef ALIGNMENTHISTOGRAMS_H
#define ALIGNMENTHISTOGRAMS_H

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

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

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

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

class TDirectory;
class TH1;
class TH1D;
class TH2D;
class TProfile;
class TTree;
class TFile;
class VldContext;
class NtpAlignmentRecord;

class AlignmentHistograms
{
   
public:
   
   AlignmentHistograms(TFile *rootfile, const VldContext &vld);

   virtual ~AlignmentHistograms();
   
   void Fill(const AlignmentStrip& astrip, const NtpAlignmentRecord* ntprec);
   
   void FillPlanePlex(const vector<AlignmentStrip> &track_strips,
                      const vector<AlignmentStrip> &cand_strips);
   
private:

   void FillTrackInfo(const NtpAlignmentRecord* ntprec);

   void FillHits(PlexStripEndId &plexid);

   void SetAxisLabels(TH1 *h, const char *xlabel=0, const char *ylabel=0);

   //map key is plane number
   map<int, TH2D *> f2hPlanePlexChargeRatio;
   map<int, TH2D *> f2hPlanePlexHits;

   //number of hits per module and strip
   map<PlexStripEndId, unsigned int> fNStripHits;
   map<PlexScintMdlId, unsigned int> fNMdlHits;

   TFile *fRootFile;
   TDirectory *fPlanePlexDirectory;
   TDirectory *fTrackDirectory;
   TDirectory *fResidualSummaryDirectory;

   int fnplanes;
   int fnstrips;
   const int fNResidualBins;
   const Double_t fResidualLimit;
   const Double_t fResidualStripBinWidth;
   const Double_t fResidualMdlBinWidth;
   const int fNResidualBinsStrip;
   const Double_t fResidualLimitStrip;
   const int fNChargeBins;
   const Double_t fChargeLimit;
   const int fNTrackChargeBins;
   const Double_t fTrackChargeLimit;
   const int fNHitsBins;
   const int fTrackResidualRMSBins;
   const Double_t fTrackResidualRMSLimit;
   const int fNCosBins;

   TH1D *fhCharge2dTracksU;
   TH1D *fhCharge2dTracksV;
   TH1D *fhCharge2dTracksDiffNorm;
   TH1D *fhNStrip2dTracksU;
   TH1D *fhNStrip2dTracksV;
   TH1D *fhNStrip2dTracksDiffNorm;
   TH1D *fhSigmaOfTPosVview;
   TH1D *fhSigmaOfTPosUview;
   TH1D *fhSigmaOfTPosDiffNorm;
   TH1D *fChargesInAllHitStrip;        
   TH1D *fChargevsResidualHitStrips;  
   TH1D *fResidualAllStrips; 
   TH1D *fResidualVStrips; 
   TH1D *fResidualUStrips; 
   TH1D *fTrackResidualRMSVview;
   TH1D *fTrackResidualRMSUview;
   TH1D *fChargeRatioOf2dTracks;
   TH1D *fhTrackHoughCosU;
   TH1D *fhTrackHoughCosV;
   TH1D *fhTrackHoughCosZ;
   TH1D *fhPlexPlaneEntries;
   TH1D *fhTrackChargeOverTotalChargeRatio;

   TH2D *f2hNumberOfTrackHitsPerStrip;        
   TH2D *f2hNumberOfTrackHitsPerModule;
   TH2D *f2hChargevsResidualHitStrips;  
   TH2D *f2hResidualvsPlane;
   TH2D *f2hMeanChargePlaneStrip;
   TH2D *f2hMeanPathCorrChargePlaneStrip;   

   //Track TProfile histograms
   TProfile *fpTPosSigmavsCosZVview;
   TProfile *fpTPosSigmavsCosZUview;
   TProfile *fpTPosSigmavsTrackChargeVview;
   TProfile *fpTPosSigmavsTrackChargeUview;
   TProfile *fpTPosSigmavsTrackNStripsVview;
   TProfile *fpTPosSigmavsTrackNStripsUview;
   TProfile *fpTPosSigmavsTrackChargeRatioVview;
   TProfile *fpTPosSigmavsTrackChargeRatioUview;
   TProfile *fp3dTrackChargevsTrackChargeRatio;
   TProfile *fp3dTrackChargevsTrackCosz;
   TProfile *fp3dTrackNStripsvsCosz;
   TProfile *fp3dTrackChargevsNStrips;
   TProfile *fp3dTrackChargeRatiovsCosz;

   //Strip TProfile histograms
   TProfile *fpResidualvsPlane;
   TProfile *fpResidualvsCosZVview;
   TProfile *fpResidualvsCosZUview;
   TProfile *fpResidualvsTPosSigmaVview;
   TProfile *fpResidualvsTPosSigmaUview;
   TProfile *fpResidualvsStripChargeVview;
   TProfile *fpResidualvsStripChargeUview;
   TProfile *fpResidualvsTrackChargeVview;
   TProfile *fpResidualvsTrackChargeUview;
   TProfile *fpResidualvsTrackNStripsVview;
   TProfile *fpResidualvsTrackNStripsUview;
   TProfile *fpResidualvsTrackChargeRatioVview;
   TProfile *fpResidualvsTrackChargeRatioUview;
   TProfile *fpMeanStripResponsevsResidual;
};

#endif //ALIGNMENTHISTOGRAMS_H

---