TFitParticleEMomDev Class Reference

#include <TFitParticleEMomDev.h>

Inheritance diagram for TFitParticleEMomDev:

Public Member Functions
TLorentzVector *	calc4Vec (const TMatrixD *params) override
TAbsFitParticle *	clone (const TString &newname="") const override
TMatrixD *	getDerivative () override
void	setIni4Vec (const TLorentzVector *pini) override
	TFitParticleEMomDev ()
	TFitParticleEMomDev (const TFitParticleEMomDev &fitParticle)
	TFitParticleEMomDev (const TString &name, const TString &title, TLorentzVector *pini, const TMatrixD *theCovMatrix)
	TFitParticleEMomDev (TLorentzVector *pini, const TMatrixD *theCovMatrix)
TMatrixD *	transform (const TLorentzVector &vec) override
	~TFitParticleEMomDev () override
Public Member Functions inherited from TAbsFitParticle
virtual void	applycorr (TMatrixD *corrMatrix)
virtual const TMatrixD *	getCovMatrix () const
virtual const TMatrixD *	getCovMatrixDeltaY ()
virtual const TMatrixD *	getCovMatrixFit () const
const TLorentzVector *	getCurr4Vec ()
TString	getInfoString ()
const TLorentzVector *	getIni4Vec ()
Int_t	getNPar () const
const TMatrixD *	getParCurr ()
const TMatrixD *	getParIni ()
virtual const TMatrixD *	getPull ()
const TVector3 *	getu1 ()
const TVector3 *	getu2 ()
const TVector3 *	getu3 ()
virtual void	print ()
virtual void	reset ()
virtual void	setCovMatrix (const TMatrixD *theCovMatrix)
virtual void	setCovMatrixFit (const TMatrixD *theCovMatrixFit)
void	setParIni (const TMatrixD *parini)
	TAbsFitParticle ()
	TAbsFitParticle (const TString &name, const TString &title)
	~TAbsFitParticle () override

Protected Member Functions
void	init (TLorentzVector *pini, const TMatrixD *theCovMatrix)
Protected Member Functions inherited from TAbsFitParticle
void	calcCovMatrixDeltaY ()

Additional Inherited Members
Protected Attributes inherited from TAbsFitParticle
TMatrixD	_covMatrix
TMatrixD	_covMatrixDeltaY
TMatrixD	_covMatrixFit
TMatrixD	_iniparameters
Int_t	_nPar
TMatrixD	_parameters
TLorentzVector	_pcurr
TLorentzVector	_pini
TMatrixD	_pull
TVector3	_u1
TVector3	_u2
TVector3	_u3

Detailed Description

Definition at line 10 of file TFitParticleEMomDev.h.

Constructor & Destructor Documentation

TFitParticleEMomDev() [1/4]

TFitParticleEMomDev::TFitParticleEMomDev

(

)

Definition at line 31 of file TFitParticleEMomDev.cc.

: TAbsFitParticle() { init(nullptr, nullptr); }

References init().

Referenced by clone().

TFitParticleEMomDev() [2/4]

TFitParticleEMomDev::TFitParticleEMomDev

(

const TFitParticleEMomDev &

fitParticle

)

Definition at line 33 of file TFitParticleEMomDev.cc.

    : TAbsFitParticle(fitParticle.GetName(), fitParticle.GetTitle()) {
  _nPar = fitParticle._nPar;
  _u1 = fitParticle._u1;
  _u2 = fitParticle._u2;
  _u3 = fitParticle._u3;
  _covMatrix.ResizeTo(fitParticle._covMatrix);
  _covMatrix = fitParticle._covMatrix;
  _iniparameters.ResizeTo(fitParticle._iniparameters);
  _iniparameters = fitParticle._iniparameters;
  _parameters.ResizeTo(fitParticle._parameters);
  _parameters = fitParticle._parameters;
  _pini = fitParticle._pini;
  _pcurr = fitParticle._pcurr;
}

References TAbsFitParticle::_covMatrix, TAbsFitParticle::_iniparameters, TAbsFitParticle::_nPar, TAbsFitParticle::_parameters, TAbsFitParticle::_pcurr, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, and TAbsFitParticle::_u3.

TFitParticleEMomDev() [3/4]

TFitParticleEMomDev::TFitParticleEMomDev	(	TLorentzVector *	pini,
		const TMatrixD *	theCovMatrix
	)

Definition at line 49 of file TFitParticleEMomDev.cc.

                                                                                           : TAbsFitParticle() {
  init(pini, theCovMatrix);
}

References init().

TFitParticleEMomDev() [4/4]

TFitParticleEMomDev::TFitParticleEMomDev	(	const TString &	name,
		const TString &	title,
		TLorentzVector *	pini,
		const TMatrixD *	theCovMatrix
	)

Definition at line 53 of file TFitParticleEMomDev.cc.

    : TAbsFitParticle(name, title) {
  init(pini, theCovMatrix);
}

References init().

~TFitParticleEMomDev()

TFitParticleEMomDev::~TFitParticleEMomDev ( )

override

Definition at line 73 of file TFitParticleEMomDev.cc.

{}

Member Function Documentation

calc4Vec()

TLorentzVector * TFitParticleEMomDev::calc4Vec ( const TMatrixD *  params )

overridevirtual

Implements TAbsFitParticle.

Definition at line 126 of file TFitParticleEMomDev.cc.

                                                                    {
  // Calculates a 4vector corresponding to the given
  // parameter values
 
  if (params == nullptr) {
    return nullptr;
  }
 
  if (params->GetNcols() != 1 || params->GetNrows() != _nPar) {
    edm::LogError("WrongMatrixSize") << GetName() << "::calc4Vec - Parameter matrix has wrong size.";
    return nullptr;
  }
 
  Double_t X = _pini.P() * (*params)(0, 0) * _u1.X() + (*params)(1, 0) * _u2.X() + (*params)(2, 0) * _u3.X();
  Double_t Y = _pini.P() * (*params)(0, 0) * _u1.Y() + (*params)(1, 0) * _u2.Y() + (*params)(2, 0) * _u3.Y();
  Double_t Z = _pini.P() * (*params)(0, 0) * _u1.Z() + (*params)(1, 0) * _u2.Z() + (*params)(2, 0) * _u3.Z();
  Double_t E = _pini.E() * (*params)(3, 0);
 
  TLorentzVector* vec = new TLorentzVector(X, Y, Z, E);
  return vec;
}

References TAbsFitParticle::_nPar, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, TAbsFitParticle::_u3, CalibrationSummaryClient_cfi::params, X, BeamSpotPI::Y, and BeamSpotPI::Z.

clone()

TAbsFitParticle * TFitParticleEMomDev::clone ( const TString &  newname = "" ) const

overridevirtual

Implements TAbsFitParticle.

Definition at line 61 of file TFitParticleEMomDev.cc.

                                                                        {
  // Returns a copy of itself
 
  TAbsFitParticle* myclone = new TFitParticleEMomDev(*this);
  if (newname.Length() > 0)
    myclone->SetName(newname);
  return myclone;
}

References TFitParticleEMomDev().

getDerivative()

TMatrixD * TFitParticleEMomDev::getDerivative ( )

overridevirtual

Implements TAbsFitParticle.

Definition at line 148 of file TFitParticleEMomDev.cc.

                                             {
  // returns derivative dP/dy with P=(p,E) and y=(r, theta, phi, d)
  // the free parameters of the fit. The columns of the matrix contain
  // (dP/dr, dP/dtheta, dP/dphi, dP/dd).
 
  TMatrixD* DerivativeMatrix = new TMatrixD(4, 4);
  (*DerivativeMatrix) *= 0.;
 
  //1st column: dP/dr
  (*DerivativeMatrix)(0, 0) = _pini.P() * _u1.X();
  (*DerivativeMatrix)(1, 0) = _pini.P() * _u1.Y();
  (*DerivativeMatrix)(2, 0) = _pini.P() * _u1.Z();
  (*DerivativeMatrix)(3, 0) = 0.;
 
  //2nd column: dP/dtheta
  (*DerivativeMatrix)(0, 1) = _u2.X();
  (*DerivativeMatrix)(1, 1) = _u2.Y();
  (*DerivativeMatrix)(2, 1) = _u2.Z();
  (*DerivativeMatrix)(3, 1) = 0.;
 
  //3rd column: dP/dphi
  (*DerivativeMatrix)(0, 2) = _u3.X();
  (*DerivativeMatrix)(1, 2) = _u3.Y();
  (*DerivativeMatrix)(2, 2) = _u3.Z();
  (*DerivativeMatrix)(3, 2) = 0.;
 
  //4th column: dP/dm
  (*DerivativeMatrix)(0, 3) = 0.;
  (*DerivativeMatrix)(1, 3) = 0.;
  (*DerivativeMatrix)(2, 3) = 0.;
  (*DerivativeMatrix)(3, 3) = _pini.E();
 
  return DerivativeMatrix;
}

References TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, and TAbsFitParticle::_u3.

init()

void TFitParticleEMomDev::init ( TLorentzVector *  pini, const TMatrixD *  theCovMatrix  )

protected

Definition at line 78 of file TFitParticleEMomDev.cc.

                                                                                 {
  _nPar = 4;
  setIni4Vec(pini);
  setCovMatrix(theCovMatrix);
}

References TAbsFitParticle::_nPar, TAbsFitParticle::setCovMatrix(), and setIni4Vec().

Referenced by TFitParticleEMomDev().

setIni4Vec()

void TFitParticleEMomDev::setIni4Vec ( const TLorentzVector *  pini )

overridevirtual

Implements TAbsFitParticle.

Definition at line 84 of file TFitParticleEMomDev.cc.

                                                               {
  // Set the initial 4vector. Will also set the
  // inital parameter values
  _iniparameters.ResizeTo(_nPar, 1);
  _parameters.ResizeTo(_nPar, 1);
 
  if (pini == nullptr) {
    _iniparameters(0, 0) = 0.;
    _iniparameters(1, 0) = 0.;
    _iniparameters(2, 0) = 0.;
    _iniparameters(3, 0) = 0.;
 
    _u1.SetXYZ(0., 0., 0.);
    _u3.SetXYZ(0., 0., 0.);
    _u2.SetXYZ(0., 0., 0.);
    _pini.SetXYZT(0., 0., 0., 0.);
    _pcurr = _pini;
 
  } else {
    _iniparameters(0, 0) = 1.;
    _iniparameters(1, 0) = 0.;
    _iniparameters(2, 0) = 0.;
    _iniparameters(3, 0) = 1.;
 
    _pini = (*pini);
    _pcurr = _pini;
 
    _u1 = pini->Vect();
    _u1 *= 1. / _u1.Mag();
 
    TVector3 uz(0., 0., 1.);
    _u3 = uz.Cross(_u1);
    _u3 *= 1. / _u3.Mag();
 
    _u2 = _u3.Cross(_u1);
    _u2 *= 1. / _u2.Mag();
  }
 
  // reset parameters
  _parameters = _iniparameters;
}

References TAbsFitParticle::_iniparameters, TAbsFitParticle::_nPar, TAbsFitParticle::_parameters, TAbsFitParticle::_pcurr, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, and TAbsFitParticle::_u3.

Referenced by init().

transform()

TMatrixD * TFitParticleEMomDev::transform ( const TLorentzVector &  vec )

overridevirtual

Implements TAbsFitParticle.

Definition at line 183 of file TFitParticleEMomDev.cc.

                                                                  {
  // Returns the parameters corresponding to the given
  // 4vector wrt. to the current base vectors u_r, u_theta, and u_phi
 
  // construct rotation matrix
  TRotation rot;
  rot.RotateAxes(_u1, _u2, _u3);
  rot.Invert();
 
  // rotate vector
  TVector3 vec3(vec.Vect());
  vec3.Transform(rot);
 
  // retrieve parameters
  TMatrixD* tparams = new TMatrixD(_nPar, 1);
  (*tparams)(0, 0) = vec3.x() / _pini.P();
  (*tparams)(1, 0) = vec3.y();
  (*tparams)(2, 0) = vec3.z();
  (*tparams)(3, 0) = vec.E() / _pini.E();
 
  return tparams;
}

References TAbsFitParticle::_nPar, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, TAbsFitParticle::_u3, and makeMuonMisalignmentScenario::rot.