CMS 3D CMS Logo

List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
SeedMultiplicityAnalyzer Class Reference

#include <myTKAnalyses/TrackRecoMonitoring/src/SeedMultiplicityAnalyzer.cc>

Inheritance diagram for SeedMultiplicityAnalyzer:
edm::EDAnalyzer edm::EDConsumerBase

Classes

class  FromTrackRefSeedFilter
 

Public Member Functions

 SeedMultiplicityAnalyzer (const edm::ParameterSet &)
 
 ~SeedMultiplicityAnalyzer () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueueglobalLuminosityBlocksQueue ()
 
SerialTaskQueueglobalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
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 analyze (const edm::Event &, const edm::EventSetup &) override
 

Private Attributes

std::vector< unsigned int > _binseta
 
std::vector< unsigned int > _binsmult
 
std::vector< TH2F * > _hbpixcluslenangle
 
std::vector< TH2F * > _hbpixclusleneta
 
std::vector< TH2F * > _hfpixcluslenangle
 
std::vector< TH2F * > _hfpixclusleneta
 
std::vector< TH1F * > _hpixelrhmult
 
std::vector< TH1F * > _hseedeta
 
std::vector< std::vector< TH2F * > > _hseedeta2D
 
std::vector< TH1F * > _hseedmult
 
std::vector< std::vector< TH2F * > > _hseedmult2D
 
std::vector< TH2F * > _hseedphieta
 
std::vector< std::string > _labels
 
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord_magFieldToken
 
std::vector< double > _maxs
 
std::vector< edm::EDGetTokenT< std::map< unsigned int, int > > > _multiplicityMapTokens
 
std::vector< unsigned int > _seedbins
 
std::vector< edm::EDGetTokenT< TrajectorySeedCollection > > _seedcollTokens
 
std::vector< FromTrackRefSeedFilter_seedfilters
 
std::vector< double > _seedmax
 
std::vector< unsigned int > _selections
 
edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord_TTRHBuilderToken
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsInputProcessBlocks ()
 
static bool wantsProcessBlocks ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
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<Transition Tr = Transition::Event>
constexpr auto esConsumes () noexcept
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag) noexcept
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 69 of file SeedMultiplicityAnalyzer.cc.

Constructor & Destructor Documentation

◆ SeedMultiplicityAnalyzer()

SeedMultiplicityAnalyzer::SeedMultiplicityAnalyzer ( const edm::ParameterSet iConfig)
explicit

Definition at line 131 of file SeedMultiplicityAnalyzer.cc.

132  : _magFieldToken(esConsumes(edm::ESInputTag{"", iConfig.getParameter<std::string>("TTRHBuilder")})),
134  _seedcollTokens(),
135  _seedbins(),
136  _seedmax(),
137  _seedfilters(),
139  _labels(),
140  _selections(),
141  _binsmult(),
142  _binseta(),
143  _maxs(),
144  _hseedmult2D(),
145  _hseedeta2D() {
146  //now do what ever initialization is needed
147 
148  //
149 
150  std::vector<edm::ParameterSet> seedCollectionConfigs =
151  iConfig.getParameter<std::vector<edm::ParameterSet>>("seedCollections");
152 
153  for (std::vector<edm::ParameterSet>::const_iterator scps = seedCollectionConfigs.begin();
154  scps != seedCollectionConfigs.end();
155  ++scps) {
156  _seedcollTokens.push_back(consumes<TrajectorySeedCollection>(scps->getParameter<edm::InputTag>("src")));
157  _seedbins.push_back(scps->getUntrackedParameter<unsigned int>("nBins", 1000));
158  _seedmax.push_back(scps->getUntrackedParameter<double>("maxValue", 100000.));
159 
160  if (scps->exists("trackFilter")) {
161  _seedfilters.push_back(
162  FromTrackRefSeedFilter(consumesCollector(), scps->getParameter<edm::ParameterSet>("trackFilter")));
163  } else {
164  _seedfilters.push_back(FromTrackRefSeedFilter());
165  }
166  }
167 
168  std::vector<edm::ParameterSet> correlationConfigs =
169  iConfig.getParameter<std::vector<edm::ParameterSet>>("multiplicityCorrelations");
170 
171  for (std::vector<edm::ParameterSet>::const_iterator ps = correlationConfigs.begin(); ps != correlationConfigs.end();
172  ++ps) {
173  _multiplicityMapTokens.push_back(
174  consumes<std::map<unsigned int, int>>(ps->getParameter<edm::InputTag>("multiplicityMap")));
175  _labels.push_back(ps->getParameter<std::string>("detLabel"));
176  _selections.push_back(ps->getParameter<unsigned int>("detSelection"));
177  _binsmult.push_back(ps->getParameter<unsigned int>("nBins"));
178  _binseta.push_back(ps->getParameter<unsigned int>("nBinsEta"));
179  _maxs.push_back(ps->getParameter<double>("maxValue"));
180  }
181 
183 
184  std::vector<unsigned int>::const_iterator nseedbins = _seedbins.begin();
185  std::vector<double>::const_iterator seedmax = _seedmax.begin();
186  std::vector<FromTrackRefSeedFilter>::const_iterator filter = _seedfilters.begin();
187 
188  for (std::vector<edm::ParameterSet>::const_iterator scps = seedCollectionConfigs.begin();
189  scps != seedCollectionConfigs.end();
190  ++scps, ++nseedbins, ++seedmax, ++filter) {
191  std::string extendedlabel = std::string(scps->getParameter<edm::InputTag>("src").encode()) + filter->suffix();
192 
193  std::string hname = extendedlabel + std::string("_mult");
194  std::string htitle = extendedlabel + std::string(" seed multiplicity");
195  _hseedmult.push_back(tfserv->make<TH1F>(
196  hname.c_str(), htitle.c_str(), *nseedbins + 1, 0.5 - *seedmax / (*nseedbins), *seedmax + 0.5));
197  _hseedmult[_hseedmult.size() - 1]->GetXaxis()->SetTitle("seeds");
198  _hseedmult[_hseedmult.size() - 1]->GetYaxis()->SetTitle("events");
199 
200  hname = extendedlabel + std::string("_eta");
201  htitle = extendedlabel + std::string(" seed pseudorapidity");
202  _hseedeta.push_back(tfserv->make<TH1F>(hname.c_str(), htitle.c_str(), 80, -4., 4.));
203  _hseedeta[_hseedeta.size() - 1]->GetXaxis()->SetTitle("#eta");
204  _hseedeta[_hseedeta.size() - 1]->GetYaxis()->SetTitle("seeds");
205 
206  hname = extendedlabel + std::string("_phieta");
207  htitle = extendedlabel + std::string(" seed phi vs pseudorapidity");
208  _hseedphieta.push_back(tfserv->make<TH2F>(hname.c_str(), htitle.c_str(), 80, -4., 4., 80, -M_PI, M_PI));
209  _hseedphieta[_hseedphieta.size() - 1]->GetXaxis()->SetTitle("#eta");
210  _hseedphieta[_hseedphieta.size() - 1]->GetYaxis()->SetTitle("#phi");
211 
212  _hseedmult2D.push_back(std::vector<TH2F*>());
213  _hseedeta2D.push_back(std::vector<TH2F*>());
214 
215  hname = extendedlabel + std::string("_npixelrh");
216  htitle = extendedlabel + std::string(" seed SiPixelRecHit multiplicity");
217  _hpixelrhmult.push_back(tfserv->make<TH1F>(hname.c_str(), htitle.c_str(), 5, -.5, 4.5));
218  _hpixelrhmult[_hpixelrhmult.size() - 1]->GetXaxis()->SetTitle("NRecHits");
219  _hpixelrhmult[_hpixelrhmult.size() - 1]->GetYaxis()->SetTitle("seeds");
220 
221  hname = extendedlabel + std::string("_bpixleneta");
222  htitle = extendedlabel + std::string(" seed BPIX cluster length vs pseudorapidity");
223  _hbpixclusleneta.push_back(tfserv->make<TH2F>(hname.c_str(), htitle.c_str(), 80, -4., 4., 40, -0.5, 39.5));
224  _hbpixclusleneta[_hbpixclusleneta.size() - 1]->GetXaxis()->SetTitle("#eta");
225  _hbpixclusleneta[_hbpixclusleneta.size() - 1]->GetYaxis()->SetTitle("length");
226 
227  hname = extendedlabel + std::string("_fpixleneta");
228  htitle = extendedlabel + std::string(" seed FPIX cluster length vs pseudorapidity");
229  _hfpixclusleneta.push_back(tfserv->make<TH2F>(hname.c_str(), htitle.c_str(), 80, -4., 4., 40, -0.5, 39.5));
230  _hfpixclusleneta[_hfpixclusleneta.size() - 1]->GetXaxis()->SetTitle("#eta");
231  _hfpixclusleneta[_hfpixclusleneta.size() - 1]->GetYaxis()->SetTitle("length");
232 
233  hname = extendedlabel + std::string("_bpixlenangle");
234  htitle = extendedlabel + std::string(" seed BPIX cluster length vs track projection");
235  _hbpixcluslenangle.push_back(tfserv->make<TH2F>(hname.c_str(), htitle.c_str(), 200, -1., 1., 40, -0.5, 39.5));
236  _hbpixcluslenangle[_hbpixcluslenangle.size() - 1]->GetXaxis()->SetTitle("projection");
237  _hbpixcluslenangle[_hbpixcluslenangle.size() - 1]->GetYaxis()->SetTitle("length");
238 
239  hname = extendedlabel + std::string("_fpixlenangle");
240  htitle = extendedlabel + std::string(" seed FPIX cluster length vs track projection");
241  _hfpixcluslenangle.push_back(tfserv->make<TH2F>(hname.c_str(), htitle.c_str(), 200, -1., 1., 40, -0.5, 39.5));
242  _hfpixcluslenangle[_hfpixcluslenangle.size() - 1]->GetXaxis()->SetTitle("projection");
243  _hfpixcluslenangle[_hfpixcluslenangle.size() - 1]->GetYaxis()->SetTitle("length");
244 
245  for (unsigned int i = 0; i < _multiplicityMapTokens.size(); ++i) {
246  std::string hname2D = extendedlabel + _labels[i];
247  hname2D += "_mult";
248  std::string htitle2D = extendedlabel + " seeds multiplicity";
249  htitle2D += " vs ";
250  htitle2D += _labels[i];
251  htitle2D += " hits";
252  _hseedmult2D[_hseedmult2D.size() - 1].push_back(tfserv->make<TH2F>(hname2D.c_str(),
253  htitle2D.c_str(),
254  _binsmult[i],
255  0.,
256  _maxs[i],
257  *nseedbins + 1,
258  0.5 - *seedmax / (*nseedbins),
259  *seedmax + 0.5));
260  _hseedmult2D[_hseedmult2D.size() - 1][_hseedmult2D[_hseedmult2D.size() - 1].size() - 1]->GetXaxis()->SetTitle(
261  "hits");
262  _hseedmult2D[_hseedmult2D.size() - 1][_hseedmult2D[_hseedmult2D.size() - 1].size() - 1]->GetYaxis()->SetTitle(
263  "seeds");
264 
265  hname2D = extendedlabel + _labels[i];
266  hname2D += "_eta";
267  htitle2D = extendedlabel + " seeds pseudorapidity";
268  htitle2D += " vs ";
269  htitle2D += _labels[i];
270  htitle2D += " hits";
271  _hseedeta2D[_hseedeta2D.size() - 1].push_back(
272  tfserv->make<TH2F>(hname2D.c_str(), htitle2D.c_str(), _binseta[i], 0., _maxs[i], 80, -4., 4.));
273  _hseedeta2D[_hseedeta2D.size() - 1][_hseedeta2D[_hseedeta2D.size() - 1].size() - 1]->GetXaxis()->SetTitle("hits");
274  _hseedeta2D[_hseedeta2D.size() - 1][_hseedeta2D[_hseedeta2D.size() - 1].size() - 1]->GetYaxis()->SetTitle("#eta");
275  }
276  }
277 }

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~SeedMultiplicityAnalyzer()

SeedMultiplicityAnalyzer::~SeedMultiplicityAnalyzer ( )
override

Definition at line 279 of file SeedMultiplicityAnalyzer.cc.

279  {
280  // do anything here that needs to be done at desctruction time
281  // (e.g. close files, deallocate resources etc.)
282 }

Member Function Documentation

◆ analyze()

void SeedMultiplicityAnalyzer::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 289 of file SeedMultiplicityAnalyzer.cc.

289  {
290  using namespace edm;
291 
292  // compute cluster multiplicities
293 
294  std::vector<int> tmpmult(_multiplicityMapTokens.size(), -1);
295  for (unsigned int i = 0; i < _multiplicityMapTokens.size(); ++i) {
297  iEvent.getByToken(_multiplicityMapTokens[i], mults);
298 
299  // check if the selection exists
300 
301  std::map<unsigned int, int>::const_iterator mult = mults->find(_selections[i]);
302 
303  if (mult != mults->end()) {
304  tmpmult[i] = mult->second;
305  } else {
306  edm::LogWarning("DetSelectionNotFound") << " DetSelection " << _selections[i] << " not found";
307  }
308  }
309 
310  // preparation for loop on seeds
311 
312  // TrajectoryStateTransform tsTransform;
313  TSCBLBuilderNoMaterial tscblBuilder; // I could have used TSCBLBuilderWithPropagator
314 
315  const auto theMF = &iSetup.getData(_magFieldToken);
316  const auto& theTTRHBuilder = iSetup.getData(_TTRHBuilderToken);
317 
318  // I need:
319  // - beamspot bs POSTPONED
320 
321  std::vector<TH1F*>::iterator histomult = _hseedmult.begin();
322  std::vector<std::vector<TH2F*>>::iterator histomult2D = _hseedmult2D.begin();
323  std::vector<TH1F*>::iterator histoeta = _hseedeta.begin();
324  std::vector<TH2F*>::iterator histophieta = _hseedphieta.begin();
325  std::vector<std::vector<TH2F*>>::iterator histoeta2D = _hseedeta2D.begin();
326  std::vector<TH1F*>::iterator hpixelrhmult = _hpixelrhmult.begin();
327  std::vector<TH2F*>::iterator histobpixleneta = _hbpixclusleneta.begin();
328  std::vector<TH2F*>::iterator histofpixleneta = _hfpixclusleneta.begin();
329  std::vector<TH2F*>::iterator histobpixlenangle = _hbpixcluslenangle.begin();
330  std::vector<TH2F*>::iterator histofpixlenangle = _hfpixcluslenangle.begin();
331  std::vector<FromTrackRefSeedFilter>::iterator filter = _seedfilters.begin();
332 
333  // loop on seed collections
334 
335  for (std::vector<edm::EDGetTokenT<TrajectorySeedCollection>>::const_iterator coll = _seedcollTokens.begin();
336  coll != _seedcollTokens.end() && histomult != _hseedmult.end() && histomult2D != _hseedmult2D.end() &&
337  histoeta != _hseedeta.end() && histoeta2D != _hseedeta2D.end() && histophieta != _hseedphieta.end() &&
338  hpixelrhmult != _hpixelrhmult.end() && histobpixleneta != _hbpixclusleneta.end() &&
339  histofpixleneta != _hfpixclusleneta.end() && histobpixlenangle != _hbpixcluslenangle.end() &&
340  histofpixlenangle != _hfpixcluslenangle.end();
341  ++coll,
342  ++histomult,
343  ++histomult2D,
344  ++histoeta,
345  ++histophieta,
346  ++histoeta2D,
347  ++hpixelrhmult,
348  ++histobpixleneta,
349  ++histofpixleneta,
350  ++histobpixlenangle,
351  ++histofpixlenangle,
352  ++filter) {
353  filter->prepareEvent(iEvent);
354 
356  iEvent.getByToken(*coll, seeds);
357 
358  /*
359  (*histomult)->Fill(seeds->size());
360 
361  for(unsigned int i=0;i<_multiplicityMaps.size();++i) {
362  if(tmpmult[i]>=0) (*histomult2D)[i]->Fill(tmpmult[i],seeds->size());
363  }
364  */
365 
366  // loop on seeds
367 
368  unsigned int nseeds = 0;
369  unsigned int iseed = 0;
370  for (TrajectorySeedCollection::const_iterator seed = seeds->begin(); seed != seeds->end(); ++seed, ++iseed) {
371  if (filter->isSelected(iseed)) {
372  ++nseeds;
373 
374  TransientTrackingRecHit::RecHitPointer recHit = theTTRHBuilder.build(&*(seed->recHits().end() - 1));
376  trajectoryStateTransform::transientState(seed->startingState(), recHit->surface(), theMF);
377  // TrajectoryStateClosestToBeamLine tsAtClosestApproachSeed = tscblBuilder(*state.freeState(),bs); // here I need them BS
378 
379  double eta = state.globalMomentum().eta();
380  double phi = state.globalMomentum().phi();
381 
382  (*histoeta)->Fill(eta);
383  (*histophieta)->Fill(eta, phi);
384 
385  for (unsigned int i = 0; i < _multiplicityMapTokens.size(); ++i) {
386  if (tmpmult[i] >= 0)
387  (*histoeta2D)[i]->Fill(tmpmult[i], eta);
388  }
389 
390  int npixelrh = 0;
391  for (auto const& hit : seed->recHits()) {
392  const SiPixelRecHit* sphit = dynamic_cast<const SiPixelRecHit*>(&hit);
393  if (sphit) {
394  ++npixelrh;
395  // compute state on recHit surface
396  TransientTrackingRecHit::RecHitPointer ttrhit = theTTRHBuilder.build(&hit);
398  trajectoryStateTransform::transientState(seed->startingState(), ttrhit->surface(), theMF);
399 
400  if (sphit->geographicalId().det() == DetId::Tracker &&
402  (*histobpixleneta)->Fill(eta, sphit->cluster()->sizeY());
403  if (tsos.isValid()) {
404  // double normdx = sin(atan2(tsos.localMomentum().x(),tsos.localMomentum().z()));
405  double normdx = tsos.localMomentum().x() / sqrt(tsos.localMomentum().x() * tsos.localMomentum().x() +
406  tsos.localMomentum().z() * tsos.localMomentum().z());
407  (*histobpixlenangle)->Fill(normdx, sphit->cluster()->sizeY());
408  }
409  } else if (sphit->geographicalId().det() == DetId::Tracker &&
411  (*histofpixleneta)->Fill(eta, sphit->cluster()->sizeX());
412  if (tsos.isValid()) {
413  // double normdy = sin(atan2(tsos.localMomentum().y(),tsos.localMomentum().z()));
414  double normdy = tsos.localMomentum().y() / sqrt(tsos.localMomentum().y() * tsos.localMomentum().y() +
415  tsos.localMomentum().z() * tsos.localMomentum().z());
416  (*histofpixlenangle)->Fill(normdy, sphit->cluster()->sizeX());
417  }
418  } else {
419  edm::LogError("InconsistentSiPixelRecHit")
420  << "SiPixelRecHit with a non-pixel DetId " << sphit->geographicalId().rawId();
421  }
422  }
423  }
424  (*hpixelrhmult)->Fill(npixelrh);
425  }
426  }
427  (*histomult)->Fill(nseeds);
428 
429  for (unsigned int i = 0; i < _multiplicityMapTokens.size(); ++i) {
430  if (tmpmult[i] >= 0)
431  (*histomult2D)[i]->Fill(tmpmult[i], nseeds);
432  }
433  }
434 }

References _hbpixcluslenangle, _hbpixclusleneta, _hfpixcluslenangle, _hfpixclusleneta, _hpixelrhmult, _hseedeta, _hseedeta2D, _hseedmult, _hseedmult2D, _hseedphieta, _magFieldToken, _multiplicityMapTokens, _seedcollTokens, _seedfilters, _selections, _TTRHBuilderToken, SiPixelRecHit::cluster(), DetId::det(), PVValHelper::eta, ALCARECOTkAlBeamHalo_cff::filter, TrackingRecHit::geographicalId(), edm::EventSetup::getData(), mps_fire::i, iEvent, iseed, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localMomentum(), VarParsing::mult, phi, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::rawId(), rpcPointValidation_cfi::recHit, fileCollector::seed, DetachedQuadStep_cff::seeds, mathSSE::sqrt(), DetId::subdetId(), DetId::Tracker, trajectoryStateTransform::transientState(), trackerHitRTTI::vector, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Member Data Documentation

◆ _binseta

std::vector<unsigned int> SeedMultiplicityAnalyzer::_binseta
private

Definition at line 106 of file SeedMultiplicityAnalyzer.cc.

◆ _binsmult

std::vector<unsigned int> SeedMultiplicityAnalyzer::_binsmult
private

Definition at line 105 of file SeedMultiplicityAnalyzer.cc.

◆ _hbpixcluslenangle

std::vector<TH2F*> SeedMultiplicityAnalyzer::_hbpixcluslenangle
private

Definition at line 114 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hbpixclusleneta

std::vector<TH2F*> SeedMultiplicityAnalyzer::_hbpixclusleneta
private

Definition at line 112 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hfpixcluslenangle

std::vector<TH2F*> SeedMultiplicityAnalyzer::_hfpixcluslenangle
private

Definition at line 115 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hfpixclusleneta

std::vector<TH2F*> SeedMultiplicityAnalyzer::_hfpixclusleneta
private

Definition at line 113 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hpixelrhmult

std::vector<TH1F*> SeedMultiplicityAnalyzer::_hpixelrhmult
private

Definition at line 111 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hseedeta

std::vector<TH1F*> SeedMultiplicityAnalyzer::_hseedeta
private

Definition at line 109 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hseedeta2D

std::vector<std::vector<TH2F*> > SeedMultiplicityAnalyzer::_hseedeta2D
private

Definition at line 117 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hseedmult

std::vector<TH1F*> SeedMultiplicityAnalyzer::_hseedmult
private

Definition at line 108 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hseedmult2D

std::vector<std::vector<TH2F*> > SeedMultiplicityAnalyzer::_hseedmult2D
private

Definition at line 116 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _hseedphieta

std::vector<TH2F*> SeedMultiplicityAnalyzer::_hseedphieta
private

Definition at line 110 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _labels

std::vector<std::string> SeedMultiplicityAnalyzer::_labels
private

◆ _magFieldToken

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> SeedMultiplicityAnalyzer::_magFieldToken
private

Definition at line 96 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _maxs

std::vector<double> SeedMultiplicityAnalyzer::_maxs
private

Definition at line 107 of file SeedMultiplicityAnalyzer.cc.

◆ _multiplicityMapTokens

std::vector<edm::EDGetTokenT<std::map<unsigned int, int> > > SeedMultiplicityAnalyzer::_multiplicityMapTokens
private

Definition at line 102 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _seedbins

std::vector<unsigned int> SeedMultiplicityAnalyzer::_seedbins
private

Definition at line 99 of file SeedMultiplicityAnalyzer.cc.

◆ _seedcollTokens

std::vector<edm::EDGetTokenT<TrajectorySeedCollection> > SeedMultiplicityAnalyzer::_seedcollTokens
private

Definition at line 98 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _seedfilters

std::vector<FromTrackRefSeedFilter> SeedMultiplicityAnalyzer::_seedfilters
private

Definition at line 101 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _seedmax

std::vector<double> SeedMultiplicityAnalyzer::_seedmax
private

Definition at line 100 of file SeedMultiplicityAnalyzer.cc.

◆ _selections

std::vector<unsigned int> SeedMultiplicityAnalyzer::_selections
private

Definition at line 104 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

◆ _TTRHBuilderToken

edm::ESGetToken<TransientTrackingRecHitBuilder, TransientRecHitRecord> SeedMultiplicityAnalyzer::_TTRHBuilderToken
private

Definition at line 97 of file SeedMultiplicityAnalyzer.cc.

Referenced by analyze().

SeedMultiplicityAnalyzer::_binsmult
std::vector< unsigned int > _binsmult
Definition: SeedMultiplicityAnalyzer.cc:105
TrajectoryStateOnSurface::localMomentum
LocalVector localMomentum() const
Definition: TrajectoryStateOnSurface.h:75
mps_fire.i
i
Definition: mps_fire.py:428
edm::ESInputTag
Definition: ESInputTag.h:87
PixelSubdetector::PixelEndcap
Definition: PixelSubdetector.h:11
SeedMultiplicityAnalyzer::_seedbins
std::vector< unsigned int > _seedbins
Definition: SeedMultiplicityAnalyzer.cc:99
PixelSubdetector::PixelBarrel
Definition: PixelSubdetector.h:11
SeedMultiplicityAnalyzer::_seedfilters
std::vector< FromTrackRefSeedFilter > _seedfilters
Definition: SeedMultiplicityAnalyzer.cc:101
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
edm::EDGetTokenT< TrajectorySeedCollection >
edm
HLT enums.
Definition: AlignableModifier.h:19
DetId::det
constexpr Detector det() const
get the detector field from this detid
Definition: DetId.h:46
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:47
SeedMultiplicityAnalyzer::_multiplicityMapTokens
std::vector< edm::EDGetTokenT< std::map< unsigned int, int > > > _multiplicityMapTokens
Definition: SeedMultiplicityAnalyzer.cc:102
SeedMultiplicityAnalyzer::_selections
std::vector< unsigned int > _selections
Definition: SeedMultiplicityAnalyzer.cc:104
TrackingRecHit::RecHitPointer
std::shared_ptr< TrackingRecHit const > RecHitPointer
Definition: TrackingRecHit.h:24
edm::Handle
Definition: AssociativeIterator.h:50
SiPixelRecHit::cluster
ClusterRef cluster() const
Definition: SiPixelRecHit.h:47
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:120
rpcPointValidation_cfi.recHit
recHit
Definition: rpcPointValidation_cfi.py:7
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
SeedMultiplicityAnalyzer::_TTRHBuilderToken
edm::ESGetToken< TransientTrackingRecHitBuilder, TransientRecHitRecord > _TTRHBuilderToken
Definition: SeedMultiplicityAnalyzer.cc:97
SeedMultiplicityAnalyzer::_hseedmult
std::vector< TH1F * > _hseedmult
Definition: SeedMultiplicityAnalyzer.cc:108
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
fileCollector.seed
seed
Definition: fileCollector.py:127
TrajectoryStateOnSurface
Definition: TrajectoryStateOnSurface.h:16
SeedMultiplicityAnalyzer::_seedcollTokens
std::vector< edm::EDGetTokenT< TrajectorySeedCollection > > _seedcollTokens
Definition: SeedMultiplicityAnalyzer.cc:98
SeedMultiplicityAnalyzer::_hfpixcluslenangle
std::vector< TH2F * > _hfpixcluslenangle
Definition: SeedMultiplicityAnalyzer.cc:115
PVValHelper::eta
Definition: PVValidationHelpers.h:70
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
SeedMultiplicityAnalyzer::_hpixelrhmult
std::vector< TH1F * > _hpixelrhmult
Definition: SeedMultiplicityAnalyzer.cc:111
SeedMultiplicityAnalyzer::_hfpixclusleneta
std::vector< TH2F * > _hfpixclusleneta
Definition: SeedMultiplicityAnalyzer.cc:113
ALCARECOTkAlBeamHalo_cff.filter
filter
Definition: ALCARECOTkAlBeamHalo_cff.py:27
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
Definition: DetId.h:48
edm::ParameterSet
Definition: ParameterSet.h:47
SeedMultiplicityAnalyzer::_magFieldToken
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > _magFieldToken
Definition: SeedMultiplicityAnalyzer.cc:96
iseed
int iseed
Definition: AMPTWrapper.h:134
DetId::Tracker
Definition: DetId.h:25
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
edm::Service< TFileService >
iEvent
int iEvent
Definition: GenABIO.cc:224
trajectoryStateTransform::transientState
TrajectoryStateOnSurface transientState(const PTrajectoryStateOnDet &ts, const Surface *surface, const MagneticField *field)
Definition: TrajectoryStateTransform.cc:35
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
trackerHitRTTI::vector
Definition: trackerHitRTTI.h:21
edm::InputTag::encode
std::string encode() const
Definition: InputTag.cc:159
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:127
DDAxes::phi
SeedMultiplicityAnalyzer::_seedmax
std::vector< double > _seedmax
Definition: SeedMultiplicityAnalyzer.cc:100
SeedMultiplicityAnalyzer::_hseedeta
std::vector< TH1F * > _hseedeta
Definition: SeedMultiplicityAnalyzer.cc:109
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
DetachedQuadStep_cff.seeds
seeds
Definition: DetachedQuadStep_cff.py:195
SeedMultiplicityAnalyzer::_binseta
std::vector< unsigned int > _binseta
Definition: SeedMultiplicityAnalyzer.cc:106
RunInfoPI::state
state
Definition: RunInfoPayloadInspectoHelper.h:16
TSCBLBuilderNoMaterial
Definition: TSCBLBuilderNoMaterial.h:13
SeedMultiplicityAnalyzer::_hseedmult2D
std::vector< std::vector< TH2F * > > _hseedmult2D
Definition: SeedMultiplicityAnalyzer.cc:116
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
SeedMultiplicityAnalyzer::_maxs
std::vector< double > _maxs
Definition: SeedMultiplicityAnalyzer.cc:107
VarParsing.mult
mult
Definition: VarParsing.py:658
edm::EDConsumerBase::esConsumes
auto esConsumes()
Definition: EDConsumerBase.h:206
SeedMultiplicityAnalyzer::_hseedeta2D
std::vector< std::vector< TH2F * > > _hseedeta2D
Definition: SeedMultiplicityAnalyzer.cc:117
SeedMultiplicityAnalyzer::_labels
std::vector< std::string > _labels
Definition: SeedMultiplicityAnalyzer.cc:103
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:159
SeedMultiplicityAnalyzer::_hbpixclusleneta
std::vector< TH2F * > _hbpixclusleneta
Definition: SeedMultiplicityAnalyzer.cc:112
SeedMultiplicityAnalyzer::_hbpixcluslenangle
std::vector< TH2F * > _hbpixcluslenangle
Definition: SeedMultiplicityAnalyzer.cc:114
edm::InputTag
Definition: InputTag.h:15
SeedMultiplicityAnalyzer::_hseedphieta
std::vector< TH2F * > _hseedphieta
Definition: SeedMultiplicityAnalyzer.cc:110
TrajectoryStateOnSurface::isValid
bool isValid() const
Definition: TrajectoryStateOnSurface.h:54
hit
Definition: SiStripHitEffFromCalibTree.cc:88