#include <TSystem.h>
#include <TVirtualMC.h>
#include <TVirtualMCApplication.h>
+#include <TDatime.h>
#include "AliAlignObj.h"
#include "AliCDBEntry.h"
#include "AliCDBManager.h"
+#include "AliGRPManager.h"
#include "AliCDBStorage.h"
#include "AliCTPRawData.h"
#include "AliCentralTrigger.h"
fLoadAlObjsListOfDets("ALL"),
fMakeSDigits("ALL"),
fMakeDigits("ALL"),
- fMakeTrigger(""),
+ fTriggerConfig(""),
fMakeDigitsFromHits(""),
fWriteRawData(""),
fRawDataFileName(""),
fSetRunNumberFromDataCalled(kFALSE),
fEmbeddingFlag(kFALSE),
fLego(NULL),
+ fKey(0),
+ fUseVertexFromCDB(0),
+ fUseMagFieldFromGRP(0),
+ fGRPWriteLocation(Form("local://%s", gSystem->pwd())),
fQADetectors("ALL"),
fQATasks("ALL"),
fRunQA(kTRUE),
fEventSpecie(AliRecoParam::kDefault),
fWriteQAExpertData(kTRUE),
+ fGeometryFile(),
fRunHLT("default"),
+ fpHLT(NULL),
fWriteGRPEntry(kTRUE)
{
// create simulation object with default parameters
SetGAliceFile("galice.root");
// for QA
- AliQAManager * qam = AliQAManager::QAManager("sim") ;
+ AliQAManager * qam = AliQAManager::QAManager(AliQAv1::kSIMMODE) ;
qam->SetActiveDetectors(fQADetectors) ;
fQATasks = Form("%d %d %d", AliQAv1::kHITS, AliQAv1::kSDIGITS, AliQAv1::kDIGITS) ;
qam->SetTasks(fQATasks) ;
if (fInitCDBCalled) return;
fInitCDBCalled = kTRUE;
- AliQAManager * qam = AliQAManager::QAManager("sim") ;
+ AliQAManager * qam = AliQAManager::QAManager(AliQAv1::kSIMMODE) ;
qam->SetActiveDetectors(fQADetectors) ;
fQATasks = Form("%d %d %d", AliQAv1::kHITS, AliQAv1::kSDIGITS, AliQAv1::kDIGITS) ;
qam->SetTasks(fQATasks) ;
// 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);
+ ULong_t key = AliCDBManager::Instance()->SetLock(1);
+ if (key) fKey = key;
}
//_____________________________________________________________________________
}
// Export ideal geometry
- if(!gAlice->IsRootGeometry()) AliGeomManager::GetGeometry()->Export("geometry.root");
+ if(!IsGeometryFromFile()) AliGeomManager::GetGeometry()->Export("geometry.root");
// Load alignment data from CDB and apply to geometry through AliGeomManager
if(fLoadAlignFromCDB){
// run the generation, simulation and digitization
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
AliSysInfo::AddStamp("Start_Run");
// Load run number and seed from environmental vars
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 (!AliGeomManager::GetGeometry()) {
// Initialize the geometry manager
AliGeomManager::LoadGeometry("geometry.root");
- AliSysInfo::AddStamp("GetGeometry");
-
-
+ AliSysInfo::AddStamp("GetGeometry");
// // Check that the consistency of symbolic names for the activated subdetectors
// // in the geometry loaded by AliGeomManager
// AliRunLoader* runLoader = LoadRun("READ");
// digits -> trigger
- if (!RunTrigger(fMakeTrigger,fMakeDigits)) {
+ if (!RunTrigger(fTriggerConfig,fMakeDigits)) {
if (fStopOnError) return kFALSE;
}
// run the generation and simulation
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
// initialize CDB storage and run number from external environment
// (either CDB manager or AliSimulation setters)
AliPDG::AddParticlesToPdgDataBase();
- gAlice->GetMCApp()->Init();
-
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();
{
// run the trigger
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
// initialize CDB storage from external environment
// (either CDB manager or AliSimulation setters),
TString trconfiguration = config;
if (trconfiguration.IsNull()) {
- if (strcmp(gAlice->GetTriggerDescriptor(),"")) {
- trconfiguration = gAlice->GetTriggerDescriptor();
+ if(!fTriggerConfig.IsNull()) {
+ trconfiguration = fTriggerConfig;
}
else
AliWarning("No trigger descriptor is specified. Loading the one that is in the CDB.");
{
// run the generation and simulation
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
// initialize CDB storage and run number from external environment
// (either CDB manager or AliSimulation setters)
//
// Initialize ALICE Simulation run
//
-
gAlice->Announce();
+ //
+ // If requested set the mag. field from the GRP entry.
+ // After this the field is loccked and cannot be changed by Config.C
+ if (fUseMagFieldFromGRP) {
+ AliGRPManager grpM;
+ grpM.ReadGRPEntry();
+ grpM.SetMagField();
+ AliInfo("Field is locked now. It cannot be changed in Config.C");
+ }
+//
+// Execute Config.C
gROOT->LoadMacro(fConfigFileName.Data());
gInterpreter->ProcessLine(gAlice->GetConfigFunction());
AliSysInfo::AddStamp("RunSimulation_Config");
+//
+// If requested obtain the vertex position and vertex sigma_z from the CDB
+// This overwrites the settings from the Config.C
+ if (fUseVertexFromCDB) {
+ Double_t vtxPos[3] = {0., 0., 0.};
+ Double_t vtxSig[3] = {0., 0., 0.};
+ AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertex");
+ if (entry) {
+ AliESDVertex* vertex = dynamic_cast<AliESDVertex*> (entry->GetObject());
+ if (vertex) {
+ if(vertex->GetXRes()>2.8) { // > pipe radius --> it's a dummy object, don't use it
+ entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertexSPD");
+ if (entry) vertex = dynamic_cast<AliESDVertex*> (entry->GetObject());
+ }
+ }
+ if (vertex) {
+ vertex->GetXYZ(vtxPos);
+ vertex->GetSigmaXYZ(vtxSig);
+ AliInfo("Overwriting Config.C vertex settings !");
+ AliInfo(Form("Vertex position from OCDB entry: x = %13.3f, y = %13.3f, z = %13.3f (sigma = %13.3f)\n",
+ vtxPos[0], vtxPos[1], vtxPos[2], vtxSig[2]));
+
+ AliGenerator *gen = gAlice->GetMCApp()->Generator();
+ gen->SetOrigin(vtxPos[0], vtxPos[1], vtxPos[2]); // vertex position
+ gen->SetSigmaZ(vtxSig[2]);
+ }
+ }
+ }
+
if(AliCDBManager::Instance()->GetRun() >= 0) {
AliRunLoader::Instance()->SetRunNumber(AliCDBManager::Instance()->GetRun());
AliRunLoader::Instance()->SetNumberOfEventsPerRun(fNEvents);
}
AliRunLoader::Instance()->CdGAFile();
-
- AliPDG::AddParticlesToPdgDataBase();
+
- gAlice->GetMCApp()->Init();
- AliSysInfo::AddStamp("RunSimulation_InitMCApp");
+ 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");
AliSysInfo::AddStamp("RunSimulation_InitLoaders");
//___________________________________________________________________________________________
- // Get the trigger descriptor string
- // Either from AliSimulation or from
- // gAlice
- if (fMakeTrigger.IsNull()) {
- if (strcmp(gAlice->GetTriggerDescriptor(),""))
- fMakeTrigger = gAlice->GetTriggerDescriptor();
- }
- else
- gAlice->SetTriggerDescriptor(fMakeTrigger.Data());
AliSysInfo::AddStamp("RunSimulation_TriggerDescriptor");
// Set run number in CDBManager
loaderName += "Loader";
AliLoader* loader = runLoader->GetLoader(loaderName);
if (!loader) {
- AliError(Form("RunSimulation", "no loader for %s found\n"
- "Number of events per file not set for %s %s",
+ AliError(Form("RunSimulation no loader for %s found\n Number of events per file not set for %s %s",
detName, typeName, detName));
continue;
}
{
// run the digitization and produce summable digits
static Int_t eventNr=0;
- AliCodeTimerAuto("") ;
+ AliCodeTimerAuto("",0) ;
// initialize CDB storage, run number, set CDB lock
InitCDB();
{
// run the digitization and produce digits from sdigits
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
// initialize CDB storage, run number, set CDB lock
InitCDB();
{
// run the digitization and produce digits from hits
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
// initialize CDB storage, run number, set CDB lock
InitCDB();
// 'selrawdata' flag can be used to enable writing of detectors raw data
// accoring to the trigger cluster.
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
AliSysInfo::AddStamp("WriteRawData_Start");
TString detStr = detectors;
// run HLT simulation on simulated DDL raw files
// and produce HLT ddl raw files to be included in date/root file
- if (IsSelected("HLT", detStr) && !fRunHLT.IsNull()) {
+ // 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;
}
if (runLoader) for (Int_t iEvent = 0;
iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
char command[256];
- sprintf(command, "rm -r raw%d", iEvent);
+ snprintf(command, 256, "rm -r raw%d", iEvent);
gSystem->Exec(command);
}
+ delete runLoader;
}
if (rootOutput) {
{
// convert the digits to raw data DDL files
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
AliRunLoader* runLoader = LoadRun("READ");
if (!runLoader) return kFALSE;
runLoader->GetEvent(iEvent);
TString baseDir = gSystem->WorkingDirectory();
char dirName[256];
- sprintf(dirName, "raw%d", iEvent);
+ snprintf(dirName, 256, "raw%d", iEvent);
gSystem->MakeDirectory(dirName);
if (!gSystem->ChangeDirectory(dirName)) {
AliError(Form("couldn't change to directory %s", dirName));
// The second argument is not empty when the user decides to write
// the detectors raw data according to the trigger cluster.
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
char* path = gSystem->Which(gSystem->Getenv("PATH"), "dateStream");
if (!path) {
char command[256];
// 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",
+ snprintf(command, 256, "dateStream -c -s -D -o %s -# %d -C -run %d",
dateFileName, runLoader->GetNumberOfEvents(),runLoader->GetHeader()->GetRun());
FILE* pipe = gSystem->OpenPipe(command, "w");
ldc += AliDAQ::NumberOfLdcs(iDet) / AliDAQ::NumberOfDdls(iDet);
char rawFileName[256];
- sprintf(rawFileName, "raw%d/%s",
+ snprintf(rawFileName, 256, "raw%d/%s",
iEvent, AliDAQ::DdlFileName(iDet,iDDL));
// check existence and size of raw data file
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",
+ snprintf(command, 256, "dateStream -c -s -D -o %s -# %d -C -run %d",
selDateFileName,selEvents,runLoader->GetHeader()->GetRun());
FILE* pipe2 = gSystem->OpenPipe(command, "w");
ldc += AliDAQ::NumberOfLdcs(iDet) / AliDAQ::NumberOfDdls(iDet);
char rawFileName[256];
- sprintf(rawFileName, "raw%d/%s",
+ snprintf(rawFileName, 256, "raw%d/%s",
iEvent, AliDAQ::DdlFileName(iDet,iDDL));
// check existence and size of raw data file
}
//_____________________________________________________________________________
-Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* esdFileName)
+Bool_t AliSimulation::ConvertRaw2SDigits(const char* rawDirectory, const char* esdFileName, Int_t N)
{
//
// Steering routine to convert raw data in directory rawDirectory/ to fake SDigits.
AliPDG::AddParticlesToPdgDataBase();
- gAlice->GetMCApp()->Init();
-
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();
AliRunLoader* runLoader = AliRunLoader::Instance();
//
// Open esd file if available
- TFile* esdFile = TFile::Open(esdFileName);
+ TFile* esdFile = 0;
TTree* treeESD = 0;
- AliESDEvent* esd = new AliESDEvent();
- esdFile->GetObject("esdTree", treeESD);
- if (treeESD) esd->ReadFromTree(treeESD);
+ AliESDEvent* esd = 0;
+ if (esdFileName && (strlen(esdFileName)>0)) {
+ esdFile = TFile::Open(esdFileName);
+ if (esdFile) {
+ esd = new AliESDEvent();
+ esdFile->GetObject("esdTree", treeESD);
+ if (treeESD) esd->ReadFromTree(treeESD);
+ }
+ }
//
// Create the RawReader
} else if (!fileName.IsNull()) {
rawReader = new AliRawReaderDate(fileName);
}
+ if (!rawReader) return (kFALSE);
+
// if (!fEquipIdMap.IsNull() && fRawReader)
// fRawReader->LoadEquipmentIdsMap(fEquipIdMap);
//
//
// Get Header
AliHeader* header = runLoader->GetHeader();
- //
- TString detStr = fMakeSDigits;
// Event loop
Int_t nev = 0;
while(kTRUE) {
if (!(rawReader->NextEvent())) break;
+ runLoader->SetEventNumber(nev);
+ runLoader->GetHeader()->Reset(rawReader->GetRunNumber(),
+ nev, nev);
+ runLoader->GetEvent(nev);
+ AliInfo(Form("We are at event %d",nev));
//
// Detector loop
+ TString detStr = fMakeSDigits;
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()));
+ AliInfo(Form("Calling Raw2SDigits for %s", det->GetName()));
det->Raw2SDigits(rawReader);
rawReader->Reset();
}
} // detectors
-
+
//
// If ESD information available obtain reconstructed vertex and store in header.
mcHeader->SetPrimaryVertex(mcV);
header->Reset(0,nev);
header->SetGenEventHeader(mcHeader);
- printf("***** Saved vertex %f %f %f \n", position[0], position[1], position[2]);
+ AliInfo(Form("***** Saved vertex %f %f %f \n", position[0], position[1], position[2]));
}
- nev++;
//
// Finish the event
runLoader->TreeE()->Fill();
- runLoader->SetNextEvent();
+ AliInfo(Form("Finished event %d",nev));
+ nev++;
+ if (N>0&&nev>=N)
+ break;
} // events
delete rawReader;
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()
{
// 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();
- // 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)( 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));
- 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;
{
// run the QA on summable hits, digits or digits
- if(!gAlice) return kFALSE;
+ //if(!gAlice) return kFALSE;
AliQAManager::QAManager()->SetRunLoader(AliRunLoader::Instance()) ;
TString detectorsw("") ;
}
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 (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 {
grpObj->SetRunType("PHYSICS");
grpObj->SetTimeStart(0);
- grpObj->SetTimeEnd(9999);
+ 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();
+ Int_t a = 0;
+ Int_t z = 0;
+
if (gen) {
- grpObj->SetBeamEnergy(gen->GetEnergyCMS());
TString projectile;
- Int_t a,z;
gen->GetProjectile(projectile,a,z);
TString target;
gen->GetTarget(target,a,z);
beamType.ReplaceAll(" ","");
if (!beamType.CompareTo("-")) {
grpObj->SetBeamType("UNKNOWN");
+ grpObj->SetBeamEnergy(gen->GetEnergyCMS()/2);
}
else {
grpObj->SetBeamType(beamType);
+ if (z != 0) {
+ grpObj->SetBeamEnergy(gen->GetEnergyCMS()/2 * a / z);
+ } else {
+ grpObj->SetBeamEnergy(gen->GetEnergyCMS()/2 );
+ }
// Heavy ion run, the event specie is set to kHighMult
fEventSpecie = AliRecoParam::kHighMult;
if ((strcmp(beamType,"p-p") == 0) ||
}
}
// CTP
- if (!fMakeTrigger.IsNull() || strcmp(gAlice->GetTriggerDescriptor(),""))
+ if (!fTriggerConfig.IsNull())
detectorPattern |= (1 << AliDAQ::DetectorID("TRG"));
// HLT
grpObj->SetDetectorMask((Int_t)detectorPattern);
grpObj->SetLHCPeriod("LHC08c");
grpObj->SetLHCState("STABLE_BEAMS");
- grpObj->SetLHCLuminosity(0,(AliGRPObject::Stats)0);
- grpObj->SetBeamIntensity(0,(AliGRPObject::Stats)0);
-
+ //
AliMagF *field = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
Float_t solenoidField = field ? TMath::Abs(field->SolenoidField()) : 0;
- Float_t factor = field ? field->Factor() : 0;
- Float_t l3current = TMath::Abs(factor)*solenoidField*30000./5.;
- grpObj->SetL3Current(l3current,(AliGRPObject::Stats)0);
-
- if (factor > 0) {
- grpObj->SetL3Polarity(0);
- grpObj->SetDipolePolarity(0);
- }
- else {
- grpObj->SetL3Polarity(1);
- grpObj->SetDipolePolarity(1);
- }
-
- if (TMath::Abs(factor) != 0)
- grpObj->SetDipoleCurrent(6000,(AliGRPObject::Stats)0);
- else
- grpObj->SetDipoleCurrent(0,(AliGRPObject::Stats)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);
+ if (field) 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();
+
+ man->SetLock(0, fKey);
+
+ AliCDBStorage* sto = man->GetStorage(fGRPWriteLocation.Data());
+
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);
+ sto->Put(grpObj,id,metadata);
+ man->SetLock(1, fKey);
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff