GEMRecHitBaseAlgo.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author M. Maggi -- INFN Bari
 */
#include "RecoLocalMuon/GEMRecHit/interface/GEMRecHitBaseAlgo.h"
#include "RecoLocalMuon/GEMRecHit/interface/GEMClusterizer.h"
#include "RecoLocalMuon/GEMRecHit/interface/GEMMaskReClusterizer.h"

#include "Geometry/GEMGeometry/interface/GEMEtaPartition.h"
#include "DataFormats/GEMDigi/interface/GEMDigiCollection.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

GEMRecHitBaseAlgo::GEMRecHitBaseAlgo(const edm::ParameterSet& config) {}

GEMRecHitBaseAlgo::~GEMRecHitBaseAlgo() {}

// Build all hits in the range associated to the layerId, at the 1st step.
edm::OwnVector<GEMRecHit> GEMRecHitBaseAlgo::reconstruct(const GEMEtaPartition& roll,
                                                         const GEMDetId& gemId,
                                                         const GEMDigiCollection::Range& digiRange,
                                                         const EtaPartitionMask& mask) {
  edm::OwnVector<GEMRecHit> result;

  GEMClusterizer clizer;
  GEMClusterContainer tcls = clizer.doAction(digiRange, mask);
  GEMMaskReClusterizer mrclizer;
  GEMClusterContainer cls = mrclizer.doAction(gemId, tcls, mask);

  for (GEMClusterContainer::const_iterator cl = cls.begin(); cl != cls.end(); cl++) {
    LocalError tmpErr;
    LocalPoint point;
    // Call the compute method
    bool OK = this->compute(roll, *cl, point, tmpErr);
    if (!OK)
      continue;

    // Build a new pair of 1D rechit
    int firstClustStrip = cl->firstStrip();
    int clusterSize = cl->clusterSize();
    GEMRecHit* recHit = new GEMRecHit(gemId, cl->bx(), firstClustStrip, clusterSize, point, tmpErr);

    result.push_back(recHit);
  }
  return result;
}