/*
+MTR DA for online
-Version 1 for MUONTRGda MUON trigger
-Working version for reading back raw data
-(Ch. Finck)
-
+Contact: Franck Manso <manso@clermont.in2p3.fr>
+Link: http://aliceinfo.cern.ch/static/Offline/dimuon/muon_html/README_mtrda.html
+Reference run: 61898
+Run Type: PHYSICS
+DA Type: MON
+Number of events needed: 1000 events
+Input files: MtgGlobalCrate.dat MtgRegionalCrate.dat MtgLocalMask.dat MtgLocalLut.dat MtgCurrent.dat DAConfig.dat
+Output Files: ExportedFiles.dat MtgGlobalCrate.dat
+Trigger types used: PHYSICS_EVENT CALIBRATION_EVENT
*/
+
+//////////////////////////////////////////////////////////////////////////////
+// Detector Algorithm for the MUON trigger configuration. //
+// //
+// Calculates masks for the global trigger input, by looking at dead //
+// channels in calibration events and at noisy channels in physics events. //
+// Transfers trigger configuration files to the File Exchange Server and //
+// writes them (if modified) into the trigger configuration data base. //
+// //
+// Authors: //
+// Christian Fink (formerly at Subatech, Nantes) //
+// Franck Manso (LPC Clermont-Ferrand, manso@clermont.in2p3.fr) //
+// Bogdan Vulpescu (LPC Clermont-Ferrand, vulpescu@clermont.in2p3.fr) //
+// //
+//////////////////////////////////////////////////////////////////////////////
+
+//#define OFFLINE
+
extern "C" {
#include <daqDA.h>
}
#include <stdio.h>
#include <stdlib.h>
-//AliRoot
+#include "AliRawReaderDate.h"
+
+#include "AliMpConstants.h"
#include "AliMUONRawStreamTrigger.h"
+#include "AliMUONRawStreamTriggerHP.h"
#include "AliMUONDarcHeader.h"
#include "AliMUONRegHeader.h"
-#include "AliMUONLocalStruct.h"
#include "AliMUONDDLTrigger.h"
-// #include "AliMUONVStore.h"
-// #include "AliMUON2DMap.h"
-// #include "AliMUONCalibParamNF.h"
-// #include "AliMpDDLStore.h"
-// #include "AliMpIntPair.h"
-#include "AliMpConstants.h"
-#include "AliRawReaderDate.h"
-
-
-//ROOT
+#include "AliMUONVStore.h"
+#include "AliMUON1DArray.h"
+#include "AliMUONTriggerIO.h"
+#include "AliMUONRegionalTriggerConfig.h"
+#include "AliMUONGlobalCrateConfig.h"
+#include "AliMUONVCalibParam.h"
+#include "AliMUONCalibParamND.h"
+#include "AliMUONCalibParamNI.h"
+#include "AliMUONLocalStruct.h"
#include "TString.h"
+#include "TSystem.h"
#include "TStopwatch.h"
-#include "TMath.h"
-#include "TTimeStamp.h"
#include "TROOT.h"
#include "TPluginManager.h"
+#include "TArrayS.h"
+/// class for DA run parameters and DA working space
+class AliDAConfig : TObject {
-// global variables
-TString command("pat");
-UInt_t runNumber = 0;
-Int_t nEvents = 0;
+public:
+ AliDAConfig() :
+ fDAConfigFileName("DAConfig.dat"),
+ fCurrentFileName("MtgCurrent.dat"),
+ fLastCurrentFileName("MtgLastCurrent.dat"),
+ fSodName(""),
+ fSodFlag(0),
+ fDAName(""),
+ fDAFlag(0),
+ fDAMode(1),
+ fGlobalFileName(""),
+ fRegionalFileName(""),
+ fLocalMaskFileName(""),
+ fLocalLutFileName(""),
+ fSignatureFileName(""),
+ fGlobalFileVersion(0),
+ fRegionalFileVersion(0),
+ fLocalMaskFileVersion(0),
+ fLocalLutFileVersion(0),
+ fSignatureFileVersion(0),
+ fGlobalFileLastVersion(0),
+ fRegionalFileLastVersion(0),
+ fLocalMaskFileLastVersion(0),
+ fLocalLutFileLastVersion(0),
+ fEventsN(0),
+ fEventsD(0),
+ fPrintLevel(0),
+ fLocalMasks(0x0),
+ fLocalMasksDA(0x0),
+ fRegionalMasks(0x0),
+ fGlobalMasks(0x0),
+ fTriggerIO(new AliMUONTriggerIO),
+ fAlgoNoisyInput(kFALSE),
+ fAlgoDeadcInput(kFALSE),
+ fThrN(0.1),
+ fThrD(0.9),
+ fThrLocN(0.1),
+ fThrLocD(0.9),
+ fMinEvents(10),
+ fSkipEvents(0),
+ fMaxEvents(65535),
+ fWithWarnings(kFALSE),
+ fUseFastDecoder(kFALSE),
+ fNLocalBoard(AliMpConstants::TotalNofLocalBoards()+1),
+ fPatternStoreN(0x0),
+ fPatternStoreD(0x0)
+ {
+ /// default constructor
+ for (Int_t ii = 0; ii < kGlobalInputs; ii++) {
+ for (Int_t il = 0; il < kGlobalInputLength; il++) {
+ fAccGlobalInputN[ii][il] = 0;
+ fAccGlobalInputD[ii][il] = 0;
+ }
+ }
+ fLocalMasks = new AliMUON1DArray(fNLocalBoard);
+ fLocalMasksDA = new AliMUON1DArray(fNLocalBoard);
+ fRegionalMasks = new AliMUONRegionalTriggerConfig();
+ fGlobalMasks = new AliMUONGlobalCrateConfig();
+ fPatternStoreN = new AliMUON1DArray(fNLocalBoard);
+ fPatternStoreD = new AliMUON1DArray(fNLocalBoard);
-//*************************************************************//
+ // Generate local trigger masks store. All masks are set to FFFF
+ for (Int_t i = 1; i <= AliMpConstants::TotalNofLocalBoards(); i++) {
+ AliMUONVCalibParam* localBoard = new AliMUONCalibParamNI(1,8,i,0,0);
+ for (Int_t x = 0; x < 2; x++) {
+ for (Int_t y = 0; y < 4; y++) {
+ Int_t index = x*4+y;
+ localBoard->SetValueAsInt(index,0,0xFFFF);
+ }
+ }
+ fLocalMasksDA->Add(localBoard);
+ }
+
+ // Create a default pattern store
+ for (Int_t i = 1; i <= AliMpConstants::TotalNofLocalBoards(); i++) {
+ AliMUONVCalibParam *patN = new AliMUONCalibParamND(2, 64, i, 0, 0.);
+ AliMUONVCalibParam *patD = new AliMUONCalibParamND(2, 64, i, 0, 0.);
+ fPatternStoreN->Add(patN);
+ fPatternStoreD->Add(patD);
+ }
-// main routine
-int main(Int_t argc, Char_t **argv)
+ }
+
+ virtual ~AliDAConfig()
+ {
+ /// destructor
+ delete fLocalMasks;
+ delete fLocalMasksDA;
+ delete fRegionalMasks;
+ delete fGlobalMasks;
+ delete fPatternStoreN;
+ delete fPatternStoreD;
+ delete fTriggerIO;
+ }
+ void PrintConfig()
+ {
+ /// print DA parameters
+ printf("DA config file name: %s \n",GetDAConfigFileName());
+ printf("Current file name: %s \n",GetCurrentFileName());
+ printf("Last current file name: %s \n",GetLastCurrentFileName());
+ printf("Global file name: %s (%d %d)\n",GetGlobalFileName(),GetGlobalFileVersion(),GetGlobalFileLastVersion());
+ printf("Regional file name: %s (%d %d)\n",GetRegionalFileName(),GetRegionalFileVersion(),GetRegionalFileLastVersion());
+ printf("Local mask file name: %s (%d %d)\n",GetLocalMaskFileName(),GetLocalMaskFileVersion(),GetLocalMaskFileLastVersion());
+ printf("Local LUT file name: %s (%d %d)\n",GetLocalLutFileName(),GetLocalLutFileVersion(),GetLocalLutFileLastVersion());
+ printf("Signature file name: %s (%d)\n",GetSignatureFileName(),GetSignatureFileVersion());
+ }
+
+ /// name of the DA configuration file from detector DB
+ const Char_t* GetDAConfigFileName() { return fDAConfigFileName.Data(); }
+ /// file with current versions of the configuration files, usually MtgCurrent.dat
+ const Char_t* GetCurrentFileName() { return fCurrentFileName.Data(); }
+ /// last known versions of the configuration files, usually MtgLastCurrent.dat
+ const Char_t* GetLastCurrentFileName() { return fLastCurrentFileName.Data(); }
+
+ /// name of the Start-of-data field in MtgCurrent.dat
+ const Char_t* GetSodName() { return fSodName.Data(); }
+ /// flag value of the Start-of-data field in MtgCurrent.dat
+ Int_t GetSodFlag() const { return fSodFlag; }
+
+ /// name of the Detector Algorithm field in MtgCurrent.dat
+ const Char_t* GetDAName() { return fDAName.Data(); }
+ /// flag value of the Detector Algorithm field in MtgCurrent.dat
+ Int_t GetDAFlag() const { return fDAFlag; }
+ /// DA active mode (GLOBAL or GLOBAL+LOCAL)
+ Int_t GetDAMode() const { return fDAMode; }
+
+ /// global crate configuration file name
+ const Char_t* GetGlobalFileName() { return fGlobalFileName.Data(); }
+ /// regional crate configuration file name
+ const Char_t* GetRegionalFileName() { return fRegionalFileName.Data(); }
+ /// masks for the local cards, file name
+ const Char_t* GetLocalMaskFileName() { return fLocalMaskFileName.Data(); }
+ /// transverse momentum Look-Up-Table, file name
+ const Char_t* GetLocalLutFileName() { return fLocalLutFileName.Data(); }
+ /// signature file name
+ const Char_t* GetSignatureFileName() { return fSignatureFileName.Data(); }
+
+ /// version of the global crate configuration in the detector DB
+ Int_t GetGlobalFileVersion() const { return fGlobalFileVersion; }
+ /// version of the regional crate configuration in the detector DB
+ Int_t GetRegionalFileVersion() const { return fRegionalFileVersion; }
+ /// version of the masks for the local cards in the detector DB
+ Int_t GetLocalMaskFileVersion() const { return fLocalMaskFileVersion; }
+ /// version of the transverse momentum Look-Up-Table in the detector DB
+ Int_t GetLocalLutFileVersion() const { return fLocalLutFileVersion; }
+ /// version of the signature file in the detector DB
+ Int_t GetSignatureFileVersion() const { return fSignatureFileVersion; }
+
+ /// last known version of the global crate configuration
+ Int_t GetGlobalFileLastVersion() const { return fGlobalFileLastVersion; }
+ /// last known version of the regional crate configuration
+ Int_t GetRegionalFileLastVersion() const { return fRegionalFileLastVersion; }
+ /// last known version of the masks for the local cards
+ Int_t GetLocalMaskFileLastVersion() const { return fLocalMaskFileLastVersion; }
+ /// last known version of the transverse momentum Look-Up-Table
+ Int_t GetLocalLutFileLastVersion() const { return fLocalLutFileLastVersion; }
+
+ /// store for the masks for the local cards (own)
+ AliMUONVStore* GetLocalMasks() const { return fLocalMasks; }
+ /// store for the DA-calculated masks for the local cards (own)
+ AliMUONVStore* GetLocalMasksDA() const { return fLocalMasksDA; }
+ /// configuration object for the regional crate (own)
+ AliMUONRegionalTriggerConfig* GetRegionalMasks() const { return fRegionalMasks; }
+ /// configuration object for the global crate (own)
+ AliMUONGlobalCrateConfig* GetGlobalMasks() const { return fGlobalMasks; }
+
+ /// read/write configurations, masks and LUT to/from online files (own)
+ AliMUONTriggerIO* GetTriggerIO() const { return fTriggerIO; }
+
+ /// store for local strips patterns (noisy strips)
+ AliMUONVStore* GetPatternStoreN() const { return fPatternStoreN; }
+ /// store for local strips patterns (dead strips)
+ AliMUONVStore* GetPatternStoreD() const { return fPatternStoreD; }
+
+ /// number of accumulated PHYSICS events for noisy channels analysis
+ Int_t GetEventsN() const { return fEventsN; }
+ /// number of accumulated CALIBRATION events for dead channels analysis
+ Int_t GetEventsD() const { return fEventsD; }
+
+ /// print verbosity of the DA
+ Int_t GetPrintLevel() const { return fPrintLevel; }
+
+ /// select PHYSICS events for noisy channels analysis
+ Bool_t GetAlgoNoisyInput() const { return fAlgoNoisyInput; }
+ /// select CALIBRATION events for dead channels analysis
+ Bool_t GetAlgoDeadcInput() const { return fAlgoDeadcInput; }
+
+ /// threshold for noisy global inputs (fraction of events)
+ Float_t GetThrN() const { return fThrN; }
+ /// threshold for dead global inputs (fraction of events)
+ Float_t GetThrD() const { return fThrD; }
+ /// threshold for noisy local strips (fraction of events)
+ Float_t GetThrLocN() const { return fThrLocN; }
+ /// threshold for dead local strips (fraction of events)
+ Float_t GetThrLocD() const { return fThrLocD; }
+
+ /// minumum nr of events for rate calculation
+ Int_t GetMinEvents() const { return fMinEvents; }
+ /// maximum number of events to analyze
+ Int_t GetMaxEvents() const { return fMaxEvents; }
+ /// number of events to skip from start
+ Int_t GetSkipEvents() const { return fSkipEvents; }
+
+ /// show warnings from the raw data decoder
+ Bool_t WithWarnings() const { return fWithWarnings; }
+ /// use the high-performance (HP) decoder
+ Bool_t UseFastDecoder() const { return fUseFastDecoder; }
+
+ /// number of global input words
+ Int_t GetGlobalInputs() const { return kGlobalInputs; }
+ /// length in bits of a global input word
+ Int_t GetGlobalInputLength() const { return kGlobalInputLength; }
+
+ /// get accumulated values for bit "ib" from global input word "ii", PHYSICS events
+ Int_t GetAccGlobalInputN(Int_t ii, Int_t ib) const { return fAccGlobalInputN[ii][ib]; }
+ /// get accumulated values for bit "ib" from global input word "ii", CALIBRATION events
+ Int_t GetAccGlobalInputD(Int_t ii, Int_t ib) const { return fAccGlobalInputD[ii][ib]; }
+
+ /// set the name of the Start-of-data field in MtgCurrent.dat
+ void SetSodName(Char_t *name) { fSodName = TString(name); }
+ /// set the flag value of the Start-of-data field in MtgCurrent.dat
+ void SetSodFlag(Int_t flag) { fSodFlag = flag; }
+
+ /// set the name of the Detector Algorithm field in MtgCurrent.dat
+ void SetDAName(Char_t *name) { fDAName = TString(name); }
+ /// set the flag value of the Detector Algorithm field in MtgCurrent.dat
+ void SetDAFlag(Int_t flag) { fDAFlag = flag; }
+ /// set DA active mode, 1 = GLOBAL (default), 2 = GLOBAL and LOCAL
+ void SetDAMode(Int_t mode) { fDAMode = mode; }
+
+ /// set the global crate configuration file name
+ void SetGlobalFileName(const Char_t *name) { fGlobalFileName = TString(name); }
+ /// set the regional crate configuration file name
+ void SetRegionalFileName(const Char_t *name) { fRegionalFileName = TString(name); }
+ /// set the masks for the local cards, file name
+ void SetLocalMaskFileName(const Char_t *name) { fLocalMaskFileName = TString(name); }
+ /// set the transverse momentum Look-Up-Table, file name
+ void SetLocalLutFileName(const Char_t *name) { fLocalLutFileName = TString(name); }
+ /// set the signature file name
+ void SetSignatureFileName(const Char_t *name) { fSignatureFileName = TString(name); }
+
+ /// set the version of the global crate configuration in the detector DB
+ void SetGlobalFileVersion(Int_t ver) { fGlobalFileVersion = ver; }
+ /// set the version of the regional crate configuration in the detector DB
+ void SetRegionalFileVersion(Int_t ver) { fRegionalFileVersion = ver; }
+ /// set the version of the masks for the local cards in the detector DB
+ void SetLocalMaskFileVersion(Int_t ver) { fLocalMaskFileVersion = ver; }
+ /// set the version of the transverse momentum Look-Up-Table in the detector DB
+ void SetLocalLutFileVersion(Int_t ver) { fLocalLutFileVersion = ver; }
+ /// set the version of the signature file in the detector DB
+ void SetSignatureFileVersion(Int_t ver) { fSignatureFileVersion = ver; }
+
+ /// set the last known version of the global crate configuration
+ void SetGlobalFileLastVersion(Int_t ver) { fGlobalFileLastVersion = ver; }
+ /// set the last known version of the regional crate configuration
+ void SetRegionalFileLastVersion(Int_t ver) { fRegionalFileLastVersion = ver; }
+ /// set the last known version of the masks for the local cards
+ void SetLocalMaskFileLastVersion(Int_t ver) { fLocalMaskFileLastVersion = ver; }
+ /// set the last known version of the transverse momentum Look-Up-Table
+ void SetLocalLutFileLastVersion(Int_t ver) { fLocalLutFileLastVersion = ver; }
+
+ /// increment the number of selected PHYSICS events
+ void IncNoiseEvent() { fEventsN++; }
+ /// increment the number of selected CALIBRATION events
+ void IncDeadcEvent() { fEventsD++; }
+
+ /// count the value of the bit "ib" of global input word "ii" (PHYSICS events)
+ void AddAccGlobalInputN(Int_t ii, Int_t ib, Int_t val) { fAccGlobalInputN[ii][ib] += val; }
+ /// count the value of the bit "ib" of global input word "ii" (CALIBRATION events)
+ void AddAccGlobalInputD(Int_t ii, Int_t ib, Int_t val) { fAccGlobalInputD[ii][ib] += val; }
+
+ /// set the print verbosity level of the DA
+ void SetPrintLevel(Int_t level) { fPrintLevel = level; }
+
+ /// select PHYSICS events for noisy channels analysis
+ void SetAlgoNoisyInput(Bool_t val) { fAlgoNoisyInput = val; }
+ /// select CALIBRATION events for dead channels analysis
+ void SetAlgoDeadcInput(Bool_t val) { fAlgoDeadcInput = val; }
+
+ /// set the threshold for noisy global inputs (fraction of events)
+ void SetThrN(Float_t val) { fThrN = val; }
+ /// set the threshold for dead global inputs (fraction of events)
+ void SetThrD(Float_t val) { fThrD = val; }
+ /// set the threshold for noisy local strips (fraction of events)
+ void SetThrLocN(Float_t val) { fThrLocN = val; }
+ /// set the threshold for dead local strips (fraction of events)
+ void SetThrLocD(Float_t val) { fThrLocD = val; }
+
+ /// set the minumum nr of events for rate calculation
+ void SetMinEvents(Int_t val) { fMinEvents = val; }
+ /// set the maximum number of events to analyze
+ void SetMaxEvents(Int_t val) { fMaxEvents = val; }
+ /// set the number of events to skip from start
+ void SetSkipEvents(Int_t val) { fSkipEvents = val; }
+
+ /// set/unset to show warnings from the raw data decoder
+ void SetWithWarnings() { fWithWarnings = kTRUE; }
+ /// set/unset the use of the high-performance (HP) decoder
+ void SetUseFastDecoder() { fUseFastDecoder = kTRUE; }
+
+ /// increment version of the global crate configuration file
+ void IncGlobalFileVersion() { fGlobalFileVersion++; }
+ /// increment version of the local mask configuration file
+ void IncLocalMaskFileVersion() { fLocalMaskFileVersion++; }
+ /// count skipped events
+ void DecSkipEvents() { fSkipEvents--; }
+
+private:
+
+ /// copy constructor, not implemented
+ AliDAConfig (const AliDAConfig& cfg);
+ /// assignment operator, not implemented
+ AliDAConfig& operator=(const AliDAConfig& cfg);
+
+ const TString fDAConfigFileName; //!< name of the DA configuration file from detector DB
+ const TString fCurrentFileName; //!< usually MtgCurrent.dat
+ const TString fLastCurrentFileName; //!< usually MtgLastCurrent.dat
+
+ TString fSodName; //!< name of the Start-of-data field in MtgCurrent.dat
+ Int_t fSodFlag; //!< flag value of the Start-of-data field in MtgCurrent.dat
+
+ TString fDAName; //!< name of the Detector Algorithm field in MtgCurrent.dat
+ Int_t fDAFlag; //!< flag value of the Detector Algorithm field in MtgCurrent.dat (enabled/disabled)
+ Int_t fDAMode; //!< DA active mode, GLOBAL or GLOBAL+LOCAL
+
+ TString fGlobalFileName; //!< global crate configuration, file name
+ TString fRegionalFileName; //!< regional crate configuration, file name
+ TString fLocalMaskFileName; //!< masks for the local cards, file name
+ TString fLocalLutFileName; //!< transverse momentum Look-Up-Table, file name
+ TString fSignatureFileName; //!< signature file name
+
+ Int_t fGlobalFileVersion; //!< version of the global crate configuration in the detector DB
+ Int_t fRegionalFileVersion; //!< version of the regional crate configuration in the detector DB
+ Int_t fLocalMaskFileVersion; //!< version of the masks for the local cards in the detector DB
+ Int_t fLocalLutFileVersion; //!< version of the transverse momentum Look-Up-Table in the detector DB
+ Int_t fSignatureFileVersion; //!< version of the signature file in the detector DB
+
+ Int_t fGlobalFileLastVersion; //!< last known version of the global crate configuration
+ Int_t fRegionalFileLastVersion; //!< last known version of the regional crate configuration
+ Int_t fLocalMaskFileLastVersion; //!< last known version of the masks for the local cards
+ Int_t fLocalLutFileLastVersion; //!< last known version of the transverse momentum Look-Up-Table
+
+ Int_t fEventsN; //!< number of accumulated PHYSICS events
+ Int_t fEventsD; //!< number of accumulated CALIBRATION events
+
+ Int_t fPrintLevel; //!< print verbosity of the DA
+
+ AliMUONVStore* fLocalMasks; //!< store for the masks for the local cards
+ AliMUONVStore* fLocalMasksDA; //!< store for the DA-calculated masks for the local cards
+ AliMUONRegionalTriggerConfig* fRegionalMasks; //!< configuration object for the regional crate
+ AliMUONGlobalCrateConfig* fGlobalMasks; //!< configuration object for the global crate
+
+ AliMUONTriggerIO *fTriggerIO; //!< read/write masks and LUT to/from online files
+
+ Bool_t fAlgoNoisyInput; //!< select PHYSICS events for noisy channels analysis
+ Bool_t fAlgoDeadcInput; //!< select CALIBRATION events for dead channels analysis
+
+ Float_t fThrN; //!< threshold for noisy global inputs (fraction of events)
+ Float_t fThrD; //!< threshold for dead global inputs (fraction of events)
+ Float_t fThrLocN; //!< threshold for noisy local strips (fraction of events)
+ Float_t fThrLocD; //!< threshold for dead local strips (fraction of events)
+ Int_t fMinEvents; //!< minumum nr of events for rate calculation
+ Int_t fSkipEvents; //!< number of events to skip from start
+ Int_t fMaxEvents; //!< maximum number of events to analyze
+ Bool_t fWithWarnings; //!< show warnings from the raw data decoder
+ Bool_t fUseFastDecoder; //!< use the high-performance (HP) decoder
+
+ const Int_t fNLocalBoard; //!< number of local boards
+
+ enum { kGlobalInputs = 4, //!< number of global input words
+ kGlobalInputLength = 32 //!< length in bits of a global input word
+ };
+
+ AliMUONVStore *fPatternStoreN; //! store for local strips patterns
+ AliMUONVStore *fPatternStoreD; //! store for local strips patterns
+
+ Int_t fAccGlobalInputN[kGlobalInputs][kGlobalInputLength]; //!< storage for global input (PHYSICS events)
+ Int_t fAccGlobalInputD[kGlobalInputs][kGlobalInputLength]; //!< storage for global input (CALIBRATION events)
+
+};
+
+//__________________________________________________________________
+Bool_t ReadDAConfig(AliDAConfig& cfg)
+{
+ /// read run parameters for the DA
+
+ char line[80];
+
+ TString file;
+ file = cfg.GetDAConfigFileName();
+ std::ifstream in(gSystem->ExpandPathName(file.Data()));
+ if (!in.good()) {
+ printf("Cannot open DA configuration file %s ; use default values.\n",file.Data());
+ return kTRUE;
+ }
+
+ TString tmp;
+ Int_t pos;
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetPrintLevel(tmp.Atoi());
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetThrN(tmp.Atof());
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetThrD(tmp.Atof());
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetMinEvents(tmp.Atoi());
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetSkipEvents(tmp.Atoi());
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetMaxEvents(tmp.Atoi());
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ if (tmp.Atoi() != 0) cfg.SetWithWarnings();
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ if (tmp.Atoi() != 0) cfg.SetUseFastDecoder();
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetThrLocN(tmp.Atof());
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetThrLocD(tmp.Atof());
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp = tmp(0,pos);
+ cfg.SetDAMode(tmp.Atoi());
+
+ return kTRUE;
+
+}
+
+//__________________________________________________________________
+void WriteLastCurrentFile(AliDAConfig& cfg, TString currentFile)
+{
+ /// write last current file
+
+ ofstream out;
+ TString file;
+ file = currentFile;
+ out.open(file.Data());
+ out << cfg.GetSodName() << " " << cfg.GetSodFlag() << endl;
+ out << cfg.GetDAName() << " " << cfg.GetDAFlag() << endl;
+
+ out << cfg.GetGlobalFileName() << " " << cfg.GetGlobalFileVersion() << endl;
+ out << cfg.GetRegionalFileName() << " " << cfg.GetRegionalFileVersion() << endl;
+ out << cfg.GetLocalMaskFileName() << " " << cfg.GetLocalMaskFileVersion() << endl;
+ out << cfg.GetLocalLutFileName() << " " << cfg.GetLocalLutFileVersion() << endl;
+ out << cfg.GetSignatureFileName() << " " << cfg.GetSignatureFileVersion() << endl;
+
+ out.close();
+}
+
+//___________________________________________________________________________________________
+Bool_t ReadCurrentFile(AliDAConfig& cfg, TString currentFile, Bool_t lastCurrentFlag = kFALSE)
+{
+ /// read last current file name and version
+
+ char line[80];
+ char name[80];
+ Int_t flag;
+
+ TString file;
+ file = currentFile;
+ std::ifstream in(gSystem->ExpandPathName(file.Data()));
+ if (!in.good()) {
+ printf("Cannot open last current file %s\n",currentFile.Data());
+ return kFALSE;
+ }
+
+ // read SOD
+ in.getline(line,80);
+ sscanf(line, "%s %d", name, &flag);
+ cfg.SetSodName(name);
+ cfg.SetSodFlag(flag);
+ if (cfg.GetPrintLevel()) printf("Sod Flag %d\n", cfg.GetSodFlag());
+
+ //read DA
+ in.getline(line,80);
+ sscanf(line, "%s %d", name, &flag);
+ cfg.SetDAName(name);
+ cfg.SetDAFlag(flag);
+ if (cfg.GetPrintLevel()) printf("DA Flag: %d\n", cfg.GetDAFlag());
+
+ // read global
+ in.getline(line,80);
+ TString tmp(line);
+ Int_t pos = tmp.First(" ");
+ TString tmp1 = tmp(0, pos);
+ cfg.SetGlobalFileName(tmp1.Data());
+
+ if (!lastCurrentFlag) {
+ cfg.SetGlobalFileVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Global File Name: %s version: %d\n",
+ cfg.GetGlobalFileName(), cfg.GetGlobalFileVersion());
+ } else {
+ cfg.SetGlobalFileLastVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Global File Name: %s last version: %d\n",
+ cfg.GetGlobalFileName(), cfg.GetGlobalFileLastVersion());
+ }
+
+ // read regional
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp1 = tmp(0, pos);
+ cfg.SetRegionalFileName(tmp1.Data());
+
+ if (!lastCurrentFlag) {
+ cfg.SetRegionalFileVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Regional File Name: %s version: %d\n",
+ cfg.GetRegionalFileName(), cfg.GetRegionalFileVersion());
+
+ } else {
+ cfg.SetRegionalFileLastVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Regional File Name: %s last version: %d\n",
+ cfg.GetRegionalFileName(), cfg.GetRegionalFileLastVersion());
+ }
+
+ // read mask
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp1 = tmp(0, pos);
+ cfg.SetLocalMaskFileName(tmp1.Data());
+
+ if (!lastCurrentFlag) {
+ cfg.SetLocalMaskFileVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Mask File Name: %s version: %d\n",
+ cfg.GetLocalMaskFileName(), cfg.GetLocalMaskFileVersion());
+ } else {
+ cfg.SetLocalMaskFileLastVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Mask File Name: %s last version: %d\n",
+ cfg.GetLocalMaskFileName(), cfg.GetLocalMaskFileLastVersion());
+ }
+ // read Lut
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp1 = tmp(0, pos);
+ cfg.SetLocalLutFileName(tmp1.Data());
+
+ if (!lastCurrentFlag) {
+ cfg.SetLocalLutFileVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Lut File Name: %s version: %d\n",
+ cfg.GetLocalLutFileName(), cfg.GetLocalLutFileVersion());
+ } else {
+ cfg.SetLocalLutFileLastVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Lut File Name: %s last version: %d\n",
+ cfg.GetLocalLutFileName(), cfg.GetLocalLutFileLastVersion());
+ }
+
+ in.getline(line,80);
+ tmp = line;
+ pos = tmp.First(" ");
+ tmp1 = tmp(0, pos);
+ cfg.SetSignatureFileName(tmp1.Data());
+ cfg.SetSignatureFileVersion(atoi(tmp(pos+1, tmp.Length()-pos).Data()));
+ if (cfg.GetPrintLevel()) printf("Lut File Name: %s version: %d\n",
+ cfg.GetSignatureFileName(), cfg.GetSignatureFileVersion());
+
+ return kTRUE;
+}
+
+//_____________
+void ReadFileNames(AliDAConfig& cfg)
{
+ /// if last current file does not exist than read current file
+
+ if (!ReadCurrentFile(cfg,cfg.GetLastCurrentFileName(), kTRUE))
+ {
+ ReadCurrentFile(cfg,cfg.GetCurrentFileName(), kTRUE);
+ }
- // needed for streamer application
- gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo",
- "*",
- "TStreamerInfo",
- "RIO",
- "TStreamerInfo()");
+ // any case read current file
+ ReadCurrentFile(cfg,cfg.GetCurrentFileName());
+
+}
+//__________________
+Bool_t ExportFiles(AliDAConfig& cfg)
+{
+ /// Export files to FES
- Int_t printLevel = 0; // Global variable defined as extern in the others .cxx files
- Int_t skipEvents = 0;
- Int_t maxEvents = 1000000;
- Char_t inputFile[256];
- TString flatOutputFile;
- TString configFile;
+ // env variables have to be set (suppose by ECS ?)
+ // setenv DATE_FES_PATH
+ // setenv DATE_RUN_NUMBER
+ // setenv DATE_ROLE_NAME
+ // setenv DATE_DETECTOR_CODE
-// option handler
+#ifdef OFFLINE
+ gSystem->Setenv("DAQDALIB_PATH", "$DATE_SITE/infoLogger");
+ gSystem->Setenv("DAQDA_TEST_DIR", "/alisoft/FES");
+#endif
- // decode the input line
- for (Int_t i = 1; i < argc; i++) // argument 0 is the executable name
- {
- Char_t* arg;
-
- arg = argv[i];
- if (arg[0] != '-') continue;
- switch (arg[1])
- {
- case 'f' :
- i++;
- sprintf(inputFile,argv[i]);
- break;
- case 'a' :
- i++;
- flatOutputFile = argv[i];
- break;
- case 'c' :
- i++;
- configFile = argv[i];
- break;
- case 'e' :
- i++;
- command = argv[i];
- break;
- case 'd' :
- i++;
- printLevel=atoi(argv[i]);
- break;
- case 's' :
- i++;
- skipEvents=atoi(argv[i]);
- break;
- case 'n' :
- i++;
- sscanf(argv[i],"%d",&maxEvents);
- break;
- case 'h' :
- i++;
- printf("\n******************* %s usage **********************",argv[0]);
- printf("\n%s -options, the available options are :",argv[0]);
- printf("\n-h help (this screen)");
- printf("\n");
- printf("\n Input");
- printf("\n-f <raw data file> (default = %s)",inputFile);
- printf("\n");
- printf("\n Output");
- printf("\n-a <Flat ASCII file> (default = %s)",flatOutputFile.Data());
- printf("\n");
- printf("\n Options");
- printf("\n-d <print level> (default = %d)",printLevel);
- printf("\n-s <skip events> (default = %d)",skipEvents);
- printf("\n-n <max events> (default = %d)",maxEvents);
- printf("\n-e <execute pattern/scaler> (default = %s)",command.Data());
-
- printf("\n\n");
- exit(-1);
- default :
- printf("%s : bad argument %s (please check %s -h)\n",argv[0],argv[i],argv[0]);
- argc = 2; exit(-1); // exit if error
- } // end of switch
- } // end of for i
-
- // set command to lower case
- command.ToLower();
-
- // decoding the events
+ // update files
+ Int_t status = 0;
+
+ Bool_t modified = kFALSE;
+ Bool_t globalExported = kFALSE;
+
+ ofstream out;
+ TString fileExp("ExportedFiles.dat");
+ TString file;
+
+ out.open(fileExp.Data());
+ if (!out.good()) {
+ printf("Failed to create file: %s\n",file.Data());
+ return kFALSE;
+ }
+
+ // check if MtgLastCurrent.dat exists
+ // if not, do initial export of all files
+ Bool_t initFES = kFALSE;
+ if (gSystem->AccessPathName("MtgLastCurrent.dat"))
+ initFES = kTRUE;
+ if (initFES) printf("Copy all configuration files to the FES.\n");
+
+ file = cfg.GetLocalMaskFileName();
+ if ((cfg.GetLocalMaskFileLastVersion() != cfg.GetLocalMaskFileVersion()) || initFES) {
+ modified = kTRUE;
+ status = daqDA_FES_storeFile(file.Data(), "LOCAL");
+ if (status) {
+ printf("Failed to export file: %s\n",cfg.GetLocalMaskFileName());
+ return kFALSE;
+ }
+ if(cfg.GetPrintLevel()) printf("Export file: %s\n",cfg.GetLocalMaskFileName());
+ out << cfg.GetLocalMaskFileName() << endl;
+ }
+
+ file = cfg.GetLocalLutFileName();
+ if ((cfg.GetLocalLutFileLastVersion() != cfg.GetLocalLutFileVersion()) || initFES) {
+ modified = kTRUE;
+ status = daqDA_FES_storeFile(file.Data(), "LUT");
+ if (status) {
+ printf("Failed to export file: %s\n",cfg.GetLocalLutFileName());
+ return kFALSE;
+ }
+ if(cfg.GetPrintLevel()) printf("Export file: %s\n",cfg.GetLocalLutFileName());
+ out << cfg.GetLocalLutFileName() << endl;
+
+ }
+
+ file = cfg.GetGlobalFileName();
+ if ((cfg.GetGlobalFileLastVersion() != cfg.GetGlobalFileVersion()) || modified || initFES) {
+ modified = kTRUE;
+ globalExported = kTRUE;
+ status = daqDA_FES_storeFile(file.Data(), "GLOBAL");
+ if (status) {
+ printf("Failed to export file: %s\n",cfg.GetGlobalFileName());
+ return kFALSE;
+ }
+ if(cfg.GetPrintLevel()) printf("Export file: %s\n",cfg.GetGlobalFileName());
+ out << cfg.GetGlobalFileName() << endl;
+ }
+
+ file = cfg.GetRegionalFileName();
+ if ( (cfg.GetRegionalFileLastVersion() != cfg.GetRegionalFileVersion()) || modified || initFES) {
+ status = daqDA_FES_storeFile(file.Data(), "REGIONAL");
+ if (status) {
+ printf("Failed to export file: %s\n",cfg.GetRegionalFileName());
+ return kFALSE;
+ }
+ if(cfg.GetPrintLevel()) printf("Export file: %s\n",cfg.GetRegionalFileName());
+ out << cfg.GetRegionalFileName() << endl;
+
+ // needed for the initialisation of the mapping
+ if (!globalExported) {
+ file = cfg.GetGlobalFileName();
+ status = daqDA_FES_storeFile(file.Data(), "GLOBAL");
+ if (status) {
+ printf("Failed to export file: %s\n",cfg.GetGlobalFileName());
+ return kFALSE;
+ }
+ if(cfg.GetPrintLevel()) printf("Export file: %s\n",cfg.GetGlobalFileName());
+ out << cfg.GetGlobalFileName() << endl;
+ }
+
+ }
+
+ out.close();
+
+ // export Exported file to FES anyway
+ status = daqDA_FES_storeFile(fileExp.Data(), "EXPORTED");
+ if (status) {
+ printf("Failed to export file: %s\n", fileExp.Data());
+ return kFALSE;
+ }
+ if(cfg.GetPrintLevel()) printf("Export file: %s\n",fileExp.Data());
+
+ // write last current file
+ WriteLastCurrentFile(cfg,cfg.GetLastCurrentFileName());
+
+ return kTRUE;
+}
+
+//__________________
+Bool_t ImportFiles(AliDAConfig& cfg)
+{
+ /// copy locally a file from daq detector config db
+ /// The current detector is identified by detector code in variable
+ /// DATE_DETECTOR_CODE. It must be defined.
+ /// If environment variable DAQDA_TEST_DIR is defined, files are copied from
+ /// DAQDA_TEST_DIR instead of the database.
+ /// The usual environment variables are not needed.
+
+ Int_t status = 0;
+
+#ifdef OFFLINE
+ gSystem->Setenv("DAQDALIB_PATH", "$DATE_SITE/db");
+#endif
+
+ status = daqDA_DB_getFile(cfg.GetDAConfigFileName(), cfg.GetDAConfigFileName());
+ if (status) {
+ printf("Failed to get DA config file from DB: %s\n",cfg.GetDAConfigFileName());
+ return kFALSE;
+ }
+
+ ReadDAConfig(cfg);
+
+ status = daqDA_DB_getFile(cfg.GetCurrentFileName(), cfg.GetCurrentFileName());
+ if (status) {
+ printf("Failed to get current config file from DB: %s\n",cfg.GetCurrentFileName());
+ return kFALSE;
+ }
+
+ ReadFileNames(cfg);
+
+ status = daqDA_DB_getFile(cfg.GetGlobalFileName(), cfg.GetGlobalFileName());
+ if (status) {
+ printf("Failed to get current config file from DB: %s\n", cfg.GetGlobalFileName());
+ return kFALSE;
+ }
+
+ status = daqDA_DB_getFile(cfg.GetRegionalFileName(), cfg.GetRegionalFileName());
+ if (status) {
+ printf("Failed to get current config file from DB: %s\n",cfg.GetRegionalFileName());
+ return kFALSE;
+ }
+
+ status = daqDA_DB_getFile(cfg.GetLocalMaskFileName(), cfg.GetLocalMaskFileName());
+ if (status) {
+ printf("Failed to get current config file from DB: %s\n",cfg.GetLocalMaskFileName());
+ return kFALSE;
+ }
+
+ status = daqDA_DB_getFile(cfg.GetLocalLutFileName(), cfg.GetLocalLutFileName());
+ if (status) {
+ printf("Failed to get current config file from DB: %s\n",cfg.GetLocalLutFileName());
+ return kFALSE;
+ }
+
+ return kTRUE;
+}
+
+//_____________
+void ReadMaskFiles(AliDAConfig& cfg)
+{
+ /// read mask files
+
+ const Char_t* localFile = cfg.GetLocalMaskFileName();
+ const Char_t* regionalFile = cfg.GetRegionalFileName();
+ const Char_t* globalFile = cfg.GetGlobalFileName();
+
+ cfg.GetTriggerIO()->ReadConfig(localFile, regionalFile, globalFile, cfg.GetLocalMasks(), cfg.GetRegionalMasks(), cfg.GetGlobalMasks());
+}
+
+//______________________________________________________________
+UInt_t GetFetMode(const AliDAConfig & cfg)
+{
+ /// FET mode = 3 to run algorithm for dead global inputs
+ /// 0x3 prepulse
+ /// 0x0 internal
+
+ return cfg.GetGlobalMasks()->GetFetRegister(3);
+
+}
+
+//______________________________________________________________
+void StoreGlobalInput(AliDAConfig& cfg, const UInt_t * const globalInput)
+{
+ /// accumulate and build statistics of global input values
- Int_t status;
- Int_t nDateEvents = 0;
+ for (Int_t ii = 0; ii < cfg.GetGlobalInputs(); ii++) {
+ for (Int_t ib = 0; ib < cfg.GetGlobalInputLength(); ib++) {
+ // lsb -> msb
+ if (cfg.GetAlgoNoisyInput())
+ cfg.AddAccGlobalInputN(ii,ib,((globalInput[ii] >> ib) & 0x1));
+ if (cfg.GetAlgoDeadcInput())
+ cfg.AddAccGlobalInputD(ii,ib,((globalInput[ii] >> ib) & 0x1));
+ }
+ }
- void* event;
+}
+
+//______________________________________________________________
+void UpdateGlobalMasks(AliDAConfig& cfg)
+{
+ /// update the global masks
+
+#ifdef OFFLINE
+ gSystem->Setenv("DAQDALIB_PATH", "$DATE_SITE/db");
+#endif
- // containers
- AliMUONDDLTrigger* ddlTrigger = 0x0;
- AliMUONDarcHeader* darcHeader = 0x0;
- AliMUONRegHeader* regHeader = 0x0;
- AliMUONLocalStruct* localStruct = 0x0;
+ Float_t rateN = 0.0, rateD = 0.0;
+ UInt_t gmask[4], omask;
+ Bool_t noise, deadc, withEvN, withEvD, updated = kFALSE;
- TStopwatch timers;
+ for (Int_t ii = 0; ii < cfg.GetGlobalInputs(); ii++) {
+ gmask[ii] = 0;
- timers.Start(kTRUE);
+ for (Int_t ib = 0; ib < cfg.GetGlobalInputLength(); ib++) {
+ // lsb -> msb
+ noise = kFALSE;
+ deadc = kFALSE;
+ withEvN = kFALSE;
+ withEvD = kFALSE;
+ if (cfg.GetEventsN() > cfg.GetMinEvents()) {
+ rateN = (Float_t)cfg.GetAccGlobalInputN(ii,ib)/(Float_t)cfg.GetEventsN();
+ noise = (rateN > cfg.GetThrN());
+ withEvN = kTRUE;
+ }
+ if (cfg.GetEventsD() > cfg.GetMinEvents()) {
+ rateD = (Float_t)cfg.GetAccGlobalInputD(ii,ib)/(Float_t)cfg.GetEventsD();
+ deadc = (rateD < cfg.GetThrD());
+ withEvD = kTRUE;
+ }
+ if (!withEvN && !withEvD) {
+ // - copy the bit from the old mask
+ gmask[ii] |= ((cfg.GetGlobalMasks()->GetGlobalMask(ii) >> ib) & 0x1) << ib;
+ if (cfg.GetPrintLevel())
+ printf("Mask not changed (just copy the old values)\n");
+ }
+ if (!withEvN && withEvD) {
+ if (!deadc) {
+ // - create a new mask, set the bit to 1
+ // not allowed!
+ //gmask[ii] |= 0x1 << ib;
+ // - copy the bit from the old mask
+ gmask[ii] |= ((cfg.GetGlobalMasks()->GetGlobalMask(ii) >> ib) & 0x1) << ib;
+ } else {
+ // - create a new mask, set the bit to 0
+ gmask[ii] |= 0x0 << ib;
+ if (cfg.GetPrintLevel())
+ printf("Found dead channel %1d:%02d (%4.2f) \n",ii,ib,rateD);
+ }
+ }
+ if (withEvN && !withEvD) {
+ if (!noise) {
+ // - create a new mask, set the bit to 1
+ // not allowed!
+ //gmask[ii] |= 0x1 << ib;
+ // - copy the bit from the old mask
+ gmask[ii] |= ((cfg.GetGlobalMasks()->GetGlobalMask(ii) >> ib) & 0x1) << ib;
+ } else {
+ // - create a new mask, set the bit to 0
+ gmask[ii] |= 0x0 << ib;
+ if (cfg.GetPrintLevel())
+ printf("Found noisy channel %1d:%02d (%4.2f) \n",ii,ib,rateN);
+ }
+ }
+ if (withEvN && withEvD) {
+ if (!noise && !deadc) {
+ // - create a new mask, set the bit to 1
+ // not allowed!
+ //gmask[ii] |= 0x1 << ib;
+ // - copy the bit from the old mask
+ gmask[ii] |= ((cfg.GetGlobalMasks()->GetGlobalMask(ii) >> ib) & 0x1) << ib;
+ } else {
+ // - create a new mask, set the bit to 0
+ gmask[ii] |= 0x0 << ib;
+ if (cfg.GetPrintLevel()) {
+ if (noise)
+ printf("Found noisy channel %1d:%02d (%4.2f) \n",ii,ib,rateN);
+ if (deadc)
+ printf("Found dead channel %1d:%02d (%4.2f) \n",ii,ib,rateD);
+ }
+ }
+ }
+ }
+ }
- // once we have a configuration file in db
- // copy locally a file from daq detector config db
- // The current detector is identified by detector code in variable
- // DATE_DETECTOR_CODE. It must be defined.
- // If environment variable DAQDA_TEST_DIR is defined, files are copied from DAQDA_TEST_DIR
- // instead of the database. The usual environment variables are not needed.
- if (!configFile.IsNull()) {
- status = daqDA_DB_getFile("myconfig", configFile.Data());
+ // check if at least one mask value has been changed from previous version
+ for (Int_t ii = 0; ii < cfg.GetGlobalInputs(); ii++) {
+ printf("Global mask [%1d] %08x \n",ii,gmask[ii]);
+ omask = cfg.GetGlobalMasks()->GetGlobalMask(ii);
+ if (gmask[ii] != omask) {
+ updated = kTRUE;
+ cfg.GetGlobalMasks()->SetGlobalMask(ii,gmask[ii]);
+ }
+ }
+
+ Int_t status = 0;
+ if (updated) {
+
+ // update version
+ cfg.IncGlobalFileVersion();
+
+ // don't change the file version ("-x.dat")
+
+ cfg.GetTriggerIO()->WriteGlobalConfig(cfg.GetGlobalFileName(),cfg.GetGlobalMasks());
+
+ // write last current file
+ WriteLastCurrentFile(cfg,cfg.GetCurrentFileName());
+
+ status = daqDA_DB_storeFile(cfg.GetGlobalFileName(), cfg.GetGlobalFileName());
if (status) {
- printf("Failed to get config file : %d\n",status);
- return -1;
+ printf("Failed to export file to DB: %s\n",cfg.GetGlobalFileName());
+ return;
}
+
+ status = daqDA_DB_storeFile(cfg.GetCurrentFileName(), cfg.GetCurrentFileName());
+ if (status) {
+ printf("Failed to export file to DB: %s\n",cfg.GetCurrentFileName());
+ return;
+ }
+
}
+
+}
+
+//______________________________________________________________
+void UpdateLocalMask(AliDAConfig& cfg, Int_t localBoardId, Int_t connector, Int_t strip)
+{
+
+ /// update local strip mask
+ AliMUONVCalibParam* localMask =
+ static_cast<AliMUONVCalibParam*>(cfg.GetLocalMasksDA()->FindObject(localBoardId));
+
+ UShort_t mask = static_cast<UShort_t>(localMask->ValueAsInt(connector,0));
+
+ mask &= ~(0x1 << strip); // set strip mask to zero
+
+ localMask->SetValueAsInt(connector, 0, mask);
+
+}
- status = monitorSetDataSource(inputFile);
- if (status) {
- cerr << "ERROR : monitorSetDataSource status (hex) = " << hex << status
- << " " << monitorDecodeError(status) << endl;
- return -1;
+//______________________________________________________________
+void MakePattern(AliDAConfig& cfg, Int_t localBoardId, const TArrayS& xPattern, const TArrayS& yPattern)
+{
+ /// calculate the hit map for each strip in x and y direction
+
+ AliMUONVCalibParam* pat = 0x0;
+
+ if (cfg.GetAlgoNoisyInput())
+ pat = static_cast<AliMUONVCalibParam*>(cfg.GetPatternStoreN()->FindObject(localBoardId));
+ if (cfg.GetAlgoDeadcInput())
+ pat = static_cast<AliMUONVCalibParam*>(cfg.GetPatternStoreD()->FindObject(localBoardId));
+
+ if (!pat) {
+ pat = new AliMUONCalibParamND(2, 64, localBoardId, 0, 0.);
+ if (cfg.GetAlgoNoisyInput())
+ cfg.GetPatternStoreN()->Add(pat);
+ if (cfg.GetAlgoDeadcInput())
+ cfg.GetPatternStoreD()->Add(pat);
}
- status = monitorDeclareMp("MUON Tracking monitoring");
- if (status) {
- cerr << "ERROR : monitorDeclareMp status (hex) = " << hex << status
- << " " << monitorDecodeError(status) << endl;
- return -1;
+
+ for (Int_t i = 0; i < 4; ++i) {
+ for (Int_t j = 0; j < 16; ++j) {
+
+ Int_t xMask = xPattern[i];
+ Int_t yMask = yPattern[i];
+
+ Int_t index = 16*i + j;
+ Double_t patOcc = 0.;
+
+ if ( (xMask >> j ) & 0x1 ) {
+ patOcc = pat->ValueAsDouble(index, 0) + 1.;
+ pat->SetValueAsDouble(index, 0, patOcc);
+ }
+ if ( (yMask >> j ) & 0x1 ) {
+ patOcc = pat->ValueAsDouble(index, 1) + 1.;
+ pat->SetValueAsDouble(index, 1, patOcc);
+ }
+ }
}
- cout << "MUONTRKda : Reading data from file " << inputFile <<endl;
+}
- while(1)
- {
- if (nEvents >= maxEvents) break;
- if (nEvents && nEvents % 100 == 0)
- cout<<"Cumulated events " << nEvents << endl;
+//______________________________________________________________
+void MakePatternStore(AliDAConfig& cfg)
+{
+ /// analyse patterns for local strips (calculate occupancy)
+
+#ifdef OFFLINE
+ gSystem->Setenv("DAQDALIB_PATH", "$DATE_SITE/db");
+#endif
- // check shutdown condition
- if (daqDA_checkShutdown())
- break;
+ AliMUONVCalibParam* pat;
+ Int_t localBoardId = 0;
+ Int_t nEventsN = 0, nEventsD = 0;
+ UShort_t strip = 0;
+ Int_t connector = 0;
+ Bool_t updated = kFALSE;
+
+ if (cfg.GetEventsN() > cfg.GetMinEvents()) {
+ nEventsN = cfg.GetEventsN();
+ }
+ if (cfg.GetEventsD() > cfg.GetMinEvents()) {
+ nEventsD = cfg.GetEventsD();
+ }
- // Skip Events if needed
- while (skipEvents) {
- status = monitorGetEventDynamic(&event);
- skipEvents--;
+ // noisy strips
+ if (nEventsN > 0) {
+
+ TIter next(cfg.GetPatternStoreN()->CreateIterator());
+
+ while ( ( pat = dynamic_cast<AliMUONVCalibParam*>(next() ) ) ) {
+
+ localBoardId = pat->ID0();
+
+ for (Int_t index = 0; index < pat->Size(); index++) {
+
+ Double_t patXOcc = pat->ValueAsDouble(index, 0)/(Double_t)nEventsN;
+ Double_t patYOcc = pat->ValueAsDouble(index, 1)/(Double_t)nEventsN;
+
+ pat->SetValueAsDouble(index, 0, patXOcc);
+ pat->SetValueAsDouble(index, 1, patYOcc);
+
+ // check for x strip
+ if (patXOcc > cfg.GetThrLocN()) {
+ strip = index % 16;
+ connector = index/16;
+ UpdateLocalMask(cfg, localBoardId, connector, strip);
+ }
+ // check for y strip
+ if (patYOcc > cfg.GetThrLocN()) {
+ strip = index % 16;
+ connector = index/16 + 4;
+ UpdateLocalMask(cfg, localBoardId, connector, strip);
+ }
+
+ }
}
+
+ }
+
+ // dead strips
+ if (nEventsD > 0) {
+
+ TIter next(cfg.GetPatternStoreD()->CreateIterator());
- // starts reading
- status = monitorGetEventDynamic(&event);
- if (status < 0) {
- cout<<"EOF found"<<endl;
- break;
+ while ( ( pat = dynamic_cast<AliMUONVCalibParam*>(next() ) ) ) {
+
+ localBoardId = pat->ID0();
+
+ for (Int_t index = 0; index < pat->Size(); index++) {
+
+ Double_t patXOcc = pat->ValueAsDouble(index, 0)/(Double_t)nEventsD;
+ Double_t patYOcc = pat->ValueAsDouble(index, 1)/(Double_t)nEventsD;
+
+ pat->SetValueAsDouble(index, 0, patXOcc);
+ pat->SetValueAsDouble(index, 1, patYOcc);
+
+ // check for x strip
+ if (patXOcc < cfg.GetThrLocD()) {
+ strip = index % 16;
+ connector = index/16;
+ UpdateLocalMask(cfg, localBoardId, connector, strip);
+ }
+ // check for y strip
+ if (patYOcc < cfg.GetThrLocD()) {
+ strip = index % 16;
+ connector = index/16 + 4;
+ UpdateLocalMask(cfg, localBoardId, connector, strip);
+ }
+
+ }
+ }
+
+ }
+
+ // make and AND with the previous version of the mask and
+ // check if the mask has changed
+ UShort_t maskDA, mask;
+ Int_t nMaskBits = AliMpConstants::TotalNofLocalBoards()*8*16;
+ Int_t nMaskBitsChanged = 0;
+ for (localBoardId = 1; localBoardId <= AliMpConstants::TotalNofLocalBoards(); localBoardId++) {
+ AliMUONVCalibParam* localMaskDA = static_cast<AliMUONVCalibParam*>(cfg.GetLocalMasksDA()->FindObject(localBoardId));
+ AliMUONVCalibParam* localMask = static_cast<AliMUONVCalibParam*>(cfg.GetLocalMasks()->FindObject(localBoardId));
+ for (connector = 0; connector < 8; connector++) {
+ maskDA = static_cast<UShort_t>(localMaskDA->ValueAsInt(connector,0));
+ mask = static_cast<UShort_t>(localMask->ValueAsInt(connector,0));
+ maskDA &= mask;
+ localMaskDA->SetValueAsInt(connector, 0, maskDA);
+ if (maskDA != mask) {
+ updated = kTRUE;
+ // calculated percentage of mask bits changed
+ for (Int_t iBit = 0; iBit < 16; iBit++) {
+ if (((maskDA >> iBit) & 0x1) != ((mask >> iBit) &0x1)) {
+ nMaskBitsChanged++;
+ }
+ }
+ }
}
+ }
- nDateEvents++;
+ printf("LOCAL mask bits changed = %5d (%7.3f %%) \n",nMaskBitsChanged,100*(Float_t)nMaskBitsChanged/(Float_t)nMaskBits);
- // decoding rawdata headers
- AliRawReader *rawReader = new AliRawReaderDate(event);
-
- Int_t eventType = rawReader->GetType();
- runNumber = rawReader->GetRunNumber();
+ Int_t status = 0;
+ if (updated) {
+
+ // update version
+ cfg.IncLocalMaskFileVersion();
+
+ // don't change the file version ("-x.dat")
+
+ cfg.GetTriggerIO()->WriteLocalMasks(cfg.GetLocalMaskFileName(),*cfg.GetLocalMasksDA());
+
+ // write last current file
+ WriteLastCurrentFile(cfg,cfg.GetCurrentFileName());
+
+ status = daqDA_DB_storeFile(cfg.GetLocalMaskFileName(), cfg.GetLocalMaskFileName());
+ if (status) {
+ printf("Failed to export file to DB: %s\n",cfg.GetLocalMaskFileName());
+ return;
+ }
+
+ status = daqDA_DB_storeFile(cfg.GetCurrentFileName(), cfg.GetCurrentFileName());
+ if (status) {
+ printf("Failed to export file to DB: %s\n",cfg.GetCurrentFileName());
+ return;
+ }
+
+ }
+
+}
+
+//*************************************************************//
+int main(Int_t argc, Char_t **argv)
+{
+ /// main routine
+
+ // needed for streamer application
+ gROOT->GetPluginManager()->AddHandler("TVirtualStreamerInfo", "*", "TStreamerInfo", "RIO", "TStreamerInfo()");
+
+ /* check that we got some arguments = list of files */
+ if (argc<2) {
+ printf("Wrong number of arguments\n");
+ return -1;
+ }
+
+ AliDAConfig cfg;
+
+ Char_t inputFile[256] = "";
+ inputFile[0] = 0;
+ if (argc > 1)
+ if (argv[1] != NULL)
+ strncpy(inputFile, argv[1], 256);
+ else {
+ printf("MUONTRGda : No input File !\n");
+ return -1;
+ }
+
+ // decoding the events
+
+ Int_t status = 0;
+ Int_t nDateEvents = 0;
+
+ void* event;
+
+ // containers
+ // old decoder
+ AliMUONDDLTrigger* ddlTrigger = 0x0;
+ AliMUONDarcHeader* darcHeader = 0x0;
+ AliMUONRegHeader* regHeader = 0x0;
+ AliMUONLocalStruct* localStruct = 0x0;
+ // new (fast) decoder
+ const AliMUONRawStreamTriggerHP::AliHeader* darcHeaderHP = 0x0;
+ const AliMUONRawStreamTriggerHP::AliLocalStruct* localStructHP = 0x0;
+
+ TArrayS xPattern(4);
+ TArrayS yPattern(4);
+ Int_t localBoardId = 0;
+
+ TStopwatch timers;
+
+ timers.Start(kTRUE);
+
+ // comment out, since we do not retrieve files from database
+ if (!ImportFiles(cfg)) {
+ printf("Import from DB failed\n");
+ printf("For local test set DAQDA_TEST_DIR to the local directory where the Mtg files are located \n");
+ return -1;
+ }
+ ReadMaskFiles(cfg);
- if (eventType != PHYSICS_EVENT)
- continue; // for the moment
+#ifdef OFFLINE
+ // the run number extracted from the file name
+ TString tmp(inputFile);
+ Int_t pos1 = tmp.First('d');
+ Int_t pos2 = tmp.Last('.');
+ Int_t len = pos2 - (pos1+3);
+ tmp = tmp(pos1+3,len);
+ gSystem->Setenv("DATE_RUN_NUMBER",tmp.Data());
+ gSystem->Exec("echo \"DATE_RUN_NUMBER = \" $DATE_RUN_NUMBER");
+#endif
- nEvents++;
- if (printLevel) printf("\nEvent # %d\n",nEvents);
+ if(!ExportFiles(cfg)) {
+ printf("ExportFiles failed\n");
+ return -1;
+ }
- // decoding MUON payload
- AliMUONRawStreamTrigger* rawStream = new AliMUONRawStreamTrigger(rawReader);
- //rawStream->SetMaxReg(1);
+ if (!cfg.GetDAFlag()) {
- // loops over DDL
- while((status = rawStream->NextDDL())) {
+ cout << "MUONTRGda: DA enable: " << cfg.GetDAFlag() << endl;
+ cout << "MUONTRGda: Print level: " << cfg.GetPrintLevel() << endl;
+
+ printf("MUONTRGda: Execution time : R:%7.2fs C:%7.2fs\n", timers.RealTime(), timers.CpuTime());
- if (printLevel) printf("iDDL %d\n", rawStream->GetDDL());
+ return status;
- ddlTrigger = rawStream->GetDDLTrigger();
- darcHeader = ddlTrigger->GetDarcHeader();
+ }
- if (printLevel) printf("Global output %x\n", (Int_t)darcHeader->GetGlobalOutput());
+ // FET is triggered by CTP
+ Bool_t modeFET3 = kTRUE;
+ if (GetFetMode(cfg) != 3) {
+ printf("FET is not in mode 3. Only PHYSICS events will be analysed (noisy channels)\n");
+ modeFET3 = kFALSE;
+ }
- // loop over regional structures
- Int_t nReg = darcHeader->GetRegHeaderEntries();
- for(Int_t iReg = 0; iReg < nReg; ++iReg){ //REG loop
+ // All 5 global cards are controlled by the Mts proxy
+ if (cfg.GetGlobalMasks()->GetGlobalCrateEnable() != 0x1F) {
+ printf("The MTS proxy does not control all global cards\n");
+ return -1;
+ }
- if (printLevel) printf("RegionalId %d\n", iReg);
+ // The global cards are ON (active on the global inputs)
+ if (!cfg.GetGlobalMasks()->GetMasksOn()) {
+ printf("Global masks are not ON\n");
+ return -1;
+ }
+
+ // make sure to catch the "rare" calib events (1 every 50s in physics)
+ const Char_t* tableSOD[] = {"ALL", "yes", "CAL", "all", NULL, NULL};
+ monitorDeclareTable(const_cast<char**>(tableSOD));
- regHeader = darcHeader->GetRegHeaderEntry(iReg);
+ status = monitorSetDataSource(inputFile);
+ if (status) {
+ cerr << "ERROR : monitorSetDataSource status (hex) = " << hex << status
+ << " " << monitorDecodeError(status) << endl;
+ return -1;
+ }
+ status = monitorDeclareMp("MUON Trigger monitoring");
+ if (status) {
+ cerr << "ERROR : monitorDeclareMp status (hex) = " << hex << status
+ << " " << monitorDecodeError(status) << endl;
+ return -1;
+ }
+
+ /* define wait event timeout - 1s max */
+ monitorSetNowait();
+ monitorSetNoWaitNetworkTimeout(1000);
+
+ cout << "MUONTRGda : Reading data from file " << inputFile <<endl;
+
+ UInt_t *globalInput = new UInt_t[4];
+ Bool_t doUpdate = kFALSE;
+ Int_t runNumber = 0;
+ Int_t nEvents = 0;
+
+ while(1)
+ {
+ if (nEvents >= cfg.GetMaxEvents()) break;
+ if (cfg.GetPrintLevel()) {
+ if (nEvents && nEvents % 1000 == 0)
+ cout<<"Cumulated events " << nEvents << endl;
+ }
+ // check shutdown condition
+ if (daqDA_checkShutdown())
+ break;
+
+ // Skip Events if needed
+ while (cfg.GetSkipEvents()) {
+ status = monitorGetEventDynamic(&event);
+ cfg.DecSkipEvents();
+ }
+
+ // starts reading
+ status = monitorGetEventDynamic(&event);
+ if (status < 0) {
+ cout << "MUONTRGda : EOF found" << endl;
+ break;
+ }
+
+ nDateEvents++;
+
+ // decoding rawdata headers
+ AliRawReader *rawReader = new AliRawReaderDate(event);
+
+ Int_t eventType = rawReader->GetType();
+ runNumber = rawReader->GetRunNumber();
+
+ // L1Swc1
+ // CALIBRATION_EVENT
+ // SYSTEM_SOFTWARE_TRIGGER_EVENT
+ // DETECTOR_SOFTWARE_TRIGGER_EVENT
+ cfg.SetAlgoNoisyInput(kFALSE);
+ cfg.SetAlgoDeadcInput(kFALSE);
+ if (eventType == PHYSICS_EVENT) {
+ cfg.SetAlgoNoisyInput(kTRUE);
+ doUpdate = kTRUE;
+ cfg.IncNoiseEvent();
+ } else if (modeFET3 && eventType == CALIBRATION_EVENT) {
+ cfg.SetAlgoDeadcInput(kTRUE);
+ doUpdate = kTRUE;
+ cfg.IncDeadcEvent();
+ } else {
+ continue;
+ }
+
+ nEvents++;
+ if (cfg.GetPrintLevel() == 2) printf("\nEvent # %d\n",nEvents);
+
+ // decoding MUON payload
+ AliMUONVRawStreamTrigger *rawStream = 0x0;
+ if (cfg.UseFastDecoder()) {
+ rawStream = new AliMUONRawStreamTriggerHP(rawReader);
+ } else {
+ rawStream = new AliMUONRawStreamTrigger(rawReader);
+ }
- // loop over local structures
- Int_t nLocal = regHeader->GetLocalEntries();
- for(Int_t iLocal = 0; iLocal < nLocal; ++iLocal) {
+ // ... without warnings from the decoder !!!
+ if (!cfg.WithWarnings())
+ rawStream->DisableWarnings();
- localStruct = regHeader->GetLocalEntry(iLocal);
+ // loops over DDL
+ while((status = rawStream->NextDDL())) {
- Int_t loStripX = (Int_t)localStruct->GetXPos();
- Int_t loStripY = (Int_t)localStruct->GetYPos();
- Int_t loDev = (Int_t)localStruct->GetXDev();
-
- if (printLevel) printf("Index %d, XPos: %d, YPos: %d Dev: %d\n",
- localStruct->GetId(), loStripX, loStripY, loDev);
+ if (cfg.GetPrintLevel() == 2) printf("iDDL %d\n", rawStream->GetDDL());
- if (printLevel) printf("X pattern %x %x %x %x, Y pattern %x %x %x %x\n",
- localStruct->GetX1(), localStruct->GetX2(),localStruct->GetX3(),localStruct->GetX4(),
- localStruct->GetY1(), localStruct->GetY2(),localStruct->GetY3(),localStruct->GetY4());
-
- } // iLocal
- } // iReg
- } // NextDDL
+ if (cfg.UseFastDecoder()) {
+ darcHeaderHP = static_cast<AliMUONRawStreamTriggerHP*>(rawStream)->GetHeaders();
+ if (cfg.GetPrintLevel() == 2) printf("Global output (fast decoder) %x\n", (Int_t)darcHeaderHP->GetGlobalOutput());
+ for (Int_t ig = 0; ig < cfg.GetGlobalInputs(); ig++) {
+ globalInput[ig] = darcHeaderHP->GetGlobalInput(ig);
+ }
+ // loop over regional structure
+ Int_t nReg = (Int_t)static_cast<AliMUONRawStreamTriggerHP*>(rawStream)->GetRegionalHeaderCount();
+ for(Int_t iReg = 0; iReg < nReg; iReg++) {
+ // loop over local structures
+ Int_t nLoc = (Int_t)static_cast<AliMUONRawStreamTriggerHP*>(rawStream)->GetLocalStructCount(iReg);
+ for(Int_t iLoc = 0; iLoc < nLoc; iLoc++) {
+ localStructHP = static_cast<AliMUONRawStreamTriggerHP*>(rawStream)->GetLocalStruct(iReg, iLoc);
+ localBoardId = cfg.GetTriggerIO()->LocalBoardId(rawStream->GetDDL(),iReg,localStructHP->GetId());
+ if (localBoardId > 0) {
+ localStructHP->GetXPattern(xPattern);
+ localStructHP->GetYPattern(yPattern);
+ MakePattern(cfg,localBoardId,xPattern,yPattern);
+ }
+ }
+ }
+ } else {
+ ddlTrigger = rawStream->GetDDLTrigger();
+ darcHeader = ddlTrigger->GetDarcHeader();
+ if (cfg.GetPrintLevel() == 2) printf("Global output %x\n", (Int_t)darcHeader->GetGlobalOutput());
+ globalInput = darcHeader->GetGlobalInput();
+ // loop over regional structure
+ Int_t nReg = darcHeader->GetRegHeaderEntries();
+ for(Int_t iReg = 0; iReg < nReg; iReg++) {
+ regHeader = darcHeader->GetRegHeaderEntry(iReg);
+ // loop over local structures
+ Int_t nLoc = regHeader->GetLocalEntries();
+ for(Int_t iLoc = 0; iLoc < nLoc; iLoc++) {
+ localStruct = regHeader->GetLocalEntry(iLoc);
+ localBoardId = cfg.GetTriggerIO()->LocalBoardId(rawStream->GetDDL(),iReg,localStruct->GetId());
+ if (localBoardId > 0) {
+ localStruct->GetXPattern(xPattern);
+ localStruct->GetYPattern(yPattern);
+ MakePattern(cfg,localBoardId,xPattern,yPattern);
+ }
+ }
+ }
+ }
+ if (rawStream->GetDDL() == 0) {
+ StoreGlobalInput(cfg,globalInput);
+ }
+
+ } // NextDDL
- delete rawReader;
- delete rawStream;
+ delete rawReader;
+ delete rawStream;
- } // while (1)
+ } // while (1)
+ // update configuration files ifrequested event types were found
+ if (doUpdate) {
+ if (cfg.GetDAMode() > 0) UpdateGlobalMasks(cfg);
+ if (cfg.GetDAMode() > 1) MakePatternStore(cfg);
+ }
- timers.Stop();
+ timers.Stop();
- cout << "MUONTRKda : Run number : " << runNumber << endl;
- cout << "MUONTRKda : Flat ASCII file generated : " << flatOutputFile << endl;
- // cout << "MUONTRKda : Histo file generated : " << histoFileName << endl;
- cout << "MUONTRKda : Nb of DATE events = " << nDateEvents << endl;
- cout << "MUONTRKda : Nb of events used = " << nEvents << endl;
+ cout << "MUONTRGda: DA enable: " << cfg.GetDAFlag() << endl;
+ cout << "MUONTRGda: Run number: " << runNumber << endl;
+ cout << "MUONTRGda: Nb of DATE events: " << nDateEvents << endl;
+ cout << "MUONTRGda: Nb of events used: " << nEvents << endl;
+ cout << "MUONTRGda: Nb of events used (noise): " << cfg.GetEventsN() << endl;
+ cout << "MUONTRGda: Nb of events used (deadc): " << cfg.GetEventsD() << endl;
+ cout << "MUONTRGda: Minumum nr of events for rate calculation: " << cfg.GetMinEvents() << endl;
+ cout << "MUONTRGda: Maximum nr of analyzed events: " << cfg.GetMaxEvents() << endl;
+ cout << "MUONTRGda: Skip events from start: " << cfg.GetSkipEvents() << endl;
+ cout << "MUONTRGda: Threshold for global noisy inputs: " << 100*cfg.GetThrN() << "%" << endl;
+ cout << "MUONTRGda: Threshold for global dead inputs: " << 100*cfg.GetThrD() << "%" << endl;
+ cout << "MUONTRGda: Threshold for local noisy inputs: " << 100*cfg.GetThrLocN() << "%" << endl;
+ cout << "MUONTRGda: Threshold for local dead inputs: " << 100*cfg.GetThrLocD() << "%" << endl;
+ cout << "MUONTRGda: Print level: " << cfg.GetPrintLevel() << endl;
+ cout << "MUONTRGda: Show decoder warnings: " << cfg.WithWarnings() << endl;
+ cout << "MUONTRGda: Use the fast decoder: " << cfg.UseFastDecoder() << endl;
+ cout << "MUONTRGda: DA mode (1=GLOBAL, 2=GLOBAL+LOCAL): " << cfg.GetDAMode() << endl;
- printf("Execution time : R:%7.2fs C:%7.2fs\n", timers.RealTime(), timers.CpuTime());
+ printf("MUONTRGda: Execution time : R:%7.2fs C:%7.2fs\n", timers.RealTime(), timers.CpuTime());
+
+ return status;
- return status;
}
+
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff