GetLumi.h

Go to the documentation of this file.

#ifndef GetLumi_H
#define GetLumi_H

// system include files
#include <memory>

// user include files
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

class LumiDetails;
class LumiSummary;

class GetLumi {
public:
  enum SQRT_S { SQRT_S_7TeV, SQRT_S_8TeV, SQRT_S_13TeV };

  static const unsigned int NUM_BX = 3564;
  static constexpr double FREQ_ORBIT = 11246.;  // Hz
  static constexpr double SECONDS_PER_LS = double(0x40000) / double(FREQ_ORBIT);

  static constexpr double INELASTIC_XSEC_7TeV = 68.0;  // mb
  static constexpr double INELASTIC_XSEC_8TeV = 69.3;  // mb
  // mb from http://inspirehep.net/record/1447965/files/FSQ-15-005-pas.pdf
  static constexpr double INELASTIC_XSEC_13TeV = 71.3;

  // from http://cmslxr.fnal.gov/source/DQM/PixelLumi/plugins/PixelLumiDQM.h
  // Using all pixel clusters:
  static constexpr double XSEC_PIXEL_CLUSTER = 10.08e-24;  //in cm^2
  static constexpr double XSEC_PIXEL_CLUSTER_UNC = 0.17e-24;

  // Excluding the inner barrel layer.
  static constexpr double rXSEC_PIXEL_CLUSTER = 9.4e-24;  //in cm^2
  static constexpr double rXSEC_PIXEL_CLUSTER_UNC = 0.119e-24;
  static constexpr double CM2_TO_NANOBARN = 1.0 / 1.e-33;
  static const unsigned int lastBunchCrossing = 3564;

  GetLumi(const edm::ParameterSet&);
  GetLumi(const edm::InputTag&, double);
  GetLumi(const edm::ParameterSet&, edm::ConsumesCollector& iC);
  GetLumi(const edm::InputTag&, double, edm::ConsumesCollector& iC);
  virtual ~GetLumi();
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

  double getRawValue(const edm::Event&);
  double getValue(const edm::Event&);

  double getRawValue(edm::LuminosityBlock const&, edm::EventSetup const&);
  double getValue(edm::LuminosityBlock const&, edm::EventSetup const&);

  double convert2PU(double, double);
  double convert2PU(double, int);

  // ----------member data ---------------------------
  edm::InputTag lumiInputTag_;
  double lumiScale_;

  edm::EDGetTokenT<LumiDetails> lumiDetailsToken_;
  edm::EDGetTokenT<LumiSummary> lumiSummaryToken_;
};
#endif