|
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
|
// -*- C++ -*-
//
// Package: FromClusterSummaryMultiplicityProducer
// Class: FromClusterSummaryMultiplicityProducer
//
/**\class FromClusterSummaryMultiplicityProducer FromClusterSummaryMultiplicityProducer.cc DPGAnalysis/SiStripTools/plugins/FromClusterSummaryMultiplicityProducer.cc
Description: EDProducer of multiplicity maps
Implementation:
*/
//
// Original Author: Andrea Venturi
// Created: Fri Dec 04 2009
//
//
// system include files
#include <memory>
#include <string>
#include <map>
#include <vector>
// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/TrackerCommon/interface/ClusterSummary.h"
//
// class decleration
//
class FromClusterSummaryMultiplicityProducer : public edm::global::EDProducer<> {
public:
explicit FromClusterSummaryMultiplicityProducer(const edm::ParameterSet&);
private:
void produce(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
// ----------member data ---------------------------
edm::EDGetTokenT<ClusterSummary> m_collectionToken;
std::vector<ClusterSummary::CMSTracker> m_subdetenums;
std::vector<int> m_subdetsel;
ClusterSummary::VariablePlacement m_subdetvar;
};
//
// constants, enums and typedefs
//
//
// static data member definitions
//
//
// constructors and destructor
//
FromClusterSummaryMultiplicityProducer::FromClusterSummaryMultiplicityProducer(const edm::ParameterSet& iConfig)
: m_collectionToken(consumes<ClusterSummary>(iConfig.getParameter<edm::InputTag>("clusterSummaryCollection"))),
m_subdetenums(),
m_subdetsel(),
m_subdetvar(ClusterSummary::NCLUSTERS) {
produces<std::map<unsigned int, int> >();
//now do what ever other initialization is needed
std::vector<edm::ParameterSet> wantedsubds(iConfig.getParameter<std::vector<edm::ParameterSet> >("wantedSubDets"));
m_subdetenums.reserve(wantedsubds.size());
m_subdetsel.reserve(wantedsubds.size());
for (std::vector<edm::ParameterSet>::iterator ps = wantedsubds.begin(); ps != wantedsubds.end(); ++ps) {
m_subdetenums.push_back((ClusterSummary::CMSTracker)ps->getParameter<int>("subDetEnum"));
m_subdetsel.push_back(ps->getParameter<int>("subDetEnum"));
}
m_subdetvar = (ClusterSummary::VariablePlacement)iConfig.getParameter<int>("varEnum");
}
//
// member functions
//
// ------------ method called to produce the data ------------
void FromClusterSummaryMultiplicityProducer::produce(edm::StreamID,
edm::Event& iEvent,
const edm::EventSetup& iSetup) const {
LogDebug("Multiplicity") << " Ready to go";
using namespace edm;
std::unique_ptr<std::map<unsigned int, int> > mults(new std::map<unsigned int, int>);
Handle<ClusterSummary> clustsumm;
iEvent.getByToken(m_collectionToken, clustsumm);
switch (m_subdetvar) {
case ClusterSummary::NCLUSTERS:
for (unsigned int iS = 0; iS < m_subdetenums.size(); ++iS)
(*mults)[m_subdetsel[iS]] = int(clustsumm->getNClus(m_subdetenums[iS]));
break;
case ClusterSummary::CLUSTERSIZE:
for (unsigned int iS = 0; iS < m_subdetenums.size(); ++iS)
(*mults)[m_subdetsel[iS]] = int(clustsumm->getClusSize(m_subdetenums[iS]));
break;
case ClusterSummary::CLUSTERCHARGE:
for (unsigned int iS = 0; iS < m_subdetenums.size(); ++iS)
(*mults)[m_subdetsel[iS]] = int(clustsumm->getClusCharge(m_subdetenums[iS]));
break;
default:
for (unsigned int iS = 0; iS < m_subdetenums.size(); ++iS)
(*mults)[m_subdetsel[iS]] = -1;
}
for (unsigned int iS = 0; iS < m_subdetenums.size(); ++iS)
LogDebug("Multiplicity") << "GetModuleLocation result: " << m_subdetenums[iS] << " "
<< clustsumm->getModuleLocation(m_subdetenums[iS]);
for (std::map<unsigned int, int>::const_iterator it = mults->begin(); it != mults->end(); ++it) {
LogDebug("Multiplicity") << " Found " << it->second << " digis/clusters in " << it->first;
}
iEvent.put(std::move(mults));
}
DEFINE_FWK_MODULE(FromClusterSummaryMultiplicityProducer);
|