|
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
128
129
130
131
132
133
|
#include "Alignment/LaserAlignment/interface/LASGlobalData.h"
#include "Alignment/LaserAlignment/interface/LASGlobalLoop.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
// Forward declarations
class TFile;
class TTree;
class RawDataConverter : public edm::one::EDAnalyzer<edm::one::WatchRuns> {
public:
explicit RawDataConverter(const edm::ParameterSet&);
~RawDataConverter() override = default;
private:
enum DigiType { ZeroSuppressed, VirginRaw, ProcessedRaw, Unknown };
void beginJob() override;
void beginRun(edm::Run const&, edm::EventSetup const&) override;
void analyze(const edm::Event&, const edm::EventSetup&) override;
void endRun(edm::Run const&, edm::EventSetup const&) override {}
void endJob() override;
void fillDetectorId(void);
void ClearData(void);
DigiType GetValidLabels(
const edm::Event&
iEvent); // Check what kind of file is being processed and get valid module and instance labels returns the type of Digis that was found
template <class T>
void GetDigis(const edm::Event&); // Copy the Digis into the local container (theData)
std::vector<std::string> theDigiModuleLabels;
std::vector<std::string> theProductInstanceLabels;
std::string CurrentModuleLabel;
std::string CurrentInstanceLabel;
TFile* theOutputFile;
TTree* theOutputTree;
LASGlobalData<std::vector<float> > theData;
int latency;
int eventnumber;
int runnumber;
int lumiBlock;
LASGlobalData<int> detectorId;
};
// Copy the Digis into the local container (theData)
// Currently this has only been implemented and tested for SiStripDigis
// SiStripRawDigis and SiStripProcessedRawDigis will need some changes to work (this is the final goal)
template <class Digitype>
void RawDataConverter::GetDigis(const edm::Event& iEvent) {
LogDebug("RawDataConverter") << "Fill ZeroSuppressed Digis into the Tree";
// Get the DetSetVector for the SiStripDigis
// This is a vector with all the modules, each module containing zero or more strips with signal (Digis)
edm::Handle<edm::DetSetVector<Digitype> > detSetVector; // Handle for holding the DetSetVector
iEvent.getByLabel(CurrentModuleLabel, CurrentInstanceLabel, detSetVector);
if (!detSetVector.isValid())
throw std::runtime_error("Could not find the Digis");
// set everything in the local container to zero
ClearData();
// Fill the Digis into the Raw Data Container
LASGlobalLoop loop; // loop helper
int det, ring, beam, disk, pos; // and its variables
// loop over TEC+- (internal) modules
det = 0;
ring = 0;
beam = 0;
disk = 0;
do {
// Find the module in the DetSetVector and get a pointer (iterator) to it
typename edm::DetSetVector<Digitype>::const_iterator theModule =
detSetVector->find(detectorId.GetTECEntry(det, ring, beam, disk));
if (theModule != detSetVector->end()) {
// loop over all the Digis in this Module
typename edm::DetSet<Digitype>::const_iterator theDigi;
for (theDigi = theModule->data.begin(); theDigi != theModule->data.end(); ++theDigi) {
// fill the number of adc counts into the local container
if (theDigi->channel() < 512)
theData.GetTECEntry(det, ring, beam, disk).at(theDigi->channel()) = theDigi->adc();
}
}
} while (loop.TECLoop(det, ring, beam, disk));
// loop TIB/TOB
det = 2;
beam = 0;
pos = 0; // <- set det = 2 (TIB)
do {
// Find the module in the DetSetVector and get a pointer (iterator) to it
typename edm::DetSetVector<Digitype>::const_iterator theModule =
detSetVector->find(detectorId.GetTIBTOBEntry(det, beam, pos));
if (theModule != detSetVector->end()) {
// loop over all the Digis in this Module
typename edm::DetSet<Digitype>::const_iterator theDigi;
for (theDigi = theModule->data.begin(); theDigi != theModule->data.end(); ++theDigi) {
// fill the number of adc counts into the local container
if (theDigi->channel() < 512)
theData.GetTIBTOBEntry(det, beam, pos).at(theDigi->channel()) = theDigi->adc();
}
}
} while (loop.TIBTOBLoop(det, beam, pos));
// loop TEC (AT)
det = 0;
beam = 0;
disk = 0;
do {
// Find the module in the DetSetVector and get a pointer (iterator) to it
typename edm::DetSetVector<Digitype>::const_iterator theModule =
detSetVector->find(detectorId.GetTEC2TECEntry(det, beam, disk));
if (theModule != detSetVector->end()) {
// loop over all the Digis in this Module
typename edm::DetSet<Digitype>::const_iterator theDigi;
for (theDigi = theModule->data.begin(); theDigi != theModule->data.end(); ++theDigi) {
// fill the number of adc counts into the local container
if (theDigi->channel() < 512)
theData.GetTEC2TECEntry(det, beam, disk).at(theDigi->channel()) = theDigi->adc();
}
}
} while (loop.TEC2TECLoop(det, beam, disk));
}
|