CMSSW Reference Manual

TFitParticleMCMomDev::TFitParticleMCMomDev ( )

Referenced by clone().

TFitParticleMCMomDev::TFitParticleMCMomDev ( const TFitParticleMCMomDev & fitParticle )

Definition at line 43 of file TFitParticleMCMomDev.cc.

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

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

TFitParticleMCMomDev::TFitParticleMCMomDev	(	TVector3 *	p,
		Double_t	M,
		const TMatrixD *	theCovMatrix
	)

Definition at line 62 of file TFitParticleMCMomDev.cc.

References init().

00063   :TAbsFitParticle()
00064 { 
00065   init(p, M, theCovMatrix);
00066 }

TFitParticleMCMomDev::TFitParticleMCMomDev	(	const TString &	name,
		const TString &	title,
		TVector3 *	p,
		Double_t	M,
		const TMatrixD *	theCovMatrix
	)

Definition at line 68 of file TFitParticleMCMomDev.cc.

References init().

00070   :TAbsFitParticle(name, title)
00071 { 
00072   init(p, M, theCovMatrix);
00073 }

TFitParticleMCMomDev::~TFitParticleMCMomDev ( ) [virtual]

Definition at line 87 of file TFitParticleMCMomDev.cc.

00087                                             {
00088 
00089 }

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

Implements TAbsFitParticle.

Definition at line 109 of file TFitParticleMCMomDev.cc.

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

00109                                                                        {
00110   // Calculates a 4vector corresponding to the given
00111   // parameter values
00112 
00113   if (params == 0) {
00114     return 0;
00115   }
00116 
00117   if ( params->GetNcols() != 1 || params->GetNrows() !=_nPar ) {
00118     edm::LogError ("WrongMatrixSize")
00119       << GetName() << "::calc4Vec - Parameter matrix has wrong size.";
00120     return 0;
00121   }
00122 
00123   Double_t X =  _pini.P() * (*params)(0,0) *_u1.X()+
00124     (*params)(1,0) * _u2.X()+
00125     (*params)(2,0) * _u3.X() ;
00126   Double_t Y = _pini.P() * (*params)(0,0) *_u1.Y()+
00127     (*params)(1,0) * _u2.Y()+
00128     (*params)(2,0) * _u3.Y() ;
00129   Double_t Z = _pini.P() * (*params)(0,0) *_u1.Z()+
00130     (*params)(1,0) * _u2.Z()+
00131     (*params)(2,0) * _u3.Z() ;
00132   Double_t E =  TMath::Sqrt(  X*X + Y*Y + Z*Z + _pini.M2() );
00133 
00134   TLorentzVector* vec = new TLorentzVector( X, Y, Z, E );
00135   return vec;
00136 
00137 }

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

Implements TAbsFitParticle.

Definition at line 75 of file TFitParticleMCMomDev.cc.

References TFitParticleMCMomDev().

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

TMatrixD * TFitParticleMCMomDev::getDerivative ( ) [virtual]

Implements TAbsFitParticle.

Definition at line 178 of file TFitParticleMCMomDev.cc.

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

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

void TFitParticleMCMomDev::init	(	TVector3 *	p,
		Double_t	M,
		const TMatrixD *	theCovMatrix
	)			`[protected]`

Definition at line 95 of file TFitParticleMCMomDev.cc.

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

Referenced by TFitParticleMCMomDev().

00095                                                                                      {
00096 
00097   _nPar = 3;
00098   setIni4Vec(p, M);
00099   _iniparameters.ResizeTo(_nPar,1);
00100   _iniparameters(0,0)=1.;
00101   _iniparameters(1,0)=0.;
00102   _iniparameters(2,0)=0.;
00103   _parameters.ResizeTo(_nPar, 1);
00104   _parameters = _iniparameters;
00105   setCovMatrix(theCovMatrix);
00106   
00107 }

void TFitParticleMCMomDev::setIni4Vec	(	const TVector3 *	p,
		Double_t	M
	)

Definition at line 148 of file TFitParticleMCMomDev.cc.

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

00148                                                                    {
00149   // Set the initial 4vector. Will also set the 
00150   // inital parameter values
00151 
00152   if ( p == 0 ) {
00153 
00154     _pini.SetXYZM( 0., 0., 0., M);
00155     _pcurr = _pini;
00156 
00157   } else {
00158 
00159     _pini.SetXYZM( p->x(), p->y(), p->z(), M);
00160     _pcurr = _pini;
00161 
00162     _u1 = (*p);
00163     _u1 *= 1./_u1.Mag();
00164 
00165     TVector3 uz(0., 0., 1.);
00166     _u3 = uz.Cross(_u1);
00167     _u3 *= 1./_u3.Mag();
00168   
00169     _u2 = _u3.Cross(_u1);
00170     _u2 *= 1./_u2.Mag();
00171   
00172   }
00173 
00174   _parameters = _iniparameters;
00175 
00176 }

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

Implements TAbsFitParticle.

Definition at line 139 of file TFitParticleMCMomDev.cc.

References muonGeometry::TVector3().

Referenced by init().

00139                                                                 {
00140   // Set the initial 4vector. Will also set the 
00141   // inital parameter values
00142 
00143   TVector3 vec( pini->Vect() );
00144   setIni4Vec( &vec, pini->M() );
00145 
00146 }

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

Implements TAbsFitParticle.

Definition at line 213 of file TFitParticleMCMomDev.cc.

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

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


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

TFitParticleMCMomDev Class Reference

Public Member Functions

Protected Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation