CMS 3D CMS Logo

TFitParticleESpher Class Reference

#include <PhysicsTools/KinFitter/interface/TFitParticleESpher.h>

Inheritance diagram for TFitParticleESpher:

TAbsFitParticle

List of all members.

Public Member Functions

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

Protected Member Functions

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

Detailed Description

Definition at line 10 of file TFitParticleESpher.h.

Constructor & Destructor Documentation

TFitParticleESpher::TFitParticleESpher (  ) 

Referenced by clone().

TFitParticleESpher::TFitParticleESpher ( const TFitParticleESpher &  fitParticle  ) 

Definition at line 36 of file TFitParticleESpher.cc.

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

00037   :TAbsFitParticle( fitParticle.GetName(), fitParticle.GetTitle() )
00038 {
00039 
00040   _nPar = fitParticle._nPar;
00041   _u1 = fitParticle._u1;
00042   _u2 = fitParticle._u2;
00043   _u3 = fitParticle._u3;
00044   _covMatrix.ResizeTo(  fitParticle._covMatrix );
00045   _covMatrix = fitParticle._covMatrix;
00046   _iniparameters.ResizeTo( fitParticle._iniparameters );
00047   _iniparameters = fitParticle._iniparameters;
00048   _parameters.ResizeTo( fitParticle._parameters );
00049   _parameters = fitParticle._parameters;
00050   _pini = fitParticle._pini;
00051   _pcurr = fitParticle._pcurr;
00052 
00053 }

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

Definition at line 55 of file TFitParticleESpher.cc.

References init().

00056   :TAbsFitParticle()  
00057 {
00058   init(pini, theCovMatrix);
00059 }

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

Definition at line 61 of file TFitParticleESpher.cc.

References init().

00063   :TAbsFitParticle(name, title)  
00064 {
00065   init(pini, theCovMatrix);
00066 }

TFitParticleESpher::~TFitParticleESpher (  )  [virtual]

Definition at line 80 of file TFitParticleESpher.cc.

00080                                         {
00081 
00082 }

Member Function Documentation

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

Implements TAbsFitParticle.

Definition at line 95 of file TFitParticleESpher.cc.

References TAbsFitParticle::_nPar, TAbsFitParticle::_pini, d, phi, r, theta, and X.

00095                                                                      {
00096   // Calculates a 4vector corresponding to the given
00097   // parameter values
00098 
00099   if (params == 0) {
00100     return 0;
00101   }
00102 
00103   if ( params->GetNcols() != 1 || params->GetNrows() !=_nPar ) {
00104     edm::LogError ("WrongMatrixSize")
00105       << GetName() << "::calc4Vec - Parameter matrix has wrong size.";
00106     return 0;
00107   }
00108 
00109   Double_t r = (*params)(0,0);
00110   Double_t theta = (*params)(1,0);
00111   Double_t phi = (*params)(2,0);
00112   Double_t d = (*params)(3,0);
00113 
00114   Double_t X = r*TMath::Cos(phi)*TMath::Sin(theta);
00115   Double_t Y = r*TMath::Sin(phi)*TMath::Sin(theta);
00116   Double_t Z = r*TMath::Cos(theta);
00117   Double_t E = d * _pini.E();
00118                 
00119   TLorentzVector* vec = new TLorentzVector( X, Y, Z, E );
00120   return vec;
00121 
00122 }

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

Implements TAbsFitParticle.

Definition at line 68 of file TFitParticleESpher.cc.

References TFitParticleESpher().

00068                                                                   {
00069   // Returns a copy of itself
00070   
00071   TAbsFitParticle* myclone = new TFitParticleESpher( *this );
00072   if ( newname.Length() > 0 ) myclone->SetName(newname);
00073   return myclone;
00074 
00075 }

TMatrixD * TFitParticleESpher::getDerivative (  )  [virtual]

Implements TAbsFitParticle.

Definition at line 174 of file TFitParticleESpher.cc.

References TAbsFitParticle::_parameters, TAbsFitParticle::_pini, phi, r, and theta.

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

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

Definition at line 87 of file TFitParticleESpher.cc.

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

Referenced by TFitParticleESpher().

00087                                                                                  {
00088 
00089   _nPar = 4;
00090   setIni4Vec(pini);
00091   setCovMatrix(theCovMatrix);
00092 
00093 }

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

Implements TAbsFitParticle.

Definition at line 124 of file TFitParticleESpher.cc.

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

Referenced by init().

00124                                                               {
00125   // Set the initial 4vector. Will also set the 
00126   // inital parameter values 
00127 
00128   if (pini == 0) {
00129 
00130     _u1.SetXYZ(0., 0., 0.);
00131     _u3.SetXYZ(0., 0., 0.);
00132     _u2.SetXYZ(0., 0., 0.);
00133     _pini.SetXYZT(0., 0., 0., 0.);
00134     _pcurr = _pini;
00135 
00136     _iniparameters.ResizeTo(_nPar,1);
00137     _iniparameters(0,0) = 0.;
00138     _iniparameters(1,0) = 0.;
00139     _iniparameters(2,0) = 0.;
00140     _iniparameters(3,0) = 1.;
00141     
00142     _parameters.ResizeTo(_nPar,1);
00143     _parameters(0,0) = 0.;
00144     _parameters(1,0) = 0.;
00145     _parameters(2,0) = 0.;   
00146     _parameters(3,0) = 1.;   
00147     
00148   } else {
00149     
00150     Double_t r = pini->P();
00151     Double_t theta = pini->Theta();
00152     Double_t phi = pini->Phi();
00153     
00154     _pini = (*pini);
00155     _pcurr = _pini;
00156     
00157     _iniparameters.ResizeTo(_nPar,1);
00158     _iniparameters(0,0) = r;
00159     _iniparameters(1,0) = theta;
00160     _iniparameters(2,0) = phi;
00161     _iniparameters(3,0) = 1.;
00162     
00163     _parameters.ResizeTo(_nPar,1);
00164     _parameters = _iniparameters;
00165 
00166     _u1.SetXYZ( TMath::Cos(phi)*TMath::Sin(theta), TMath::Sin(phi)*TMath::Sin(theta), TMath::Cos(theta) );
00167     _u2.SetXYZ( TMath::Cos(phi)*TMath::Cos(theta), TMath::Sin(phi)*TMath::Cos(theta), -1.*TMath::Sin(theta) );
00168     _u3.SetXYZ( -1.*TMath::Sin(phi), TMath::Cos(phi), 0. );
00169 
00170   }
00171 
00172 }

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

Implements TAbsFitParticle.

Definition at line 215 of file TFitParticleESpher.cc.

References TAbsFitParticle::_nPar, and TAbsFitParticle::_pini.

00215                                                                  {
00216   // Returns the parameters corresponding to the given 
00217   // 4vector
00218 
00219   // retrieve parameters
00220   TMatrixD* tparams = new TMatrixD( _nPar, 1 );
00221   (*tparams)(0,0) = vec.P();
00222   (*tparams)(1,0) = vec.Theta();
00223   (*tparams)(2,0) = vec.Phi();
00224   (*tparams)(3,0) = vec.E()/_pini.E();
00225 
00226   return tparams;
00227 
00228 }

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:33:23 2009 for CMSSW by  doxygen 1.5.4