AlignmentStripHistograms.cxx

Go to the documentation of this file.

//
//
// AlignmentStripHistograms
//
// Package: DetectorAlignment
//
// Main alignment histograms class
//
// Contact: rustem@fnal.gov
//
// Created on: Tue Jul 26 11:40:05 2004
//

//Local package headers
#include "AlignmentStripHistograms.h"
#include "AlignmentStrip.h"

//MINOS headers
#include "MessageService/MsgService.h"
#include "MessageService/MsgFormat.h"

//ROOT headers
#include "TH1.h"
#include "TProfile.h"
#include "TDirectory.h"

//C++ headers
#include <cassert>
#include <cmath>
#include <iostream>
#include <string>
#include <sstream>

using namespace std;

CVSID("$Id: AlignmentStripHistograms.cxx,v 1.4 2009/02/28 21:46:10 gmieg Exp $");

//__________________________________________________________________________________________________
AlignmentStripHistograms::AlignmentStripHistograms(TDirectory *ResidualDir,
                                                   TDirectory *ChargeDir,
                                                   TDirectory *AttenDir,
                                                   TDirectory *OccupancyDir,
                                                   TDirectory *RotResDir,
                                                   const AlignmentStrip &astrip)
   :fPlane(astrip.plane),
    fStrip(astrip.strip),
    fPlexSEId(astrip.plexseid),
    fLength(astrip.length),
    fPigTail(astrip.pigtail),
    fhCharge(0),
    fhChargePathCorr(0), 
    fhResidual(0), 
    fhOccupancy(0),
    fpAtten(0),
    fpAttenPigCorr(0),
    fpRotationalResidual(0),
    fNChargeBins(200),
    fNResidualBins(60),
    fChargeLimit(5000.0),
    fResidualLimit(0.12),
    fOccupancyBinWidth(0.04),
    fAttenuationBinWidth(0.4),
    fRotationalResidualBinWidth(0.5)
{

   MSG("Align", Msg::kVerbose) << "Constructor AlignmentStripHistograms() " << endl;

   if(!ResidualDir || !ChargeDir || !AttenDir || !OccupancyDir || !RotResDir)
   {
      MSG("Align", Msg::kFatal) << "AlignmentStripHistograms ctor: TDirectory pointer is zero." << endl;
      abort();
   }

   fhCharge             = InitChargeHistogram(ChargeDir);
   fhChargePathCorr     = InitPathCorrChargeHistogram(ChargeDir);
   fhResidual           = InitResidualHistogram(ResidualDir);
   fhOccupancy          = InitOccupancyHistogram(OccupancyDir);
   fpAtten              = InitAttenHistogram(AttenDir);
   fpAttenPigCorr       = InitPigtailAttenHistogram(AttenDir);
   fpRotationalResidual = InitRotResidualHistogram(RotResDir);
}


//__________________________________________________________________________________________________
void AlignmentStripHistograms::FillHistograms(const AlignmentStrip& astrip)
{
   fhCharge         -> Fill(astrip.charge);
   fhChargePathCorr -> Fill(astrip.charge_pc);
   fhResidual       -> Fill(astrip.residual);
   fhOccupancy      -> Fill(astrip.lhitpos);
   
   if(!astrip.goodhit)
      return;
   
   fpRotationalResidual -> Fill(astrip.lhitpos+astrip.lposrelmdl, astrip.residual);

   if(astrip.wlsbypass > 0.001)
      return;

   fpAtten        -> Fill(fabs(astrip.lhitpos - astrip.length/2.0), astrip.charge_pc);
   fpAttenPigCorr -> Fill(astrip.pigtail + fabs(astrip.lhitpos - astrip.length/2.0), astrip.charge_pc);  
}

//__________________________________________________________________________________________________
TH1D* AlignmentStripHistograms::InitChargeHistogram(TDirectory *dir)
{
   char name[100], title[100];
   sprintf(name,"ChargeInPlane%03dStrip%03d", fPlane, fStrip);
   sprintf(title,"Charge plane %03d, strip %03d", fPlane, fStrip);
   MSG("Align", Msg::kVerbose) << "Creating histogram " << name << endl;
   TH1D *h = new TH1D(name, title, fNChargeBins, 0, fChargeLimit);
   h -> SetDirectory(dir);
   h -> GetXaxis() -> SetTitle("Raw ADC");
   h -> GetXaxis() -> CenterTitle();      
   return h;
}

//__________________________________________________________________________________________________
TH1D* AlignmentStripHistograms::InitPathCorrChargeHistogram(TDirectory *dir)
{
   char name[100], title[100];
   sprintf(name,"PathCorrChargeInPlane%03dStrip%03d", fPlane, fStrip);
   sprintf(title,"Path corrected charge plane %03d, strip %03d", fPlane, fStrip);
   MSG("Align", Msg::kVerbose) << "Creating histogram " << name << endl;
   TH1D *h = new TH1D(name, title, fNChargeBins, 0, fChargeLimit);
   h -> SetDirectory(dir);
   h -> GetXaxis() -> SetTitle("Raw ADC");
   h -> GetXaxis() -> CenterTitle();      
   return h;
}

//__________________________________________________________________________________________________
TH1D* AlignmentStripHistograms::InitResidualHistogram(TDirectory *dir)
{
   char name[100], title[100];
   sprintf(name,"ResidualInPlane%03dStrip%03d", fPlane, fStrip);
   sprintf(title,"Residual in hit strip, plane %03d, strip %03d", fPlane, fStrip);
   TH1D *h = new TH1D(name, title, 2*fNResidualBins, -fResidualLimit, fResidualLimit);
   h -> SetDirectory(dir);
   h -> GetXaxis() -> SetTitle("residual (m)");
   h -> GetXaxis() -> CenterTitle();
   return h;
}


//__________________________________________________________________________________________________
TH1D* AlignmentStripHistograms::InitOccupancyHistogram(TDirectory *dir)
{
   int nbin = int(fLength/fAttenuationBinWidth);
   char name[100], title[100];
   sprintf(name,"OccupancyInPlane%03dStrip%03d", fPlane, fStrip);
   sprintf(title,"Longitudinal occupancy, plane %03d, strip %03d", fPlane, fStrip);
   TH1D *h = new TH1D(name, title, nbin, -fLength/2.0, fLength/2.0);
   h -> SetDirectory(dir);
   h -> GetXaxis() -> SetTitle("residual (m)");
   h -> GetXaxis() -> CenterTitle();
   return h;
}

//__________________________________________________________________________________________________
TProfile* AlignmentStripHistograms::InitAttenHistogram(TDirectory *dir)
{       
   char name[100], title[100];
   int nbin = int(fLength/fAttenuationBinWidth);
   sprintf(name,"ProfileStripAtten%03d_%03d", fPlane, fStrip);
   sprintf(title,"Strip attenuation plane %03d and strip %03d", fPlane, fStrip);
   MSG("Align", Msg::kVerbose) << "Creating " << name << " histogram with "<< nbin << " bins"<<endl;
   TProfile *p = new TProfile(name, title, nbin, 0.0, fLength);
   p -> SetDirectory(dir);
   p -> GetXaxis() -> SetTitle("wls fiber length (m)");
   p -> GetXaxis() -> CenterTitle();
   p -> GetYaxis() -> SetTitle("charge (ADC)");
   p -> GetYaxis() -> CenterTitle();
   return p;
}   

//__________________________________________________________________________________________________
TProfile* AlignmentStripHistograms::InitPigtailAttenHistogram(TDirectory *dir)
{       
   char name[100], title[100];
   const int nbin = int(floor(fLength/fAttenuationBinWidth));
   sprintf(name,"ProfileStripPigAtten%03d_%03d", fPlane, fStrip);
   sprintf(title,"Strip attenuation plane %03d and strip %03d", fPlane, fStrip);
   MSG("Align", Msg::kVerbose) << "Creating " << name << " histogram with "<< nbin << " bins"<<endl;
   Double_t bins[10];
   bins[0] = 0.0;
   bins[1] = fPigTail;
   for(int i = 1; i <= 8; ++i)
   {
      double bin = i;
      bins[i+1] = bin*fAttenuationBinWidth + fPigTail;
   }
   TProfile *p = new TProfile(name, title, 9, bins);
   p -> SetDirectory(dir);
   p -> GetXaxis() -> SetTitle("wls fiber length (m)");
   p -> GetXaxis() -> CenterTitle();
   p -> GetYaxis() -> SetTitle("charge (ADC)");
   p -> GetYaxis() -> CenterTitle();
   return p;
}   

//__________________________________________________________________________________________________
TProfile* AlignmentStripHistograms::InitRotResidualHistogram(TDirectory *dir)
{
   char name[100], title[100];   
   sprintf(name,"RotResidualInPlane%03dStrip%03d", fPlane, fStrip);
   sprintf(title,"Profile rotational residuals plane %03d strip %03d", fPlane, fStrip);
   unsigned int nbin = int(floor(fLength/fRotationalResidualBinWidth));
   TProfile *p = new TProfile(name, title, nbin, -fLength/2.0, fLength/2.0);
   p -> SetDirectory(dir);
   p -> GetXaxis() -> SetTitle("position (m)");
   p -> GetXaxis() -> CenterTitle();
   p -> GetYaxis() -> SetTitle("residual (m)");
   p -> GetYaxis() -> CenterTitle();   
   return p;
}


//____________________________________________________________________________________________
AlignmentStripHistograms::~AlignmentStripHistograms()
{
   MSG("Align", Msg::kVerbose) << "Destructor ~AlignmentStripHistograms()..." << endl;
}

---