CMS 3D CMS Logo

Public Member Functions | Protected Member Functions

TFitParticleMCCart Class Reference

#include <TFitParticleMCCart.h>

Inheritance diagram for TFitParticleMCCart:
TAbsFitParticle TNamed

List of all members.

Public Member Functions

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

Protected Member Functions

void init (TVector3 *p, Double_t M, const TMatrixD *theCovMatrix)

Detailed Description

Definition at line 9 of file TFitParticleMCCart.h.


Constructor & Destructor Documentation

TFitParticleMCCart::TFitParticleMCCart ( )

Definition at line 27 of file TFitParticleMCCart.cc.

References init().

Referenced by clone().

  :TAbsFitParticle()
{ 
  init( 0, 0., 0);
}
TFitParticleMCCart::TFitParticleMCCart ( const TFitParticleMCCart fitParticle)

Definition at line 33 of file TFitParticleMCCart.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;

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

Definition at line 52 of file TFitParticleMCCart.cc.

References init().

  :TAbsFitParticle()
{ 
  init(p, M, theCovMatrix);
}
TFitParticleMCCart::TFitParticleMCCart ( const TString &  name,
const TString &  title,
TVector3 *  p,
Double_t  M,
const TMatrixD *  theCovMatrix 
)

Definition at line 58 of file TFitParticleMCCart.cc.

References init().

  :TAbsFitParticle(name, title)
{ 
  init(p, M, theCovMatrix);
}
TFitParticleMCCart::~TFitParticleMCCart ( ) [virtual]

Definition at line 77 of file TFitParticleMCCart.cc.

                                        {

}

Member Function Documentation

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

Implements TAbsFitParticle.

Definition at line 94 of file TFitParticleMCCart.cc.

References TAbsFitParticle::_nPar, TAbsFitParticle::_pini, X, and Gflash::Z.

                                                                     {
  // Calculates a 4vector corresponding to the given
  // parameter values

  if (params == 0) {
    return 0;
  }

  if ( params->GetNcols() != 1 || params->GetNrows() !=_nPar ) {
    edm::LogError ("WrongMatrixSize")
      << GetName() << "::calc4Vec - Parameter matrix has wrong size.";
    return 0;
  }
  
  Double_t X = (*params)(0,0);
  Double_t Y = (*params)(1,0); 
  Double_t Z = (*params)(2,0); 
  Double_t E = TMath::Sqrt(  X*X + Y*Y + Z*Z + _pini.M2() );

  TLorentzVector* vec = new TLorentzVector( X, Y, Z, E );
  return vec;

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

Implements TAbsFitParticle.

Definition at line 65 of file TFitParticleMCCart.cc.

References TFitParticleMCCart().

                                                                  {
  // Returns a copy of itself

  TAbsFitParticle* myclone = new TFitParticleMCCart( *this );
  if ( newname.Length() > 0 ) myclone->SetName(newname);
  return myclone;

}
TMatrixD * TFitParticleMCCart::getDerivative ( ) [virtual]

Implements TAbsFitParticle.

Definition at line 163 of file TFitParticleMCCart.cc.

References TAbsFitParticle::_parameters, and TAbsFitParticle::_pcurr.

                                            {
  // returns derivative dP/dy with P=(p,E) and y=(px, py, pz) 
  // the free parameters of the fit. The columns of the matrix contain 
  // (dP/dpx, dP/dpy, ...).

  TMatrixD* DerivativeMatrix = new TMatrixD(4,3);
  (*DerivativeMatrix) *= 0.;

  //1st column: dP/dx
  (*DerivativeMatrix)(0,0) =  1.;
  (*DerivativeMatrix)(1,0) =  0.;
  (*DerivativeMatrix)(2,0) =  0.;
  (*DerivativeMatrix)(3,0) = _parameters(0,0)/_pcurr.E();

  //2nd column: dP/dy
  (*DerivativeMatrix)(0,1) = 0.;
  (*DerivativeMatrix)(1,1) = 1.;
  (*DerivativeMatrix)(2,1) = 0.;
  (*DerivativeMatrix)(3,1) = _parameters(1,0)/_pcurr.E();

  //3rd column: dP/dz
  (*DerivativeMatrix)(0,2) = 0.;
  (*DerivativeMatrix)(1,2) = 0.;
  (*DerivativeMatrix)(2,2) = 1.;
  (*DerivativeMatrix)(3,2) = _parameters(2,0)/_pcurr.E();

  return DerivativeMatrix;  

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

Definition at line 85 of file TFitParticleMCCart.cc.

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

Referenced by TFitParticleMCCart().

                                                                                   {

  _nPar = 3;
  setIni4Vec(p, M);
  setCovMatrix(theCovMatrix);
  
}
void TFitParticleMCCart::setIni4Vec ( const TLorentzVector *  pini) [virtual]

Implements TAbsFitParticle.

Definition at line 118 of file TFitParticleMCCart.cc.

Referenced by init().

                                                              {
  // Set the initial 4vector. Will also set the 
  // inital parameter values 

  TVector3 vec( pini->Vect() );
  setIni4Vec( &vec, pini->M() );

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

Definition at line 127 of file TFitParticleMCCart.cc.

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

                                                                 {
  // Set the initial 4vector. Will also set the 
  // inital parameter values

  if ( p == 0 ) {

    _iniparameters.ResizeTo(_nPar,1);
    _iniparameters(0,0) = 0.;
    _iniparameters(1,0) = 0.;
    _iniparameters(2,0) = 0.;
    _parameters.ResizeTo(_nPar,1);
    _parameters = _iniparameters;

    _pini.SetXYZM( 0., 0., 0., M);
    _pcurr = _pini;

  } else {

    _iniparameters.ResizeTo(_nPar,1);
    _iniparameters(0,0) = p->x();
    _iniparameters(1,0) = p->y();
    _iniparameters(2,0) = p->z();
    _parameters.ResizeTo(_nPar,1);
    _parameters = _iniparameters;

    _pini.SetXYZM( p->x(), p->y(), p->z(), M);
    _pcurr = _pini;

    _u1.SetXYZ( 1., 0., 0. );
    _u2.SetXYZ( 0., 1., 0. );
    _u3.SetXYZ( 0., 0., 1. );
  
  }

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

Implements TAbsFitParticle.

Definition at line 193 of file TFitParticleMCCart.cc.

References TAbsFitParticle::_nPar.

                                                                 {
  // Returns the parameters corresponding to the given 
  // 4vector

  TVector3 vec3( vec.Vect() );

  // retrieve parameters
  TMatrixD* tparams = new TMatrixD( _nPar, 1 );
  (*tparams)(0,0) = vec3.x();
  (*tparams)(1,0) = vec3.y();
  (*tparams)(2,0) = vec3.z();

  return tparams;

}