EcalUncalibRecHitPhase2WeightsProducer Class Reference

Inheritance diagram for EcalUncalibRecHitPhase2WeightsProducer:

Public Member Functions
	EcalUncalibRecHitPhase2WeightsProducer (const edm::ParameterSet &ps)
void	produce (edm::Event &evt, const edm::EventSetup &es) override
Public Member Functions inherited from edm::stream::EDProducer<>
	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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
const edm::EDGetTokenT< EBDigiCollectionPh2 >	ebDigiCollectionToken_
const edm::EDPutTokenT< EBUncalibratedRecHitCollection >	ebUncalibRecHitCollectionToken_
const float	tFall_
const float	tRise_
const std::vector< double >	weights_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
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

Detailed Description

Definition at line 15 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

Constructor & Destructor Documentation

EcalUncalibRecHitPhase2WeightsProducer()

EcalUncalibRecHitPhase2WeightsProducer::EcalUncalibRecHitPhase2WeightsProducer ( const edm::ParameterSet &  ps )

explicit

Definition at line 60 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

    : tRise_(ps.getParameter<double>("tRise")),
      tFall_(ps.getParameter<double>("tFall")),
      weights_(ps.getParameter<std::vector<double>>("weights")),
      ebDigiCollectionToken_(consumes<EBDigiCollectionPh2>(ps.getParameter<edm::InputTag>("BarrelDigis"))),
      ebUncalibRecHitCollectionToken_(
          produces<EBUncalibratedRecHitCollection>(ps.getParameter<std::string>("EBhitCollection"))) {}

Member Function Documentation

fillDescriptions()

void EcalUncalibRecHitPhase2WeightsProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 30 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, HLT_2022v14_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                        {
  edm::ParameterSetDescription desc;

  desc.add<std::string>("EBhitCollection", "EcalUncalibRecHitsEB");
  desc.add<double>("tRise", 0.2);
  desc.add<double>("tFall", 2.);
  // The below weights values should be kept up to date with GPU version of this module
  desc.add<std::vector<double>>("weights",
                                {-0.121016,
                                 -0.119899,
                                 -0.120923,
                                 -0.0848959,
                                 0.261041,
                                 0.509881,
                                 0.373591,
                                 0.134899,
                                 -0.0233605,
                                 -0.0913195,
                                 -0.112452,
                                 -0.118596,
                                 -0.121737,
                                 -0.121737,
                                 -0.121737,
                                 -0.121737});

  desc.add<edm::InputTag>("BarrelDigis", edm::InputTag("simEcalUnsuppressedDigis", ""));

  descriptions.addWithDefaultLabel(desc);
}

produce()

void EcalUncalibRecHitPhase2WeightsProducer::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

override

Definition at line 68 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

References EcalLiteDTUSample::adc(), edm::DataFrameContainer::begin(), ebDigiCollectionToken_, ebUncalibRecHitCollectionToken_, edm::DataFrameContainer::end(), diffTwoXMLs::g1, EcalLiteDTUSample::gainId(), ecalPh2::gains, edm::Event::get(), EBDataFrame::id(), EcalUncalibratedRecHit::kHasSwitchToGain1, eostools::move(), PresampleTask_cfi::nSamples, edm::Event::put(), ecalGpuTask_cfi::sample, EcalUncalibratedRecHit::setAmplitudeError(), EcalUncalibratedRecHit::setFlagBit(), EcalUncalibratedRecHit::setJitterError(), EcalDataFrame::size(), and weights_.

                                                                                           {
  // retrieve digis
  const EBDigiCollectionPh2* pdigis = &evt.get(ebDigiCollectionToken_);

  // prepare output
  auto ebUncalibRechits = std::make_unique<EBUncalibratedRecHitCollection>();

  for (auto itdg = pdigis->begin(); itdg != pdigis->end(); ++itdg) {
    EBDataFrame digi(*itdg);
    EcalDataFrame_Ph2 dataFrame(*itdg);
    DetId detId(digi.id());

    bool g1 = false;

    std::vector<float> timetrace;
    std::vector<float> adctrace;

    int nSamples = digi.size();

    timetrace.reserve(nSamples);
    adctrace.reserve(nSamples);

    float amp = 0;
    float gratio;

    for (int sample = 0; sample < nSamples; ++sample) {
      EcalLiteDTUSample thisSample = dataFrame[sample];
      gratio = ecalPh2::gains[thisSample.gainId()];
      adctrace.push_back(thisSample.adc() * gratio);

      amp = amp + adctrace[sample] * weights_[sample];

      if (thisSample.gainId() == 1)
        g1 = true;

      timetrace.push_back(sample);

    }  // loop on samples

    float amp_e = 1;
    float t0_e = 0;

    EcalUncalibratedRecHit rhit(detId, amp, 0., 0., 0., 0);  // rhit(detIt, amp, pedestal, t0, chi2, flags)
    rhit.setAmplitudeError(amp_e);
    rhit.setJitterError(t0_e);
    if (g1)
      rhit.setFlagBit(EcalUncalibratedRecHit::kHasSwitchToGain1);

    ebUncalibRechits->push_back(rhit);

  }  // loop on digis

  evt.put(ebUncalibRecHitCollectionToken_, std::move(ebUncalibRechits));
}

Member Data Documentation

ebDigiCollectionToken_

const edm::EDGetTokenT<EBDigiCollectionPh2> EcalUncalibRecHitPhase2WeightsProducer::ebDigiCollectionToken_

private

Definition at line 26 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

Referenced by produce().

ebUncalibRecHitCollectionToken_

const edm::EDPutTokenT<EBUncalibratedRecHitCollection> EcalUncalibRecHitPhase2WeightsProducer::ebUncalibRecHitCollectionToken_

private

Definition at line 27 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

Referenced by produce().

tFall_

const float EcalUncalibRecHitPhase2WeightsProducer::tFall_

private

Definition at line 23 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

tRise_

const float EcalUncalibRecHitPhase2WeightsProducer::tRise_

private

Definition at line 22 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

weights_

const std::vector<double> EcalUncalibRecHitPhase2WeightsProducer::weights_

private

Definition at line 24 of file EcalUncalibRecHitPhase2WeightsProducer.cc.

Referenced by produce().