|
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
|
#ifndef DTCalibValidation_H
#define DTCalibValidation_H
/** \class DTCalibValidation
* Analysis on DT residuals to validate the kFactor
*
*
* \author G. Mila - INFN Torino
*/
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "DataFormats/MuonDetId/interface/DTSuperLayerId.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "DataFormats/DTRecHit/interface/DTRecHitCollection.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment2DCollection.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include <string>
#include <map>
#include <vector>
// To remove into CMSSW versions before 20X
// To add into CMSSW versions before 20X
//class DaqMonitorBEInterface;
class DTGeometry;
class DTChamber;
// FR class DTCalibValidation: public edm::EDAnalyzer{
class DTCalibValidation : public DQMEDAnalyzer {
public:
/// Constructor
DTCalibValidation(const edm::ParameterSet& pset);
/// Destructor
~DTCalibValidation() override;
/// BeginRun
void dqmBeginRun(const edm::Run& r, const edm::EventSetup& c) override;
// Operations
void analyze(const edm::Event& event, const edm::EventSetup& setup) override;
protected:
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
private:
// Switch for verbosity
//bool debug;
edm::ParameterSet parameters;
int wrongSegment;
int rightSegment;
int nevent;
// the analysis type
bool detailedAnalysis;
// the geometry
edm::ESGetToken<DTGeometry, MuonGeometryRecord> muonGeomToken_;
const DTGeometry* dtGeom;
// Lable of 1D rechits in the event
edm::EDGetTokenT<DTRecHitCollection> recHits1DToken_;
// Lable of 2D segments in the event
edm::EDGetTokenT<DTRecSegment2DCollection> segment2DToken_;
// Lable of 4D segments in the event
edm::EDGetTokenT<DTRecSegment4DCollection> segment4DToken_;
// Return a map between DTRecHit1DPair and wireId
std::map<DTWireId, std::vector<DTRecHit1DPair> > map1DRecHitsPerWire(const DTRecHitCollection* dt1DRecHitPairs);
// Return a map between DTRecHit1D and wireId
std::map<DTWireId, std::vector<DTRecHit1D> > map1DRecHitsPerWire(const DTRecSegment2DCollection* segment2Ds);
// Return a map between DTRecHit1D and wireId
std::map<DTWireId, std::vector<DTRecHit1D> > map1DRecHitsPerWire(const DTRecSegment4DCollection* segment4Ds);
template <typename type>
const type* findBestRecHit(const DTLayer* layer,
DTWireId wireId,
const std::vector<type>& recHits,
const float simHitDist);
// Compute the distance from wire (cm) of a hits in a DTRecHit1DPair
float recHitDistFromWire(const DTRecHit1DPair& hitPair, const DTLayer* layer);
// Compute the distance from wire (cm) of a hits in a DTRecHit1D
float recHitDistFromWire(const DTRecHit1D& recHit, const DTLayer* layer);
// Compute the position with respect to the wire (cm) of a hits in a DTRecHit1DPair
float recHitPosition(
const DTRecHit1DPair& hitPair, const DTLayer* layer, const DTChamber* chamber, float segmPos, int sl);
// Compute the position with respect to the wire (cm) of a hits in a DTRecHit1D
float recHitPosition(const DTRecHit1D& recHit, const DTLayer* layer, const DTChamber* chamber, float segmPos, int sl);
// Does the real job
template <typename type>
void compute(const DTGeometry* dtGeom,
const DTRecSegment4D& segment,
const std::map<DTWireId, std::vector<type> >& recHitsPerWire,
int step);
// Book a set of histograms for a give chamber
void bookHistos(DTSuperLayerId slId, int step);
// Fill a set of histograms for a give chamber
void fillHistos(DTSuperLayerId slId,
float distance,
float residualOnDistance,
float position,
float residualOnPosition,
int step);
std::map<std::pair<DTSuperLayerId, int>, std::vector<MonitorElement*> > histosPerSL;
};
#endif
/* Local Variables: */
/* show-trailing-whitespace: t */
/* truncate-lines: t */
/* End: */
|