#include <CalibTracker/SiPixelGainCalibrationAnalysis/interface/SiPixelGainCalibrationAnalysis.h>

Inheritance diagram for SiPixelGainCalibrationAnalysis:

Public Member Functions
bool	checkCorrectCalibrationType () override

bool	doFits (uint32_t detid, std::vector< SiPixelCalibDigi >::const_iterator ipix) override

void	doSetup (const edm::ParameterSet &)

	SiPixelGainCalibrationAnalysis (const edm::ParameterSet &iConfig)

	~SiPixelGainCalibrationAnalysis () override

Public Member Functions inherited from SiPixelOfflineCalibAnalysisBase
void	addTF1ToDQMMonitoringElement (MonitorElement ele, TF1 func)

MonitorElement *	bookDQMHistogram1D (uint32_t detid, std::string name, std::string title, int nchX, double lowX, double highX)

MonitorElement *	bookDQMHistogram1D (uint32_t detid, std::string name, std::string title, int nchX, float *xbinsize)

MonitorElement *	bookDQMHistogram2D (uint32_t detid, std::string name, std::string title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)

MonitorElement *	bookDQMHistoPlaquetteSummary2D (uint32_t detid, std::string name, std::string title)

std::vector< uint32_t > &	getRunNumbers ()

bool	setDQMDirectory (std::string dirName)

bool	setDQMDirectory (uint32_t detID)

	SiPixelOfflineCalibAnalysisBase ()

	SiPixelOfflineCalibAnalysisBase (const edm::ParameterSet &)

std::string	translateDetIdToString (uint32_t detid)

	~SiPixelOfflineCalibAnalysisBase () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
std::vector< float >	CalculateAveragePerColumn (uint32_t detid, std::string label)

void	calibrationEnd () override

void	calibrationSetup (const edm::EventSetup &iSetup) override

void	fillDatabase ()

void	newDetID (uint32_t detid) override

void	printSummary ()

Private Attributes
bool	appendMode_

bool	bookBIGCalibPayload_

std::map< uint32_t, std::map< std::string, MonitorElement * > >	bookkeeper_

std::map< uint32_t, std::map< std::string, MonitorElement * > >	bookkeeper_pixels_

double	chi2ProbThreshold_

double	chi2Threshold_

edm::ParameterSet	conf_

uint32_t	currentDetID_

bool	filldb_

std::string	fitfunction_

TF1 *	func_

float	gainhi_

float	gainlow_

TGraphErrors *	graph_

bool	ignoreMode_

std::vector< uint32_t >	listofdetids_

double	maxChi2InHist_

double	maxGainInHist_

uint16_t	min_nentries_

int	nfitparameters_

float	pedhi_

float	pedlow_

double	plateau_max_slope_

std::string	recordName_

bool	reject_badpoints_

double	reject_badpoints_frac_

bool	reject_first_point_

bool	reject_plateaupoints_

bool	reject_single_entries_

bool	saveALLHistograms_

bool	savePixelHists_

double	scalarVcalHigh_VcalLow_

int *	statusNumbers_

bool	sum_gain_cols_

bool	sum_ped_cols_

std::ofstream	summary_

bool	useVcalHigh_

bool	writeSummary_

Additional Inherited Members
Public Types inherited from SiPixelOfflineCalibAnalysisBase
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from SiPixelOfflineCalibAnalysisBase
static const std::vector< short > *	getVcalValues ()

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Static Public Attributes inherited from SiPixelOfflineCalibAnalysisBase
static TF1 *	fitFunction_ = nullptr

Protected Member Functions inherited from SiPixelOfflineCalibAnalysisBase
uint32_t &	EventNumber ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Protected Attributes inherited from SiPixelOfflineCalibAnalysisBase
edm::ESHandle< SiPixelCalibConfiguration >	calib_

std::string	calibrationMode_

edm::ESHandle< TrackerGeometry >	geom_

short	nTriggers_

edm::ESHandle< SiPixelFedCablingMap >	theCablingMap_

Static Protected Attributes inherited from SiPixelOfflineCalibAnalysisBase
static std::vector< short >	vCalValues_

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 47 of file SiPixelGainCalibrationAnalysis.h.

Constructor & Destructor Documentation

◆ SiPixelGainCalibrationAnalysis()

SiPixelGainCalibrationAnalysis::SiPixelGainCalibrationAnalysis ( const edm::ParameterSet & iConfig )

explicit

Definition at line 34 of file SiPixelGainCalibrationAnalysis.cc.

     : SiPixelOfflineCalibAnalysisBase(iConfig),
       conf_(iConfig),
       bookkeeper_(),
       bookkeeper_pixels_(),
       nfitparameters_(iConfig.getUntrackedParameter<int>("numberOfFitParameters", 2)),
       fitfunction_(iConfig.getUntrackedParameter<std::string>("fitFunctionRootFormula", "pol1")),
       listofdetids_(conf_.getUntrackedParameter<std::vector<uint32_t> >("listOfDetIDs")),
       ignoreMode_(conf_.getUntrackedParameter<bool>("ignoreMode", false)),
       reject_plateaupoints_(iConfig.getUntrackedParameter<bool>("suppressPlateauInFit", true)),
       reject_single_entries_(iConfig.getUntrackedParameter<bool>("suppressPointsWithOneEntryOrLess", true)),
       plateau_max_slope_(iConfig.getUntrackedParameter<double>("plateauSlopeMax", 1.0)),
       reject_first_point_(iConfig.getUntrackedParameter<bool>("rejectVCalZero", true)),
       reject_badpoints_frac_(iConfig.getUntrackedParameter<double>("suppressZeroAndPlateausInFitFrac", 0)),
       bookBIGCalibPayload_(iConfig.getUntrackedParameter<bool>("saveFullPayloads", false)),
       savePixelHists_(iConfig.getUntrackedParameter<bool>("savePixelLevelHists", false)),
       chi2Threshold_(iConfig.getUntrackedParameter<double>("minChi2NDFforHistSave", 10)),
       chi2ProbThreshold_(iConfig.getUntrackedParameter<double>("minChi2ProbforHistSave", 0.05)),
       maxGainInHist_(iConfig.getUntrackedParameter<double>("maxGainInHist", 10)),
       maxChi2InHist_(iConfig.getUntrackedParameter<double>("maxChi2InHist", 25)),
       saveALLHistograms_(iConfig.getUntrackedParameter<bool>("saveAllHistograms", false)),
  
       filldb_(iConfig.getUntrackedParameter<bool>("writeDatabase", false)),
       writeSummary_(iConfig.getUntrackedParameter<bool>("writeSummary", true)),
       recordName_(conf_.getParameter<std::string>("record")),
  
       appendMode_(conf_.getUntrackedParameter<bool>("appendMode", true)),
       /*theGainCalibrationDbInput_(0),
   theGainCalibrationDbInputOffline_(0),
   theGainCalibrationDbInputHLT_(0),
   theGainCalibrationDbInputService_(iConfig),*/
       gainlow_(10.),
       gainhi_(0.),
       pedlow_(255.),
       pedhi_(0.),
       useVcalHigh_(conf_.getParameter<bool>("useVCALHIGH")),
       scalarVcalHigh_VcalLow_(conf_.getParameter<double>("vcalHighToLowConversionFac")) {
   if (reject_single_entries_)
     min_nentries_ = 1;
   else
     min_nentries_ = 0;
   ::putenv((char *)"CORAL_AUTH_USER=me");
   ::putenv((char *)"CORAL_AUTH_PASSWORD=test");
   edm::LogInfo("SiPixelGainCalibrationAnalysis") << "now using fit function " << fitfunction_ << ", which has "
                                                  << nfitparameters_ << " free parameters. " << std::endl;
   func_ = new TF1("func", fitfunction_.c_str(), 0, 256 * scalarVcalHigh_VcalLow_);
   graph_ = new TGraphErrors();
   currentDetID_ = 0;
   summary_.open("SummaryPerDetID.txt");
   statusNumbers_ = new int[10];
   for (int ii = 0; ii < 10; ii++)
     statusNumbers_[ii] = 0;
 }

References currentDetID_, fitfunction_, func_, graph_, cuy::ii, min_nentries_, nfitparameters_, reject_single_entries_, scalarVcalHigh_VcalLow_, statusNumbers_, summary_, and tools::TF1.

◆ ~SiPixelGainCalibrationAnalysis()

SiPixelGainCalibrationAnalysis::~SiPixelGainCalibrationAnalysis ( )

override

Definition at line 88 of file SiPixelGainCalibrationAnalysis.cc.

88 {}

Member Function Documentation

◆ CalculateAveragePerColumn()

std::vector< float > SiPixelGainCalibrationAnalysis::CalculateAveragePerColumn	(	uint32_t	detid,
		std::string	label
	)

private

Definition at line 93 of file SiPixelGainCalibrationAnalysis.cc.

                                                                                                           {
   std::vector<float> result;
   int ncols = bookkeeper_[detid][label]->getNbinsX();
   int nrows = bookkeeper_[detid][label]->getNbinsY();
   for (int icol = 1; icol <= ncols; ++icol) {
     float val = 0;
     float ntimes = 0;
     for (int irow = 1; irow <= nrows; ++irow) {
       val += bookkeeper_[detid][label]->getBinContent(icol, irow);
       ntimes++;
     }
     val /= ntimes;
     result.push_back(val);
   }
   return result;
 }

References bookkeeper_, label, hgcalPlots::ncols, mps_fire::result, and heppy_batch::val.

Referenced by printSummary().

◆ calibrationEnd()

void SiPixelGainCalibrationAnalysis::calibrationEnd ( )

overrideprivatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 167 of file SiPixelGainCalibrationAnalysis.cc.

                                                     {
   if (writeSummary_)
     printSummary();
  
   // this is where we loop over all histograms and save the database objects
   if (filldb_)
     fillDatabase();
 }

References fillDatabase(), filldb_, printSummary(), and writeSummary_.

◆ calibrationSetup()

void SiPixelGainCalibrationAnalysis::calibrationSetup ( const edm::EventSetup & iSetup )

overrideprivatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 125 of file SiPixelGainCalibrationAnalysis.cc.

125 {}

◆ checkCorrectCalibrationType()

bool SiPixelGainCalibrationAnalysis::checkCorrectCalibrationType ( )

overridevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 110 of file SiPixelGainCalibrationAnalysis.cc.

                                                                  {
   if (calibrationMode_ == "GainCalibration")
     return true;
   else if (ignoreMode_ == true)
     return true;
   else if (calibrationMode_ == "unknown") {
     edm::LogInfo("SiPixelGainCalibrationAnalysis")
         << "calibration mode is: " << calibrationMode_ << ", continuing anyway...";
     return true;
   } else {
     //    edm::LogError("SiPixelGainCalibrationAnalysis") << "unknown calibration mode for Gain calibration, should be \"Gain\" and is \"" << calibrationMode_ << "\"";
   }
   return false;
 }

References SiPixelOfflineCalibAnalysisBase::calibrationMode_, and ignoreMode_.

◆ doFits()

bool SiPixelGainCalibrationAnalysis::doFits	(	uint32_t	detid,
		std::vector< SiPixelCalibDigi >::const_iterator	ipix
	)

overridevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 185 of file SiPixelGainCalibrationAnalysis.cc.

                                                                                                           {
   float lowmeanval = 255;
   float highmeanval = 0;
   bool makehistopersistent = saveALLHistograms_;
   std::vector<uint32_t>::const_iterator detidfinder = find(listofdetids_.begin(), listofdetids_.end(), detid);
   if (detidfinder != listofdetids_.end())
     makehistopersistent = true;
   // first, fill the input arrays to the TLinearFitter.
   double xvals[257];
   double yvals[256];
   double yerrvals[256];
   double xvalsall[257];
   double yvalsall[256];
   double yerrvalsall[256];
   int npoints = 0;
   int nallpoints = 0;
   bool use_point = true;
   int status = 0;
   statusNumbers_[9]++;
  
   bookkeeper_[detid]["status_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, 0);
   if (writeSummary_ && detid != currentDetID_) {
     currentDetID_ = detid;
     summary_ << endl << "DetId_" << currentDetID_ << endl;
   }
  
   for (uint32_t ii = 0; ii < ipix->getnpoints() && ii < 200; ii++) {
     //    std::cout << ipix->getsum(ii) << " " << ipix->getnentries(ii) << " " << ipix->getsumsquares(ii) << std::endl;
     nallpoints++;
     use_point = true;
     if (useVcalHigh_) {
       xvalsall[ii] = vCalValues_[ii] * scalarVcalHigh_VcalLow_;
     } else
       xvalsall[ii] = vCalValues_[ii];
     yerrvalsall[ii] = yvalsall[ii] = 0;
  
     if (ipix->getnentries(ii) > min_nentries_) {
       yvalsall[ii] = ipix->getsum(ii) / (float)ipix->getnentries(ii);
       yerrvalsall[ii] = ipix->getsumsquares(ii) / (float)(ipix->getnentries(ii));
       yerrvalsall[ii] -= pow(yvalsall[ii], 2);
       yerrvalsall[ii] = sqrt(yerrvalsall[ii]) / sqrt(ipix->getnentries(ii));
  
       if (yvalsall[ii] < lowmeanval)
         lowmeanval = yvalsall[ii];
       if (yvalsall[ii] > highmeanval)
         highmeanval = yvalsall[ii];
     }
   }
  
   // calculate plateau value from last 4 entries
   double plateauval = 0;
   bool noPlateau = false;
   if (nallpoints >= 4) {
     for (int ii = nallpoints - 1; ii > nallpoints - 5; --ii)
       plateauval += yvalsall[ii];
     plateauval /= 4;
     for (int ii = nallpoints - 1; ii > nallpoints - 5; --ii) {
       if (fabs(yvalsall[ii] - plateauval) > 5) {
         plateauval = 255;
         noPlateau = true;
         continue;
       }
     }
  
     int NbofPointsInPlateau = 0;
     for (int ii = 0; ii < nallpoints; ++ii)
       if (fabs(yvalsall[ii] - plateauval) < 10 || yvalsall[ii] == 0)
         NbofPointsInPlateau++;
     //summary_<<"row_"<<ipix->row()<<" col_"<<ipix->col()<<"   "<<plateauval<<"  "<<NbofPointsInPlateau<<"  "<<nallpoints<<endl;
     if (NbofPointsInPlateau >= (nallpoints - 2)) {
       status = -2;
       bookkeeper_[detid]["status_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, status);
       if (writeSummary_) {
         summary_ << "row_" << ipix->row() << " col_" << ipix->col() << " status_" << status << endl;
         statusNumbers_[abs(status)]++;
       }
       return false;
     }
   } else
     plateauval = 255;
  
   double maxgoodvalinfit = plateauval * (1. - reject_badpoints_frac_);
   npoints = 0;
   for (int ii = 0; ii < nallpoints; ++ii) {
     // now selecting the appropriate points for the fit.
     use_point = true;
     if (reject_first_point_ && xvalsall[ii] < 0.1)
       use_point = false;
     if (ipix->getnentries(ii) <= min_nentries_ && reject_single_entries_)
       use_point = false;
     if (ipix->getnentries(ii) == 0 && reject_badpoints_)
       use_point = false;
     if (yvalsall[ii] > maxgoodvalinfit && !noPlateau)
       use_point = false;
     if (ii > 1 && fabs(yvalsall[ii] - yvalsall[ii - 1]) < 5. && yvalsall[ii] > 0.8 * maxgoodvalinfit &&
         reject_plateaupoints_) {
       use_point = false;
       break;
     }
  
     if (use_point) {
       xvals[npoints] = xvalsall[ii];
       yvals[npoints] = yvalsall[ii];
       yerrvals[npoints] = yerrvalsall[ii];
       npoints++;
     }
   }
  
   float chi2, slope, intercept, prob, slopeerror, intercepterror;
   prob = chi2 = -1;
   slope = intercept = slopeerror = intercepterror = 0;
  
   // now check on number of points. If bad just start taking the first 4:
  
   if (npoints < 4) {
     npoints = 0;
     for (int ii = 0; ii < nallpoints && npoints < 4 && yvalsall[ii] < plateauval * 0.97; ++ii) {
       if (yvalsall[ii] > 0) {
         if (ii > 0 && yvalsall[ii] - yvalsall[ii - 1] < 0.1)
           continue;
         xvals[npoints] = xvalsall[ii];
         yvals[npoints] = yvalsall[ii];
         yerrvals[npoints] = yerrvalsall[ii];
         npoints++;
       }
     }
   }
   if (npoints < 2) {
     status = -7;
     bookkeeper_[detid]["status_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, status);
     if (writeSummary_) {
       summary_ << "row_" << ipix->row() << " col_" << ipix->col() << " status_" << status << endl;
       statusNumbers_[abs(status)]++;
     }
     std::ostringstream pixelinfo;
     pixelinfo << "GainCurve_row_" << ipix->row() << "_col_" << ipix->col();
     std::string tempname = translateDetIdToString(detid);
     tempname += "_";
     tempname += pixelinfo.str();
     setDQMDirectory(detid);
     bookkeeper_pixels_[detid][pixelinfo.str()] = bookDQMHistogram1D(
         detid, pixelinfo.str(), tempname, 105 * nallpoints, xvalsall[0], xvalsall[nallpoints - 1] * 1.05);
     for (int ii = 0; ii < nallpoints; ++ii)
       bookkeeper_pixels_[detid][pixelinfo.str()]->Fill(xvalsall[ii], yvalsall[ii]);
     return false;
   }
  
   //  std::cout << "starting fit!" << std::endl;
   graph_->Set(npoints);
  
   func_->SetParameter(0, 50.);
   func_->SetParameter(1, 0.25);
   for (int ipointtemp = 0; ipointtemp < npoints; ++ipointtemp) {
     graph_->SetPoint(ipointtemp, xvals[ipointtemp], yvals[ipointtemp]);
     graph_->SetPointError(ipointtemp, 0, yerrvals[ipointtemp]);
   }
   Int_t tempresult = graph_->Fit(func_, "FQ0N");
   slope = func_->GetParameter(1);
   slopeerror = func_->GetParError(1);
   intercept = func_->GetParameter(0);
   intercepterror = func_->GetParError(0);
   chi2 = func_->GetChisquare() / ((float)npoints - func_->GetNpar());
   prob = TMath::Prob(func_->GetChisquare(), npoints - func_->GetNpar());
   size_t ntimes = 0;
   while ((isnan(slope) || isnan(intercept)) && ntimes < 10) {
     ntimes++;
     makehistopersistent = true;
     //    std::cout << slope << " " << intercept << " " << prob << std::endl;
     edm::LogWarning("SiPixelGainCalibrationAnalysis") << "impossible to fit values, try " << ntimes << ": ";
     for (int ii = 0; ii < npoints; ++ii) {
       edm::LogWarning("SiPixelGainCalibrationAnalysis")
           << "vcal " << xvals[ii] << " response: " << yvals[ii] << "+/-" << yerrvals[ii] << std::endl;
     }
     tempresult = graph_->Fit(func_, "FQ0NW");
     slope = func_->GetParameter(1);
     slopeerror = func_->GetParError(1);
     intercept = func_->GetParameter(0);
     intercepterror = func_->GetParError(0);
     chi2 = func_->GetChisquare() / ((float)npoints - func_->GetNpar());
     prob = TMath::Prob(func_->GetChisquare(), npoints - func_->GetNpar());
   }
  
   if (tempresult == 0)
     status = 1;
   else
     status = 0;
   if (slope != 0)
     slope = 1. / slope;
   if (isnan(slope) || isnan(intercept)) {
     status = -6;
     bookkeeper_[detid]["status_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, status);
     if (writeSummary_) {
       summary_ << "row_" << ipix->row() << " col_" << ipix->col() << " status_" << status << endl;
       statusNumbers_[abs(status)]++;
     }
     //return false;
   }
   if (chi2 > chi2Threshold_ && chi2Threshold_ >= 0)
     status = 5;
   if (prob < chi2ProbThreshold_)
     status = 5;
   if (noPlateau)
     status = 3;
   if (nallpoints < 4)
     status = -7;
   if (TMath::Abs(slope > maxGainInHist_) || slope < 0)
     status = -8;
   if (status != 1)
     makehistopersistent = true;
   statusNumbers_[abs(status)]++;
  
   if (slope < gainlow_)
     gainlow_ = slope;
   if (slope > gainhi_)
     gainhi_ = slope;
   if (intercept > pedhi_)
     pedhi_ = intercept;
   if (intercept < pedlow_)
     pedlow_ = intercept;
   bookkeeper_[detid]["gain_1d"]->Fill(slope);
   if (slope > maxGainInHist_) {
     makehistopersistent = true;
     edm::LogWarning("SiPixelGainCalibration")
         << "For DETID " << detid << "pixel row,col " << ipix->row() << "," << ipix->col() << " Gain was measured to be "
         << slope << " which is outside the range of the summary plot (" << maxGainInHist_ << ") !!!! " << std::endl;
   }
   bookkeeper_[detid]["dynamicrange_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, highmeanval - lowmeanval);
   bookkeeper_[detid]["plateau_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, highmeanval);
   bookkeeper_[detid]["gain_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, slope);
   bookkeeper_[detid]["errorgain_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, slopeerror);
   bookkeeper_[detid]["ped_1d"]->Fill(intercept);
   bookkeeper_[detid]["ped_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, intercept);
   bookkeeper_[detid]["errorped_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, intercepterror);
   bookkeeper_[detid]["chi2_1d"]->Fill(chi2);
   bookkeeper_[detid]["chi2_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, chi2);
   bookkeeper_[detid]["prob_1d"]->Fill(prob);
   bookkeeper_[detid]["prob_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, prob);
   bookkeeper_[detid]["lowpoint_1d"]->Fill(xvals[0]);
   bookkeeper_[detid]["lowpoint_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, xvals[0]);
   bookkeeper_[detid]["highpoint_1d"]->Fill(xvals[npoints - 1]);
   bookkeeper_[detid]["highpoint_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, xvals[npoints - 1]);
   bookkeeper_[detid]["nfitpoints_1d"]->Fill(npoints);
   bookkeeper_[detid]["endpoint_1d"]->Fill((255 - intercept) * slope);
   bookkeeper_[detid]["status_2d"]->setBinContent(ipix->col() + 1, ipix->row() + 1, status);
  
   if (!savePixelHists_)
     return true;
   if (detidfinder == listofdetids_.end() && !listofdetids_.empty())
     return true;
   if (makehistopersistent) {
     std::ostringstream pixelinfo;
     pixelinfo << "GainCurve_row_" << ipix->row() << "_col_" << ipix->col();
     std::string tempname = translateDetIdToString(detid);
     tempname += "_";
     tempname += pixelinfo.str();
  
     // and book the histo
     // fill the last value of the vcal array...
  
     setDQMDirectory(detid);
     bookkeeper_pixels_[detid][pixelinfo.str()] = bookDQMHistogram1D(
         detid, pixelinfo.str(), tempname, 105 * nallpoints, xvalsall[0], xvalsall[nallpoints - 1] * 1.05);
  
     edm::LogInfo("SiPixelGainCalibrationAnalysis")
         << "now saving histogram for pixel " << tempname << ", gain = " << slope << ", pedestal = " << intercept
         << ", chi2/NDF=" << chi2 << "(prob:" << prob << "), fit status " << status;
     for (int ii = 0; ii < nallpoints; ++ii) {
       //      std::cout << xvalsall[ii]<<","<<yvalsall[ii]<< " " << tempfloats[ii] << std::endl;
       bookkeeper_pixels_[detid][pixelinfo.str()]->Fill(xvalsall[ii], yvalsall[ii]);
     }
  
     //    addTF1ToDQMMonitoringElement(bookkeeper_pixels_[detid][pixelinfo.str()],func_);
  
     if (writeSummary_) {
       summary_ << "row_" << ipix->row() << " col_" << ipix->col();
       summary_ << " status_" << status;
       summary_ << endl;
  
       //std::cout<<detid<<"  "<<"row " <<ipix->row()<<" col "<<ipix->col()<<"  "<<status<<"  "<<chi2<<"  "<<prob<<"  "<<npoints<<"  "<<xvals[0]<<"  "<<xvals[npoints-1]<<"  "<<plateauval<<std::endl;
     }
   }
   return true;
 }

References funct::abs(), Abs(), SiPixelOfflineCalibAnalysisBase::bookDQMHistogram1D(), bookkeeper_, bookkeeper_pixels_, hltPixelTracks_cff::chi2, chi2ProbThreshold_, chi2Threshold_, currentDetID_, spr::find(), dqmMemoryStats::float, func_, gainhi_, gainlow_, graph_, cuy::ii, CommonMethods::isnan(), listofdetids_, maxGainInHist_, min_nentries_, npoints, pedhi_, pedlow_, funct::pow(), TtFullHadEvtBuilder_cfi::prob, reject_badpoints_, reject_badpoints_frac_, reject_first_point_, reject_plateaupoints_, reject_single_entries_, saveALLHistograms_, savePixelHists_, scalarVcalHigh_VcalLow_, SiPixelOfflineCalibAnalysisBase::setDQMDirectory(), slope, mathSSE::sqrt(), mps_update::status, statusNumbers_, AlCaHLTBitMon_QueryRunRegistry::string, summary_, SiPixelOfflineCalibAnalysisBase::translateDetIdToString(), useVcalHigh_, SiPixelOfflineCalibAnalysisBase::vCalValues_, and writeSummary_.

◆ doSetup()

void SiPixelGainCalibrationAnalysis::doSetup ( const edm::ParameterSet & )

◆ fillDatabase()

void SiPixelGainCalibrationAnalysis::fillDatabase ( )

private

Definition at line 176 of file SiPixelGainCalibrationAnalysis.cc.

                                                   {
   // only create when necessary.
   // process the minimum and maximum gain & ped values...
   edm::LogError("SiPixelGainCalibration::fillDatabase()")
       << "PLEASE do not fill the database directly from the gain calibration analyzer. This function is currently "
          "disabled and no DB payloads will be produced!"
       << std::endl;
 }

Referenced by calibrationEnd().

◆ newDetID()

void SiPixelGainCalibrationAnalysis::newDetID ( uint32_t detid )

overrideprivatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 469 of file SiPixelGainCalibrationAnalysis.cc.

                                                             {
   setDQMDirectory(detid);
   std::string tempname = translateDetIdToString(detid);
   bookkeeper_[detid]["gain_1d"] = bookDQMHistogram1D(detid, "Gain1d", "gain for " + tempname, 100, 0., maxGainInHist_);
   bookkeeper_[detid]["gain_2d"] = bookDQMHistoPlaquetteSummary2D(detid, "Gain2d", "gain for " + tempname);
   bookkeeper_[detid]["errorgain_2d"] =
       bookDQMHistoPlaquetteSummary2D(detid, "ErrorGain2d", "error on gain for " + tempname);
   bookkeeper_[detid]["ped_1d"] = bookDQMHistogram1D(detid, "Pedestal1d", "pedestal for " + tempname, 256, 0., 256.0);
   bookkeeper_[detid]["ped_2d"] = bookDQMHistoPlaquetteSummary2D(detid, "Pedestal2d", "pedestal for " + tempname);
   bookkeeper_[detid]["errorped_2d"] =
       bookDQMHistoPlaquetteSummary2D(detid, "ErrorPedestal2d", "error on pedestal for " + tempname);
   bookkeeper_[detid]["chi2_1d"] =
       bookDQMHistogram1D(detid, "GainChi2NDF1d", "#chi^{2}/NDOF for " + tempname, 100, 0., maxChi2InHist_);
   bookkeeper_[detid]["chi2_2d"] =
       bookDQMHistoPlaquetteSummary2D(detid, "GainChi2NDF2d", "#chi^{2}/NDOF for " + tempname);
   bookkeeper_[detid]["prob_1d"] =
       bookDQMHistogram1D(detid, "GainChi2Prob1d", "P(#chi^{2},NDOF) for " + tempname, 100, 0., 1.0);
   bookkeeper_[detid]["prob_2d"] =
       bookDQMHistoPlaquetteSummary2D(detid, "GainChi2Prob2d", "P(#chi^{2},NDOF) for " + tempname);
   bookkeeper_[detid]["status_2d"] =
       bookDQMHistoPlaquetteSummary2D(detid, "GainFitResult2d", "Fit result for " + tempname);
   bookkeeper_[detid]["endpoint_1d"] = bookDQMHistogram1D(
       detid, "GainEndPoint1d", "point where fit meets ADC=255 for " + tempname, 256, 0., 256. * scalarVcalHigh_VcalLow_);
   bookkeeper_[detid]["lowpoint_1d"] = bookDQMHistogram1D(
       detid, "GainLowPoint1d", "lowest fit point for " + tempname, 256, 0., 256. * scalarVcalHigh_VcalLow_);
   bookkeeper_[detid]["highpoint_1d"] = bookDQMHistogram1D(
       detid, "GainHighPoint1d", "highest fit point for " + tempname, 256, 0., 256. * scalarVcalHigh_VcalLow_);
   bookkeeper_[detid]["nfitpoints_1d"] =
       bookDQMHistogram1D(detid, "GainNPoints1d", "number of fit point for " + tempname, 20, 0., 20);
   bookkeeper_[detid]["dynamicrange_2d"] = bookDQMHistoPlaquetteSummary2D(
       detid, "GainDynamicRange2d", "Difference lowest and highest points on gain curve for " + tempname);
   bookkeeper_[detid]["lowpoint_2d"] =
       bookDQMHistoPlaquetteSummary2D(detid, "GainLowPoint2d", "lowest fit point for " + tempname);
   bookkeeper_[detid]["highpoint_2d"] =
       bookDQMHistoPlaquetteSummary2D(detid, "GainHighPoint2d", "highest fit point for " + tempname);
   bookkeeper_[detid]["plateau_2d"] =
       bookDQMHistoPlaquetteSummary2D(detid, "GainSaturate2d", "Highest points on gain curve for " + tempname);
 }

References SiPixelOfflineCalibAnalysisBase::bookDQMHistogram1D(), SiPixelOfflineCalibAnalysisBase::bookDQMHistoPlaquetteSummary2D(), bookkeeper_, maxChi2InHist_, maxGainInHist_, scalarVcalHigh_VcalLow_, SiPixelOfflineCalibAnalysisBase::setDQMDirectory(), AlCaHLTBitMon_QueryRunRegistry::string, and SiPixelOfflineCalibAnalysisBase::translateDetIdToString().

◆ printSummary()

void SiPixelGainCalibrationAnalysis::printSummary ( )

private

Definition at line 127 of file SiPixelGainCalibrationAnalysis.cc.

                                                   {
   uint32_t detid = 0;
   for (std::map<uint32_t, std::map<std::string, MonitorElement *> >::const_iterator idet = bookkeeper_.begin();
        idet != bookkeeper_.end();
        ++idet) {
     if (detid == idet->first)
       continue;  // only do things once per detid
     detid = idet->first;
     std::vector<float> gainvec = CalculateAveragePerColumn(detid, "gain_2d");
     std::vector<float> pedvec = CalculateAveragePerColumn(detid, "ped_2d");
     std::vector<float> chi2vec = CalculateAveragePerColumn(detid, "chi2_2d");
     std::ostringstream summarytext;
  
     summarytext << "Summary for det ID " << detid << "(" << translateDetIdToString(detid) << ")\n";
     summarytext << "\t Following: values per column: column #, gain, pedestal, chi2\n";
     for (uint32_t i = 0; i < gainvec.size(); i++)
       summarytext << "\t " << i << " \t" << gainvec[i] << " \t" << pedvec[i] << " \t" << chi2vec[i] << "\n";
     summarytext << "\t list of pixels with high chi2 (chi2> " << chi2Threshold_ << "): \n";
  
     for (std::map<std::string, MonitorElement *>::const_iterator ipix = bookkeeper_pixels_[detid].begin();
          ipix != bookkeeper_pixels_[detid].end();
          ++ipix)
       summarytext << "\t " << ipix->first << "\n";
     edm::LogInfo("SiPixelGainCalibrationAnalysis") << summarytext.str() << std::endl;
   }
   if (summary_.is_open()) {
     summary_.close();
     summary_.open("Summary.txt");
     summary_ << "Total Number of Pixel computed :" << statusNumbers_[9] << endl;
     summary_ << "Number of pixel tagged with status :" << endl;
     for (int ii = 0; ii < 9; ii++)
       summary_ << ii << " -> " << statusNumbers_[ii] << " ~ "
                << double(statusNumbers_[ii]) / double(statusNumbers_[9]) * 100. << " %" << endl;
  
     summary_.close();
   }
 }

References bookkeeper_, bookkeeper_pixels_, CalculateAveragePerColumn(), chi2Threshold_, mps_fire::i, cuy::ii, genParticles_cff::map, statusNumbers_, summary_, and SiPixelOfflineCalibAnalysisBase::translateDetIdToString().

Referenced by calibrationEnd().

Member Data Documentation

◆ appendMode_

bool SiPixelGainCalibrationAnalysis::appendMode_

private

Definition at line 101 of file SiPixelGainCalibrationAnalysis.h.

◆ bookBIGCalibPayload_

bool SiPixelGainCalibrationAnalysis::bookBIGCalibPayload_

private

Definition at line 87 of file SiPixelGainCalibrationAnalysis.h.

◆ bookkeeper_

std::map<uint32_t, std::map<std::string, MonitorElement *> > SiPixelGainCalibrationAnalysis::bookkeeper_

private

Definition at line 68 of file SiPixelGainCalibrationAnalysis.h.

Referenced by CalculateAveragePerColumn(), doFits(), newDetID(), and printSummary().

◆ bookkeeper_pixels_

std::map<uint32_t, std::map<std::string, MonitorElement *> > SiPixelGainCalibrationAnalysis::bookkeeper_pixels_

private

Definition at line 69 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and printSummary().

◆ chi2ProbThreshold_

double SiPixelGainCalibrationAnalysis::chi2ProbThreshold_

private

Definition at line 90 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ chi2Threshold_

double SiPixelGainCalibrationAnalysis::chi2Threshold_

private

Definition at line 89 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and printSummary().

◆ conf_

edm::ParameterSet SiPixelGainCalibrationAnalysis::conf_

private

Definition at line 66 of file SiPixelGainCalibrationAnalysis.h.

◆ currentDetID_

uint32_t SiPixelGainCalibrationAnalysis::currentDetID_

private

Definition at line 118 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

◆ filldb_

bool SiPixelGainCalibrationAnalysis::filldb_

private

Definition at line 96 of file SiPixelGainCalibrationAnalysis.h.

Referenced by calibrationEnd().

◆ fitfunction_

std::string SiPixelGainCalibrationAnalysis::fitfunction_

private

Definition at line 73 of file SiPixelGainCalibrationAnalysis.h.

Referenced by SiPixelGainCalibrationAnalysis().

◆ func_

TF1* SiPixelGainCalibrationAnalysis::func_

private

Definition at line 74 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

◆ gainhi_

float SiPixelGainCalibrationAnalysis::gainhi_

private

Definition at line 109 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ gainlow_

float SiPixelGainCalibrationAnalysis::gainlow_

private

Definition at line 108 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ graph_

TGraphErrors* SiPixelGainCalibrationAnalysis::graph_

private

Definition at line 75 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

◆ ignoreMode_

bool SiPixelGainCalibrationAnalysis::ignoreMode_

private

Definition at line 78 of file SiPixelGainCalibrationAnalysis.h.

Referenced by checkCorrectCalibrationType().

◆ listofdetids_

std::vector<uint32_t> SiPixelGainCalibrationAnalysis::listofdetids_

private

Definition at line 77 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ maxChi2InHist_

double SiPixelGainCalibrationAnalysis::maxChi2InHist_

private

Definition at line 92 of file SiPixelGainCalibrationAnalysis.h.

Referenced by newDetID().

◆ maxGainInHist_

double SiPixelGainCalibrationAnalysis::maxGainInHist_

private

Definition at line 91 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and newDetID().

◆ min_nentries_

uint16_t SiPixelGainCalibrationAnalysis::min_nentries_

private

Definition at line 112 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

◆ nfitparameters_

int SiPixelGainCalibrationAnalysis::nfitparameters_

private

Definition at line 72 of file SiPixelGainCalibrationAnalysis.h.

Referenced by SiPixelGainCalibrationAnalysis().

◆ pedhi_

float SiPixelGainCalibrationAnalysis::pedhi_

private

Definition at line 111 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ pedlow_

float SiPixelGainCalibrationAnalysis::pedlow_

private

Definition at line 110 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ plateau_max_slope_

double SiPixelGainCalibrationAnalysis::plateau_max_slope_

private

Definition at line 84 of file SiPixelGainCalibrationAnalysis.h.

◆ recordName_

std::string SiPixelGainCalibrationAnalysis::recordName_

private

Definition at line 100 of file SiPixelGainCalibrationAnalysis.h.

◆ reject_badpoints_

bool SiPixelGainCalibrationAnalysis::reject_badpoints_

private

Definition at line 81 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ reject_badpoints_frac_

double SiPixelGainCalibrationAnalysis::reject_badpoints_frac_

private

Definition at line 86 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ reject_first_point_

bool SiPixelGainCalibrationAnalysis::reject_first_point_

private

Definition at line 85 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ reject_plateaupoints_

bool SiPixelGainCalibrationAnalysis::reject_plateaupoints_

private

Definition at line 82 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ reject_single_entries_

bool SiPixelGainCalibrationAnalysis::reject_single_entries_

private

Definition at line 83 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

◆ saveALLHistograms_

bool SiPixelGainCalibrationAnalysis::saveALLHistograms_

private

Definition at line 93 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ savePixelHists_

bool SiPixelGainCalibrationAnalysis::savePixelHists_

private

Definition at line 88 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ scalarVcalHigh_VcalLow_

double SiPixelGainCalibrationAnalysis::scalarVcalHigh_VcalLow_

private

Definition at line 114 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), newDetID(), and SiPixelGainCalibrationAnalysis().

◆ statusNumbers_

int* SiPixelGainCalibrationAnalysis::statusNumbers_

private

Definition at line 119 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), printSummary(), and SiPixelGainCalibrationAnalysis().

◆ sum_gain_cols_

bool SiPixelGainCalibrationAnalysis::sum_gain_cols_

private

Definition at line 95 of file SiPixelGainCalibrationAnalysis.h.

◆ sum_ped_cols_

bool SiPixelGainCalibrationAnalysis::sum_ped_cols_

private

Definition at line 94 of file SiPixelGainCalibrationAnalysis.h.

◆ summary_

std::ofstream SiPixelGainCalibrationAnalysis::summary_

private

Definition at line 117 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), printSummary(), and SiPixelGainCalibrationAnalysis().

◆ useVcalHigh_

bool SiPixelGainCalibrationAnalysis::useVcalHigh_

private

Definition at line 113 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

◆ writeSummary_

bool SiPixelGainCalibrationAnalysis::writeSummary_

private

Definition at line 97 of file SiPixelGainCalibrationAnalysis.h.

Referenced by calibrationEnd(), and doFits().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SiPixelGainCalibrationAnalysis()

◆ ~SiPixelGainCalibrationAnalysis()

Member Function Documentation

◆ CalculateAveragePerColumn()

◆ calibrationEnd()

◆ calibrationSetup()

◆ checkCorrectCalibrationType()

◆ doFits()

◆ doSetup()

◆ fillDatabase()

◆ newDetID()

◆ printSummary()

Member Data Documentation

◆ appendMode_

◆ bookBIGCalibPayload_

◆ bookkeeper_

◆ bookkeeper_pixels_

◆ chi2ProbThreshold_

◆ chi2Threshold_

◆ conf_

◆ currentDetID_

◆ filldb_

◆ fitfunction_

◆ func_

◆ gainhi_

◆ gainlow_

◆ graph_

◆ ignoreMode_

◆ listofdetids_

◆ maxChi2InHist_

◆ maxGainInHist_

◆ min_nentries_

◆ nfitparameters_

◆ pedhi_

◆ pedlow_

◆ plateau_max_slope_

◆ recordName_

◆ reject_badpoints_

◆ reject_badpoints_frac_

◆ reject_first_point_

◆ reject_plateaupoints_

◆ reject_single_entries_

◆ saveALLHistograms_

◆ savePixelHists_

◆ scalarVcalHigh_VcalLow_

◆ statusNumbers_

◆ sum_gain_cols_

◆ sum_ped_cols_

◆ summary_

◆ useVcalHigh_

◆ writeSummary_