Correct sign for calculated b_y (outside measured region).
[u/mrichter/AliRoot.git] / TPC / AliTPCclusterInfo.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
16 //-------------------------------------------------------
17 //          Implementation of the TPC cluster debug information
18 //
19 //   Origin: Marian Ivanov   Marian.Ivanov@cern.ch
20 //   
21 //   Additional cluster information to monitor clustering performance
22 //   and to extract a features of detector response
23 //   Information attached to the default cluster
24 //   ONLY in DEBUG MODE
25 //    
26 //-------------------------------------------------------
27
28 /* $Id$ */
29
30 #include "AliTPCclusterInfo.h"
31 #include "AliLog.h"
32
33 ClassImp(AliTPCclusterInfo)
34
35
36 AliTPCclusterInfo::AliTPCclusterInfo():
37   fNPads(0),
38   fNTimeBins(0),
39   fNBins(0),
40   fGraph(0)
41 {
42   //
43   // default constructor
44   //  
45   for (Int_t i=0; i<25;i++){
46     fMatrix[i] = i; 
47   }
48 }
49
50 AliTPCclusterInfo::AliTPCclusterInfo(const  AliTPCclusterInfo & info):
51   TObject(info),
52   fNPads(info.fNPads),
53   fNTimeBins(info.fNTimeBins),
54   fNBins(info.fNBins),
55   fGraph(0)    
56 {
57   //
58   // copy constructor
59   //
60   //  AliInfo("Copy constructor\n");
61   for (Int_t i=0; i<25;i++){
62     fMatrix[i] = info.fMatrix[i]; 
63   }
64   if (info.fGraph) fGraph = new Float_t[fNBins];
65   for (Int_t i=0;i<fNBins; i++){
66     fGraph[i] = info.fGraph[i];
67   }
68   
69 }
70
71
72 AliTPCclusterInfo::AliTPCclusterInfo(Bool_t extend):
73   fNPads(0),
74   fNTimeBins(0),
75   fNBins(0),
76   fGraph(0)
77 {
78   //
79   // allocate dummy graph - neccessary for IO part to use automatic branching
80   //  
81   for (Int_t i=0; i<25;i++){
82     fMatrix[i] = i; 
83   }
84   if (extend){
85     fNBins = 1;
86     fGraph  = new Float_t[1];
87     fGraph[0]=-1;
88   }
89 }
90
91 AliTPCclusterInfo::AliTPCclusterInfo(Float_t *matrix, Int_t nbins, Float_t* graph){
92   //
93   // constructor of the info
94   //
95   for (Int_t i=0;i<25;i++){
96     fMatrix[i]=matrix[i];
97   }
98   fNPads=0;
99   fNTimeBins=0;
100   Int_t center = 5+5+2;
101   for (Int_t i=-2; i<=2;i++) if (matrix[center+i]>0) fNTimeBins++;
102   for (Int_t i=-2; i<=2;i++) if (matrix[center+i*5]>0) fNPads++;
103   fNBins = nbins;
104   fGraph = 0;
105   if (fNBins>0) {
106     fGraph = new Float_t[fNBins];
107     for (Int_t i=0;i<fNBins; i++){
108       fGraph[i] = graph[i];
109     }
110   }
111 }
112
113 UChar_t AliTPCclusterInfo::GetNPads(Float_t threshold) const { 
114   //
115   //
116   Int_t nPads=0;
117   Int_t center = 5+5+2;
118   for (Int_t i=-2; i<=2;i++) if (fMatrix[center+i*5]>threshold) nPads++;
119   return nPads;
120 }
121
122 UChar_t AliTPCclusterInfo::GetNTimeBins(Float_t threshold) const { 
123   //
124   //
125   //
126   Int_t nTimeBins=0;
127   Int_t center = 5+5+2;
128   for (Int_t i=-2; i<=2;i++) if (fMatrix[center+i]>0) nTimeBins++;
129   return nTimeBins;
130 }
131
132
133
134
135 AliTPCclusterInfo::~AliTPCclusterInfo(){
136   //
137   // destructor
138   //
139   if (fGraph)  delete [] fGraph;
140 }
141