Bumping the version number
[u/mrichter/AliRoot.git] / MUON / MUONTRKOCCda.cxx
CommitLineData
2825ba79 1/*
2 MCH DA for online occupancy
3
4 Contact: Laurent Aphecetche <laurent.aphecetche@subatech.in2p3.fr>, Jean-Luc Charvet <jean-luc.charvet@cea.fr>, Alberto Baldisseri <alberto.baldisseri@cea.fr>
5 Link:
6 Run Type: PHYSICS STANDALONE
7 DA Type: MON
8 Number of events needed: all (or at least as much as possible...)
be9b25ed 9 Input Files: 09000094301009.10.raw
2825ba79 10 Output Files: mch.occupancy, to be exported to the DAQ FXS
11 Trigger types used: PHYSICS_EVENT
12*/
13
14/**************************************************************************
15 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
16 * *
17 * Author: The ALICE Off-line Project. *
18 * Contributors are mentioned in the code where appropriate. *
19 * *
20 * Permission to use, copy, modify and distribute this software and its *
21 * documentation strictly for non-commercial purposes is hereby granted *
22 * without fee, provided that the above copyright notice appears in all *
23 * copies and that both the copyright notice and this permission notice *
24 * appear in the supporting documentation. The authors make no claims *
25 * about the suitability of this software for any purpose. It is *
26 * provided "as is" without express or implied warranty. *
27 **************************************************************************/
28
29///
30/// MUON TRACKER DA to compute the hit count at manu level.
31///
32/// In the end, this DA produces an ASCII file containing
33/// the hit count of all the manus (that were seen in the data flow)
34/// of the MUON Tracker (and the number of seen events, so we can
35/// later on compute the occupancy)
36///
37/// $Id$
38
39#include "AliMUON2DMap.h"
40#include "AliMUONCalibParamNI.h"
41#include "AliMUONRawStreamTrackerHP.h"
42#include "AliMpConstants.h"
43#include "AliRawEventHeaderBase.h"
134d9649 44#include "AliRawReaderDate.h"
2825ba79 45#include "Riostream.h"
46#include "TPluginManager.h"
47#include "TROOT.h"
48#include "TString.h"
49#include "TTimeStamp.h"
50#include "TStopwatch.h"
51#include "daqDA.h"
52#include "event.h"
134d9649 53#include "monitor.h"
65cd97ed 54#include "signal.h"
2825ba79 55
56#ifdef ALI_AMORE
57#include <AmoreDA.h>
5431405e 58#include "TObjString.h"
59#include "TSystem.h"
60#include <sstream>
2825ba79 61#endif
62
63const char* OUTPUT_FILE = "mch.occupancy";
29e40e09 64const char* DAVERSION = "MUONTRKOCCda v1.5 ($Id$)";
2825ba79 65
66//______________________________________________________________________________
67void Add(AliMUONVStore& destStore, const AliMUONVStore& srcStore)
68{
69 /// Add all elements from srcStore to destStore
70 /// Each element of srcStore is supposed to be an AliMUONCalibParamNI,
71 /// with ID0=busPatchId and ID1=manuId
72
73 TIter next(srcStore.CreateIterator());
74 AliMUONVCalibParam* source;
75
76 while ( ( source = static_cast<AliMUONVCalibParam*>(next()) ) )
77 {
78 AliMUONCalibParamNI* dest = static_cast<AliMUONCalibParamNI*>(destStore.FindObject(source->ID0(),source->ID1()));
79 if (!dest)
80 {
81 dest = static_cast<AliMUONCalibParamNI*>(source->Clone());
82 destStore.Add(dest);
83 }
84 else
85 {
86 for ( Int_t i = 0; i < source->Size(); ++i )
87 {
88 for ( Int_t j = 0; j < source->Dimension(); ++j )
89 {
be9b25ed 90 dest->SetValueAsIntFast(i,j,dest->ValueAsIntFast(i,j)+source->ValueAsIntFast(i,j));
2825ba79 91 }
92 }
93 }
94 }
95}
96
97//______________________________________________________________________________
98void GenerateOutputFile(const AliMUONVStore& store, ostream& out,
99 Int_t runNumber, Int_t nevents)
100{
101 /// Write the channel hit count (grouped by manu) in the output file.
102
103 TIter next(store.CreateIterator());
104 AliMUONVCalibParam* manu;
105
106 out << "//===========================================================================" << endl;
107 out << "// Hit counter file calculated by " << __FILE__ << endl;
108 out << "//===========================================================================" << endl;
109 out << "//" << endl;
110 out << "// * Run Number : " << runNumber << endl;
111 out << "// * File Creation Date : " << TTimeStamp().AsString("l") << endl;
112 out << "//---------------------------------------------------------------------------" << endl;
113 out << "// BP MANU SUM_N NEVENTS" << endl;
114 out << "//---------------------------------------------------------------------------" << endl;
115
116 while ( ( manu = static_cast<AliMUONVCalibParam*>(next()) ) )
117 {
118 Int_t sum(0);
119// Int_t nevents(0);
120
121 for ( Int_t i = 0; i < manu->Size(); ++i )
122 {
123 sum += manu->ValueAsInt(i);
124 // nevents = TMath::Max(nevents,manu->ValueAsInt(i,1));
125 // nevents = TMath::Max(nevents,manu->ValueAsInt(i,1));
126 }
127
128 out << Form("%5d %5d %10d %10d",manu->ID0(),manu->ID1(),sum,nevents) << endl;
129 }
130}
131
132//______________________________________________________________________________
133int main(int argc, char **argv)
134{
135 /// Main method.
136 ///
137 /// We loop over all physics events.
138 /// For each event we store the channels that were hit for that event.
139 /// If the event is good, we then increment the list of channels hit for
140 /// the whole run.
141 /// We delete the store for a single event and move to next event.
142 ///
143 /// In the end we output an ASCII file with the necessary information
144 /// to compute the occupancy later on, i.e. the number of times channels
145 /// were seen per manu, and the number of events.
146 ///
147
a9d45f70 148 signal(SIGSEGV,SIG_DFL); // to be able to get core dumps...
149
2825ba79 150 TStopwatch timers;
151 timers.Start(kTRUE);
152
153 ios::sync_with_stdio();
154
155 cout << "Running " << DAVERSION << endl;
156
157 if ( argc < 2 )
158 {
159 cout << "Wrong number of arguments" << endl;
160 cout << "Usage : " << argv[0] << " datasource1 [datasource2] ..." << endl;
161 return -1;
162 }
163
164 // needed for streamer application
165 gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
166 "*",
167 "TStreamerInfo",
168 "RIO",
169 "TStreamerInfo()");
170
171 Int_t numberOfEvents(0);
172 Int_t numberOfPhysicsEvent(0);
173 Int_t numberOfBadEvents(0);
174 Int_t numberOfUsedEvents(0);
175
176 AliMUON2DMap oneEventData(kTRUE);
177 AliMUON2DMap accumulatedData(kTRUE);
178
179 UInt_t runNumber(0);
180
181 for ( Int_t i = 1; i < argc; ++i )
182 {
134d9649 183 int status;
184 AliRawReaderDate* rawReader(0x0);
185
186// define data source :
033da0c5 187 status=monitorSetDataSource(argv[i]);
134d9649 188 if (status!=0)
189 {
190 printf("MCH Occupancy DA ERROR: monitorSetDataSource() failed: %s\n", monitorDecodeError(status));
191 return -1;
192 }
193
194// Declare monitoring program
195 status=monitorDeclareMp("MUON_TRK_OCC");
196 if (status!=0)
197 {
198 printf("MCH Occupancy DA ERROR: monitorDeclareMp() failed: %s\n", monitorDecodeError(status));
199 return -1;
200 }
201// Define wait event timeout - 1s max
202 monitorSetNowait();
203 monitorSetNoWaitNetworkTimeout(1000);
204
205 for(;;)
2825ba79 206 {
63108d64 207 struct eventHeaderStruct *event(0x0);
134d9649 208 eventTypeType eventT;
209
210 status=monitorGetEventDynamic((void **)&event);
211 if (status!=0)
212 {
be9b25ed 213 printf("MCH Occupancy DA ERROR: %s\n", monitorDecodeError(status));
214 delete event;
215 break;
134d9649 216 }
217
2825ba79 218 /* check shutdown condition */
be9b25ed 219 if (daqDA_checkShutdown())
2825ba79 220 {
be9b25ed 221 delete event;
222 break;
2825ba79 223 }
be9b25ed 224
134d9649 225 /* retry if got no event */
226 if (event==NULL) continue;
227
228 ++numberOfEvents;
229
230 eventT=event->eventType;
be9b25ed 231 if ((eventT == END_OF_RUN)||(eventT == END_OF_RUN_FILES))
232 {
233 delete event;
234 break;
235 }
236 if (eventT != PHYSICS_EVENT)
237 {
238 delete event;
239 continue;
240 }
134d9649 241
2825ba79 242 ++numberOfPhysicsEvent;
243
134d9649 244 rawReader = new AliRawReaderDate((void*)event);
245
2825ba79 246 if ( rawReader->GetRunNumber() != runNumber )
247 {
248 if ( runNumber != 0 )
249 {
250 cout << "Uh oh. That's bad... Changing of run number ???" << endl;
be9b25ed 251 delete event;
252 delete rawReader;
2825ba79 253 return -9999;
254 }
255 runNumber = rawReader->GetRunNumber();
256 }
257
258 AliMUONRawStreamTrackerHP stream(rawReader);
259
260 stream.DisableWarnings();
2825ba79 261
262 oneEventData.Clear();
263
264 Int_t buspatchId;
265 UShort_t manuId;
266 UChar_t manuChannel;
267 UShort_t adc;
268
269 stream.First();
134d9649 270
2825ba79 271 while ( stream.Next(buspatchId,manuId,manuChannel,adc,kTRUE) )
272 {
273 AliMUONVCalibParam* one = static_cast<AliMUONVCalibParam*>(oneEventData.FindObject(buspatchId,manuId));
274
275 if (!one)
276 {
277 one = new AliMUONCalibParamNI(1,AliMpConstants::ManuNofChannels(),buspatchId,manuId);
278 oneEventData.Add(one);
279 }
280
281 one->SetValueAsInt(manuChannel,0,one->ValueAsInt(manuChannel,0)+1);
282 }
283
284 Bool_t badEvent = stream.HasPaddingError() || stream.HasGlitchError();
285
286 if ( !badEvent )
287 {
288 ++numberOfUsedEvents;
289 Add(accumulatedData,oneEventData);
290 }
291 else
292 {
293 ++numberOfBadEvents;
294 }
be9b25ed 295
be9b25ed 296 delete event;
297
298 delete rawReader;
299 }
2825ba79 300 }
301
302
303 cout << Form("%12d events processed : %12d physics %d used ones %d bad ones",
304 numberOfEvents,numberOfPhysicsEvent,numberOfUsedEvents,numberOfBadEvents) << endl;
305
306 ofstream fout(OUTPUT_FILE);
307
308 GenerateOutputFile(accumulatedData,fout,runNumber,numberOfUsedEvents);
309
310 fout.close();
311
312#ifdef ALI_AMORE
313
4b50c4a8 314 if ( numberOfUsedEvents ) // do not update the AMORE pool with an empty object...
2825ba79 315 {
4b50c4a8 316 // Send occupancy store (as a big string) to the AMORE DB
317 amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender);
318
319 ostringstream str;
320
321 GenerateOutputFile(accumulatedData,str,runNumber,numberOfUsedEvents);
322
323 TObjString occupancyAsString(str.str().c_str());
324
325 Int_t status = amoreDA.Send("Occupancy",&occupancyAsString);
326 if ( status )
327 {
328 cerr << "ERROR : Failed to write occupancies in the AMORE database : " << status << endl;
329 }
330 }
2825ba79 331
332#endif
333
334 /* store the result file on FXS */
335 if (daqDA_FES_storeFile(OUTPUT_FILE,"OCCUPANCY")) return -9;
336
337 timers.Stop();
338 printf("\nExecution time : R:%7.2fs C:%7.2fs\n", timers.RealTime(), timers.CpuTime());
339
340 return 0;
341}