|
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
|
#ifndef DTChamberEfficiencyTask_H
#define DTChamberEfficiencyTask_H
/** \class DTChamberEfficiencyTask
* DQM Analysis of 4D DT segments, it produces plots about: <br>
* - single chamber efficiency
* All histos are produced per Chamber
*
* Class based on the code written by S. Lacaprara :
* RecoLocalMuon / DTSegment / test / DTEffAnalyzer.h
*
* \author G. Mila - INFN Torino
*/
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "DataFormats/MuonDetId/interface/DTLayerId.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMOneEDAnalyzer.h"
#include <string>
#include <map>
#include <vector>
class DTChamberEfficiencyTask : public DQMOneEDAnalyzer<edm::one::WatchLuminosityBlocks> {
public:
/// Constructor
DTChamberEfficiencyTask(const edm::ParameterSet& pset);
/// Destructor
~DTChamberEfficiencyTask() override;
/// BeginRun
void dqmBeginRun(const edm::Run& run, const edm::EventSetup& setup) override;
/// To reset the MEs
void beginLuminosityBlock(edm::LuminosityBlock const& lumiSeg, edm::EventSetup const& context) override;
void endLuminosityBlock(edm::LuminosityBlock const& lumiSeg, edm::EventSetup const& context) final {}
// Operations
void analyze(const edm::Event& event, const edm::EventSetup& setup) override;
protected:
// Book the histograms
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
private:
const DTRecSegment4D& getBestSegment(const DTRecSegment4DCollection::range& segs) const;
const DTRecSegment4D* getBestSegment(const DTRecSegment4D* s1, const DTRecSegment4D* s2) const;
bool isGoodSegment(const DTRecSegment4D& seg) const;
LocalPoint interpolate(const DTRecSegment4D& seg1, const DTRecSegment4D& seg3, const DTChamberId& MB2) const;
void bookHistos(DQMStore::IBooker& ibooker, DTChamberId chId);
// Switch for verbosity
bool debug;
// The running mode
bool onlineMonitor;
// The analysis mode
bool detailedAnalysis;
// Lable of 4D segments in the event
edm::EDGetTokenT<DTRecSegment4DCollection> recHits4DToken_;
edm::ParameterSet parameters;
std::map<DTChamberId, std::vector<MonitorElement*> > histosPerCh;
unsigned int theMinHitsSegment;
double theMinChi2NormSegment;
double theMinCloseDist;
//Load geometry
edm::ESGetToken<DTGeometry, MuonGeometryRecord> muonGeomToken_;
const DTGeometry* dtGeom;
edm::Handle<DTRecSegment4DCollection> segs;
};
#endif
/* Local Variables: */
/* show-trailing-whitespace: t */
/* truncate-lines: t */
/* End: */
|