|
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
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
|
#ifndef DTCalibration_DTCalibrationMap_H
#define DTCalibration_DTCalibrationMap_H
/** \class DTCalibrationMap
* Allow saving and retrieving of calibration constants to/from txt file.
* This was originally provided for backward compatibility with the ORCA MuBarDigiParameters file.
* Can be used to save an arbitrary number of constants with the
* needed granularity and to retrieve them back using the wireId.
* Current field allocation: fields for each key are allocated to:
* --First block: legacy descriptors--
* [0] ttrig
* [1] sigma_ttrig (obsolete)
* [2] kfactor (obsolete)
* [3] vdrift
* [4] sigma_vdrift (obsolete, was formerly hacked to include reoslution)
* [5] t0
* [6] t0rms
* [7] noisy or dead flag
* [8-9] left for future usage
* --Second block (optional): free fields
* [10] Encoded information of free fields: (1000*version)+(100*type)+(number of fields); type is: ttrig=0, vdrift=1, uncertainties=3
* [11-end] free fields
*
* \author G. Cerminara - INFN Torino
*/
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include <string>
#include <map>
#include <vector>
namespace edm {
class ParameterSet;
}
class DTCalibrationMap {
public:
/// Constructor
DTCalibrationMap(const edm::ParameterSet& pset);
/// Destructor
virtual ~DTCalibrationMap();
// Operations
/// Return the t_trig (ns) for a particular wire
float tTrig(DTWireId wireId) const;
/// Return the sigma of the t_trig (ns) for a particular wire
float sigma_tTrig(DTWireId wireId) const;
/// Return the kfactor for a particular wire
float kFactor(DTWireId wireId) const;
/// Return the mean drift velocity for a particular wire (cm/ns)
float meanVDrift(DTWireId wireId) const;
/// Return the sigma of the mean drift velocity for a particular wire (cm/ns)
float sigma_meanVDrift(DTWireId wireId) const;
typedef std::vector<float> CalibConsts;
typedef DTWireId Key;
typedef std::map<Key, CalibConsts>::const_iterator const_iterator;
// Clean the map
void cleanTheConsts() { theMap.clear(); }
// Get a particular number (field) between all the calibration
// constants available for a particluar wire
float getField(DTWireId wireId, int field) const;
// Get from the map the calibration constants for a particular wire
const CalibConsts* getConsts(DTWireId wireId) const;
// Add to the map the calibration consts for a given key
void addCell(Key wireId, const CalibConsts& calibConst);
// Write the calibration consts to a file
void writeConsts(const std::string& outputFileName) const;
// Get a key to read calibration constants for a particular wire
// with the given granularity
Key getKey(DTWireId wireId) const;
const_iterator keyAndConsts_begin() const { return theMap.begin(); }
const_iterator keyAndConsts_end() const { return theMap.end(); }
protected:
private:
// Specify the granularity for the calibration constants
enum CalibGranularity { byChamber, bySL, byLayer, byWire };
CalibGranularity theGranularity;
// Read the calibration consts from a file
void readConsts(const std::string& inputFileName);
// Check the consistency of a given key with the selected granularity
bool checkGranularity(Key aKey) const;
// The number of fields (calibration numbers) to be read from file
unsigned int nFields;
// The name of the file containing the calibration constants
std::string calibConstFileName;
// Define the granularity to be used for t0
std::string calibConstGranularity;
// The map between the Key and the calibration constants
std::map<Key, CalibConsts> theMap;
};
#endif
|