|
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
|
/**_________________________________________________________________
class: AlcaPCCProducer.cc
authors:Sam Higginbotham (shigginb@cern.ch) and Chris Palmer (capalmer@cern.ch)
________________________________________________________________**/
// C++ standard
#include <string>
// CMS
#include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"
#include "DataFormats/Luminosity/interface/PixelClusterCounts.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/LuminosityBlock.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
//The class
class AlcaPCCProducer
: public edm::one::EDProducer<edm::EndLuminosityBlockProducer, edm::LuminosityBlockCache<reco::PixelClusterCounts>> {
public:
explicit AlcaPCCProducer(const edm::ParameterSet&);
private:
std::shared_ptr<reco::PixelClusterCounts> globalBeginLuminosityBlock(edm::LuminosityBlock const& lumiSeg,
edm::EventSetup const& iSetup) const final;
void endLuminosityBlockProduce(edm::LuminosityBlock& lumiSeg, const edm::EventSetup& iSetup) final;
void globalEndLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) final {}
void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) final;
edm::EDGetTokenT<edmNew::DetSetVector<SiPixelCluster>> pixelToken_;
edm::EDPutTokenT<reco::PixelClusterCounts> putToken_;
};
//--------------------------------------------------------------------------------------------------
AlcaPCCProducer::AlcaPCCProducer(const edm::ParameterSet& iConfig)
: pixelToken_(consumes(iConfig.getParameter<edm::InputTag>("pixelClusterLabel"))),
//specifies the trigger Rand or ZeroBias
putToken_(produces<reco::PixelClusterCounts, edm::Transition::EndLuminosityBlock>(
iConfig.getUntrackedParameter<std::string>("trigstring", "alcaPCC"))) {}
//--------------------------------------------------------------------------------------------------
void AlcaPCCProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
unsigned int bx = iEvent.bunchCrossing();
//std::cout<<"The Bunch Crossing"<<bx<<std::endl;
auto* pccOb = luminosityBlockCache(iEvent.getLuminosityBlock().index());
pccOb->eventCounter(bx);
//Looping over the clusters and adding the counts up
const edmNew::DetSetVector<SiPixelCluster>& clustColl = iEvent.get(pixelToken_);
// ----------------------------------------------------------------------
// -- Clusters without tracks
for (auto const& mod : clustColl) {
if (mod.empty()) {
continue;
}
DetId detId = mod.id();
//--The following will be used when we make a theshold for the clusters.
//--Keeping this for features that may be implemented later.
// -- clusters on this det
//edmNew::DetSet<SiPixelCluster>::const_iterator di;
//int nClusterCount=0;
//for (di = mod.begin(); di != mod.end(); ++di) {
// nClusterCount++;
//}
int nCluster = mod.size();
pccOb->increment(detId(), bx, nCluster);
}
}
//--------------------------------------------------------------------------------------------------
std::shared_ptr<reco::PixelClusterCounts> AlcaPCCProducer::globalBeginLuminosityBlock(
edm::LuminosityBlock const& lumiSeg, edm::EventSetup const& iSetup) const {
//New PCC object at the beginning of each lumi section
return std::make_shared<reco::PixelClusterCounts>();
}
//--------------------------------------------------------------------------------------------------
void AlcaPCCProducer::endLuminosityBlockProduce(edm::LuminosityBlock& lumiSeg, const edm::EventSetup& iSetup) {
//Saving the PCC object
lumiSeg.emplace(putToken_, std::move(*luminosityBlockCache(lumiSeg.index())));
}
DEFINE_FWK_MODULE(AlcaPCCProducer);
|