/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/RecoMuon/Navigation/src/MuonNavigableLayer.cc

Go to the documentation of this file.

#include "RecoMuon/Navigation/interface/MuonNavigableLayer.h"

/* Collaborating Class Header */
#include "TrackingTools/DetLayers/interface/DetLayer.h"
#include "TrackingTools/DetLayers/interface/BarrelDetLayer.h"
#include "RecoMuon/Navigation/interface/MuonDetLayerMap.h"
#include "RecoMuon/Navigation/interface/MuonEtaRange.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
/* C++ Headers */
#include <algorithm>


using namespace std;

extern float calculateEta(float r, float z)  {
  if ( z > 0 ) return -log((tan(atan(r/z)/2.)));
  return log(-(tan(atan(r/z)/2.)));

}

MuonEtaRange MuonNavigableLayer::trackingRange(const FreeTrajectoryState& fts) const
{  
  float z = fts.position().z();
  float r = fts.position().perp();
  float eta;
  if ( z>0 ) eta = -log((tan(atan(r/z)/2.)));
  else eta = log(-(tan(atan(r/z)/2.)));

  double theta = atan(r/z);

  double spread = 5.0*sqrt(fts.curvilinearError().matrix()(2,2))/fabs(sin(theta));  //5*sigma(eta)

  //C.L.: this spread could be too large to use.
  // convert it to a smaller one by assuming a virtual radius
  // that transforms the error on angle to error on z axis.
  // not accurate, but works!

  double eta_max = 0;

  if ( z > 0 ) eta_max = calculateEta(r, z+spread); 
  else eta_max = calculateEta(r, z-spread); 

  spread = std::min(0.07, fabs(eta_max-eta));

  MuonEtaRange range(eta+spread,eta-spread);

  spread = 0.07; 
  // special treatment for special geometry in overlap region
  
  if ( eta > 1.0 && eta < 1.1 )  range = MuonEtaRange(eta+3.0*spread,eta-spread);
  if ( eta < -1.0 && eta > -1.1 ) range = MuonEtaRange(eta+spread,eta-3.0*spread);

  return range;
}

bool MuonNavigableLayer::isInsideOut(const FreeTrajectoryState& fts) const {
  
  return (fts.position().basicVector().dot(fts.momentum().basicVector())>0);
  
}