// Setup a trackrep which uses geane
  CbmGeanePro* gePro=new CbmGeanePro();

  TVector3 pos(0.1,0.1,0.1);
  TVector3 mom(0.1,0.1,1.);
  TVector3 poserr=TVector3(0.1,0.1,0.1);
  TVector3 momerr=0.2*mom;

  DetPlane initplane(TVector3(pos),TVector3(1.,0.,0.),TVector3(0.,1.,0.));
  double q=-1;
  
 CbmTrackParP 	par(pos,mom,poserr,momerr,q,
	            initplane.getO(),initplane.getU(),initplane.getV());
  par.Print();

======================================================
Position : (0.1, 0.1, 0.1)
Slopes : dx/dz = 0.0008, dy/dz = 0.0004
q/p = -0.99
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

//      Tranform error matrix 
//      FROM   MASTER                  VARIABLES (px, py,pz, x, y, z)
//      TO     SD (transverse or local system) 
//                                     VARIABLES (1/p, v', w', v, w)
//
//      Method: the MARS system is considered as a detector system SD1 with
//              y-z as the detector plane. Hence eq (79) of the 
//              report CMS 2006/001 is used to go from canonical to SD1 variables.
//              Then, the SD1 to SD routine is used.
//              In this way the track length variation and the magnetic field
//              effects are correctly taken into account.
// 
//      Authors: A. Rotondi and A. Fontana (July 2007)

TVector3 positionsd = util.FromMARSToSDCoord(TVector3(fX, fY, fZ), forigin, fiver, fjver, fkver);
   fU = positionsd.X(); 
   fV = positionsd.Y(); 
   fW = positionsd.Z(); 

DetPlane plane(TVector3(0.0,0.0,0.0),TVector3(1.,0.,0.),TVector3(0.,1.,0.));

-   Double_t PD[3], RD[6][6], H[3], CH, SP1, PC[3], RC[3];
+   Double_t PD[3], RD[6][6], H[3], CH, SP1, PC[3], RC[15];

memset(PD,0,sizeof(PD));