/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_7/src/RecoLocalMuon/GEMRecHit/src/GEMRecHitStandardAlgo.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  $Date: 2013/04/24 17:16:35 $
 *  $Revision: 1.1 $
 *  \author M. Maggi -- INFN
 */

#include "GEMCluster.h"
#include "RecoLocalMuon/GEMRecHit/src/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.);
  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;
}