+
+// OLD function, with roads automatically calculated instead from being parametrized
+// kept because it would be a better solution,
+// if one can really find the right values.
+// //__________________________________________________________________________
+// void AliMUONSegment::UpdateFromStationTrackParam(AliMUONTrackParam *TrackParam, Double_t MCSfactor, Double_t Dz1, Double_t Dz2)
+// {
+// // Fill data members with values calculated from the array of track parameters
+// // pointed to by "TrackParam" (index = 0 and 1 for first and second chambers
+// // of the station, respectively).
+// // Multiple Coulomb scattering is taking into account with "MCSfactor"
+// // corresponding to one chamber,
+// // with one chamber for the coordinate, two chambers for the angle,
+// // due to the arrangement in stations.
+// // Resolution coming from:
+// // coordinate in closest station at "Dz1",
+// // slope between closest stations, with "Dz2" interval between them,
+// // extrapolation over "Dz" from closest station.
+// // When called, "fBendingCoorReso2" and "fNonBendingCoorReso2"
+// // are assumed to be filled
+// // with the variance on bending and non bending coordinates.
+// AliMUONTrackParam *param0;
+// Double_t cReso2, sReso2;
+// param0 = &(TrackParam[0]);
+// // Bending plane
+// fBendingCoor = param0->GetBendingCoor(); // coordinate
+// fBendingSlope = param0->GetBendingSlope(); // slope
+// cReso2 = fBendingCoorReso2;
+// sReso2 = 2.0 * cReso2 / Dz2 / Dz2;
+// fBendingCoorReso2 = cReso2 + (sReso2 + MCSfactor) * Dz1 * Dz1;
+// fBendingSlopeReso2 = sReso2 + 2.0 * MCSfactor;
+// // Non bending plane
+// fNonBendingCoor = param0->GetNonBendingCoor(); // coordinate
+// fNonBendingSlope = param0->GetNonBendingSlope(); // slope
+// cReso2 = fNonBendingCoorReso2;
+// sReso2 = 2.0 * cReso2 / Dz2 / Dz2;
+// fNonBendingCoorReso2 = cReso2 + (sReso2 + MCSfactor) * Dz1 * Dz1;
+// fNonBendingSlopeReso2 = sReso2 + 2.0 * MCSfactor;
+// return;
+// }