// this feature is not available for all detectors and that merging is not //
// possible, when digits are created directly from hits. //
// //
-// Backgound events can be merged by calling //
+// Background events can be merged by calling //
// //
// sim.MergeWith("background/galice.root", 2); //
// //
// The first argument is the file name of the background galice file. The //
// second argument is the number of signal events per background event. //
-// The default value for this is 1. MergeWith can be called several times //
-// to merge more than two event streams. It is assumed that the sdigits //
-// were already produced for the background events. //
+// By default this number is calculated from the number of available //
+// background events. MergeWith can be called several times to merge more //
+// than two event streams. It is assumed that the sdigits were already //
+// produced for the background events. //
// //
-// The methods RunSimulation, RunSDigitization, RunDigitization and //
-// RunHitsDigitization can be used to run only parts of the full simulation //
-// chain. //
+// The output of raw data can be switched on by calling //
+// //
+// sim.SetWriteRawData("MUON"); // write raw data for MUON //
+// //
+// The default output format of the raw data are DDL files. They are //
+// converted to a DATE file, if a file name is given as second argument. //
+// For this conversion the program "dateStream" is required. If the file //
+// name has the extension ".root", the DATE file is converted to a root //
+// file. The program "alimdc" is used for this purpose. For the conversion //
+// to DATE and root format the two conversion programs have to be installed. //
+// Only the raw data in the final format is kept if the third argument is //
+// kTRUE. //
+// //
+// The methods RunSimulation, RunSDigitization, RunDigitization, //
+// RunHitsDigitization and WriteRawData can be used to run only parts of //
+// the full simulation chain. The creation of raw data DDL files and their //
+// conversion to the DATE or root format can be run directly by calling //
+// the methods WriteRawFiles, ConvertRawFilesToDate and ConvertDateToRoot. //
+// //
+// The default number of events per file, which is usually set in the //
+// config file, can be changed for individual detectors and data types //
+// by calling //
+// //
+// sim.SetEventsPerFile("PHOS", "Reconstructed Points", 3); //
+// //
+// The first argument is the detector, the second one the data type and the //
+// last one the number of events per file. Valid data types are "Hits", //
+// "Summable Digits", "Digits", "Reconstructed Points" and "Tracks". //
+// The number of events per file has to be set before the simulation of //
+// hits. Otherwise it has no effect. //
// //
///////////////////////////////////////////////////////////////////////////////
-
-#include "AliSimulation.h"
-#include "AliRunLoader.h"
-#include "AliRun.h"
-#include "AliModule.h"
+#include <TGeoManager.h>
+#include <TObjString.h>
+#include <TStopwatch.h>
+#include <TSystem.h>
+#include <TFile.h>
+
+#include "AliCDBStorage.h"
+#include "AliCDBEntry.h"
+#include "AliCDBManager.h"
+#include "AliAlignObj.h"
+#include "AliCentralTrigger.h"
+#include "AliDAQConfig.h"
+#include "AliDigitizer.h"
#include "AliGenerator.h"
+#include "AliLog.h"
+#include "AliModule.h"
+#include "AliRun.h"
#include "AliRunDigitizer.h"
-#include "AliDigitizer.h"
-#include <TObjString.h>
-
+#include "AliRunLoader.h"
+#include "AliSimulation.h"
+#include "AliVertexGenFile.h"
+#include "AliCentralTrigger.h"
ClassImp(AliSimulation)
fRunSimulation(kTRUE),
fMakeSDigits("ALL"),
fMakeDigits("ALL"),
+ fMakeTrigger(""),
fMakeDigitsFromHits(""),
+ fWriteRawData(""),
+ fRawDataFileName(""),
+ fDeleteIntermediateFiles(kFALSE),
fStopOnError(kFALSE),
fNEvents(1),
fConfigFileName(configFileName),
fGAliceFileName("galice.root"),
+ fEventsPerFile(),
fBkgrdFileNames(NULL),
- fRegionOfInterest(kTRUE)
+ fAlignObjArray(NULL),
+ fUseBkgrdVertex(kTRUE),
+ fRegionOfInterest(kFALSE)
{
// create simulation object with default parameters
+ SetGAliceFile("galice.root");
}
//_____________________________________________________________________________
fRunSimulation(sim.fRunSimulation),
fMakeSDigits(sim.fMakeSDigits),
fMakeDigits(sim.fMakeDigits),
+ fMakeTrigger(sim.fMakeTrigger),
fMakeDigitsFromHits(sim.fMakeDigitsFromHits),
+ fWriteRawData(sim.fWriteRawData),
+ fRawDataFileName(""),
+ fDeleteIntermediateFiles(kFALSE),
fStopOnError(sim.fStopOnError),
fNEvents(sim.fNEvents),
fConfigFileName(sim.fConfigFileName),
fGAliceFileName(sim.fGAliceFileName),
+ fEventsPerFile(),
fBkgrdFileNames(NULL),
+ fAlignObjArray(NULL),
+ fUseBkgrdVertex(sim.fUseBkgrdVertex),
fRegionOfInterest(sim.fRegionOfInterest)
{
// copy constructor
+ for (Int_t i = 0; i < sim.fEventsPerFile.GetEntriesFast(); i++) {
+ if (!sim.fEventsPerFile[i]) continue;
+ fEventsPerFile.Add(sim.fEventsPerFile[i]->Clone());
+ }
+
fBkgrdFileNames = new TObjArray;
for (Int_t i = 0; i < sim.fBkgrdFileNames->GetEntriesFast(); i++) {
if (!sim.fBkgrdFileNames->At(i)) continue;
{
// clean up
+ fEventsPerFile.Delete();
+
if (fBkgrdFileNames) {
fBkgrdFileNames->Delete();
delete fBkgrdFileNames;
fConfigFileName = fileName;
}
+//_____________________________________________________________________________
+void AliSimulation::SetGAliceFile(const char* fileName)
+{
+// set the name of the galice file
+// the path is converted to an absolute one if it is relative
+
+ fGAliceFileName = fileName;
+ if (!gSystem->IsAbsoluteFileName(fGAliceFileName)) {
+ char* absFileName = gSystem->ConcatFileName(gSystem->WorkingDirectory(),
+ fGAliceFileName);
+ fGAliceFileName = absFileName;
+ delete[] absFileName;
+ }
+
+ AliDebug(2, Form("galice file name set to %s", fileName));
+}
+
+//_____________________________________________________________________________
+void AliSimulation::SetEventsPerFile(const char* detector, const char* type,
+ Int_t nEvents)
+{
+// set the number of events per file for the given detector and data type
+// ("Hits", "Summable Digits", "Digits", "Reconstructed Points" or "Tracks")
+
+ TNamed* obj = new TNamed(detector, type);
+ obj->SetUniqueID(nEvents);
+ fEventsPerFile.Add(obj);
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::ApplyAlignObjsToGeom(const char* fileName, const char* clArrayName)
+{
+ // read collection of alignment objects (AliAlignObj derived) saved
+ // in the TClonesArray ClArrayName in the file fileName and apply
+ // them to the TGeo geometry passed as argument
+ //
+
+ TFile* inFile = TFile::Open(fileName,"READ");
+ if (!inFile || !inFile->IsOpen()) {
+ AliErrorClass(Form("Could not open file %s !",fileName));
+ return kFALSE;
+ }
+
+ TClonesArray* alObjArray = ((TClonesArray*) inFile->Get(clArrayName));
+ inFile->Close();
+ if (!alObjArray) {
+ AliErrorClass(Form("Could not get array (%s) from file (%s) !",clArrayName,fileName));
+ return kFALSE;
+ }
+
+ return gAlice->ApplyAlignObjsToGeom(alObjArray);
+
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::ApplyAlignObjsToGeom(AliCDBParam* param, AliCDBId& Id)
+{
+ // read collection of alignment objects (AliAlignObj derived) saved
+ // in the TClonesArray ClArrayName in the AliCDBEntry identified by
+ // param (to get the AliCDBStorage) and Id; apply the alignment objects
+ // to the TGeo geometry passed as argument
+ //
+
+ AliCDBStorage* storage = AliCDBManager::Instance()->GetStorage(param);
+ AliCDBEntry* entry = storage->Get(Id);
+ TClonesArray* AlObjArray = ((TClonesArray*) entry->GetObject());
+
+ return gAlice->ApplyAlignObjsToGeom(AlObjArray);
+
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::ApplyAlignObjsToGeom(const char* uri, const char* path, Int_t runnum, Int_t version, Int_t sversion)
+{
+ // read collection of alignment objects (AliAlignObj derived) saved
+ // in the TClonesArray ClArrayName in the AliCDBEntry identified by
+ // param (to get the AliCDBStorage) and Id; apply the alignment objects
+ // to the TGeo geometry passed as argument
+ //
+
+ AliCDBParam* param = AliCDBManager::Instance()->CreateParameter(uri);
+ AliCDBId id(path, runnum, runnum, version, sversion);
+
+ return ApplyAlignObjsToGeom(param, id);
+
+}
+
//_____________________________________________________________________________
void AliSimulation::MergeWith(const char* fileName, Int_t nSignalPerBkgrd)
{
{
// run the generation, simulation and digitization
+ // First check if we have any CDB storage set, because it is used
+ // to retrieve the calibration and alignment constants
+
+ AliCDBManager* man = AliCDBManager::Instance();
+ if (!man->IsDefaultStorageSet())
+ {
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ AliWarning("No default CDB storage set, so I will use $ALICE_ROOT");
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ man->SetDefaultStorage("local://$ALICE_ROOT");
+ }
+
if (nEvents > 0) fNEvents = nEvents;
// generation and simulation -> hits
// hits -> digits
if (!fMakeDigitsFromHits.IsNull()) {
if (fBkgrdFileNames && (fBkgrdFileNames->GetEntriesFast() > 0)) {
- Warning("Run", "Merging and direct creation of digits from hits "
- "was selected for some detectors. "
- "No merging will be done for the following detectors: %s",
- fMakeDigitsFromHits.Data());
+ AliWarning(Form("Merging and direct creation of digits from hits "
+ "was selected for some detectors. "
+ "No merging will be done for the following detectors: %s",
+ fMakeDigitsFromHits.Data()));
}
if (!RunHitsDigitization(fMakeDigitsFromHits)) {
if (fStopOnError) return kFALSE;
}
}
+ // digits -> trigger
+ if (!fMakeTrigger.IsNull()) {
+ if (!RunTrigger(fMakeTrigger)) {
+ if (fStopOnError) return kFALSE;
+ }
+ }
+
+ // digits -> raw data
+ if (!fWriteRawData.IsNull()) {
+ if (!WriteRawData(fWriteRawData, fRawDataFileName,
+ fDeleteIntermediateFiles)) {
+ if (fStopOnError) return kFALSE;
+ }
+ }
+
return kTRUE;
}
+//_____________________________________________________________________________
+Bool_t AliSimulation::RunTrigger(const char* descriptors)
+{
+ // run the trigger
+
+ TStopwatch stopwatch;
+ stopwatch.Start();
+
+ AliRunLoader* runLoader = LoadRun("READ");
+ if (!runLoader) return kFALSE;
+ TString des = descriptors;
+ // Load Descriptors
+ AliCentralTrigger* aCTP = new AliCentralTrigger( des );
+
+ // digits -> trigger
+ if( !aCTP->RunTrigger( runLoader ) ) {
+ if (fStopOnError) {
+ delete aCTP;
+ return kFALSE;
+ }
+ }
+
+/*
+ // Process each event
+ for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
+ AliInfo(Form("processing event %d", iEvent));
+ runLoader->GetEvent(iEvent);
+
+ TObjArray* detArray = runLoader->GetAliRun()->Detectors();
+ for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
+ AliModule* det = (AliModule*) detArray->At(iDet);
+ if (!det || !det->IsActive()) continue;
+ if (IsSelected(det->GetName(), detStr)) {
+ AliInfo(Form("triggering from digits for %s", det->GetName()));
+
+ // AliLoader* loader = fLoader[iDet];
+ // loader->LoadDigits("read");
+ // TTree* digitsTree = loader->TreeD();
+ // det->Trigger( digitsTree );
+ // or
+ AliTriggerDetector* tdet = det->CreateTriggerDetector();
+ TObjArray* detInp = dtrg->GetTriggerInputs();
+ for( Int_t i=0; i<detInp->GetEntriesFast(); i++ )
+ fInputs.AddLast( detInp->At(i) );
+
+ AliInfo(Form("Execution time for %s: R:%.2fs C:%.2fs",
+ det->GetName(),stopwatchDet.RealTime(),stopwatchDet.CpuTime()));
+ }
+ }
+
+ if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
+ AliError(Form("the following detectors were not found: %s",
+ detStr.Data()));
+ if (fStopOnError) {
+ delete centralTP;
+ return kFALSE;
+ }
+ }
+
+ // Check trigger conditions
+ centralTP->TriggerConditions();
+
+ // Write trigger ????
+ centralTP->Write();
+
+ } */
+
+ AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
+ stopwatch.RealTime(),stopwatch.CpuTime()));
+
+ delete aCTP;
+ delete runLoader;
+
+ return kTRUE;
+}
+
+
+
//_____________________________________________________________________________
Bool_t AliSimulation::RunSimulation(Int_t nEvents)
{
stopwatch.Start();
if (!gAlice) {
- Error("RunSimulation", "no gAlice object. Restart aliroot and try again.");
+ AliError("no gAlice object. Restart aliroot and try again.");
return kFALSE;
}
if (gAlice->Modules()->GetEntries() > 0) {
- Error("RunSimulation",
- "gAlice was already run. Restart aliroot and try again.");
+ AliError("gAlice was already run. Restart aliroot and try again.");
return kFALSE;
}
- Info("RunSimulation", "initializing gAlice with config file %s",
- fConfigFileName.Data());
- gAlice->Init(fConfigFileName.Data());
+ AliInfo(Form("initializing gAlice with config file %s",
+ fConfigFileName.Data()));
+ StdoutToAliInfo(StderrToAliError(
+ gAlice->Init(fConfigFileName.Data());
+ ););
+
+ // Check if the array with alignment objects was
+ // provided by the user. If yes, apply the objects
+ // to the present TGeo geometry
+ if (fAlignObjArray) {
+ if (gGeoManager && gGeoManager->IsClosed()) {
+ if (gAlice->ApplyAlignObjsToGeom(fAlignObjArray) == kFALSE) {
+ AliError("The application of misalignment failed! Restart aliroot and try again. ");
+ return kFALSE;
+ }
+ }
+ else {
+ AliError("Can't apply the misalignment! gGeoManager doesn't exist or it is still opened!");
+ return kFALSE;
+ }
+ }
+
+ // Export TGeo geometry
+ if (gGeoManager) gGeoManager->Export("geometry.root");
+
AliRunLoader* runLoader = gAlice->GetRunLoader();
if (!runLoader) {
- Error("RunSimulation", "gAlice has no run loader object. "
- "Check your config file: %s", fConfigFileName.Data());
+ AliError(Form("gAlice has no run loader object. "
+ "Check your config file: %s", fConfigFileName.Data()));
return kFALSE;
}
- fGAliceFileName = runLoader->GetFileName();
+ SetGAliceFile(runLoader->GetFileName());
+
+ if (!gAlice->Generator()) {
+ AliError(Form("gAlice has no generator object. "
+ "Check your config file: %s", fConfigFileName.Data()));
+ return kFALSE;
+ }
+ if (nEvents <= 0) nEvents = fNEvents;
+
+ // get vertex from background file in case of merging
+ if (fUseBkgrdVertex &&
+ fBkgrdFileNames && (fBkgrdFileNames->GetEntriesFast() > 0)) {
+ Int_t signalPerBkgrd = GetNSignalPerBkgrd(nEvents);
+ const char* fileName = ((TObjString*)
+ (fBkgrdFileNames->At(0)))->GetName();
+ AliInfo(Form("The vertex will be taken from the background "
+ "file %s with nSignalPerBackground = %d",
+ fileName, signalPerBkgrd));
+ AliVertexGenFile* vtxGen = new AliVertexGenFile(fileName, signalPerBkgrd);
+ gAlice->Generator()->SetVertexGenerator(vtxGen);
+ }
if (!fRunSimulation) {
- if (!gAlice->Generator()) {
- Error("RunSimulation", "gAlice has no generator object. "
- "Check your config file: %s", fConfigFileName.Data());
- return kFALSE;
- }
gAlice->Generator()->SetTrackingFlag(0);
}
- Info("RunSimulation", "running gAlice");
- if (nEvents <= 0) nEvents = fNEvents;
- gAlice->Run(nEvents);
+ // set the number of events per file for given detectors and data types
+ for (Int_t i = 0; i < fEventsPerFile.GetEntriesFast(); i++) {
+ if (!fEventsPerFile[i]) continue;
+ const char* detName = fEventsPerFile[i]->GetName();
+ const char* typeName = fEventsPerFile[i]->GetTitle();
+ TString loaderName(detName);
+ loaderName += "Loader";
+ AliLoader* loader = runLoader->GetLoader(loaderName);
+ if (!loader) {
+ AliError(Form("RunSimulation", "no loader for %s found\n"
+ "Number of events per file not set for %s %s",
+ detName, typeName, detName));
+ continue;
+ }
+ AliDataLoader* dataLoader =
+ loader->GetDataLoader(typeName);
+ if (!dataLoader) {
+ AliError(Form("no data loader for %s found\n"
+ "Number of events per file not set for %s %s",
+ typeName, detName, typeName));
+ continue;
+ }
+ dataLoader->SetNumberOfEventsPerFile(fEventsPerFile[i]->GetUniqueID());
+ AliDebug(1, Form("number of events per file set to %d for %s %s",
+ fEventsPerFile[i]->GetUniqueID(), detName, typeName));
+ }
+
+ AliInfo("running gAlice");
+ StdoutToAliInfo(StderrToAliError(
+ gAlice->Run(nEvents);
+ ););
delete runLoader;
- Info("RunSimulation", "execution time:");
- stopwatch.Print();
+ AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
+ stopwatch.RealTime(),stopwatch.CpuTime()));
return kTRUE;
}
AliModule* det = (AliModule*) detArray->At(iDet);
if (!det || !det->IsActive()) continue;
if (IsSelected(det->GetName(), detStr)) {
- Info("RunSDigitization", "creating summable digits for %s",
- det->GetName());
+ AliInfo(Form("creating summable digits for %s", det->GetName()));
+ TStopwatch stopwatchDet;
+ stopwatchDet.Start();
det->Hits2SDigits();
+ AliInfo(Form("Execution time for %s: R:%.2fs C:%.2fs",
+ det->GetName(),stopwatchDet.RealTime(),stopwatchDet.CpuTime()));
}
}
if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
- Error("RunSDigitization", "the following detectors were not found: %s",
- detStr.Data());
+ AliError(Form("the following detectors were not found: %s",
+ detStr.Data()));
if (fStopOnError) return kFALSE;
}
delete runLoader;
- Info("RunSDigitization", "execution time:");
- stopwatch.Print();
+ AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
+ stopwatch.RealTime(),stopwatch.CpuTime()));
return kTRUE;
}
gAlice = NULL;
Int_t nStreams = 1;
- Int_t signalPerBkgrd = 1;
- if (fBkgrdFileNames) {
- nStreams = fBkgrdFileNames->GetEntriesFast() + 1;
- if (nStreams > 1) signalPerBkgrd = fBkgrdFileNames->At(0)->GetUniqueID();
- }
+ if (fBkgrdFileNames) nStreams = fBkgrdFileNames->GetEntriesFast() + 1;
+ Int_t signalPerBkgrd = GetNSignalPerBkgrd();
AliRunDigitizer* manager = new AliRunDigitizer(nStreams, signalPerBkgrd);
manager->SetInputStream(0, fGAliceFileName.Data());
for (Int_t iStream = 1; iStream < nStreams; iStream++) {
!IsSelected(det->GetName(), detExcl)) {
AliDigitizer* digitizer = det->CreateDigitizer(manager);
if (!digitizer) {
- Error("RunDigitization", "no digitizer for %s", det->GetName());
+ AliError(Form("no digitizer for %s", det->GetName()));
if (fStopOnError) return kFALSE;
} else {
digitizer->SetRegionOfInterest(fRegionOfInterest);
}
if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
- Error("RunDigitization", "the following detectors were not found: %s",
- detStr.Data());
+ AliError(Form("the following detectors were not found: %s",
+ detStr.Data()));
if (fStopOnError) return kFALSE;
}
if (!manager->GetListOfTasks()->IsEmpty()) {
- Info("RunDigitization", "executing digitization");
+ AliInfo("executing digitization");
manager->Exec("");
}
delete manager;
- Info("RunDigitization", "execution time:");
- stopwatch.Print();
-
+ AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
+ stopwatch.RealTime(),stopwatch.CpuTime()));
+
return kTRUE;
}
TStopwatch stopwatch;
stopwatch.Start();
- AliRunLoader* runLoader = LoadRun();
+ AliRunLoader* runLoader = LoadRun("READ");
if (!runLoader) return kFALSE;
TString detStr = detectors;
AliModule* det = (AliModule*) detArray->At(iDet);
if (!det || !det->IsActive()) continue;
if (IsSelected(det->GetName(), detStr)) {
- Info("RunHitsDigitization", "creating digits from hits for %s",
- det->GetName());
+ AliInfo(Form("creating digits from hits for %s", det->GetName()));
det->Hits2Digits();
}
}
if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
- Error("RunHitsDigitization", "the following detectors were not found: %s",
- detStr.Data());
+ AliError(Form("the following detectors were not found: %s",
+ detStr.Data()));
if (fStopOnError) return kFALSE;
}
delete runLoader;
+ //PH Temporary fix to avoid interference with the PHOS loder/getter
+ //PH The problem has to be solved in more general way 09/06/05
- Info("RunHitsDigitization", "execution time:");
- stopwatch.Print();
+ AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
+ stopwatch.RealTime(),stopwatch.CpuTime()));
return kTRUE;
}
+//_____________________________________________________________________________
+Bool_t AliSimulation::WriteRawData(const char* detectors,
+ const char* fileName,
+ Bool_t deleteIntermediateFiles)
+{
+// convert the digits to raw data
+// First DDL raw data files for the given detectors are created.
+// If a file name is given, the DDL files are then converted to a DATE file.
+// If deleteIntermediateFiles is true, the DDL raw files are deleted
+// afterwards.
+// If the file name has the extension ".root", the DATE file is converted
+// to a root file.
+// If deleteIntermediateFiles is true, the DATE file is deleted afterwards.
+
+ TStopwatch stopwatch;
+ stopwatch.Start();
+
+ if (!WriteRawFiles(detectors)) {
+ if (fStopOnError) return kFALSE;
+ }
+
+ TString dateFileName(fileName);
+ if (!dateFileName.IsNull()) {
+ Bool_t rootOutput = dateFileName.EndsWith(".root");
+ if (rootOutput) dateFileName += ".date";
+ if (!ConvertRawFilesToDate(dateFileName)) {
+ if (fStopOnError) return kFALSE;
+ }
+ if (deleteIntermediateFiles) {
+ AliRunLoader* runLoader = LoadRun("READ");
+ if (runLoader) for (Int_t iEvent = 0;
+ iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
+ char command[256];
+ sprintf(command, "rm -r raw%d", iEvent);
+ gSystem->Exec(command);
+ }
+ }
+
+ if (rootOutput) {
+ if (!ConvertDateToRoot(dateFileName, fileName)) {
+ if (fStopOnError) return kFALSE;
+ }
+ if (deleteIntermediateFiles) {
+ gSystem->Unlink(dateFileName);
+ }
+ }
+ }
+
+ AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
+ stopwatch.RealTime(),stopwatch.CpuTime()));
+
+ return kTRUE;
+}
//_____________________________________________________________________________
-AliRunLoader* AliSimulation::LoadRun() const
+Bool_t AliSimulation::WriteRawFiles(const char* detectors)
+{
+// convert the digits to raw data DDL files
+
+ AliRunLoader* runLoader = LoadRun("READ");
+ if (!runLoader) return kFALSE;
+
+ // write raw data to DDL files
+ for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
+ AliInfo(Form("processing event %d", iEvent));
+ runLoader->GetEvent(iEvent);
+ TString baseDir = gSystem->WorkingDirectory();
+ char dirName[256];
+ sprintf(dirName, "raw%d", iEvent);
+ gSystem->MakeDirectory(dirName);
+ if (!gSystem->ChangeDirectory(dirName)) {
+ AliError(Form("couldn't change to directory %s", dirName));
+ if (fStopOnError) return kFALSE; else continue;
+ }
+
+ TString detStr = detectors;
+ TObjArray* detArray = runLoader->GetAliRun()->Detectors();
+ for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
+ AliModule* det = (AliModule*) detArray->At(iDet);
+ if (!det || !det->IsActive()) continue;
+ if (IsSelected(det->GetName(), detStr)) {
+ AliInfo(Form("creating raw data from digits for %s", det->GetName()));
+ det->Digits2Raw();
+ }
+ }
+
+ gSystem->ChangeDirectory(baseDir);
+ if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
+ AliError(Form("the following detectors were not found: %s",
+ detStr.Data()));
+ if (fStopOnError) return kFALSE;
+ }
+ }
+
+ delete runLoader;
+ return kTRUE;
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::ConvertRawFilesToDate(const char* dateFileName)
+{
+// convert raw data DDL files to a DATE file with the program "dateStream"
+
+ char* path = gSystem->Which(gSystem->Getenv("PATH"), "dateStream");
+ if (!path) {
+ AliError("the program dateStream was not found");
+ if (fStopOnError) return kFALSE;
+ } else {
+ delete[] path;
+ }
+
+ AliRunLoader* runLoader = LoadRun("READ");
+ if (!runLoader) return kFALSE;
+
+ AliInfo(Form("converting raw data DDL files to DATE file %s", dateFileName));
+ char command[256];
+ sprintf(command, "dateStream -o %s -# %d -C",
+ dateFileName, runLoader->GetNumberOfEvents());
+ FILE* pipe = gSystem->OpenPipe(command, "w");
+
+ for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
+ fprintf(pipe, "GDC\n");
+ Float_t ldc = 0;
+ Int_t prevLDC = -1;
+
+ // loop over detectors and DDLs
+ for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
+ for (Int_t iDDL = 0; iDDL < kDetectorDDLs[iDet]; iDDL++) {
+
+ Int_t ddlID = 0x100*iDet + iDDL;
+ Int_t ldcID = Int_t(ldc + 0.0001);
+ ldc += kDetectorLDCs[iDet] / kDetectorDDLs[iDet];
+
+ char rawFileName[256];
+ sprintf(rawFileName, "raw%d/%s_%d.ddl",
+ iEvent, kDetectors[iDet], ddlID);
+
+ // check existence and size of raw data file
+ FILE* file = fopen(rawFileName, "rb");
+ if (!file) continue;
+ fseek(file, 0, SEEK_END);
+ unsigned long size = ftell(file);
+ fclose(file);
+ if (!size) continue;
+
+ if (ldcID != prevLDC) {
+ fprintf(pipe, " LDC Id %d\n", ldcID);
+ prevLDC = ldcID;
+ }
+ fprintf(pipe, " Equipment Id %d Payload %s\n", ddlID, rawFileName);
+ }
+ }
+ }
+
+ Int_t result = gSystem->ClosePipe(pipe);
+
+ delete runLoader;
+ return (result == 0);
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::ConvertDateToRoot(const char* dateFileName,
+ const char* rootFileName)
+{
+// convert a DATE file to a root file with the program "alimdc"
+
+ // ALIMDC setup
+ const Int_t kDBSize = 1000000000;
+ const Int_t kTagDBSize = 1000000000;
+ const Bool_t kFilter = kFALSE;
+ const Int_t kCompression = 1;
+
+ char* path = gSystem->Which(gSystem->Getenv("PATH"), "alimdc");
+ if (!path) {
+ AliError("the program alimdc was not found");
+ if (fStopOnError) return kFALSE;
+ } else {
+ delete[] path;
+ }
+
+ AliInfo(Form("converting DATE file %s to root file %s",
+ dateFileName, rootFileName));
+
+ gSystem->Exec("rm -rf /tmp/mdc1");
+ gSystem->Exec("rm -rf /tmp/mdc2");
+
+ gSystem->Exec("mkdir /tmp/mdc1");
+ gSystem->Exec("mkdir /tmp/mdc2");
+
+ char command[256];
+ sprintf(command, "alimdc %d %d %d %d %s",
+ kDBSize, kTagDBSize, kFilter, kCompression, dateFileName);
+ Int_t result = gSystem->Exec(command);
+ sprintf(command, "mv /tmp/mdc1/*.root %s", rootFileName);
+ gSystem->Exec(command);
+ gSystem->Exec("rm -rf /tmp/mdc1");
+ gSystem->Exec("rm -rf /tmp/mdc2");
+
+ return (result == 0);
+}
+
+
+//_____________________________________________________________________________
+AliRunLoader* AliSimulation::LoadRun(const char* mode) const
{
// delete existing run loaders, open a new one and load gAlice
while (AliRunLoader::GetRunLoader()) delete AliRunLoader::GetRunLoader();
AliRunLoader* runLoader =
AliRunLoader::Open(fGAliceFileName.Data(),
- AliConfig::fgkDefaultEventFolderName, "UPDATE");
+ AliConfig::GetDefaultEventFolderName(), mode);
if (!runLoader) {
- Error("LoadRun", "no run loader found in file %s",
- fGAliceFileName.Data());
+ AliError(Form("no run loader found in file %s", fGAliceFileName.Data()));
return NULL;
}
runLoader->LoadgAlice();
gAlice = runLoader->GetAliRun();
if (!gAlice) {
- Error("LoadRun", "no gAlice object found in file %s",
- fGAliceFileName.Data());
+ AliError(Form("no gAlice object found in file %s",
+ fGAliceFileName.Data()));
return NULL;
}
return runLoader;
}
+//_____________________________________________________________________________
+Int_t AliSimulation::GetNSignalPerBkgrd(Int_t nEvents) const
+{
+// get or calculate the number of signal events per background event
+
+ if (!fBkgrdFileNames) return 1;
+ Int_t nBkgrdFiles = fBkgrdFileNames->GetEntriesFast();
+ if (nBkgrdFiles == 0) return 1;
+
+ // get the number of signal events
+ if (nEvents <= 0) {
+ AliRunLoader* runLoader =
+ AliRunLoader::Open(fGAliceFileName.Data(), "SIGNAL");
+ if (!runLoader) return 1;
+ nEvents = runLoader->GetNumberOfEvents();
+ delete runLoader;
+ }
+
+ Int_t result = 0;
+ for (Int_t iBkgrdFile = 0; iBkgrdFile < nBkgrdFiles; iBkgrdFile++) {
+ // get the number of background events
+ const char* fileName = ((TObjString*)
+ (fBkgrdFileNames->At(iBkgrdFile)))->GetName();
+ AliRunLoader* runLoader =
+ AliRunLoader::Open(fileName, "BKGRD");
+ if (!runLoader) continue;
+ Int_t nBkgrdEvents = runLoader->GetNumberOfEvents();
+ delete runLoader;
+
+ // get or calculate the number of signal per background events
+ Int_t nSignalPerBkgrd = fBkgrdFileNames->At(iBkgrdFile)->GetUniqueID();
+ if (nSignalPerBkgrd <= 0) {
+ nSignalPerBkgrd = (nEvents-1) / nBkgrdEvents + 1;
+ } else if (result && (result != nSignalPerBkgrd)) {
+ AliInfo(Form("the number of signal events per background event "
+ "will be changed from %d to %d for stream %d",
+ nSignalPerBkgrd, result, iBkgrdFile+1));
+ nSignalPerBkgrd = result;
+ }
+
+ if (!result) result = nSignalPerBkgrd;
+ if (nSignalPerBkgrd * nBkgrdEvents < nEvents) {
+ AliWarning(Form("not enough background events (%d) for %d signal events "
+ "using %d signal per background events for stream %d",
+ nBkgrdEvents, nEvents, nSignalPerBkgrd, iBkgrdFile+1));
+ }
+ }
+
+ return result;
+}
+
//_____________________________________________________________________________
Bool_t AliSimulation::IsSelected(TString detName, TString& detectors) const
{