GetLumi Class Reference

#include <GetLumi.h>

Public Types
enum	SQRT_S { SQRT_S_7TeV, SQRT_S_8TeV, SQRT_S_13TeV }

Public Member Functions
double	convert2PU (double, double)
double	convert2PU (double, int)
	GetLumi (const edm::ParameterSet &)
	GetLumi (const edm::InputTag &, double)
	GetLumi (const edm::ParameterSet &, edm::ConsumesCollector &iC)
	GetLumi (const edm::InputTag &, double, edm::ConsumesCollector &iC)
double	getRawValue (const edm::Event &)
double	getRawValue (edm::LuminosityBlock const &, edm::EventSetup const &)
double	getValue (const edm::Event &)
double	getValue (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual	~GetLumi ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Public Attributes
edm::EDGetTokenT< LumiDetails >	lumiDetailsToken_
edm::InputTag	lumiInputTag_
double	lumiScale_
edm::EDGetTokenT< LumiSummary >	lumiSummaryToken_

Static Public Attributes
static double	CM2_TO_NANOBARN = 1.0/1.e-33
static double	FREQ_ORBIT = 11246.
static double	INELASTIC_XSEC_13TeV = 71.3
static double	INELASTIC_XSEC_7TeV = 68.0
static double	INELASTIC_XSEC_8TeV = 69.3
static const unsigned int	lastBunchCrossing = 3564
static const unsigned int	NUM_BX = 3564
static double	rXSEC_PIXEL_CLUSTER = 9.4e-24
static double	rXSEC_PIXEL_CLUSTER_UNC = 0.119e-24
static double	SECONDS_PER_LS = double(0x40000)/double(FREQ_ORBIT)
static double	XSEC_PIXEL_CLUSTER = 10.08e-24
static double	XSEC_PIXEL_CLUSTER_UNC = 0.17e-24

Detailed Description

Definition at line 21 of file GetLumi.h.

Member Enumeration Documentation

enum GetLumi::SQRT_S

Enumerator
SQRT_S_7TeV
SQRT_S_8TeV
SQRT_S_13TeV

Definition at line 25 of file GetLumi.h.

             {
    SQRT_S_7TeV,
    SQRT_S_8TeV,
    SQRT_S_13TeV
  };

Constructor & Destructor Documentation

GetLumi::GetLumi

(

const edm::ParameterSet &

iConfig

)

Definition at line 14 of file GetLumi.cc.

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

GetLumi::GetLumi	(	const edm::InputTag &	lumiInputTag,
		double	lumiScale
	)

Definition at line 20 of file GetLumi.cc.

  : lumiInputTag_ ( lumiInputTag )
  , lumiScale_    ( lumiScale    )
{
}

GetLumi::GetLumi	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &	iC
	)

Definition at line 26 of file GetLumi.cc.

References edm::ConsumesCollector::consumes(), edm::InLumi, lumiDetailsToken_, lumiInputTag_, and lumiSummaryToken_.

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

Definition at line 32 of file GetLumi.cc.

References edm::ConsumesCollector::consumes(), edm::InLumi, lumiDetailsToken_, lumiInputTag_, and lumiSummaryToken_.

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

GetLumi::~GetLumi ( )

virtual

Definition at line 38 of file GetLumi.cc.

{
}

Member Function Documentation

double GetLumi::convert2PU	(	double	instLumi,
		double	inelastic_xSec = GetLumi::INELASTIC_XSEC_8TeV
	)

Definition at line 106 of file GetLumi.cc.

References FREQ_ORBIT.

Referenced by convert2PU().

{

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

Definition at line 114 of file GetLumi.cc.

References convert2PU(), INELASTIC_XSEC_7TeV, INELASTIC_XSEC_8TeV, SQRT_S_7TeV, and 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);

}

static void GetLumi::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

double GetLumi::getRawValue

(

const edm::Event &

iEvent

)

Definition at line 43 of file GetLumi.cc.

References edm::EventBase::bunchCrossing(), edm::LuminosityBlock::getByToken(), edm::Event::getLuminosityBlock(), edm::HandleBase::isValid(), LumiDetails::kOCC1, fjr2json::lumi, and lumiSummaryToken_.

Referenced by getValue().

{

  // 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::getRawValue	(	edm::LuminosityBlock const &	lumiBlock,
		edm::EventSetup const &	eSetup
	)

Definition at line 71 of file GetLumi.cc.

References gather_cfg::cout, edm::LuminosityBlock::getByToken(), edm::HandleBase::isValid(), fjr2json::lumi, and lumiSummaryToken_.

{
  
  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

(

const edm::Event &

iEvent

)

Definition at line 64 of file GetLumi.cc.

References getRawValue(), and lumiScale_.

Referenced by VertexMonitor::analyze(), TrackingMonitor::analyze(), and LogMessageMonitor::analyze().

{
  //    bxlumi = lumi->lumiValue(LumiDetails::kOCC1,iEvent.bunchCrossing())*6.37;
  return getRawValue(iEvent)*lumiScale_;
}

double GetLumi::getValue	(	edm::LuminosityBlock const &	lumiBlock,
		edm::EventSetup const &	eSetup
	)

Definition at line 98 of file GetLumi.cc.

References getRawValue(), and lumiScale_.

{
  return getRawValue(lumiBlock,eSetup)*lumiScale_;
}

Member Data Documentation

double GetLumi::CM2_TO_NANOBARN = 1.0/1.e-33

static

Definition at line 47 of file GetLumi.h.

Referenced by LumiMonitor::analyze().

double GetLumi::FREQ_ORBIT = 11246.

static

Definition at line 32 of file GetLumi.h.

Referenced by convert2PU().

double GetLumi::INELASTIC_XSEC_13TeV = 71.3

static

Definition at line 37 of file GetLumi.h.

double GetLumi::INELASTIC_XSEC_7TeV = 68.0

static

Definition at line 35 of file GetLumi.h.

Referenced by convert2PU().

double GetLumi::INELASTIC_XSEC_8TeV = 69.3

static

Definition at line 36 of file GetLumi.h.

Referenced by convert2PU().

const unsigned int GetLumi::lastBunchCrossing = 3564

static

Definition at line 48 of file GetLumi.h.

edm::EDGetTokenT<LumiDetails> GetLumi::lumiDetailsToken_

Definition at line 71 of file GetLumi.h.

Referenced by GetLumi().

edm::InputTag GetLumi::lumiInputTag_

Definition at line 68 of file GetLumi.h.

Referenced by GetLumi().

double GetLumi::lumiScale_

Definition at line 69 of file GetLumi.h.

Referenced by getValue().

edm::EDGetTokenT<LumiSummary> GetLumi::lumiSummaryToken_

Definition at line 72 of file GetLumi.h.

Referenced by GetLumi(), and getRawValue().

const unsigned int GetLumi::NUM_BX = 3564

static

Definition at line 31 of file GetLumi.h.

double GetLumi::rXSEC_PIXEL_CLUSTER = 9.4e-24

static

Definition at line 45 of file GetLumi.h.

double GetLumi::rXSEC_PIXEL_CLUSTER_UNC = 0.119e-24

static

Definition at line 46 of file GetLumi.h.

double GetLumi::SECONDS_PER_LS = double(0x40000)/double(FREQ_ORBIT)

static

Definition at line 33 of file GetLumi.h.

double GetLumi::XSEC_PIXEL_CLUSTER = 10.08e-24

static

Definition at line 41 of file GetLumi.h.

double GetLumi::XSEC_PIXEL_CLUSTER_UNC = 0.17e-24

static

Definition at line 42 of file GetLumi.h.