|
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
|
/*
* See header file for a description of this class.
*
* \author S. Bolognesi
*/
#include "CalibMuon/DTCalibration/plugins/DTTPDeadWriter.h"
#include "CalibMuon/DTCalibration/interface/DTCalibDBUtils.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "CondFormats/DTObjects/interface/DTT0.h"
#include "CondFormats/DataRecord/interface/DTT0Rcd.h"
#include "CondFormats/DataRecord/interface/DTDeadFlagRcd.h"
#include "CondFormats/DTObjects/interface/DTDeadFlag.h"
/* C++ Headers */
#include <vector>
#include <set>
#include <iostream>
#include <fstream>
#include <string>
#include <sstream>
#include "TFile.h"
#include "TH1.h"
using namespace std;
using namespace edm;
// Constructor
DTTPDeadWriter::DTTPDeadWriter(const ParameterSet& pset) : dtGeomToken_(esConsumes<edm::Transition::BeginRun>()) {
// get selected debug option
debug = pset.getUntrackedParameter<bool>("debug", false);
t0Token_ = esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", pset.getParameter<string>("debug")));
// Create the object to be written to DB
tpDeadList = new DTDeadFlag();
if (debug)
cout << "[DTTPDeadWriter]Constructor called!" << endl;
}
// Destructor
DTTPDeadWriter::~DTTPDeadWriter() {
if (debug)
cout << "[DTTPDeadWriter]Destructor called!" << endl;
}
void DTTPDeadWriter::beginRun(const edm::Run&, const EventSetup& setup) {
// Get the t0 map
tZeroMap = &setup.getData(t0Token_);
// Get the muon Geometry
muonGeom = setup.getHandle(dtGeomToken_);
}
// Do the job
void DTTPDeadWriter::analyze(const Event& event, const EventSetup& eventSetup) {
set<DTLayerId> analyzedLayers;
//Loop on tzero map
for (DTT0::const_iterator tzero = tZeroMap->begin(); tzero != tZeroMap->end(); ++tzero) {
//Consider what layers have been already considered
// @@@ NEW DTT0 FORMAT
// DTLayerId layerId = (DTWireId((*tzero).first.wheelId,
// (*tzero).first.stationId,
// (*tzero).first.sectorId,
// (*tzero).first.slId,
// (*tzero).first.layerId,
// (*tzero).first.cellId)).layerId();
int channelId = tzero->channelId;
if (channelId == 0)
continue;
DTLayerId layerId = (DTWireId(channelId)).layerId();
// @@@ NEW DTT0 END
if (analyzedLayers.find(layerId) == analyzedLayers.end()) {
analyzedLayers.insert(layerId);
//Take the layer topology
const DTTopology& dtTopo = muonGeom->layer(layerId)->specificTopology();
const int firstWire = dtTopo.firstChannel();
//const int lastWire = dtTopo.lastChannel();
const int nWires = muonGeom->layer(layerId)->specificTopology().channels();
//Loop on wires
for (int wire = firstWire; wire <= nWires; wire++) {
DTWireId wireId(layerId, wire);
float t0 = 0;
float t0rms = 0;
tZeroMap->get(wireId, t0, t0rms, DTTimeUnits::ns);
//If no t0 stored then is a tp dead channel
if (!t0) {
tpDeadList->setCellDead_TP(wireId, true);
cout << "Wire id " << wireId << " is TP dead" << endl;
}
}
}
}
}
// Write objects to DB
void DTTPDeadWriter::endJob() {
if (debug)
cout << "[DTTPDeadWriter]Writing ttrig object to DB!" << endl;
// FIXME: to be read from cfg?
string deadRecord = "DTDeadFlagRcd";
// Write the object to DB
DTCalibDBUtils::writeToDB(deadRecord, tpDeadList);
}
|