|
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
|
#ifndef DTNoiseCalibration_H
#define DTNoiseCalibration_H
/*
* \file DTNoiseCalibration.h
*
* \author G. Mila - INFN Torino
* A. Vilela Pereira - INFN Torino
*
*/
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "CondFormats/DataRecord/interface/DTTtrigRcd.h"
#include "DataFormats/DTDigi/interface/DTDigiCollection.h"
#include <string>
#include <vector>
#include <map>
#include <ctime>
class DTGeometry;
class DTChamberId;
class DTSuperLayerId;
class DTLayerId;
class DTWireId;
class DTTtrig;
class TFile;
class TH2F;
class TH1F;
class DTNoiseCalibration : public edm::one::EDAnalyzer<edm::one::WatchRuns> {
public:
/// Constructor
DTNoiseCalibration(const edm::ParameterSet& ps);
/// Destructor
~DTNoiseCalibration() override;
void beginJob() override;
void beginRun(const edm::Run& run, const edm::EventSetup& setup) override;
void analyze(const edm::Event& e, const edm::EventSetup& c) override;
void endRun(const edm::Run& run, const edm::EventSetup& setup) override {}
void endJob() override;
private:
std::string getChannelName(const DTWireId&) const;
// Get the name of the layer
std::string getLayerName(const DTLayerId&) const;
// Get the name of the superLayer
std::string getSuperLayerName(const DTSuperLayerId&) const;
// Get the name of the chamber
std::string getChamberName(const DTChamberId&) const;
const edm::EDGetTokenT<DTDigiCollection> digiToken_;
const bool useTimeWindow_;
const double triggerWidth_;
const int timeWindowOffset_;
const double maximumNoiseRate_;
const bool useAbsoluteRate_;
bool readDB_;
int defaultTtrig_;
std::vector<DTWireId> wireIdWithHisto_;
unsigned int lumiMax_;
int nevents_;
//int counter;
time_t runBeginTime_;
time_t runEndTime_;
// Get the DT Geometry
edm::ESHandle<DTGeometry> dtGeom_;
const edm::ESGetToken<DTGeometry, MuonGeometryRecord> dtGeomToken_;
// tTrig map
edm::ESHandle<DTTtrig> tTrigMap_;
const edm::ESGetToken<DTTtrig, DTTtrigRcd> ttrigToken_;
TFile* rootFile_;
// TDC digi distribution
TH1F* hTDCTriggerWidth_;
// Map of the occupancy histograms by layer
std::map<DTLayerId, TH1F*> theHistoOccupancyMap_;
// Map of occupancy by lumi by wire
std::map<DTWireId, TH1F*> theHistoOccupancyVsLumiMap_;
// Map of occupancy by lumi by chamber
std::map<DTChamberId, TH1F*> chamberOccupancyVsLumiMap_;
// Map of occupancy by time by chamber
std::map<DTChamberId, TH1F*> chamberOccupancyVsTimeMap_;
};
#endif
|