GetLumi.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author:  Mia Tosi,40 3-B32,+41227671609 
 */

#include <iostream>
#include "DQM/TrackingMonitor/interface/GetLumi.h"

#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "DataFormats/Luminosity/interface/LumiDetails.h"
#include "DataFormats/Luminosity/interface/LumiSummary.h"

GetLumi::GetLumi(const edm::ParameterSet& iConfig)
    : lumiInputTag_(iConfig.getParameter<edm::InputTag>("lumi")),
      lumiScale_(iConfig.getParameter<double>("lumiScale")) {}

GetLumi::GetLumi(const edm::InputTag& lumiInputTag, double lumiScale)
    : lumiInputTag_(lumiInputTag), lumiScale_(lumiScale) {}

GetLumi::GetLumi(const edm::ParameterSet& iConfig, edm::ConsumesCollector& iC) : GetLumi(iConfig) {
  lumiDetailsToken_ = iC.consumes<LumiDetails>(lumiInputTag_);
  lumiSummaryToken_ = iC.consumes<LumiSummary, edm::InLumi>(lumiInputTag_);
}

GetLumi::GetLumi(const edm::InputTag& lumiInputTag, double lumiScale, edm::ConsumesCollector& iC)
    : GetLumi(lumiInputTag, lumiScale) {
  lumiDetailsToken_ = iC.consumes<LumiDetails>(lumiInputTag_);
  lumiSummaryToken_ = iC.consumes<LumiSummary, edm::InLumi>(lumiInputTag_);
}

GetLumi::~GetLumi() {}

double GetLumi::getRawValue(const edm::Event& iEvent) {
  // taken from
  // DPGAnalysis/SiStripTools/src/DigiLumiCorrHistogramMaker.cc
  // the scale factor 6.37 should follow the lumi prescriptions
  edm::Handle<LumiDetails> lumi;
  iEvent.getLuminosityBlock().getByToken(lumiSummaryToken_, lumi);

  double bxlumi = -1.;
  if (lumi->isValid()) {
    bxlumi = lumi->lumiValue(LumiDetails::kOCC1, iEvent.bunchCrossing());
  }

  return bxlumi;
}

double GetLumi::getValue(const edm::Event& iEvent) {
  //    bxlumi = lumi->lumiValue(LumiDetails::kOCC1,iEvent.bunchCrossing())*6.37;
  return getRawValue(iEvent) * lumiScale_;
}

double GetLumi::getRawValue(edm::LuminosityBlock const& lumiBlock, edm::EventSetup const& eSetup) {
  double lumi = -1.;
  double intDelLumi = -1.;

  //  size_t LS = lumiBlock.luminosityBlockAuxiliary().luminosityBlock();
  // accumulate HF data at every LS as it is closed.
  // note: lumi unit from DIPLumiSummary and Detail is microbarns
  edm::Handle<LumiSummary> lumiSummary_;
  lumiBlock.getByToken(lumiSummaryToken_, lumiSummary_);
  if (lumiSummary_->isValid()) {
    lumi = lumiSummary_->avgInsDelLumi();
    intDelLumi = lumiSummary_->intgDelLumi();
    std::cout << "Luminosity in this Lumi Section " << lumi << " --> " << intDelLumi << std::endl;
  } else {
    std::cout << "No valid data found!" << std::endl;
  }

  return lumi;
}

double GetLumi::getValue(edm::LuminosityBlock const& lumiBlock, edm::EventSetup const& eSetup) {
  return getRawValue(lumiBlock, eSetup) * lumiScale_;
}

double GetLumi::convert2PU(double instLumi,
                           double inelastic_xSec = GetLumi::INELASTIC_XSEC_8TeV)  // inelastic_xSec in mb
{
  // from https://cmswbm.web.cern.ch/cmswbm/images/pileup.png
  return instLumi * inelastic_xSec / FREQ_ORBIT;
}

double GetLumi::convert2PU(double instLumi, int sqrt_s = GetLumi::SQRT_S_8TeV) {
  double inelastic_xSec = 0.;

  switch (sqrt_s) {
    case GetLumi::SQRT_S_7TeV:
      inelastic_xSec = GetLumi::INELASTIC_XSEC_7TeV;
      break;
    case GetLumi::SQRT_S_8TeV:
      inelastic_xSec = GetLumi::INELASTIC_XSEC_8TeV;
      break;
    default:
      break;
  }

  return convert2PU(instLumi, inelastic_xSec);
}