|
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
|
/*
* See header file for a description of this class.
*
* $Date: 2012/03/21 13:48:46 $
* $Revision: 1.2 $
* \author A. Vilela Pereira
*/
#include "DTT0FillChamberFromDB.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "CondFormats/DTObjects/interface/DTT0.h"
#include "CondFormats/DataRecord/interface/DTT0Rcd.h"
#include <string>
#include <sstream>
using namespace std;
using namespace edm;
namespace dtCalibration {
DTT0FillChamberFromDB::DTT0FillChamberFromDB(const ParameterSet& pset, edm::ConsumesCollector cc)
: chamberRef_(pset.getParameter<string>("chamberId")),
t0Token_(cc.esConsumes<edm::Transition::BeginRun>()),
t0RefToken_(
cc.esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", pset.getParameter<string>("dbLabelRef")))) {
//DTChamberId chosenChamberId;
if (!chamberRef_.empty() && chamberRef_ != "None") {
stringstream linestr;
int selWheel, selStation, selSector;
linestr << chamberRef_;
linestr >> selWheel >> selStation >> selSector;
chosenChamberId_ = DTChamberId(selWheel, selStation, selSector);
LogVerbatim("Calibration") << "[DTT0FillChamberFromDB] Chosen chamber: " << chosenChamberId_ << endl;
}
//FIXME: Check if chosen chamber is valid.
}
DTT0FillChamberFromDB::~DTT0FillChamberFromDB() {}
void DTT0FillChamberFromDB::setES(const EventSetup& setup) {
// Get t0 record from DB
ESHandle<DTT0> t0H = setup.getHandle(t0Token_);
t0Map_ = &*t0H;
LogVerbatim("Calibration") << "[DTT0FillChamberFromDB] T0 version: " << t0H->version();
// Get reference t0 DB
ESHandle<DTT0> t0RefH = setup.getHandle(t0RefToken_);
t0MapRef_ = &*t0RefH;
LogVerbatim("Calibration") << "[DTT0FillChamberFromDB] Reference T0 version: " << t0RefH->version();
}
DTT0Data DTT0FillChamberFromDB::correction(const DTWireId& wireId) {
// If wire belongs to chosen chamber, use t0 value from reference DB
// Otherwise use value from default DB
DTChamberId chamberId = wireId.layerId().superlayerId().chamberId();
if (!chamberRef_.empty() && chamberRef_ != "None" && chamberId == chosenChamberId_) {
// Access reference DB
float t0MeanRef, t0RMSRef;
int statusRef = t0MapRef_->get(wireId, t0MeanRef, t0RMSRef, DTTimeUnits::counts);
if (!statusRef) {
return DTT0Data(t0MeanRef, t0RMSRef);
} else {
//...
throw cms::Exception("[DTT0FillChamberFromDB]")
<< "Could not find t0 entry in reference DB for" << wireId << endl;
}
} else {
// Access default DB
float t0Mean, t0RMS;
int status = t0Map_->get(wireId, t0Mean, t0RMS, DTTimeUnits::counts);
if (!status) {
return DTT0Data(t0Mean, t0RMS);
} else {
//...
throw cms::Exception("[DTT0FillChamberFromDB]") << "Could not find t0 entry in DB for" << wireId << endl;
}
}
}
} // namespace dtCalibration
|