|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
/*
* See header file for a description of this class.
*
* \author: Mia Tosi,40 3-B32,+41227671609
*/
#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() = default;
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();
edm::LogInfo("GetLumi") << "Luminosity in this Lumi Section " << lumi << " --> " << intDelLumi << std::endl;
} else {
edm::LogWarning("GetLumi") << "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);
}
|