EcalDigisValidation Class Reference

#include <EcalDigisValidation.h>

Inheritance diagram for EcalDigisValidation:

Public Member Functions
void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override
	EcalDigisValidation (const edm::ParameterSet &ps)
	Constructor. More...
	~EcalDigisValidation () override
	Destructor. More...
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
	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

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &c) override
	Analyze. More...
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Types
typedef std::map< uint32_t, float, std::less< uint32_t > >	MapType

Private Member Functions
void	checkCalibrations (edm::EventSetup const &c)

Private Attributes
double	barrelADCtoGeV_
edm::EDGetTokenT< CrossingFrame< PCaloHit > >	crossingFramePCaloHitEBToken_
edm::EDGetTokenT< CrossingFrame< PCaloHit > >	crossingFramePCaloHitEEToken_
edm::EDGetTokenT< CrossingFrame< PCaloHit > >	crossingFramePCaloHitESToken_
edm::EDGetTokenT< EBDigiCollection >	EBdigiCollectionToken_
edm::EDGetTokenT< EEDigiCollection >	EEdigiCollectionToken_
double	endcapADCtoGeV_
edm::EDGetTokenT< ESDigiCollection >	ESdigiCollectionToken_
edm::EDGetTokenT< edm::SimTrackContainer >	g4TkInfoToken_
edm::EDGetTokenT< edm::SimVertexContainer >	g4VtxInfoToken_
std::map< int, double, std::less< int > >	gainConv_
edm::EDGetTokenT< edm::HepMCProduct >	HepMCToken_
MonitorElement *	meEBDigiSimRatio_
MonitorElement *	meEBDigiSimRatiogt100ADC_
MonitorElement *	meEBDigiSimRatiogt10ADC_
MonitorElement *	meEEDigiSimRatio_
MonitorElement *	meEEDigiSimRatiogt100ADC_
MonitorElement *	meEEDigiSimRatiogt20ADC_
MonitorElement *	meGunEnergy_
MonitorElement *	meGunEta_
MonitorElement *	meGunPhi_
std::string	outputFile_
bool	verbose_

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 47 of file EcalDigisValidation.h.

Member Typedef Documentation

MapType

typedef std::map<uint32_t, float, std::less<uint32_t> > EcalDigisValidation::MapType

private

Definition at line 48 of file EcalDigisValidation.h.

Constructor & Destructor Documentation

EcalDigisValidation()

EcalDigisValidation::EcalDigisValidation

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 14 of file EcalDigisValidation.cc.

    : HepMCToken_(consumes<edm::HepMCProduct>(edm::InputTag(ps.getParameter<std::string>("moduleLabelMC")))),
      g4TkInfoToken_(consumes<edm::SimTrackContainer>(edm::InputTag(ps.getParameter<std::string>("moduleLabelG4")))),
      g4VtxInfoToken_(consumes<edm::SimVertexContainer>(edm::InputTag(ps.getParameter<std::string>("moduleLabelG4")))),
      EBdigiCollectionToken_(consumes<EBDigiCollection>(ps.getParameter<edm::InputTag>("EBdigiCollection"))),
      EEdigiCollectionToken_(consumes<EEDigiCollection>(ps.getParameter<edm::InputTag>("EEdigiCollection"))),
      ESdigiCollectionToken_(consumes<ESDigiCollection>(ps.getParameter<edm::InputTag>("ESdigiCollection"))),
      crossingFramePCaloHitEBToken_(consumes<CrossingFrame<PCaloHit> >(edm::InputTag(
          std::string("mix"), ps.getParameter<std::string>("moduleLabelG4") + std::string("EcalHitsEB")))),
      crossingFramePCaloHitEEToken_(consumes<CrossingFrame<PCaloHit> >(edm::InputTag(
          std::string("mix"), ps.getParameter<std::string>("moduleLabelG4") + std::string("EcalHitsEE")))),
      crossingFramePCaloHitESToken_(consumes<CrossingFrame<PCaloHit> >(edm::InputTag(
          std::string("mix"), ps.getParameter<std::string>("moduleLabelG4") + std::string("EcalHitsES")))) {
  // DQM ROOT output
  outputFile_ = ps.getUntrackedParameter<std::string>("outputFile", "");
 
  if (!outputFile_.empty()) {
    edm::LogInfo("OutputInfo") << " Ecal Digi Task histograms will be saved to '" << outputFile_.c_str() << "'";
  } else {
    edm::LogInfo("OutputInfo") << " Ecal Digi Task histograms will NOT be saved";
  }
 
  // verbosity switch
  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
 
  gainConv_[1] = 1.;
  gainConv_[2] = 2.;
  gainConv_[3] = 12.;
  gainConv_[0] = 12.;  // saturated channels
  barrelADCtoGeV_ = 0.035;
  endcapADCtoGeV_ = 0.06;
 
  meGunEnergy_ = nullptr;
  meGunEta_ = nullptr;
  meGunPhi_ = nullptr;
 
  meEBDigiSimRatio_ = nullptr;
  meEEDigiSimRatio_ = nullptr;
 
  meEBDigiSimRatiogt10ADC_ = nullptr;
  meEEDigiSimRatiogt20ADC_ = nullptr;
 
  meEBDigiSimRatiogt100ADC_ = nullptr;
  meEEDigiSimRatiogt100ADC_ = nullptr;
}

References barrelADCtoGeV_, endcapADCtoGeV_, gainConv_, edm::ParameterSet::getUntrackedParameter(), meEBDigiSimRatio_, meEBDigiSimRatiogt100ADC_, meEBDigiSimRatiogt10ADC_, meEEDigiSimRatio_, meEEDigiSimRatiogt100ADC_, meEEDigiSimRatiogt20ADC_, meGunEnergy_, meGunEta_, meGunPhi_, outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

~EcalDigisValidation()

EcalDigisValidation::~EcalDigisValidation ( )

override

Destructor.

Definition at line 60 of file EcalDigisValidation.cc.

{}

Member Function Documentation

analyze()

void EcalDigisValidation::analyze ( edm::Event const &  e, edm::EventSetup const &  c  )

overrideprotectedvirtual

Analyze.

Reimplemented from DQMEDAnalyzer.

Definition at line 97 of file EcalDigisValidation.cc.

                                                                           {
  edm::LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
 
  std::vector<SimTrack> theSimTracks;
  std::vector<SimVertex> theSimVertexes;
 
  edm::Handle<edm::HepMCProduct> MCEvt;
  edm::Handle<edm::SimTrackContainer> SimTk;
  edm::Handle<edm::SimVertexContainer> SimVtx;
  edm::Handle<CrossingFrame<PCaloHit> > crossingFrame;
  edm::Handle<EBDigiCollection> EcalDigiEB;
  edm::Handle<EEDigiCollection> EcalDigiEE;
  edm::Handle<ESDigiCollection> EcalDigiES;
 
  bool skipMC = false;
  e.getByToken(HepMCToken_, MCEvt);
  if (!MCEvt.isValid()) {
    skipMC = true;
  }
  e.getByToken(g4TkInfoToken_, SimTk);
  e.getByToken(g4VtxInfoToken_, SimVtx);
 
  const EBDigiCollection* EBdigis = nullptr;
  const EEDigiCollection* EEdigis = nullptr;
  const ESDigiCollection* ESdigis = nullptr;
 
  bool isBarrel = true;
  e.getByToken(EBdigiCollectionToken_, EcalDigiEB);
  if (EcalDigiEB.isValid()) {
    EBdigis = EcalDigiEB.product();
    LogDebug("DigiInfo") << "total # EBdigis: " << EBdigis->size();
    if (EBdigis->empty())
      isBarrel = false;
  } else {
    isBarrel = false;
  }
 
  bool isEndcap = true;
  e.getByToken(EEdigiCollectionToken_, EcalDigiEE);
  if (EcalDigiEE.isValid()) {
    EEdigis = EcalDigiEE.product();
    LogDebug("DigiInfo") << "total # EEdigis: " << EEdigis->size();
    if (EEdigis->empty())
      isEndcap = false;
  } else {
    isEndcap = false;
  }
 
  bool isPreshower = true;
  e.getByToken(ESdigiCollectionToken_, EcalDigiES);
  if (EcalDigiES.isValid()) {
    ESdigis = EcalDigiES.product();
    LogDebug("DigiInfo") << "total # ESdigis: " << ESdigis->size();
    if (ESdigis->empty())
      isPreshower = false;
  } else {
    isPreshower = false;
  }
 
  theSimTracks.insert(theSimTracks.end(), SimTk->begin(), SimTk->end());
  theSimVertexes.insert(theSimVertexes.end(), SimVtx->begin(), SimVtx->end());
 
  if (!skipMC) {
    double theGunEnergy = 0.;
    for (HepMC::GenEvent::particle_const_iterator p = MCEvt->GetEvent()->particles_begin();
         p != MCEvt->GetEvent()->particles_end();
         ++p) {
      theGunEnergy = (*p)->momentum().e();
      double htheta = (*p)->momentum().theta();
      double heta = -log(tan(htheta * 0.5));
      double hphi = (*p)->momentum().phi();
      hphi = (hphi >= 0) ? hphi : hphi + 2 * M_PI;
      hphi = hphi / M_PI * 180.;
      LogDebug("EventInfo") << "Particle gun type form MC = " << abs((*p)->pdg_id()) << "\n"
                            << "Energy = " << (*p)->momentum().e() << " Eta = " << heta << " Phi = " << hphi;
 
      if (meGunEnergy_)
        meGunEnergy_->Fill(theGunEnergy);
      if (meGunEta_)
        meGunEta_->Fill(heta);
      if (meGunPhi_)
        meGunPhi_->Fill(hphi);
    }
  }
 
  int nvtx = 0;
  for (std::vector<SimVertex>::iterator isimvtx = theSimVertexes.begin(); isimvtx != theSimVertexes.end(); ++isimvtx) {
    LogDebug("EventInfo") << " Vertex index = " << nvtx << " event Id = " << isimvtx->eventId().rawId() << "\n"
                          << " vertex dump: " << *isimvtx;
    ++nvtx;
  }
 
  int ntrk = 0;
  for (std::vector<SimTrack>::iterator isimtrk = theSimTracks.begin(); isimtrk != theSimTracks.end(); ++isimtrk) {
    LogDebug("EventInfo") << " Track index = " << ntrk << " track Id = " << isimtrk->trackId()
                          << " event Id = " << isimtrk->eventId().rawId() << "\n"
                          << " track dump: " << *isimtrk;
    ++ntrk;
  }
 
  // BARREL
 
  // loop over simHits
 
  if (isBarrel) {
    e.getByToken(crossingFramePCaloHitEBToken_, crossingFrame);
    const MixCollection<PCaloHit> barrelHits(crossingFrame.product());
 
    MapType ebSimMap;
    for (auto const& iHit : barrelHits) {
      EBDetId ebid = EBDetId(iHit.id());
 
      LogDebug("HitInfo") << " CaloHit " << iHit.getName() << "\n"
                          << " DetID = " << iHit.id() << " EBDetId = " << ebid.ieta() << " " << ebid.iphi() << "\n"
                          << " Time = " << iHit.time() << " Event id. = " << iHit.eventId().rawId() << "\n"
                          << " Track Id = " << iHit.geantTrackId() << "\n"
                          << " Energy = " << iHit.energy();
 
      uint32_t crystid = ebid.rawId();
      ebSimMap[crystid] += iHit.energy();
    }
 
    // loop over Digis
 
    const EBDigiCollection* barrelDigi = EcalDigiEB.product();
 
    std::vector<double> ebAnalogSignal;
    std::vector<double> ebADCCounts;
    std::vector<double> ebADCGains;
    ebAnalogSignal.reserve(EBDataFrame::MAXSAMPLES);
    ebADCCounts.reserve(EBDataFrame::MAXSAMPLES);
    ebADCGains.reserve(EBDataFrame::MAXSAMPLES);
 
    for (unsigned int digis = 0; digis < EcalDigiEB->size(); ++digis) {
      EBDataFrame ebdf = (*barrelDigi)[digis];
      int nrSamples = ebdf.size();
 
      EBDetId ebid = ebdf.id();
 
      double Emax = 0.;
      int Pmax = 0;
      double pedestalPreSample = 0.;
      double pedestalPreSampleAnalog = 0.;
 
      for (int sample = 0; sample < nrSamples; ++sample) {
        ebAnalogSignal[sample] = 0.;
        ebADCCounts[sample] = 0.;
        ebADCGains[sample] = -1.;
      }
 
      for (int sample = 0; sample < nrSamples; ++sample) {
        EcalMGPASample mySample = ebdf[sample];
 
        ebADCCounts[sample] = (mySample.adc());
        ebADCGains[sample] = (mySample.gainId());
        ebAnalogSignal[sample] = (ebADCCounts[sample] * gainConv_[(int)ebADCGains[sample]] * barrelADCtoGeV_);
        if (Emax < ebAnalogSignal[sample]) {
          Emax = ebAnalogSignal[sample];
          Pmax = sample;
        }
        if (sample < 3) {
          pedestalPreSample += ebADCCounts[sample];
          pedestalPreSampleAnalog += ebADCCounts[sample] * gainConv_[(int)ebADCGains[sample]] * barrelADCtoGeV_;
        }
        LogDebug("DigiInfo") << "EB sample " << sample << " ADC counts = " << ebADCCounts[sample]
                             << " Gain Id = " << ebADCGains[sample] << " Analog eq = " << ebAnalogSignal[sample];
      }
 
      pedestalPreSample /= 3.;
      pedestalPreSampleAnalog /= 3.;
      double Erec = Emax - pedestalPreSampleAnalog * gainConv_[(int)ebADCGains[Pmax]];
 
      if (ebSimMap[ebid.rawId()] != 0.) {
        LogDebug("DigiInfo") << " Digi / Hit = " << Erec << " / " << ebSimMap[ebid.rawId()] << " gainConv "
                             << gainConv_[(int)ebADCGains[Pmax]];
        if (meEBDigiSimRatio_)
          meEBDigiSimRatio_->Fill(Erec / ebSimMap[ebid.rawId()]);
        if (Erec > 10. * barrelADCtoGeV_ && meEBDigiSimRatiogt10ADC_)
          meEBDigiSimRatiogt10ADC_->Fill(Erec / ebSimMap[ebid.rawId()]);
        if (Erec > 100. * barrelADCtoGeV_ && meEBDigiSimRatiogt100ADC_)
          meEBDigiSimRatiogt100ADC_->Fill(Erec / ebSimMap[ebid.rawId()]);
      }
    }
  }
 
  // ENDCAP
 
  // loop over simHits
 
  if (isEndcap) {
    e.getByToken(crossingFramePCaloHitEEToken_, crossingFrame);
    const MixCollection<PCaloHit> endcapHits(crossingFrame.product());
 
    MapType eeSimMap;
    for (auto const& iHit : endcapHits) {
      EEDetId eeid = EEDetId(iHit.id());
 
      LogDebug("HitInfo") << " CaloHit " << iHit.getName() << "\n"
                          << " DetID = " << iHit.id() << " EEDetId side = " << eeid.zside() << " = " << eeid.ix() << " "
                          << eeid.iy() << "\n"
                          << " Time = " << iHit.time() << " Event id. = " << iHit.eventId().rawId() << "\n"
                          << " Track Id = " << iHit.geantTrackId() << "\n"
                          << " Energy = " << iHit.energy();
 
      uint32_t crystid = eeid.rawId();
      eeSimMap[crystid] += iHit.energy();
    }
 
    // loop over Digis
 
    const EEDigiCollection* endcapDigi = EcalDigiEE.product();
 
    std::vector<double> eeAnalogSignal;
    std::vector<double> eeADCCounts;
    std::vector<double> eeADCGains;
    eeAnalogSignal.reserve(EEDataFrame::MAXSAMPLES);
    eeADCCounts.reserve(EEDataFrame::MAXSAMPLES);
    eeADCGains.reserve(EEDataFrame::MAXSAMPLES);
 
    for (unsigned int digis = 0; digis < EcalDigiEE->size(); ++digis) {
      EEDataFrame eedf = (*endcapDigi)[digis];
      int nrSamples = eedf.size();
 
      EEDetId eeid = eedf.id();
 
      double Emax = 0.;
      int Pmax = 0;
      double pedestalPreSample = 0.;
      double pedestalPreSampleAnalog = 0.;
 
      for (int sample = 0; sample < nrSamples; ++sample) {
        eeAnalogSignal[sample] = 0.;
        eeADCCounts[sample] = 0.;
        eeADCGains[sample] = -1.;
      }
 
      for (int sample = 0; sample < nrSamples; ++sample) {
        EcalMGPASample mySample = eedf[sample];
 
        eeADCCounts[sample] = (mySample.adc());
        eeADCGains[sample] = (mySample.gainId());
        eeAnalogSignal[sample] = (eeADCCounts[sample] * gainConv_[(int)eeADCGains[sample]] * endcapADCtoGeV_);
        if (Emax < eeAnalogSignal[sample]) {
          Emax = eeAnalogSignal[sample];
          Pmax = sample;
        }
        if (sample < 3) {
          pedestalPreSample += eeADCCounts[sample];
          pedestalPreSampleAnalog += eeADCCounts[sample] * gainConv_[(int)eeADCGains[sample]] * endcapADCtoGeV_;
        }
        LogDebug("DigiInfo") << "EE sample " << sample << " ADC counts = " << eeADCCounts[sample]
                             << " Gain Id = " << eeADCGains[sample] << " Analog eq = " << eeAnalogSignal[sample];
      }
      pedestalPreSample /= 3.;
      pedestalPreSampleAnalog /= 3.;
      double Erec = Emax - pedestalPreSampleAnalog * gainConv_[(int)eeADCGains[Pmax]];
 
      if (eeSimMap[eeid.rawId()] != 0.) {
        LogDebug("DigiInfo") << " Digi / Hit = " << Erec << " / " << eeSimMap[eeid.rawId()] << " gainConv "
                             << gainConv_[(int)eeADCGains[Pmax]];
        if (meEEDigiSimRatio_)
          meEEDigiSimRatio_->Fill(Erec / eeSimMap[eeid.rawId()]);
        if (Erec > 20. * endcapADCtoGeV_ && meEEDigiSimRatiogt20ADC_)
          meEEDigiSimRatiogt20ADC_->Fill(Erec / eeSimMap[eeid.rawId()]);
        if (Erec > 100. * endcapADCtoGeV_ && meEEDigiSimRatiogt100ADC_)
          meEEDigiSimRatiogt100ADC_->Fill(Erec / eeSimMap[eeid.rawId()]);
      }
    }
  }
 
  if (isPreshower) {
    e.getByToken(crossingFramePCaloHitESToken_, crossingFrame);
    const MixCollection<PCaloHit> preshowerHits(crossingFrame.product());
    for (auto const& iHit : preshowerHits) {
      ESDetId esid = ESDetId(iHit.id());
 
      LogDebug("HitInfo") << " CaloHit " << iHit.getName() << "\n"
                          << " DetID = " << iHit.id() << "ESDetId: z side " << esid.zside() << "  plane "
                          << esid.plane() << esid.six() << ',' << esid.siy() << ':' << esid.strip() << "\n"
                          << " Time = " << iHit.time() << " Event id. = " << iHit.eventId().rawId() << "\n"
                          << " Track Id = " << iHit.geantTrackId() << "\n"
                          << " Energy = " << iHit.energy();
    }
  }
}

References funct::abs(), barrelADCtoGeV_, HLT_2018_cff::barrelHits, crossingFramePCaloHitEBToken_, crossingFramePCaloHitEEToken_, crossingFramePCaloHitESToken_, MillePedeFileConverter_cfg::e, EBdigiCollectionToken_, EEdigiCollectionToken_, edm::DataFrameContainer::empty(), endcapADCtoGeV_, HLT_2018_cff::endcapHits, ESdigiCollectionToken_, dqm::impl::MonitorElement::Fill(), g4TkInfoToken_, g4VtxInfoToken_, gainConv_, edm::HepMCProduct::GetEvent(), HepMCToken_, EEDataFrame::id(), EBDataFrame::id(), EBDetId::ieta(), createfilelist::int, EBDetId::iphi(), PixelPluginsPhase0_cfi::isBarrel, GeomDetEnumerators::isEndcap(), edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), dqm-mbProfile::log, LogDebug, M_PI, EcalDataFrame::MAXSAMPLES, meEBDigiSimRatio_, meEBDigiSimRatiogt100ADC_, meEBDigiSimRatiogt10ADC_, meEEDigiSimRatio_, meEEDigiSimRatiogt100ADC_, meEEDigiSimRatiogt20ADC_, meGunEnergy_, meGunEta_, meGunPhi_, ecalTB2006H4_GenSimDigiReco_cfg::mySample, AlCaHLTBitMon_ParallelJobs::p, ESDetId::plane(), dataAnalyzerFineBiningParameters_cff::Pmax, gedPhotons_cfi::preshowerHits, edm::Handle< T >::product(), DetId::rawId(), edm::DataFrameContainer::reserve(), simplePhotonAnalyzer_cfi::sample, ESDetId::six(), ESDetId::siy(), EcalDataFrame::size(), edm::DataFrameContainer::size(), ESDetId::strip(), funct::tan(), ESDetId::zside(), and EEDetId::zside().

bookHistograms()

void EcalDigisValidation::bookHistograms ( DQMStore::IBooker &  i, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 64 of file EcalDigisValidation.cc.

                                                                                                      {
  Char_t histo[200];
 
  ibooker.setCurrentFolder("EcalDigisV/EcalDigiTask");
 
  sprintf(histo, "EcalDigiTask Gun Momentum");
  meGunEnergy_ = ibooker.book1D(histo, histo, 100, 0., 1000.);
 
  sprintf(histo, "EcalDigiTask Gun Eta");
  meGunEta_ = ibooker.book1D(histo, histo, 700, -3.5, 3.5);
 
  sprintf(histo, "EcalDigiTask Gun Phi");
  meGunPhi_ = ibooker.book1D(histo, histo, 360, 0., 360.);
 
  sprintf(histo, "EcalDigiTask Barrel maximum Digi over Sim ratio");
  meEBDigiSimRatio_ = ibooker.book1D(histo, histo, 100, 0., 2.);
 
  sprintf(histo, "EcalDigiTask Endcap maximum Digi over Sim ratio");
  meEEDigiSimRatio_ = ibooker.book1D(histo, histo, 100, 0., 2.);
 
  sprintf(histo, "EcalDigiTask Barrel maximum Digi over Sim ratio gt 10 ADC");
  meEBDigiSimRatiogt10ADC_ = ibooker.book1D(histo, histo, 100, 0., 2.);
 
  sprintf(histo, "EcalDigiTask Endcap maximum Digi over Sim ratio gt 20 ADC");
  meEEDigiSimRatiogt20ADC_ = ibooker.book1D(histo, histo, 100, 0., 2.);
 
  sprintf(histo, "EcalDigiTask Barrel maximum Digi over Sim ratio gt 100 ADC");
  meEBDigiSimRatiogt100ADC_ = ibooker.book1D(histo, histo, 100, 0., 2.);
 
  sprintf(histo, "EcalDigiTask Endcap maximum Digi over Sim ratio gt 100 ADC");
  meEEDigiSimRatiogt100ADC_ = ibooker.book1D(histo, histo, 100, 0., 2.);
}

References dqm::implementation::IBooker::book1D(), timingPdfMaker::histo, meEBDigiSimRatio_, meEBDigiSimRatiogt100ADC_, meEBDigiSimRatiogt10ADC_, meEEDigiSimRatio_, meEEDigiSimRatiogt100ADC_, meEEDigiSimRatiogt20ADC_, meGunEnergy_, meGunEta_, meGunPhi_, and dqm::implementation::NavigatorBase::setCurrentFolder().

checkCalibrations()

void EcalDigisValidation::checkCalibrations ( edm::EventSetup const &  c )

private

Definition at line 383 of file EcalDigisValidation.cc.

                                                                           {
  // ADC -> GeV Scale
  edm::ESHandle<EcalADCToGeVConstant> pAgc;
  eventSetup.get<EcalADCToGeVConstantRcd>().get(pAgc);
  const EcalADCToGeVConstant* agc = pAgc.product();
 
  EcalMGPAGainRatio* defaultRatios = new EcalMGPAGainRatio();
 
  gainConv_[1] = 1.;
  gainConv_[2] = defaultRatios->gain12Over6();
  gainConv_[3] = gainConv_[2] * (defaultRatios->gain6Over1());
  gainConv_[0] = gainConv_[2] * (defaultRatios->gain6Over1());  // saturated channels
 
  LogDebug("EcalDigi") << " Gains conversions: "
                       << "\n"
                       << " g1 = " << gainConv_[1] << "\n"
                       << " g2 = " << gainConv_[2] << "\n"
                       << " g3 = " << gainConv_[3];
  LogDebug("EcalDigi") << " Gains conversions: "
                       << "\n"
                       << " saturation = " << gainConv_[0];
 
  delete defaultRatios;
 
  const double barrelADCtoGeV_ = agc->getEBValue();
  LogDebug("EcalDigi") << " Barrel GeV/ADC = " << barrelADCtoGeV_;
  const double endcapADCtoGeV_ = agc->getEEValue();
  LogDebug("EcalDigi") << " Endcap GeV/ADC = " << endcapADCtoGeV_;
}

References barrelADCtoGeV_, endcapADCtoGeV_, EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), gainConv_, edm::EventSetup::get(), get, EcalADCToGeVConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), LogDebug, and edm::ESHandle< T >::product().

Referenced by dqmBeginRun().

dqmBeginRun()

void EcalDigisValidation::dqmBeginRun ( edm::Run const &  , edm::EventSetup const &  c  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 62 of file EcalDigisValidation.cc.

{ checkCalibrations(c); }

References HltBtagPostValidation_cff::c, and checkCalibrations().

Member Data Documentation

barrelADCtoGeV_

double EcalDigisValidation::barrelADCtoGeV_

private

Definition at line 84 of file EcalDigisValidation.h.

Referenced by analyze(), checkCalibrations(), and EcalDigisValidation().

crossingFramePCaloHitEBToken_

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalDigisValidation::crossingFramePCaloHitEBToken_

private

Definition at line 79 of file EcalDigisValidation.h.

Referenced by analyze().

crossingFramePCaloHitEEToken_

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalDigisValidation::crossingFramePCaloHitEEToken_

private

Definition at line 79 of file EcalDigisValidation.h.

Referenced by analyze().

crossingFramePCaloHitESToken_

edm::EDGetTokenT<CrossingFrame<PCaloHit> > EcalDigisValidation::crossingFramePCaloHitESToken_

private

Definition at line 79 of file EcalDigisValidation.h.

Referenced by analyze().

EBdigiCollectionToken_

edm::EDGetTokenT<EBDigiCollection> EcalDigisValidation::EBdigiCollectionToken_

private

Definition at line 75 of file EcalDigisValidation.h.

Referenced by analyze().

EEdigiCollectionToken_

edm::EDGetTokenT<EEDigiCollection> EcalDigisValidation::EEdigiCollectionToken_

private

Definition at line 76 of file EcalDigisValidation.h.

Referenced by analyze().

endcapADCtoGeV_

double EcalDigisValidation::endcapADCtoGeV_

private

Definition at line 85 of file EcalDigisValidation.h.

Referenced by analyze(), checkCalibrations(), and EcalDigisValidation().

ESdigiCollectionToken_

edm::EDGetTokenT<ESDigiCollection> EcalDigisValidation::ESdigiCollectionToken_

private

Definition at line 77 of file EcalDigisValidation.h.

Referenced by analyze().

g4TkInfoToken_

edm::EDGetTokenT<edm::SimTrackContainer> EcalDigisValidation::g4TkInfoToken_

private

Definition at line 72 of file EcalDigisValidation.h.

Referenced by analyze().

g4VtxInfoToken_

edm::EDGetTokenT<edm::SimVertexContainer> EcalDigisValidation::g4VtxInfoToken_

private

Definition at line 73 of file EcalDigisValidation.h.

Referenced by analyze().

gainConv_

std::map<int, double, std::less<int> > EcalDigisValidation::gainConv_

private

Definition at line 82 of file EcalDigisValidation.h.

Referenced by analyze(), checkCalibrations(), and EcalDigisValidation().

HepMCToken_

edm::EDGetTokenT<edm::HepMCProduct> EcalDigisValidation::HepMCToken_

private

Definition at line 71 of file EcalDigisValidation.h.

Referenced by analyze().

meEBDigiSimRatio_

MonitorElement* EcalDigisValidation::meEBDigiSimRatio_

private

Definition at line 91 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

meEBDigiSimRatiogt100ADC_

MonitorElement* EcalDigisValidation::meEBDigiSimRatiogt100ADC_

private

Definition at line 97 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

meEBDigiSimRatiogt10ADC_

MonitorElement* EcalDigisValidation::meEBDigiSimRatiogt10ADC_

private

Definition at line 94 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

meEEDigiSimRatio_

MonitorElement* EcalDigisValidation::meEEDigiSimRatio_

private

Definition at line 92 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

meEEDigiSimRatiogt100ADC_

MonitorElement* EcalDigisValidation::meEEDigiSimRatiogt100ADC_

private

Definition at line 98 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

meEEDigiSimRatiogt20ADC_

MonitorElement* EcalDigisValidation::meEEDigiSimRatiogt20ADC_

private

Definition at line 95 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

meGunEnergy_

MonitorElement* EcalDigisValidation::meGunEnergy_

private

Definition at line 87 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

meGunEta_

MonitorElement* EcalDigisValidation::meGunEta_

private

Definition at line 88 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

meGunPhi_

MonitorElement* EcalDigisValidation::meGunPhi_

private

Definition at line 89 of file EcalDigisValidation.h.

Referenced by analyze(), bookHistograms(), and EcalDigisValidation().

outputFile_

std::string EcalDigisValidation::outputFile_

private

Definition at line 69 of file EcalDigisValidation.h.

Referenced by EcalDigisValidation().

verbose_

bool EcalDigisValidation::verbose_

private

Definition at line 67 of file EcalDigisValidation.h.

Referenced by EcalDigisValidation().