|
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
|
#ifndef PhotonEnergyCalibratorRun2_h
#define PhotonEnergyCalibratorRun2_h
#include <TRandom.h>
#include "EgammaAnalysis/ElectronTools/interface/EnergyScaleCorrection_class.h"
#include "EgammaAnalysis/ElectronTools/interface/SimplePhoton.h"
#include "FWCore/Utilities/interface/StreamID.h"
#include <vector>
class PhotonEnergyCalibratorRun2 {
public:
// dummy constructor for persistence
PhotonEnergyCalibratorRun2() {}
// further configuration will be added when we will learn how it will work
PhotonEnergyCalibratorRun2(bool isMC, bool synchronization, std::string correctionFile);
~PhotonEnergyCalibratorRun2();
/// Initialize with a random number generator (if not done, it will use the CMSSW service)
/// Caller code owns the TRandom.
void initPrivateRng(TRandom *rnd);
/// Correct this electron.
/// StreamID is needed when used with CMSSW Random Number Generator
void calibrate(SimplePhoton &photon, edm::StreamID const &id = edm::StreamID::invalidStreamID()) const;
/// Correct this electron.
/// StreamID is needed when used with CMSSW Random Number Generator
void calibrate(reco::Photon &photon,
unsigned int runNumber,
edm::StreamID const &id = edm::StreamID::invalidStreamID()) const;
protected:
// whatever data will be needed
bool isMC_, synchronization_;
TRandom *rng_;
std::vector<double> smearings_;
std::vector<double> scales_;
/// Return a number distributed as a unit gaussian, drawn from the private RNG if initPrivateRng was called,
/// or from the CMSSW RandomNumberGenerator service
/// If synchronization is set to true, it returns a fixed number (1.0)
double gauss(edm::StreamID const &id) const;
EnergyScaleCorrection_class _correctionRetriever;
};
#endif
|