|
|
|
#include <EBDataCertificationTask.h>
Definition at line 19 of file EBDataCertificationTask.h.
| EBDataCertificationTask::EBDataCertificationTask | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 28 of file EBDataCertificationTask.cc.
References cloneME_, dqmStore_, enableCleanup_, edm::ParameterSet::getUntrackedParameter(), hDAQ_, hDCS_, hDQM_, hFrontendByLumi_, hIntegrityByLumi_, hSynchronizationByLumi_, i, meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, mergeRuns_, cmsCodeRules::cppFunctionSkipper::operator, and prefixME_.
{
dqmStore_ = edm::Service<DQMStore>().operator->();
// cloneME switch
cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
meEBDataCertificationSummary_ = 0;
meEBDataCertificationSummaryMap_ = 0;
for (int i = 0; i < 36; i++) {
meEBDataCertification_[i] = 0;
}
hDQM_ = 0;
hDAQ_ = 0;
hDCS_ = 0;
hIntegrityByLumi_ = 0;
hFrontendByLumi_ = 0;
hSynchronizationByLumi_ = 0;
}
| EBDataCertificationTask::~EBDataCertificationTask | ( | ) | [virtual] |
| void EBDataCertificationTask::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
| void EBDataCertificationTask::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 60 of file EBDataCertificationTask.cc.
References DQMStore::book2D(), DQMStore::bookFloat(), dqmStore_, MonitorElement::Fill(), i, meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, mergeVDriftHistosByStation::name, prefixME_, Numbers::sEB(), MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().
{
std::string name;
if ( dqmStore_ ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
name = "CertificationSummary";
meEBDataCertificationSummary_ = dqmStore_->bookFloat(name);
meEBDataCertificationSummary_->Fill(-1.0);
name = "CertificationSummaryMap";
meEBDataCertificationSummaryMap_ = dqmStore_->book2D(name, name, 72, 0., 72., 34, 0., 34.);
meEBDataCertificationSummaryMap_->setAxisTitle("jphi", 1);
meEBDataCertificationSummaryMap_->setAxisTitle("jeta", 2);
dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/CertificationContents");
for (int i = 0; i < 36; i++) {
name = "EcalBarrel_" + Numbers::sEB(i+1);
meEBDataCertification_[i] = dqmStore_->bookFloat(name);
meEBDataCertification_[i]->Fill(-1.0);
}
}
}
| void EBDataCertificationTask::beginLuminosityBlock | ( | const edm::LuminosityBlock & | lumiBlock, |
| const edm::EventSetup & | iSetup | ||
| ) | [protected, virtual] |
BeginLuminosityBlock.
Reimplemented from edm::EDAnalyzer.
Definition at line 95 of file EBDataCertificationTask.cc.
{
}
| void EBDataCertificationTask::beginRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
BeginRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 241 of file EBDataCertificationTask.cc.
References mergeRuns_, and reset().
{
if ( ! mergeRuns_ ) this->reset();
}
| void EBDataCertificationTask::cleanup | ( | void | ) | [protected] |
Cleanup.
Definition at line 337 of file EBDataCertificationTask.cc.
References cloneME_, dqmStore_, edm::getName(), MonitorElement::getName(), hDAQ_, hDCS_, hDQM_, hFrontendByLumi_, hIntegrityByLumi_, hSynchronizationByLumi_, i, meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().
Referenced by endJob().
{
if ( cloneME_ ) {
if( hDQM_ ) delete hDQM_;
if( hDAQ_ ) delete hDAQ_;
if( hDCS_ ) delete hDCS_;
if( hIntegrityByLumi_ ) delete hIntegrityByLumi_;
if( hFrontendByLumi_ ) delete hFrontendByLumi_;
if( hSynchronizationByLumi_ ) delete hSynchronizationByLumi_;
}
hDQM_ = 0;
hDAQ_ = 0;
hDCS_ = 0;
hIntegrityByLumi_ = 0;
hFrontendByLumi_ = 0;
hSynchronizationByLumi_ = 0;
if ( dqmStore_ ) {
dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
if ( meEBDataCertificationSummary_ ) dqmStore_->removeElement( meEBDataCertificationSummary_->getName() );
if ( meEBDataCertificationSummaryMap_ ) dqmStore_->removeElement( meEBDataCertificationSummaryMap_->getName() );
dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/CertificationContents");
for (int i = 0; i < 36; i++) {
if ( meEBDataCertification_[i] ) dqmStore_->removeElement( meEBDataCertification_[i]->getName() );
}
}
}
| void EBDataCertificationTask::endJob | ( | void | ) | [protected, virtual] |
EndJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 89 of file EBDataCertificationTask.cc.
References cleanup(), and enableCleanup_.
{
if ( enableCleanup_ ) this->cleanup();
}
| void EBDataCertificationTask::endLuminosityBlock | ( | const edm::LuminosityBlock & | lumiBlock, |
| const edm::EventSetup & | iSetup | ||
| ) | [protected, virtual] |
EndLuminosityBlock.
Reimplemented from edm::EDAnalyzer.
Definition at line 99 of file EBDataCertificationTask.cc.
References abs, cloneME_, dqmStore_, benchmark_cfg::errors, MonitorElement::Fill(), DQMStore::get(), hDAQ_, hDCS_, hFrontendByLumi_, hIntegrityByLumi_, hSynchronizationByLumi_, i, ecalpyutils::ism(), meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, min, prefixME_, reset(), Numbers::sEB(), and MonitorElement::setBinContent().
{
this->reset();
MonitorElement* me;
// evaluate the DQM quality of observables checked by lumi
float DQMVal[36];
for (int i = 0; i < 36; i++) {
DQMVal[i] = -1.;
}
me = dqmStore_->get(prefixME_ + "/EBIntegrityTask/EBIT weighted integrity errors by lumi");
hIntegrityByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hIntegrityByLumi_ );
me = dqmStore_->get(prefixME_ + "/EBStatusFlagsTask/FEStatus/EBSFT weighted frontend errors by lumi");
hFrontendByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hFrontendByLumi_ );
me = dqmStore_->get(prefixME_ + "/EBRawDataTask/EBRDT FE synchronization errors by lumi");
hSynchronizationByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hSynchronizationByLumi_ );
if( hIntegrityByLumi_ && hFrontendByLumi_ && hSynchronizationByLumi_ ) {
float integrityErrSum = 0.;
float integrityQual = 1.0;
float frontendErrSum = 0.;
float frontendQual = 1.0;
float synchronizationErrSum = 0.;
float synchronizationQual = 1.0;
for ( int i=0; i<36; i++) {
float ismIntegrityQual = 1.0;
if( hIntegrityByLumi_->GetBinContent(0) > 0 ) {
float errors = hIntegrityByLumi_->GetBinContent(i+1);
ismIntegrityQual = 1.0 - errors/hIntegrityByLumi_->GetBinContent(0);
integrityErrSum += errors;
}
float ismFrontendQual = 1.0;
if( hFrontendByLumi_->GetBinContent(0) > 0 ) {
float errors = hFrontendByLumi_->GetBinContent(i+1);
ismFrontendQual = 1.0 - errors/hFrontendByLumi_->GetBinContent(0);
frontendErrSum += errors;
}
float ismSynchronizationQual = 1.0;
if( hSynchronizationByLumi_->GetBinContent(0) > 0 ) {
float errors = hSynchronizationByLumi_->GetBinContent(i+1);
ismSynchronizationQual = 1.0 - errors/hSynchronizationByLumi_->GetBinContent(0);
synchronizationErrSum += errors;
}
float minVal= std::min(ismIntegrityQual,ismFrontendQual);
DQMVal[i] = std::min(minVal,ismSynchronizationQual);
}
if( hIntegrityByLumi_->GetBinContent(0) > 0 ) integrityQual = 1.0 - integrityErrSum/hIntegrityByLumi_->GetBinContent(0)/36.;
if( hFrontendByLumi_->GetBinContent(0) > 0 ) frontendQual = 1.0 - frontendErrSum/hFrontendByLumi_->GetBinContent(0)/36.;
if( hSynchronizationByLumi_->GetBinContent(0) > 0 ) synchronizationQual = 1.0 - synchronizationErrSum/hSynchronizationByLumi_->GetBinContent(0)/36.;
float minVal = std::min(integrityQual,frontendQual);
float totDQMVal = std::min(minVal,synchronizationQual);
me = dqmStore_->get((prefixME_ + "/EventInfo/reportSummary"));
if( me ) me->Fill(totDQMVal);
for ( int i=0; i<36; i++) {
me = dqmStore_->get(prefixME_ + "/EventInfo/reportSummaryContents/EcalBarrel_" + Numbers::sEB(i+1) ) ;
if( me ) me->Fill(DQMVal[i]);
me = dqmStore_->get(prefixME_ + "/EventInfo/reportSummaryMap");
if( me ) {
for ( int iett = 0; iett < 34; iett++ ) {
for ( int iptt = 0; iptt < 72; iptt++ ) {
int ism = ( iett<17 ) ? iptt/4+1 : 18+iptt/4+1;
if( i == (ism-1) ) me->setBinContent(iptt+1, iett+1, DQMVal[ism-1]);
}
}
}
}
}
// now combine reduced DQM with DCS and DAQ
me = dqmStore_->get(prefixME_ + "/EventInfo/DAQSummaryMap");
hDAQ_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDAQ_ );
me = dqmStore_->get(prefixME_ + "/EventInfo/DCSSummaryMap");
hDCS_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDCS_ );
float sumCert = 0.;
float sumCertEB[36];
int nValidChannels = 0;
int nValidChannelsEB[36];
for (int i = 0; i < 36; i++) {
sumCertEB[i] = 0.;
nValidChannelsEB[i] = 0;
}
for ( int iett = 0; iett < 34; iett++ ) {
for ( int iptt = 0; iptt < 72; iptt++ ) {
int ism = ( iett<17 ) ? iptt/4+1 : 18+iptt/4+1;
float xvalDQM = DQMVal[ism-1];
float xvalDAQ, xvalDCS;
xvalDAQ = xvalDCS = -1.;
float xcert = -1.;
if ( hDAQ_ ) xvalDAQ = hDAQ_->GetBinContent( iptt+1, iett+1 );
if ( hDCS_ ) xvalDCS = hDCS_->GetBinContent( iptt+1, iett+1 );
if ( xvalDQM == -1 || ( xvalDAQ == -1 && xvalDCS == -1 ) ) {
// problems: DQM empty or DAQ and DCS not available
xcert = 0.0;
} else {
// do not consider the white value of DAQ and DCS (problems with DB)
xcert = std::abs(xvalDQM) * std::abs(xvalDAQ) * std::abs(xvalDCS);
}
if ( meEBDataCertificationSummaryMap_ ) meEBDataCertificationSummaryMap_->setBinContent( iptt+1, iett+1, xcert );
sumCertEB[ism-1] += xcert;
nValidChannelsEB[ism-1]++;
sumCert += xcert;
nValidChannels++;
}
}
if( meEBDataCertificationSummary_ ) {
if( nValidChannels>0 ) meEBDataCertificationSummary_->Fill( sumCert/nValidChannels );
else meEBDataCertificationSummary_->Fill( 0.0 );
}
for (int i = 0; i < 36; i++) {
if( meEBDataCertification_[i] ) {
if( nValidChannelsEB[i]>0 ) meEBDataCertification_[i]->Fill( sumCertEB[i]/nValidChannelsEB[i] );
else meEBDataCertification_[i]->Fill( 0.0 );
}
}
}
| void EBDataCertificationTask::endRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
EndRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 247 of file EBDataCertificationTask.cc.
References abs, cloneME_, dqmStore_, MonitorElement::Fill(), DQMStore::get(), hDAQ_, hDCS_, hDQM_, i, ecalpyutils::ism(), meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, prefixME_, reset(), and MonitorElement::setBinContent().
{
this->reset();
MonitorElement* me;
me = dqmStore_->get(prefixME_ + "/EventInfo/reportSummaryMap");
hDQM_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDQM_ );
me = dqmStore_->get(prefixME_ + "/EventInfo/DAQSummaryMap");
hDAQ_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDAQ_ );
me = dqmStore_->get(prefixME_ + "/EventInfo/DCSSummaryMap");
hDCS_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDCS_ );
float sumCert = 0.;
float sumCertEB[36];
int nValidChannels = 0;
int nValidChannelsEB[36];
for (int i = 0; i < 36; i++) {
sumCertEB[i] = 0.;
nValidChannelsEB[i] = 0;
}
for ( int iett = 0; iett < 34; iett++ ) {
for ( int iptt = 0; iptt < 72; iptt++ ) {
float xvalDQM, xvalDAQ, xvalDCS;
xvalDQM = xvalDAQ = xvalDCS = -1.;
float xcert = -1;
if ( hDQM_ ) xvalDQM = hDQM_->GetBinContent( iptt+1, iett+1 );
if ( hDAQ_ ) xvalDAQ = hDAQ_->GetBinContent( iptt+1, iett+1 );
if ( hDCS_ ) xvalDCS = hDCS_->GetBinContent( iptt+1, iett+1 );
if ( xvalDQM == -1 || ( xvalDAQ == -1 && xvalDCS == -1 ) ) {
// problems: DQM empty or DAQ and DCS not available
xcert = 0.0;
} else {
// do not consider the white value of DAQ and DCS (problems with DB)
xcert = std::abs(xvalDQM) * std::abs(xvalDAQ) * std::abs(xvalDCS);
}
if ( meEBDataCertificationSummaryMap_ ) meEBDataCertificationSummaryMap_->setBinContent( iptt+1, iett+1, xcert );
int ism = ( iett<17 ) ? iptt/4+1 : 18+iptt/4+1;
sumCertEB[ism-1] += xcert;
nValidChannelsEB[ism-1]++;
sumCert += xcert;
nValidChannels++;
}
}
if( meEBDataCertificationSummary_ ) {
if( nValidChannels>0 ) {
meEBDataCertificationSummary_->Fill( sumCert/nValidChannels );
} else {
meEBDataCertificationSummary_->Fill( 0.0 );
}
}
for (int i = 0; i < 36; i++) {
if( meEBDataCertification_[i] ) {
if( nValidChannelsEB[i]>0 ) {
meEBDataCertification_[i]->Fill( sumCertEB[i]/nValidChannelsEB[i] );
} else {
meEBDataCertification_[i]->Fill( 0.0 );
}
}
}
}
| void EBDataCertificationTask::reset | ( | void | ) | [protected] |
Reset.
Definition at line 324 of file EBDataCertificationTask.cc.
References i, meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, and MonitorElement::Reset().
Referenced by beginRun(), endLuminosityBlock(), and endRun().
{
if ( meEBDataCertificationSummary_ ) meEBDataCertificationSummary_->Reset();
for (int i = 0; i < 36; i++) {
if ( meEBDataCertification_[i] ) meEBDataCertification_[i]->Reset();
}
if ( meEBDataCertificationSummaryMap_ ) meEBDataCertificationSummaryMap_->Reset();
}
bool EBDataCertificationTask::cloneME_ [private] |
Definition at line 60 of file EBDataCertificationTask.h.
Referenced by cleanup(), EBDataCertificationTask(), endLuminosityBlock(), and endRun().
DQMStore* EBDataCertificationTask::dqmStore_ [private] |
Definition at line 62 of file EBDataCertificationTask.h.
Referenced by beginJob(), cleanup(), EBDataCertificationTask(), endLuminosityBlock(), and endRun().
bool EBDataCertificationTask::enableCleanup_ [private] |
Definition at line 66 of file EBDataCertificationTask.h.
Referenced by EBDataCertificationTask(), and endJob().
TH2F* EBDataCertificationTask::hDAQ_ [private] |
Definition at line 71 of file EBDataCertificationTask.h.
Referenced by cleanup(), EBDataCertificationTask(), endLuminosityBlock(), and endRun().
TH2F* EBDataCertificationTask::hDCS_ [private] |
Definition at line 72 of file EBDataCertificationTask.h.
Referenced by cleanup(), EBDataCertificationTask(), endLuminosityBlock(), and endRun().
TH2F* EBDataCertificationTask::hDQM_ [private] |
Definition at line 70 of file EBDataCertificationTask.h.
Referenced by cleanup(), EBDataCertificationTask(), and endRun().
TH1F* EBDataCertificationTask::hFrontendByLumi_ [private] |
Definition at line 74 of file EBDataCertificationTask.h.
Referenced by cleanup(), EBDataCertificationTask(), and endLuminosityBlock().
TH1F* EBDataCertificationTask::hIntegrityByLumi_ [private] |
Definition at line 73 of file EBDataCertificationTask.h.
Referenced by cleanup(), EBDataCertificationTask(), and endLuminosityBlock().
TH1F* EBDataCertificationTask::hSynchronizationByLumi_ [private] |
Definition at line 75 of file EBDataCertificationTask.h.
Referenced by cleanup(), EBDataCertificationTask(), and endLuminosityBlock().
MonitorElement* EBDataCertificationTask::meEBDataCertification_[36] [private] |
Definition at line 78 of file EBDataCertificationTask.h.
Referenced by beginJob(), cleanup(), EBDataCertificationTask(), endLuminosityBlock(), endRun(), and reset().
Definition at line 77 of file EBDataCertificationTask.h.
Referenced by beginJob(), cleanup(), EBDataCertificationTask(), endLuminosityBlock(), endRun(), and reset().
Definition at line 79 of file EBDataCertificationTask.h.
Referenced by beginJob(), cleanup(), EBDataCertificationTask(), endLuminosityBlock(), endRun(), and reset().
Definition at line 80 of file EBDataCertificationTask.h.
Definition at line 81 of file EBDataCertificationTask.h.
bool EBDataCertificationTask::mergeRuns_ [private] |
Definition at line 68 of file EBDataCertificationTask.h.
Referenced by beginRun(), and EBDataCertificationTask().
std::string EBDataCertificationTask::prefixME_ [private] |
Definition at line 64 of file EBDataCertificationTask.h.
Referenced by beginJob(), cleanup(), EBDataCertificationTask(), endLuminosityBlock(), and endRun().
1.7.3