PixelResolutionHistograms Class Reference

#include <PixelResolutionHistograms.h>

Public Member Functions
int	Fill (double dx, double dy, double cotalpha, double cotbeta, int qbin, int nxpix, int nypix)
const SimpleHistogramGenerator *	getGeneratorX (double cotalpha, double cotbeta, int qbin, bool singlex)
const SimpleHistogramGenerator *	getGeneratorY (double cotalpha, double cotbeta, int qbin, bool singley)
	PixelResolutionHistograms (std::string filename, std::string rootdir, std::string descTitle, unsigned int detType, double cotbetaBinWidth, double cotbetaLowEdge, int cotbetaBins, double cotalphaBinWidth, double cotalphaLowEdge, int cotalphaBins)
	PixelResolutionHistograms (std::string filename, std::string rootdir="", int detType=-1, bool ignore_multi=false, bool ignore_single=false, bool ignore_qBin=false)
int	status ()
virtual	~PixelResolutionHistograms ()

Private Attributes
TH2F *	binningHisto_
TAxis *	cotalphaAxis_
int	cotalphaBins_
double	cotalphaBinWidth_
double	cotalphaLowEdge_
TAxis *	cotbetaAxis_
int	cotbetaBins_
double	cotbetaBinWidth_
double	cotbetaLowEdge_
unsigned int	detType_
std::unique_ptr< TFile >	file_
SimpleHistogramGenerator *	qbinGen_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX]
TH1F *	qbinHist_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX]
int	qbins_
int	qbinWidth_
SimpleHistogramGenerator *	resMultiPixelXGen_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX][QBIN_HIST_MAX]
TH1F *	resMultiPixelXHist_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX][QBIN_HIST_MAX]
SimpleHistogramGenerator *	resMultiPixelYGen_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX][QBIN_HIST_MAX]
TH1F *	resMultiPixelYHist_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX][QBIN_HIST_MAX]
SimpleHistogramGenerator *	resSinglePixelXGen_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX]
TH1F *	resSinglePixelXHist_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX]
SimpleHistogramGenerator *	resSinglePixelYGen_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX]
TH1F *	resSinglePixelYHist_ [COTBETA_HIST_MAX][COTALPHA_HIST_MAX]
int	status_
bool	weOwnHistograms_

Detailed Description

Definition at line 22 of file PixelResolutionHistograms.h.

Constructor & Destructor Documentation

PixelResolutionHistograms() [1/2]

PixelResolutionHistograms::PixelResolutionHistograms	(	std::string	filename,
		std::string	rootdir,
		std::string	descTitle,
		unsigned int	detType,
		double	cotbetaBinWidth,
		double	cotbetaLowEdge,
		int	cotbetaBins,
		double	cotalphaBinWidth,
		double	cotalphaLowEdge,
		int	cotalphaBins
	)

Definition at line 62 of file PixelResolutionHistograms.cc.

References binningHisto_, cotalphaBins_, cotalphaBinWidth_, cotalphaLowEdge_, cotbetaBins_, cotbetaBinWidth_, cotbetaLowEdge_, detType_, file_, corrVsCorr::filename, timingPdfMaker::histo, cuy::ii, findQualityFiles::jj, GetRecoTauVFromDQM_MC_cff::kk, qbinHist_, qbins_, resMultiPixelXHist_, resMultiPixelYHist_, resSinglePixelXHist_, resSinglePixelYHist_, and runGCPTkAlMap::title.

    : weOwnHistograms_(true),                                           // we'll be making some histos
      detType_(detType),
      cotbetaBinWidth_(cotbetaBinWidth),
      cotbetaLowEdge_(cotbetaLowEdge),
      cotbetaBins_(cotbetaBins),
      cotalphaBinWidth_(cotalphaBinWidth),
      cotalphaLowEdge_(cotalphaLowEdge),
      cotalphaBins_(cotalphaBins),
      qbinWidth_(1),
      qbins_(4),
      binningHisto_(nullptr),
      resMultiPixelXHist_(),
      resSinglePixelXHist_(),  // all to nullptr
      resMultiPixelYHist_(),
      resSinglePixelYHist_(),  // all to nullptr
      qbinHist_(),             // all to nullptr
      file_(nullptr),
      status_(0),
      resMultiPixelXGen_(),
      resSinglePixelXGen_(),
      resMultiPixelYGen_(),
      resSinglePixelYGen_(),
      qbinGen_() {
  file_ = std::make_unique<TFile>(filename.c_str(), "RECREATE");
  //Resolution binning
  // const double   cotbetaBinWidth = 1.0;
  // const double   cotbetaLowEdge  = -11.5 ;
  // const int  cotbetaBins = 23;
  // const double   cotalphaBinWidth    = 0.08 ;
  // const double   cotalphaLowEdge = -0.36 ;
  // const int  cotalphaBins    = 9;
  // const int  qbinWidth   = 1;
  // const int  qbins       = 4;

  // Dummy 2D histogram to store binning:
  binningHisto_ = new TH2F("ResHistoBinning",
                           descTitle.c_str(),
                           cotbetaBins,
                           cotbetaLowEdge,
                           cotbetaLowEdge + cotbetaBins * cotbetaBinWidth,
                           cotalphaBins,
                           cotalphaLowEdge,
                           cotalphaLowEdge + cotalphaBins * cotalphaBinWidth);

  // Store detType in the underflow bin
  binningHisto_->SetBinContent(0, 0, detType_);

  // All other histograms:
  Char_t histo[200];
  Char_t title[200];
  //
  //--- Histograms for clusters with multiple pixels hit in a given direction.
  //
  for (int ii = 0; ii < cotbetaBins_; ii++) {
    for (int jj = 0; jj < cotalphaBins_; jj++) {
      for (int kk = 0; kk < qbins_; kk++) {
        //
        sprintf(histo, "hx%d1%02d%d%d", detType_, ii + 1, jj + 1, kk + 1);  //information of bits of histogram names
        //--- First bit 1/0 barrel/forward, second 1/0 multi/single, cotbeta, cotalpha, qbins
        sprintf(title,
                "cotbeta %.1f-%.1f cotalpha %.2f-%.2f qbin %d npixel>1 X",
                cotbetaLowEdge_ + ii * cotbetaBinWidth_,
                cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
                cotalphaLowEdge_ + jj * cotalphaBinWidth_,
                cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_,
                kk + 1);
        //
        resMultiPixelXHist_[ii][jj][kk] = new TH1F(histo, title, 1000, -0.05, 0.05);

        sprintf(histo, "hy%d1%02d%d%d", detType_, ii + 1, jj + 1, kk + 1);
        sprintf(title,
                "cotbeta %.1f-%.1f cotalpha %.2f-%.2f qbin %d npixel>1 Y",
                cotbetaLowEdge_ + ii * cotbetaBinWidth_,
                cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
                cotalphaLowEdge_ + jj * cotalphaBinWidth_,
                cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_,
                kk + 1);
        //
        resMultiPixelYHist_[ii][jj][kk] = new TH1F(histo, title, 1000, -0.05, 0.05);
      }
    }
  }

  //
  //--- Histograms for clusters where only a single pixel was hit in a given direction.
  //
  for (int ii = 0; ii < cotbetaBins_; ii++) {
    for (int jj = 0; jj < cotalphaBins_; jj++) {
      sprintf(histo, "hx%d0%02d%d", detType_, ii + 1, jj + 1);  //information of bits of histogram names
      //first bit 1/0 barrel/forward, second 1/0 multi/single, cotbeta, cotalpha
      sprintf(title,
              "cotbeta %.1f-%.1f cotalpha %.2f-%.2f npixel=1 X",
              cotbetaLowEdge_ + ii * cotbetaBinWidth_,
              cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
              cotalphaLowEdge_ + jj * cotalphaBinWidth_,
              cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_);
      //
      resSinglePixelXHist_[ii][jj] = new TH1F(histo, title, 1000, -0.05, 0.05);

      sprintf(histo, "hy%d0%02d%d", detType_, ii + 1, jj + 1);
      sprintf(title,
              "cotbeta %.1f-%.1f cotalpha %.2f-%.2f npixel=1 Y",
              cotbetaLowEdge_ + ii * cotbetaBinWidth_,
              cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
              cotalphaLowEdge_ + jj * cotalphaBinWidth_,
              cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_);
      //
      resSinglePixelYHist_[ii][jj] = new TH1F(histo, title, 1000, -0.05, 0.05);

      sprintf(histo, "hqbin%d%02d%d", detType_, ii + 1, jj + 1);
      sprintf(title,
              "cotbeta %.1f-%.1f cotalpha %.2f-%.2f qbin",
              cotbetaLowEdge_ + ii * cotbetaBinWidth_,
              cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
              cotalphaLowEdge_ + jj * cotalphaBinWidth_,
              cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_);
      //
      qbinHist_[ii][jj] = new TH1F(histo, title, 4, -0.49, 3.51);
    }
  }
}

PixelResolutionHistograms() [2/2]

PixelResolutionHistograms::PixelResolutionHistograms	(	std::string	filename,
		std::string	rootdir = "",
		int	detType = -1,
		bool	ignore_multi = false,
		bool	ignore_single = false,
		bool	ignore_qBin = false
	)

Definition at line 204 of file PixelResolutionHistograms.cc.

References binningHisto_, cotalphaAxis_, cotalphaBins_, cotalphaBinWidth_, cotalphaLowEdge_, cotbetaAxis_, cotbetaBins_, cotbetaBinWidth_, cotbetaLowEdge_, detType_, file_, corrVsCorr::filename, timingPdfMaker::histo, cuy::ii, findQualityFiles::jj, GetRecoTauVFromDQM_MC_cff::kk, LOGDEBUG, LOGERROR, LOGINFO, qbinGen_, qbinHist_, qbins_, resMultiPixelXGen_, resMultiPixelXHist_, resMultiPixelYGen_, resMultiPixelYHist_, resSinglePixelXGen_, resSinglePixelXHist_, resSinglePixelYGen_, resSinglePixelYHist_, status_, AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::title.

    : weOwnHistograms_(false),  // resolution histograms are owned by the ROOT file
      detType_(-1),
      cotbetaBinWidth_(0),
      cotbetaLowEdge_(0),
      cotbetaBins_(0),
      cotalphaBinWidth_(0),
      cotalphaLowEdge_(0),
      cotalphaBins_(0),
      qbinWidth_(1),
      qbins_(4),
      binningHisto_(nullptr),
      resMultiPixelXHist_(),
      resSinglePixelXHist_(),  // all to nullptr
      resMultiPixelYHist_(),
      resSinglePixelYHist_(),  // all to nullptr
      qbinHist_(),             // all to nullptr
      file_(nullptr),
      status_(0),
      resMultiPixelXGen_(),
      resSinglePixelXGen_(),
      resMultiPixelYGen_(),
      resSinglePixelYGen_(),
      qbinGen_() {
  Char_t histo[200];        // the name of the histogram
  Char_t title[200];        // histo title, for debugging and sanity checking (compare inside file)
  TH1F* tmphist = nullptr;  // cache for histo pointer

  //--- Open the file for reading.
  file_ = std::make_unique<TFile>(filename.c_str(), "READ");
  if (!file_) {
    status_ = 1;
    LOGERROR << "PixelResolutionHistograms:: Error, file " << filename << " not found.";
    return;  // PixelTemplateSmearerBase will throw an exception upon our return.
  }

  //--- The dummy 2D histogram with the binning of cot\beta and cot\alpha:
  binningHisto_ = (TH2F*)file_->Get(Form("%s%s", rootdir.c_str(), "ResHistoBinning"));
  if (!binningHisto_) {
    status_ = 11;
    LOGERROR << "PixelResolutionHistograms:: Error, binning histogrram ResHistoBinning not found.";
    return;  // PixelTemplateSmearerBase will throw an exception upon our return.
  }

  if (detType == -1) {
    //--- Fish out detType from the underflow bin:
    detType_ = binningHisto_->GetBinContent(0, 0);
  } else {
    detType_ = detType;  // constructor's argument overrides what's in ResHistoBinning histogram.
  }

  //--- Now we fill the binning variables:
  cotbetaAxis_ = binningHisto_->GetXaxis();
  cotbetaBinWidth_ = binningHisto_->GetXaxis()->GetBinWidth(1);  // assume all same width
  cotbetaLowEdge_ = binningHisto_->GetXaxis()->GetXmin();        // low edge of the first bin
  cotbetaBins_ = binningHisto_->GetXaxis()->GetNbins();
  cotalphaAxis_ = binningHisto_->GetYaxis();
  cotalphaBinWidth_ = binningHisto_->GetYaxis()->GetBinWidth(1);  // assume all same width;
  cotalphaLowEdge_ = binningHisto_->GetYaxis()->GetXmin();        // low edge of the first bin;
  cotalphaBins_ = binningHisto_->GetYaxis()->GetNbins();

  if (!ignore_multi) {
    //
    //--- Histograms for clusters with multiple pixels hit in a given direction.
    //
    for (int ii = 0; ii < cotbetaBins_; ii++) {
      for (int jj = 0; jj < cotalphaBins_; jj++) {
        for (int kk = 0; kk < qbins_; kk++) {
          //
          sprintf(histo, "hx%d1%02d%d%d", detType_, ii + 1, jj + 1, kk + 1);  //information of bits of histogram names
          //--- First bit 1/0 barrel/forward, second 1/0 multi/single, cotbeta, cotalpha, qbins
          sprintf(title,
                  "cotbeta %.1f-%.1f cotalpha %.2f-%.2f qbin %d npixel>1 X",
                  cotbetaLowEdge_ + ii * cotbetaBinWidth_,
                  cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
                  cotalphaLowEdge_ + jj * cotalphaBinWidth_,
                  cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_,
                  kk + 1);
          //
          tmphist = (TH1F*)file_->Get(Form("%s%s", rootdir.c_str(), histo));
          if (!tmphist) {
            status_ = 2;
            LOGERROR << "Failed to find histogram=" << std::string(histo);
            return;
          }
          LOGDEBUG << "Found histo " << std::string(histo) << " with title = " << std::string(tmphist->GetTitle())
                   << std::endl;
          if (tmphist->GetEntries() < 5) {
            LOGINFO << "Histogram " << std::string(histo) << " has only " << tmphist->GetEntries()
                    << " entries. Trouble ahead." << std::endl;
          }
          resMultiPixelXHist_[ii][jj][kk] = tmphist;
          resMultiPixelXGen_[ii][jj][kk] = new SimpleHistogramGenerator(tmphist);

          sprintf(histo, "hy%d1%02d%d%d", detType_, ii + 1, jj + 1, kk + 1);
          sprintf(title,
                  "cotbeta %.1f-%.1f cotalpha %.2f-%.2f qbin %d npixel>1 Y",
                  cotbetaLowEdge_ + ii * cotbetaBinWidth_,
                  cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
                  cotalphaLowEdge_ + jj * cotalphaBinWidth_,
                  cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_,
                  kk + 1);
          //
          tmphist = (TH1F*)file_->Get(Form("%s%s", rootdir.c_str(), histo));
          if (!tmphist) {
            status_ = 3;
            LOGERROR << "Failed to find histogram=" << std::string(histo);
            return;
          }
          LOGDEBUG << "Found histo " << std::string(histo) << " with title = " << std::string(tmphist->GetTitle())
                   << std::endl;
          if (tmphist->GetEntries() < 5) {
            LOGINFO << "Histogram " << std::string(histo) << " has only " << tmphist->GetEntries()
                    << " entries. Trouble ahead." << std::endl;
          }
          resMultiPixelYHist_[ii][jj][kk] = tmphist;
          resMultiPixelYGen_[ii][jj][kk] = new SimpleHistogramGenerator(tmphist);
        }
      }
    }
    //
  }  // if (not ignore multi)

  //
  //--- Histograms for clusters where only a single pixel was hit in a given direction.
  //
  for (int ii = 0; ii < cotbetaBins_; ii++) {
    for (int jj = 0; jj < cotalphaBins_; jj++) {
      //--- Single pixel, along X.
      //
      sprintf(histo, "hx%d0%02d%d", detType_, ii + 1, jj + 1);  //information of bits of histogram names
      //--- First bit 1/0 barrel/forward, second 1/0 multi/single, cotbeta, cotalpha
      sprintf(title,
              "cotbeta %.1f-%.1f cotalpha %.2f-%.2f npixel=1 X",
              cotbetaLowEdge_ + ii * cotbetaBinWidth_,
              cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
              cotalphaLowEdge_ + jj * cotalphaBinWidth_,
              cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_);
      //
      tmphist = (TH1F*)file_->Get(Form("%s%s", rootdir.c_str(), histo));
      if (!tmphist) {
        if (!ignore_single) {
          LOGERROR << "Failed to find histogram=" << std::string(histo);
          status_ = 4;
          return;
        }
      } else {
        LOGDEBUG << "Found histo " << std::string(histo) << " with title = " << std::string(tmphist->GetTitle())
                 << std::endl;
        LOGDEBUG << "Found histo with title = " << std::string(tmphist->GetTitle()) << std::endl;
        if (tmphist->GetEntries() < 5) {
          LOGINFO << "Histogram " << std::string(histo) << " has only " << tmphist->GetEntries()
                  << " entries. Trouble ahead." << std::endl;
        }
        resSinglePixelXHist_[ii][jj] = tmphist;
        resSinglePixelXGen_[ii][jj] = new SimpleHistogramGenerator(tmphist);
      }

      //--- Single pixel, along Y.
      //
      sprintf(histo, "hy%d0%02d%d", detType_, ii + 1, jj + 1);
      sprintf(title,
              "cotbeta %.1f-%.1f cotalpha %.2f-%.2f npixel=1 Y",
              cotbetaLowEdge_ + ii * cotbetaBinWidth_,
              cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
              cotalphaLowEdge_ + jj * cotalphaBinWidth_,
              cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_);
      //
      tmphist = (TH1F*)file_->Get(Form("%s%s", rootdir.c_str(), histo));
      if (!tmphist) {
        if (!ignore_single) {
          LOGERROR << "Failed to find histogram=" << std::string(histo);
          status_ = 5;
          return;
        }
      } else {
        LOGDEBUG << "Found histo " << std::string(histo) << " with title = " << std::string(tmphist->GetTitle())
                 << std::endl;
        if (tmphist->GetEntries() < 5) {
          LOGINFO << "Histogram " << std::string(histo) << " has only " << tmphist->GetEntries()
                  << " entries. Trouble ahead." << std::endl;
        }
        resSinglePixelYHist_[ii][jj] = tmphist;
        resSinglePixelYGen_[ii][jj] = new SimpleHistogramGenerator(tmphist);
      }

      //--- qBin distribution, for this (cotbeta, cotalpha) bin.
      //
      sprintf(histo, "hqbin%d%02d%d", detType_, ii + 1, jj + 1);
      sprintf(title,
              "cotbeta %.1f-%.1f cotalpha %.2f-%.2f qbin",
              cotbetaLowEdge_ + ii * cotbetaBinWidth_,
              cotbetaLowEdge_ + (ii + 1) * cotbetaBinWidth_,
              cotalphaLowEdge_ + jj * cotalphaBinWidth_,
              cotalphaLowEdge_ + (jj + 1) * cotalphaBinWidth_);
      //
      tmphist = (TH1F*)file_->Get(Form("%s%s", rootdir.c_str(), histo));
      if (!tmphist) {
        if (!ignore_qBin) {
          LOGERROR << "Failed to find histogram=" << std::string(histo);
          status_ = 6;
          return;
        }
      } else {
        LOGDEBUG << "Found histo " << std::string(histo) << " with title = " << std::string(tmphist->GetTitle())
                 << std::endl;
        if (tmphist->GetEntries() < 5) {
          LOGINFO << "Histogram " << std::string(histo) << " has only " << tmphist->GetEntries()
                  << " entries. Trouble ahead." << std::endl;
        }
        qbinHist_[ii][jj] = tmphist;
        qbinGen_[ii][jj] = new SimpleHistogramGenerator(tmphist);
      }
    }
  }
}

~PixelResolutionHistograms()

PixelResolutionHistograms::~PixelResolutionHistograms ( )

virtual

delete file_ ; // no need to delete if unique_ptr<> file_ = 0;

Definition at line 431 of file PixelResolutionHistograms.cc.

References cotalphaBins_, cotbetaBins_, file_, cuy::ii, findQualityFiles::jj, GetRecoTauVFromDQM_MC_cff::kk, LOGINFO, qbinGen_, qbins_, resMultiPixelXGen_, resMultiPixelYGen_, resSinglePixelXGen_, resSinglePixelYGen_, and weOwnHistograms_.

                                                      {
  //--- Delete histograms, but only if we own them. If
  //--- they came from a file, let them be.
  //
  if (!weOwnHistograms_) {
    //--- Read the histograms from the TFile, the file will take care of them.
    file_->Close();
  } else {
    //--- We made the histograms, so first write them inthe output ROOT file and close it.
    LOGINFO << "PixelResHistoStore: Writing the histograms to the output file. "  // << filename
            << std::endl;
    file_->Write();
    file_->Close();

    // ROOT file has the ownership, and once the file is closed,
    // all of these guys are deleted.  So, we don't need to do anything.
  }  // else

  //--- Delete FastSim generators. (It's safe to delete a nullptr.)
  for (int ii = 0; ii < cotbetaBins_; ii++) {
    for (int jj = 0; jj < cotalphaBins_; jj++) {
      for (int kk = 0; kk < qbins_; kk++) {
        delete resMultiPixelXGen_[ii][jj][kk];
        delete resMultiPixelYGen_[ii][jj][kk];
      }
    }
  }
  for (int ii = 0; ii < cotbetaBins_; ii++) {
    for (int jj = 0; jj < cotalphaBins_; jj++) {
      delete resSinglePixelXGen_[ii][jj];
      delete resSinglePixelYGen_[ii][jj];
      delete qbinGen_[ii][jj];
    }
  }
}

Member Function Documentation

Fill()

int PixelResolutionHistograms::Fill	(	double	dx,
		double	dy,
		double	cotalpha,
		double	cotbeta,
		int	qbin,
		int	nxpix,
		int	nypix
	)

Definition at line 473 of file PixelResolutionHistograms.cc.

References cmtomicron, cotalphaBinWidth_, cotalphaLowEdge_, cotbetaBins_, cotbetaBinWidth_, cotbetaLowEdge_, PVValHelper::dx, PVValHelper::dy, createfilelist::int, qbinHist_, resMultiPixelXHist_, resMultiPixelYHist_, resSinglePixelXHist_, and resSinglePixelYHist_.

                                                                                           {
  int icotalpha, icotbeta, iqbin;
  icotalpha = (int)floor((cotalpha - cotalphaLowEdge_) / cotalphaBinWidth_);
  icotbeta = (int)floor((cotbeta - cotbetaLowEdge_) / cotbetaBinWidth_);
  iqbin = qbin > 2 ? 3 : qbin;
  if (icotalpha >= 0 && icotalpha < cotalphaBins_ && icotbeta >= 0 && icotbeta < cotbetaBins_) {
    qbinHist_[icotbeta][icotalpha]->Fill((double)iqbin);
    if (nxpix == 1)
      resSinglePixelXHist_[icotbeta][icotalpha]->Fill(dx / cmtomicron);
    else
      resMultiPixelXHist_[icotbeta][icotalpha][iqbin]->Fill(dx / cmtomicron);
    if (nypix == 1)
      resSinglePixelYHist_[icotbeta][icotalpha]->Fill(dy / cmtomicron);
    else
      resMultiPixelYHist_[icotbeta][icotalpha][iqbin]->Fill(dy / cmtomicron);
  }

  return 0;
}

getGeneratorX()

const SimpleHistogramGenerator * PixelResolutionHistograms::getGeneratorX	(	double	cotalpha,
		double	cotbeta,
		int	qbin,
		bool	singlex
	)

Definition at line 500 of file PixelResolutionHistograms.cc.

References cotalphaBins_, cotalphaBinWidth_, cotalphaLowEdge_, cotbetaBins_, cotbetaBinWidth_, cotbetaLowEdge_, createfilelist::int, resMultiPixelXGen_, resSinglePixelXGen_, and single().

                                                                                      {
  int icotalpha, icotbeta, iqbin;
  icotalpha = (int)floor((cotalpha - cotalphaLowEdge_) / cotalphaBinWidth_);
  icotbeta = (int)floor((cotbeta - cotbetaLowEdge_) / cotbetaBinWidth_);
  iqbin = qbin > 2 ? 3 : qbin;  // if (qbin>2) then = 3, else return qbin
  //
  //if( icotalpha >= 0 && icotalpha < cotalphaBins_ && icotbeta >= 0 && icotbeta < cotbetaBins_ ) {

  if (icotalpha < 0)
    icotalpha = 0;
  if (icotalpha >= cotalphaBins_)
    icotalpha = cotalphaBins_ - 1;

  if (icotbeta < 0)
    icotbeta = 0;
  if (icotbeta >= cotbetaBins_)
    icotbeta = cotbetaBins_ - 1;

  // At this point we are sure to return *some bin* from the 3D histogram

  if (single)
    return resSinglePixelXGen_[icotbeta][icotalpha];
  else
    return resMultiPixelXGen_[icotbeta][icotalpha][iqbin];

  // }
  //else
  //return nullptr;
}

getGeneratorY()

const SimpleHistogramGenerator * PixelResolutionHistograms::getGeneratorY	(	double	cotalpha,
		double	cotbeta,
		int	qbin,
		bool	singley
	)

Definition at line 539 of file PixelResolutionHistograms.cc.

References cotalphaBins_, cotalphaBinWidth_, cotalphaLowEdge_, cotbetaBins_, cotbetaBinWidth_, cotbetaLowEdge_, createfilelist::int, resMultiPixelYGen_, resSinglePixelYGen_, and single().

                                                                                      {
  int icotalpha, icotbeta, iqbin;
  icotalpha = (int)floor((cotalpha - cotalphaLowEdge_) / cotalphaBinWidth_);
  icotbeta = (int)floor((cotbeta - cotbetaLowEdge_) / cotbetaBinWidth_);
  iqbin = qbin > 2 ? 3 : qbin;  // if (qbin>2) then = 3, else return qbin
  //
  //if( icotalpha >= 0 && icotalpha < cotalphaBins_ && icotbeta >= 0 && icotbeta < cotbetaBins_ ) {

  if (icotalpha < 0)
    icotalpha = 0;
  if (icotalpha >= cotalphaBins_)
    icotalpha = cotalphaBins_ - 1;

  if (icotbeta < 0)
    icotbeta = 0;
  if (icotbeta >= cotbetaBins_)
    icotbeta = cotbetaBins_ - 1;

  // At this point we are sure to return *some bin* from the 3D histogram

  if (single)
    return resSinglePixelYGen_[icotbeta][icotalpha];
  else
    return resMultiPixelYGen_[icotbeta][icotalpha][iqbin];

  //}
  //else
  //return nullptr;
}

status()

int PixelResolutionHistograms::status ( void  )

inline

Definition at line 58 of file PixelResolutionHistograms.h.

References status_.

{ return status_; }

Member Data Documentation

binningHisto_

TH2F* PixelResolutionHistograms::binningHisto_

private

Definition at line 92 of file PixelResolutionHistograms.h.

Referenced by PixelResolutionHistograms().

cotalphaAxis_

TAxis* PixelResolutionHistograms::cotalphaAxis_

private

Definition at line 94 of file PixelResolutionHistograms.h.

Referenced by PixelResolutionHistograms().

cotalphaBins_

int PixelResolutionHistograms::cotalphaBins_

private

Definition at line 87 of file PixelResolutionHistograms.h.

Referenced by getGeneratorX(), getGeneratorY(), PixelResolutionHistograms(), and ~PixelResolutionHistograms().

cotalphaBinWidth_

double PixelResolutionHistograms::cotalphaBinWidth_

private

Definition at line 85 of file PixelResolutionHistograms.h.

Referenced by Fill(), getGeneratorX(), getGeneratorY(), and PixelResolutionHistograms().

cotalphaLowEdge_

double PixelResolutionHistograms::cotalphaLowEdge_

private

Definition at line 86 of file PixelResolutionHistograms.h.

Referenced by Fill(), getGeneratorX(), getGeneratorY(), and PixelResolutionHistograms().

cotbetaAxis_

TAxis* PixelResolutionHistograms::cotbetaAxis_

private

Definition at line 93 of file PixelResolutionHistograms.h.

Referenced by PixelResolutionHistograms().

cotbetaBins_

int PixelResolutionHistograms::cotbetaBins_

private

Definition at line 84 of file PixelResolutionHistograms.h.

Referenced by Fill(), getGeneratorX(), getGeneratorY(), PixelResolutionHistograms(), and ~PixelResolutionHistograms().

cotbetaBinWidth_

double PixelResolutionHistograms::cotbetaBinWidth_

private

Definition at line 82 of file PixelResolutionHistograms.h.

Referenced by Fill(), getGeneratorX(), getGeneratorY(), and PixelResolutionHistograms().

cotbetaLowEdge_

double PixelResolutionHistograms::cotbetaLowEdge_

private

Definition at line 83 of file PixelResolutionHistograms.h.

Referenced by Fill(), getGeneratorX(), getGeneratorY(), and PixelResolutionHistograms().

detType_

unsigned int PixelResolutionHistograms::detType_

private

Definition at line 79 of file PixelResolutionHistograms.h.

Referenced by PixelResolutionHistograms().

file_

std::unique_ptr<TFile> PixelResolutionHistograms::file_

private

Definition at line 115 of file PixelResolutionHistograms.h.

Referenced by PixelResolutionHistograms(), and ~PixelResolutionHistograms().

qbinGen_

SimpleHistogramGenerator* PixelResolutionHistograms::qbinGen_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX]

private

Definition at line 125 of file PixelResolutionHistograms.h.

Referenced by PixelResolutionHistograms(), and ~PixelResolutionHistograms().

qbinHist_

TH1F* PixelResolutionHistograms::qbinHist_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX]

private

Definition at line 112 of file PixelResolutionHistograms.h.

Referenced by Fill(), and PixelResolutionHistograms().

qbins_

int PixelResolutionHistograms::qbins_

private

Definition at line 89 of file PixelResolutionHistograms.h.

Referenced by PixelResolutionHistograms(), and ~PixelResolutionHistograms().

qbinWidth_

int PixelResolutionHistograms::qbinWidth_

private

Definition at line 88 of file PixelResolutionHistograms.h.

resMultiPixelXGen_

SimpleHistogramGenerator* PixelResolutionHistograms::resMultiPixelXGen_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX][QBIN_HIST_MAX]

private

Definition at line 121 of file PixelResolutionHistograms.h.

Referenced by getGeneratorX(), PixelResolutionHistograms(), and ~PixelResolutionHistograms().

resMultiPixelXHist_

TH1F* PixelResolutionHistograms::resMultiPixelXHist_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX][QBIN_HIST_MAX]

private

Definition at line 108 of file PixelResolutionHistograms.h.

Referenced by Fill(), and PixelResolutionHistograms().

resMultiPixelYGen_

SimpleHistogramGenerator* PixelResolutionHistograms::resMultiPixelYGen_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX][QBIN_HIST_MAX]

private

Definition at line 123 of file PixelResolutionHistograms.h.

Referenced by getGeneratorY(), PixelResolutionHistograms(), and ~PixelResolutionHistograms().

resMultiPixelYHist_

TH1F* PixelResolutionHistograms::resMultiPixelYHist_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX][QBIN_HIST_MAX]

private

Definition at line 110 of file PixelResolutionHistograms.h.

Referenced by Fill(), and PixelResolutionHistograms().

resSinglePixelXGen_

SimpleHistogramGenerator* PixelResolutionHistograms::resSinglePixelXGen_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX]

private

Definition at line 122 of file PixelResolutionHistograms.h.

Referenced by getGeneratorX(), PixelResolutionHistograms(), and ~PixelResolutionHistograms().

resSinglePixelXHist_

TH1F* PixelResolutionHistograms::resSinglePixelXHist_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX]

private

Definition at line 109 of file PixelResolutionHistograms.h.

Referenced by Fill(), and PixelResolutionHistograms().

resSinglePixelYGen_

SimpleHistogramGenerator* PixelResolutionHistograms::resSinglePixelYGen_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX]

private

Definition at line 124 of file PixelResolutionHistograms.h.

Referenced by getGeneratorY(), PixelResolutionHistograms(), and ~PixelResolutionHistograms().

resSinglePixelYHist_

TH1F* PixelResolutionHistograms::resSinglePixelYHist_[COTBETA_HIST_MAX][COTALPHA_HIST_MAX]

private

Definition at line 111 of file PixelResolutionHistograms.h.

Referenced by Fill(), and PixelResolutionHistograms().

status_

int PixelResolutionHistograms::status_

private

Definition at line 118 of file PixelResolutionHistograms.h.

Referenced by PixelResolutionHistograms(), and status().

weOwnHistograms_

bool PixelResolutionHistograms::weOwnHistograms_

private

Definition at line 76 of file PixelResolutionHistograms.h.

Referenced by ~PixelResolutionHistograms().