Phase1L1TJetCalibrator Class Reference

#include <Phase1L1TJetCalibrator.cc>

Inheritance diagram for Phase1L1TJetCalibrator:

Public Member Functions
	Phase1L1TJetCalibrator (const edm::ParameterSet &)
	~Phase1L1TJetCalibrator () 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 Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::vector< std::vector < double > >	_jetCalibrationFactorsPtBins
std::vector< double >	absEtaBinning_
std::vector< edm::ParameterSet >	calibration_
edm::EDGetTokenT< std::vector < reco::CaloJet > >	inputCollectionTag_
std::vector< std::vector < double > >	jetCalibrationFactorsBinnedInEta_
size_t	nBinsEta_
std::string	outputCollectionName_

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

Description: Applies calibrations to reco::calojets

INPUT PARAMETERS ***
inputCollectionTag, InputTag, collection of reco calojet to calibrate
absEtaBinning, vdouble with eta bins, allows for non-homogeneous binning
calibration, VPSet with calibration factors
outputCollectionName, string, output collection tag name

Definition at line 43 of file Phase1L1TJetCalibrator.cc.

Constructor & Destructor Documentation

Phase1L1TJetCalibrator::Phase1L1TJetCalibrator ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 76 of file Phase1L1TJetCalibrator.cc.

References _jetCalibrationFactorsPtBins, calibration_, jetCalibrationFactorsBinnedInEta_, outputCollectionName_, and TrackValidation_cff::pset.

    : inputCollectionTag_(edm::EDGetTokenT<std::vector<reco::CaloJet>>(
          consumes<std::vector<reco::CaloJet>>(iConfig.getParameter<edm::InputTag>("inputCollectionTag")))),
      absEtaBinning_(iConfig.getParameter<std::vector<double>>("absEtaBinning")),
      nBinsEta_(absEtaBinning_.size() - 1),
      calibration_(iConfig.getParameter<std::vector<edm::ParameterSet>>("calibration")),
      outputCollectionName_(iConfig.getParameter<std::string>("outputCollectionName"))

{
  for (const auto& pset : calibration_) {
    _jetCalibrationFactorsPtBins.emplace_back(pset.getParameter<std::vector<double>>("l1tPtBins"));
    jetCalibrationFactorsBinnedInEta_.emplace_back(pset.getParameter<std::vector<double>>("l1tCalibrationFactors"));
  }

  produces<std::vector<reco::CaloJet>>(outputCollectionName_).setBranchAlias(outputCollectionName_);
}

Phase1L1TJetCalibrator::~Phase1L1TJetCalibrator ( )

override

Definition at line 93 of file Phase1L1TJetCalibrator.cc.

{}

Member Function Documentation

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

static

Definition at line 169 of file Phase1L1TJetCalibrator.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addVPSet(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                        {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("inputCollectionTag",
                          edm::InputTag("Phase1L1TJetProducer", "UncalibratedPhase1L1TJetFromPfCandidates"));
  desc.add<std::vector<double>>("absEtaBinning");
  std::vector<edm::ParameterSet> vDefaults;
  edm::ParameterSetDescription validator;
  validator.add<double>("etaMax");
  validator.add<double>("etaMin");
  validator.add<std::vector<double>>("l1tCalibrationFactors");
  validator.add<std::vector<double>>("l1tPtBins");
  desc.addVPSet("calibration", validator, vDefaults);
  desc.add<std::string>("outputCollectionName", "Phase1L1TJetFromPfCandidates");
  descriptions.add("Phase1L1TJetCalibrator", desc);
}

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

overrideprivate

Definition at line 100 of file Phase1L1TJetCalibrator.cc.

References _jetCalibrationFactorsPtBins, absEtaBinning_, PVValHelper::eta, etaBin(), edm::Event::getByToken(), inputCollectionTag_, jetCalibrationFactorsBinnedInEta_, LogDebug, eostools::move(), outputCollectionName_, DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), edm::Event::put(), and cuda_std::upper_bound().

                                                                                  {
  edm::Handle<std::vector<reco::CaloJet>> inputCollectionHandle;
  iEvent.getByToken(inputCollectionTag_, inputCollectionHandle);

  auto calibratedCollectionPtr = std::make_unique<std::vector<reco::CaloJet>>();

  // for each candidate:
  // 1 get pt and eta
  // 2 run a dicotomic search on the eta vector to find the eta index
  // 3 fetch the corresponding calibration elements
  // 4 run a dicotomic search on the pt binning to find the pt index
  // 5 fetch the calibration factor
  // 6 update the candidate pt by applying the calibration factor
  // 7 store calibrated candidate in a new collection

  calibratedCollectionPtr->reserve(inputCollectionHandle->size());

  for (const auto& candidate : *inputCollectionHandle) {
    // 1
    float pt = candidate.pt();
    float eta = candidate.eta();

    //2
    auto etaBin = upper_bound(absEtaBinning_.begin(), absEtaBinning_.end(), fabs(eta));
    int etaIndex = etaBin - absEtaBinning_.begin() - 1;

    //3
    const std::vector<double>& l1tPtBins = _jetCalibrationFactorsPtBins[etaIndex];
    const std::vector<double>& l1tCalibrationFactors = jetCalibrationFactorsBinnedInEta_[etaIndex];

    //4
    auto ptBin = upper_bound(l1tPtBins.begin(), l1tPtBins.end(), pt);
    int ptBinIndex = ptBin - l1tPtBins.begin() - 1;

    //5
    const double& l1tCalibrationFactor = l1tCalibrationFactors[ptBinIndex];

    //6 and 7
    reco::Candidate::PolarLorentzVector candidateP4(candidate.polarP4());
    candidateP4.SetPt(candidateP4.pt() * l1tCalibrationFactor);
    calibratedCollectionPtr->emplace_back(candidate);
    calibratedCollectionPtr->back().setP4(candidateP4);

#ifdef DEBUG
    if (newCandidate->pt() < 0) {
      LogDebug("Phase1L1TJetCalibrator") << "######################" << std::endl;
      LogDebug("Phase1L1TJetCalibrator") << "PRE-CALIBRATION " << std::endl;
      LogDebug("Phase1L1TJetCalibrator") << "\t Jet properties (pt, eta, phi, pile-up): " << candidate.pt() << "\t"
                                         << candidate.eta() << "\t" LogDebug("Phase1L1TJetCalibrator")
                                         << candidate.phi() << "\t" << candidate.pileup() << std::endl;
      LogDebug("Phase1L1TJetCalibrator") << "CALIBRATION " << std::endl;
      LogDebug("Phase1L1TJetCalibrator") << "\t Using eta - pt - factor " << *etaBin << " - " << *ptBin << " - "
                                         << l1tCalibrationFactor LogDebug("Phase1L1TJetCalibrator") << std::endl;
      LogDebug("Phase1L1TJetCalibrator") << "POST-CALIBRATION " << std::endl;
      LogDebug("Phase1L1TJetCalibrator") << "\t Jet properties (pt, eta, phi, pile-up): " << newCandidate->pt() << "\t"
                                         << newCandidate->eta() << "\t" << newCandidate->phi() << "\t"
                                         << newCandidate->pileup() << std::endl;
    }
#endif
  }

  // finally, sort the collection by pt
  std::sort(calibratedCollectionPtr->begin(),
            calibratedCollectionPtr->end(),
            [](const reco::CaloJet& jet1, const reco::CaloJet& jet2) { return jet1.pt() > jet2.pt(); });

  iEvent.put(std::move(calibratedCollectionPtr), outputCollectionName_);
}

Member Data Documentation

std::vector<std::vector<double> > Phase1L1TJetCalibrator::_jetCalibrationFactorsPtBins

private

Definition at line 60 of file Phase1L1TJetCalibrator.cc.

Referenced by Phase1L1TJetCalibrator(), and produce().

std::vector<double> Phase1L1TJetCalibrator::absEtaBinning_

private

Definition at line 54 of file Phase1L1TJetCalibrator.cc.

Referenced by produce().

std::vector<edm::ParameterSet> Phase1L1TJetCalibrator::calibration_

private

Definition at line 56 of file Phase1L1TJetCalibrator.cc.

Referenced by Phase1L1TJetCalibrator().

edm::EDGetTokenT<std::vector<reco::CaloJet> > Phase1L1TJetCalibrator::inputCollectionTag_

private

Definition at line 53 of file Phase1L1TJetCalibrator.cc.

Referenced by produce().

std::vector<std::vector<double> > Phase1L1TJetCalibrator::jetCalibrationFactorsBinnedInEta_

private

Definition at line 59 of file Phase1L1TJetCalibrator.cc.

Referenced by Phase1L1TJetCalibrator(), and produce().

size_t Phase1L1TJetCalibrator::nBinsEta_

private

Definition at line 55 of file Phase1L1TJetCalibrator.cc.

std::string Phase1L1TJetCalibrator::outputCollectionName_

private

Definition at line 57 of file Phase1L1TJetCalibrator.cc.

Referenced by Phase1L1TJetCalibrator(), and produce().