#include "AliGenEventHeader.h"
#include "AliMC.h"
#include "AliHLTSimulation.h"
-#include "AliQA.h"
#include "AliQADataMakerSteer.h"
+#include "AliSysInfo.h"
ClassImp(AliSimulation)
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),
fWriteRawData(""),
fRawDataFileName(""),
fDeleteIntermediateFiles(kFALSE),
+ fWriteSelRawData(kFALSE),
fStopOnError(kFALSE),
fNEvents(1),
fAlignObjArray(NULL),
fUseBkgrdVertex(kTRUE),
fRegionOfInterest(kFALSE),
- fCDBUri(cdbUri),
- fRemoteCDBUri(""),
+ fCDBUri(""),
fSpecCDBUri(),
+ fRun(-1),
+ fSeed(0),
+ fInitCDBCalled(kFALSE),
+ fInitRunNumberCalled(kFALSE),
+ fSetRunNumberFromDataCalled(kFALSE),
fEmbeddingFlag(kFALSE),
+ fRunQA(kTRUE),
fRunHLT("default")
{
// create simulation object with default parameters
fWriteRawData(sim.fWriteRawData),
fRawDataFileName(""),
fDeleteIntermediateFiles(kFALSE),
+ fWriteSelRawData(kFALSE),
fStopOnError(sim.fStopOnError),
fNEvents(sim.fNEvents),
fUseBkgrdVertex(sim.fUseBkgrdVertex),
fRegionOfInterest(sim.fRegionOfInterest),
fCDBUri(sim.fCDBUri),
- fRemoteCDBUri(sim.fRemoteCDBUri),
fSpecCDBUri(),
+ fRun(-1),
+ fSeed(0),
+ fInitCDBCalled(sim.fInitCDBCalled),
+ fInitRunNumberCalled(sim.fInitRunNumberCalled),
+ fSetRunNumberFromDataCalled(sim.fSetRunNumberFromDataCalled),
fEmbeddingFlag(sim.fEmbeddingFlag),
+ fRunQA(kTRUE),
fRunHLT(sim.fRunHLT)
{
// copy constructor
}
//_____________________________________________________________________________
-void AliSimulation::InitCDBStorage()
+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())
{
AliWarning("Default CDB storage has been already set !");
AliWarning(Form("Ignoring the default storage declared in AliSimulation: %s",fCDBUri.Data()));
AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
- fCDBUri = "";
+ fCDBUri = man->GetDefaultStorage()->GetURI();
}
else {
- AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
- AliDebug(2, Form("Default CDB storage is set to: %s",fCDBUri.Data()));
- AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
+ 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);
}
- // Remote storage (the Grid storage) is used if it is activated
- // and if the object is not found in the default storage
- // OBSOLETE: Removed
- // if (man->IsRemoteStorageSet())
- // {
- // AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
- // AliWarning("Remote CDB storage has been already set !");
- // AliWarning(Form("Ignoring the remote storage declared in AliSimulation: %s",fRemoteCDBUri.Data()));
- // AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
- // fRemoteCDBUri = "";
- // }
- // else {
- // AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
- // AliDebug(2, Form("Remote CDB storage is set to: %s",fRemoteCDBUri.Data()));
- // AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
- // man->SetRemoteStorage(fRemoteCDBUri);
- // }
-
// Now activate the detector specific CDB storage locations
for (Int_t i = 0; i < fSpecCDBUri.GetEntriesFast(); i++) {
TObject* obj = fSpecCDBUri[i];
AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
man->SetSpecificStorage(obj->GetName(), obj->GetTitle());
}
- man->Print();
+
}
//_____________________________________________________________________________
-void AliSimulation::SetDefaultStorage(const char* uri) {
-// Store the desired default CDB storage location
-// Activate it later within the Run() method
+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);
- fCDBUri = uri;
+ man->Print();
}
//_____________________________________________________________________________
-void AliSimulation::SetRemoteStorage(const char* uri) {
-// Store the desired remote CDB storage location
+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::SetDefaultStorage(const char* uri) {
+// Store the desired default CDB storage location
// Activate it later within the Run() method
-// Remote storage (the Grid storage) is used if it is activated
-// and if the object is not found in the default storage (the local cache)
- fRemoteCDBUri = uri;
+ fCDBUri = uri;
}
}
+//_____________________________________________________________________________
+void AliSimulation::SetRunNumber(Int_t run)
+{
+// sets run number
+// Activate it later within the Run() method
+
+ fRun = run;
+}
+
+//_____________________________________________________________________________
+void AliSimulation::SetSeed(Int_t seed)
+{
+// sets seed number
+// Activate it later within the Run() method
+
+ fSeed = seed;
+}
+
+//_____________________________________________________________________________
+Bool_t AliSimulation::SetRunNumberFromData()
+{
+ // 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;
+}
+
//_____________________________________________________________________________
void AliSimulation::SetConfigFile(const char* fileName)
{
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");
return kTRUE;
}
-
-//_____________________________________________________________________________
-Bool_t AliSimulation::SetRunNumber()
-{
- // Set the CDB manager run number
- // The run number is retrieved from gAlice
-
- 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;
-}
-
//_____________________________________________________________________________
void AliSimulation::MergeWith(const char* fileName, Int_t nSignalPerBkgrd)
{
// run the generation, simulation and digitization
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("")
- InitCDBStorage();
+ // 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;
}
-
-//QA
- AliQADataMakerSteer qas ;
- qas.Run(AliQA::kHITS);
-
- // Set run number in CDBManager (if it is not already set in RunSimulation)
- if (!SetRunNumber()) 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()) {
if(!MisalignGeometry()) if (fStopOnError) return kFALSE;
}
+
// hits -> summable digits
+ AliSysInfo::AddStamp("Start_sdigitization");
if (!fMakeSDigits.IsNull()) {
if (!RunSDigitization(fMakeSDigits)) if (fStopOnError) return kFALSE;
+
}
+ AliSysInfo::AddStamp("Stop_sdigitization");
- //QA
- qas.Reset() ;
- qas.Run(AliQA::kSDIGITS);
-
- // summable digits -> digits
+ 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)) {
if (fStopOnError) return kFALSE;
}
}
-
- //QA
- qas.Reset() ;
- qas.Run(AliQA::kDIGITS);
+
+
// 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
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( "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;
}
}
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
return kFALSE;
}
SetGAliceFile(runLoader->GetFileName());
-
+
// Misalign geometry
#if ROOT_VERSION_CODE < 331527
AliGeomManager::SetGeometry(gGeoManager);
}
AliInfo("running gAlice");
+ AliSysInfo::AddStamp("Start_simulation");
StdoutToAliInfo(StderrToAliError(
gAlice->Run(nEvents);
););
-
+ AliSysInfo::AddStamp("Stop_simulation");
delete runLoader;
-
return kTRUE;
}
Bool_t AliSimulation::RunSDigitization(const char* detectors)
{
// run the digitization and produce summable digits
-
+ static Int_t eventNr=0;
AliCodeTimerAuto("")
+ // initialize CDB storage, run number, set CDB lock
+ InitCDB();
+ if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE;
+ SetCDBLock();
+
AliRunLoader* runLoader = LoadRun();
if (!runLoader) return kFALSE;
if (IsSelected(det->GetName(), detStr)) {
AliInfo(Form("creating summable digits for %s", det->GetName()));
AliCodeTimerAuto(Form("creating summable digits for %s", det->GetName()));
-
det->Hits2SDigits();
+ AliSysInfo::AddStamp(Form("Digit_%s_%d",det->GetName(),eventNr), 0,1, eventNr);
}
}
detStr.Data()));
if (fStopOnError) return kFALSE;
}
-
+ eventNr++;
delete runLoader;
return kTRUE;
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;
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;
//_____________________________________________________________________________
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.
AliCodeTimerAuto("")
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);
+ }
+ }
}
}
}
//_____________________________________________________________________________
-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("")
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];
// Note the option -s. It is used in order to avoid
// the generation of SOR/EOR events.
- sprintf(command, "dateStream -s -D -o %s -# %d -C",
- dateFileName, runLoader->GetNumberOfEvents());
+ 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 < AliDAQ::kNDetectors; iDet++) {
for (Int_t iDDL = 0; iDDL < AliDAQ::NumberOfDdls(iDet); iDDL++) {
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));
}
//_____________________________________________________________________________
return NULL;
}
runLoader->LoadgAlice();
+ runLoader->LoadHeader();
gAlice = runLoader->GetAliRun();
if (!gAlice) {
AliError(Form("no gAlice object found in file %s",
return result;
}
+//_____________________________________________________________________________
Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* esdFileName)
{
//
StdoutToAliInfo(StderrToAliError(gAlice->Init(fConfigFileName.Data());););
//
// Initialize CDB
- InitCDBStorage();
- AliCDBManager* man = AliCDBManager::Instance();
- man->SetRun(0); // Should this come from rawdata header ?
+ InitCDB();
+ //AliCDBManager* man = AliCDBManager::Instance();
+ //man->SetRun(0); // Should this come from rawdata header ?
Int_t iDet;
//
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);
}
// init the HLT simulation
- if (fRunHLT.CompareTo("default")==0) fRunHLT="";
+ 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, fRunHLT.Data())))<0) {
+ 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
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 ;
+ detectorsw = qas.Run("ALL", AliQA::kHITS) ;
+// qas.Reset() ;
+ detectorsw += qas.Run(fMakeSDigits.Data(), AliQA::kSDIGITS) ;
+// qas.Reset() ;
+ detectorsw += qas.Run(fMakeDigits.Data(), AliQA::kDIGITS) ;
+// qas.Reset() ;
+ detectorsw += qas.Run(fMakeDigitsFromHits.Data(), AliQA::kDIGITS) ;
+
+ if ( detectorsw.IsNull() )
+ rv = kFALSE ;
+ return rv ;
+}
+
+//_____________________________________________________________________________
+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