|
|
|
#include <CalibTracker/SiPixelGainCalibrationAnalysis/interface/SiPixelGainCalibrationAnalysis.h>
Description: <one line="" class="" summary>="">
Implementation: <Notes on="" implementation>="">
Definition at line 49 of file SiPixelGainCalibrationAnalysis.h.
| SiPixelGainCalibrationAnalysis::SiPixelGainCalibrationAnalysis | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 35 of file SiPixelGainCalibrationAnalysis.cc.
References currentDetID_, fitfunction_, func_, graph_, 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 87 of file SiPixelGainCalibrationAnalysis.cc.
{
}
| std::vector< float > SiPixelGainCalibrationAnalysis::CalculateAveragePerColumn | ( | uint32_t | detid, |
| std::string | label | ||
| ) | [private] |
Definition at line 94 of file SiPixelGainCalibrationAnalysis.cc.
References bookkeeper_, cond::rpcobgas::detid, label, and query::result.
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 | ( | ) | [private, virtual] |
Reimplemented from SiPixelOfflineCalibAnalysisBase.
Definition at line 174 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 | ) | [private, virtual] |
Reimplemented from SiPixelOfflineCalibAnalysisBase.
Definition at line 128 of file SiPixelGainCalibrationAnalysis.cc.
{
}
| bool SiPixelGainCalibrationAnalysis::checkCorrectCalibrationType | ( | ) | [virtual] |
Reimplemented from SiPixelOfflineCalibAnalysisBase.
Definition at line 112 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 191 of file SiPixelGainCalibrationAnalysis.cc.
References abs, SiPixelOfflineCalibAnalysisBase::bookDQMHistogram1D(), bookkeeper_, bookkeeper_pixels_, ExpressReco_HICollisions_FallBack::chi2, chi2ProbThreshold_, chi2Threshold_, currentDetID_, cond::rpcobgas::detid, spr::find(), func_, gainhi_, gainlow_, graph_, edm::detail::isnan(), listofdetids_, maxGainInHist_, min_nentries_, 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(), ntuplemaker::status, statusNumbers_, 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 183 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 | ) | [private, virtual] |
Reimplemented from SiPixelOfflineCalibAnalysisBase.
Definition at line 472 of file SiPixelGainCalibrationAnalysis.cc.
References SiPixelOfflineCalibAnalysisBase::bookDQMHistogram1D(), SiPixelOfflineCalibAnalysisBase::bookDQMHistoPlaquetteSummary2D(), bookkeeper_, cond::rpcobgas::detid, maxChi2InHist_, maxGainInHist_, scalarVcalHigh_VcalLow_, SiPixelOfflineCalibAnalysisBase::setDQMDirectory(), 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 133 of file SiPixelGainCalibrationAnalysis.cc.
References begin, bookkeeper_, bookkeeper_pixels_, CalculateAveragePerColumn(), chi2Threshold_, cond::rpcobgas::detid, i, Association::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();
}
}
bool SiPixelGainCalibrationAnalysis::appendMode_ [private] |
Definition at line 104 of file SiPixelGainCalibrationAnalysis.h.
bool SiPixelGainCalibrationAnalysis::bookBIGCalibPayload_ [private] |
Definition at line 90 of file SiPixelGainCalibrationAnalysis.h.
std::map<uint32_t,std::map<std::string, MonitorElement *> > SiPixelGainCalibrationAnalysis::bookkeeper_ [private] |
Definition at line 71 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 72 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), and printSummary().
double SiPixelGainCalibrationAnalysis::chi2ProbThreshold_ [private] |
Definition at line 93 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
double SiPixelGainCalibrationAnalysis::chi2Threshold_ [private] |
Definition at line 92 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), and printSummary().
Definition at line 69 of file SiPixelGainCalibrationAnalysis.h.
uint32_t SiPixelGainCalibrationAnalysis::currentDetID_ [private] |
Definition at line 121 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), and SiPixelGainCalibrationAnalysis().
bool SiPixelGainCalibrationAnalysis::filldb_ [private] |
Definition at line 99 of file SiPixelGainCalibrationAnalysis.h.
Referenced by calibrationEnd().
std::string SiPixelGainCalibrationAnalysis::fitfunction_ [private] |
Definition at line 76 of file SiPixelGainCalibrationAnalysis.h.
Referenced by SiPixelGainCalibrationAnalysis().
TF1* SiPixelGainCalibrationAnalysis::func_ [private] |
Definition at line 77 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), and SiPixelGainCalibrationAnalysis().
float SiPixelGainCalibrationAnalysis::gainhi_ [private] |
Definition at line 112 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
float SiPixelGainCalibrationAnalysis::gainlow_ [private] |
Definition at line 111 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
TGraphErrors* SiPixelGainCalibrationAnalysis::graph_ [private] |
Definition at line 78 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), and SiPixelGainCalibrationAnalysis().
bool SiPixelGainCalibrationAnalysis::ignoreMode_ [private] |
Definition at line 81 of file SiPixelGainCalibrationAnalysis.h.
Referenced by checkCorrectCalibrationType().
std::vector<uint32_t> SiPixelGainCalibrationAnalysis::listofdetids_ [private] |
Definition at line 80 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
double SiPixelGainCalibrationAnalysis::maxChi2InHist_ [private] |
Definition at line 95 of file SiPixelGainCalibrationAnalysis.h.
Referenced by newDetID().
double SiPixelGainCalibrationAnalysis::maxGainInHist_ [private] |
Definition at line 94 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), and newDetID().
uint16_t SiPixelGainCalibrationAnalysis::min_nentries_ [private] |
Definition at line 115 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), and SiPixelGainCalibrationAnalysis().
int SiPixelGainCalibrationAnalysis::nfitparameters_ [private] |
Definition at line 75 of file SiPixelGainCalibrationAnalysis.h.
Referenced by SiPixelGainCalibrationAnalysis().
float SiPixelGainCalibrationAnalysis::pedhi_ [private] |
Definition at line 114 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
float SiPixelGainCalibrationAnalysis::pedlow_ [private] |
Definition at line 113 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
double SiPixelGainCalibrationAnalysis::plateau_max_slope_ [private] |
Definition at line 87 of file SiPixelGainCalibrationAnalysis.h.
std::string SiPixelGainCalibrationAnalysis::recordName_ [private] |
Definition at line 103 of file SiPixelGainCalibrationAnalysis.h.
bool SiPixelGainCalibrationAnalysis::reject_badpoints_ [private] |
Definition at line 84 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
double SiPixelGainCalibrationAnalysis::reject_badpoints_frac_ [private] |
Definition at line 89 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
bool SiPixelGainCalibrationAnalysis::reject_first_point_ [private] |
Definition at line 88 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
bool SiPixelGainCalibrationAnalysis::reject_plateaupoints_ [private] |
Definition at line 85 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
bool SiPixelGainCalibrationAnalysis::reject_single_entries_ [private] |
Definition at line 86 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), and SiPixelGainCalibrationAnalysis().
bool SiPixelGainCalibrationAnalysis::saveALLHistograms_ [private] |
Definition at line 96 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
bool SiPixelGainCalibrationAnalysis::savePixelHists_ [private] |
Definition at line 91 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
double SiPixelGainCalibrationAnalysis::scalarVcalHigh_VcalLow_ [private] |
Definition at line 117 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), newDetID(), and SiPixelGainCalibrationAnalysis().
int* SiPixelGainCalibrationAnalysis::statusNumbers_ [private] |
Definition at line 122 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), printSummary(), and SiPixelGainCalibrationAnalysis().
bool SiPixelGainCalibrationAnalysis::sum_gain_cols_ [private] |
Definition at line 98 of file SiPixelGainCalibrationAnalysis.h.
bool SiPixelGainCalibrationAnalysis::sum_ped_cols_ [private] |
Definition at line 97 of file SiPixelGainCalibrationAnalysis.h.
ofstream SiPixelGainCalibrationAnalysis::summary_ [private] |
Definition at line 120 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits(), printSummary(), and SiPixelGainCalibrationAnalysis().
bool SiPixelGainCalibrationAnalysis::useVcalHigh_ [private] |
Definition at line 116 of file SiPixelGainCalibrationAnalysis.h.
Referenced by doFits().
bool SiPixelGainCalibrationAnalysis::writeSummary_ [private] |
Definition at line 100 of file SiPixelGainCalibrationAnalysis.h.
Referenced by calibrationEnd(), and doFits().
1.7.3