|
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
|
/*
* See header file for a description of this class.
*
* \author A. Vilela Pereira
*/
#include "DTTTrigT0SegCorrection.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/MuonDetId/interface/DTSuperLayerId.h"
#include "CondFormats/DTObjects/interface/DTTtrig.h"
#include "CondFormats/DataRecord/interface/DTTtrigRcd.h"
#include "TFile.h"
#include "TH1F.h"
#include <string>
#include <sstream>
using namespace std;
using namespace edm;
namespace dtCalibration {
DTTTrigT0SegCorrection::DTTTrigT0SegCorrection(const ParameterSet& pset, edm::ConsumesCollector cc) {
string t0SegRootFile = pset.getParameter<string>("t0SegRootFile");
rootFile_ = new TFile(t0SegRootFile.c_str(), "READ");
ttrigToken_ =
cc.esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", pset.getUntrackedParameter<string>("dbLabel")));
}
DTTTrigT0SegCorrection::~DTTTrigT0SegCorrection() { delete rootFile_; }
void DTTTrigT0SegCorrection::setES(const EventSetup& setup) {
// Get tTrig record from DB
ESHandle<DTTtrig> tTrig = setup.getHandle(ttrigToken_);
tTrigMap_ = &*tTrig;
}
DTTTrigData DTTTrigT0SegCorrection::correction(const DTSuperLayerId& slId) {
float tTrigMean, tTrigSigma, kFactor;
int status = tTrigMap_->get(slId, tTrigMean, tTrigSigma, kFactor, DTTimeUnits::ns);
if (status != 0)
throw cms::Exception("[DTTTrigT0SegCorrection]") << "Could not find tTrig entry in DB for" << slId << endl;
const TH1F* t0SegHisto = getHisto(slId);
double corrMean = tTrigMean;
double corrSigma = tTrigSigma;
//FIXME: can we fit the t0seg histo? How do we remove the peak at 0?;
double corrKFact = (kFactor * tTrigSigma + t0SegHisto->GetMean()) / tTrigSigma;
return DTTTrigData(corrMean, corrSigma, corrKFact);
}
const TH1F* DTTTrigT0SegCorrection::getHisto(const DTSuperLayerId& slId) {
string histoName = getHistoName(slId);
TH1F* histo = static_cast<TH1F*>(rootFile_->Get(histoName.c_str()));
if (!histo)
throw cms::Exception("[DTTTrigT0SegCorrection]") << "t0-seg histogram not found:" << histoName << endl;
return histo;
}
string DTTTrigT0SegCorrection::getHistoName(const DTSuperLayerId& slId) {
DTChamberId chId = slId.chamberId();
// Compose the chamber name
std::string wheel = std::to_string(chId.wheel());
std::string station = std::to_string(chId.station());
std::string sector = std::to_string(chId.sector());
string chHistoName = "_W" + wheel + "_St" + station + "_Sec" + sector;
return (slId.superLayer() != 2) ? ("hRPhiSegT0" + chHistoName) : ("hRZSegT0" + chHistoName);
}
} // namespace dtCalibration
|