CMS 3D CMS Logo
List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
ShallowClustersProducer Class Reference

#include <ShallowClustersProducer.h>

Inheritance diagram for ShallowClustersProducer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Classes

struct  moduleVars
 
struct  NearDigis
 

Public Member Functions

 ShallowClustersProducer (const edm::ParameterSet &)
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
SerialTaskQueue * globalLuminosityBlocksQueue ()
 
SerialTaskQueue * globalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDProducer () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void produce (edm::Event &, const edm::EventSetup &) override
 

Private Attributes

std::string Prefix
 
edm::InputTag theClustersLabel
 
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > > theClustersToken_
 
edm::EDGetTokenT< edm::DetSetVector< SiStripProcessedRawDigi > > theDigisToken_
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 17 of file ShallowClustersProducer.h.

Constructor & Destructor Documentation

ShallowClustersProducer::ShallowClustersProducer ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 12 of file ShallowClustersProducer.cc.

References edm::ParameterSet::getParameter(), Prefix, produce(), theClustersToken_, and theDigisToken_.

13  : Prefix(iConfig.getParameter<std::string>("Prefix") )
14 {
15  produces <std::vector<unsigned> > ( Prefix + "number" );
16  produces <std::vector<unsigned> > ( Prefix + "width" );
17  produces <std::vector<float> > ( Prefix + "variance" );
18  produces <std::vector<float> > ( Prefix + "barystrip" );
19  produces <std::vector<float> > ( Prefix + "middlestrip" );
20  produces <std::vector<unsigned> > ( Prefix + "charge" );
21  produces <std::vector<float> > ( Prefix + "noise" );
22  produces <std::vector<float> > ( Prefix + "ston" );
23  produces <std::vector<unsigned> > ( Prefix + "seedstrip" );
24  produces <std::vector<unsigned> > ( Prefix + "seedindex" );
25  produces <std::vector<unsigned> > ( Prefix + "seedcharge" );
26  produces <std::vector<float> > ( Prefix + "seednoise" );
27  produces <std::vector<float> > ( Prefix + "seedgain" );
28  produces <std::vector<unsigned> > ( Prefix + "qualityisbad" );
29 
30  produces <std::vector<float> > ( Prefix + "rawchargeC" );
31  produces <std::vector<float> > ( Prefix + "rawchargeL" );
32  produces <std::vector<float> > ( Prefix + "rawchargeR" );
33  produces <std::vector<float> > ( Prefix + "rawchargeLL" );
34  produces <std::vector<float> > ( Prefix + "rawchargeRR" );
35  produces <std::vector<float> > ( Prefix + "eta" );
36  produces <std::vector<float> > ( Prefix + "foldedeta" );
37  produces <std::vector<float> > ( Prefix + "etaX" );
38  produces <std::vector<float> > ( Prefix + "etaasymm" );
39  produces <std::vector<float> > ( Prefix + "outsideasymm");
40  produces <std::vector<float> > ( Prefix + "neweta");
41  produces <std::vector<float> > ( Prefix + "newetaerr");
42 
43  produces <std::vector<unsigned> > ( Prefix + "detid" );
44  produces <std::vector<int> > ( Prefix + "subdetid" );
45  produces <std::vector<int> > ( Prefix + "module" );
46  produces <std::vector<int> > ( Prefix + "side" );
47  produces <std::vector<int> > ( Prefix + "layerwheel" );
48  produces <std::vector<int> > ( Prefix + "stringringrod" );
49  produces <std::vector<int> > ( Prefix + "petal" );
50  produces <std::vector<int> > ( Prefix + "stereo" );
51 
52  theClustersToken_ = consumes<edmNew::DetSetVector<SiStripCluster> > (iConfig.getParameter<edm::InputTag>("Clusters"));
53  theDigisToken_ = consumes<edm::DetSetVector<SiStripProcessedRawDigi> > (edm::InputTag("siStripProcessedRawDigis", ""));
54 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
ShallowClustersProducer::theClustersToken_
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > > theClustersToken_
Definition: ShallowClustersProducer.h:45
ShallowClustersProducer::theDigisToken_
edm::EDGetTokenT< edm::DetSetVector< SiStripProcessedRawDigi > > theDigisToken_
Definition: ShallowClustersProducer.h:46

Member Function Documentation

void ShallowClustersProducer::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivate

Definition at line 57 of file ShallowClustersProducer.cc.

References edmNew::DetSetVector< T >::begin(), SiStripClusterInfo::charge(), ALCARECOTkAlJpsiMuMu_cff::charge, fastPrimaryVertexProducer_cfi::clusters, edmNew::DetSetVector< T >::end(), ShallowClustersProducer::NearDigis::eta(), PVValHelper::eta, ShallowClustersProducer::NearDigis::etaasymm(), ShallowClustersProducer::NearDigis::etaX(), ShallowClustersProducer::NearDigis::first, SiStripClusterInfo::firstStrip(), edm::EventSetup::get(), edm::Event::getByToken(), edmNew::DetSetVector< T >::id(), info(), SiStripClusterInfo::IsAnythingBad(), edm::HandleBase::isValid(), ShallowClustersProducer::NearDigis::last, ShallowClustersProducer::moduleVars::layerwheel, ShallowClustersProducer::NearDigis::left, ShallowClustersProducer::NearDigis::Lleft, ShallowClustersProducer::NearDigis::max, SiStripClusterInfo::maxCharge(), SiStripClusterInfo::maxIndex(), SiStripClusterInfo::maxStrip(), ShallowClustersProducer::moduleVars::module, eostools::move(), ShallowClustersProducer::NearDigis::NearDigis(), SiStripClusterInfo::noiseRescaledByGain(), ShallowClustersProducer::NearDigis::outsideasymm(), ShallowClustersProducer::moduleVars::petal, funct::pow(), Prefix, edm::ESHandle< T >::product(), edm::Event::put(), ShallowClustersProducer::NearDigis::right, ShallowClustersProducer::NearDigis::Rright, ShallowClustersProducer::moduleVars::side, SiStripClusterInfo::signalOverNoise(), mathSSE::sqrt(), ShallowClustersProducer::moduleVars::stereo, ShallowClustersProducer::moduleVars::stringringrod, SiStripClusterInfo::stripGains(), SiStripClusterInfo::stripNoisesRescaledByGain(), ShallowClustersProducer::moduleVars::subdetid, theClustersToken_, theDigisToken_, SiStripClusterInfo::variance(), ApeEstimator_cff::width, and SiStripClusterInfo::width().

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), HTMLExport.HTMLExportStatic::export(), and ShallowClustersProducer().

57  {
58  //Retrieve tracker topology from geometry
59  edm::ESHandle<TrackerTopology> tTopoHandle;
60  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
61  const TrackerTopology* const tTopo = tTopoHandle.product();
62 
63  auto number = std::make_unique<std::vector<unsigned>>(7,0);
64  auto width = std::make_unique<std::vector<unsigned>>();
65  auto variance = std::make_unique<std::vector<float>>();
66  auto barystrip = std::make_unique<std::vector<float>>();
67  auto middlestrip = std::make_unique<std::vector<float>>();
68  auto charge = std::make_unique<std::vector<unsigned>>();
69  auto noise = std::make_unique<std::vector<float>>();
70  auto ston = std::make_unique<std::vector<float>>();
71  auto seedstrip = std::make_unique<std::vector<unsigned>>();
72  auto seedindex = std::make_unique<std::vector<unsigned>>();
73  auto seedcharge = std::make_unique<std::vector<unsigned>>();
74  auto seednoise = std::make_unique<std::vector<float>>();
75  auto seedgain = std::make_unique<std::vector<float>>();
76  auto qualityisbad = std::make_unique<std::vector<unsigned>>();
77 
78  auto rawchargeC = std::make_unique<std::vector<float>>();
79  auto rawchargeL = std::make_unique<std::vector<float>>();
80  auto rawchargeR = std::make_unique<std::vector<float>>();
81  auto rawchargeLL = std::make_unique<std::vector<float>>();
82  auto rawchargeRR = std::make_unique<std::vector<float>>();
83  auto etaX = std::make_unique<std::vector<float>>();
84  auto eta = std::make_unique<std::vector<float>>();
85  auto foldedeta = std::make_unique<std::vector<float>>();
86  auto etaasymm = std::make_unique<std::vector<float>>();
87  auto outsideasymm = std::make_unique<std::vector<float>>();
88  auto neweta = std::make_unique<std::vector<float>>();
89  auto newetaerr = std::make_unique<std::vector<float>>();
90 
91  auto detid = std::make_unique<std::vector<unsigned>>();
92  auto subdetid = std::make_unique<std::vector<int>>();
93  auto side = std::make_unique<std::vector<int>>();
94  auto module = std::make_unique<std::vector<int>>();
95  auto layerwheel = std::make_unique<std::vector<int>>();
96  auto stringringrod = std::make_unique<std::vector<int>>();
97  auto petal = std::make_unique<std::vector<int>>();
98  auto stereo = std::make_unique<std::vector<int>>();
99 
100  edm::Handle<edmNew::DetSetVector<SiStripCluster> > clusters;
101  // iEvent.getByLabel(theClustersLabel, clusters);
102  iEvent.getByToken(theClustersToken_, clusters);
103 
104  edm::Handle<edm::DetSetVector<SiStripProcessedRawDigi> > rawProcessedDigis;
105  // iEvent.getByLabel("siStripProcessedRawDigis", "", rawProcessedDigis);
106  iEvent.getByToken(theDigisToken_,rawProcessedDigis);
107 
108  edmNew::DetSetVector<SiStripCluster>::const_iterator itClusters=clusters->begin();
109  for(;itClusters!=clusters->end();++itClusters){
110  uint32_t id = itClusters->id();
111  const moduleVars moduleV(id, tTopo);
112  for(edmNew::DetSet<SiStripCluster>::const_iterator cluster=itClusters->begin(); cluster!=itClusters->end();++cluster){
113 
114  const SiStripClusterInfo info(*cluster, iSetup, id);
115  const NearDigis digis = rawProcessedDigis.isValid() ? NearDigis(info, *rawProcessedDigis) : NearDigis(info);
116 
117  (number->at(0))++;
118  (number->at(moduleV.subdetid))++;
119  width->push_back( cluster->amplitudes().size() );
120  barystrip->push_back( cluster->barycenter() );
121  variance->push_back( info.variance() );
122  middlestrip->push_back( info.firstStrip() + info.width()/2.0 );
123  charge->push_back( info.charge() );
124  noise->push_back( info.noiseRescaledByGain() );
125  ston->push_back( info.signalOverNoise() );
126  seedstrip->push_back( info.maxStrip() );
127  seedindex->push_back( info.maxIndex() );
128  seedcharge->push_back( info.maxCharge() );
129  seednoise->push_back( info.stripNoisesRescaledByGain().at(info.maxIndex()) );
130  seedgain->push_back( info.stripGains().at(info.maxIndex()) );
131  qualityisbad->push_back( info.IsAnythingBad() );
132 
133  rawchargeC->push_back( digis.max );
134  rawchargeL->push_back( digis.left );
135  rawchargeR->push_back( digis.right );
136  rawchargeLL->push_back( digis.Lleft );
137  rawchargeRR->push_back( digis.Rright );
138  etaX->push_back( digis.etaX() );
139  eta->push_back( digis.eta() );
140  etaasymm->push_back( digis.etaasymm() );
141  outsideasymm->push_back( digis.outsideasymm() );
142  neweta->push_back( (digis.last-digis.first)/info.charge() );
143  newetaerr->push_back( (sqrt(digis.last+digis.first))/pow(info.charge(),1.5) );
144 
145  detid->push_back( id );
146  subdetid->push_back( moduleV.subdetid );
147  side->push_back( moduleV.side );
148  module->push_back( moduleV.module );
149  layerwheel->push_back( moduleV.layerwheel );
150  stringringrod->push_back( moduleV.stringringrod );
151  petal->push_back( moduleV.petal );
152  stereo->push_back( moduleV.stereo );
153  }
154  }
155 
156  iEvent.put(std::move(number), Prefix + "number" );
157  iEvent.put(std::move(width), Prefix + "width" );
158  iEvent.put(std::move(variance), Prefix + "variance" );
159  iEvent.put(std::move(barystrip), Prefix + "barystrip" );
160  iEvent.put(std::move(middlestrip), Prefix + "middlestrip" );
161  iEvent.put(std::move(charge), Prefix + "charge" );
162  iEvent.put(std::move(noise), Prefix + "noise" );
163  iEvent.put(std::move(ston), Prefix + "ston" );
164  iEvent.put(std::move(seedstrip), Prefix + "seedstrip" );
165  iEvent.put(std::move(seedindex), Prefix + "seedindex" );
166  iEvent.put(std::move(seedcharge), Prefix + "seedcharge" );
167  iEvent.put(std::move(seednoise), Prefix + "seednoise" );
168  iEvent.put(std::move(seedgain), Prefix + "seedgain" );
169  iEvent.put(std::move(qualityisbad),Prefix + "qualityisbad" );
170 
171  iEvent.put(std::move(rawchargeC), Prefix + "rawchargeC" );
172  iEvent.put(std::move(rawchargeL), Prefix + "rawchargeL" );
173  iEvent.put(std::move(rawchargeR), Prefix + "rawchargeR" );
174  iEvent.put(std::move(rawchargeLL), Prefix + "rawchargeLL" );
175  iEvent.put(std::move(rawchargeRR), Prefix + "rawchargeRR" );
176  iEvent.put(std::move(etaX), Prefix + "etaX" );
177  iEvent.put(std::move(eta), Prefix + "eta" );
178  iEvent.put(std::move(foldedeta), Prefix + "foldedeta" );
179  iEvent.put(std::move(etaasymm), Prefix + "etaasymm" );
180  iEvent.put(std::move(outsideasymm),Prefix + "outsideasymm" );
181  iEvent.put(std::move(neweta), Prefix + "neweta" );
182  iEvent.put(std::move(newetaerr), Prefix + "newetaerr" );
183 
184  iEvent.put(std::move(detid), Prefix + "detid" );
185  iEvent.put(std::move(subdetid), Prefix + "subdetid" );
186  iEvent.put(std::move(module), Prefix + "module" );
187  iEvent.put(std::move(side), Prefix + "side" );
188  iEvent.put(std::move(layerwheel), Prefix + "layerwheel" );
189  iEvent.put(std::move(stringringrod),Prefix + "stringringrod" );
190  iEvent.put(std::move(petal), Prefix + "petal" );
191  iEvent.put(std::move(stereo), Prefix + "stereo" );
192 
193 }
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:220
edmNew::DetSetVector::end
const_iterator end(bool update=false) const
Definition: DetSetVectorNew.h:592
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:170
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:125
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
edmNew::DetSet::const_iterator
data_type const * const_iterator
Definition: DetSetNew.h:30
edmNew::DetSetVector::id
id_type id(size_t cell) const
Definition: DetSetVectorNew.h:630
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
fastPrimaryVertexProducer_cfi.clusters
clusters
Definition: fastPrimaryVertexProducer_cfi.py:5
ApeEstimator_cff.width
width
Definition: ApeEstimator_cff.py:24
ShallowClustersProducer::theClustersToken_
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > > theClustersToken_
Definition: ShallowClustersProducer.h:45
ShallowClustersProducer::theDigisToken_
edm::EDGetTokenT< edm::DetSetVector< SiStripProcessedRawDigi > > theDigisToken_
Definition: ShallowClustersProducer.h:46
edm::EventSetup::get
T get() const
Definition: EventSetup.h:71
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
module
Definition: vlib.h:208
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
eostools.move
def move(src, dest)
Definition: eostools.py:511
edmNew::DetSetVector::begin
const_iterator begin(bool update=false) const
Definition: DetSetVectorNew.h:586

Member Data Documentation

std::string ShallowClustersProducer::Prefix
private

Definition at line 26 of file ShallowClustersProducer.h.

Referenced by produce(), and ShallowClustersProducer().

edm::InputTag ShallowClustersProducer::theClustersLabel
private

Definition at line 25 of file ShallowClustersProducer.h.

edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster> > ShallowClustersProducer::theClustersToken_
private

Definition at line 45 of file ShallowClustersProducer.h.

Referenced by produce(), and ShallowClustersProducer().

edm::EDGetTokenT<edm::DetSetVector<SiStripProcessedRawDigi> > ShallowClustersProducer::theDigisToken_
private

Definition at line 46 of file ShallowClustersProducer.h.

Referenced by produce(), and ShallowClustersProducer().

Generated for CMSSW Reference Manual by doxygen 1.8.11