Removed completelly high noise at channel > 700 (Marian)
[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   fNPads(0),
93   fNTimeBins(0),
94   fNBins(0),
95   fGraph(0)
96 {
97   //
98   // constructor of the info
99   //
100   for (Int_t i=0;i<25;i++){
101     fMatrix[i]=matrix[i];
102   }
103   fNPads=0;
104   fNTimeBins=0;
105   Int_t center = 5+5+2;
106   for (Int_t i=-2; i<=2;i++) if (matrix[center+i]>0) fNTimeBins++;
107   for (Int_t i=-2; i<=2;i++) if (matrix[center+i*5]>0) fNPads++;
108   fNBins = nbins;
109   fGraph = 0;
110   if (fNBins>0) {
111     fGraph = new Float_t[fNBins];
112     for (Int_t i=0;i<fNBins; i++){
113       fGraph[i] = graph[i];
114     }
115   }
116 }
117
118 AliTPCclusterInfo& AliTPCclusterInfo::operator=(const AliTPCclusterInfo& info){
119   //
120   // assignment operator
121   //
122   if (this != &info) {
123     new (this) AliTPCclusterInfo(info);
124   }
125   return *this;
126 }
127
128
129 UChar_t AliTPCclusterInfo::GetNPads(Float_t threshold) const { 
130   //
131   //
132   Int_t nPads=0;
133   Int_t center = 5+5+2;
134   for (Int_t i=-2; i<=2;i++) if (fMatrix[center+i*5]>threshold) nPads++;
135   return nPads;
136 }
137
138 UChar_t AliTPCclusterInfo::GetNTimeBins(Float_t threshold) const { 
139   //
140   //
141   //
142   Int_t nTimeBins=0;
143   Int_t center = 5+5+2;
144   for (Int_t i=-2; i<=2;i++) if (fMatrix[center+i]>threshold) nTimeBins++;
145   return nTimeBins;
146 }
147
148
149
150
151 AliTPCclusterInfo::~AliTPCclusterInfo(){
152   //
153   // destructor
154   //
155   if (fGraph)  delete [] fGraph;
156 }
157