]>
Commit | Line | Data |
---|---|---|
10757ee9 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | //////////////////////////////////////////////////////////////////////////// | |
17 | // | |
18 | // === Class for fitting properties specific to pad regions === | |
19 | // | |
20 | // For each pad size region a separate TLinearFitter object is assigned. | |
21 | // Commonly used functions such as getting the center of a pad size | |
22 | // region or visualization functions are provided. Also, choosing the | |
23 | // segment and pad type is made easy due to different methods that | |
24 | // can calculate these informations from other coordinates. | |
25 | // | |
26 | //////////////////////////////////////////////////////////////////////////// | |
27 | ||
28 | #include "AliTPCFitPad.h" | |
29 | ||
30 | ClassImp(AliTPCFitPad) | |
31 | ||
10757ee9 | 32 | |
f1c2a4a3 | 33 | |
10757ee9 | 34 | |
35 | AliTPCFitPad::AliTPCFitPad(Int_t ndim, const char* formula, Option_t* opt) : | |
36 | AliTPCCalPadRegion("", ""), | |
37 | fNdim(ndim), | |
38 | fFormula(formula), | |
39 | fOpt(opt) | |
40 | { | |
41 | // | |
42 | // Constructor. The parameters are used for generating new TLinearFitter | |
43 | // objects and are described in its documentation. | |
44 | // | |
45 | } | |
46 | ||
f1c2a4a3 | 47 | AliTPCFitPad& AliTPCFitPad::operator=(const AliTPCFitPad& rhs) |
48 | { | |
49 | // | |
50 | // Assignment operator. | |
51 | // | |
52 | ||
53 | if (this != &rhs) { | |
54 | AliTPCCalPadRegion::operator=(rhs); | |
55 | fNdim = rhs.fNdim; | |
56 | fFormula = rhs.fFormula; | |
57 | fOpt = rhs.fOpt; | |
58 | } | |
59 | return *this; | |
60 | } | |
61 | ||
62 | AliTPCFitPad::~AliTPCFitPad() { | |
10757ee9 | 63 | // |
f1c2a4a3 | 64 | // Destructor. |
10757ee9 | 65 | // |
66 | ||
f1c2a4a3 | 67 | Delete(); |
10757ee9 | 68 | } |
69 | ||
f1c2a4a3 | 70 | |
71 | ||
10757ee9 | 72 | void AliTPCFitPad::Add(AliTPCFitPad* fit) { |
73 | // | |
74 | // Adds another AliTPCFitPad object to this object. The formula should be the | |
75 | // same, though it won't be checked! | |
76 | // | |
77 | ||
78 | for (UInt_t iSegment = 0; iSegment < GetNSegments(); iSegment++) { | |
79 | for (UInt_t iPadType = 0; iPadType < GetNPadTypes(); iPadType++) { | |
80 | TLinearFitter* fitter = fit->GetFitterSimple(iSegment, iPadType); | |
81 | // parameter workaround == kTRUE because it is not possible to add another | |
82 | // TLinearFitter object to a "virgin" one. Thus a dummy data point is added | |
83 | // and cleared again immediately afterwards. Due to a broken TLinearFitter | |
84 | // copy constructor this is a necessary workaround. | |
85 | if (fitter) { | |
86 | cout << "TLinearFitter::Add called for " << iSegment << ", " << iPadType << endl; | |
87 | GetFitter(iSegment, iPadType, kTRUE)->Add(fitter); | |
88 | } | |
89 | } | |
90 | } | |
91 | } | |
92 | ||
93 | TLinearFitter* AliTPCFitPad::GetFitterSimple(UInt_t segment, UInt_t padType) { | |
94 | // | |
95 | // This method returns the fitter corresponding to segment and pad type. | |
96 | // In contrast to GetFitter() no fitter will be created, if it does | |
97 | // not exist, but a null pointer is returned. | |
98 | // | |
99 | ||
100 | return (TLinearFitter*)(GetObject(segment, padType)); | |
101 | } | |
102 | ||
103 | TLinearFitter* AliTPCFitPad::GetFitter(UInt_t segment, UInt_t padType, Bool_t workaround) { | |
104 | // | |
105 | // This method returns the fitter corresponding | |
106 | // to segment and pad type. | |
107 | // If the fitter doesn't exist yet, it will be created on the fly | |
108 | // according to the parameters passed to the constructor. | |
109 | // | |
110 | // The workaround parameter should always be kFALSE. (It is only used by the Add method.) | |
111 | // | |
112 | ||
113 | TLinearFitter* fitter = GetFitterSimple(segment, padType); | |
114 | if (fitter == 0) { | |
115 | SetObject(new TLinearFitter(fNdim, fFormula, fOpt), segment, padType); | |
116 | fitter = (TLinearFitter*)(GetObject(segment, padType)); | |
117 | if (workaround) { | |
118 | Double_t x[1000]; | |
119 | for (Int_t i = 0; i < fNdim; i++) x[i] = 3.141592; | |
120 | fitter->AddPoint(x, 31.41592); | |
121 | fitter->ClearPoints(); | |
122 | //cout << "workaround called for " << segment << ", " << padType << endl; | |
123 | } | |
124 | } | |
125 | return fitter; | |
126 | } | |
127 | ||
128 | Int_t AliTPCFitPad::Evaluate(Bool_t robust, Double_t frac) { | |
129 | // | |
130 | // Evaluates all fitters. Returns 0 if successful, 1 in case of an error. | |
131 | // If the robust option is set to kTRUE a robust fit is performed with frac as | |
132 | // the minimal fraction of good points (see TLinearFitter::EvalRobust for details). | |
133 | // Beware: Robust fitting is much slower! | |
134 | // | |
135 | ||
136 | Int_t returnCode = 0; | |
137 | for (UInt_t iSegment = 0; iSegment < GetNSegments(); iSegment++) { | |
138 | for (UInt_t iPadType = 0; iPadType < GetNPadTypes(); iPadType++) { | |
139 | if (TLinearFitter* fitter = GetFitterSimple(iSegment, iPadType)) { | |
140 | Int_t status = robust ? fitter->EvalRobust(frac) : fitter->Eval(); | |
141 | if (status != 0) { | |
142 | returnCode = 1; | |
143 | Error("Evaluate", "Error in evaluation of fitter in segment %d, pad region %d", iSegment, iPadType); | |
144 | } | |
145 | } | |
146 | } | |
147 | } | |
148 | return returnCode; | |
149 | } |