GEMRecHitStandardAlgo.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author M. Maggi -- INFN
 */
#include "RecoLocalMuon/GEMRecHit/interface/GEMCluster.h"
#include "RecoLocalMuon/GEMRecHit/plugins/GEMRecHitStandardAlgo.h"
#include "DataFormats/MuonDetId/interface/GEMDetId.h"
#include "Geometry/GEMGeometry/interface/GEMEtaPartition.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Utilities/interface/Exception.h"

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

GEMRecHitStandardAlgo::~GEMRecHitStandardAlgo()
{}

void GEMRecHitStandardAlgo::setES(const edm::EventSetup& setup) {}

// First Step
bool GEMRecHitStandardAlgo::compute(const GEMEtaPartition& roll,
                    const GEMCluster& cluster,
                    LocalPoint& Point,
                    LocalError& error)  const
{
  // Get Average Strip position
  float fstrip = (roll.centreOfStrip(cluster.firstStrip())).x();
  float lstrip = (roll.centreOfStrip(cluster.lastStrip())).x();
  float centreOfCluster = (fstrip + lstrip)/2;

  LocalPoint loctemp2(centreOfCluster,0.,0.);
 
  Point = loctemp2;
  error = roll.localError((cluster.firstStrip()+cluster.lastStrip())/2., cluster.clusterSize());
  return true;
}

bool GEMRecHitStandardAlgo::compute(const GEMEtaPartition& roll,
                    const GEMCluster& cl,
                    const float& angle,
                    const GlobalPoint& globPos, 
                    LocalPoint& Point,
                    LocalError& error)  const
{
  // Glob Pos and angle not used so far...
  if (globPos.z()<0){ } // Fake use to avoid warnings
  if (angle<0.){ }      // Fake use to avoid warnings
  this->compute(roll,cl,Point,error);
  return true;
}