DTClusterer.h

Go to the documentation of this file.

#ifndef DTCLUSTERER_H
#define DTCLUSTERER_H

/* Base Class Headers */
#include "FWCore/Framework/interface/EDProducer.h"

namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
}

/* Collaborating Class Declarations */
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/DTRecHit/interface/DTSLRecCluster.h"
#include "DataFormats/DTRecHit/interface/DTRecHit1DPair.h"
#include "DataFormats/DTRecHit/interface/DTRecHitCollection.h"
class DTSuperLayer;

/* C++ Headers */
#include <vector>
#include <utility>

/* ====================================================================== */

/* Class DTClusterer Interface */

class DTClusterer : public edm::EDProducer {

  public:

    /* Constructor */ 
    DTClusterer(const edm::ParameterSet&) ;

    /* Destructor */ 
    ~DTClusterer() override ;

    /* Operations */ 
    void produce(edm::Event& event, const edm::EventSetup& setup) override;

  private:
    // build clusters from hits
    std::vector<DTSLRecCluster> buildClusters(const DTSuperLayer* sl,
                                              std::vector<DTRecHit1DPair>& pairs);

    std::vector<std::pair<float, DTRecHit1DPair> > initHits(const DTSuperLayer* sl,
                                                            std::vector<DTRecHit1DPair>& pairs);

    unsigned int differentLayers(std::vector<DTRecHit1DPair>& hits);

  private:
    // to sort hits by x
    struct sortClusterByX {
      bool operator()(const std::pair<float, DTRecHit1DPair>& lhs, 
                      const std::pair<float, DTRecHit1DPair>& rhs) {
        return lhs.first < rhs.first;
      }
    };

  private:
    // Switch on verbosity
    bool debug;

    unsigned int theMinHits; // min number of hits to build a cluster
    unsigned int theMinLayers; // min number of layers to build a cluster
    edm::EDGetTokenT<DTRecHitCollection> recHits1DToken_;
  protected:

};
#endif // DTCLUSTERER_H