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	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () 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 ecalMGPA::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.