ecaldqm::MESetTrend Class Reference

#include <MESetTrend.h>

Inheritance diagram for ecaldqm::MESetTrend:

Public Member Functions
void	book ()
void	fill (DetId const &, double, double _wy=1., double _w=1.)
void	fill (unsigned, double, double _wy=1., double _w=1.)
void	fill (double, double _wy=1., double _w=1.)
double	getBinContent (DetId const &, int _bin=0) const
double	getBinContent (unsigned, int _bin=0) const
double	getBinEntries (DetId const &, int _bin=0) const
double	getBinEntries (unsigned, int _bin=0) const
double	getBinError (DetId const &, int _bin=0) const
double	getBinError (unsigned, int _bin=0) const
bool	getMinutely () const
time_t	getTimeZero () const
	MESetTrend (std::string const &, MEData const &, bool _readOnly=false)
void	setBinContent (DetId const &, double, double _err=0.)
void	setBinContent (unsigned, double, double _err=0.)
void	setBinEntries (DetId const &, double)
void	setBinEntries (unsigned, double)
void	setMinutely (bool _minutely)
void	setTimeZero (time_t _t0)
	~MESetTrend ()
Public Member Functions inherited from ecaldqm::MESetEcal
std::vector< std::string >	generateNames () const
	MESetEcal (std::string const &, MEData const &, int logicalDimensions_, bool _readOnly=false)
void	reset (double _content=0., double _err=0., double _entries=0.)
bool	retrieve () const
	~MESetEcal ()
Public Member Functions inherited from ecaldqm::MESet
virtual void	clear () const
virtual void	fill (EcalElectronicsId const &, double _wx=1., double _wy=1., double _w=1.)
virtual double	getBinContent (EcalElectronicsId const &, int _bin=0) const
virtual double	getBinEntries (EcalElectronicsId const &, int _bin=0) const
virtual double	getBinError (EcalElectronicsId const &, int _bin=0) const
BinService::BinningType	getBinType () const
std::string const &	getDir () const
virtual MonitorElement const *	getME (unsigned _offset) const
std::string const &	getName () const
BinService::ObjectType	getObjType () const
bool	isActive () const
	MESet (std::string const &, MEData const &, bool _readOnly=false)
void	name (std::map< std::string, std::string > const &) const
virtual void	resetAll (double _content=0., double _err=0., double _entries=0.)
virtual void	setAxisTitle (std::string const &, int _axis=1)
virtual void	setBinContent (EcalElectronicsId const &, double, double _err=0.)
virtual void	setBinEntries (EcalElectronicsId const &, double)
virtual void	setBinLabel (unsigned, int, std::string const &, int _axis=1)
void	setDir (std::string const &_dir)
void	setName (std::string const &_name)
virtual	~MESet ()

Private Member Functions
bool	shift_ (time_t)

Private Attributes
bool	minutely_
time_t	t0_
time_t	tLow_

Additional Inherited Members
Protected Member Functions inherited from ecaldqm::MESetEcal
virtual void	fill_ (double)
virtual void	find_ (uint32_t) const
Protected Member Functions inherited from ecaldqm::MESet
virtual void	fill_ (unsigned, int, double)
virtual void	fill_ (unsigned, double, double, double)
virtual double	getBinContent_ (unsigned, int) const
virtual double	getBinEntries_ (unsigned, int) const
virtual double	getBinError_ (unsigned, int) const
virtual void	setBinContent_ (unsigned, int, double, double)
virtual void	setBinEntries_ (unsigned, int, double)
Protected Attributes inherited from ecaldqm::MESetEcal
std::pair< unsigned, std::vector< int > >	cache_
uint32_t	cacheId_
const unsigned	logicalDimensions_
Protected Attributes inherited from ecaldqm::MESet
bool	active_
MEData const *	data_
std::string	dir_
std::vector< MonitorElement * >	mes_
std::string	name_
bool	readOnly_
Static Protected Attributes inherited from ecaldqm::MESet
static BinService const *	binService_
static DQMStore *	dqmStore_

Detailed Description

Definition at line 11 of file MESetTrend.h.

Constructor & Destructor Documentation

ecaldqm::MESetTrend::MESetTrend	(	std::string const &	_fullpath,
		MEData const &	_data,
		bool	_readOnly = false
	)

Definition at line 7 of file MESetTrend.cc.

References EcalDQMBinningService::AxisSpecs::edges, edm::hlt::Exception, and ecaldqm::MEData::xaxis.

                                                                                         :
    MESetEcal(_fullpath, _data, 1, _readOnly),
    t0_(0),
    minutely_(false),
    tLow_(0)
  {
    if(!_data.xaxis || _data.xaxis->edges)
      throw cms::Exception("InvalidConfiguration") << "MESetTrend";
  }

ecaldqm::MESetTrend::~MESetTrend

(

)

Definition at line 17 of file MESetTrend.cc.

  {
  }

Member Function Documentation

void ecaldqm::MESetTrend::book ( )

virtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 22 of file MESetTrend.cc.

References ecaldqm::MESetEcal::book(), ecaldqm::MEData::btype, edm::conversion(), ecaldqm::MESet::data_, EcalDQMBinningService::AxisSpecs::high, ecaldqm::MEData::kind, EcalDQMBinningService::AxisSpecs::low, ecaldqm::MESet::mes_, minutely_, EcalDQMBinningService::AxisSpecs::nbins, ecaldqm::MEData::otype, ecaldqm::MEData::pathName, t0_, groupFilesInBlocks::temp, tLow_, create_public_lumi_plots::width, ecaldqm::MEData::xaxis, ecaldqm::MEData::yaxis, and ecaldqm::MEData::zaxis.

  {
    int conversion(minutely_ ? 60 : 3600);
    time_t width((data_->xaxis->high - data_->xaxis->low) * conversion);

    tLow_ = t0_;

    MEData const* temp(data_);
    BinService::AxisSpecs xaxis(*temp->xaxis);
    xaxis.low = tLow_;
    xaxis.high = tLow_ + width;

    data_ = new MEData(temp->pathName, temp->otype, temp->btype, temp->kind, &xaxis, temp->yaxis, temp->zaxis);

    MESetEcal::book();

    delete data_;
    data_ = temp;

    // if yaxis was variable bin size, xaxis will be booked as variable too

    for(unsigned iME(0); iME < mes_.size(); iME++){
      TAxis* axis(mes_[iME]->getTH1()->GetXaxis());
      if(axis->IsVariableBinSize())
    axis->Set(data_->xaxis->nbins, data_->xaxis->low, data_->xaxis->high);
    }
  }

void ecaldqm::MESetTrend::fill ( DetId const &  _id, double  _t, double  _wy = 1., double  _w = 1.  )

virtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 51 of file MESetTrend.cc.

References ecaldqm::MESet::active_, ecaldqm::MESet::binService_, ecaldqm::MESet::data_, ecaldqm::MESet::fill_(), EcalDQMBinningService::findOffset(), evf::evtn::offset(), ecaldqm::MEData::otype, and shift_().

  {
    if(!active_) return;

    unsigned offset(binService_->findOffset(data_->otype, _id));

    if(shift_(time_t(_t)))
      MESet::fill_(offset, _t, _wy, _w);
  }

void ecaldqm::MESetTrend::fill ( unsigned  _dcctccid, double  _t, double  _wy = 1., double  _w = 1.  )

virtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 62 of file MESetTrend.cc.

References ecaldqm::MESet::active_, ecaldqm::MESet::binService_, ecaldqm::MEData::btype, ecaldqm::MESet::data_, ecaldqm::MESet::fill_(), EcalDQMBinningService::findOffset(), evf::evtn::offset(), ecaldqm::MEData::otype, and shift_().

  {
    if(!active_) return;

    unsigned offset(binService_->findOffset(data_->otype, data_->btype, _dcctccid));

    if(shift_(time_t(_t)))
      MESet::fill_(offset, _t, _wy, _w);
  }

void ecaldqm::MESetTrend::fill ( double  _t, double  _wy = 1., double  _w = 1.  )

virtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 73 of file MESetTrend.cc.

References ecaldqm::MESet::active_, ecaldqm::MESet::fill_(), ecaldqm::MESet::mes_, and shift_().

  {
    if(!active_) return;
    if(mes_.size() != 1)
      throw cms::Exception("InvalidCall") << "MESet type incompatible" << std::endl;

    if(shift_(time_t(_t)))
      MESet::fill_(0, _t, _wy, _w);
  }

double ecaldqm::MESetTrend::getBinContent ( DetId const &  , int  _bin = 0  ) const

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 29 of file MESetTrend.h.

{ return 0.; }

double ecaldqm::MESetTrend::getBinContent ( unsigned  , int  _bin = 0  ) const

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 30 of file MESetTrend.h.

{ return 0.; }

double ecaldqm::MESetTrend::getBinEntries ( DetId const &  , int  _bin = 0  ) const

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 35 of file MESetTrend.h.

{ return 0.; }

double ecaldqm::MESetTrend::getBinEntries ( unsigned  , int  _bin = 0  ) const

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 36 of file MESetTrend.h.

{ return 0.; }

double ecaldqm::MESetTrend::getBinError ( DetId const &  , int  _bin = 0  ) const

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 32 of file MESetTrend.h.

{ return 0.; }

double ecaldqm::MESetTrend::getBinError ( unsigned  , int  _bin = 0  ) const

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 33 of file MESetTrend.h.

{ return 0.; }

bool ecaldqm::MESetTrend::getMinutely ( ) const

inline

Definition at line 41 of file MESetTrend.h.

References minutely_.

{ return minutely_; }

time_t ecaldqm::MESetTrend::getTimeZero ( ) const

inline

Definition at line 39 of file MESetTrend.h.

References t0_.

{ return t0_; }

void ecaldqm::MESetTrend::setBinContent ( DetId const &  , double  , double  _err = 0.  )

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 23 of file MESetTrend.h.

{}

void ecaldqm::MESetTrend::setBinContent ( unsigned  , double  , double  _err = 0.  )

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 24 of file MESetTrend.h.

{}

void ecaldqm::MESetTrend::setBinEntries ( DetId const &  , double    )

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 26 of file MESetTrend.h.

{}

void ecaldqm::MESetTrend::setBinEntries ( unsigned  , double    )

inlinevirtual

Reimplemented from ecaldqm::MESetEcal.

Definition at line 27 of file MESetTrend.h.

{}

void ecaldqm::MESetTrend::setMinutely ( bool  _minutely )

inline

Definition at line 40 of file MESetTrend.h.

References minutely_.

{ minutely_ = _minutely; }

void ecaldqm::MESetTrend::setTimeZero ( time_t  _t0 )

inline

Definition at line 38 of file MESetTrend.h.

References t0_.

{ t0_ = _t0; }

bool ecaldqm::MESetTrend::shift_ ( time_t  _t )

private

Definition at line 84 of file MESetTrend.cc.

References edm::conversion(), ecaldqm::MESet::data_, alignCSCRings::dest, MonitorElement::DQM_KIND_TH1F, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TPROFILE, end, python.tagInventory::entries, MonitorElement::getBinContent(), MonitorElement::getBinEntries(), MonitorElement::getBinError(), MonitorElement::getNbinsY(), MonitorElement::getTH1(), EcalDQMBinningService::AxisSpecs::high, ecaldqm::MEData::kind, EcalDQMBinningService::AxisSpecs::low, ecaldqm::MESet::mes_, minutely_, EcalDQMBinningService::AxisSpecs::nbins, plotscripts::rms(), MonitorElement::setBinContent(), MonitorElement::setBinEntries(), MonitorElement::setBinError(), MonitorElement::setEntries(), mathSSE::sqrt(), dqm_diff::start, relval_parameters_module::step, tLow_, create_public_lumi_plots::width, and ecaldqm::MEData::xaxis.

Referenced by fill().

  {
    int conversion(minutely_ ? 60 : 3600);
    time_t width((data_->xaxis->high - data_->xaxis->low) * conversion);

    time_t tHigh(tLow_ + width);
    int nbinsX(data_->xaxis->nbins);
    MonitorElement::Kind kind(data_->kind);

    int dtPerBin(width / nbinsX);
    int dBin(0);

    if(_t >= tLow_ && _t < tHigh)
      return true;
    else if(_t >= tHigh)
      dBin = (_t - tHigh) / dtPerBin + 1;
    else if(_t < tLow_){
      int maxBin(0);

      for(unsigned iME(0); iME < mes_.size(); iME++){
    MonitorElement* me(mes_[iME]);

    bool filled(false);
    int iMax(nbinsX + 1);
    while(--iMax > 0 && !filled){
      switch(kind){
      case  MonitorElement::DQM_KIND_TH1F:
        if(me->getBinContent(iMax) != 0) filled = true;
        break;
      case MonitorElement::DQM_KIND_TPROFILE:
        if(me->getBinEntries(iMax) != 0) filled = true;
        break;
      case MonitorElement::DQM_KIND_TH2F:
        for(int iy(1); iy <= me->getNbinsY(); iy++)
          if(me->getBinContent(iMax, iy) != 0) filled = true;
        break;
      default:
        return false;
      }
    }

    if(iMax > maxBin) maxBin = iMax;
      }

      if(_t < tLow_ - (nbinsX - maxBin) * dtPerBin) return false;

      dBin = (_t - dtPerBin - tLow_) / dtPerBin;
    }

    int start(dBin > 0 ? dBin + 1 : nbinsX + dBin);
    int end(dBin > 0 ? nbinsX + 1 : 0);
    int step(dBin > 0 ? 1 : -1);

    tLow_ += dBin * dtPerBin;
    tHigh += dBin * dtPerBin;

    for(unsigned iME(0); iME < mes_.size(); iME++){
      MonitorElement* me(mes_[iME]);

      me->setEntries(0.);

      double entries(0.);

      for(int ix(start); (dBin > 0 ? (ix < end) : (ix > end)); ix += step){
    switch(kind){
    case  MonitorElement::DQM_KIND_TH1F:
      entries += me->getBinContent(ix);
      me->setBinContent(ix - dBin, me->getBinContent(ix));
      me->setBinError(ix - dBin, me->getBinError(ix));
      me->setBinContent(ix, 0.);
      me->setBinError(ix, 0.);
      break;
    case MonitorElement::DQM_KIND_TPROFILE:
      entries += me->getBinEntries(ix);
      me->setBinEntries(ix - dBin, me->getBinEntries(ix));
      me->setBinContent(ix - dBin, me->getBinContent(ix) * me->getBinEntries(ix));
      if(me->getBinEntries(ix) > 0){
        double rms(me->getBinError(ix) * std::sqrt(me->getBinEntries(ix)));
        double sumw2((rms * rms + me->getBinContent(ix) * me->getBinContent(ix)) * me->getBinEntries(ix));
        me->setBinError(ix - dBin, std::sqrt(sumw2));
      }
      me->setBinEntries(ix, 0.);
      me->setBinContent(ix, 0.);
      me->setBinError(ix, 0.);
      break;
    case MonitorElement::DQM_KIND_TH2F:
      for(int iy(1); iy <= me->getNbinsY(); iy++){
        int orig(me->getTH1()->GetBin(ix, iy));
        int dest(me->getTH1()->GetBin(ix - dBin, iy));
        entries += me->getBinContent(orig);
        me->setBinContent(dest, me->getBinContent(orig));
        me->setBinError(dest, me->getBinError(orig));
        me->setBinContent(orig, 0.);
        me->setBinError(orig, 0.);
      }
      break;
    default:
      break;
    }
      }

      me->setEntries(entries);
      me->getTH1()->GetXaxis()->SetLimits(tLow_, tHigh);
    }

    return true;
  }

Member Data Documentation

bool ecaldqm::MESetTrend::minutely_

private

Definition at line 47 of file MESetTrend.h.

Referenced by book(), getMinutely(), setMinutely(), and shift_().

time_t ecaldqm::MESetTrend::t0_

private

Definition at line 46 of file MESetTrend.h.

Referenced by book(), getTimeZero(), and setTimeZero().

time_t ecaldqm::MESetTrend::tLow_

private

Definition at line 49 of file MESetTrend.h.

Referenced by book(), and shift_().