2 MCH DA for online occupancy
4 Contact: Laurent Aphecetche <laurent.aphecetche@subatech.in2p3.fr>, Jean-Luc Charvet <jean-luc.charvet@cea.fr>, Alberto Baldisseri <alberto.baldisseri@cea.fr>
6 Run Type: PHYSICS STANDALONE
8 Number of events needed: all (or at least as much as possible...)
9 Input Files: 09000094301009.10.raw
10 Output Files: mch.occupancy, to be exported to the DAQ FXS
11 Trigger types used: PHYSICS_EVENT
14 /**************************************************************************
15 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
17 * Author: The ALICE Off-line Project. *
18 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
30 /// MUON TRACKER DA to compute the hit count at manu level.
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)
39 #include "AliMUON2DMap.h"
40 #include "AliMUONCalibParamNI.h"
41 #include "AliMUONRawStreamTrackerHP.h"
42 #include "AliMpConstants.h"
43 #include "AliRawEventHeaderBase.h"
44 #include "AliRawReaderDate.h"
45 #include "Riostream.h"
46 #include "TPluginManager.h"
49 #include "TTimeStamp.h"
50 #include "TStopwatch.h"
57 #include "TObjString.h"
62 const char* OUTPUT_FILE = "mch.occupancy";
63 const char* DAVERSION = "MUONTRKOCCda v1.2 ($Id$)";
65 //______________________________________________________________________________
66 void Add(AliMUONVStore& destStore, const AliMUONVStore& srcStore)
68 /// Add all elements from srcStore to destStore
69 /// Each element of srcStore is supposed to be an AliMUONCalibParamNI,
70 /// with ID0=busPatchId and ID1=manuId
72 TIter next(srcStore.CreateIterator());
73 AliMUONVCalibParam* source;
75 while ( ( source = static_cast<AliMUONVCalibParam*>(next()) ) )
77 AliMUONCalibParamNI* dest = static_cast<AliMUONCalibParamNI*>(destStore.FindObject(source->ID0(),source->ID1()));
80 dest = static_cast<AliMUONCalibParamNI*>(source->Clone());
85 for ( Int_t i = 0; i < source->Size(); ++i )
87 for ( Int_t j = 0; j < source->Dimension(); ++j )
89 dest->SetValueAsIntFast(i,j,dest->ValueAsIntFast(i,j)+source->ValueAsIntFast(i,j));
96 //______________________________________________________________________________
97 void GenerateOutputFile(const AliMUONVStore& store, ostream& out,
98 Int_t runNumber, Int_t nevents)
100 /// Write the channel hit count (grouped by manu) in the output file.
102 TIter next(store.CreateIterator());
103 AliMUONVCalibParam* manu;
105 out << "//===========================================================================" << endl;
106 out << "// Hit counter file calculated by " << __FILE__ << endl;
107 out << "//===========================================================================" << endl;
109 out << "// * Run Number : " << runNumber << endl;
110 out << "// * File Creation Date : " << TTimeStamp().AsString("l") << endl;
111 out << "//---------------------------------------------------------------------------" << endl;
112 out << "// BP MANU SUM_N NEVENTS" << endl;
113 out << "//---------------------------------------------------------------------------" << endl;
115 while ( ( manu = static_cast<AliMUONVCalibParam*>(next()) ) )
120 for ( Int_t i = 0; i < manu->Size(); ++i )
122 sum += manu->ValueAsInt(i);
123 // nevents = TMath::Max(nevents,manu->ValueAsInt(i,1));
124 // nevents = TMath::Max(nevents,manu->ValueAsInt(i,1));
127 out << Form("%5d %5d %10d %10d",manu->ID0(),manu->ID1(),sum,nevents) << endl;
131 //______________________________________________________________________________
132 int main(int argc, char **argv)
136 /// We loop over all physics events.
137 /// For each event we store the channels that were hit for that event.
138 /// If the event is good, we then increment the list of channels hit for
140 /// We delete the store for a single event and move to next event.
142 /// In the end we output an ASCII file with the necessary information
143 /// to compute the occupancy later on, i.e. the number of times channels
144 /// were seen per manu, and the number of events.
150 ios::sync_with_stdio();
152 cout << "Running " << DAVERSION << endl;
156 cout << "Wrong number of arguments" << endl;
157 cout << "Usage : " << argv[0] << " datasource1 [datasource2] ..." << endl;
161 // needed for streamer application
162 gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
168 Int_t numberOfEvents(0);
169 Int_t numberOfPhysicsEvent(0);
170 Int_t numberOfBadEvents(0);
171 Int_t numberOfUsedEvents(0);
173 AliMUON2DMap oneEventData(kTRUE);
174 AliMUON2DMap accumulatedData(kTRUE);
178 for ( Int_t i = 1; i < argc; ++i )
181 AliRawReaderDate* rawReader(0x0);
183 // define data source :
184 status=monitorSetDataSource(argv[i]);
187 printf("MCH Occupancy DA ERROR: monitorSetDataSource() failed: %s\n", monitorDecodeError(status));
191 // Declare monitoring program
192 status=monitorDeclareMp("MUON_TRK_OCC");
195 printf("MCH Occupancy DA ERROR: monitorDeclareMp() failed: %s\n", monitorDecodeError(status));
198 // Define wait event timeout - 1s max
200 monitorSetNoWaitNetworkTimeout(1000);
204 struct eventHeaderStruct *event;
205 eventTypeType eventT;
207 status=monitorGetEventDynamic((void **)&event);
210 printf("MCH Occupancy DA ERROR: %s\n", monitorDecodeError(status));
215 /* check shutdown condition */
216 if (daqDA_checkShutdown())
222 /* retry if got no event */
223 if (event==NULL) continue;
227 eventT=event->eventType;
228 if ((eventT == END_OF_RUN)||(eventT == END_OF_RUN_FILES))
233 if (eventT != PHYSICS_EVENT)
239 ++numberOfPhysicsEvent;
241 rawReader = new AliRawReaderDate((void*)event);
243 if ( rawReader->GetRunNumber() != runNumber )
245 if ( runNumber != 0 )
247 cout << "Uh oh. That's bad... Changing of run number ???" << endl;
252 runNumber = rawReader->GetRunNumber();
255 AliMUONRawStreamTrackerHP stream(rawReader);
257 stream.DisableWarnings();
259 oneEventData.Clear();
268 while ( stream.Next(buspatchId,manuId,manuChannel,adc,kTRUE) )
270 AliMUONVCalibParam* one = static_cast<AliMUONVCalibParam*>(oneEventData.FindObject(buspatchId,manuId));
274 one = new AliMUONCalibParamNI(1,AliMpConstants::ManuNofChannels(),buspatchId,manuId);
275 oneEventData.Add(one);
278 one->SetValueAsInt(manuChannel,0,one->ValueAsInt(manuChannel,0)+1);
281 Bool_t badEvent = stream.HasPaddingError() || stream.HasGlitchError();
285 ++numberOfUsedEvents;
286 Add(accumulatedData,oneEventData);
300 cout << Form("%12d events processed : %12d physics %d used ones %d bad ones",
301 numberOfEvents,numberOfPhysicsEvent,numberOfUsedEvents,numberOfBadEvents) << endl;
303 ofstream fout(OUTPUT_FILE);
305 GenerateOutputFile(accumulatedData,fout,runNumber,numberOfUsedEvents);
311 // Send occupancy store (as a big string) to the AMORE DB
312 amore::da::AmoreDA amoreDA(amore::da::AmoreDA::kSender);
316 GenerateOutputFile(accumulatedData,str,runNumber,numberOfUsedEvents);
318 TObjString occupancyAsString(str.str().c_str());
320 Int_t status = amoreDA.Send("Occupancy",&occupancyAsString);
323 cerr << "ERROR : Failed to write occupancies in the AMORE database : " << status << endl;
328 /* store the result file on FXS */
329 if (daqDA_FES_storeFile(OUTPUT_FILE,"OCCUPANCY")) return -9;
332 printf("\nExecution time : R:%7.2fs C:%7.2fs\n", timers.RealTime(), timers.CpuTime());