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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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<>
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

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::Phase1L1TJetCalibrator ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 76 of file Phase1L1TJetCalibrator.cc.

    : 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_);
}

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

~Phase1L1TJetCalibrator()

Phase1L1TJetCalibrator::~Phase1L1TJetCalibrator ( )

override

Definition at line 93 of file Phase1L1TJetCalibrator.cc.

{}

Member Function Documentation

fillDescriptions()

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

static

Definition at line 169 of file Phase1L1TJetCalibrator.cc.

                                                                                        {
  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);
}

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addVPSet(), HLT_2018_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

produce()

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

overrideprivate

Definition at line 100 of file Phase1L1TJetCalibrator.cc.

                                                                                  {
  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_);
}

References _jetCalibrationFactorsPtBins, absEtaBinning_, PVValHelper::eta, etaBin(), iEvent, inputCollectionTag_, jetCalibrationFactorsBinnedInEta_, Phase1L1TJetCalibrator_cfi::l1tCalibrationFactors, Phase1L1TJetCalibrator_cfi::l1tPtBins, LogDebug, eostools::move(), outputCollectionName_, DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), muonRecoAnalyzer_cfi::ptBin, and cuda_std::upper_bound().

Member Data Documentation

_jetCalibrationFactorsPtBins

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

private

Definition at line 60 of file Phase1L1TJetCalibrator.cc.

Referenced by Phase1L1TJetCalibrator(), and produce().

absEtaBinning_

std::vector<double> Phase1L1TJetCalibrator::absEtaBinning_

private

Definition at line 54 of file Phase1L1TJetCalibrator.cc.

Referenced by produce().

calibration_

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

private

Definition at line 56 of file Phase1L1TJetCalibrator.cc.

Referenced by Phase1L1TJetCalibrator().

inputCollectionTag_

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

private

Definition at line 53 of file Phase1L1TJetCalibrator.cc.

Referenced by produce().

jetCalibrationFactorsBinnedInEta_

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

private

Definition at line 59 of file Phase1L1TJetCalibrator.cc.

Referenced by Phase1L1TJetCalibrator(), and produce().

nBinsEta_

size_t Phase1L1TJetCalibrator::nBinsEta_

private

Definition at line 55 of file Phase1L1TJetCalibrator.cc.

outputCollectionName_

std::string Phase1L1TJetCalibrator::outputCollectionName_

private

Definition at line 57 of file Phase1L1TJetCalibrator.cc.

Referenced by Phase1L1TJetCalibrator(), and produce().