L1TZDCProducer Class Reference

Inheritance diagram for L1TZDCProducer:

Public Member Functions
	L1TZDCProducer (const edm::ParameterSet &ps)
	~L1TZDCProducer () override=default
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 Member Functions
void	beginRun (edm::Run const &, edm::EventSetup const &) override
void	produce (edm::Event &, const edm::EventSetup &) override
int	zdcLUTIndexHelper (int iDetPos, int iBXPos)

Private Attributes
int	bxFirst_
int	bxLast_
edm::ESGetToken< CaloParams, L1TCaloParamsRcd >	candidateToken_
edm::EDPutTokenT< EtSumBxCollection >	etToken_
std::unique_ptr< CaloParamsHelper >	params_
int	sampleToCenterBX_
int	scaleFactor_
edm::EDGetTokenT< QIE10DigiCollection >	zdcToken_

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 62 of file L1TZDCProducer.cc.

Constructor & Destructor Documentation

L1TZDCProducer()

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

explicit

Definition at line 97 of file L1TZDCProducer.cc.

References edm::ParameterSet::getParameter().

                                                        {
  // register what you produce
  etToken_ = produces<EtSumBxCollection>();

  // register what you consume and keep token for later access:
  zdcToken_ = consumes<QIE10DigiCollection>(ps.getParameter<edm::InputTag>("zdcDigis"));
  candidateToken_ = esConsumes<CaloParams, L1TCaloParamsRcd, edm::Transition::BeginRun>();

  bxFirst_ = ps.getParameter<int>("bxFirst");
  bxLast_ = ps.getParameter<int>("bxLast");
  sampleToCenterBX_ = ps.getParameter<int>("sampleToCenterBX");
}

~L1TZDCProducer()

L1TZDCProducer::~L1TZDCProducer ( )

overridedefault

Member Function Documentation

beginRun()

void L1TZDCProducer::beginRun ( edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overrideprivate

Definition at line 237 of file L1TZDCProducer.cc.

References edm::EventSetup::getHandle(), and edm::ESHandle< T >::product().

                                                                             {
  edm::ESHandle<CaloParams> candidateHandle = iSetup.getHandle(candidateToken_);

  params_ = std::make_unique<CaloParamsHelper>(*candidateHandle.product());
  scaleFactor_ = params_->zdcLUT()->data(0);  //First position is the integer scaling factor
  edm::LogInfo("L1TZDCProducer") << "SCALE FACTOR FOR LUT: " << scaleFactor_ << std::endl;
}

fillDescriptions()

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

static

Definition at line 249 of file L1TZDCProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("zdcDigis", edm::InputTag("hcalDigis", "ZDC"));
  desc.add<int>("bxFirst", -2);
  desc.add<int>("bxLast", 2);
  desc.add<int>("sampleToCenterBX", 4);
  descriptions.add("l1tZDCProducer", desc);
}

produce()

void L1TZDCProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 111 of file L1TZDCProducer.cc.

References gpuClustering::adc, edm::DataFrameContainer::begin(), HcalZDCDetId::channel(), edm::DataFrameContainer::end(), l1t::CaloTools::etSumP4Demux(), amptDefault_cfi::frame, iEvent, edm::HandleBase::isValid(), L1Analysis::kZDCM, L1Analysis::kZDCP, LogDebug, eostools::move(), PresampleTask_cfi::nSamples, BXVector< T >::push_back(), QIE10DataFrame::samples(), HcalZDCDetId::section(), hgcalPlots::section, l1t::L1Candidate::setHwEta(), l1t::L1Candidate::setHwPhi(), l1t::L1Candidate::setHwPt(), l1t::EtSum::setType(), HcalZDCDetId::zside(), and ecaldqm::zside().

                                                                          {
  using namespace edm;
  using namespace l1t;

  LogDebug("L1TZDCProducer") << "L1TZDCProducer::produce function called..." << std::endl;

  // reduced collection to be emplaced in output
  EtSumBxCollection etsumsReduced(0, bxFirst_, bxLast_);

  //inputs
  Handle<QIE10DigiCollection> zdcDigiCollection;
  iEvent.getByToken(zdcToken_, zdcDigiCollection);

  //Produce ZDC EtSums IF input zdcDigiCollection is valid
  if (zdcDigiCollection.isValid()) {
    //In lieu of bxFirst, bxLast, use the number of the timeslice samples
    const QIE10DataFrame& frametest = (*zdcDigiCollection)[0];
    int nSamples = frametest.samples();

    //Check if the sample requested as center bx is valid
    int sampleToCenterBXChecked = sampleToCenterBX_;

    if (sampleToCenterBXChecked < 0) {
      edm::LogWarning("L1TZDCProducer") << "sampleToCenterBX LT 0; Set bx 0 to sample 0 (minimum allowed)" << std::endl;
      sampleToCenterBXChecked = 0;
    } else if (sampleToCenterBXChecked >= nSamples) {
      edm::LogWarning("L1TZDCProducer")
          << "sampleToCenterBX GE nsamples; Set bx 0 to sample nsamples-1 (maximum allowed)" << std::endl;
      sampleToCenterBXChecked = nSamples - 1;
    }

    //rawadc[detector index][time slices]
    unsigned short rawadc[18][10];

    // the loop below loops over all the elements of the QIE10DigiCollection. Each entry corresponds to one channel
    for (QIE10DigiCollection::const_iterator it = zdcDigiCollection->begin(); it != zdcDigiCollection->end(); it++) {
      const QIE10DataFrame& frame(*it);
      HcalZDCDetId cell = frame.id();
      int zside = cell.zside();
      int section = cell.section();
      int channel = cell.channel();

      if (zside != -1 && zside != 1)
        continue;
      if (section != 1 && section != 2)
        continue;
      if (section == 1 && (channel < 1 || channel > 5))
        continue;
      if (section == 2 && (channel < 1 || channel > 4))
        continue;

      int ihitid = (zside == 1 ? 9 : 0) + (section == 2 ? 5 : 0) + (channel - 1);
      //the loop below iterates over the time slices
      for (int iTS = 0; iTS < nSamples; iTS++) {
        unsigned short adc = (unsigned short)frame[iTS].adc();
        rawadc[ihitid][iTS] = adc;
      }  // end of loop over iTS
    }    //end of loop over channels

    for (int ibx = 0; ibx < nSamples; ibx++) {
      double cEMP = 0, cEMM = 0, cHDP = 0, cHDM = 0;
      double sumcEMP = 0, sumcEMM = 0, sumcHDP = 0, sumcHDM = 0;
      //idet=0-4 correpond to the EM channels
      for (int idet = 0; idet < 5; idet++) {
        unsigned short EMP = rawadc[idet + 9][ibx];
        unsigned short EMM = rawadc[idet][ibx];

        int cEMP_LUTIndex = zdcLUTIndexHelper(idet + 9, (int)EMP);
        int cEMM_LUTIndex = zdcLUTIndexHelper(idet, (int)EMM);

        cEMP = ((double)params_->zdcLUT()->data(cEMP_LUTIndex)) / ((double)scaleFactor_);
        cEMM = ((double)params_->zdcLUT()->data(cEMM_LUTIndex)) / ((double)scaleFactor_);

        sumcEMP = sumcEMP + cEMP;
        sumcEMM = sumcEMM + cEMM;
      }
      //idet=5-8 correspond to HAD channels
      for (int idet = 5; idet < 9; idet++) {
        unsigned short HDP = rawadc[idet + 9][ibx];
        unsigned short HDM = rawadc[idet][ibx];

        int cHDP_LUTIndex = zdcLUTIndexHelper(idet + 9, (int)HDP);
        int cHDM_LUTIndex = zdcLUTIndexHelper(idet, (int)HDM);

        cHDP = ((double)params_->zdcLUT()->data(cHDP_LUTIndex)) / ((double)scaleFactor_);
        cHDM = ((double)params_->zdcLUT()->data(cHDM_LUTIndex)) / ((double)scaleFactor_);

        sumcHDP = sumcHDP + cHDP;
        sumcHDM = sumcHDM + cHDM;
      }
      double sumM = sumcEMM + sumcHDM;
      double sumP = sumcEMP + sumcHDP;

      if (ibx == 4) {
        edm::LogInfo("L1TZDCProducer") << ", sumM= " << sumM << std::endl;
        edm::LogInfo("L1TZDCProducer") << ", sumP= " << sumP << std::endl;
      }
      l1t::EtSum tempEtM = l1t::EtSum();
      tempEtM.setHwPt(sumM);
      tempEtM.setHwEta(-1.);
      tempEtM.setHwPhi(0.);
      tempEtM.setType(EtSum::EtSumType::kZDCM);

      l1t::EtSum tempEtP = l1t::EtSum();
      tempEtP.setHwPt(sumP);
      tempEtP.setHwEta(1.);
      tempEtP.setHwPhi(0.);
      tempEtP.setType(EtSum::EtSumType::kZDCP);

      if (ibx >= sampleToCenterBXChecked + bxFirst_ && ibx <= sampleToCenterBXChecked + bxLast_) {
        etsumsReduced.push_back(ibx - sampleToCenterBXChecked, CaloTools::etSumP4Demux(tempEtP));
        etsumsReduced.push_back(ibx - sampleToCenterBXChecked, CaloTools::etSumP4Demux(tempEtM));
      }
    }  // end of loop over bunch crossings
  }    // end if(zdcDigiCollection.isValid())
  else {
    // If the collection is not valid issue a warning before putting an empty collection
    edm::LogWarning("L1TZDCProducer") << "zdcDigis not valid; return empty ZDC Et Sum BXCollection" << std::endl;
  }

  // Emplace even if !zdcDigiCollection.isValid()
  // Output in this case will be an empty collection
  iEvent.emplace(etToken_, std::move(etsumsReduced));
}

zdcLUTIndexHelper()

int L1TZDCProducer::zdcLUTIndexHelper ( int  iDetPos, int  iBXPos  )

private

Definition at line 246 of file L1TZDCProducer.cc.

{ return 1 + iDetPos * 256 + iBxPos; }

Member Data Documentation

bxFirst_

int L1TZDCProducer::bxFirst_

private

Definition at line 85 of file L1TZDCProducer.cc.

bxLast_

int L1TZDCProducer::bxLast_

private

Definition at line 86 of file L1TZDCProducer.cc.

candidateToken_

edm::ESGetToken<CaloParams, L1TCaloParamsRcd> L1TZDCProducer::candidateToken_

private

Definition at line 82 of file L1TZDCProducer.cc.

etToken_

edm::EDPutTokenT<EtSumBxCollection> L1TZDCProducer::etToken_

private

Definition at line 90 of file L1TZDCProducer.cc.

params_

std::unique_ptr<CaloParamsHelper> L1TZDCProducer::params_

private

Definition at line 93 of file L1TZDCProducer.cc.

sampleToCenterBX_

int L1TZDCProducer::sampleToCenterBX_

private

Definition at line 87 of file L1TZDCProducer.cc.

scaleFactor_

int L1TZDCProducer::scaleFactor_

private

Definition at line 94 of file L1TZDCProducer.cc.

zdcToken_

edm::EDGetTokenT<QIE10DigiCollection> L1TZDCProducer::zdcToken_

private

Definition at line 81 of file L1TZDCProducer.cc.