class AliMillepede: public TObject
\class AliMillepede Detector independent alignment class This modified C++ version is based on a C++ translation of Millepede used for LHCb Vertex Detector alignment (lhcb-2005-101), available here http://isscvs.cern.ch/cgi-bin/cvsweb.cgi/Alignment/AlignmentTools/src/?cvsroot=lhcb The original millepede fortran package is available at: http://www.desy.de/~blobel/wwwmille.html \author Javier Castillo
Function Members (Methods)
public:
| AliMillepede() | |
| AliMillepede(const AliMillepede&) | |
| virtual | ~AliMillepede() |
| void | TObject::AbstractMethod(const char* method) const |
| virtual void | TObject::AppendPad(Option_t* option = "") |
| virtual void | TObject::Browse(TBrowser* b) |
| static TClass* | Class() |
| virtual const char* | TObject::ClassName() const |
| virtual void | TObject::Clear(Option_t* = "") |
| virtual TObject* | TObject::Clone(const char* newname = "") const |
| virtual Int_t | TObject::Compare(const TObject* obj) const |
| virtual void | TObject::Copy(TObject& object) const |
| virtual void | TObject::Delete(Option_t* option = "")MENU |
| virtual Int_t | TObject::DistancetoPrimitive(Int_t px, Int_t py) |
| virtual void | TObject::Draw(Option_t* option = "") |
| virtual void | TObject::DrawClass() constMENU |
| virtual TObject* | TObject::DrawClone(Option_t* option = "") constMENU |
| virtual void | TObject::Dump() constMENU |
| virtual void | TObject::Error(const char* method, const char* msgfmt) const |
| virtual void | TObject::Execute(const char* method, const char* params, Int_t* error = 0) |
| virtual void | TObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0) |
| virtual void | TObject::ExecuteEvent(Int_t event, Int_t px, Int_t py) |
| virtual void | TObject::Fatal(const char* method, const char* msgfmt) const |
| virtual TObject* | TObject::FindObject(const char* name) const |
| virtual TObject* | TObject::FindObject(const TObject* obj) const |
| virtual Option_t* | TObject::GetDrawOption() const |
| static Long_t | TObject::GetDtorOnly() |
| virtual const char* | TObject::GetIconName() const |
| virtual const char* | TObject::GetName() const |
| virtual Int_t | GetNLocalEquations() const |
| virtual char* | TObject::GetObjectInfo(Int_t px, Int_t py) const |
| static Bool_t | TObject::GetObjectStat() |
| virtual Option_t* | TObject::GetOption() const |
| virtual Double_t | GetParError(int iPar) const |
| virtual const char* | TObject::GetTitle() const |
| virtual UInt_t | TObject::GetUniqueID() const |
| virtual Int_t | GlobalFit(double* par, double* error, double* pull) |
| virtual Bool_t | TObject::HandleTimer(TTimer* timer) |
| virtual ULong_t | TObject::Hash() const |
| virtual void | TObject::Info(const char* method, const char* msgfmt) const |
| virtual Bool_t | TObject::InheritsFrom(const char* classname) const |
| virtual Bool_t | TObject::InheritsFrom(const TClass* cl) const |
| virtual Int_t | InitMille(int nglo, int nloc, int nstd, double lResCut, double lResCutInit) |
| virtual void | TObject::Inspect() constMENU |
| void | TObject::InvertBit(UInt_t f) |
| virtual TClass* | IsA() const |
| virtual Bool_t | TObject::IsEqual(const TObject* obj) const |
| virtual Bool_t | TObject::IsFolder() const |
| Bool_t | TObject::IsOnHeap() const |
| virtual Bool_t | TObject::IsSortable() const |
| Bool_t | TObject::IsZombie() const |
| virtual Int_t | LocalFit(int n, double* localParams, Bool_t bSingleFit) |
| virtual void | TObject::ls(Option_t* option = "") const |
| void | TObject::MayNotUse(const char* method) const |
| virtual Bool_t | TObject::Notify() |
| void | TObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const |
| static void | TObject::operator delete(void* ptr) |
| static void | TObject::operator delete(void* ptr, void* vp) |
| static void | TObject::operator delete[](void* ptr) |
| static void | TObject::operator delete[](void* ptr, void* vp) |
| void* | TObject::operator new(size_t sz) |
| void* | TObject::operator new(size_t sz, void* vp) |
| void* | TObject::operator new[](size_t sz) |
| void* | TObject::operator new[](size_t sz, void* vp) |
| AliMillepede& | operator=(const AliMillepede&) |
| virtual void | TObject::Paint(Option_t* option = "") |
| virtual void | TObject::Pop() |
| virtual void | TObject::Print(Option_t* option = "") const |
| virtual Int_t | PrintGlobalParameters() const |
| virtual Int_t | TObject::Read(const char* name) |
| virtual void | TObject::RecursiveRemove(TObject* obj) |
| void | TObject::ResetBit(UInt_t f) |
| virtual void | TObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU |
| virtual void | TObject::SavePrimitive(basic_ostream<char,char_traits<char> >& out, Option_t* option = "") |
| void | TObject::SetBit(UInt_t f) |
| void | TObject::SetBit(UInt_t f, Bool_t set) |
| virtual void | TObject::SetDrawOption(Option_t* option = "")MENU |
| static void | TObject::SetDtorOnly(void* obj) |
| virtual Int_t | SetGlobalConstraint(double* dercs, double rhs) |
| virtual Int_t | SetGlobalParameter(int iPar, double vPar) |
| virtual Int_t | SetGlobalParameters(double* param) |
| virtual Int_t | SetIterations(double cutfac) |
| virtual Int_t | SetLocalEquation(double* dergb, double* derlc, double rmeas, double sigma) |
| virtual void | SetNLocalEquations(int value) |
| virtual Int_t | SetNonLinear(int index) |
| static void | TObject::SetObjectStat(Bool_t stat) |
| virtual Int_t | SetParSigma(int iPar, double sigma) |
| virtual void | TObject::SetUniqueID(UInt_t uid) |
| virtual void | ShowMembers(TMemberInspector&) |
| virtual void | Streamer(TBuffer&) |
| void | StreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b) |
| virtual void | TObject::SysError(const char* method, const char* msgfmt) const |
| Bool_t | TObject::TestBit(UInt_t f) const |
| Int_t | TObject::TestBits(UInt_t f) const |
| virtual void | TObject::UseCurrentStyle() |
| virtual void | TObject::Warning(const char* method, const char* msgfmt) const |
| virtual Int_t | TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) |
| virtual Int_t | TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const |
protected:
| virtual void | TObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
| void | TObject::MakeZombie() |
private:
| double | Chi2DoFLim(int n, int nd) |
| int | SpAVAt(double v[][fgkMaxLocalPar], double a[][fgkMaxLocalPar], double w[][fgkMaxGlobalPar], int n, int m) |
| int | SpAX(double a[][fgkMaxLocalPar], double* x, double* y, int n, int m) |
| int | SpmInv(double v[][fgkMaxGloPC], double* b, int n) |
Data Members
public:
| enum TObject::EStatusBits { | kCanDelete | |
| kMustCleanup | ||
| kObjInCanvas | ||
| kIsReferenced | ||
| kHasUUID | ||
| kCannotPick | ||
| kNoContextMenu | ||
| kInvalidObject | ||
| }; | ||
| enum TObject::[unnamed] { | kIsOnHeap | |
| kNotDeleted | ||
| kZombie | ||
| kBitMask | ||
| kSingleKey | ||
| kOverwrite | ||
| kWriteDelete | ||
| }; |
private:
| int | fCGLRow2Glo[1000] | /< Row in "used" g*l matrix to global parameter |
| double | fChi2CutFactor | /< Cut factor for chi2 cut to accept local fit |
| double | fChi2CutRef | /< Reference cut for chi2 cut to accept local fit |
| double | fDeltaPar[1000] | /< Variation of global parameters |
| TArrayD | fDerivAllEqs | /< derivative in all loc. eq. storage for iterations |
| TArrayD | fDerivLocEq | /< Table of local equation derivatives wrt. parameter |
| double | fDiagCGlo[1020] | /< Initial diagonal elements of C global matrix |
| int | fGlo2CGLRow[1000] | /< Global parameter to row in "used" g*l matrix |
| TArrayI | fIndexAllEqs | /< Index in all loc. eq. storage for iterations |
| TArrayI | fIndexLocEq | /< Table of parameter indexes in local equation |
| double | fInitPar[1000] | /< Initial global parameters |
| Bool_t | fIsNonLinear[1000] | /< Flag for non linear parameters |
| int | fIter | /< Current iteration |
| double | fLagMult[20] | /< Lagrange multipliers of constrained equations |
| TArrayI | fLocEqPlace | /< Loc. Eq. position in AllEqs storage |
| double | fMatCGlo[1020][1020] | /< Matrix C global |
| double | fMatCGloCorr[1000][1000] | /< Correction of matrix C global |
| double | fMatCGloLoc[1000][20] | /< Rectangular matrix C g*l |
| double | fMatCLoc[20][20] | /< Matrix C local |
| double | fMatDerConstr[20][1000] | /< Constrained derivatives |
| int | fMaxIter | /< Maximum number of iterations |
| Int_t | fNDerivAllEqs | /< Number of entries in fDerivAllEqs |
| Int_t | fNDerivLocEq | /< Number of entries in fDerivLocEq |
| int | fNGlobalConstraints | /< Number of constraint equations |
| int | fNGlobalPar | /< Number of global parameters |
| Int_t | fNIndexAllEqs | /< Number of entries in fIndexAllEqs |
| Int_t | fNIndexLocEq | /< Number of entries in fIndexLocEq |
| Int_t | fNLocEqPlace | /< Number of entries in fLocEqPlace |
| int | fNLocalEquations | /< Number of local equations |
| int | fNLocalFits | /< Number of local fits |
| int | fNLocalFitsRejected | /< Number of local fits rejected |
| int | fNLocalPar | /< Number of local parameters |
| int | fNStdDev | /< Number of standard deviations for chi2 cut |
| double | fResCut | /< Cut in residual for other iterartiona |
| double | fResCutInit | /< Cut in residual for first iterartion |
| double | fSigmaPar[1000] | /< Sigma of allowed variation of global parameter |
| double | fVecBGlo[1020] | /< Vector B global (parameters) |
| double | fVecBGloCorr[1000] | /< Correction of vector B global |
| double | fVecBLoc[20] | /< Vector B local (parameters) |
| static const int | fgkMaxGloCsts | /< Max. number of constraint equations |
| static const int | fgkMaxGloPC | /< fgkMaxGlobalPar+fgkMaxGloCsts |
| static const int | fgkMaxGlobalPar | /< Max. number of global parameters |
| static const int | fgkMaxLocalPar | /< Max. number of local parameters |
Class Charts
| Inheritance Chart: | ||||||||
|
Function documentation
Int_t InitMille(int nglo, int nloc, int nstd, double lResCut, double lResCutInit)
Initialization of millepede
Int_t SetParSigma(int iPar, double sigma)
Define a range [-sigma;sigma] where iPar is encourage to vary
Int_t SetLocalEquation(double* dergb, double* derlc, double rmeas, double sigma)
Write one local equation
Int_t LocalFit(int n, double* localParams, Bool_t bSingleFit)
Perform local parameters fit once all the local equations have been set
Int_t GlobalFit(double* par, double* error, double* pull)
perform global parameters fit once all the local equations have been fitted
int SpmInv(double v[][fgkMaxGloPC], double* b, int n)
Obtain solution of a system of linear equations with symmetric matrix and the inverse (using 'singular-value friendly' GAUSS pivot)
Int_t SpAVAt(double v[][fgkMaxLocalPar], double a[][fgkMaxLocalPar], double w[][fgkMaxGlobalPar], int n, int m)
multiply symmetric N-by-N matrix from the left with general M-by-N matrix and from the right with the transposed of the same general matrix to form symmetric M-by-M matrix.
Int_t SpAX(double a[][fgkMaxLocalPar], double* x, double* y, int n, int m)
multiply general M-by-N matrix A and N-vector X