#include <ClusterChecker.h>

Public Member Functions
	ClusterChecker (const edm::ParameterSet &conf, edm::ConsumesCollector &iC)

	ClusterChecker (const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)

	ClusterChecker ()=delete

size_t	tooManyClusters (const edm::Event &e) const

	~ClusterChecker ()

Static Public Member Functions
static void	fillDescriptions (edm::ParameterSetDescription &description)

Private Attributes
edm::InputTag	clusterCollectionInputTag_

bool	doACheck_

unsigned int	ignoreDetsAboveNClusters_

unsigned int	maxNrOfCosmicClusters_

unsigned int	maxNrOfPixelClusters_

edm::InputTag	pixelClusterCollectionInputTag_

StringCutObjectSelector< reco::utils::ClusterTotals >	selector_

edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > >	token_pc

edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > >	token_sc

Detailed Description

Definition at line 30 of file ClusterChecker.h.

Constructor & Destructor Documentation

◆ ClusterChecker() [1/3]

ClusterChecker::ClusterChecker	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	iC
	)

Definition at line 12 of file ClusterChecker.cc.

References clusterCollectionInputTag_, edm::ConsumesCollector::consumes(), doACheck_, edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), ignoreDetsAboveNClusters_, maxNrOfCosmicClusters_, maxNrOfPixelClusters_, pixelClusterCollectionInputTag_, token_pc, and token_sc.

     : doACheck_(conf.getParameter<bool>("doClusterCheck")),
       selector_(conf.getParameter<bool>("doClusterCheck") && conf.existsAs<std::string>("cut")
                     ? conf.getParameter<std::string>("cut")
                     : "") {
   if (doACheck_) {
     clusterCollectionInputTag_ = conf.getParameter<edm::InputTag>("ClusterCollectionLabel");
     pixelClusterCollectionInputTag_ = conf.getParameter<edm::InputTag>("PixelClusterCollectionLabel");
     token_sc = iC.consumes<edmNew::DetSetVector<SiStripCluster> >(clusterCollectionInputTag_);
     token_pc = iC.consumes<edmNew::DetSetVector<SiPixelCluster> >(pixelClusterCollectionInputTag_);
     maxNrOfCosmicClusters_ = conf.getParameter<unsigned int>("MaxNumberOfCosmicClusters");
     maxNrOfPixelClusters_ = conf.getParameter<unsigned int>("MaxNumberOfPixelClusters");
     if (conf.existsAs<uint32_t>("DontCountDetsAboveNClusters")) {
       ignoreDetsAboveNClusters_ = conf.getParameter<uint32_t>("DontCountDetsAboveNClusters");
     } else {
       ignoreDetsAboveNClusters_ = 0;
     }
   }
 }

◆ ClusterChecker() [2/3]

ClusterChecker::ClusterChecker	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	iC
	)

Definition at line 10 of file ClusterChecker.cc.

10 : ClusterChecker(conf, iC) {}

ClusterChecker::ClusterChecker

ClusterChecker()=delete

◆ ClusterChecker() [3/3]

ClusterChecker::ClusterChecker ( )

delete

◆ ~ClusterChecker()

ClusterChecker::~ClusterChecker ( )

Definition at line 42 of file ClusterChecker.cc.

42 {}

Member Function Documentation

◆ fillDescriptions()

void ClusterChecker::fillDescriptions ( edm::ParameterSetDescription & description )

static

Definition at line 32 of file ClusterChecker.cc.

References submitPVResolutionJobs::desc, HLT_2022v14_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ClusterCheckerEDProducer::fillDescriptions().

                                                                       {
   desc.add<bool>("doClusterCheck", true);
   desc.add<unsigned>("MaxNumberOfCosmicClusters", 400000);
   desc.add<edm::InputTag>("ClusterCollectionLabel", edm::InputTag("siStripClusters"));
   desc.add<unsigned>("MaxNumberOfPixelClusters", 40000);
   desc.add<edm::InputTag>("PixelClusterCollectionLabel", edm::InputTag("siPixelClusters"));
   desc.add<std::string>("cut",
                         "strip < 400000 && pixel < 40000 && (strip < 50000 + 10*pixel) && (pixel < 5000 + 0.1*strip)");
 }

◆ tooManyClusters()

size_t ClusterChecker::tooManyClusters ( const edm::Event & e ) const

Definition at line 44 of file ClusterChecker.cc.

References doACheck_, MillePedeFileConverter_cfg::e, edm::HandleBase::failedToGet(), ignoreDetsAboveNClusters_, input, maxNrOfCosmicClusters_, maxNrOfPixelClusters_, reco::utils::ClusterTotals::pixel, pixelClusterCollectionInputTag_, reco::utils::ClusterTotals::pixeldets, selector_, edmNew::DetSetVector< T >::size(), reco::utils::ClusterTotals::strip, reco::utils::ClusterTotals::stripdets, token_pc, and token_sc.

Referenced by PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::analyze(), PhotonConversionTrajectorySeedProducerFromSingleLegAlgo::find(), ClusterCheckerEDProducer::produce(), SeedGeneratorFromRegionHitsEDProducer::produce(), CosmicSeedGenerator::produce(), SimpleCosmicBONSeeder::produce(), and CtfSpecialSeedGenerator::produce().

                                                               {
   if (!doACheck_)
     return 0;
 
   // get special input for cosmic cluster multiplicity filter
   edm::Handle<edmNew::DetSetVector<SiStripCluster> > clusterDSV;
   e.getByToken(token_sc, clusterDSV);
   reco::utils::ClusterTotals totals;
   if (!clusterDSV.failedToGet()) {
     const edmNew::DetSetVector<SiStripCluster>& input = *clusterDSV;
 
     if (ignoreDetsAboveNClusters_ == 0) {
       totals.strip = input.dataSize();
       totals.stripdets = input.size();
     } else {
       //loop over detectors
       totals.strip = 0;
       totals.stripdets = 0;
       edmNew::DetSetVector<SiStripCluster>::const_iterator DSViter = input.begin(), DSViter_end = input.end();
       for (; DSViter != DSViter_end; DSViter++) {
         size_t siz = DSViter->size();
         if (siz > ignoreDetsAboveNClusters_)
           continue;
         totals.strip += siz;
         totals.stripdets++;
       }
     }
   }
   if (totals.strip > int(maxNrOfCosmicClusters_))
     return totals.strip;
 
   // get special input for pixel cluster multiplicity filter
   edm::Handle<edmNew::DetSetVector<SiPixelCluster> > pixelClusterDSV;
   e.getByToken(token_pc, pixelClusterDSV);
   if (!pixelClusterDSV.failedToGet()) {
     const edmNew::DetSetVector<SiPixelCluster>& input = *pixelClusterDSV;
 
     if (ignoreDetsAboveNClusters_ == 0) {
       totals.pixel = input.dataSize();
       totals.pixeldets = input.size();
     } else {
       //loop over detectors
       totals.pixel = 0;
       totals.pixeldets = 0;
       edmNew::DetSetVector<SiPixelCluster>::const_iterator DSViter = input.begin(), DSViter_end = input.end();
       for (; DSViter != DSViter_end; DSViter++) {
         size_t siz = DSViter->size();
         if (siz > ignoreDetsAboveNClusters_)
           continue;
         totals.pixel += siz;
         totals.pixeldets++;
       }
     }
   } else {
     //say something's wrong.
     edm::LogError("ClusterChecker")
         << "could not get any SiPixel cluster collections of type edm::DetSetVector<SiPixelCluster>  with label: "
         << pixelClusterCollectionInputTag_;
     totals.pixel = 999999;
   }
   if (totals.pixel > int(maxNrOfPixelClusters_))
     return totals.pixel;
 
   if (!selector_(totals))
     return totals.strip;
   return 0;
 }