X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=STEER%2FAliSimulation.cxx;h=84b6514e35f6dd9015e56306d2367d60feba4bf4;hb=5cb6ddd48d25c3aed7454e04854230334823c92b;hp=e6e718e561d618706ccb8e5cfcb84bbfc3701345;hpb=b856a511a4b888291cfab97fc6ddac72c07a745a;p=u%2Fmrichter%2FAliRoot.git diff --git a/STEER/AliSimulation.cxx b/STEER/AliSimulation.cxx index e6e718e561d..84b6514e35f 100644 --- a/STEER/AliSimulation.cxx +++ b/STEER/AliSimulation.cxx @@ -105,43 +105,58 @@ // // /////////////////////////////////////////////////////////////////////////////// -#include +#include +#include +#include #include #include -#include +#include #include -#include +#include +#include +#include -#include "AliCDBStorage.h" +#include "AliAlignObj.h" #include "AliCDBEntry.h" #include "AliCDBManager.h" -#include "AliGeomManager.h" -#include "AliAlignObj.h" +#include "AliCDBStorage.h" +#include "AliCTPRawData.h" #include "AliCentralTrigger.h" +#include "AliCentralTrigger.h" +#include "AliCodeTimer.h" #include "AliDAQ.h" #include "AliDigitizer.h" +#include "AliESDEvent.h" +#include "AliGRPObject.h" +#include "AliGenEventHeader.h" #include "AliGenerator.h" +#include "AliGeomManager.h" +#include "AliHLTSimulation.h" +#include "AliHeader.h" +#include "AliLego.h" +#include "AliLegoGenerator.h" #include "AliLog.h" +#include "AliMC.h" +#include "AliMagF.h" #include "AliModule.h" +#include "AliPDG.h" +#include "AliRawReaderDate.h" +#include "AliRawReaderFile.h" +#include "AliRawReaderRoot.h" #include "AliRun.h" #include "AliRunDigitizer.h" #include "AliRunLoader.h" #include "AliSimulation.h" +#include "AliSysInfo.h" #include "AliVertexGenFile.h" -#include "AliCentralTrigger.h" -#include "AliCTPRawData.h" -#include "AliRawReaderFile.h" -#include "AliESD.h" -#include "AliHeader.h" -#include "AliGenEventHeader.h" -#include "AliMC.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), @@ -156,8 +171,8 @@ AliSimulation::AliSimulation(const char* configFileName, const char* cdbUri, fWriteRawData(""), fRawDataFileName(""), fDeleteIntermediateFiles(kFALSE), + fWriteSelRawData(kFALSE), fStopOnError(kFALSE), - fNEvents(1), fConfigFileName(configFileName), fGAliceFileName("galice.root"), @@ -166,71 +181,34 @@ AliSimulation::AliSimulation(const char* configFileName, const char* cdbUri, fAlignObjArray(NULL), fUseBkgrdVertex(kTRUE), fRegionOfInterest(kFALSE), - fCDBUri(cdbUri), + fCDBUri(""), + fQARefUri(""), fSpecCDBUri(), - fEmbeddingFlag(kFALSE) + fRun(-1), + fSeed(0), + fInitCDBCalled(kFALSE), + fInitRunNumberCalled(kFALSE), + fSetRunNumberFromDataCalled(kFALSE), + fEmbeddingFlag(kFALSE), + fLego(NULL), + fQADetectors("ALL"), + fQATasks("ALL"), + fRunQA(kTRUE), + fEventSpecie(AliRecoParam::kDefault), + fWriteQAExpertData(kTRUE), + fRunHLT("default"), + fpHLT(NULL), + fWriteGRPEntry(kTRUE) { // create simulation object with default parameters fgInstance = this; SetGAliceFile("galice.root"); -} - -//_____________________________________________________________________________ -AliSimulation::AliSimulation(const AliSimulation& sim) : - TNamed(sim), - - fRunGeneration(sim.fRunGeneration), - fRunSimulation(sim.fRunSimulation), - fLoadAlignFromCDB(sim.fLoadAlignFromCDB), - fLoadAlObjsListOfDets(sim.fLoadAlObjsListOfDets), - 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), - fCDBUri(sim.fCDBUri), - fSpecCDBUri(), - fEmbeddingFlag(sim.fEmbeddingFlag) -{ -// 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; - 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; -} - -//_____________________________________________________________________________ -AliSimulation& AliSimulation::operator = (const AliSimulation& sim) -{ -// assignment operator - - this->~AliSimulation(); - new(this) AliSimulation(sim); - return *this; + +// for QA + AliQAManager * qam = AliQAManager::QAManager(AliQAv1::kSIMMODE) ; + qam->SetActiveDetectors(fQADetectors) ; + fQATasks = Form("%d %d %d", AliQAv1::kHITS, AliQAv1::kSDIGITS, AliQAv1::kDIGITS) ; + qam->SetTasks(fQATasks) ; } //_____________________________________________________________________________ @@ -249,6 +227,10 @@ AliSimulation::~AliSimulation() fSpecCDBUri.Delete(); if (fgInstance==this) fgInstance = 0; + + AliQAManager::QAManager()->ShowQA() ; + AliQAManager::Destroy() ; + AliCodeTimer::Instance()->Print(); } @@ -261,25 +243,77 @@ void AliSimulation::SetNumberOfEvents(Int_t nEvents) } //_____________________________________________________________________________ -void AliSimulation::InitCDBStorage() +void AliSimulation::InitQA() +{ + // 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; + + AliQAManager * qam = AliQAManager::QAManager(AliQAv1::kSIMMODE) ; + qam->SetActiveDetectors(fQADetectors) ; + fQATasks = Form("%d %d %d", AliQAv1::kHITS, AliQAv1::kSDIGITS, AliQAv1::kDIGITS) ; + qam->SetTasks(fQATasks) ; + if (fWriteQAExpertData) + qam->SetWriteExpert() ; + + if (qam->IsDefaultStorageSet()) { + AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); + AliWarning("Default QA reference storage has been already set !"); + AliWarning(Form("Ignoring the default storage declared in AliSimulation: %s",fQARefUri.Data())); + AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); + fQARefUri = qam->GetDefaultStorage()->GetURI(); + } else { + if (fQARefUri.Length() > 0) { + AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); + AliDebug(2, Form("Default QA reference storage is set to: %s", fQARefUri.Data())); + AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); + } else { + fQARefUri="local://$ALICE_ROOT/QARef"; + AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); + AliWarning("Default QA reference storage not yet set !!!!"); + AliWarning(Form("Setting it now to: %s", fQARefUri.Data())); + AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); + } + qam->SetDefaultStorage(fQARefUri); + } +} + +//_____________________________________________________________________________ +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("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); AliWarning("Default CDB storage has been already set !"); - AliWarning(Form("Ignoring the default storage declared in AliReconstruction: %s",fCDBUri.Data())); + 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/OCDB"; + AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); + AliWarning("Default CDB storage not yet set !!!!"); + AliWarning(Form("Setting it now to: %s", fCDBUri.Data())); + AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); + + } man->SetDefaultStorage(fCDBUri); } @@ -292,16 +326,60 @@ void AliSimulation::InitCDBStorage() AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"); man->SetSpecificStorage(obj->GetName(), obj->GetTitle()); } + +} + +//_____________________________________________________________________________ +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, Form("Setting CDB run number to: %d",fRun)); + } else { + fRun=0; + AliWarning(Form("Run number not yet set !!!! Setting it now to: %d", + fRun)); + } + man->SetRun(fRun); + man->Print(); + } //_____________________________________________________________________________ -void AliSimulation::SetDefaultStorage(const char* uri) { -// Store the desired default CDB storage location -// Activate it later within the Run() method +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 + fCDBUri = uri; + +} +//_____________________________________________________________________________ +void AliSimulation::SetQARefDefaultStorage(const char* uri) { + // Store the desired default CDB storage location + // Activate it later within the Run() method + + fQARefUri = uri; + AliQAv1::SetQARefStorage(fQARefUri.Data()) ; } //_____________________________________________________________________________ @@ -321,6 +399,67 @@ void AliSimulation::SetSpecificStorage(const char* calibType, const char* 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->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) { @@ -369,15 +508,25 @@ Bool_t AliSimulation::MisalignGeometry(AliRunLoader *runLoader) // 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 (gGeoManager) gGeoManager->Export("geometry.root"); + if(!gAlice->IsRootGeometry()) AliGeomManager::GetGeometry()->Export("geometry.root"); // Load alignment data from CDB and apply to geometry through AliGeomManager if(fLoadAlignFromCDB){ @@ -390,29 +539,23 @@ Bool_t AliSimulation::MisalignGeometry(AliRunLoader *runLoader) 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 += " "; + //add det to list of dets to be aligned from CDB + loadAlObjsListOfDets += det->GetName(); + loadAlObjsListOfDets += " "; } } // end loop over detectors - (AliGeomManager::Instance())->ApplyAlignObjsFromCDB(loadAlObjsListOfDets.Data()); + 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 (gGeoManager && gGeoManager->IsClosed()) { - if ((AliGeomManager::Instance())->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 (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; } } } @@ -432,28 +575,6 @@ Bool_t AliSimulation::MisalignGeometry(AliRunLoader *runLoader) if (delRunLoader) delete runLoader; - // Update the TGeoPhysicalNodes - gGeoManager->RefreshPhysicalNodes(); - - 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; } @@ -480,38 +601,95 @@ Bool_t AliSimulation::Run(Int_t nEvents) { // run the generation, simulation and digitization - InitCDBStorage(); + + AliCodeTimerAuto("",0) + AliSysInfo::AddStamp("Start_Run"); + + // Load run number and seed from environmental vars + ProcessEnvironmentVars(); + AliSysInfo::AddStamp("ProcessEnvironmentVars"); + gRandom->SetSeed(fSeed); + if (nEvents > 0) fNEvents = nEvents; + // create and setup the HLT instance + if (!fRunHLT.IsNull() && !CreateHLT()) { + if (fStopOnError) return kFALSE; + // disable HLT + fRunHLT=""; + } + // 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; + AliSysInfo::AddStamp("RunSimulation"); + + // initialize CDB storage from external environment + // (either CDB manager or AliSimulation setters), + // if not already done in RunSimulation() + InitCDB(); + AliSysInfo::AddStamp("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 (!gGeoManager) { - TGeoManager::Import("geometry.root"); - if (!gGeoManager) if (fStopOnError) return kFALSE; - // Initialize the geometry manager (if not already done) + if (!AliGeomManager::GetGeometry()) { + // Initialize the geometry manager + AliGeomManager::LoadGeometry("geometry.root"); + AliSysInfo::AddStamp("GetGeometry"); + + +// // 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; } + AliSysInfo::AddStamp("MissalignGeometry"); + // 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)) { @@ -525,60 +703,237 @@ Bool_t AliSimulation::Run(Int_t nEvents) } } + AliSysInfo::AddStamp("Hits2Digits"); + + // digits -> trigger - if (!RunTrigger(fMakeTrigger)) { + if (!RunTrigger(fMakeTrigger,fMakeDigits)) { if (fStopOnError) return kFALSE; } + AliSysInfo::AddStamp("RunTrigger"); + + // digits -> raw data if (!fWriteRawData.IsNull()) { if (!WriteRawData(fWriteRawData, fRawDataFileName, - fDeleteIntermediateFiles)) { + fDeleteIntermediateFiles,fWriteSelRawData)) { + if (fStopOnError) return kFALSE; + } + } + + AliSysInfo::AddStamp("WriteRaw"); + + // run HLT simulation on simulated digit data if raw data is not + // simulated, otherwise its called as part of WriteRawData + if (!fRunHLT.IsNull() && fWriteRawData.IsNull()) { + if (!RunHLT()) { if (fStopOnError) return kFALSE; } } + AliSysInfo::AddStamp("RunHLT"); + + //QA + if (fRunQA) { + Bool_t rv = RunQA() ; + if (!rv) + if (fStopOnError) + return kFALSE ; + } + + AliSysInfo::AddStamp("RunQA"); + + // Cleanup of CDB manager: cache and active storages! + AliCDBManager::Instance()->ClearCache(); + + return kTRUE; +} + +//_______________________________________________________________________ +Bool_t AliSimulation::RunLego(const char *setup, Int_t nc1, Float_t c1min, + Float_t c1max,Int_t nc2,Float_t c2min,Float_t c2max, + Float_t rmin,Float_t rmax,Float_t zmax, AliLegoGenerator* gener, Int_t nev) +{ + // + // Generates lego plots of: + // - radiation length map phi vs theta + // - radiation length map phi vs eta + // - interaction length map + // - g/cm2 length map + // + // ntheta bins in theta, eta + // themin minimum angle in theta (degrees) + // themax maximum angle in theta (degrees) + // nphi bins in phi + // phimin minimum angle in phi (degrees) + // phimax maximum angle in phi (degrees) + // rmin minimum radius + // rmax maximum radius + // + // + // The number of events generated = ntheta*nphi + // run input parameters in macro setup (default="Config.C") + // + // Use macro "lego.C" to visualize the 3 lego plots in spherical coordinates + //Begin_Html + /* + + */ + //End_Html + //Begin_Html + /* + + */ + //End_Html + //Begin_Html + /* + + */ + //End_Html + // + +// run the generation and simulation + + AliCodeTimerAuto("",0) + + // 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; + } + 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())); + + // Number of events + if (nev == -1) nev = nc1 * nc2; + + // check if initialisation has been done + // If runloader has been initialized, set the number of events per file to nc1 * nc2 + + // Set new generator + if (!gener) gener = new AliLegoGenerator(); + // + // Configure Generator + + gener->SetRadiusRange(rmin, rmax); + gener->SetZMax(zmax); + gener->SetCoor1Range(nc1, c1min, c1max); + gener->SetCoor2Range(nc2, c2min, c2max); + + + //Create Lego object + fLego = new AliLego("lego",gener); + + //__________________________________________________________________________ + + gAlice->Announce(); + + gROOT->LoadMacro(setup); + gInterpreter->ProcessLine(gAlice->GetConfigFunction()); + + if(AliCDBManager::Instance()->GetRun() >= 0) { + SetRunNumber(AliCDBManager::Instance()->GetRun()); + } else { + AliWarning("Run number not initialized!!"); + } + + AliRunLoader::Instance()->CdGAFile(); + + AliPDG::AddParticlesToPdgDataBase(); + + gMC->SetMagField(TGeoGlobalMagField::Instance()->GetField()); + + gAlice->GetMCApp()->Init(); + + + //Must be here because some MCs (G4) adds detectors here and not in Config.C + gAlice->InitLoaders(); + AliRunLoader::Instance()->MakeTree("E"); + + // + // Save stuff at the beginning of the file to avoid file corruption + AliRunLoader::Instance()->CdGAFile(); + gAlice->Write(); + + //Save current generator + AliGenerator *gen=gAlice->GetMCApp()->Generator(); + gAlice->GetMCApp()->ResetGenerator(gener); + //Prepare MC for Lego Run + gMC->InitLego(); + + //Run Lego Object + + + AliRunLoader::Instance()->SetNumberOfEventsPerFile(nev); + gMC->ProcessRun(nev); + + // End of this run, close files + FinishRun(); + // Restore current generator + gAlice->GetMCApp()->ResetGenerator(gen); + // Delete Lego Object + delete fLego; + 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("",0) + // 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 (gAlice->GetTriggerDescriptor() != "") { - des = gAlice->GetTriggerDescriptor(); - } - else { - AliWarning("No trigger descriptor is specified. Skipping the trigger simulation..."); - return kTRUE; + if (trconfiguration.IsNull()) { + if (strcmp(gAlice->GetTriggerDescriptor(),"")) { + 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; } } - AliInfo(Form("Execution time: R:%.2fs C:%.2fs", - stopwatch.RealTime(),stopwatch.CpuTime())); - delete runLoader; return kTRUE; @@ -601,9 +956,17 @@ Bool_t AliSimulation::RunSimulation(Int_t nEvents) { // run the generation and simulation - TStopwatch stopwatch; - stopwatch.Start(); + AliCodeTimerAuto("",0) + // initialize CDB storage and run number from external environment + // (either CDB manager or AliSimulation setters) + AliSysInfo::AddStamp("RunSimulation_Begin"); + InitCDB(); + AliSysInfo::AddStamp("RunSimulation_InitCDB"); + InitRunNumber(); + SetCDBLock(); + AliSysInfo::AddStamp("RunSimulation_SetCDBLock"); + if (!gAlice) { AliError("no gAlice object. Restart aliroot and try again."); return kFALSE; @@ -615,39 +978,91 @@ Bool_t AliSimulation::RunSimulation(Int_t nEvents) AliInfo(Form("initializing gAlice with config file %s", fConfigFileName.Data())); - StdoutToAliInfo(StderrToAliError( - gAlice->Init(fConfigFileName.Data()); - );); + // + // Initialize ALICE Simulation run + // + + gAlice->Announce(); + + gROOT->LoadMacro(fConfigFileName.Data()); + gInterpreter->ProcessLine(gAlice->GetConfigFunction()); + AliSysInfo::AddStamp("RunSimulation_Config"); + + if(AliCDBManager::Instance()->GetRun() >= 0) { + AliRunLoader::Instance()->SetRunNumber(AliCDBManager::Instance()->GetRun()); + AliRunLoader::Instance()->SetNumberOfEventsPerRun(fNEvents); + } else { + AliWarning("Run number not initialized!!"); + } + + AliRunLoader::Instance()->CdGAFile(); + + AliPDG::AddParticlesToPdgDataBase(); + + gMC->SetMagField(TGeoGlobalMagField::Instance()->GetField()); + AliSysInfo::AddStamp("RunSimulation_GetField"); + + gAlice->GetMCApp()->Init(); + AliSysInfo::AddStamp("RunSimulation_InitMCApp"); + + //Must be here because some MCs (G4) adds detectors here and not in Config.C + gAlice->InitLoaders(); + AliRunLoader::Instance()->MakeTree("E"); + AliRunLoader::Instance()->LoadKinematics("RECREATE"); + AliRunLoader::Instance()->LoadTrackRefs("RECREATE"); + AliRunLoader::Instance()->LoadHits("all","RECREATE"); + // + // Save stuff at the beginning of the file to avoid file corruption + AliRunLoader::Instance()->CdGAFile(); + gAlice->Write(); + gAlice->SetEventNrInRun(-1); //important - we start Begin event from increasing current number in run + AliSysInfo::AddStamp("RunSimulation_InitLoaders"); + //___________________________________________________________________________________________ + // Get the trigger descriptor string // Either from AliSimulation or from // gAlice if (fMakeTrigger.IsNull()) { - if (gAlice->GetTriggerDescriptor() != "") + if (strcmp(gAlice->GetTriggerDescriptor(),"")) fMakeTrigger = gAlice->GetTriggerDescriptor(); } else gAlice->SetTriggerDescriptor(fMakeTrigger.Data()); + AliSysInfo::AddStamp("RunSimulation_TriggerDescriptor"); // Set run number in CDBManager - AliCDBManager::Instance()->SetRun(gAlice->GetRunNumber()); AliInfo(Form("Run number: %d",AliCDBManager::Instance()->GetRun())); - AliRunLoader* runLoader = gAlice->GetRunLoader(); + AliRunLoader* runLoader = AliRunLoader::Instance(); if (!runLoader) { AliError(Form("gAlice has no run loader object. " "Check your config file: %s", fConfigFileName.Data())); return kFALSE; } SetGAliceFile(runLoader->GetFileName()); - + // Misalign geometry -// MisalignGeometry(runLoader); - - // Export (mis)aligned geometry - if (gGeoManager) gGeoManager->Export("misaligned_geometry.root"); +#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); + AliSysInfo::AddStamp("RunSimulation_MisalignGeometry"); +#endif -// AliRunLoader* runLoader = gAlice->GetRunLoader(); +// AliRunLoader* runLoader = AliRunLoader::Instance(); // if (!runLoader) { // AliError(Form("gAlice has no run loader object. " // "Check your config file: %s", fConfigFileName.Data())); @@ -655,11 +1070,17 @@ Bool_t AliSimulation::RunSimulation(Int_t nEvents) // } // SetGAliceFile(runLoader->GetFileName()); - if (!gAlice->Generator()) { + if (!gAlice->GetMCApp()->Generator()) { AliError(Form("gAlice has no generator object. " "Check your config file: %s", fConfigFileName.Data())); return kFALSE; } + + // Write GRP entry corresponding to the setting found in Cofig.C + if (fWriteGRPEntry) + WriteGRPEntry(); + AliSysInfo::AddStamp("RunSimulation_WriteGRP"); + if (nEvents <= 0) nEvents = fNEvents; // get vertex from background file in case of merging @@ -672,11 +1093,11 @@ Bool_t AliSimulation::RunSimulation(Int_t nEvents) "file %s with nSignalPerBackground = %d", fileName, signalPerBkgrd)); AliVertexGenFile* vtxGen = new AliVertexGenFile(fileName, signalPerBkgrd); - gAlice->Generator()->SetVertexGenerator(vtxGen); + gAlice->GetMCApp()->Generator()->SetVertexGenerator(vtxGen); } if (!fRunSimulation) { - gAlice->Generator()->SetTrackingFlag(0); + gAlice->GetMCApp()->Generator()->SetTrackingFlag(0); } // set the number of events per file for given detectors and data types @@ -707,14 +1128,18 @@ Bool_t AliSimulation::RunSimulation(Int_t nEvents) } AliInfo("running gAlice"); - StdoutToAliInfo(StderrToAliError( - gAlice->Run(nEvents); - );); + AliSysInfo::AddStamp("Start_ProcessRun"); - delete runLoader; + // Create the Root Tree with one branch per detector + //Hits moved to begin event -> now we are crating separate tree for each event + + gMC->ProcessRun(nEvents); - AliInfo(Form("Execution time: R:%.2fs C:%.2fs", - stopwatch.RealTime(),stopwatch.CpuTime())); + // End of this run, close files + if(nEvents>0) FinishRun(); + + AliSysInfo::AddStamp("Stop_ProcessRun"); + delete runLoader; return kTRUE; } @@ -723,10 +1148,14 @@ Bool_t AliSimulation::RunSimulation(Int_t nEvents) Bool_t AliSimulation::RunSDigitization(const char* detectors) { // run the digitization and produce summable digits + static Int_t eventNr=0; + AliCodeTimerAuto("",0) ; - 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; @@ -737,11 +1166,10 @@ Bool_t AliSimulation::RunSDigitization(const char* detectors) if (!det || !det->IsActive()) continue; if (IsSelected(det->GetName(), detStr)) { AliInfo(Form("creating summable digits for %s", det->GetName())); - TStopwatch stopwatchDet; - stopwatchDet.Start(); + AliCodeTimerStart(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())); + AliCodeTimerStop(Form("creating summable digits for %s", det->GetName())); + AliSysInfo::AddStamp(Form("Digit_%s_%d",det->GetName(),eventNr), 0,1, eventNr); } } @@ -750,12 +1178,9 @@ Bool_t AliSimulation::RunSDigitization(const char* detectors) detStr.Data())); if (fStopOnError) return kFALSE; } - + eventNr++; delete runLoader; - AliInfo(Form("Execution time: R:%.2fs C:%.2fs", - stopwatch.RealTime(),stopwatch.CpuTime())); - return kTRUE; } @@ -766,11 +1191,15 @@ Bool_t AliSimulation::RunDigitization(const char* detectors, { // run the digitization and produce digits from sdigits - TStopwatch stopwatch; - stopwatch.Start(); + AliCodeTimerAuto("",0) - while (AliRunLoader::GetRunLoader()) delete AliRunLoader::GetRunLoader(); - if (gAlice) delete gAlice; + // initialize CDB storage, run number, set CDB lock + InitCDB(); + if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE; + SetCDBLock(); + + delete AliRunLoader::Instance(); + delete gAlice; gAlice = NULL; Int_t nStreams = 1; @@ -820,9 +1249,6 @@ Bool_t AliSimulation::RunDigitization(const char* detectors, delete manager; - AliInfo(Form("Execution time: R:%.2fs C:%.2fs", - stopwatch.RealTime(),stopwatch.CpuTime())); - return kTRUE; } @@ -831,9 +1257,13 @@ Bool_t AliSimulation::RunHitsDigitization(const char* detectors) { // run the digitization and produce digits from hits - TStopwatch stopwatch; - stopwatch.Start(); + AliCodeTimerAuto("",0) + // initialize CDB storage, run number, set CDB lock + InitCDB(); + if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE; + SetCDBLock(); + AliRunLoader* runLoader = LoadRun("READ"); if (!runLoader) return kFALSE; @@ -854,20 +1284,14 @@ Bool_t AliSimulation::RunHitsDigitization(const char* detectors) 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 - - 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. @@ -877,21 +1301,46 @@ Bool_t AliSimulation::WriteRawData(const char* detectors, // 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(); - - if (!WriteRawFiles(detectors)) { + AliCodeTimerAuto("",0) + AliSysInfo::AddStamp("WriteRawData_Start"); + + TString detStr = detectors; + if (!WriteRawFiles(detStr.Data())) { if (fStopOnError) return kFALSE; } + AliSysInfo::AddStamp("WriteRawFiles"); + + // run HLT simulation on simulated DDL raw files + // and produce HLT ddl raw files to be included in date/root file + // bugfix 2009-06-26: the decision whether to write HLT raw data + // is taken in RunHLT. Here HLT always needs to be run in order to + // create HLT digits, unless its switched off. This is due to the + // special placement of the HLT between the generation of DDL files + // and conversion to DATE/Root file. + detStr.ReplaceAll("HLT", ""); + if (!fRunHLT.IsNull()) { + if (!RunHLT()) { + if (fStopOnError) return kFALSE; + } + } + AliSysInfo::AddStamp("WriteRawData_RunHLT"); TString dateFileName(fileName); 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; } + AliSysInfo::AddStamp("ConvertRawFilesToDate"); if (deleteIntermediateFiles) { AliRunLoader* runLoader = LoadRun("READ"); if (runLoader) for (Int_t iEvent = 0; @@ -906,15 +1355,23 @@ Bool_t AliSimulation::WriteRawData(const char* detectors, if (!ConvertDateToRoot(dateFileName, fileName)) { if (fStopOnError) return kFALSE; } + AliSysInfo::AddStamp("ConvertDateToRoot"); 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; } @@ -923,6 +1380,8 @@ Bool_t AliSimulation::WriteRawFiles(const char* detectors) { // convert the digits to raw data DDL files + AliCodeTimerAuto("",0) + AliRunLoader* runLoader = LoadRun("READ"); if (!runLoader) return kFALSE; @@ -939,7 +1398,15 @@ Bool_t AliSimulation::WriteRawFiles(const char* detectors) if (fStopOnError) return kFALSE; else continue; } + ofstream runNbFile(Form("run%u",runLoader->GetHeader()->GetRun())); + runNbFile.close(); + TString detStr = detectors; + if (IsSelected("HLT", detStr)) { + // Do nothing. "HLT" will be removed from detStr and HLT raw + // data files are generated in RunHLT. + } + TObjArray* detArray = runLoader->GetAliRun()->Detectors(); for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) { AliModule* det = (AliModule*) detArray->At(iDet); @@ -962,14 +1429,20 @@ Bool_t AliSimulation::WriteRawFiles(const char* detectors) } 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("",0) + char* path = gSystem->Which(gSystem->Getenv("PATH"), "dateStream"); if (!path) { AliError("the program dateStream was not found"); @@ -982,15 +1455,43 @@ Bool_t AliSimulation::ConvertRawFilesToDate(const char* dateFileName) 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"); + if (!pipe) { + AliError(Form("Cannot execute command: %s",command)); + return kFALSE; + } + + Int_t selEvents = 0; for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) { - fprintf(pipe, "GDC\n"); + + UInt_t detectorPattern = 0; + runLoader->GetEvent(iEvent); + if (!runLoader->LoadTrigger()) { + AliCentralTrigger *aCTP = runLoader->GetTrigger(); + detectorPattern = aCTP->GetClusterMask(); + // Check if the event was triggered by CTP + if (selrawdata) { + if (aCTP->GetClassMask()) selEvents++; + } + } + else { + AliWarning("No trigger can be loaded! Some fields in the event header will be empty !"); + if (selrawdata) { + AliWarning("No trigger can be loaded! Writing of selected raw data is abandoned !"); + selrawdata = kFALSE; + } + } + + fprintf(pipe, "GDC DetectorPattern %u\n", detectorPattern); Float_t ldc = 0; Int_t prevLDC = -1; @@ -1025,8 +1526,72 @@ Bool_t AliSimulation::ConvertRawFilesToDate(const char* dateFileName) 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 + UInt_t detectorPattern = 0; + TString detClust; + runLoader->GetEvent(iEvent); + if (!runLoader->LoadTrigger()) { + AliCentralTrigger *aCTP = runLoader->GetTrigger(); + if (aCTP->GetClassMask() == 0) continue; + detectorPattern = aCTP->GetClusterMask(); + detClust = AliDAQ::ListOfTriggeredDetectors(detectorPattern); + AliInfo(Form("List of detectors to be read out: %s",detClust.Data())); + } + + fprintf(pipe2, "GDC DetectorPattern %u\n", detectorPattern); + 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)); } //_____________________________________________________________________________ @@ -1039,7 +1604,7 @@ Bool_t AliSimulation::ConvertDateToRoot(const char* dateFileName, const Int_t kDBSize = 2000000000; const Int_t kTagDBSize = 1000000000; const Bool_t kFilter = kFALSE; - const Int_t kCompression = 0; + const Int_t kCompression = 1; char* path = gSystem->Which(gSystem->Getenv("PATH"), "alimdc"); if (!path) { @@ -1088,7 +1653,7 @@ AliRunLoader* AliSimulation::LoadRun(const char* mode) const { // delete existing run loaders, open a new one and load gAlice - while (AliRunLoader::GetRunLoader()) delete AliRunLoader::GetRunLoader(); + delete AliRunLoader::Instance(); AliRunLoader* runLoader = AliRunLoader::Open(fGAliceFileName.Data(), AliConfig::GetDefaultEventFolderName(), mode); @@ -1097,6 +1662,7 @@ AliRunLoader* AliSimulation::LoadRun(const char* mode) const return NULL; } runLoader->LoadgAlice(); + runLoader->LoadHeader(); gAlice = runLoader->GetAliRun(); if (!gAlice) { AliError(Form("no gAlice object found in file %s", @@ -1191,6 +1757,7 @@ Bool_t AliSimulation::IsSelected(TString detName, TString& detectors) const return result; } +//_____________________________________________________________________________ Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* esdFileName) { // @@ -1200,45 +1767,77 @@ Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* e // // 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; - } + 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())); + + gAlice->Announce(); + + gROOT->LoadMacro(fConfigFileName.Data()); + gInterpreter->ProcessLine(gAlice->GetConfigFunction()); + + if(AliCDBManager::Instance()->GetRun() >= 0) { + SetRunNumber(AliCDBManager::Instance()->GetRun()); + } else { + AliWarning("Run number not initialized!!"); + } + + AliRunLoader::Instance()->CdGAFile(); - AliInfo(Form("initializing gAlice with config file %s",fConfigFileName.Data())); - StdoutToAliInfo(StderrToAliError(gAlice->Init(fConfigFileName.Data()););); + AliPDG::AddParticlesToPdgDataBase(); + + gMC->SetMagField(TGeoGlobalMagField::Instance()->GetField()); + + gAlice->GetMCApp()->Init(); + + //Must be here because some MCs (G4) adds detectors here and not in Config.C + gAlice->InitLoaders(); + AliRunLoader::Instance()->MakeTree("E"); + AliRunLoader::Instance()->LoadKinematics("RECREATE"); + AliRunLoader::Instance()->LoadTrackRefs("RECREATE"); + AliRunLoader::Instance()->LoadHits("all","RECREATE"); + // + // Save stuff at the beginning of the file to avoid file corruption + AliRunLoader::Instance()->CdGAFile(); + gAlice->Write(); // // 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; // // Get the runloader - AliRunLoader* runLoader = gAlice->GetRunLoader(); + AliRunLoader* runLoader = AliRunLoader::Instance(); // // 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); - } - } + AliESDEvent* esd = new AliESDEvent(); + esdFile->GetObject("esdTree", treeESD); + if (treeESD) esd->ReadFromTree(treeESD); + // // Create the RawReader - AliRawReaderFile* rawReader = new AliRawReaderFile(rawDirectory); + 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); + } +// if (!fEquipIdMap.IsNull() && fRawReader) +// fRawReader->LoadEquipmentIdsMap(fEquipIdMap); // // Get list of detectors TObjArray* detArray = runLoader->GetAliRun()->Detectors(); @@ -1246,6 +1845,7 @@ Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* e // Get Header AliHeader* header = runLoader->GetHeader(); // + TString detStr = fMakeSDigits; // Event loop Int_t nev = 0; while(kTRUE) { @@ -1254,14 +1854,18 @@ Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* e // Detector loop for (iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) { AliModule* det = (AliModule*) detArray->At(iDet); - AliInfo(Form("Calling Raw2SDigits for %s\n", det->GetName())); - det->Raw2SDigits(rawReader); - rawReader->Reset(); + 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) { + if (treeESD) { treeESD->GetEvent(nev); const AliESDVertex* esdVertex = esd->GetPrimaryVertex(); Double_t position[3]; @@ -1291,3 +1895,419 @@ Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* e return kTRUE; } + +//_____________________________________________________________________________ +void AliSimulation::FinishRun() +{ + // + // Called at the end of the run. + // + + if(IsLegoRun()) + { + AliDebug(1, "Finish Lego"); + AliRunLoader::Instance()->CdGAFile(); + fLego->FinishRun(); + } + + // Clean detector information + TIter next(gAlice->Modules()); + AliModule *detector; + while((detector = dynamic_cast(next()))) { + AliDebug(2, Form("%s->FinishRun()", detector->GetName())); + detector->FinishRun(); + } + + AliDebug(1, "AliRunLoader::Instance()->WriteHeader(OVERWRITE)"); + AliRunLoader::Instance()->WriteHeader("OVERWRITE"); + + // Write AliRun info and all detectors parameters + AliRunLoader::Instance()->CdGAFile(); + gAlice->Write(0,TObject::kOverwrite);//write AliRun + AliRunLoader::Instance()->Write(0,TObject::kOverwrite);//write RunLoader itself + + if(gAlice->GetMCApp()) gAlice->GetMCApp()->FinishRun(); + AliRunLoader::Instance()->Synchronize(); +} + +//_____________________________________________________________________________ +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::CreateHLT() +{ + // Init the HLT simulation. + // The function loads the library and creates the instance of AliHLTSimulation. + // the main reason for the decoupled creation is to set the transient OCDB + // objects before the OCDB is locked + + // 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) { + AliWarning(Form("%s version does not match: compiled for version %d, loaded %d", ALIHLTSIMULATION_LIBRARY, ALIHLTSIMULATION_LIBRARY_VERSION, fctVersion())); + } + + // print compile info + typedef void (*CompileInfo)( const char*& date, const char*& time); + CompileInfo fctInfo=(CompileInfo)gSystem->DynFindSymbol(ALIHLTSIMULATION_LIBRARY, "CompileInfo"); + if (fctInfo) { + const char* date=""; + const 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)); + if (fctCreate==NULL || (fpHLT=(fctCreate()))==NULL) { + AliError(Form("can not create instance of HLT simulation (creator %p)", fctCreate)); + return kFALSE; + } + + TString specObjects; + for (Int_t i = 0; i < fSpecCDBUri.GetEntriesFast(); i++) { + if (specObjects.Length()>0) specObjects+=" "; + specObjects+=fSpecCDBUri[i]->GetName(); + } + + AliHLTSimulationSetup_t fctSetup=(AliHLTSimulationSetup_t)(gSystem->DynFindSymbol(ALIHLTSIMULATION_LIBRARY, ALIHLTSIMULATION_SETUP)); + if (fctSetup==NULL || fctSetup(fpHLT, this, specObjects.Data())<0) { + AliWarning(Form("failed to setup HLT simulation (function %p)", fctSetup)); + } + + return kTRUE; +} + +//_____________________________________________________________________________ +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-current/ + // http://web.ift.uib.no/~kjeks/doc/alice-hlt-current/classAliHLTModuleAgent.html + // + // The libraries to be loaded can be specified as an option. + // 
+  // AliSimulation sim;
+  // sim.SetRunHLT("libAliHLTSample.so");
+  //
+ // will only load libAliHLTSample.so + + // Other available options: + // \li loglevel=level
+ // logging level for this processing + // \li alilog=off + // disable redirection of log messages to AliLog class + // \li config=macro + // configuration macro + // \li chains=configuration + // comma separated list of configurations to be run during simulation + // \li rawfile=file + // source for the RawReader to be created, the default is ./ if + // raw data is simulated + + int iResult=0; + + if (!fpHLT && !CreateHLT()) { + return kFALSE; + } + AliHLTSimulation* pHLT=fpHLT; + + AliRunLoader* pRunLoader = LoadRun("READ"); + if (!pRunLoader) return kFALSE; + + // initialize CDB storage, run number, set CDB lock + // thats for the case of running HLT simulation without all the other steps + // multiple calls are handled by the function, so we can just call + InitCDB(); + if (!SetRunNumberFromData()) if (fStopOnError) return kFALSE; + SetCDBLock(); + + // init the HLT simulation + TString options; + if (fRunHLT.CompareTo("default")!=0) options=fRunHLT; + TString detStr = fWriteRawData; + if (!IsSelected("HLT", detStr)) { + options+=" writerawfiles="; + } else { + options+=" writerawfiles=HLT"; + } + + if (!detStr.IsNull() && !options.Contains("rawfile=")) { + // as a matter of fact, HLT will run reconstruction and needs the RawReader + // in order to get detector data. By default, RawReaderFile is used to read + // the already simulated ddl files. Date and Root files from the raw data + // are generated after the HLT simulation. + options+=" rawfile=./"; + } + + 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 + + if(!gAlice) return kFALSE; + AliQAManager::QAManager()->SetRunLoader(AliRunLoader::Instance()) ; + + TString detectorsw("") ; + Bool_t rv = kTRUE ; + AliQAManager::QAManager()->SetEventSpecie(fEventSpecie) ; + detectorsw = AliQAManager::QAManager()->Run(fQADetectors.Data()) ; + 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") ){ + TString tmp = Form("%s %s", fMakeDigits.Data(), fMakeDigitsFromHits.Data()) ; + Int_t minus = fQADetectors.Last('-') ; + TString toKeep = Form("%s %s", fMakeDigits.Data(), fMakeDigitsFromHits.Data()) ; + TString toRemove("") ; + while (minus >= 0) { + toRemove = fQADetectors(minus+1, fQADetectors.Length()) ; + toRemove = toRemove.Strip() ; + toKeep.ReplaceAll(toRemove, "") ; + fQADetectors.ReplaceAll(Form("-%s", toRemove.Data()), "") ; + minus = fQADetectors.Last('-') ; + } + fQADetectors = toKeep ; + } + fQATasks = detAndAction(colon+1, detAndAction.Sizeof() ) ; + if (fQATasks.Contains("ALL") ) { + fQATasks = Form("%d %d %d", AliQAv1::kHITS, AliQAv1::kSDIGITS, AliQAv1::kDIGITS) ; + } else { + fQATasks.ToUpper() ; + TString tempo("") ; + if ( fQATasks.Contains("HIT") ) + tempo = Form("%d ", AliQAv1::kHITS) ; + if ( fQATasks.Contains("SDIGIT") ) + tempo += Form("%d ", AliQAv1::kSDIGITS) ; + if ( fQATasks.Contains("DIGIT") ) + tempo += Form("%d ", AliQAv1::kDIGITS) ; + fQATasks = tempo ; + if (fQATasks.IsNull()) { + AliInfo("No QA requested\n") ; + fRunQA = kFALSE ; + return kTRUE ; + } + } + TString tempo(fQATasks) ; + tempo.ReplaceAll(Form("%d", AliQAv1::kHITS), AliQAv1::GetTaskName(AliQAv1::kHITS)) ; + tempo.ReplaceAll(Form("%d", AliQAv1::kSDIGITS), AliQAv1::GetTaskName(AliQAv1::kSDIGITS)) ; + tempo.ReplaceAll(Form("%d", AliQAv1::kDIGITS), AliQAv1::GetTaskName(AliQAv1::kDIGITS)) ; + AliInfo( Form("QA will be done on \"%s\" for \"%s\"\n", fQADetectors.Data(), tempo.Data()) ) ; + fRunQA = kTRUE ; + AliQAManager::QAManager()->SetActiveDetectors(fQADetectors) ; + AliQAManager::QAManager()->SetTasks(fQATasks) ; + for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++) + AliQAManager::QAManager()->SetWriteExpert(AliQAv1::DETECTORINDEX_t(det)) ; + + 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)); +} + +//--------------------------------------------------------------------- +void AliSimulation::WriteGRPEntry() +{ + // Get the necessary information from galice (generator, trigger etc) and + // write a GRP entry corresponding to the settings in the Config.C used + // note that Hall probes and Cavern and Surface Atmos pressures are not simulated. + + + AliInfo("Writing global run parameters entry into the OCDB"); + + AliGRPObject* grpObj = new AliGRPObject(); + + grpObj->SetRunType("PHYSICS"); + grpObj->SetTimeStart(0); + TDatime curtime; + grpObj->SetTimeStart(0); + grpObj->SetTimeEnd(curtime.Convert()); + grpObj->SetBeamEnergyIsSqrtSHalfGeV(); // new format of GRP: store sqrt(s)/2 in GeV + + const AliGenerator *gen = gAlice->GetMCApp()->Generator(); + if (gen) { + grpObj->SetBeamEnergy(gen->GetEnergyCMS()/2); + TString projectile; + Int_t a,z; + gen->GetProjectile(projectile,a,z); + TString target; + gen->GetTarget(target,a,z); + TString beamType = projectile + "-" + target; + beamType.ReplaceAll(" ",""); + if (!beamType.CompareTo("-")) { + grpObj->SetBeamType("UNKNOWN"); + } + else { + grpObj->SetBeamType(beamType); + // Heavy ion run, the event specie is set to kHighMult + fEventSpecie = AliRecoParam::kHighMult; + if ((strcmp(beamType,"p-p") == 0) || + (strcmp(beamType,"p-") == 0) || + (strcmp(beamType,"-p") == 0) || + (strcmp(beamType,"P-P") == 0) || + (strcmp(beamType,"P-") == 0) || + (strcmp(beamType,"-P") == 0)) { + // Proton run, the event specie is set to kLowMult + fEventSpecie = AliRecoParam::kLowMult; + } + } + } else { + AliWarning("Unknown beam type and energy! Setting energy to 0"); + grpObj->SetBeamEnergy(0); + grpObj->SetBeamType("UNKNOWN"); + } + + UInt_t detectorPattern = 0; + Int_t nDets = 0; + TObjArray *detArray = gAlice->Detectors(); + for (Int_t iDet = 0; iDet < AliDAQ::kNDetectors-1; iDet++) { + if (detArray->FindObject(AliDAQ::OfflineModuleName(iDet))) { + detectorPattern |= (1 << iDet); + nDets++; + } + } + // CTP + if (!fMakeTrigger.IsNull() || strcmp(gAlice->GetTriggerDescriptor(),"")) + detectorPattern |= (1 << AliDAQ::DetectorID("TRG")); + + // HLT + if (!fRunHLT.IsNull()) + detectorPattern |= (1 << AliDAQ::kHLTId); + + grpObj->SetNumberOfDetectors((Char_t)nDets); + grpObj->SetDetectorMask((Int_t)detectorPattern); + grpObj->SetLHCPeriod("LHC08c"); + grpObj->SetLHCState("STABLE_BEAMS"); + // + AliMagF *field = (AliMagF*)TGeoGlobalMagField::Instance()->GetField(); + Float_t solenoidField = field ? TMath::Abs(field->SolenoidField()) : 0; + + Float_t factorSol = field ? field->GetFactorSol() : 0; + Float_t currentSol = TMath::Abs(factorSol)>1E-6 ? + TMath::Nint(TMath::Abs(solenoidField/factorSol))/5.*30000.*TMath::Abs(factorSol) : 0; + // + Float_t factorDip = field ? field->GetFactorDip() : 0; + Float_t currentDip = 6000.*TMath::Abs(factorDip); + // + grpObj->SetL3Current(currentSol,(AliGRPObject::Stats)0); + grpObj->SetDipoleCurrent(currentDip,(AliGRPObject::Stats)0); + grpObj->SetL3Polarity(factorSol>0 ? 0:1); + grpObj->SetDipolePolarity(factorDip>0 ? 0:1); + grpObj->SetUniformBMap(field->IsUniform()); // for special MC with k5kGUniform map + grpObj->SetPolarityConventionLHC(); // LHC convention +/+ current -> -/- field main components + // + grpObj->SetCavernTemperature(0,(AliGRPObject::Stats)0); + + //grpMap->Add(new TObjString("fCavernPressure"),new TObjString("0")); ---> not inserted in simulation with the new object, since it is now an AliDCSSensor + + // Now store the entry in OCDB + AliCDBManager* man = AliCDBManager::Instance(); + + AliCDBId id("GRP/GRP/Data", man->GetRun(), man->GetRun(), 1, 1); + AliCDBMetaData *metadata= new AliCDBMetaData(); + + metadata->SetResponsible("alice-off@cern.ch"); + metadata->SetComment("Automatically produced GRP entry for Monte Carlo"); + + man->Put(grpObj,id,metadata); +} + +