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