#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.

       : 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 -----------------------------------------------------------------------------------
   }

References csclayermap_, gemetapartmap_, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~GEMCSCSegFit()

virtual GEMCSCSegFit::~GEMCSCSegFit ( )

inlinevirtual

Definition at line 110 of file GEMCSCSegFit.h.

110 {}

Member Function Documentation

◆ chi2()

double GEMCSCSegFit::chi2 ( void ) const

inline

Definition at line 134 of file GEMCSCSegFit.h.

134 { return chi2_; }

References chi2_.

Referenced by GEMCSCSegAlgoRR::buildSegments().

◆ 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.

                                                   {
     // 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
   }

References relativeConstraints::chamber, csclayermap_, and Exception.

Referenced by fit2(), and fitlsq().

◆ csclayer()

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

inline

Definition at line 148 of file GEMCSCSegFit.h.

                                               {
     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;
     }
   }

References csclayermap_, and Exception.

◆ derivativeMatrix()

SMatrix16by4 GEMCSCSegFit::derivativeMatrix ( void )

protected

◆ fit()

void GEMCSCSegFit::fit ( void )

Definition at line 19 of file GEMCSCSegFit.cc.

                            {
   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!!";
   }
 }

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

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

◆ fit2()

void GEMCSCSegFit::fit2 ( void )

private

Definition at line 44 of file GEMCSCSegFit.cc.

                             {
   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;
 }

References chi2_, MuonSubdetId::CSC, cscchamber(), d1, PVValHelper::dz, fitdone_, MuonSubdetId::GEM, gemetapartition(), hits_, intercept_, CSCChamber::layer(), CSCDetId::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().

◆ fitdone()

bool GEMCSCSegFit::fitdone ( ) const

inline

Definition at line 171 of file GEMCSCSegFit.h.

171 { return fitdone_; }

References fitdone_.

Referenced by fit().

◆ fitlsq()

void GEMCSCSegFit::fitlsq ( void )

private

Definition at line 133 of file GEMCSCSegFit.cc.

References TtFullHadDaughter::B, MuonSubdetId::CSC, cscchamber(), ztail::d, fitdone_, MuonSubdetId::GEM, gemetapartition(), runTauDisplay::gp, hits_, CSCChamber::id(), intercept_, CSCChamber::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.

                                                             {
     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;
     }
   }

References Exception, and gemetapartmap_.

Referenced by fit2(), and fitlsq().

◆ hits()

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

inline

Definition at line 131 of file GEMCSCSegFit.h.

131 { return hits_; }

References hits_.

◆ intercept()

LocalPoint GEMCSCSegFit::intercept ( ) const

inline

Definition at line 136 of file GEMCSCSegFit.h.

136 { return intercept_; }

References intercept_.

Referenced by GEMCSCSegAlgoRR::buildSegments().

◆ localdir()

LocalVector GEMCSCSegFit::localdir ( ) const

inline

Definition at line 137 of file GEMCSCSegFit.h.

137 { return localdir_; }

References localdir_.

Referenced by GEMCSCSegAlgoRR::buildSegments().

◆ ndof()

int GEMCSCSegFit::ndof ( void ) const

inline

Definition at line 135 of file GEMCSCSegFit.h.

135 { return ndof_; }

References ndof_.

◆ nhits()

size_t GEMCSCSegFit::nhits ( void ) const

inline

Definition at line 133 of file GEMCSCSegFit.h.

133 { return hits_.size(); }

References hits_.

Referenced by fit().

◆ 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.

                                           {
     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();
     }
   }

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

Referenced by fit2(), and fitlsq().

◆ scaleXError()

double GEMCSCSegFit::scaleXError ( void ) const

inline

Definition at line 132 of file GEMCSCSegFit.h.

132 { return scaleXError_; }

References 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.

119 { scaleXError_ = factor; }

References DQMScaleToClient_cfi::factor, and scaleXError_.

◆ 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().

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Member Typedef Documentation

◆ SMatrix16by4

◆ SMatrix4

◆ SMatrixSym16

◆ SMatrixSym2

◆ SMatrixSym4

◆ SVector4

Constructor & Destructor Documentation

◆ GEMCSCSegFit()

◆ ~GEMCSCSegFit()

Member Function Documentation

◆ chi2()

◆ covarianceMatrix()

◆ cscchamber()

◆ csclayer()

◆ derivativeMatrix()

◆ fit()

◆ fit2()

◆ fitdone()

◆ fitlsq()

◆ flipErrors()

◆ gemetapartition()

◆ hits()

◆ intercept()

◆ localdir()

◆ ndof()

◆ nhits()

◆ Rdev()

◆ refcscchamber()

◆ scaleXError()

◆ setChi2()

◆ setOutFromIP()

◆ setScaleXError()

◆ weightMatrix()

◆ xdev()

◆ xfit()

◆ ydev()

◆ yfit()

Member Data Documentation

◆ chi2_

◆ csclayermap_

◆ fitdone_

◆ gemetapartmap_

◆ hits_

◆ intercept_

◆ localdir_

◆ ndof_

◆ refcscchamber_

◆ refid_

◆ scaleXError_

◆ uslope_

◆ vslope_