adding map include file to solve compilation issue (bug https://savannah.cern.ch...
[u/mrichter/AliRoot.git] / HLT / ITS / AliHLTITSQAComponent.cxx
1 // $Id$
2 //**************************************************************************
3 //* This file is property of and copyright by the ALICE HLT Project        * 
4 //* ALICE Experiment at CERN, All rights reserved.                         *
5 //*                                                                        *
6 //* Primary Authors: Piergiorgio Cerello <cerello@to.infn.it>              *
7 //*                  for The ALICE HLT Project.                            *
8 //*                                                                        *
9 //* Permission to use, copy, modify and distribute this software and its   *
10 //* documentation strictly for non-commercial purposes is hereby granted   *
11 //* without fee, provided that the above copyright notice appears in all   *
12 //* copies and that both the copyright notice and this permission notice   *
13 //* appear in the supporting documentation. The authors make no claims     *
14 //* about the suitability of this software for any purpose. It is          *
15 //* provided "as is" without express or implied warranty.                  *
16 //**************************************************************************
17
18 /// @file   AliHLTITSQHistoComponent.cxx
19 /// @author Piergiorgio Cerello cerello@to.infn.it
20 /// @date   2009-07-03
21 /// @brief  Interface component to the ITS QA
22
23 #if __GNUC__>= 3
24 using namespace std;
25 #endif
26
27 #include "AliHLTITSQAComponent.h"
28 #include "AliHLTITSClusterDataFormat.h"
29 #include "AliCDBEntry.h"
30 #include "AliCDBManager.h"
31 #include "AliITSRecPoint.h"
32 #include <TFile.h>
33 #include <TString.h>
34 #include "TObjString.h"
35 #include "TObjArray.h"
36
37
38 //#include <stdlib.h>
39 //#include <cerrno>
40
41 /** ROOT macro for the implementation of ROOT specific class methods */
42 ClassImp(AliHLTITSQAComponent)
43
44 AliHLTITSQAComponent::AliHLTITSQAComponent()
45 :
46 fAliITSQADataMakerRec(NULL)
47 {
48   // see header file for class documentation
49   // or
50   // refer to README to build package
51   // or
52   // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
53
54 }
55
56 AliHLTITSQAComponent::~AliHLTITSQAComponent()
57 {
58   // see header file for class documentation
59 }
60
61 // Public functions to implement AliHLTComponent's interface.
62 // These functions are required for the registration process
63
64 const char* AliHLTITSQAComponent::GetComponentID()
65 {
66   // see header file for class documentation
67   
68   return "ITSClusterQA";
69 }
70
71 void AliHLTITSQAComponent::GetInputDataTypes(AliHLTComponentDataTypeList& list)
72 {
73   // see header file for class documentation
74   list.clear();
75   list.push_back( kAliHLTDataTypeTObjArray );
76 }
77
78 AliHLTComponentDataType AliHLTITSQAComponent::GetOutputDataType()
79 {
80   // see header file for class documentation
81   return kAliHLTDataTypeHistogram;
82
83 }
84
85 void AliHLTITSQAComponent::GetOutputDataSize( unsigned long& constBase, double& inputMultiplier )
86 {
87   // see header file for class documentation
88   // XXX TODO: Find more realistic values.
89   constBase = 80000;
90   inputMultiplier = 10;
91 }
92
93 AliHLTComponent* AliHLTITSQAComponent::Spawn()
94 {
95   // see header file for class documentation
96   return new AliHLTITSQAComponent;
97 }
98
99 int AliHLTITSQAComponent::DoInit( int argc, const char** argv )
100 {
101 // add AliITSQADataMakerRec constructor
102   fAliITSQADataMakerRec = new AliITSQADataMakerRec();
103   int iResult=0;
104   TString configuration="";
105   TString argument="";
106   for (int i=0; i<argc && iResult>=0; i++) {
107     argument=argv[i];
108     if (!configuration.IsNull()) configuration+=" ";
109     configuration+=argument;
110   }
111   
112   if (!configuration.IsNull()) {
113     iResult=Configure(configuration.Data());
114   }  
115
116   return iResult; 
117 }
118   
119 int AliHLTITSQAComponent::DoDeinit()
120 {
121   // see header file for class documentation
122 // add AliITSQADataMakerRec destruction
123   if(fAliITSQADataMakerRec) delete fAliITSQADataMakerRec;
124   return 0;
125 }
126
127 int AliHLTITSQAComponent::DoEvent(const AliHLTComponentEventData& /*evtData*/, AliHLTComponentTriggerData& /*trigData*/)
128 {
129   
130   int TotalSpacePoint = 0;
131   
132   const AliHLTComponentBlockData* iter = NULL;
133   
134   if(!IsDataEvent())
135     return 0;
136   
137   static Int_t rp = 1; 
138   // Check id histograms already created for this Event Specie
139   if ( rp ) { fAliITSQADataMakerRec->InitRecPoints(); rp = 0; }
140   for ( iter = GetFirstInputBlock(kAliHLTDataTypeClusters|kAliHLTDataOriginITSSSD); iter != NULL; iter = GetNextInputBlock() ) {
141     
142     const AliHLTITSClusterData* clusterData = (const AliHLTITSClusterData*) iter->fPtr;
143     Int_t nSpacepoint = (Int_t) clusterData->fSpacePointCnt;
144     TotalSpacePoint += nSpacepoint;
145     AliHLTITSSpacePointData *clusters = (AliHLTITSSpacePointData*) clusterData->fSpacePoints;
146
147     for(int i=0;i<nSpacepoint;i++){
148       Int_t lab[4]={0,0,0,0};
149       Float_t hit[6]={0,0,0,0,0,0};
150       Int_t info[3]={0,0,0};
151       
152       lab[0]=clusters[i].fTracks[0];
153       lab[1]=clusters[i].fTracks[1];
154       lab[2]=clusters[i].fTracks[2];
155       lab[3]=clusters[i].fIndex;
156       hit[0]=clusters[i].fY;
157       hit[1]=clusters[i].fZ;
158       hit[2]=clusters[i].fSigmaY2;
159       hit[3]=clusters[i].fSigmaZ2;
160       hit[4]=clusters[i].fQ;
161       hit[5]=clusters[i].fSigmaYZ;
162       info[0]=clusters[i].fNy;
163       info[1]=clusters[i].fNz;
164       info[2]=clusters[i].fLayer;
165
166       AliITSRecPoint recpoint(lab,hit,info);
167       fAliITSQADataMakerRec->FillRecPoint(recpoint);
168     }
169   }
170         for(Int_t i=0; i<6; i++) {
171                 //if(fPlotCharge){
172                 AliHLTUInt32_t specification = 0x0;
173                 PushBack( (TObject*) fAliITSQADataMakerRec->GetITSGlobalHisto(i),kAliHLTDataTypeHistogram,specification);
174                 //}
175         }
176         HLTInfo("ITSClusterHisto found %d Total Spacepoints", TotalSpacePoint);
177   
178   return 0;
179 }
180
181 int AliHLTITSQAComponent::Configure(const char* arguments)
182 {
183   
184   int iResult=0;
185   
186   if (!arguments) return iResult;
187   
188   TString allArgs=arguments;
189   TString argument;
190   
191   TObjArray* pTokens=allArgs.Tokenize(" ");
192   
193   if (pTokens) {
194     for (int i=0; i<pTokens->GetEntries() && iResult>=0; i++) {
195       argument=((TObjString*)pTokens->At(i))->GetString();
196       if (argument.IsNull()) continue;
197 /*      
198       if (argument.CompareTo("-plot-all")==0) {
199         HLTInfo("Ploting all historgams");
200         fPlotXY = kTRUE;
201         fPlotXYZ = kTRUE;
202         fPlotCharge = kTRUE;
203         continue;
204       }
205       else if (argument.CompareTo("-plot-xy")==0) {
206         HLTInfo("Ploting Global XY");
207         fPlotXY = kTRUE;
208         continue;
209       }
210
211       else if (argument.CompareTo("-plot-xyz")==0) {
212         HLTInfo("Ploting Global XYZ");
213         //fPlotXYZ = kTRUE;
214         continue;
215       }
216       else if (argument.CompareTo("-plot-charge")==0) {
217         HLTInfo("Ploting charge of clusters");
218         fPlotCharge = kTRUE;
219         continue;
220       }
221 */
222       else {
223         HLTError("unknown argument %s", argument.Data());
224         iResult=-EINVAL;
225         break;
226       }
227     }
228     delete pTokens;
229   }
230   
231   return iResult;
232 }
233
234 int AliHLTITSQAComponent::Reconfigure(const char* cdbEntry, const char* chainId)
235 {
236   // see header file for class documentation
237   int iResult=0;
238   
239   const char* path="HLT/ConfigITS/HistoComponent";
240   const char* defaultNotify="";
241   if (cdbEntry) {
242     path=cdbEntry;
243     defaultNotify=" (default)";
244   }
245   if (path) {
246     HLTInfo("reconfigure from entry %s%s, chain id %s", path, defaultNotify,(chainId!=NULL && chainId[0]!=0)?chainId:"<none>");
247     AliCDBEntry *pEntry = AliCDBManager::Instance()->Get(path);
248     if (pEntry) {
249       TObjString* pString=dynamic_cast<TObjString*>(pEntry->GetObject());
250       if (pString) {
251         HLTInfo("received configuration object string: \'%s\'", pString->GetString().Data());
252         iResult=Configure(pString->GetString().Data());
253       } else {
254         HLTError("configuration object \"%s\" has wrong type, required TObjString", path);
255       }
256     } else {
257       HLTError("can not fetch object \"%s\" from CDB", path);
258     }
259   }
260   
261   return iResult;
262 }