|
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
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
|
#ifndef DTVDriftCalibration_H
#define DTVDriftCalibration_H
/** \class DTVDriftCalibration
* No description available.
*
* \author M. Giunta
*/
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
#include "CalibMuon/DTCalibration/interface/vDriftHistos.h"
#include "CalibMuon/DTCalibration/interface/DTTMax.h"
#include "CalibMuon/DTCalibration/interface/DTSegmentSelector.h"
#include "DTCalibrationMap.h"
#include <string>
#include <vector>
namespace edm {
class ParameterSet;
class Event;
class EventSetup;
} // namespace edm
class TFile;
class DTMeanTimerFitter;
class DTGeometry;
class MuonGeometryRecord;
class DTVDriftCalibration : public edm::one::EDAnalyzer<> {
public:
/// Constructor
DTVDriftCalibration(const edm::ParameterSet& pset);
/// Destructor
~DTVDriftCalibration() override;
// Operations
void analyze(const edm::Event& event, const edm::EventSetup& eventSetup) override;
void endJob() override;
protected:
private:
std::unique_ptr<DTSegmentSelector> select_;
// The class containing TMax information
typedef DTTMax::TMax TMax;
dtcalibration::Histograms histograms_;
// class to create/manage histos for each partition (SL)
class cellInfo {
public:
cellInfo(const TString& name) { histos = new hTMaxCell(name); }
~cellInfo() { delete histos; }
void add(const std::vector<const TMax*>& tMaxes);
void update() { addedCells.clear(); }
hTMaxCell* getHists() { return histos; }
private:
cellInfo() {}
cellInfo(const cellInfo&) {}
std::vector<dttmaxenums::TMaxCells> addedCells;
hTMaxCell* histos;
};
TH1F* hChi2;
h2DSegm* h2DSegmRZ;
h2DSegm* h2DSegmRPhi;
h4DSegm* h4DSegmAllCh;
// Divide cellInfo by given granularity (to be implemented)
// DTVDriftCalibration::cellInfo* partition(const DTWireId& wireId);
// Specify the granularity for the TMax histograms
enum TMaxGranularity { byChamber, bySL, byPartition };
TMaxGranularity theGranularity;
// The label used to retrieve 4D segments from the event
edm::EDGetTokenT<DTRecSegment4DCollection> theRecHits4DToken;
const edm::ESGetToken<DTGeometry, MuonGeometryRecord> theDTGeomToken;
// Debug flag
bool debug;
// The label used to retrieve digis from the event
std::string digiLabel;
// The file which will contain the tMax histograms
TFile* theFile;
// The fitter
std::unique_ptr<DTMeanTimerFitter> theFitter;
// Perform the vDrift and t0 evaluation or just fill the
// tMaxHists (if you read the dataset in different jobs)
bool findVDriftAndT0;
// The name of the output text file
std::string theVDriftOutputFile;
// Map of wires and cellInfo with coarse granularity
std::map<DTWireId, cellInfo*> theWireIdAndCellMap;
// Switch for checking of noisy channels
//bool checkNoisyChannels;
// The module for t0 subtraction
std::unique_ptr<DTTTrigBaseSync> theSync; //FIXME: should be const
// parameter set for DTCalibrationMap constructor
edm::ParameterSet theCalibFilePar;
// Maximum value for the 4D Segment chi2
//double theMaxChi2;
// Maximum incident angle for Phi Seg
//double theMaxPhiAngle;
// Maximum incident angle for Theta Seg
//double theMaxZAngle;
// Choose the chamber you want to calibrate
std::string theCalibChamber;
// which format to be created
bool writeLegacyVDriftDB;
};
#endif
|