#include <TFitParticleMomDev.h>

Inheritance diagram for TFitParticleMomDev:

Public Member Functions
TLorentzVector *	calc4Vec (const TMatrixD *params) override

TAbsFitParticle *	clone (const TString &newname=TString("")) const override

TMatrixD *	getDerivative () override

void	setIni4Vec (const TLorentzVector *pini) override

	TFitParticleMomDev ()

	TFitParticleMomDev (const TFitParticleMomDev &fitParticle)

	TFitParticleMomDev (TLorentzVector pini, const TMatrixD theCovMatrix)

	TFitParticleMomDev (const TString &name, const TString &title, TLorentzVector pini, const TMatrixD theCovMatrix)

TMatrixD *	transform (const TLorentzVector &vec) override

	~TFitParticleMomDev () 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 TFitParticleMomDev.h.

Constructor & Destructor Documentation

TFitParticleMomDev::TFitParticleMomDev ( )

Definition at line 32 of file TFitParticleMomDev.cc.

References init().

Referenced by clone().

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

TAbsFitParticle::TAbsFitParticle

TAbsFitParticle()

Definition: TAbsFitParticle.cc:18

TFitParticleMomDev::init

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

Definition: TFitParticleMomDev.cc:79

TFitParticleMomDev::TFitParticleMomDev ( const TFitParticleMomDev & fitParticle )

Definition at line 34 of file TFitParticleMomDev.cc.

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

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

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

Definition at line 50 of file TFitParticleMomDev.cc.

References init().

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

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

Definition at line 54 of file TFitParticleMomDev.cc.

References init().

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

TFitParticleMomDev::~TFitParticleMomDev ( )

override

Definition at line 74 of file TFitParticleMomDev.cc.

74 {}

Member Function Documentation

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

overridevirtual

Implements TAbsFitParticle.

Definition at line 122 of file TFitParticleMomDev.cc.

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

                                                                    {
   // 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 = TMath::Sqrt(X * X + Y * Y + Z * Z + (*params)(3, 0) * (*params)(3, 0) * _pini.M2());
 
   TLorentzVector* vec = new TLorentzVector(X, Y, Z, E);
   return vec;
 }

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

overridevirtual

Implements TAbsFitParticle.

Definition at line 62 of file TFitParticleMomDev.cc.

References TFitParticleMomDev().

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

TMatrixD * TFitParticleMomDev::getDerivative ( )

overridevirtual

Implements TAbsFitParticle.

Definition at line 144 of file TFitParticleMomDev.cc.

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

                                             {
   // 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.M() * _pini.M() * _parameters(3, 0) / _pcurr.E();
 
   return DerivativeMatrix;
 }

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

protected

Definition at line 79 of file TFitParticleMomDev.cc.

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

Referenced by TFitParticleMomDev().

                                                                                 {
   _nPar = 4;
   _iniparameters.ResizeTo(_nPar, 1);
   _iniparameters(0, 0) = 1.;
   _iniparameters(1, 0) = 0.;
   _iniparameters(2, 0) = 0.;
   _iniparameters(3, 0) = 1.;
   _parameters.ResizeTo(_nPar, 1);
   _parameters = _iniparameters;
   setIni4Vec(pini);
   setCovMatrix(theCovMatrix);
 }

void TFitParticleMomDev::setIni4Vec ( const TLorentzVector * pini )

overridevirtual

Implements TAbsFitParticle.

Definition at line 92 of file TFitParticleMomDev.cc.

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

Referenced by init().

                                                               {
   // Set the initial 4vector. Will also set the
   // inital parameter values
 
   if (pini == nullptr) {
     _u1.SetXYZ(0., 0., 0.);
     _u3.SetXYZ(0., 0., 0.);
     _u2.SetXYZ(0., 0., 0.);
     _pini.SetXYZT(0., 0., 0., 0.);
     _pcurr = _pini;
 
   } else {
     _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;
 }

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

overridevirtual

Implements TAbsFitParticle.

Definition at line 179 of file TFitParticleMomDev.cc.

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

                                                                  {
   // 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(0) / _pini.P();
   (*tparams)(1, 0) = vec3(1);
   (*tparams)(2, 0) = vec3(2);
   (*tparams)(3, 0) = vec.M() / _pini.M();
 
   return tparams;
 }

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation