// //
///////////////////////////////////////////////////////////////////////////////
+#include <TVirtualMCApplication.h>
#include <TGeoManager.h>
#include <TObjString.h>
-#include <TStopwatch.h>
#include <TSystem.h>
#include <TFile.h>
+#include "AliCodeTimer.h"
#include "AliCDBStorage.h"
#include "AliCDBEntry.h"
#include "AliCDBManager.h"
+#include "AliGeomManager.h"
#include "AliAlignObj.h"
#include "AliCentralTrigger.h"
-#include "AliDAQConfig.h"
+#include "AliDAQ.h"
#include "AliDigitizer.h"
#include "AliGenerator.h"
#include "AliLog.h"
#include "AliVertexGenFile.h"
#include "AliCentralTrigger.h"
#include "AliCTPRawData.h"
+#include "AliRawReaderFile.h"
+#include "AliRawReaderRoot.h"
+#include "AliRawReaderDate.h"
+#include "AliESD.h"
+#include "AliHeader.h"
+#include "AliGenEventHeader.h"
+#include "AliMC.h"
+#include "AliHLTSimulation.h"
+#include "AliQADataMakerSteer.h"
+#include "AliSysInfo.h"
ClassImp(AliSimulation)
+AliSimulation *AliSimulation::fgInstance = 0;
+const char* AliSimulation::fgkDetectorName[AliSimulation::fgkNDetectors] = {"ITS", "TPC", "TRD", "TOF", "PHOS", "HMPID", "EMCAL", "MUON", "FMD", "ZDC", "PMD", "T0", "VZERO", "ACORDE", "HLT"};
//_____________________________________________________________________________
-AliSimulation::AliSimulation(const char* configFileName, const char* cdbUri,
+AliSimulation::AliSimulation(const char* configFileName,
const char* name, const char* title) :
TNamed(name, title),
fRunGeneration(kTRUE),
fRunSimulation(kTRUE),
fLoadAlignFromCDB(kTRUE),
- fLoadAlignData("ALL"),
+ fLoadAlObjsListOfDets("ALL"),
fMakeSDigits("ALL"),
fMakeDigits("ALL"),
fMakeTrigger(""),
fWriteRawData(""),
fRawDataFileName(""),
fDeleteIntermediateFiles(kFALSE),
+ fWriteSelRawData(kFALSE),
fStopOnError(kFALSE),
fNEvents(1),
fAlignObjArray(NULL),
fUseBkgrdVertex(kTRUE),
fRegionOfInterest(kFALSE),
- fCDBUri(cdbUri)
+ fCDBUri(""),
+ fSpecCDBUri(),
+ fRun(-1),
+ fSeed(0),
+ fInitCDBCalled(kFALSE),
+ fInitRunNumberCalled(kFALSE),
+ fSetRunNumberFromDataCalled(kFALSE),
+ fEmbeddingFlag(kFALSE),
+ fQADetectors("ALL"),
+ fQATasks("ALL"),
+ fRunQA(kTRUE),
+ fRunHLT("default")
{
// create simulation object with default parameters
-
+ fgInstance = this;
SetGAliceFile("galice.root");
+
+// for QA
+ for (Int_t iDet = 0; iDet < fgkNDetectors; iDet++)
+ fQACycles[iDet] = 999999;
}
//_____________________________________________________________________________
fRunGeneration(sim.fRunGeneration),
fRunSimulation(sim.fRunSimulation),
fLoadAlignFromCDB(sim.fLoadAlignFromCDB),
- fLoadAlignData(sim.fLoadAlignData),
+ fLoadAlObjsListOfDets(sim.fLoadAlObjsListOfDets),
fMakeSDigits(sim.fMakeSDigits),
fMakeDigits(sim.fMakeDigits),
fMakeTrigger(sim.fMakeTrigger),
fWriteRawData(sim.fWriteRawData),
fRawDataFileName(""),
fDeleteIntermediateFiles(kFALSE),
+ fWriteSelRawData(kFALSE),
fStopOnError(sim.fStopOnError),
fNEvents(sim.fNEvents),
fAlignObjArray(NULL),
fUseBkgrdVertex(sim.fUseBkgrdVertex),
fRegionOfInterest(sim.fRegionOfInterest),
- fCDBUri(sim.fCDBUri)
+ fCDBUri(sim.fCDBUri),
+ fSpecCDBUri(),
+ fRun(-1),
+ fSeed(0),
+ fInitCDBCalled(sim.fInitCDBCalled),
+ fInitRunNumberCalled(sim.fInitRunNumberCalled),
+ fSetRunNumberFromDataCalled(sim.fSetRunNumberFromDataCalled),
+ fEmbeddingFlag(sim.fEmbeddingFlag),
+ fQADetectors(sim.fQADetectors),
+ fQATasks(sim.fQATasks),
+ fRunQA(sim.fRunQA),
+ fRunHLT(sim.fRunHLT)
{
// copy constructor
if (!sim.fBkgrdFileNames->At(i)) continue;
fBkgrdFileNames->Add(sim.fBkgrdFileNames->At(i)->Clone());
}
+
+ for (Int_t i = 0; i < sim.fSpecCDBUri.GetEntriesFast(); i++) {
+ if (sim.fSpecCDBUri[i]) fSpecCDBUri.Add(sim.fSpecCDBUri[i]->Clone());
+ }
+ fgInstance = this;
+
+// for QA
+ for (Int_t iDet = 0; iDet < fgkNDetectors; iDet++)
+ fQACycles[iDet] = sim.fQACycles[iDet];
}
//_____________________________________________________________________________
fBkgrdFileNames->Delete();
delete fBkgrdFileNames;
}
+
+ fSpecCDBUri.Delete();
+ if (fgInstance==this) fgInstance = 0;
+
+ AliCodeTimer::Instance()->Print();
}
}
//_____________________________________________________________________________
-void AliSimulation::InitCDBStorage(const char* uri)
+void AliSimulation::InitCDB()
{
// activate a default CDB storage
// First check if we have any CDB storage set, because it is used
// to retrieve the calibration and alignment constants
+ if (fInitCDBCalled) return;
+ fInitCDBCalled = kTRUE;
+
AliCDBManager* man = AliCDBManager::Instance();
- if (!man->IsDefaultStorageSet())
+ if (man->IsDefaultStorageSet())
{
- AliWarningClass("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
- AliWarningClass("Default CDB storage not yet set");
- AliWarningClass(Form("Using default storage declared in AliSimulation: %s",uri));
- AliWarningClass("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
- SetDefaultStorage(uri);
- }
-
-}
-
-//_____________________________________________________________________________
-void AliSimulation::SetDefaultStorage(const char* uri) {
-// activate a default CDB storage
-
- AliCDBManager::Instance()->SetDefaultStorage(uri);
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ AliWarning("Default CDB storage has been already set !");
+ AliWarning(Form("Ignoring the default storage declared in AliSimulation: %s",fCDBUri.Data()));
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ fCDBUri = man->GetDefaultStorage()->GetURI();
+ }
+ else {
+ if (fCDBUri.Length() > 0)
+ {
+ AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ AliDebug(2, Form("Default CDB storage is set to: %s", fCDBUri.Data()));
+ AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ } else {
+ fCDBUri="local://$ALICE_ROOT";
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ AliWarning("Default CDB storage not yet set !!!!");
+ AliWarning(Form("Setting it now to: %s", fCDBUri.Data()));
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+
+ }
+ man->SetDefaultStorage(fCDBUri);
+ }
+ // Now activate the detector specific CDB storage locations
+ for (Int_t i = 0; i < fSpecCDBUri.GetEntriesFast(); i++) {
+ TObject* obj = fSpecCDBUri[i];
+ if (!obj) continue;
+ AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ AliDebug(2, Form("Specific CDB storage for %s is set to: %s",obj->GetName(),obj->GetTitle()));
+ AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ man->SetSpecificStorage(obj->GetName(), obj->GetTitle());
+ }
+
}
//_____________________________________________________________________________
-void AliSimulation::SetSpecificStorage(const char* detName, const char* uri) {
-// activate a detector-specific CDB storage
+void AliSimulation::InitRunNumber(){
+// check run number. If not set, set it to 0 !!!!
+
+ if (fInitRunNumberCalled) return;
+ fInitRunNumberCalled = kTRUE;
+
+ AliCDBManager* man = AliCDBManager::Instance();
+ if (man->GetRun() >= 0)
+ {
+ AliFatal(Form("Run number cannot be set in AliCDBManager before start of simulation: "
+ "Use external variable DC_RUN or AliSimulation::SetRun()!"));
+ }
+
+ if(fRun >= 0) {
+ AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ AliDebug(2, Form("Setting CDB run number to: %d",fRun));
+ AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ } else {
+ fRun=0;
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ AliWarning("Run number not yet set !!!!");
+ AliWarning(Form("Setting it now to: %d", fRun));
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+
+ }
+ man->SetRun(fRun);
- AliCDBManager::Instance()->SetSpecificStorage(detName, uri);
+ man->Print();
}
//_____________________________________________________________________________
-void AliSimulation::SetConfigFile(const char* fileName)
-{
-// set the name of the config file
-
- fConfigFileName = fileName;
+void AliSimulation::SetCDBLock() {
+ // Set CDB lock: from now on it is forbidden to reset the run number
+ // or the default storage or to activate any further storage!
+
+ AliCDBManager::Instance()->SetLock(1);
}
//_____________________________________________________________________________
-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
+void AliSimulation::SetDefaultStorage(const char* uri) {
+// Store the desired default CDB storage location
+// Activate it later within the Run() method
- fGAliceFileName = fileName;
- if (!gSystem->IsAbsoluteFileName(fGAliceFileName)) {
- char* absFileName = gSystem->ConcatFileName(gSystem->WorkingDirectory(),
- fGAliceFileName);
- fGAliceFileName = absFileName;
- delete[] absFileName;
- }
+ fCDBUri = uri;
- 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")
+void AliSimulation::SetSpecificStorage(const char* calibType, const char* uri) {
+// Store a detector-specific CDB storage location
+// Activate it later within the Run() method
+
+ AliCDBPath aPath(calibType);
+ if(!aPath.IsValid()){
+ AliError(Form("Not a valid path: %s", calibType));
+ return;
+ }
+
+ TObject* obj = fSpecCDBUri.FindObject(calibType);
+ if (obj) fSpecCDBUri.Remove(obj);
+ fSpecCDBUri.Add(new TNamed(calibType, uri));
- TNamed* obj = new TNamed(detector, type);
- obj->SetUniqueID(nEvents);
- fEventsPerFile.Add(obj);
}
//_____________________________________________________________________________
-Bool_t AliSimulation::ApplyAlignObjsToGeom(TObjArray* alObjArray)
+void AliSimulation::SetRunNumber(Int_t run)
{
- // Read collection of alignment objects (AliAlignObj derived) saved
- // in the TClonesArray ClArrayName and apply them to the geometry
- // manager singleton.
- //
- alObjArray->Sort();
- Int_t nvols = alObjArray->GetEntriesFast();
-
- Bool_t flag = kTRUE;
-
- for(Int_t j=0; j<nvols; j++)
- {
- AliAlignObj* alobj = (AliAlignObj*) alObjArray->UncheckedAt(j);
- if (alobj->ApplyToGeometry() == kFALSE) flag = kFALSE;
- }
-
- if (AliDebugLevelClass() >= 1) {
- gGeoManager->CheckOverlaps(20);
- TObjArray* ovexlist = gGeoManager->GetListOfOverlaps();
- if(ovexlist->GetEntriesFast()){
- AliError("The application of alignment objects to the geometry caused huge overlaps/extrusions!");
- }
- }
-
- return flag;
+// sets run number
+// Activate it later within the Run() method
+ fRun = run;
}
//_____________________________________________________________________________
-Bool_t AliSimulation::ApplyAlignObjsToGeom(const char* fileName, const char* clArrayName)
+void AliSimulation::SetSeed(Int_t seed)
{
- // 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 ApplyAlignObjsToGeom(alObjArray);
+// sets seed number
+// Activate it later within the Run() method
+ fSeed = seed;
}
//_____________________________________________________________________________
-Bool_t AliSimulation::ApplyAlignObjsToGeom(AliCDBParam* param, AliCDBId& Id)
+Bool_t AliSimulation::SetRunNumberFromData()
{
- // 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 ApplyAlignObjsToGeom(AlObjArray);
+ // Set the CDB manager run number
+ // The run number is retrieved from gAlice
+ if (fSetRunNumberFromDataCalled) return kTRUE;
+ fSetRunNumberFromDataCalled = kTRUE;
+
+ AliCDBManager* man = AliCDBManager::Instance();
+ Int_t runData = -1, runCDB = -1;
+
+ AliRunLoader* runLoader = LoadRun("READ");
+ if (!runLoader) return kFALSE;
+ else {
+ runData = runLoader->GetAliRun()->GetHeader()->GetRun();
+ delete runLoader;
+ }
+
+ runCDB = man->GetRun();
+ if(runCDB >= 0) {
+ if (runCDB != runData) {
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ AliWarning(Form("A run number was previously set in AliCDBManager: %d !", runCDB));
+ AliWarning(Form("It will be replaced with the run number got from run header: %d !", runData));
+ AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ }
+
+ }
+
+ man->SetRun(runData);
+ fRun = runData;
+
+ if(man->GetRun() < 0) {
+ AliError("Run number not properly initalized!");
+ return kFALSE;
+ }
+
+ man->Print();
+
+ return kTRUE;
}
//_____________________________________________________________________________
-Bool_t AliSimulation::ApplyAlignObjsToGeom(const char* uri, const char* path, Int_t runnum, Int_t version, Int_t sversion)
+void AliSimulation::SetConfigFile(const char* fileName)
{
- // 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);
+// set the name of the config file
+ fConfigFileName = fileName;
}
//_____________________________________________________________________________
-Bool_t AliSimulation::ApplyAlignObjsToGeom(const char* detName, Int_t runnum, Int_t version, Int_t sversion)
+void AliSimulation::SetGAliceFile(const char* fileName)
{
- // 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
- //
-
- InitCDBStorage("local://$ALICE_ROOT");
- AliCDBPath path(detName,"Align","Data");
- AliCDBEntry* entry = AliCDBManager::Instance()->Get(path.GetPath(),runnum,version,sversion);
+// set the name of the galice file
+// the path is converted to an absolute one if it is relative
- if(!entry) return kFALSE;
- TClonesArray* AlObjArray = ((TClonesArray*) entry->GetObject());
+ fGAliceFileName = fileName;
+ if (!gSystem->IsAbsoluteFileName(fGAliceFileName)) {
+ char* absFileName = gSystem->ConcatFileName(gSystem->WorkingDirectory(),
+ fGAliceFileName);
+ fGAliceFileName = absFileName;
+ delete[] absFileName;
+ }
- return ApplyAlignObjsToGeom(AlObjArray);
+ AliDebug(2, Form("galice file name set to %s", fileName));
}
//_____________________________________________________________________________
-Bool_t AliSimulation::SetAlignObjArraySingleDet(const char* detName)
+void AliSimulation::SetEventsPerFile(const char* detector, const char* type,
+ Int_t nEvents)
{
- // Fills array of single detector's alignable objects from CDB
-
- AliDebug(2, Form("Loading alignment data for detector: %s",detName));
-
- AliCDBEntry *entry;
-
- AliCDBPath path(detName,"Align","Data");
-
- entry=AliCDBManager::Instance()->Get(path.GetPath());
- if(!entry){
- AliDebug(2,Form("Couldn't load alignment data for detector %s",detName));
- return kFALSE;
- }
- entry->SetOwner(1);
- TClonesArray *alignArray = (TClonesArray*) entry->GetObject();
- alignArray->SetOwner(0);
- AliDebug(2,Form("Found %d alignment objects for %s",
- alignArray->GetEntries(),detName));
-
- AliAlignObj *alignObj=0;
- TIter iter(alignArray);
-
- // loop over align objects in detector
- while( ( alignObj=(AliAlignObj *) iter.Next() ) ){
- fAlignObjArray->Add(alignObj);
- }
- // delete entry --- Don't delete, it is cached!
-
- AliDebug(2, Form("fAlignObjArray entries: %d",fAlignObjArray->GetEntries() ));
- return kTRUE;
+// 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);
}
//_____________________________________________________________________________
// then applied to the TGeo geometry.
// Finally an overlaps check is performed.
+ if (!AliGeomManager::GetGeometry() || !AliGeomManager::GetGeometry()->IsClosed()) {
+ AliError("Can't apply the misalignment! Geometry is not loaded or it is still opened!");
+ return kFALSE;
+ }
+
+ // initialize CDB storage, run number, set CDB lock
+ InitCDB();
+// if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE;
+ SetCDBLock();
+
Bool_t delRunLoader = kFALSE;
if (!runLoader) {
runLoader = LoadRun("READ");
if (!runLoader) return kFALSE;
delRunLoader = kTRUE;
}
+
+ // Export ideal geometry
+ if(!gAlice->IsRootGeometry()) AliGeomManager::GetGeometry()->Export("geometry.root");
- // Load alignment data from CDB and fill fAlignObjArray
+ // Load alignment data from CDB and apply to geometry through AliGeomManager
if(fLoadAlignFromCDB){
- if(!fAlignObjArray) fAlignObjArray = new TObjArray();
-
- //fAlignObjArray->RemoveAll();
- fAlignObjArray->Clear();
- fAlignObjArray->SetOwner(0);
-
- TString detStr = fLoadAlignData;
- TString dataNotLoaded="";
- TString dataLoaded="";
-
- 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)) {
- if(!SetAlignObjArraySingleDet(det->GetName())){
- dataNotLoaded += det->GetName();
- dataNotLoaded += " ";
- } else {
- dataLoaded += det->GetName();
- dataLoaded += " ";
- }
- }
- } // end loop over detectors
-
- if ((detStr.CompareTo("ALL") == 0)) detStr = "";
- dataNotLoaded += detStr;
- AliInfo(Form("Alignment data loaded for: %s",
- dataLoaded.Data()));
- AliInfo(Form("Didn't/couldn't load alignment data for: %s",
- dataNotLoaded.Data()));
- } // fLoadAlignFromCDB flag
-
- // 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 (ApplyAlignObjsToGeom(fAlignObjArray) == kFALSE) {
- AliError("The misalignment of one or more volumes failed!"
- "Compare the list of simulated detectors and the list of detector alignment data!");
- if (delRunLoader) delete runLoader;
- return kFALSE;
+
+ TString detStr = fLoadAlObjsListOfDets;
+ TString loadAlObjsListOfDets = "";
+
+ 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)) {
+ //add det to list of dets to be aligned from CDB
+ loadAlObjsListOfDets += det->GetName();
+ loadAlObjsListOfDets += " ";
+ }
+ } // end loop over detectors
+ loadAlObjsListOfDets.Prepend("GRP "); //add alignment objects for non-sensitive modules
+ AliGeomManager::ApplyAlignObjsFromCDB(loadAlObjsListOfDets.Data());
+ }else{
+ // 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 (AliGeomManager::ApplyAlignObjsToGeom(*fAlignObjArray) == kFALSE) {
+ AliError("The misalignment of one or more volumes failed!"
+ "Compare the list of simulated detectors and the list of detector alignment data!");
+ if (delRunLoader) delete runLoader;
+ return kFALSE;
}
- }
- else {
- AliError("Can't apply the misalignment! gGeoManager doesn't exist or it is still opened!");
- if (delRunLoader) delete runLoader;
- return kFALSE;
}
}
- if (delRunLoader) delete runLoader;
+ // Update the internal geometry of modules (ITS needs it)
+ TString detStr = fLoadAlObjsListOfDets;
+ TObjArray* detArray = runLoader->GetAliRun()->Detectors();
+ for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
- return kTRUE;
-}
+ AliModule* det = (AliModule*) detArray->At(iDet);
+ if (!det || !det->IsActive()) continue;
+ if (IsSelected(det->GetName(), detStr)) {
+ det->UpdateInternalGeometry();
+ }
+ } // end loop over detectors
-//_____________________________________________________________________________
-Bool_t AliSimulation::SetRunNumber()
-{
- // Set the CDB manager run number
- // The run number is retrieved from gAlice
+ if (delRunLoader) delete runLoader;
- if(AliCDBManager::Instance()->GetRun() < 0) {
- AliRunLoader* runLoader = LoadRun("READ");
- if (!runLoader) return kFALSE;
- else {
- AliCDBManager::Instance()->SetRun(runLoader->GetAliRun()->GetRunNumber());
- AliInfo(Form("Run number: %d",AliCDBManager::Instance()->GetRun()));
- delete runLoader;
- }
- }
return kTRUE;
}
fBkgrdFileNames->Add(fileNameStr);
}
+void AliSimulation::EmbedInto(const char* fileName, Int_t nSignalPerBkgrd)
+{
+// add a file with background events for embeddin
+ MergeWith(fileName, nSignalPerBkgrd);
+ fEmbeddingFlag = kTRUE;
+}
//_____________________________________________________________________________
Bool_t AliSimulation::Run(Int_t nEvents)
{
// run the generation, simulation and digitization
- InitCDBStorage(fCDBUri);
+
+ AliCodeTimerAuto("")
+
+ // Load run number and seed from environmental vars
+ ProcessEnvironmentVars();
+ gRandom->SetSeed(fSeed);
+
if (nEvents > 0) fNEvents = nEvents;
// generation and simulation -> hits
if (fRunGeneration) {
if (!RunSimulation()) if (fStopOnError) return kFALSE;
}
-
- // Set run number in CDBManager (if it is not already set in RunSimulation)
- if (!SetRunNumber()) if (fStopOnError) return kFALSE;
-
- // Load and misalign the geometry
- if (!gGeoManager) {
- TGeoManager::Import("geometry.root");
- if (!gGeoManager) if (fStopOnError) return kFALSE;
- if (!MisalignGeometry()) if (fStopOnError) return kFALSE;
+
+ // initialize CDB storage from external environment
+ // (either CDB manager or AliSimulation setters),
+ // if not already done in RunSimulation()
+ InitCDB();
+
+ // Set run number in CDBManager from data
+ // From this point on the run number must be always loaded from data!
+ if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE;
+
+ // Set CDB lock: from now on it is forbidden to reset the run number
+ // or the default storage or to activate any further storage!
+ SetCDBLock();
+
+ // If RunSimulation was not called, load the geometry and misalign it
+ if (!AliGeomManager::GetGeometry()) {
+ // Initialize the geometry manager
+ AliGeomManager::LoadGeometry("geometry.root");
+
+// // Check that the consistency of symbolic names for the activated subdetectors
+// // in the geometry loaded by AliGeomManager
+// AliRunLoader* runLoader = LoadRun("READ");
+// if (!runLoader) return kFALSE;
+//
+// TString detsToBeChecked = "";
+// 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;
+// detsToBeChecked += det->GetName();
+// detsToBeChecked += " ";
+// } // end loop over detectors
+// if(!AliGeomManager::CheckSymNamesLUT(detsToBeChecked.Data()))
+ if(!AliGeomManager::CheckSymNamesLUT("ALL"))
+ AliFatalClass("Current loaded geometry differs in the definition of symbolic names!");
+
+ if (!AliGeomManager::GetGeometry()) if (fStopOnError) return kFALSE;
+ // Misalign geometry
+ if(!MisalignGeometry()) if (fStopOnError) return kFALSE;
}
+
// hits -> summable digits
+ AliSysInfo::AddStamp("Start_sdigitization");
if (!fMakeSDigits.IsNull()) {
if (!RunSDigitization(fMakeSDigits)) if (fStopOnError) return kFALSE;
+
}
-
- // summable digits -> digits
+ AliSysInfo::AddStamp("Stop_sdigitization");
+
+ AliSysInfo::AddStamp("Start_digitization");
+ // summable digits -> digits
if (!fMakeDigits.IsNull()) {
if (!RunDigitization(fMakeDigits, fMakeDigitsFromHits)) {
if (fStopOnError) return kFALSE;
}
- }
+ }
+ AliSysInfo::AddStamp("Stop_digitization");
+
+
// hits -> digits
if (!fMakeDigitsFromHits.IsNull()) {
if (fBkgrdFileNames && (fBkgrdFileNames->GetEntriesFast() > 0)) {
}
}
+
+
// digits -> trigger
- if (!RunTrigger(fMakeTrigger)) {
+ if (!RunTrigger(fMakeTrigger,fMakeDigits)) {
if (fStopOnError) return kFALSE;
}
+
+
// digits -> raw data
if (!fWriteRawData.IsNull()) {
if (!WriteRawData(fWriteRawData, fRawDataFileName,
- fDeleteIntermediateFiles)) {
+ fDeleteIntermediateFiles,fWriteSelRawData)) {
if (fStopOnError) return kFALSE;
}
}
+
+
+ // run HLT simulation
+ if (!fRunHLT.IsNull()) {
+ if (!RunHLT()) {
+ if (fStopOnError) return kFALSE;
+ }
+ }
+
+ //QA
+ if (fRunQA) {
+ Bool_t rv = RunQA() ;
+ if (!rv)
+ if (fStopOnError)
+ return kFALSE ;
+ }
+
+ // Cleanup of CDB manager: cache and active storages!
+ AliCDBManager::Instance()->ClearCache();
+
return kTRUE;
}
//_____________________________________________________________________________
-Bool_t AliSimulation::RunTrigger(const char* descriptors)
+Bool_t AliSimulation::RunTrigger(const char* config, const char* detectors)
{
// run the trigger
- TStopwatch stopwatch;
- stopwatch.Start();
+ AliCodeTimerAuto("")
+ // initialize CDB storage from external environment
+ // (either CDB manager or AliSimulation setters),
+ // if not already done in RunSimulation()
+ InitCDB();
+
+ // Set run number in CDBManager from data
+ // From this point on the run number must be always loaded from data!
+ if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE;
+
+ // Set CDB lock: from now on it is forbidden to reset the run number
+ // or the default storage or to activate any further storage!
+ SetCDBLock();
+
AliRunLoader* runLoader = LoadRun("READ");
if (!runLoader) return kFALSE;
- TString des = descriptors;
+ TString trconfiguration = config;
- if (des.IsNull()) {
+ if (trconfiguration.IsNull()) {
if (gAlice->GetTriggerDescriptor() != "") {
- des = gAlice->GetTriggerDescriptor();
- }
- else {
- AliWarning("No trigger descriptor is specified. Skipping the trigger simulation...");
- return kTRUE;
+ trconfiguration = gAlice->GetTriggerDescriptor();
}
+ else
+ AliWarning("No trigger descriptor is specified. Loading the one that is in the CDB.");
}
- runLoader->MakeTree( "CT" );
+ runLoader->MakeTree( "GG" );
AliCentralTrigger* aCTP = runLoader->GetTrigger();
- // Load Descriptors
- aCTP->LoadDescriptor( des );
+ // Load Configuration
+ if (!aCTP->LoadConfiguration( trconfiguration ))
+ return kFALSE;
- // digits -> trigger
- if( !aCTP->RunTrigger( runLoader ) ) {
+ // digits -> trigger
+ if( !aCTP->RunTrigger( runLoader , detectors ) ) {
if (fStopOnError) {
- // delete aCTP;
- return kFALSE;
+ // delete aCTP;
+ return kFALSE;
}
}
- AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
- stopwatch.RealTime(),stopwatch.CpuTime()));
-
delete runLoader;
return kTRUE;
{
// run the generation and simulation
- TStopwatch stopwatch;
- stopwatch.Start();
+ AliCodeTimerAuto("")
+ // initialize CDB storage and run number from external environment
+ // (either CDB manager or AliSimulation setters)
+ InitCDB();
+ InitRunNumber();
+ SetCDBLock();
+
if (!gAlice) {
AliError("no gAlice object. Restart aliroot and try again.");
return kFALSE;
StdoutToAliInfo(StderrToAliError(
gAlice->Init(fConfigFileName.Data());
););
-
+
// Get the trigger descriptor string
// Either from AliSimulation or from
// gAlice
gAlice->SetTriggerDescriptor(fMakeTrigger.Data());
// Set run number in CDBManager
- AliCDBManager::Instance()->SetRun(gAlice->GetRunNumber());
AliInfo(Form("Run number: %d",AliCDBManager::Instance()->GetRun()));
AliRunLoader* runLoader = gAlice->GetRunLoader();
return kFALSE;
}
SetGAliceFile(runLoader->GetFileName());
-
- // Export ideal geometry
- if (gGeoManager) gGeoManager->Export("geometry.root");
-
+
// Misalign geometry
-// if (!MisalignGeometry(runLoader)) {
-// delete runLoader;
-// return kFALSE;
-// }
+#if ROOT_VERSION_CODE < 331527
+ AliGeomManager::SetGeometry(gGeoManager);
+
+ // Check that the consistency of symbolic names for the activated subdetectors
+ // in the geometry loaded by AliGeomManager
+ TString detsToBeChecked = "";
+ 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;
+ detsToBeChecked += det->GetName();
+ detsToBeChecked += " ";
+ } // end loop over detectors
+ if(!AliGeomManager::CheckSymNamesLUT(detsToBeChecked.Data()))
+ AliFatalClass("Current loaded geometry differs in the definition of symbolic names!");
MisalignGeometry(runLoader);
-
- // Temporary fix by A.Gheata
- // Could be removed with the next Root version (>5.11)
- if (gGeoManager) {
- TIter next(gGeoManager->GetListOfVolumes());
- TGeoVolume *vol;
- while ((vol = (TGeoVolume *)next())) {
- if (vol->GetVoxels()) {
- if (vol->GetVoxels()->NeedRebuild()) {
- vol->GetVoxels()->Voxelize();
- vol->FindOverlaps();
- }
- }
- }
- }
-
- // Export (mis)aligned geometry
- if (gGeoManager) gGeoManager->Export("misaligned_geometry.root");
+#endif
// AliRunLoader* runLoader = gAlice->GetRunLoader();
// if (!runLoader) {
}
AliInfo("running gAlice");
+ AliSysInfo::AddStamp("Start_simulation");
StdoutToAliInfo(StderrToAliError(
gAlice->Run(nEvents);
););
-
+ AliSysInfo::AddStamp("Stop_simulation");
delete runLoader;
- AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
- stopwatch.RealTime(),stopwatch.CpuTime()));
-
return kTRUE;
}
Bool_t AliSimulation::RunSDigitization(const char* detectors)
{
// run the digitization and produce summable digits
+ static Int_t eventNr=0;
+ AliCodeTimerAuto("")
- TStopwatch stopwatch;
- stopwatch.Start();
-
+ // initialize CDB storage, run number, set CDB lock
+ InitCDB();
+ if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE;
+ SetCDBLock();
+
AliRunLoader* runLoader = LoadRun();
if (!runLoader) return kFALSE;
if (!det || !det->IsActive()) continue;
if (IsSelected(det->GetName(), detStr)) {
AliInfo(Form("creating summable digits for %s", det->GetName()));
- TStopwatch stopwatchDet;
- stopwatchDet.Start();
+ AliCodeTimerAuto(Form("creating summable digits for %s", det->GetName()));
det->Hits2SDigits();
- AliInfo(Form("Execution time for %s: R:%.2fs C:%.2fs",
- det->GetName(),stopwatchDet.RealTime(),stopwatchDet.CpuTime()));
+ AliSysInfo::AddStamp(Form("Digit_%s_%d",det->GetName(),eventNr), 0,1, eventNr);
}
}
detStr.Data()));
if (fStopOnError) return kFALSE;
}
-
+ eventNr++;
delete runLoader;
- AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
- stopwatch.RealTime(),stopwatch.CpuTime()));
-
return kTRUE;
}
{
// run the digitization and produce digits from sdigits
- TStopwatch stopwatch;
- stopwatch.Start();
+ AliCodeTimerAuto("")
+ // initialize CDB storage, run number, set CDB lock
+ InitCDB();
+ if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE;
+ SetCDBLock();
+
while (AliRunLoader::GetRunLoader()) delete AliRunLoader::GetRunLoader();
if (gAlice) delete gAlice;
gAlice = NULL;
if (fBkgrdFileNames) nStreams = fBkgrdFileNames->GetEntriesFast() + 1;
Int_t signalPerBkgrd = GetNSignalPerBkgrd();
AliRunDigitizer* manager = new AliRunDigitizer(nStreams, signalPerBkgrd);
+ // manager->SetEmbeddingFlag(fEmbeddingFlag);
manager->SetInputStream(0, fGAliceFileName.Data());
for (Int_t iStream = 1; iStream < nStreams; iStream++) {
const char* fileName = ((TObjString*)
if (IsSelected(det->GetName(), detStr) &&
!IsSelected(det->GetName(), detExcl)) {
AliDigitizer* digitizer = det->CreateDigitizer(manager);
+
if (!digitizer) {
AliError(Form("no digitizer for %s", det->GetName()));
if (fStopOnError) return kFALSE;
delete manager;
- AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
- stopwatch.RealTime(),stopwatch.CpuTime()));
-
return kTRUE;
}
{
// run the digitization and produce digits from hits
- TStopwatch stopwatch;
- stopwatch.Start();
+ AliCodeTimerAuto("")
+ // initialize CDB storage, run number, set CDB lock
+ InitCDB();
+ if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE;
+ SetCDBLock();
+
AliRunLoader* runLoader = LoadRun("READ");
if (!runLoader) return kFALSE;
//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
- 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)
+ Bool_t deleteIntermediateFiles,
+ Bool_t selrawdata)
{
// convert the digits to raw data
// First DDL raw data files for the given detectors are created.
// 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.
+// 'selrawdata' flag can be used to enable writing of detectors raw data
+// accoring to the trigger cluster.
- TStopwatch stopwatch;
- stopwatch.Start();
+ AliCodeTimerAuto("")
if (!WriteRawFiles(detectors)) {
if (fStopOnError) return kFALSE;
if (!dateFileName.IsNull()) {
Bool_t rootOutput = dateFileName.EndsWith(".root");
if (rootOutput) dateFileName += ".date";
- if (!ConvertRawFilesToDate(dateFileName)) {
+ TString selDateFileName;
+ if (selrawdata) {
+ selDateFileName = "selected.";
+ selDateFileName+= dateFileName;
+ }
+ if (!ConvertRawFilesToDate(dateFileName,selDateFileName)) {
if (fStopOnError) return kFALSE;
}
if (deleteIntermediateFiles) {
if (deleteIntermediateFiles) {
gSystem->Unlink(dateFileName);
}
+ if (selrawdata) {
+ TString selFileName = "selected.";
+ selFileName += fileName;
+ if (!ConvertDateToRoot(selDateFileName, selFileName)) {
+ if (fStopOnError) return kFALSE;
+ }
+ if (deleteIntermediateFiles) {
+ gSystem->Unlink(selDateFileName);
+ }
+ }
}
}
- AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
- stopwatch.RealTime(),stopwatch.CpuTime()));
-
return kTRUE;
}
{
// convert the digits to raw data DDL files
+ AliCodeTimerAuto("")
+
AliRunLoader* runLoader = LoadRun("READ");
if (!runLoader) return kFALSE;
}
delete runLoader;
+
return kTRUE;
}
//_____________________________________________________________________________
-Bool_t AliSimulation::ConvertRawFilesToDate(const char* dateFileName)
+Bool_t AliSimulation::ConvertRawFilesToDate(const char* dateFileName,
+ const char* selDateFileName)
{
// convert raw data DDL files to a DATE file with the program "dateStream"
+// The second argument is not empty when the user decides to write
+// the detectors raw data according to the trigger cluster.
+ AliCodeTimerAuto("")
+
char* path = gSystem->Which(gSystem->Getenv("PATH"), "dateStream");
if (!path) {
AliError("the program dateStream was not found");
if (!runLoader) return kFALSE;
AliInfo(Form("converting raw data DDL files to DATE file %s", dateFileName));
+ Bool_t selrawdata = kFALSE;
+ if (strcmp(selDateFileName,"") != 0) selrawdata = kTRUE;
+
char command[256];
- sprintf(command, "dateStream -D -o %s -# %d -C",
- dateFileName, runLoader->GetNumberOfEvents());
+ // Note the option -s. It is used in order to avoid
+ // the generation of SOR/EOR events.
+ sprintf(command, "dateStream -c -s -D -o %s -# %d -C -run %d",
+ dateFileName, runLoader->GetNumberOfEvents(),runLoader->GetHeader()->GetRun());
FILE* pipe = gSystem->OpenPipe(command, "w");
+ Int_t selEvents = 0;
for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
fprintf(pipe, "GDC\n");
Float_t ldc = 0;
Int_t prevLDC = -1;
+ if (selrawdata) {
+ // Check if the event was triggered by CTP
+ runLoader->GetEvent(iEvent);
+ if (!runLoader->LoadTrigger()) {
+ AliCentralTrigger *aCTP = runLoader->GetTrigger();
+ if (aCTP->GetClassMask()) selEvents++;
+ }
+ else {
+ AliWarning("No trigger can be loaded! Writing of selected raw data is abandoned !");
+ selrawdata = kFALSE;
+ }
+ }
+
// loop over detectors and DDLs
- for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
- for (Int_t iDDL = 0; iDDL < kDetectorDDLs[iDet]; iDDL++) {
+ for (Int_t iDet = 0; iDet < AliDAQ::kNDetectors; iDet++) {
+ for (Int_t iDDL = 0; iDDL < AliDAQ::NumberOfDdls(iDet); iDDL++) {
- Int_t ddlID = 0x100*iDet + iDDL;
+ Int_t ddlID = AliDAQ::DdlID(iDet,iDDL);
Int_t ldcID = Int_t(ldc + 0.0001);
- ldc += kDetectorLDCs[iDet] / kDetectorDDLs[iDet];
+ ldc += AliDAQ::NumberOfLdcs(iDet) / AliDAQ::NumberOfDdls(iDet);
char rawFileName[256];
- sprintf(rawFileName, "raw%d/%s_%d.ddl",
- iEvent, kDetectors[iDet], ddlID);
+ sprintf(rawFileName, "raw%d/%s",
+ iEvent, AliDAQ::DdlFileName(iDet,iDDL));
// check existence and size of raw data file
FILE* file = fopen(rawFileName, "rb");
Int_t result = gSystem->ClosePipe(pipe);
+ if (!(selrawdata && selEvents > 0)) {
+ delete runLoader;
+ return (result == 0);
+ }
+
+ AliInfo(Form("converting selected by trigger cluster raw data DDL files to DATE file %s", selDateFileName));
+
+ sprintf(command, "dateStream -c -s -D -o %s -# %d -C -run %d",
+ selDateFileName,selEvents,runLoader->GetHeader()->GetRun());
+ FILE* pipe2 = gSystem->OpenPipe(command, "w");
+
+ for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
+
+ // Get the trigger decision and cluster
+ TString detClust;
+ runLoader->GetEvent(iEvent);
+ if (!runLoader->LoadTrigger()) {
+ AliCentralTrigger *aCTP = runLoader->GetTrigger();
+ if (aCTP->GetClassMask() == 0) continue;
+ detClust = aCTP->GetTriggeredDetectors();
+ AliInfo(Form("List of detectors to be read out: %s",detClust.Data()));
+ }
+
+ fprintf(pipe2, "GDC\n");
+ Float_t ldc = 0;
+ Int_t prevLDC = -1;
+
+ // loop over detectors and DDLs
+ for (Int_t iDet = 0; iDet < AliDAQ::kNDetectors; iDet++) {
+ // Write only raw data from detectors that
+ // are contained in the trigger cluster(s)
+ if (!IsSelected(AliDAQ::DetectorName(iDet),detClust)) continue;
+
+ for (Int_t iDDL = 0; iDDL < AliDAQ::NumberOfDdls(iDet); iDDL++) {
+
+ Int_t ddlID = AliDAQ::DdlID(iDet,iDDL);
+ Int_t ldcID = Int_t(ldc + 0.0001);
+ ldc += AliDAQ::NumberOfLdcs(iDet) / AliDAQ::NumberOfDdls(iDet);
+
+ char rawFileName[256];
+ sprintf(rawFileName, "raw%d/%s",
+ iEvent, AliDAQ::DdlFileName(iDet,iDDL));
+
+ // 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(pipe2, " LDC Id %d\n", ldcID);
+ prevLDC = ldcID;
+ }
+ fprintf(pipe2, " Equipment Id %d Payload %s\n", ddlID, rawFileName);
+ }
+ }
+ }
+
+ Int_t result2 = gSystem->ClosePipe(pipe2);
+
delete runLoader;
- return (result == 0);
+ return ((result == 0) && (result2 == 0));
}
//_____________________________________________________________________________
// convert a DATE file to a root file with the program "alimdc"
// ALIMDC setup
- const Int_t kDBSize = 1000000000;
+ const Int_t kDBSize = 2000000000;
const Int_t kTagDBSize = 1000000000;
const Bool_t kFilter = kFALSE;
const Int_t kCompression = 1;
const char* rawDBFS[2] = { "/tmp/mdc1", "/tmp/mdc2" };
const char* tagDBFS = "/tmp/mdc1/tags";
- const char* runDBFS = "/tmp/mdc1/meta";
// User defined file system locations
if (gSystem->Getenv("ALIMDC_RAWDB1"))
rawDBFS[1] = gSystem->Getenv("ALIMDC_RAWDB2");
if (gSystem->Getenv("ALIMDC_TAGDB"))
tagDBFS = gSystem->Getenv("ALIMDC_TAGDB");
- if (gSystem->Getenv("ALIMDC_RUNDB"))
- runDBFS = gSystem->Getenv("ALIMDC_RUNDB");
gSystem->Exec(Form("rm -rf %s",rawDBFS[0]));
gSystem->Exec(Form("rm -rf %s",rawDBFS[1]));
gSystem->Exec(Form("rm -rf %s",tagDBFS));
- gSystem->Exec(Form("rm -rf %s",runDBFS));
gSystem->Exec(Form("mkdir %s",rawDBFS[0]));
gSystem->Exec(Form("mkdir %s",rawDBFS[1]));
gSystem->Exec(Form("mkdir %s",tagDBFS));
- gSystem->Exec(Form("mkdir %s",runDBFS));
Int_t result = gSystem->Exec(Form("alimdc %d %d %d %d %s",
kDBSize, kTagDBSize, kFilter, kCompression, dateFileName));
gSystem->Exec(Form("rm -rf %s",rawDBFS[0]));
gSystem->Exec(Form("rm -rf %s",rawDBFS[1]));
gSystem->Exec(Form("rm -rf %s",tagDBFS));
- gSystem->Exec(Form("rm -rf %s",runDBFS));
return (result == 0);
}
return NULL;
}
runLoader->LoadgAlice();
+ runLoader->LoadHeader();
gAlice = runLoader->GetAliRun();
if (!gAlice) {
AliError(Form("no gAlice object found in file %s",
// get the number of signal events
if (nEvents <= 0) {
AliRunLoader* runLoader =
- AliRunLoader::Open(fGAliceFileName.Data(), "SIGNAL");
+ AliRunLoader::Open(fGAliceFileName.Data(), "SIGNAL");
if (!runLoader) return 1;
+
nEvents = runLoader->GetNumberOfEvents();
delete runLoader;
}
// get the number of background events
const char* fileName = ((TObjString*)
(fBkgrdFileNames->At(iBkgrdFile)))->GetName();
- AliRunLoader* runLoader =
+ 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) {
return result;
}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* esdFileName)
+{
+//
+// Steering routine to convert raw data in directory rawDirectory/ to fake SDigits.
+// These can be used for embedding of MC tracks into RAW data using the standard
+// merging procedure.
+//
+// If an ESD file is given the reconstructed vertex is taken from it and stored in the event header.
+//
+ if (!gAlice) {
+ AliError("no gAlice object. Restart aliroot and try again.");
+ return kFALSE;
+ }
+ if (gAlice->Modules()->GetEntries() > 0) {
+ AliError("gAlice was already run. Restart aliroot and try again.");
+ return kFALSE;
+ }
+
+ AliInfo(Form("initializing gAlice with config file %s",fConfigFileName.Data()));
+ StdoutToAliInfo(StderrToAliError(gAlice->Init(fConfigFileName.Data());););
+//
+// Initialize CDB
+ InitCDB();
+ //AliCDBManager* man = AliCDBManager::Instance();
+ //man->SetRun(0); // Should this come from rawdata header ?
+
+ Int_t iDet;
+ //
+ // Get the runloader
+ AliRunLoader* runLoader = gAlice->GetRunLoader();
+ //
+ // Open esd file if available
+ TFile* esdFile = TFile::Open(esdFileName);
+ Bool_t esdOK = (esdFile != 0);
+ AliESD* esd = new AliESD;
+ TTree* treeESD = 0;
+ if (esdOK) {
+ treeESD = (TTree*) esdFile->Get("esdTree");
+ if (!treeESD) {
+ AliWarning("No ESD tree found");
+ esdOK = kFALSE;
+ } else {
+ treeESD->SetBranchAddress("ESD", &esd);
+ }
+ }
+ //
+ // Create the RawReader
+ TString fileName(rawDirectory);
+ AliRawReader* rawReader = 0x0;
+ if (fileName.EndsWith("/")) {
+ rawReader = new AliRawReaderFile(fileName);
+ } else if (fileName.EndsWith(".root")) {
+ rawReader = new AliRawReaderRoot(fileName);
+ } else if (!fileName.IsNull()) {
+ rawReader = new AliRawReaderDate(fileName);
+ rawReader->SelectEvents(7);
+ }
+// if (!fEquipIdMap.IsNull() && fRawReader)
+// fRawReader->LoadEquipmentIdsMap(fEquipIdMap);
+ //
+ // Get list of detectors
+ TObjArray* detArray = runLoader->GetAliRun()->Detectors();
+ //
+ // Get Header
+ AliHeader* header = runLoader->GetHeader();
+ //
+ TString detStr = fMakeSDigits;
+ // Event loop
+ Int_t nev = 0;
+ while(kTRUE) {
+ if (!(rawReader->NextEvent())) break;
+ //
+ // Detector loop
+ for (iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
+ AliModule* det = (AliModule*) detArray->At(iDet);
+ if (!det || !det->IsActive()) continue;
+ if (IsSelected(det->GetName(), detStr)) {
+ AliInfo(Form("Calling Raw2SDigits for %s\n", det->GetName()));
+ det->Raw2SDigits(rawReader);
+ rawReader->Reset();
+ }
+ } // detectors
+
+
+ //
+ // If ESD information available obtain reconstructed vertex and store in header.
+ if (esdOK) {
+ treeESD->GetEvent(nev);
+ const AliESDVertex* esdVertex = esd->GetPrimaryVertex();
+ Double_t position[3];
+ esdVertex->GetXYZ(position);
+ AliGenEventHeader* mcHeader = new AliGenEventHeader("ESD");
+ TArrayF mcV;
+ mcV.Set(3);
+ for (Int_t i = 0; i < 3; i++) mcV[i] = position[i];
+ mcHeader->SetPrimaryVertex(mcV);
+ header->Reset(0,nev);
+ header->SetGenEventHeader(mcHeader);
+ printf("***** Saved vertex %f %f %f \n", position[0], position[1], position[2]);
+ }
+ nev++;
+//
+// Finish the event
+ runLoader->TreeE()->Fill();
+ runLoader->SetNextEvent();
+ } // events
+
+ delete rawReader;
+//
+// Finish the run
+ runLoader->CdGAFile();
+ runLoader->WriteHeader("OVERWRITE");
+ runLoader->WriteRunLoader();
+
+ return kTRUE;
+}
+
+//_____________________________________________________________________________
+Int_t AliSimulation::GetDetIndex(const char* detector)
+{
+ // return the detector index corresponding to detector
+ Int_t index = -1 ;
+ for (index = 0; index < fgkNDetectors ; index++) {
+ if ( strcmp(detector, fgkDetectorName[index]) == 0 )
+ break ;
+ }
+ return index ;
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::RunHLT()
+{
+ // Run the HLT simulation
+ // HLT simulation is implemented in HLT/sim/AliHLTSimulation
+ // Disabled if fRunHLT is empty, default vaule is "default".
+ // AliSimulation::SetRunHLT can be used to set the options for HLT simulation
+ // The default simulation depends on the HLT component libraries and their
+ // corresponding agents which define components and chains to run. See
+ // http://web.ift.uib.no/~kjeks/doc/alice-hlt/
+ // http://web.ift.uib.no/~kjeks/doc/alice-hlt/classAliHLTModuleAgent.html
+ //
+ // The libraries to be loaded can be specified as an option.
+ // <pre>
+ // AliSimulation sim;
+ // sim.SetRunHLT("libAliHLTSample.so");
+ // </pre>
+ // will only load <tt>libAliHLTSample.so</tt>
+
+ // Other available options:
+ // \li loglevel=<i>level</i> <br>
+ // logging level for this processing
+ // \li alilog=off
+ // disable redirection of log messages to AliLog class
+ // \li config=<i>macro</i>
+ // configuration macro
+ // \li localrec=<i>configuration</i>
+ // comma separated list of configurations to be run during simulation
+
+ int iResult=0;
+ AliRunLoader* pRunLoader = LoadRun("READ");
+ if (!pRunLoader) return kFALSE;
+
+ // initialize CDB storage, run number, set CDB lock
+ InitCDB();
+ if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE;
+ SetCDBLock();
+
+ // load the library dynamically
+ gSystem->Load(ALIHLTSIMULATION_LIBRARY);
+
+ // check for the library version
+ AliHLTSimulationGetLibraryVersion_t fctVersion=(AliHLTSimulationGetLibraryVersion_t)(gSystem->DynFindSymbol(ALIHLTSIMULATION_LIBRARY, ALIHLTSIMULATION_GET_LIBRARY_VERSION));
+ if (!fctVersion) {
+ AliError(Form("can not load library %s", ALIHLTSIMULATION_LIBRARY));
+ return kFALSE;
+ }
+ if (fctVersion()!= ALIHLTSIMULATION_LIBRARY_VERSION) {
+ AliError(Form("%s version does not match: compiled for version %d, loaded %d", ALIHLTSIMULATION_LIBRARY, ALIHLTSIMULATION_LIBRARY_VERSION, fctVersion()));
+ return kFALSE;
+ }
+
+ // print compile info
+ typedef void (*CompileInfo)( char*& date, char*& time);
+ CompileInfo fctInfo=(CompileInfo)gSystem->DynFindSymbol(ALIHLTSIMULATION_LIBRARY, "CompileInfo");
+ if (fctInfo) {
+ char* date="";
+ char* time="";
+ (*fctInfo)(date, time);
+ if (!date) date="unknown";
+ if (!time) time="unknown";
+ AliInfo(Form("%s build on %s (%s)", ALIHLTSIMULATION_LIBRARY, date, time));
+ } else {
+ AliInfo(Form("no build info available for %s", ALIHLTSIMULATION_LIBRARY));
+ }
+
+ // create instance of the HLT simulation
+ AliHLTSimulationCreateInstance_t fctCreate=(AliHLTSimulationCreateInstance_t)(gSystem->DynFindSymbol(ALIHLTSIMULATION_LIBRARY, ALIHLTSIMULATION_CREATE_INSTANCE));
+ AliHLTSimulation* pHLT=NULL;
+ if (fctCreate==NULL || (pHLT=(fctCreate()))==NULL) {
+ AliError(Form("can not create instance of HLT simulation (creator %p)", fctCreate));
+ return kFALSE;
+ }
+
+ // init the HLT simulation
+ TString options;
+ if (fRunHLT.CompareTo("default")!=0) options=fRunHLT;
+ if (!IsSelected("HLT", fWriteRawData)) {
+ options+=" writerawfiles=";
+ } else {
+ options+=" writerawfiles=HLT";
+ }
+ AliHLTSimulationInit_t fctInit=(AliHLTSimulationInit_t)(gSystem->DynFindSymbol(ALIHLTSIMULATION_LIBRARY, ALIHLTSIMULATION_INIT));
+ if (fctInit==NULL || (iResult=(fctInit(pHLT, pRunLoader, options.Data())))<0) {
+ AliError(Form("can not init HLT simulation: error %d (init %p)", iResult, fctInit));
+ } else {
+ // run the HLT simulation
+ AliHLTSimulationRun_t fctRun=(AliHLTSimulationRun_t)(gSystem->DynFindSymbol(ALIHLTSIMULATION_LIBRARY, ALIHLTSIMULATION_RUN));
+ if (fctRun==NULL || (iResult=(fctRun(pHLT, pRunLoader)))<0) {
+ AliError(Form("can not run HLT simulation: error %d (run %p)", iResult, fctRun));
+ }
+ }
+
+ // delete the instance
+ AliHLTSimulationDeleteInstance_t fctDelete=(AliHLTSimulationDeleteInstance_t)(gSystem->DynFindSymbol(ALIHLTSIMULATION_LIBRARY, ALIHLTSIMULATION_DELETE_INSTANCE));
+ if (fctDelete==NULL || fctDelete(pHLT)<0) {
+ AliError(Form("can not delete instance of HLT simulation (creator %p)", fctDelete));
+ }
+ pHLT=NULL;
+
+ return iResult>=0?kTRUE:kFALSE;
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::RunQA()
+{
+ // run the QA on summable hits, digits or digits
+
+ AliQADataMakerSteer qas ;
+ qas.SetRunLoader(gAlice->GetRunLoader()) ;
+
+ TString detectorsw("") ;
+ Bool_t rv = kTRUE ;
+ if (fQATasks.Contains(Form("%d", AliQA::kHITS)))
+ detectorsw = qas.Run(fQADetectors.Data(), AliQA::kHITS) ;
+// qas.Reset() ;
+ if (fQATasks.Contains(Form("%d", AliQA::kSDIGITS)))
+ detectorsw += qas.Run(fQADetectors.Data(), AliQA::kSDIGITS) ;
+// qas.Reset() ;
+ if (fQATasks.Contains(Form("%d", AliQA::kDIGITS)))
+ detectorsw += qas.Run(fQADetectors.Data(), AliQA::kDIGITS) ;
+
+ if ( detectorsw.IsNull() )
+ rv = kFALSE ;
+ return rv ;
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::SetRunQA(TString detAndAction)
+{
+ // Allows to run QA for a selected set of detectors
+ // and a selected set of tasks among HITS, SDIGITS and DIGITS
+ // all selected detectors run the same selected tasks
+
+ if (!detAndAction.Contains(":")) {
+ AliError( Form("%s is a wrong syntax, use \"DetectorList:ActionList\" \n", detAndAction.Data()) ) ;
+ fRunQA = kFALSE ;
+ return kFALSE ;
+ }
+ Int_t colon = detAndAction.Index(":") ;
+ fQADetectors = detAndAction(0, colon) ;
+ if (fQADetectors.Contains("ALL") )
+ fQADetectors = Form("%s %s", fMakeDigits.Data(), fMakeDigitsFromHits.Data()) ;
+ fQATasks = detAndAction(colon+1, detAndAction.Sizeof() ) ;
+ if (fQATasks.Contains("ALL") ) {
+ fQATasks = Form("%d %d %d", AliQA::kHITS, AliQA::kSDIGITS, AliQA::kDIGITS) ;
+ } else {
+ fQATasks.ToUpper() ;
+ TString tempo("") ;
+ if ( fQATasks.Contains("HIT") )
+ tempo = Form("%d ", AliQA::kHITS) ;
+ if ( fQATasks.Contains("SDIGIT") )
+ tempo += Form("%d ", AliQA::kSDIGITS) ;
+ if ( fQATasks.Contains("DIGIT") )
+ tempo += Form("%d ", AliQA::kDIGITS) ;
+ fQATasks = tempo ;
+ if (fQATasks.IsNull()) {
+ AliInfo("No QA requested\n") ;
+ fRunQA = kFALSE ;
+ return kTRUE ;
+ }
+ }
+ TString tempo(fQATasks) ;
+ tempo.ReplaceAll(Form("%d", AliQA::kHITS), AliQA::GetTaskName(AliQA::kHITS)) ;
+ tempo.ReplaceAll(Form("%d", AliQA::kSDIGITS), AliQA::GetTaskName(AliQA::kSDIGITS)) ;
+ tempo.ReplaceAll(Form("%d", AliQA::kDIGITS), AliQA::GetTaskName(AliQA::kDIGITS)) ;
+ AliInfo( Form("QA will be done on \"%s\" for \"%s\"\n", fQADetectors.Data(), tempo.Data()) ) ;
+ fRunQA = kTRUE ;
+ return kTRUE;
+}
+
+//_____________________________________________________________________________
+void AliSimulation::ProcessEnvironmentVars()
+{
+// Extract run number and random generator seed from env variables
+
+ AliInfo("Processing environment variables");
+
+ // Random Number seed
+
+ // first check that seed is not already set
+ if (fSeed == 0) {
+ if (gSystem->Getenv("CONFIG_SEED")) {
+ fSeed = atoi(gSystem->Getenv("CONFIG_SEED"));
+ }
+ } else {
+ if (gSystem->Getenv("CONFIG_SEED")) {
+ AliInfo(Form("Seed for random number generation already set (%d)"
+ ": CONFIG_SEED variable ignored!", fSeed));
+ }
+ }
+
+ AliInfo(Form("Seed for random number generation = %d ", fSeed));
+
+ // Run Number
+
+ // first check that run number is not already set
+ if(fRun < 0) {
+ if (gSystem->Getenv("DC_RUN")) {
+ fRun = atoi(gSystem->Getenv("DC_RUN"));
+ }
+ } else {
+ if (gSystem->Getenv("DC_RUN")) {
+ AliInfo(Form("Run number already set (%d): DC_RUN variable ignored!", fRun));
+ }
+ }
+
+ AliInfo(Form("Run number = %d", fRun));
+}
+
+
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff