ESTimingTask Class Reference

#include <ESTimingTask.h>

Inheritance diagram for ESTimingTask:

Public Member Functions
	ESTimingTask (const edm::ParameterSet &ps)
	~ESTimingTask () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	set (const edm::EventSetup &es)

Private Attributes
edm::EDGetTokenT < ESDigiCollection >	digilabel_
int	eCount_
edm::ESHandle< ESGain >	esgain_
edm::ESGetToken< ESGain, ESGainRcd >	esgainToken_
TF1 *	fit_
MonitorElement *	h2DTiming_
TH1F *	htESM_
TH1F *	htESP_
MonitorElement *	hTiming_ [2][2]
Double_t	n_
std::string	prefixME_
int	runNum_
Double_t	wc_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 21 of file ESTimingTask.h.

Constructor & Destructor Documentation

ESTimingTask::ESTimingTask

(

const edm::ParameterSet &

ps

)

Definition at line 39 of file ESTimingTask.cc.

References DeDxTools::esConsumes(), fitf(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, and dqmiolumiharvest::j.

                                                    {
  digilabel_ = consumes<ESDigiCollection>(ps.getParameter<InputTag>("DigiLabel"));
  prefixME_ = ps.getUntrackedParameter<string>("prefixME", "EcalPreshower");
  esgainToken_ = esConsumes();
  eCount_ = 0;

  fit_ = new TF1("fitShape", fitf, -200, 200, 4);
  fit_->SetParameters(50, 10, 0, 0);

  //Histogram init
  for (int i = 0; i < 2; ++i)
    for (int j = 0; j < 2; ++j)
      hTiming_[i][j] = nullptr;

  htESP_ = new TH1F("htESP", "Timing ES+", 81, -20.5, 20.5);
  htESM_ = new TH1F("htESM", "Timing ES-", 81, -20.5, 20.5);
}

ESTimingTask::~ESTimingTask ( )

override

Definition at line 77 of file ESTimingTask.cc.

                            {
  delete fit_;
  delete htESP_;
  delete htESM_;
}

Member Function Documentation

void ESTimingTask::analyze ( const edm::Event &  e, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 83 of file ESTimingTask.cc.

References ESSample::adc(), gpuClustering::adc, funct::exp(), edm::Event::getByToken(), mps_fire::i, ESDataFrame::id(), edm::EventBase::id(), dqmiolumiharvest::j, isotrackApplyRegressor::k, log, edm::EventID::run(), ESDataFrame::sample(), ESDataFrame::size(), mps_update::status, FrontierCondition_GT_autoExpress_cfi::t0, ESDetId::zside(), and ecaldqm::zside().

                                                                         {
  set(iSetup);

  runNum_ = e.id().run();
  eCount_++;

  htESP_->Reset();
  htESM_->Reset();

  //Digis
  int zside, plane, ix, iy, is;
  double adc[3];
  //  double para[10];
  //double tx[3] = {-5., 20., 45.};
  Handle<ESDigiCollection> digis;
  if (e.getByToken(digilabel_, digis)) {
    for (ESDigiCollection::const_iterator digiItr = digis->begin(); digiItr != digis->end(); ++digiItr) {
      ESDataFrame dataframe = (*digiItr);
      ESDetId id = dataframe.id();

      zside = id.zside();
      plane = id.plane();
      ix = id.six();
      iy = id.siy();
      is = id.strip();

      //if (zside==1 && plane==1 && ix==15 && iy==6) continue;
      if (zside == 1 && plane == 1 && ix == 7 && iy == 28)
        continue;
      if (zside == 1 && plane == 1 && ix == 24 && iy == 9 && is == 21)
        continue;
      if (zside == -1 && plane == 2 && ix == 35 && iy == 17 && is == 23)
        continue;

      int i = (zside == 1) ? 0 : 1;
      int j = plane - 1;

      for (int k = 0; k < dataframe.size(); ++k)
        adc[k] = dataframe.sample(k).adc();

      double status = 0;
      if (adc[1] < 200)
        status = 1;
      if (fabs(adc[0]) > 10)
        status = 1;
      if (adc[1] < 0 || adc[2] < 0)
        status = 1;
      if (adc[0] > adc[1] || adc[0] > adc[2])
        status = 1;
      if (adc[2] > adc[1])
        status = 1;

      if (int(status) == 0) {
        double A1 = adc[1];
        double A2 = adc[2];
        double DeltaT = 25.;
        double aaa = (A2 > 0 && A1 > 0) ? log(A2 / A1) / n_ : 20.;  // if A1=0, t0=20
        double bbb = wc_ / n_ * DeltaT;
        double ccc = exp(aaa + bbb);

        double t0 = (2. - ccc) / (1. - ccc) * DeltaT - 5;
        hTiming_[i][j]->Fill(t0);
        //cout<<"t0 : "<<t0<<endl;
        /*
    TGraph *gr = new TGraph(3, tx, adc);
    fit_->SetParameters(50, 10, wc_, n_);
    fit_->FixParameter(2, wc_);
    fit_->FixParameter(3, n_);
    fit_->Print();
    gr->Fit("fitShape", "MQ");
    fit_->GetParameters(para); 
    delete gr;
    //hTiming_[i][j]->Fill(para[1]);
    */
        //cout<<"ADC : "<<zside<<" "<<plane<<" "<<ix<<" "<<iy<<" "<<is<<" "<<adc[0]<<" "<<adc[1]<<" "<<adc[2]<<" "<<para[1]<<" "<<wc_<<" "<<n_<<endl;

        if (zside == 1)
          htESP_->Fill(t0);
        else if (zside == -1)
          htESM_->Fill(t0);
      }
    }
  } else {
    LogWarning("ESTimingTask") << "DigiCollection not available";
  }

  if (htESP_->GetEntries() > 0 && htESM_->GetEntries() > 0)
    h2DTiming_->Fill(htESM_->GetMean(), htESP_->GetMean());
}

void ESTimingTask::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  , edm::EventSetup const &    )

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 57 of file ESTimingTask.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), timingPdfMaker::histo, mps_fire::i, dqmiolumiharvest::j, dqm::impl::MonitorElement::setAxisTitle(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                         {
  iBooker.setCurrentFolder(prefixME_ + "/ESTimingTask");

  //Booking Histograms
  //Notice: Change ESRenderPlugin under DQM/RenderPlugins/src if you change this histogram name.
  char histo[200];
  for (int i = 0; i < 2; ++i)
    for (int j = 0; j < 2; ++j) {
      int iz = (i == 0) ? 1 : -1;
      sprintf(histo, "ES Timing Z %d P %d", iz, j + 1);
      hTiming_[i][j] = iBooker.book1D(histo, histo, 81, -20.5, 20.5);
      hTiming_[i][j]->setAxisTitle("ES Timing (ns)", 1);
    }

  sprintf(histo, "ES 2D Timing");
  h2DTiming_ = iBooker.book2D(histo, histo, 81, -20.5, 20.5, 81, -20.5, 20.5);
  h2DTiming_->setAxisTitle("ES- Timing (ns)", 1);
  h2DTiming_->setAxisTitle("ES+ Timing (ns)", 2);
}

void ESTimingTask::set ( const edm::EventSetup &  es )

private

Definition at line 173 of file ESTimingTask.cc.

References edm::EventSetup::getData(), and ESGain::getESGain().

                                              {
  const ESGain* gain = &es.getData(esgainToken_);

  int ESGain = (int)gain->getESGain();

  if (ESGain == 1) {  // LG
    wc_ = 0.0837264;
    n_ = 2.016;
  } else {  // HG
    wc_ = 0.07291;
    n_ = 1.798;
  }

  //cout<<"gain : "<<ESGain<<endl;
  //cout<<wc_<<" "<<n_<<endl;
}

Member Data Documentation

edm::EDGetTokenT<ESDigiCollection> ESTimingTask::digilabel_

private

Definition at line 32 of file ESTimingTask.h.

int ESTimingTask::eCount_

private

Definition at line 45 of file ESTimingTask.h.

edm::ESHandle<ESGain> ESTimingTask::esgain_

private

Definition at line 39 of file ESTimingTask.h.

edm::ESGetToken<ESGain, ESGainRcd> ESTimingTask::esgainToken_

private

Definition at line 38 of file ESTimingTask.h.

TF1* ESTimingTask::fit_

private

Definition at line 41 of file ESTimingTask.h.

MonitorElement* ESTimingTask::h2DTiming_

private

Definition at line 36 of file ESTimingTask.h.

TH1F* ESTimingTask::htESM_

private

Definition at line 43 of file ESTimingTask.h.

TH1F* ESTimingTask::htESP_

private

Definition at line 42 of file ESTimingTask.h.

MonitorElement* ESTimingTask::hTiming_[2][2]

private

Definition at line 35 of file ESTimingTask.h.

Double_t ESTimingTask::n_

private

Definition at line 46 of file ESTimingTask.h.

std::string ESTimingTask::prefixME_

private

Definition at line 33 of file ESTimingTask.h.

int ESTimingTask::runNum_

private

Definition at line 45 of file ESTimingTask.h.

Double_t ESTimingTask::wc_

private

Definition at line 46 of file ESTimingTask.h.