|
|
|
#include <CalibrationSummaryFactory.h>
Protected Member Functions | |
| void | extract (Iterator) |
| void | format () |
Definition at line 6 of file CalibrationSummaryFactory.h.
| void CalibrationSummaryFactory::extract | ( | Iterator | iter | ) | [protected, virtual] |
Reimplemented from SummaryPlotFactory< CommissioningAnalysis * >.
Definition at line 12 of file CalibrationSummaryFactory.cc.
References CalibrationAnalysis::amplitude(), CalibrationAnalysis::amplitudeMax(), CalibrationAnalysis::amplitudeMean(), CalibrationAnalysis::amplitudeMin(), CalibrationAnalysis::amplitudeSpread(), combineCards::bins, sistrip::CALIBRATION_AMPLITUDE, sistrip::CALIBRATION_AMPLITUDE_ALLSTRIPS, sistrip::CALIBRATION_AMPLITUDE_MAX, sistrip::CALIBRATION_AMPLITUDE_MIN, sistrip::CALIBRATION_CHI2, sistrip::CALIBRATION_CHI2_ALLSTRIPS, sistrip::CALIBRATION_CHI2_MAX, sistrip::CALIBRATION_CHI2_MIN, sistrip::CALIBRATION_RISETIME, sistrip::CALIBRATION_RISETIME_ALLSTRIPS, sistrip::CALIBRATION_RISETIME_MAX, sistrip::CALIBRATION_RISETIME_MIN, sistrip::CALIBRATION_SMEARING, sistrip::CALIBRATION_SMEARING_ALLSTRIPS, sistrip::CALIBRATION_SMEARING_MAX, sistrip::CALIBRATION_SMEARING_MIN, sistrip::CALIBRATION_TAIL, sistrip::CALIBRATION_TAIL_ALLSTRIPS, sistrip::CALIBRATION_TAIL_MAX, sistrip::CALIBRATION_TAIL_MIN, sistrip::CALIBRATION_TIMECONSTANT, sistrip::CALIBRATION_TIMECONSTANT_ALLSTRIPS, sistrip::CALIBRATION_TIMECONSTANT_MAX, sistrip::CALIBRATION_TIMECONSTANT_MIN, SiStripFecKey::ccuAddr(), SiStripFecKey::ccuChan(), CalibrationAnalysis::chi2(), CalibrationAnalysis::chi2Max(), CalibrationAnalysis::chi2Mean(), CalibrationAnalysis::chi2Min(), CalibrationAnalysis::chi2Spread(), error, SiStripFecKey::fecCrate(), SiStripFecKey::fecRing(), SiStripFecKey::fecSlot(), SummaryGenerator::fillMap(), SummaryPlotFactoryBase::generator_, SummaryPlotFactoryBase::gran_, i, SiStripFecKey::i2cAddr(), sistrip::invalid_, SiStripKey::key(), SummaryPlotFactoryBase::level_, SiStripFecKey::lldChan(), sistrip::mlSummaryPlots_, SummaryPlotFactoryBase::mon_, SiStripEnumsAndStrings::monitorable(), CalibrationAnalysis::riseTime(), CalibrationAnalysis::riseTimeMax(), CalibrationAnalysis::riseTimeMean(), CalibrationAnalysis::riseTimeMin(), CalibrationAnalysis::riseTimeSpread(), CalibrationAnalysis::smearing(), CalibrationAnalysis::smearingMax(), CalibrationAnalysis::smearingMean(), CalibrationAnalysis::smearingMin(), CalibrationAnalysis::smearingSpread(), CalibrationAnalysis::tail(), CalibrationAnalysis::tailMax(), CalibrationAnalysis::tailMean(), CalibrationAnalysis::tailMin(), CalibrationAnalysis::tailSpread(), groupFilesInBlocks::temp, CalibrationAnalysis::timeConstant(), CalibrationAnalysis::timeConstantMax(), CalibrationAnalysis::timeConstantMean(), CalibrationAnalysis::timeConstantMin(), CalibrationAnalysis::timeConstantSpread(), and relativeConstraints::value.
{
CalibrationAnalysis* anal = dynamic_cast<CalibrationAnalysis*>( iter->second );
if ( !anal ) { return; }
std::vector<float> temp(128, 1. * sistrip::invalid_ );
std::vector< std::vector<float> > value( 2, temp );
std::vector< std::vector<float> > error( 2, temp );
std::vector< std::vector<float> > amplitude( 2, temp );
std::vector< std::vector<float> > tail( 2, temp );
std::vector< std::vector<float> > riseTime( 2, temp );
std::vector< std::vector<float> > timeConstant( 2, temp );
std::vector< std::vector<float> > smearing( 2, temp );
std::vector< std::vector<float> > chi2( 2, temp );
amplitude[0] = anal->amplitude()[0];
amplitude[1] = anal->amplitude()[1];
tail[0] = anal->tail()[0];
tail[1] = anal->tail()[1];
riseTime[0] = anal->riseTime()[0];
riseTime[1] = anal->riseTime()[1];
timeConstant[0] = anal->timeConstant()[0];
timeConstant[1] = anal->timeConstant()[1];
smearing[0] = anal->smearing()[0];
smearing[1] = anal->smearing()[1];
chi2[0] = anal->chi2()[0];
chi2[1] = anal->chi2()[1];
SiStripFecKey lldKey = SiStripFecKey(iter->first);
uint32_t key1 = SiStripFecKey(lldKey.fecCrate(),
lldKey.fecSlot(),
lldKey.fecRing(),
lldKey.ccuAddr(),
lldKey.ccuChan(),
lldKey.lldChan(),
lldKey.i2cAddr(lldKey.lldChan(),true)).key();
uint32_t key2 = SiStripFecKey(lldKey.fecCrate(),
lldKey.fecSlot(),
lldKey.fecRing(),
lldKey.ccuAddr(),
lldKey.ccuChan(),
lldKey.lldChan(),
lldKey.i2cAddr(lldKey.lldChan(),false)).key();
bool all_strips = false;
bool with_error = false;
if ( mon_ == sistrip::CALIBRATION_AMPLITUDE_ALLSTRIPS) {
all_strips = true;
uint16_t bins = amplitude[amplitude[0].size() < amplitude[1].size() ? 1 : 0].size();
for ( uint16_t i = 0; i < bins; i++ ) {
value[0][i] = amplitude[0][i]/10.;
value[1][i] = amplitude[1][i]/10.;
}
} else if ( mon_ == sistrip::CALIBRATION_AMPLITUDE) {
with_error = true;
value[0][0] = anal->amplitudeMean()[0]/10.;
value[1][0] = anal->amplitudeMean()[1]/10.;
error[0][0] = anal->amplitudeSpread()[0]/10.;
error[1][0] = anal->amplitudeSpread()[1]/10.;
} else if ( mon_ == sistrip::CALIBRATION_AMPLITUDE_MIN) {
value[0][0] = anal->amplitudeMin()[0]/10.;
value[1][0] = anal->amplitudeMin()[1]/10.;
} else if ( mon_ == sistrip::CALIBRATION_AMPLITUDE_MAX) {
value[0][0] = anal->amplitudeMax()[0]/10.;
value[1][0] = anal->amplitudeMax()[1]/10.;
} else if ( mon_ == sistrip::CALIBRATION_TAIL_ALLSTRIPS) {
all_strips = true;
uint16_t bins = tail[tail[0].size() < tail[1].size() ? 1 : 0].size();
for ( uint16_t i = 0; i < bins; i++ ) {
value[0][i] = tail[0][i];
value[1][i] = tail[1][i];
}
} else if ( mon_ == sistrip::CALIBRATION_TAIL) {
with_error = true;
value[0][0] = anal->tailMean()[0];
value[1][0] = anal->tailMean()[1];
error[0][0] = anal->tailSpread()[0];
error[1][0] = anal->tailSpread()[1];
} else if ( mon_ == sistrip::CALIBRATION_TAIL_MIN) {
value[0][0] = anal->tailMin()[0];
value[1][0] = anal->tailMin()[1];
} else if ( mon_ == sistrip::CALIBRATION_TAIL_MAX) {
value[0][0] = anal->tailMax()[0];
value[1][0] = anal->tailMax()[1];
} else if ( mon_ == sistrip::CALIBRATION_RISETIME_ALLSTRIPS) {
all_strips = true;
uint16_t bins = riseTime[riseTime[0].size() < riseTime[1].size() ? 1 : 0].size();
for ( uint16_t i = 0; i < bins; i++ ) {
value[0][i] = riseTime[0][i];
value[1][i] = riseTime[1][i];
}
} else if ( mon_ == sistrip::CALIBRATION_RISETIME) {
with_error = true;
value[0][0] = anal->riseTimeMean()[0];
value[1][0] = anal->riseTimeMean()[1];
error[0][0] = anal->riseTimeSpread()[0];
error[1][0] = anal->riseTimeSpread()[1];
} else if ( mon_ == sistrip::CALIBRATION_RISETIME_MIN) {
value[0][0] = anal->riseTimeMin()[0];
value[1][0] = anal->riseTimeMin()[1];
} else if ( mon_ == sistrip::CALIBRATION_RISETIME_MAX) {
value[0][0] = anal->riseTimeMax()[0];
value[1][0] = anal->riseTimeMax()[1];
} else if ( mon_ == sistrip::CALIBRATION_TIMECONSTANT_ALLSTRIPS) {
all_strips = true;
uint16_t bins = timeConstant[timeConstant[0].size() < timeConstant[1].size() ? 1 : 0].size();
for ( uint16_t i = 0; i < bins; i++ ) {
value[0][i] = timeConstant[0][i];
value[1][i] = timeConstant[1][i];
}
} else if ( mon_ == sistrip::CALIBRATION_TIMECONSTANT) {
with_error = true;
value[0][0] = anal->timeConstantMean()[0];
value[1][0] = anal->timeConstantMean()[1];
error[0][0] = anal->timeConstantSpread()[0];
error[1][0] = anal->timeConstantSpread()[1];
} else if ( mon_ == sistrip::CALIBRATION_TIMECONSTANT_MIN) {
value[0][0] = anal->timeConstantMin()[0];
value[1][0] = anal->timeConstantMin()[1];
} else if ( mon_ == sistrip::CALIBRATION_TIMECONSTANT_MAX) {
value[0][0] = anal->timeConstantMax()[0];
value[1][0] = anal->timeConstantMax()[1];
} else if ( mon_ == sistrip::CALIBRATION_SMEARING_ALLSTRIPS) {
all_strips = true;
uint16_t bins = smearing[smearing[0].size() < smearing[1].size() ? 1 : 0].size();
for ( uint16_t i = 0; i < bins; i++ ) {
value[0][i] = smearing[0][i];
value[1][i] = smearing[1][i];
}
} else if ( mon_ == sistrip::CALIBRATION_SMEARING) {
with_error = true;
value[0][0] = anal->smearingMean()[0];
value[1][0] = anal->smearingMean()[1];
error[0][0] = anal->smearingSpread()[0];
error[1][0] = anal->smearingSpread()[1];
} else if ( mon_ == sistrip::CALIBRATION_SMEARING_MIN) {
value[0][0] = anal->smearingMin()[0];
value[1][0] = anal->smearingMin()[1];
} else if ( mon_ == sistrip::CALIBRATION_SMEARING_MAX) {
value[0][0] = anal->smearingMax()[0];
value[1][0] = anal->smearingMax()[1];
} else if ( mon_ == sistrip::CALIBRATION_CHI2_ALLSTRIPS) {
all_strips = true;
uint16_t bins = chi2[chi2[0].size() < chi2[1].size() ? 1 : 0].size();
for ( uint16_t i = 0; i < bins; i++ ) {
value[0][i] = chi2[0][i];
value[1][i] = chi2[1][i];
}
} else if ( mon_ == sistrip::CALIBRATION_CHI2) {
with_error = true;
value[0][0] = anal->chi2Mean()[0]/100.;
value[1][0] = anal->chi2Mean()[1]/100.;
error[0][0] = anal->chi2Spread()[0]/100.;
error[1][0] = anal->chi2Spread()[1]/100.;
} else if ( mon_ == sistrip::CALIBRATION_CHI2_MIN) {
value[0][0] = anal->chi2Min()[0]/100.;
value[1][0] = anal->chi2Min()[1]/100.;
} else if ( mon_ == sistrip::CALIBRATION_CHI2_MAX) {
value[0][0] = anal->chi2Max()[0]/100.;
value[1][0] = anal->chi2Max()[1]/100.;
} else {
edm::LogWarning(mlSummaryPlots_)
<< "[SummaryPlotFactory::" << __func__ << "]"
<< " Unexpected monitorable: "
<< SiStripEnumsAndStrings::monitorable( SummaryPlotFactoryBase::mon_ );
return;
}
if ( !all_strips ) {
if( !with_error) {
SummaryPlotFactoryBase::generator_->fillMap( SummaryPlotFactoryBase::level_,
SummaryPlotFactoryBase::gran_,
key1,
value[0][0] );
SummaryPlotFactoryBase::generator_->fillMap( SummaryPlotFactoryBase::level_,
SummaryPlotFactoryBase::gran_,
key2,
value[1][0] );
} else {
SummaryPlotFactoryBase::generator_->fillMap( SummaryPlotFactoryBase::level_,
SummaryPlotFactoryBase::gran_,
key1,
value[0][0],
error[0][0]);
SummaryPlotFactoryBase::generator_->fillMap( SummaryPlotFactoryBase::level_,
SummaryPlotFactoryBase::gran_,
key2,
value[1][0],
error[1][0]);
}
} else {
for ( uint16_t istr = 0; istr < value[0].size(); istr++ ) {
SummaryPlotFactoryBase::generator_->fillMap( SummaryPlotFactoryBase::level_,
SummaryPlotFactoryBase::gran_,
key1,
value[0][istr] );
}
for ( uint16_t istr = 0; istr < value[1].size(); istr++ ) {
SummaryPlotFactoryBase::generator_->fillMap( SummaryPlotFactoryBase::level_,
SummaryPlotFactoryBase::gran_,
key2,
value[1][istr] );
}
}
}
| void CalibrationSummaryFactory::format | ( | ) | [protected, virtual] |
Reimplemented from SummaryPlotFactory< CommissioningAnalysis * >.
Definition at line 223 of file CalibrationSummaryFactory.cc.
References sistrip::CALIBRATION_AMPLITUDE, sistrip::CALIBRATION_CHI2, sistrip::CALIBRATION_RISETIME, sistrip::CALIBRATION_SMEARING, sistrip::CALIBRATION_TAIL, sistrip::CALIBRATION_TIMECONSTANT, sistrip::mlSummaryPlots_, SummaryPlotFactoryBase::mon_, and SiStripEnumsAndStrings::monitorable().
{
// Histogram formatting
if ( mon_ == sistrip::CALIBRATION_AMPLITUDE ) {
generator_->axisLabel( "Amplitude (ADC*Nevt/10.)" );
} else if ( mon_ == sistrip::CALIBRATION_TAIL ) {
generator_->axisLabel( "Tail (%)" );
} else if ( mon_ == sistrip::CALIBRATION_RISETIME ) {
generator_->axisLabel( "Rise time (ns)" );
} else if ( mon_ == sistrip::CALIBRATION_TIMECONSTANT ) {
generator_->axisLabel( "Time constant (ns)" );
} else if ( mon_ == sistrip::CALIBRATION_SMEARING ) {
generator_->axisLabel( "Smearing (ns)" );
} else if ( mon_ == sistrip::CALIBRATION_CHI2 ) {
generator_->axisLabel( "Chi2/100." );
} else {
edm::LogWarning(mlSummaryPlots_)
<< "[SummaryPlotFactory::" << __func__ << "]"
<< " Unexpected SummaryHisto value:"
<< SiStripEnumsAndStrings::monitorable( SummaryPlotFactoryBase::mon_ ) ;
}
}
1.7.3