]> git.uio.no Git - u/mrichter/AliRoot.git/blob - FMD/AliFMDCalibGain.cxx
Starting a collection of QA/Comparison macros
[u/mrichter/AliRoot.git] / FMD / AliFMDCalibGain.cxx
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 /* $Id$ */
16 /** @file    AliFMDCalibGain.cxx
17     @author  Christian Holm Christensen <cholm@nbi.dk>
18     @date    Sun Mar 26 18:30:02 2006
19     @brief   Per strip gain calibration 
20 */
21 //____________________________________________________________________
22 //                                                                          
23 // Gain value and width for each strip in the FMD. 
24 // Foo 
25 // Bar 
26 // Baz
27 // Gnus
28 //
29 #include "AliFMDCalibGain.h"    // ALIFMDCALIBGAIN_H
30 #include <iostream>
31 #include <TString.h>
32 #include <AliLog.h>
33 #include "AliFMDDebug.h"
34
35 #include "AliFMDBoolMap.h"
36
37
38 //____________________________________________________________________
39 ClassImp(AliFMDCalibGain)
40 #if 0
41   ; // This is here to keep Emacs for indenting the next line
42 #endif
43
44 //____________________________________________________________________
45 AliFMDCalibGain::AliFMDCalibGain()
46   : fValue(0), // nDet == 0 mean 51200 slots
47     fThreshold(-1.)
48 {
49   // CTOR
50   fValue.Reset(-1.);
51   fThreshold = -1.;
52 }
53
54 //____________________________________________________________________
55 AliFMDCalibGain::AliFMDCalibGain(const AliFMDCalibGain& o)
56   : TObject(o), 
57     fValue(o.fValue), 
58     fThreshold(o.fThreshold)
59 {
60   // Copy CTOR 
61 }
62
63 //____________________________________________________________________
64 AliFMDCalibGain&
65 AliFMDCalibGain::operator=(const AliFMDCalibGain& o)
66 {
67   // Assignment operator 
68   if (&o == this) return *this; 
69   fValue     = o.fValue;
70   fThreshold = o.fThreshold;
71   return (*this);
72 }
73
74 //____________________________________________________________________
75 void
76 AliFMDCalibGain::Set(UShort_t det, Char_t ring, UShort_t sec, 
77                      UShort_t str, Float_t val)
78 {
79   // Set the value for a strip 
80   if (fValue.CheckIndex(det, ring, sec, str) < 0) return;
81   fValue(det, ring, sec, str) = val;
82 }
83
84 //____________________________________________________________________
85 Float_t
86 AliFMDCalibGain::Value(UShort_t det, Char_t ring, UShort_t sec, 
87                        UShort_t str)
88 {
89   // Get the value for a strip 
90   return fValue(det, ring, sec, str);
91 }
92
93 //____________________________________________________________________
94 namespace {
95   struct MakeDead : public AliFMDMap::ForOne
96   {
97     MakeDead(AliFMDBoolMap* dead, Float_t min, Float_t max) 
98       : fDead(dead), fMin(min), fMax(max), fCount(0)
99     {}
100     MakeDead(const MakeDead& other) 
101       : AliFMDMap::ForOne(other),
102         fDead(other.fDead), fMin(other.fMin), fMax(other.fMax), 
103         fCount(other.fCount)
104     {}
105     MakeDead& operator=(const MakeDead& other) 
106     { 
107       if (&other == this) return *this; 
108       fDead   = other.fDead;
109       fMin    = other.fMin;
110       fMax    = other.fMax;
111       fCount  = other.fCount;
112       return *this;
113     }
114     Bool_t operator()(UShort_t d, Char_t r, UShort_t s, UShort_t t, Float_t v)
115     {
116       AliDebugGeneral("AliFMDCalibGain::MakeDeadMap", 100, 
117                       Form("Checking if gain of FMD%d%c[%2d,%3d]=%f "
118                            "is out of range [%f,%f]", 
119                            d, r, s, t, v, fMin, fMax));
120       if (v > fMax || v < fMin) {
121         fDead->operator()(d,r,s,t) = kTRUE;
122         fCount++;
123       }
124       return kTRUE;
125     }
126     Bool_t operator()(UShort_t, Char_t, UShort_t, UShort_t, Int_t) 
127     { return kFALSE; }
128     Bool_t operator()(UShort_t, Char_t, UShort_t, UShort_t, UShort_t) 
129     { return kFALSE; }
130     Bool_t operator()(UShort_t, Char_t, UShort_t, UShort_t, Bool_t)
131     { return kFALSE; }
132     AliFMDBoolMap* fDead;
133     Float_t        fMin;
134     Float_t        fMax;
135     Int_t          fCount;
136   };
137 }
138
139 //____________________________________________________________________
140 AliFMDBoolMap*
141 AliFMDCalibGain::MakeDeadMap(Float_t min, Float_t max, 
142                              AliFMDBoolMap* dead) const
143 {
144   if (!dead) { 
145     dead = new AliFMDBoolMap(0,0,0,0);
146     dead->Reset(kFALSE);
147   }
148   MakeDead dm(dead, min, max);
149   fValue.ForEach(dm);
150   AliFMDDebug(1, ("Found a total of %d dead channels", dm.fCount));
151   return dead;
152 }
153
154 //____________________________________________________________________
155 Bool_t
156 AliFMDCalibGain::ReadFromFile(std::istream& in)
157 {
158   //Get header (how long is it ?)
159   TString header;
160   header.ReadLine(in);
161   header.ToLower();
162   if(!header.Contains("gains")) {
163     AliError("File does not contain gains!");
164     return kFALSE;;
165   }
166
167   // Read column headers
168   header.ReadLine(in);
169   
170   int lineno  = 2;
171   // Read until EOF 
172   while(in.peek()!=EOF) {
173     if(in.bad()) { 
174       AliError(Form("Bad read at line %d of input", lineno));
175       break;
176     }
177     UShort_t det, sec, strip;
178     Char_t ring;
179     
180     Float_t gain,error,  chi2ndf;
181     Char_t c[6];
182     
183     in >> det      >> c[0] 
184        >> ring     >> c[1]
185        >> sec      >> c[2]
186        >> strip    >> c[3]
187        >> gain     >> c[4]
188        >> error    >> c[5]
189        >> chi2ndf;
190     lineno++;
191     Set(det,ring,sec,strip,gain);
192   }
193   return kTRUE;
194 }
195 //____________________________________________________________________
196 //
197 // EOF
198 //