1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
19 /// \file MUONOfflineShift.C
20 /// \brief Macro to be used to check raw data during MUON offline shifts.
24 /// MUONOfflineShift("path_to_raw_file","basename of output file"); > log
26 /// (the redirection to an output log file is recommended as the output from
27 /// this macro might be quite long...)
29 /// This will read the raw data and process it several times, varying what's done :
30 /// only decoding, decoding + zero-suppression, etc... (TBE)
32 /// Two outputs files will be created :
34 /// - basename.root, containing AliMUONVTrackerData objects that can then
35 /// be displayed using the mchview program (using mchview --use basename.root)
37 /// - basename.log, containing (for the moment) the occupancy numbers of the various detection elements
39 /// \author Laurent Aphecetche
41 #if !defined(__CINT__) || defined(__MAKECINT__)
43 #include "AliCDBManager.h"
44 #include "AliCodeTimer.h"
45 #include "AliMUONPainterRegistry.h"
46 #include "AliMUONTrackerCalibratedDataMaker.h"
47 #include "AliMUONTrackerRawDataMaker.h"
48 #include "AliMUONVTrackerData.h"
50 #include "AliMpConstants.h"
51 #include "AliMpDEIterator.h"
52 #include "AliRawReader.h"
53 #include "AliRawReaderDate.h"
54 #include "AliRawReaderRoot.h"
55 #include <Riostream.h>
58 #include <TStopwatch.h>
64 //______________________________________________________________________________
65 Int_t DataMakerReading(const char* input,
67 const char* cdbPath="",
68 const char* calibMode="",
69 Bool_t histogram=kFALSE,
71 Double_t xmax = 4096.0)
73 /// Run over the data and calibrate it if so required (if cdbPath != "")
74 /// calibMode can be :
75 /// - NOGAIN : only zero-suppression will be done
76 /// - GAINCONSTANTCAPA : zero-suppression + gain, but with a single capa value for all channels
77 /// - GAIN : zero-suppression + gain w/ individual capacitance per channel.
79 TString fileName(gSystem->ExpandPathName(input));
81 AliRawReader* rawReader(0x0);
83 // check extention to choose the rawdata file format
84 if (fileName.EndsWith(".root"))
86 rawReader = new AliRawReaderRoot(fileName);
88 else if (!fileName.IsNull())
90 rawReader = new AliRawReaderDate(fileName); // DATE file
93 if (!rawReader) return 0;
95 AliMUONVTrackerDataMaker* dm(0x0);
97 if ( strlen(cdbPath) > 0 )
99 dm = new AliMUONTrackerCalibratedDataMaker(rawReader,cdbPath,calibMode,histogram,xmin,xmax);
103 dm = new AliMUONTrackerRawDataMaker(rawReader,kTRUE);
106 AliMUONPainterRegistry::Instance()->Register(dm);
111 dm->SetRunning(kTRUE);
113 while (dm->NextEvent())
123 //______________________________________________________________________________
124 void Print(const char* method, TStopwatch& timer, Int_t n)
126 /// Get the timing for a given method
128 cout << Form("%20s %10d events. Total CPU time %7.2f seconds.",
129 method,n,timer.CpuTime());
131 Double_t cpu = timer.CpuTime()/n;
132 Double_t real = timer.RealTime()/n;
134 cout << Form(" ms real/event = %7.2f ms CPU/event = %7.2f",real*1E3,cpu*1E3)
138 //______________________________________________________________________________
139 void Occupancy(ostream& outfile)
141 /// Write occupancy numbers to output text file
143 outfile << "-----------------------------------------------------" << endl;
144 outfile << "Occupancy numbers" << endl;
145 outfile << "-----------------------------------------------------" << endl;
147 const Int_t occIndex = 2;
149 AliMUONPainterRegistry* reg = AliMUONPainterRegistry::Instance();
151 Int_t nofDataSources = reg->NumberOfDataSources();
153 outfile << Form("%11s|"," ");
155 for ( Int_t ids = 0; ids < nofDataSources; ++ ids )
157 AliMUONVTrackerData* data = reg->DataSource(ids);
158 outfile << Form(" %13s |",data->GetName());
163 for ( Int_t chamberId = 0; chamberId < AliMpConstants::NofTrackingChambers(); ++chamberId )
165 Bool_t nonZero(kFALSE);
166 for ( Int_t ids = 0; ids < nofDataSources && nonZero == kFALSE; ++ ids )
168 if ( reg->DataSource(ids)->Chamber(chamberId,occIndex) ) nonZero = kTRUE;
171 if ( !nonZero ) continue;
173 outfile << Form("Chamber %2d |",chamberId);
174 for ( Int_t ids = 0; ids < nofDataSources; ++ ids )
176 AliMUONVTrackerData* data = reg->DataSource(ids);
177 outfile << Form(" %7.2f %% |",100.0*data->Chamber(chamberId,occIndex));
185 Int_t detElemId = it.CurrentDEId();
186 Bool_t nonZero(kFALSE);
187 for ( Int_t ids = 0; ids < nofDataSources && nonZero == kFALSE; ++ ids )
189 AliMUONVTrackerData* data = reg->DataSource(ids);
190 if ( data->DetectionElement(detElemId,occIndex) > 0 )
198 outfile << Form(" DE %04d |",detElemId);
199 for ( Int_t ids = 0; ids < nofDataSources; ++ ids )
201 AliMUONVTrackerData* data = reg->DataSource(ids);
202 outfile << Form(" %7.2f %% |",100.0*data->DetectionElement(detElemId,occIndex));
211 //______________________________________________________________________________
212 void MUONOfflineShift(const char* input="alien:///alice/data/2008/LHC08a/000021931/raw/08000021931001.50.root",
213 const char* outputBase="21931.001.50",
214 const char* ocdbPath="alien://folder=/alice/data/2008/LHC08a/OCDB")
216 /// Entry point of the macro.
217 /// Example of syntax for an input file (from alien)
219 /// alien::///alice/data/2007/LHC07w/000014493/raw/07000014493001.10.root
221 /// and for an OCDB path :
223 /// alien://folder=/alice/data/2007/LHC07w/OCDB
226 TGrid::Connect("alien://");
228 AliCodeTimer::Instance()->Reset();
230 AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
231 AliCDBManager::Instance()->SetRun(0);
232 AliMpCDB::LoadDDLStore();
239 Int_t n1 = DataMakerReading(input,timer1);
241 Int_t n2 = DataMakerReading(input,timer2,ocdbPath,"NOGAIN",kTRUE);
243 Int_t n3 = DataMakerReading(input,timer3,ocdbPath,"GAINCONSTANTCAPA",kTRUE);
245 Int_t n4 = DataMakerReading(input,timer4,ocdbPath,"GAIN",kTRUE);
247 Print("DataMakerReading(HRAW)",timer1,n1);
248 Print("DataMakerReading(HCALZ)",timer2,n2);
249 Print("DataMakerReading(HCALG)",timer3,n3);
250 Print("DataMakerReading(HCALC)",timer4,n4);
252 AliMUONPainterRegistry* reg = AliMUONPainterRegistry::Instance();
254 TFile f(gSystem->ExpandPathName(Form("%s.root",outputBase)),"RECREATE");
255 ofstream out(gSystem->ExpandPathName(Form("%s.log",outputBase)));
259 for ( Int_t i = 0; i < reg->NumberOfDataSources(); ++i )
261 AliMUONVTrackerData* data = reg->DataSource(i);