|
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
|
// -*- C++ -*-
//
// Package: SiStripChannelChargeFilter
// Class : TECClusterFilter
//
//
// Original Author: sfricke
#include "EventFilter/SiStripChannelChargeFilter/interface/TECClusterFilter.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
using namespace std;
namespace cms {
TECClusterFilter::TECClusterFilter(const edm::ParameterSet& ps) {
//
ModulesToBeExcluded.clear();
ModulesToBeExcluded = ps.getParameter<std::vector<unsigned> >("ModulesToBeExcluded");
edm::LogInfo("TECClusterFilter") << "Clusters from " << ModulesToBeExcluded.size()
<< " modules will be ignored in the filter:";
for (std::vector<uint32_t>::const_iterator imod = ModulesToBeExcluded.begin(); imod != ModulesToBeExcluded.end();
imod++) {
edm::LogInfo("TECClusterFilter") << *imod;
}
//
ChargeThresholdTEC = ps.getParameter<int>("ChargeThresholdTEC");
edm::LogInfo("TECClusterFilter") << "ChargeThresholdTEC" << ChargeThresholdTEC;
minNrOfTECClusters = ps.getParameter<int>("MinNrOfTECClusters");
edm::LogInfo("TECClusterFilter") << "MinNrOfTECClusters" << minNrOfTECClusters;
clusterProducer = ps.getParameter<string>("ClusterProducer");
edm::LogInfo("TECClusterFilter") << "ClusterProducer" << clusterProducer;
// also put decision in the event
produces<int>();
}
bool TECClusterFilter::filter(edm::Event& e, edm::EventSetup const& c) {
edm::Handle<edm::DetSetVector<SiStripCluster> > h; //get SiStripCluster
e.getByLabel(clusterProducer, h);
bool decision = false; // default value, only accept if set true in this loop
unsigned int nr_clusters_above_threshold = 0;
for (edm::DetSetVector<SiStripCluster>::const_iterator it = h->begin(); it != h->end(); it++) {
DetId thedetId = DetId(it->detId());
bool exclude_this_detid = false;
for (vector<SiStripCluster>::const_iterator vit = (it->data).begin(); vit != (it->data).end(); vit++) {
for (std::vector<uint32_t>::const_iterator imod = ModulesToBeExcluded.begin();
imod != ModulesToBeExcluded.end();
imod++) {
if (*imod == thedetId.rawId())
exclude_this_detid = true;
} // found in exclusion list
if ((!exclude_this_detid) &&
(thedetId.subdetId() == StripSubdetector::TEC)) // if not excluded and if TEC module
{ // calculate sum of amplitudes
unsigned int amplclus = 0;
// int amplclus=0;
for (auto ia = vit->amplitudes().begin(); ia != vit->amplitudes().end(); ia++) {
if ((*ia) > 0)
amplclus += (*ia);
} // why should this be negative?
if (amplclus > ChargeThresholdTEC)
nr_clusters_above_threshold++;
}
}
}
if (nr_clusters_above_threshold >= minNrOfTECClusters)
decision = true;
e.put(std::make_unique<int>(decision));
return decision;
}
} // namespace cms
|