#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 (const edm::ParameterSet &)

	SiPixelOfflineCalibAnalysisBase ()

std::string	translateDetIdToString (uint32_t detid)

	~SiPixelOfflineCalibAnalysisBase () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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::one::EDAnalyzerBase
typedef EDAnalyzerBase	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::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

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...

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

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

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

void	resetItemsToGetFrom (BranchType iType)

Protected Attributes inherited from SiPixelOfflineCalibAnalysisBase
edm::ESGetToken < SiPixelFedCablingMap, SiPixelFedCablingMapRcd >	cablingMapToken_

edm::ESHandle < SiPixelCalibConfiguration >	calib_

std::string	calibrationMode_

edm::ESGetToken < SiPixelCalibConfiguration, SiPixelCalibConfigurationRcd >	calibToken_

edm::ESGetToken < SiPixelCalibConfiguration, SiPixelCalibConfigurationRcd >	calibTokenBeginRun_

edm::ESHandle< TrackerGeometry >	geom_

short	nTriggers_

edm::ESHandle < SiPixelFedCablingMap >	theCablingMap_

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord >	trackerGeomToken_

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 46 of file SiPixelGainCalibrationAnalysis.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 32 of file SiPixelGainCalibrationAnalysis.cc.

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

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

SiPixelGainCalibrationAnalysis::~SiPixelGainCalibrationAnalysis ( )

override

Definition at line 86 of file SiPixelGainCalibrationAnalysis.cc.

86 {}

Member Function Documentation

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

private

Definition at line 91 of file SiPixelGainCalibrationAnalysis.cc.

References bookkeeper_, label, mps_fire::result, and hgcalPerformanceValidation::val.

Referenced by printSummary().

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

void SiPixelGainCalibrationAnalysis::calibrationEnd ( )

overrideprivatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 165 of file SiPixelGainCalibrationAnalysis.cc.

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

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

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

overrideprivatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 123 of file SiPixelGainCalibrationAnalysis.cc.

123 {}

bool SiPixelGainCalibrationAnalysis::checkCorrectCalibrationType ( )

overridevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 108 of file SiPixelGainCalibrationAnalysis.cc.

References SiPixelOfflineCalibAnalysisBase::calibrationMode_, and ignoreMode_.

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

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

overridevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 183 of file SiPixelGainCalibrationAnalysis.cc.

References funct::abs(), Abs(), SiPixelOfflineCalibAnalysisBase::bookDQMHistogram1D(), bookkeeper_, bookkeeper_pixels_, HLT_FULL_cff::chi2, chi2ProbThreshold_, chi2Threshold_, currentDetID_, spr::find(), func_, gainhi_, gainlow_, graph_, cuy::ii, CommonMethods::isnan(), listofdetids_, maxGainInHist_, min_nentries_, npoints, pedhi_, pedlow_, funct::pow(), 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_.

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

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

void SiPixelGainCalibrationAnalysis::fillDatabase ( )

private

Definition at line 174 of file SiPixelGainCalibrationAnalysis.cc.

Referenced by calibrationEnd().

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

void SiPixelGainCalibrationAnalysis::newDetID ( uint32_t detid )

overrideprivatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 467 of file SiPixelGainCalibrationAnalysis.cc.

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

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

void SiPixelGainCalibrationAnalysis::printSummary ( )

private

Definition at line 125 of file SiPixelGainCalibrationAnalysis.cc.

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

Referenced by calibrationEnd().

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

Member Data Documentation

bool SiPixelGainCalibrationAnalysis::appendMode_

private

Definition at line 100 of file SiPixelGainCalibrationAnalysis.h.

bool SiPixelGainCalibrationAnalysis::bookBIGCalibPayload_

private

Definition at line 86 of file SiPixelGainCalibrationAnalysis.h.

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

private

Definition at line 67 of file SiPixelGainCalibrationAnalysis.h.

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

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

private

Definition at line 68 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and printSummary().

double SiPixelGainCalibrationAnalysis::chi2ProbThreshold_

private

Definition at line 89 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::chi2Threshold_

private

Definition at line 88 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and printSummary().

edm::ParameterSet SiPixelGainCalibrationAnalysis::conf_

private

Definition at line 65 of file SiPixelGainCalibrationAnalysis.h.

uint32_t SiPixelGainCalibrationAnalysis::currentDetID_

private

Definition at line 117 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

bool SiPixelGainCalibrationAnalysis::filldb_

private

Definition at line 95 of file SiPixelGainCalibrationAnalysis.h.

Referenced by calibrationEnd().

std::string SiPixelGainCalibrationAnalysis::fitfunction_

private

Definition at line 72 of file SiPixelGainCalibrationAnalysis.h.

Referenced by SiPixelGainCalibrationAnalysis().

TF1* SiPixelGainCalibrationAnalysis::func_

private

Definition at line 73 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

float SiPixelGainCalibrationAnalysis::gainhi_

private

Definition at line 108 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

float SiPixelGainCalibrationAnalysis::gainlow_

private

Definition at line 107 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

TGraphErrors* SiPixelGainCalibrationAnalysis::graph_

private

Definition at line 74 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

bool SiPixelGainCalibrationAnalysis::ignoreMode_

private

Definition at line 77 of file SiPixelGainCalibrationAnalysis.h.

Referenced by checkCorrectCalibrationType().

std::vector<uint32_t> SiPixelGainCalibrationAnalysis::listofdetids_

private

Definition at line 76 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::maxChi2InHist_

private

Definition at line 91 of file SiPixelGainCalibrationAnalysis.h.

Referenced by newDetID().

double SiPixelGainCalibrationAnalysis::maxGainInHist_

private

Definition at line 90 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and newDetID().

uint16_t SiPixelGainCalibrationAnalysis::min_nentries_

private

Definition at line 111 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

int SiPixelGainCalibrationAnalysis::nfitparameters_

private

Definition at line 71 of file SiPixelGainCalibrationAnalysis.h.

Referenced by SiPixelGainCalibrationAnalysis().

float SiPixelGainCalibrationAnalysis::pedhi_

private

Definition at line 110 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

float SiPixelGainCalibrationAnalysis::pedlow_

private

Definition at line 109 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::plateau_max_slope_

private

Definition at line 83 of file SiPixelGainCalibrationAnalysis.h.

std::string SiPixelGainCalibrationAnalysis::recordName_

private

Definition at line 99 of file SiPixelGainCalibrationAnalysis.h.

bool SiPixelGainCalibrationAnalysis::reject_badpoints_

private

Definition at line 80 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::reject_badpoints_frac_

private

Definition at line 85 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::reject_first_point_

private

Definition at line 84 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::reject_plateaupoints_

private

Definition at line 81 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::reject_single_entries_

private

Definition at line 82 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

bool SiPixelGainCalibrationAnalysis::saveALLHistograms_

private

Definition at line 92 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::savePixelHists_

private

Definition at line 87 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::scalarVcalHigh_VcalLow_

private

Definition at line 113 of file SiPixelGainCalibrationAnalysis.h.

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

int* SiPixelGainCalibrationAnalysis::statusNumbers_

private

Definition at line 118 of file SiPixelGainCalibrationAnalysis.h.

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

bool SiPixelGainCalibrationAnalysis::sum_gain_cols_

private

Definition at line 94 of file SiPixelGainCalibrationAnalysis.h.

bool SiPixelGainCalibrationAnalysis::sum_ped_cols_

private

Definition at line 93 of file SiPixelGainCalibrationAnalysis.h.

std::ofstream SiPixelGainCalibrationAnalysis::summary_

private

Definition at line 116 of file SiPixelGainCalibrationAnalysis.h.

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

bool SiPixelGainCalibrationAnalysis::useVcalHigh_

private

Definition at line 112 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::writeSummary_

private

Definition at line 96 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

Member Function Documentation

Member Data Documentation