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

Inheritance diagram for SiPixelGainCalibrationAnalysis:

Public Member Functions
virtual bool	checkCorrectCalibrationType ()

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

void	doSetup (const edm::ParameterSet &)

	SiPixelGainCalibrationAnalysis (const edm::ParameterSet &iConfig)

	~SiPixelGainCalibrationAnalysis ()

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 ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 &	itemsToGetFromEvent () 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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

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

virtual	~EDConsumerBase () noexcept(false)

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

virtual void	calibrationEnd ()

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

void	fillDatabase ()

virtual void	newDetID (uint32_t detid)

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 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 Public Attributes inherited from SiPixelOfflineCalibAnalysisBase
static TF1 *	fitFunction_ = NULL

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)

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

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

Constructor & Destructor Documentation

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

explicit

Definition at line 34 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 ( )

Definition at line 86 of file SiPixelGainCalibrationAnalysis.cc.

87 {

88 }

Member Function Documentation

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

private

Definition at line 93 of file SiPixelGainCalibrationAnalysis.cc.

References bookkeeper_, diffTwoXMLs::label, trackingPlots::ncols, mps_fire::result, and heppy_batch::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 ( )

privatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 173 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 )

privatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 127 of file SiPixelGainCalibrationAnalysis.cc.

128 {

129 }

bool SiPixelGainCalibrationAnalysis::checkCorrectCalibrationType ( )

virtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 111 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
	)

virtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 190 of file SiPixelGainCalibrationAnalysis.cc.

References funct::abs(), Abs(), SiPixelOfflineCalibAnalysisBase::bookDQMHistogram1D(), bookkeeper_, bookkeeper_pixels_, HiEvtPlane_cfi::chi2, chi2ProbThreshold_, chi2Threshold_, currentDetID_, spr::find(), objects.autophobj::float, func_, gainhi_, gainlow_, graph_, cuy::ii, edm::detail::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_.

 {
   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=0;
   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=1;
         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_.size()!=0)
     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 182 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 )

privatevirtual

Reimplemented from SiPixelOfflineCalibAnalysisBase.

Definition at line 471 of file SiPixelGainCalibrationAnalysis.cc.

References SiPixelOfflineCalibAnalysisBase::bookDQMHistogram1D(), SiPixelOfflineCalibAnalysisBase::bookDQMHistoPlaquetteSummary2D(), bookkeeper_, DEFINE_FWK_MODULE, 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 132 of file SiPixelGainCalibrationAnalysis.cc.

References begin, bookkeeper_, bookkeeper_pixels_, CalculateAveragePerColumn(), chi2Threshold_, mps_fire::i, cuy::ii, genParticles_cff::map, 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 103 of file SiPixelGainCalibrationAnalysis.h.

bool SiPixelGainCalibrationAnalysis::bookBIGCalibPayload_

private

Definition at line 89 of file SiPixelGainCalibrationAnalysis.h.

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

private

Definition at line 70 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 71 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and printSummary().

double SiPixelGainCalibrationAnalysis::chi2ProbThreshold_

private

Definition at line 92 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::chi2Threshold_

private

Definition at line 91 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and printSummary().

edm::ParameterSet SiPixelGainCalibrationAnalysis::conf_

private

Definition at line 68 of file SiPixelGainCalibrationAnalysis.h.

uint32_t SiPixelGainCalibrationAnalysis::currentDetID_

private

Definition at line 120 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

bool SiPixelGainCalibrationAnalysis::filldb_

private

Definition at line 98 of file SiPixelGainCalibrationAnalysis.h.

Referenced by calibrationEnd().

std::string SiPixelGainCalibrationAnalysis::fitfunction_

private

Definition at line 75 of file SiPixelGainCalibrationAnalysis.h.

Referenced by SiPixelGainCalibrationAnalysis().

TF1* SiPixelGainCalibrationAnalysis::func_

private

Definition at line 76 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

float SiPixelGainCalibrationAnalysis::gainhi_

private

Definition at line 111 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

float SiPixelGainCalibrationAnalysis::gainlow_

private

Definition at line 110 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

TGraphErrors* SiPixelGainCalibrationAnalysis::graph_

private

Definition at line 77 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

bool SiPixelGainCalibrationAnalysis::ignoreMode_

private

Definition at line 80 of file SiPixelGainCalibrationAnalysis.h.

Referenced by checkCorrectCalibrationType().

std::vector<uint32_t> SiPixelGainCalibrationAnalysis::listofdetids_

private

Definition at line 79 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::maxChi2InHist_

private

Definition at line 94 of file SiPixelGainCalibrationAnalysis.h.

Referenced by newDetID().

double SiPixelGainCalibrationAnalysis::maxGainInHist_

private

Definition at line 93 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and newDetID().

uint16_t SiPixelGainCalibrationAnalysis::min_nentries_

private

Definition at line 114 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

int SiPixelGainCalibrationAnalysis::nfitparameters_

private

Definition at line 74 of file SiPixelGainCalibrationAnalysis.h.

Referenced by SiPixelGainCalibrationAnalysis().

float SiPixelGainCalibrationAnalysis::pedhi_

private

Definition at line 113 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

float SiPixelGainCalibrationAnalysis::pedlow_

private

Definition at line 112 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::plateau_max_slope_

private

Definition at line 86 of file SiPixelGainCalibrationAnalysis.h.

std::string SiPixelGainCalibrationAnalysis::recordName_

private

Definition at line 102 of file SiPixelGainCalibrationAnalysis.h.

bool SiPixelGainCalibrationAnalysis::reject_badpoints_

private

Definition at line 83 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::reject_badpoints_frac_

private

Definition at line 88 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::reject_first_point_

private

Definition at line 87 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::reject_plateaupoints_

private

Definition at line 84 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::reject_single_entries_

private

Definition at line 85 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits(), and SiPixelGainCalibrationAnalysis().

bool SiPixelGainCalibrationAnalysis::saveALLHistograms_

private

Definition at line 95 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::savePixelHists_

private

Definition at line 90 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

double SiPixelGainCalibrationAnalysis::scalarVcalHigh_VcalLow_

private

Definition at line 116 of file SiPixelGainCalibrationAnalysis.h.

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

int* SiPixelGainCalibrationAnalysis::statusNumbers_

private

Definition at line 121 of file SiPixelGainCalibrationAnalysis.h.

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

bool SiPixelGainCalibrationAnalysis::sum_gain_cols_

private

Definition at line 97 of file SiPixelGainCalibrationAnalysis.h.

bool SiPixelGainCalibrationAnalysis::sum_ped_cols_

private

Definition at line 96 of file SiPixelGainCalibrationAnalysis.h.

std::ofstream SiPixelGainCalibrationAnalysis::summary_

private

Definition at line 119 of file SiPixelGainCalibrationAnalysis.h.

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

bool SiPixelGainCalibrationAnalysis::useVcalHigh_

private

Definition at line 115 of file SiPixelGainCalibrationAnalysis.h.

Referenced by doFits().

bool SiPixelGainCalibrationAnalysis::writeSummary_

private

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