#include <ClusterChecker.h>

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

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

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

	~ClusterChecker ()

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

Private Member Functions
	ClusterChecker ()

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 24 of file ClusterChecker.h.

Constructor & Destructor Documentation

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

Definition at line 15 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::ClusterChecker	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	iC
	)

Definition at line 10 of file ClusterChecker.cc.

                                :
   ClusterChecker(conf, iC)
 {}

ClusterChecker::~ClusterChecker ( )

Definition at line 47 of file ClusterChecker.cc.

48 {

49 }

ClusterChecker::ClusterChecker ( )

private

Member Function Documentation

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

static

Definition at line 37 of file ClusterChecker.cc.

References edm::ParameterSetDescription::add(), 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)");
 }

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

Definition at line 51 of file ClusterChecker.cc.

References edmNew::DetSetVector< T >::begin(), edmNew::DetSetVector< T >::dataSize(), doACheck_, edmNew::DetSetVector< T >::end(), edm::HandleBase::failedToGet(), edm::Event::getByToken(), 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;
 }