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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

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_
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 >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	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::ESTimingTask

(

const edm::ParameterSet &

ps

)

Definition at line 39 of file ESTimingTask.cc.

                                                    {
  digilabel_ = consumes<ESDigiCollection>(ps.getParameter<InputTag>("DigiLabel"));
  prefixME_ = ps.getUntrackedParameter<string>("prefixME", "EcalPreshower");
 
  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);
}

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

~ESTimingTask()

ESTimingTask::~ESTimingTask ( )

override

Definition at line 77 of file ESTimingTask.cc.

                            {
  delete fit_;
  delete htESP_;
  delete htESM_;
}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 83 of file ESTimingTask.cc.

                                                                         {
  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());
}

References ecalLiteDTU::adc(), ESSample::adc(), edm::DataFrameContainer::begin(), MillePedeFileConverter_cfg::e, edm::DataFrameContainer::end(), JetChargeProducer_cfi::exp, mps_fire::i, ESDataFrame::id(), dqmiolumiharvest::j, dqmdumpme::k, dqm-mbProfile::log, ESDataFrame::sample(), ESDataFrame::size(), mps_update::status, FrontierCondition_GT_autoExpress_cfi::t0, and ecaldqm::zside().

bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 57 of file ESTimingTask.cc.

                                                                                         {
  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);
}

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().

set()

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

private

Definition at line 173 of file ESTimingTask.cc.

                                              {
  es.get<ESGainRcd>().get(esgain_);
  const ESGain* gain = esgain_.product();
 
  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;
}

References PedestalClient_cfi::gain, edm::EventSetup::get(), get, and createfilelist::int.

Member Data Documentation

digilabel_

edm::EDGetTokenT<ESDigiCollection> ESTimingTask::digilabel_

private

Definition at line 32 of file ESTimingTask.h.

eCount_

int ESTimingTask::eCount_

private

Definition at line 44 of file ESTimingTask.h.

esgain_

edm::ESHandle<ESGain> ESTimingTask::esgain_

private

Definition at line 38 of file ESTimingTask.h.

fit_

TF1* ESTimingTask::fit_

private

Definition at line 40 of file ESTimingTask.h.

h2DTiming_

MonitorElement* ESTimingTask::h2DTiming_

private

Definition at line 36 of file ESTimingTask.h.

htESM_

TH1F* ESTimingTask::htESM_

private

Definition at line 42 of file ESTimingTask.h.

htESP_

TH1F* ESTimingTask::htESP_

private

Definition at line 41 of file ESTimingTask.h.

hTiming_

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

private

Definition at line 35 of file ESTimingTask.h.

n_

Double_t ESTimingTask::n_

private

Definition at line 45 of file ESTimingTask.h.

prefixME_

std::string ESTimingTask::prefixME_

private

Definition at line 33 of file ESTimingTask.h.

runNum_

int ESTimingTask::runNum_

private

Definition at line 44 of file ESTimingTask.h.

wc_

Double_t ESTimingTask::wc_

private

Definition at line 45 of file ESTimingTask.h.