class AliITSAlignMille2: public TObject

 main constructor that takes input from configuration file

 Destructor

 read new record from config file

 check the correctness of the record

 return 0 if success
        1 if error in module index or voluid

 build the hieararhy of the modules to align

 set the current supermodule
 new meaning

 set minimum number of points in specific detector or layer
 where = LAYER or DETECTOR
 ndet = detector number: 1-6 for LAYER and 1-3 for DETECTOR (SPD=1, SDD=2, SSD=3)
 updw = 1 for Y>0, -1 for Y<0, 0 if not specified
 nreqpts = minimum number of points of that type

 index from symname

 index from volume ID

 volume ID from symname
 works for sensitive volumes only

 volume ID from index

 initialize ideal geometry

 Load the global deltas from this file. The local gaussian constraints on some modules
 will be defined with respect to the deltas from this reference file, converted to local
 delta format. Note: conversion to local format requires reloading the geometry!

 perform global initialization

 Constrain equation defined by par to value

 Initialize global parameters with par array

 Parameter iPar is encourage to vary in [-value;value].
 If value == 0, parameter is fixed

 Reset the derivative vectors

 apply prealignment to ideal geometry

 quality factors from prealignment

 create a new AliTrackPointArray keeping only defined modules
 move points according to a given prealignment, if any
 sort alitrackpoints w.r.t. global Y direction, if cosmics

 sort alitrackpoints w.r.t. global Y direction

 get current layer id

 initialize geometry parameters for a given detector
 for current cluster (fCluster)
 fGlobalInitParam[] is set as:
    [tx,ty,tz,psi,theta,phi]
    (old was [tx,ty,tz,theta,psi,phi] ROOT's angles...)
 *** At the moment: using Raffalele's angles definition ***

 return 0 if success
 If module is found but has no parameters to vary, return 1

 print current status

 return pointer to a defined supermodule
 return NULL if error

 return pointer to a defined supermodule
 return NULL if error

 return pointer to a defined/contained supermodule
 return NULL otherwise

 get delta from prealignment for given volume

 get delta with respect to which the constraint is declared

 Initialize Riemann Fitter for current track
 return kFALSE if error

 initialize local parameters with different methods
 for current track (fTrack)

 checks if supermodule with this symname is defined and return the internal index
 return -1 if not.

 checks if module with this symname is defined and return the internal index
 return -1 if not.

 checks if supermodule 'voluid' is defined and return the internal index
 return -1 if not.

 checks if the sensitive module 'voluid' is contained inside a supermodule
 and return the internal index of the last identified supermodule
 return -1 if error
 IMPORTANT: always start from the end to start from the sensors

 checks if the sensitive module 'voluid' is contained inside a supermodule
 and return the internal index of the last identified supermodule
 return -1 if error
 IMPORTANT: always start from the end to start from the sensors

 checks if AliTrackPoints belongs to defined modules
 return number of good poins
 return 0 if not enough points

 Process track; Loop over hits and set local equations
 here 'track' is a AliTrackPointArray
 return 0 if success;

 Fit the track with selected constraints

 calculate track intersection point in local coordinates
 according with a given set of parameters (local(4) and global(6))
 and fill fPintLoc/Glo
    local are:   pgx0, pgz0, ugx, ugz   OR   riemann fitters pars
    global are:  tx,ty,tz,psi,theta,phi (Raff's delta angles in deg.)
 return 0 if success

 calculate numerically (ROOT's style) the derivatives for
 local X intersection  and local Z intersection
 parlist: local  (islpar=kTRUE)  pgx0, pgz0, ugx0, ugz0  OR riemann's params
          global (islpar=kFALSE) tx, ty, tz, psi, theta, phi (Raf's angles in deg)
 return 0 if success

 Define local equation for current cluster in X and Z coor.
 and store them to memory
 return -1 in case of failure to build some equation
         0 if no free global parameters were found but local eq is built
         1 if both local and global eqs are built

 store first intersection point

 Define local equation for current cluster in X Y and Z coor.
 and store them to memory
 return -1 in case of failure to build some equation
         0 if no free global parameters were found but local eq is built
         1 if both local and global eqs are built

 Set local equations with data stored in m
 return 0 if success

 Call global fit; Global parameters are stored in parameters

 Print global parameters

 load definitions of supermodules from a root file
 return 0 if success

 require that sum of modifications for the childs of this module is = val, i.e.
 the internal corrections moves the module as a whole by fixed value (0 by default).
 pattern is the bit pattern for the parameters to constrain

 require that median of the modifications for the childs of this module is = val, i.e.
 the internal corrections moves the module as a whole by fixed value (0 by default)
 module the outliers.
 pattern is the bit pattern for the parameters to constrain
 The difference between the mean and the median will be transfered to the parent

 require that median of the modifications for the supermodules which have no parents is = val, i.e.
 the corrections moves the whole setup by fixed value (0 by default) modulo the outliers.
 pattern is the bit pattern for the parameters to constrain

 require that median of the modifications for the supermodules which have no parents is = val, i.e.
 the corrections moves the whole setup by fixed value (0 by default) modulo the outliers.
 pattern is the bit pattern for the parameters to constrain

 apply constraint on parameters in the local frame

 apply the constraint on the local corrections of a list of modules

 apply constriants which cannot be imposed after the fit

 apply constraints which can be imposed after the fit

 require that sum of modifications for the childs of this module is = val, i.e.
 the internal corrections moves the module as a whole by fixed value (0 by default).
 pattern is the bit pattern for the parameters to constrain

 require that median of the modifications for the supermodules which have no parents is = val, i.e.
 the corrections moves the whole setup by fixed value (0 by default) modulo the outliers.
 pattern is the bit pattern for the parameters to constrain

 require that median or mean of the modifications for the childs of this module is = val, i.e.
 the internal corrections moves the module as a whole by fixed value (0 by default)
 module the outliers.
 pattern is the bit pattern for the parameters to constrain
 The difference between the mean and the median will be transfered to the parent

 require that median or mean of modifications for the supermodules which have no parents is = val, i.e.
 the corrections moves the whole setup by fixed value (0 by default).
 pattern is the bit pattern for the parameters to constrain

 check if par of the module participates in some constraint, and set the flag for their types

 check if parameter par is varied for this module or its children up to the level depth

 check if parameter par is varied and is not subjected to gaussian constraint for the children up to the level depth

 obtain drift time corrected for t0

 obtain corrected drift speed

 are there fixed modules with no parent (normally in such a case
 the constraints on the orphans should not be applied

 convert params in local frame to global one

 convert params in global frame to local one

 set Bz value

 extract calibration information used for TrackPointArray creation from run info

 extract the path for specific CDB path from user info. If it is the same as already loaded, set corresponing bit

 load SDD response

 load VDrift object

 Load SDD correction map

 Load SDD correction map for prealignment from current CDB

 load vertex constraint

 load ITS geom deltas

 build arrays for the fast access to sensor matrices from their sensor ID

 build arrays for the fast access to sensor original matrices (used for production)

 suppress constraint

 constrain q and pT of the helical fit of the track (should be set before process.track)

 constrain charge and curvature of the helical fit of the track (should be set before process.track)

 Create \Deltas for every explicitly or implicitly (via non-alignable volumes) varied
 or prealigned module.
 If the module has inded J in the hierarchy of alignable volumes (0 - the top, most
 coarse level), then its Delta is expressed via MP2 \deltas (in global frame) and
 prealignment \DeltaP's as:
 \Delta_J = Y X Y^-1
 where X = \delta_J * \DeltaP_J
 Y = Prod_{K=0,J-1} \delta_K
 Note that \delta_L accounts not only for its own correction but also of all non-alignable
 modules in the hierarchy chain from L up to the closest alignable:
 while (parent && !parent->IsAlignable()) {
   \delta_L->MultiplyLeft( \delta_parent );
   parent = parent->GetParent();
 }

 create object with SDD repsonse (t0 and vdrift corrections) accounting for
 eventual precalibration

 if there was a precalibration provided, copy it to new arrray

 Use provided UserInfo to
 load the initial calib parameters (geometry, SDD response...)
 Can be used if set of data was processed with different calibration

 Load the initial calib parameters (geometry, SDD response...)
 Can be used if set of data was processed with different calibration

 calculate the locid row of the jacobian for transformation of the local coordinate to global at current point

 impose equality of Left/Right sides VDrift correction for SDD

 extract the drift information from SDD track point

 creates dummy module for vertex constraint

 return object from the OCDB

 set vertex for constraint

 convert prealignment deltas from old geometry to new one
 NOTE: the target geometry must be loaded at time this method is called

 NOTE: This method can be ONLY used when as a prealignment deltas those used for the production
 of trackpoints (e.g. extracted from the UserInfo).
 The prealignment deltas provided by user via config file must be already converted to target geometry:
 this can be done externally using the macro ConvertDeltas.C

 delta_j_new = delta_j_old * Xj_old * Xj_new^-1
 where X = Prod{delta_i,i=j-1:0} M_j
 with j - the level of the alignable volume in the hierarchy, M - corresponding ideal matrix
 Note that delta_j * Xj is equal to final (misaligned) matrix of corresponding geometry, G_j.
 Since this method is used ONLY in the case where the prealignment deltas are equal to production deltas,
 we have already loaded G_j_old in the fConvAlgMatOld (filled in the CacheMatricesOrig)
 Hence, delta_j_new = G_j_old * Xj_new^-1

 Used only for the deltas conversion from one geometry to another
 Sort the matrices according to hiearachy (coarse -> fine)

 Used only for the deltas conversion from one geometry to another
 True if matJ is higher in hierarchy than

 find the delta for given module

 configuration methods

 geometry stuffs

 fitting methods

 methods for point unbiasing (via scaling its inverted cov.matrix)

 debug stuffs

 pepo270809

 endpepo270809

 pepo
 flag for AliITSAlignMille compatibility

 modified existing methods

 old methods recovered

 moved from private to public

 millepede methods