|
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
|
/*
* See header file for a description of this class.
*
* \author A. Vilela Pereira
*/
#include "DTVDriftMeanTimer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "CalibMuon/DTCalibration/interface/DTMeanTimerFitter.h"
#include "CalibMuon/DTCalibration/interface/DTCalibDBUtils.h"
#include <string>
#include <vector>
#include "TFile.h"
#include "TString.h"
using namespace std;
using namespace edm;
namespace dtCalibration {
DTVDriftMeanTimer::DTVDriftMeanTimer(const ParameterSet& pset, edm::ConsumesCollector cc) {
string rootFileName = pset.getParameter<string>("rootFileName");
rootFile_ = new TFile(rootFileName.c_str(), "READ");
fitter_ = new DTMeanTimerFitter(rootFile_);
bool debug = pset.getUntrackedParameter<bool>("debug", false);
if (debug)
fitter_->setVerbosity(1);
}
DTVDriftMeanTimer::~DTVDriftMeanTimer() {
rootFile_->Close();
delete fitter_;
}
void DTVDriftMeanTimer::setES(const edm::EventSetup& setup) {}
DTVDriftData DTVDriftMeanTimer::compute(DTSuperLayerId const& slId) {
// Evaluate v_drift and sigma from the TMax histograms
DTWireId wireId(slId, 0, 0);
TString N = (((((TString) "TMax" + (long)wireId.wheel()) + (long)wireId.station()) + (long)wireId.sector()) +
(long)wireId.superLayer());
vector<float> vDriftAndReso = fitter_->evaluateVDriftAndReso(N);
// Don't write the constants for the SL if the vdrift was not computed
if (vDriftAndReso.front() == -1)
throw cms::Exception("DTCalibration") << "Could not compute valid vDrift value for SL " << slId << endl;
return DTVDriftData(vDriftAndReso[0], vDriftAndReso[1]);
}
} // namespace dtCalibration
|