CMSSW Reference Manual

TFitParticleEMomDev::TFitParticleEMomDev ( )

Referenced by clone().

TFitParticleEMomDev::TFitParticleEMomDev ( const TFitParticleEMomDev & fitParticle )

Definition at line 42 of file TFitParticleEMomDev.cc.

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

00043   :TAbsFitParticle( fitParticle.GetName(), fitParticle.GetTitle() )
00044 {
00045 
00046   _nPar = fitParticle._nPar;
00047   _u1 = fitParticle._u1;
00048   _u2 = fitParticle._u2;
00049   _u3 = fitParticle._u3;
00050   _covMatrix.ResizeTo(  fitParticle._covMatrix );
00051   _covMatrix = fitParticle._covMatrix;
00052   _iniparameters.ResizeTo( fitParticle._iniparameters );
00053   _iniparameters = fitParticle._iniparameters;
00054   _parameters.ResizeTo( fitParticle._parameters );
00055   _parameters = fitParticle._parameters;
00056   _pini = fitParticle._pini;
00057   _pcurr = fitParticle._pcurr;
00058 
00059 }

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

Definition at line 61 of file TFitParticleEMomDev.cc.

References init().

00062   :TAbsFitParticle()  
00063 {
00064   init(pini, theCovMatrix);
00065 }

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

Definition at line 67 of file TFitParticleEMomDev.cc.

References init().

00069   :TAbsFitParticle(name, title)  
00070 {
00071   init(pini, theCovMatrix);
00072 }

TFitParticleEMomDev::~TFitParticleEMomDev ( ) [virtual]

Definition at line 86 of file TFitParticleEMomDev.cc.

00086                                           {
00087 
00088 }

TLorentzVector * TFitParticleEMomDev::calc4Vec ( const TMatrixD * params ) [virtual]

Implements TAbsFitParticle.

Definition at line 145 of file TFitParticleEMomDev.cc.

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

00145                                                                       {
00146   // Calculates a 4vector corresponding to the given
00147   // parameter values
00148 
00149   if (params == 0) {
00150     return 0;
00151   }
00152 
00153   if ( params->GetNcols() != 1 || params->GetNrows() !=_nPar ) {
00154     edm::LogError ("WrongMatrixSize")
00155       << GetName() << "::calc4Vec - Parameter matrix has wrong size.";
00156     return 0;
00157   }
00158 
00159   Double_t X = _pini.P() * (*params)(0,0) *_u1.X() +
00160     (*params)(1,0) * _u2.X()+
00161     (*params)(2,0) * _u3.X();
00162   Double_t Y = _pini.P() * (*params)(0,0) *_u1.Y() +
00163     (*params)(1,0) * _u2.Y()+
00164     (*params)(2,0) * _u3.Y();
00165   Double_t Z = _pini.P() * (*params)(0,0) *_u1.Z() +
00166     (*params)(1,0) * _u2.Z()+
00167     (*params)(2,0) * _u3.Z();
00168   Double_t E = _pini.E()*(*params)(3,0);
00169 
00170   TLorentzVector* vec = new TLorentzVector( X, Y, Z, E );
00171   return vec;
00172 
00173 }

TAbsFitParticle * TFitParticleEMomDev::clone ( TString newname = "" ) const [virtual]

Implements TAbsFitParticle.

Definition at line 74 of file TFitParticleEMomDev.cc.

References TFitParticleEMomDev().

00074                                                                    {
00075   // Returns a copy of itself
00076   
00077   TAbsFitParticle* myclone = new TFitParticleEMomDev( *this );
00078   if ( newname.Length() > 0 ) myclone->SetName(newname);
00079   return myclone;
00080 
00081 }

TMatrixD * TFitParticleEMomDev::getDerivative ( ) [virtual]

Implements TAbsFitParticle.

Definition at line 175 of file TFitParticleEMomDev.cc.

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

00175                                              {
00176   // returns derivative dP/dy with P=(p,E) and y=(r, theta, phi, d) 
00177   // the free parameters of the fit. The columns of the matrix contain 
00178   // (dP/dr, dP/dtheta, dP/dphi, dP/dd).
00179 
00180   TMatrixD* DerivativeMatrix = new TMatrixD(4,4);
00181   (*DerivativeMatrix) *= 0.;
00182 
00183   //1st column: dP/dr
00184   (*DerivativeMatrix)(0,0)=_pini.P()*_u1.X();
00185   (*DerivativeMatrix)(1,0)=_pini.P()*_u1.Y();
00186   (*DerivativeMatrix)(2,0)=_pini.P()*_u1.Z();
00187   (*DerivativeMatrix)(3,0)=0.;
00188 
00189   //2nd column: dP/dtheta
00190   (*DerivativeMatrix)(0,1)=_u2.X();
00191   (*DerivativeMatrix)(1,1)=_u2.Y();
00192   (*DerivativeMatrix)(2,1)=_u2.Z();
00193   (*DerivativeMatrix)(3,1)=0.;
00194 
00195    //3rd column: dP/dphi
00196   (*DerivativeMatrix)(0,2)=_u3.X();
00197   (*DerivativeMatrix)(1,2)=_u3.Y();
00198   (*DerivativeMatrix)(2,2)=_u3.Z();
00199   (*DerivativeMatrix)(3,2)=0.;
00200 
00201    //4th column: dP/dm
00202   (*DerivativeMatrix)(0,3)=0.;
00203   (*DerivativeMatrix)(1,3)=0.;
00204   (*DerivativeMatrix)(2,3)=0.;
00205   (*DerivativeMatrix)(3,3)=_pini.E();
00206 
00207   return DerivativeMatrix;  
00208 }

void TFitParticleEMomDev::init	(	TLorentzVector *	pini,
		const TMatrixD *	theCovMatrix
	)			`[protected]`

Definition at line 93 of file TFitParticleEMomDev.cc.

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

Referenced by TFitParticleEMomDev().

00093                                                                                   {
00094 
00095   _nPar = 4;
00096   setIni4Vec(pini);
00097   setCovMatrix(theCovMatrix);
00098 
00099 }

void TFitParticleEMomDev::setIni4Vec ( const TLorentzVector * pini ) [virtual]

Implements TAbsFitParticle.

Definition at line 101 of file TFitParticleEMomDev.cc.

References TAbsFitParticle::_iniparameters, TAbsFitParticle::_nPar, TAbsFitParticle::_parameters, TAbsFitParticle::_pcurr, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, TAbsFitParticle::_u3, and muonGeometry::TVector3().

Referenced by init().

00101                                                                {
00102   // Set the initial 4vector. Will also set the 
00103   // inital parameter values
00104   _iniparameters.ResizeTo(_nPar,1);
00105   _parameters.ResizeTo(_nPar,1);
00106 
00107   if (pini == 0) {
00108     _iniparameters(0,0)=0.;
00109     _iniparameters(1,0)=0.;
00110     _iniparameters(2,0)=0.;
00111     _iniparameters(3,0)=0.;
00112 
00113     _u1.SetXYZ(0., 0., 0.);
00114     _u3.SetXYZ(0., 0., 0.);
00115     _u2.SetXYZ(0., 0., 0.);
00116     _pini.SetXYZT(0., 0., 0., 0.);
00117     _pcurr = _pini;
00118     
00119   } else {
00120     _iniparameters(0,0)=1.;
00121     _iniparameters(1,0)=0.;
00122     _iniparameters(2,0)=0.;
00123     _iniparameters(3,0)=1.;
00124 
00125     _pini = (*pini);
00126     _pcurr = _pini;
00127 
00128     _u1 = pini->Vect();
00129     _u1 *= 1./_u1.Mag();
00130     
00131     TVector3 uz(0., 0., 1.);
00132     _u3 = uz.Cross(_u1);
00133     _u3 *= 1./_u3.Mag();
00134     
00135     _u2 = _u3.Cross(_u1);
00136     _u2 *= 1./_u2.Mag();
00137 
00138   }
00139 
00140   // reset parameters
00141   _parameters = _iniparameters;
00142 
00143 }

TMatrixD * TFitParticleEMomDev::transform ( const TLorentzVector & vec ) [virtual]

Implements TAbsFitParticle.

Definition at line 210 of file TFitParticleEMomDev.cc.

References TAbsFitParticle::_nPar, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, TAbsFitParticle::_u3, rot, and muonGeometry::TVector3().

00210                                                                   {
00211   // Returns the parameters corresponding to the given 
00212   // 4vector wrt. to the current base vectors u_r, u_theta, and u_phi
00213 
00214   // construct rotation matrix
00215   TRotation rot;
00216   rot.RotateAxes( _u1, _u2, _u3 );
00217   rot.Invert();
00218 
00219   // rotate vector
00220   TVector3 vec3( vec.Vect() );
00221   vec3.Transform( rot );
00222 
00223   // retrieve parameters
00224   TMatrixD* tparams = new TMatrixD( _nPar, 1 );
00225   (*tparams)(0,0) = vec3.x()/_pini.P();
00226   (*tparams)(1,0) = vec3.y();
00227   (*tparams)(2,0) = vec3.z();
00228   (*tparams)(3,0) = vec.E()/_pini.E();
00229 
00230   return tparams;
00231 
00232 }


Public Member Functions
virtual TLorentzVector *	calc4Vec (const TMatrixD *params)
virtual TAbsFitParticle *	clone (TString newname="") const
virtual TMatrixD *	getDerivative ()
virtual void	setIni4Vec (const TLorentzVector *pini)
	TFitParticleEMomDev (const TString &name, const TString &title, TLorentzVector pini, const TMatrixD theCovMatrix)
	TFitParticleEMomDev (TLorentzVector pini, const TMatrixD theCovMatrix)
	TFitParticleEMomDev (const TFitParticleEMomDev &fitParticle)
	TFitParticleEMomDev ()
virtual TMatrixD *	transform (const TLorentzVector &vec)
virtual	~TFitParticleEMomDev ()
Protected Member Functions
void	init (TLorentzVector pini, const TMatrixD theCovMatrix)

TFitParticleEMomDev Class Reference

Public Member Functions

Protected Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation