GEMCSCSegFit Class Reference

#include <GEMCSCSegFit.h>

Public Types
typedef ROOT::Math::SMatrix< double, 16, 4 >	SMatrix16by4
typedef ROOT::Math::SMatrix< double, 4 >	SMatrix4
typedef ROOT::Math::SMatrix< double, 16, 16, ROOT::Math::MatRepSym< double, 16 > >	SMatrixSym16
typedef ROOT::Math::SMatrix< double, 2, 2, ROOT::Math::MatRepSym< double, 2 > >	SMatrixSym2
typedef ROOT::Math::SMatrix< double, 4, 4, ROOT::Math::MatRepSym< double, 4 > >	SMatrixSym4
typedef ROOT::Math::SVector< double, 4 >	SVector4

Public Member Functions
double	chi2 (void) const
AlgebraicSymMatrix	covarianceMatrix (void)
const CSCChamber *	cscchamber (uint32_t id) const
const CSCLayer *	csclayer (uint32_t id) const
void	fit (void)
bool	fitdone () const
	GEMCSCSegFit (std::map< uint32_t, const CSCLayer *> csclayermap, std::map< uint32_t, const GEMEtaPartition *> gemrollmap, const std::vector< const TrackingRecHit *> hits)
const GEMEtaPartition *	gemetapartition (uint32_t id) const
std::vector< const TrackingRecHit * >	hits (void) const
LocalPoint	intercept () const
LocalVector	localdir () const
int	ndof (void) const
size_t	nhits (void) const
float	Rdev (float x, float y, float z) const
const CSCChamber *	refcscchamber () const
double	scaleXError (void) const
void	setScaleXError (double factor)
float	xdev (float x, float z) const
float	xfit (float z) const
float	ydev (float y, float z) const
float	yfit (float z) const
virtual	~GEMCSCSegFit ()

Protected Member Functions
SMatrix16by4	derivativeMatrix (void)
AlgebraicSymMatrix	flipErrors (const SMatrixSym4 &)
void	setOutFromIP (void)
SMatrixSym16	weightMatrix (void)

Protected Attributes
double	chi2_
std::map< uint32_t, const CSCLayer * >	csclayermap_
bool	fitdone_
std::map< uint32_t, const GEMEtaPartition * >	gemetapartmap_
std::vector< const TrackingRecHit * >	hits_
LocalPoint	intercept_
LocalVector	localdir_
int	ndof_
const CSCChamber *	refcscchamber_
uint32_t	refid_
double	scaleXError_
float	uslope_
float	vslope_

Private Member Functions
void	fit2 (void)
void	fitlsq (void)
void	setChi2 (void)

Detailed Description

Definition at line 47 of file GEMCSCSegFit.h.

Member Typedef Documentation

SMatrix16by4

typedef ROOT::Math::SMatrix<double, 16, 4> GEMCSCSegFit::SMatrix16by4

Definition at line 55 of file GEMCSCSegFit.h.

SMatrix4

typedef ROOT::Math::SMatrix<double, 4> GEMCSCSegFit::SMatrix4

Definition at line 58 of file GEMCSCSegFit.h.

SMatrixSym16

typedef ROOT::Math::SMatrix<double, 16, 16, ROOT::Math::MatRepSym<double, 16> > GEMCSCSegFit::SMatrixSym16

Definition at line 52 of file GEMCSCSegFit.h.

SMatrixSym2

typedef ROOT::Math::SMatrix<double, 2, 2, ROOT::Math::MatRepSym<double, 2> > GEMCSCSegFit::SMatrixSym2

Definition at line 62 of file GEMCSCSegFit.h.

SMatrixSym4

typedef ROOT::Math::SMatrix<double, 4, 4, ROOT::Math::MatRepSym<double, 4> > GEMCSCSegFit::SMatrixSym4

Definition at line 59 of file GEMCSCSegFit.h.

SVector4

typedef ROOT::Math::SVector<double, 4> GEMCSCSegFit::SVector4

Definition at line 65 of file GEMCSCSegFit.h.

Constructor & Destructor Documentation

GEMCSCSegFit()

GEMCSCSegFit::GEMCSCSegFit ( std::map< uint32_t, const CSCLayer *>  csclayermap, std::map< uint32_t, const GEMEtaPartition *>  gemrollmap, const std::vector< const TrackingRecHit *>  hits  )

inline

Definition at line 70 of file GEMCSCSegFit.h.

References csclayermap_, gemetapartmap_, and AlCaHLTBitMon_QueryRunRegistry::string.

      : csclayermap_(csclayermap),
        gemetapartmap_(gemrollmap),
        hits_(hits),
        scaleXError_(1.0),
        refid_(csclayermap_.begin()->first),
        fitdone_(false) {
    // --- LogDebug info about reading of CSC Layer map and GEM Eta Partition map -----------------------
    edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit::ctor] cached the csclayermap and the gemrollmap";

    // --- LogDebug for CSC Layer map -------------------------------------------------------------------
    std::stringstream csclayermapss;
    csclayermapss << "[GEMCSCSegFit::ctor] :: csclayermap :: elements [" << std::endl;
    for (std::map<uint32_t, const CSCLayer*>::const_iterator mapIt = csclayermap_.begin(); mapIt != csclayermap_.end();
         ++mapIt) {
      csclayermapss << "[CSC DetId " << mapIt->first << " =" << CSCDetId(mapIt->first) << ", CSC Layer "
                    << mapIt->second << " =" << (mapIt->second)->id() << "]," << std::endl;
    }
    csclayermapss << "]" << std::endl;
    std::string csclayermapstr = csclayermapss.str();
    edm::LogVerbatim("GEMCSCSegFit") << csclayermapstr;
    // --- End LogDebug -----------------------------------------------------------------------------------

    // --- LogDebug for GEM Eta Partition map ------------------------------------------------------------
    std::stringstream gemetapartmapss;
    gemetapartmapss << "[GEMCSCSegFit::ctor] :: gemetapartmap :: elements [" << std::endl;
    for (std::map<uint32_t, const GEMEtaPartition*>::const_iterator mapIt = gemetapartmap_.begin();
         mapIt != gemetapartmap_.end();
         ++mapIt) {
      gemetapartmapss << "[GEM DetId " << mapIt->first << " =" << GEMDetId(mapIt->first) << ", GEM EtaPart "
                      << mapIt->second << "]," << std::endl;
    }
    gemetapartmapss << "]" << std::endl;
    std::string gemetapartmapstr = gemetapartmapss.str();
    edm::LogVerbatim("GEMCSCSegFit") << gemetapartmapstr;
    // --- End LogDebug -----------------------------------------------------------------------------------
  }

~GEMCSCSegFit()

virtual GEMCSCSegFit::~GEMCSCSegFit ( )

inlinevirtual

Definition at line 110 of file GEMCSCSegFit.h.

{}

Member Function Documentation

chi2()

double GEMCSCSegFit::chi2 ( void  ) const

inline

Definition at line 134 of file GEMCSCSegFit.h.

References chi2_.

Referenced by GEMCSCSegAlgoRR::buildSegments().

{ return chi2_; }

covarianceMatrix()

AlgebraicSymMatrix GEMCSCSegFit::covarianceMatrix

(

void

)

Referenced by GEMCSCSegAlgoRR::buildSegments().

cscchamber()

const CSCChamber* GEMCSCSegFit::cscchamber ( uint32_t  id ) const

inline

Definition at line 139 of file GEMCSCSegFit.h.

References relativeConstraints::chamber, csclayermap_, and Exception.

Referenced by fit2(), and fitlsq().

                                                  {
    // Test whether id is found
    if (csclayermap_.find(id) == csclayermap_.end()) {  // id is not found
      throw cms::Exception("InvalidDetId") << "[GEMCSCSegFit] Failed to find CSCChamber in CSCLayerMap" << std::endl;
    }       // chamber is not found and exception is thrown
    else {  // id is found
      return (csclayermap_.find(id)->second)->chamber();
    }  // chamber found and returned
  }

csclayer()

const CSCLayer* GEMCSCSegFit::csclayer ( uint32_t  id ) const

inline

Definition at line 148 of file GEMCSCSegFit.h.

References csclayermap_, and Exception.

                                              {
    if (csclayermap_.find(id) == csclayermap_.end()) {
      throw cms::Exception("InvalidDetId") << "[GEMCSCSegFit] Failed to find CSCLayer in CSCLayerMap" << std::endl;
    } else {
      return csclayermap_.find(id)->second;
    }
  }

derivativeMatrix()

SMatrix16by4 GEMCSCSegFit::derivativeMatrix ( void  )

protected

fit()

void GEMCSCSegFit::fit

(

void

)

Definition at line 19 of file GEMCSCSegFit.cc.

References fit2(), fitdone(), fitlsq(), dqmiodumpmetadata::n, and nhits().

Referenced by GEMCSCSegAlgoRR::buildSegments(), and trackingPlots.Iteration::modules().

                           {
  edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit::fit] - start the fitting fun :: nhits = " << nhits();
  if (fitdone())
    return;  // don't redo fit unnecessarily
  short n = nhits();
  switch (n) {
    case 1:
      edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit::fit] - cannot fit just 1 hit!!";
      break;
    case 2:
      fit2();
      break;
    case 3:
    case 4:
    case 5:
    case 6:
    case 7:
    case 8:
      fitlsq();
      break;
    default:
      edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit::fit] - cannot fit more than 8 hits!!";
  }
}

fit2()

void GEMCSCSegFit::fit2 ( void  )

private

Definition at line 44 of file GEMCSCSegFit.cc.

References chi2_, MuonSubdetId::CSC, cscchamber(), d1, PVValHelper::dz, fitdone_, MuonSubdetId::GEM, gemetapartition(), hits_, intercept_, CSCChamber::layer(), CSCDetId::layer(), phase1PixelTopology::layer, GEMDetId::layer(), TrackingRecHit::localPosition(), ndof_, refcscchamber(), setOutFromIP(), DetId::subdetId(), GeomDet::toGlobal(), GeomDet::toLocal(), uslope_, vslope_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by fit().

                            {
  edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit::fit2] - start fit2()";
  // Just join the two points
  // Equation of straight line between (x1, y1) and (x2, y2) in xy-plane is
  //       y = mx + c
  // with m = (y2-y1)/(x2-x1)
  // and  c = (y1*x2-x2*y1)/(x2-x1)

  // 1) Check whether hits are on the same layer
  // -------------------------------------------
  std::vector<const TrackingRecHit*>::const_iterator ih = hits_.begin();
  // layer numbering: GEM: (1,2) CSC (3,4,5,6,7,8)
  int il1 = 0, il2 = 0;
  DetId d1 = DetId((*ih)->rawId());
  // check whether first hit is GEM or CSC
  if (d1.subdetId() == MuonSubdetId::GEM) {
    il1 = GEMDetId(d1).layer();
  } else if (d1.subdetId() == MuonSubdetId::CSC) {
    il1 = CSCDetId(d1).layer() + 2;
  } else {
    edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit:fit2] - TrackingRecHit is not in GEM or CSC subdetector";
  }
  const TrackingRecHit& h1 = (**ih);
  ++ih;
  DetId d2 = DetId((*ih)->rawId());
  // check whether second hit is GEM or CSC
  if (d2.subdetId() == MuonSubdetId::GEM) {
    il2 = GEMDetId(d2).layer();
  } else if (d2.subdetId() == MuonSubdetId::CSC) {
    il2 = GEMDetId(d2).layer() + 2;
  } else {
    edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit:fit2] - TrackingRecHit is not in GEM or CSC subdetector";
  }
  const TrackingRecHit& h2 = (**ih);
  // Skip if on same layer, but should be impossible :)
  if (il1 == il2) {
    edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit:fit2] - 2 hits on same layer!!";
    return;
  }

  // 2) Global Positions of hit 1 and 2 and
  //    Local  Positions of hit 1 and 2 w.r.t. reference CSC Chamber Frame
  // ---------------------------------------------------------------------
  GlobalPoint h1glopos, h2glopos;
  // global position hit 1
  if (d1.subdetId() == MuonSubdetId::GEM) {
    const GEMEtaPartition* roll1 = gemetapartition(GEMDetId(d1));
    h1glopos = roll1->toGlobal(h1.localPosition());
  } else if (d1.subdetId() == MuonSubdetId::CSC) {
    const CSCLayer* layer1 = cscchamber(CSCDetId(d1))->layer(CSCDetId(d1).layer());
    h1glopos = layer1->toGlobal(h1.localPosition());
  } else {
    edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit:fit2] - TrackingRecHit is not in GEM or CSC subdetector";
  }
  // global position hit 2
  if (d2.subdetId() == MuonSubdetId::GEM) {
    const GEMEtaPartition* roll2 = gemetapartition(GEMDetId(d2));
    h2glopos = roll2->toGlobal(h2.localPosition());
  } else if (d2.subdetId() == MuonSubdetId::CSC) {
    const CSCLayer* layer2 = cscchamber(CSCDetId(d2))->layer(CSCDetId(d2).layer());
    h2glopos = layer2->toGlobal(h2.localPosition());
  } else {
    edm::LogVerbatim("GEMCSCSegFit") << "[GEMCSCSegFit:fit2] - TrackingRecHit is not in GEM or CSC subdetector";
  }
  // local positions hit 1 and 2 w.r.t. ref CSC Chamber Frame
  // We want hit wrt chamber (and local z will be != 0)
  LocalPoint h1pos = refcscchamber()->toLocal(h1glopos);
  LocalPoint h2pos = refcscchamber()->toLocal(h2glopos);

  // 3) Now make straight line between the two points in local coords
  // ----------------------------------------------------------------
  float dz = h2pos.z() - h1pos.z();
  if (dz != 0) {
    uslope_ = (h2pos.x() - h1pos.x()) / dz;
    vslope_ = (h2pos.y() - h1pos.y()) / dz;
  }
  float uintercept = (h1pos.x() * h2pos.z() - h2pos.x() * h1pos.z()) / dz;
  float vintercept = (h1pos.y() * h2pos.z() - h2pos.y() * h1pos.z()) / dz;
  intercept_ = LocalPoint(uintercept, vintercept, 0.);

  setOutFromIP();

  //@@ NOT SURE WHAT IS SENSIBLE FOR THESE...
  chi2_ = 0.;
  ndof_ = 0;

  fitdone_ = true;
}

fitdone()

bool GEMCSCSegFit::fitdone ( ) const

inline

Definition at line 171 of file GEMCSCSegFit.h.

References fitdone_.

Referenced by fit().

{ return fitdone_; }

fitlsq()

void GEMCSCSegFit::fitlsq ( void  )

private

Definition at line 133 of file GEMCSCSegFit.cc.

References B, MuonSubdetId::CSC, cscchamber(), ztail::d, fitdone_, MuonSubdetId::GEM, gemetapartition(), runTauDisplay::gp, hits_, CSCChamber::id(), intercept_, CSCChamber::layer(), phase1PixelTopology::layer, convertSQLiteXML::ok, AlCaHLTBitMon_ParallelJobs::p, DetId::rawId(), refcscchamber(), setChi2(), setOutFromIP(), AlCaHLTBitMon_QueryRunRegistry::string, GeomDet::toGlobal(), GeomDet::toLocal(), uslope_, findQualityFiles::v, vslope_, and z.

Referenced by fit().

flipErrors()

AlgebraicSymMatrix GEMCSCSegFit::flipErrors ( const SMatrixSym4 &  )

protected

gemetapartition()

const GEMEtaPartition* GEMCSCSegFit::gemetapartition ( uint32_t  id ) const

inline

Definition at line 155 of file GEMCSCSegFit.h.

References Exception, and gemetapartmap_.

Referenced by fit2(), and fitlsq().

                                                            {
    if (gemetapartmap_.find(id) == gemetapartmap_.end()) {
      throw cms::Exception("InvalidDetId")
          << "[GEMCSCSegFit] Failed to find GEMEtaPartition in GEMEtaPartMap" << std::endl;
    } else {
      return gemetapartmap_.find(id)->second;
    }
  }

hits()

std::vector<const TrackingRecHit*> GEMCSCSegFit::hits ( void  ) const

inline

Definition at line 131 of file GEMCSCSegFit.h.

References hits_.

{ return hits_; }

intercept()

LocalPoint GEMCSCSegFit::intercept ( ) const

inline

Definition at line 136 of file GEMCSCSegFit.h.

References intercept_.

Referenced by GEMCSCSegAlgoRR::buildSegments().

{ return intercept_; }

localdir()

LocalVector GEMCSCSegFit::localdir ( ) const

inline

Definition at line 137 of file GEMCSCSegFit.h.

References localdir_.

Referenced by GEMCSCSegAlgoRR::buildSegments().

{ return localdir_; }

ndof()

int GEMCSCSegFit::ndof ( void  ) const

inline

Definition at line 135 of file GEMCSCSegFit.h.

References ndof_.

{ return ndof_; }

nhits()

size_t GEMCSCSegFit::nhits ( void  ) const

inline

Definition at line 133 of file GEMCSCSegFit.h.

References hits_.

Referenced by fit().

{ return hits_.size(); }

Rdev()

float GEMCSCSegFit::Rdev	(	float	x,
		float	y,
		float	z
	)		const

refcscchamber()

const CSCChamber* GEMCSCSegFit::refcscchamber ( ) const

inline

Definition at line 163 of file GEMCSCSegFit.h.

References relativeConstraints::chamber, csclayermap_, Exception, and refid_.

Referenced by fit2(), and fitlsq().

                                          {
    if (csclayermap_.find(refid_) == csclayermap_.end()) {
      throw cms::Exception("InvalidDetId")
          << "[GEMCSCSegFit] Failed to find Reference CSCChamber in CSCLayerMap" << std::endl;
    } else {
      return (csclayermap_.find(refid_)->second)->chamber();
    }
  }

scaleXError()

double GEMCSCSegFit::scaleXError ( void  ) const

inline

Definition at line 132 of file GEMCSCSegFit.h.

References scaleXError_.

{ return scaleXError_; }

setChi2()

void GEMCSCSegFit::setChi2 ( void  )

private

Referenced by fitlsq().

setOutFromIP()

void GEMCSCSegFit::setOutFromIP ( void  )

protected

Referenced by fit2(), and fitlsq().

setScaleXError()

void GEMCSCSegFit::setScaleXError ( double  factor )

inline

Definition at line 119 of file GEMCSCSegFit.h.

References DQMScaleToClient_cfi::factor, and scaleXError_.

{ scaleXError_ = factor; }

weightMatrix()

SMatrixSym16 GEMCSCSegFit::weightMatrix ( void  )

protected

xdev()

float GEMCSCSegFit::xdev	(	float	x,
		float	z
	)		const

xfit()

float GEMCSCSegFit::xfit

(

float

z

)

const

ydev()

float GEMCSCSegFit::ydev	(	float	y,
		float	z
	)		const

yfit()

float GEMCSCSegFit::yfit

(

float

z

)

const

Member Data Documentation

chi2_

double GEMCSCSegFit::chi2_

protected

Definition at line 201 of file GEMCSCSegFit.h.

Referenced by chi2(), and fit2().

csclayermap_

std::map<uint32_t, const CSCLayer*> GEMCSCSegFit::csclayermap_

protected

Definition at line 192 of file GEMCSCSegFit.h.

Referenced by cscchamber(), csclayer(), GEMCSCSegFit(), and refcscchamber().

fitdone_

bool GEMCSCSegFit::fitdone_

protected

Definition at line 205 of file GEMCSCSegFit.h.

Referenced by fit2(), fitdone(), and fitlsq().

gemetapartmap_

std::map<uint32_t, const GEMEtaPartition*> GEMCSCSegFit::gemetapartmap_

protected

Definition at line 193 of file GEMCSCSegFit.h.

Referenced by GEMCSCSegFit(), and gemetapartition().

hits_

std::vector<const TrackingRecHit*> GEMCSCSegFit::hits_

protected

Definition at line 196 of file GEMCSCSegFit.h.

Referenced by fit2(), fitlsq(), hits(), and nhits().

intercept_

LocalPoint GEMCSCSegFit::intercept_

protected

Definition at line 199 of file GEMCSCSegFit.h.

Referenced by fit2(), fitlsq(), and intercept().

localdir_

LocalVector GEMCSCSegFit::localdir_

protected

Definition at line 200 of file GEMCSCSegFit.h.

Referenced by localdir().

ndof_

int GEMCSCSegFit::ndof_

protected

Definition at line 202 of file GEMCSCSegFit.h.

Referenced by fit2(), and ndof().

refcscchamber_

const CSCChamber* GEMCSCSegFit::refcscchamber_

protected

Definition at line 194 of file GEMCSCSegFit.h.

refid_

uint32_t GEMCSCSegFit::refid_

protected

Definition at line 204 of file GEMCSCSegFit.h.

Referenced by refcscchamber().

scaleXError_

double GEMCSCSegFit::scaleXError_

protected

Definition at line 203 of file GEMCSCSegFit.h.

Referenced by scaleXError(), and setScaleXError().

uslope_

float GEMCSCSegFit::uslope_

protected

Definition at line 197 of file GEMCSCSegFit.h.

Referenced by fit2(), and fitlsq().

vslope_

float GEMCSCSegFit::vslope_

protected

Definition at line 198 of file GEMCSCSegFit.h.

Referenced by fit2(), and fitlsq().