#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 21 of file PixelResolutionHistograms.h.

Constructor & Destructor Documentation

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 67 of file PixelResolutionHistograms.cc.

References binningHisto_, cotalphaBins_, cotalphaBinWidth_, cotalphaLowEdge_, cotbetaBins_, cotbetaBinWidth_, cotbetaLowEdge_, detType_, file_, trackerHits::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::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 194 of file PixelResolutionHistograms.cc.

References binningHisto_, cotalphaAxis_, cotalphaBins_, cotalphaBinWidth_, cotalphaLowEdge_, cotbetaAxis_, cotbetaBins_, cotbetaBinWidth_, cotbetaLowEdge_, detType_, file_, trackerHits::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 ( )

virtual

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

Definition at line 418 of file PixelResolutionHistograms.cc.

References cotalphaBins_, cotbetaBins_, file_, Fill(), 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

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_, getGeneratorX(), createfilelist::int, qbinHist_, resMultiPixelXHist_, resMultiPixelYHist_, resSinglePixelXHist_, and resSinglePixelYHist_.

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

 {
   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;  
 }

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

Definition at line 504 of file PixelResolutionHistograms.cc.

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

Referenced by Fill(), and status().

 {
   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;
 }

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

Definition at line 543 of file PixelResolutionHistograms.cc.

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

Referenced by getGeneratorX(), and status().

 {
   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;
 }