|
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
|
#ifndef CalibMuon_DTT0Calibration_H
#define CalibMuon_DTT0Calibration_H
/** \class DTT0Calibration
* Analyzer class computes the mean and RMS of t0 from pulses.
* Those values are written in the DB with cell granularity. The
* mean value for each channel is normalized to a reference time common to all the sector.
* The t0 of wires in odd layers are corrected for the relative difference between
* odd and even layers
*
* \author S. Bolognesi - INFN Torino
*/
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/DTDigi/interface/DTDigiCollection.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include <string>
#include <vector>
#include <map>
#include "TFile.h"
#include "TH1I.h"
#include "TH1D.h"
#include "TSpectrum.h"
//class TFile;
//class TH1I;
//class TH1D;
//class TSpectrum;
class DTT0;
class DTT0Calibration : public edm::one::EDAnalyzer<> {
public:
/// Constructor
DTT0Calibration(const edm::ParameterSet& pset);
/// Destructor
~DTT0Calibration() override;
// Operations
/// Fill the maps with t0 (by channel)
void analyze(const edm::Event& event, const edm::EventSetup& eventSetup) override;
/// Compute the mean and the RMS of the t0 from the maps and write them to the DB with channel granularity
void endJob() override;
protected:
private:
// Generate the histo name
std::string getHistoName(const DTWireId& wId) const;
std::string getHistoName(const DTLayerId& lId) const;
// Debug flag
bool debug;
// The label used to retrieve digis from the event
edm::EDGetTokenT<DTDigiCollection> digiToken;
// The root file which contain the histos per layer
TFile theFile;
//The event counter
unsigned int nevents;
//Number of events to be used for the t0 per layer histos
unsigned int eventsForLayerT0;
//Number of events to be used for the t0 reference per wire
unsigned int eventsForWireT0;
//Acceptance of t0 w.r.t. reference peak
double tpPeakWidth;
//Acceptance of t0 w.r.t. reference peak
double tpPeakWidthPerLayer;
//Digi's will be rejected if too far from TP peak
unsigned int rejectDigiFromPeak;
//The wheels,sector to be calibrated (default All)
std::string theCalibWheel;
int selWheel;
std::string theCalibSector;
int selSector;
// Map of the histos and graph by layer
std::map<DTLayerId, TH1I> theHistoLayerMap;
// Histogram containing position of all peaks
TH1D hLayerPeaks;
TSpectrum spectrum;
//Layer with histos for each wire
std::vector<DTWireId> wireIdWithHistos;
std::vector<DTLayerId> layerIdWithWireHistos;
//Maps with t0, sigma, number of digi per wire
std::map<DTWireId, double> theAbsoluteT0PerWire;
std::map<DTWireId, double> theRelativeT0PerWire;
std::map<DTWireId, double> theSigmaT0PerWire;
std::map<DTWireId, int> nDigiPerWire;
std::map<DTWireId, int> nDigiPerWire_ref;
std::map<DTWireId, double> mK;
std::map<DTWireId, double> mK_ref;
std::map<DTWireId, double> qK;
//Map with histo per wire for the chosen layer
std::map<DTWireId, TH1I> theHistoWireMap;
//Map with mean and RMS of t0 per layer
std::map<std::string, double> theT0LayerMap;
std::map<std::string, double> theSigmaT0LayerMap;
std::map<DTLayerId, double> theTPPeakMap;
//Ref. t0 per chamber
std::map<DTChamberId, double> theSumT0ByChamber;
std::map<DTChamberId, int> theCountT0ByChamber;
std::map<DTChamberId, double> theSigmaT0ByChamber;
std::map<DTChamberId, double> theMeanT0ByChamber;
std::map<DTChamberId, double> theRefT0ByChamber;
//DTGeometry used to loop on the SL in the endJob
edm::ESHandle<DTGeometry> dtGeom;
const edm::ESGetToken<DTGeometry, MuonGeometryRecord> dtGeomToken_;
};
#endif
|