#include <ClusterChecker.h>

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

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

	~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_

Detailed Description

Definition at line 19 of file ClusterChecker.h.

Constructor & Destructor Documentation

ClusterChecker::ClusterChecker ( const edm::ParameterSet & conf )

Definition at line 13 of file ClusterChecker.cc.

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

                                                           :
     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");
         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 ( )

Definition at line 33 of file ClusterChecker.cc.

34 {

35 }

Member Function Documentation

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

Definition at line 37 of file ClusterChecker.cc.

References edmNew::DetSetVector< T >::begin(), clusterCollectionInputTag_, edmNew::DetSetVector< T >::dataSize(), doACheck_, edmNew::DetSetVector< T >::end(), edm::HandleBase::failedToGet(), edm::Event::getByLabel(), ignoreDetsAboveNClusters_, LaserDQM_cfg::input, maxNrOfCosmicClusters_, maxNrOfPixelClusters_, reco::utils::ClusterTotals::pixel, pixelClusterCollectionInputTag_, reco::utils::ClusterTotals::pixeldets, selector_, edmNew::DetSetVector< T >::size(), reco::utils::ClusterTotals::strip, and reco::utils::ClusterTotals::stripdets.

Referenced by PhotonConversionTrajectorySeedProducerFromQuadrupletsAlgo::analyze(), PhotonConversionTrajectorySeedProducerFromSingleLegAlgo::analyze(), SeedGeneratorFromRegionHitsEDProducer::produce(), CosmicSeedGenerator::produce(), RoadSearchSeedFinder::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.getByLabel(clusterCollectionInputTag_, 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++;
             }
         }
     }
     else{
         edm::Handle<edm::LazyGetter<SiStripCluster> > lazyGH;
         e.getByLabel(clusterCollectionInputTag_, lazyGH);
         totals.stripdets = 0; // don't know how to count this online
         if (!lazyGH.failedToGet()){
             totals.strip = lazyGH->size();
         }else{
             //say something's wrong.
             edm::LogError("ClusterChecker")<<"could not get any SiStrip cluster collections of type edm::DetSetVector<SiStripCluster> or edm::LazyGetter<SiStripCluster, with label: "<<clusterCollectionInputTag_;
             totals.strip = 999999;
         }
     }
     if (totals.strip > int(maxNrOfCosmicClusters_)) return totals.strip;
 
     // get special input for pixel cluster multiplicity filter
     edm::Handle<edmNew::DetSetVector<SiPixelCluster> > pixelClusterDSV;
     e.getByLabel(pixelClusterCollectionInputTag_, 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;
 }