CMS 3D CMS Logo

TFitParticleEtEtaPhi Class Reference

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

Inheritance diagram for TFitParticleEtEtaPhi:

TAbsFitParticle

List of all members.

Public Member Functions

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

Protected Member Functions

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


Detailed Description

Definition at line 10 of file TFitParticleEtEtaPhi.h.


Constructor & Destructor Documentation

TFitParticleEtEtaPhi::TFitParticleEtEtaPhi (  ) 

Referenced by clone().

TFitParticleEtEtaPhi::TFitParticleEtEtaPhi ( const TFitParticleEtEtaPhi fitParticle  ) 

Definition at line 37 of file TFitParticleEtEtaPhi.cc.

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

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

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

Definition at line 56 of file TFitParticleEtEtaPhi.cc.

References init().

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

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

Definition at line 62 of file TFitParticleEtEtaPhi.cc.

References init().

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

TFitParticleEtEtaPhi::~TFitParticleEtEtaPhi (  )  [virtual]

Definition at line 81 of file TFitParticleEtEtaPhi.cc.

00081                                             {
00082 
00083 }


Member Function Documentation

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

Implements TAbsFitParticle.

Definition at line 96 of file TFitParticleEtEtaPhi.cc.

References TAbsFitParticle::_nPar, eta, phi, and X.

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

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

Implements TAbsFitParticle.

Definition at line 69 of file TFitParticleEtEtaPhi.cc.

References TFitParticleEtEtaPhi().

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

TMatrixD * TFitParticleEtEtaPhi::getDerivative (  )  [virtual]

Implements TAbsFitParticle.

Definition at line 171 of file TFitParticleEtEtaPhi.cc.

References TAbsFitParticle::_parameters, eta, and phi.

00171                                               {
00172   // returns derivative dP/dy with P=(p,E) and y=(et, eta, phi) 
00173   // the free parameters of the fit. The columns of the matrix contain 
00174   // (dP/d(et), dP/d(eta), dP/d(phi)).
00175 
00176   TMatrixD* DerivativeMatrix = new TMatrixD(4,3);
00177   (*DerivativeMatrix) *= 0.;
00178 
00179   Double_t et = _parameters(0,0);
00180   Double_t eta = _parameters(1,0);
00181   Double_t phi = _parameters(2,0);
00182 
00183   //1st column: dP/d(et)
00184   (*DerivativeMatrix)(0,0) = TMath::Cos(phi);
00185   (*DerivativeMatrix)(1,0) = TMath::Sin(phi);
00186   (*DerivativeMatrix)(2,0) = TMath::SinH(eta);
00187   (*DerivativeMatrix)(3,0) = TMath::CosH(eta);
00188 
00189   //2nd column: dP/d(eta)
00190   (*DerivativeMatrix)(0,1) = 0.;
00191   (*DerivativeMatrix)(1,1) = 0.;
00192   (*DerivativeMatrix)(2,1) = et*TMath::CosH(eta);
00193   (*DerivativeMatrix)(3,1) = et*TMath::SinH(eta);
00194 
00195    //3rd column: dP/d(phi)
00196   (*DerivativeMatrix)(0,2) = -1.*et*TMath::Sin(phi);
00197   (*DerivativeMatrix)(1,2) = et*TMath::Cos(phi);
00198   (*DerivativeMatrix)(2,2) = 0.;
00199   (*DerivativeMatrix)(3,2) = 0.;
00200 
00201   return DerivativeMatrix;
00202 
00203 }

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

Definition at line 88 of file TFitParticleEtEtaPhi.cc.

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

Referenced by TFitParticleEtEtaPhi().

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

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

Implements TAbsFitParticle.

Definition at line 124 of file TFitParticleEtEtaPhi.cc.

References TAbsFitParticle::_iniparameters, TAbsFitParticle::_nPar, TAbsFitParticle::_parameters, TAbsFitParticle::_pcurr, TAbsFitParticle::_pini, TAbsFitParticle::_u1, TAbsFitParticle::_u2, TAbsFitParticle::_u3, eta, phi, and funct::sin().

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     
00141     _parameters.ResizeTo(_nPar,1);
00142     _parameters(0,0) = 0.;
00143     _parameters(1,0) = 0.;
00144     _parameters(2,0) = 0.;   
00145     
00146   } else {
00147     
00148     Double_t et = pini->E()*fabs(sin(pini->Theta()));
00149     Double_t eta = pini->Eta();
00150     Double_t phi = pini->Phi();
00151     
00152     _pini = (*pini);
00153     _pcurr = _pini;
00154     
00155     _iniparameters.ResizeTo(_nPar,1);
00156     _iniparameters(0,0) = et;
00157     _iniparameters(1,0) = eta;
00158     _iniparameters(2,0) = phi;
00159     
00160     _parameters.ResizeTo(_nPar,1);
00161     _parameters = _iniparameters;
00162 
00163     _u1.SetXYZ( TMath::Cos(phi), TMath::Sin(phi), 0.); // the base vector of Et
00164     _u2.SetXYZ( -1.*TMath::Cos(phi)*TMath::TanH(eta), -1.*TMath::Sin(phi)*TMath::TanH(eta), 1./TMath::CosH(eta) );// the base vector of Eta ( same as the base vector for Theta)
00165     _u3.SetXYZ( -1.*TMath::Sin(phi), TMath::Cos(phi), 0. );// the base vector of Phi
00166 
00167   }
00168 
00169 }

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

Implements TAbsFitParticle.

Definition at line 205 of file TFitParticleEtEtaPhi.cc.

References TAbsFitParticle::_nPar, and funct::sin().

00205                                                                    {
00206   // Returns the parameters corresponding to the given 
00207   // 4vector
00208 
00209   // retrieve parameters
00210   TMatrixD* tparams = new TMatrixD( _nPar, 1 );
00211   (*tparams)(0,0) = vec.E()*fabs(sin(vec.Theta()));
00212   (*tparams)(1,0) = vec.Eta();
00213   (*tparams)(2,0) = vec.Phi();
00214 
00215   return tparams;
00216 
00217 }


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