]> git.uio.no Git - u/mrichter/AliRoot.git/blame - STEER/AliSimulation.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
New detector numbering scheme (common for DAQ/HLT/Offline). All the subdetectors...
[u/mrichter/AliRoot.git] / STEER / AliSimulation.cxx
CommitLineData
85a5290f 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/* $Id$ */
17
18///////////////////////////////////////////////////////////////////////////////
19// //
20// class for running generation, simulation and digitization //
21// //
22// Hits, sdigits and digits are created for all detectors by typing: //
23// //
24// AliSimulation sim; //
25// sim.Run(); //
26// //
27// The Run method returns kTRUE in case of successful execution. //
28// The number of events can be given as argument to the Run method or it //
29// can be set by //
30// //
31// sim.SetNumberOfEvents(n); //
32// //
95601147 33// The name of the configuration file can be passed as argument to the //
34// AliSimulation constructor or can be specified by //
85a5290f 35// //
36// sim.SetConfigFile("..."); //
37// //
38// The generation of particles and the simulation of detector hits can be //
39// switched on or off by //
40// //
41// sim.SetRunGeneration(kTRUE); // generation of primary particles //
42// sim.SetRunSimulation(kFALSE); // but no tracking //
43// //
44// For which detectors sdigits and digits will be created, can be steered //
45// by //
46// //
47// sim.SetMakeSDigits("ALL"); // make sdigits for all detectors //
48// sim.SetMakeDigits("ITS TPC"); // make digits only for ITS and TPC //
49// //
50// The argument is a (case sensitive) string with the names of the //
51// detectors separated by a space. An empty string ("") can be used to //
52// disable the creation of sdigits or digits. The special string "ALL" //
53// selects all available detectors. This is the default. //
54// //
55// The creation of digits from hits instead of from sdigits can be selected //
56// by //
57// //
58// sim.SetMakeDigitsFromHits("TRD"); //
59// //
60// The argument is again a string with the selected detectors. Be aware that //
61// this feature is not available for all detectors and that merging is not //
62// possible, when digits are created directly from hits. //
63// //
05526d44 64// Background events can be merged by calling //
85a5290f 65// //
66// sim.MergeWith("background/galice.root", 2); //
67// //
68// The first argument is the file name of the background galice file. The //
69// second argument is the number of signal events per background event. //
05526d44 70// By default this number is calculated from the number of available //
71// background events. MergeWith can be called several times to merge more //
72// than two event streams. It is assumed that the sdigits were already //
73// produced for the background events. //
85a5290f 74// //
0421c3d1 75// The output of raw data can be switched on by calling //
76// //
77// sim.SetWriteRawData("MUON"); // write raw data for MUON //
78// //
06c7e392 79// The default output format of the raw data are DDL files. They are //
80// converted to a DATE file, if a file name is given as second argument. //
81// For this conversion the program "dateStream" is required. If the file //
82// name has the extension ".root", the DATE file is converted to a root //
83// file. The program "alimdc" is used for this purpose. For the conversion //
84// to DATE and root format the two conversion programs have to be installed. //
85// Only the raw data in the final format is kept if the third argument is //
86// kTRUE. //
87// //
0421c3d1 88// The methods RunSimulation, RunSDigitization, RunDigitization, //
89// RunHitsDigitization and WriteRawData can be used to run only parts of //
06c7e392 90// the full simulation chain. The creation of raw data DDL files and their //
91// conversion to the DATE or root format can be run directly by calling //
92// the methods WriteRawFiles, ConvertRawFilesToDate and ConvertDateToRoot. //
95601147 93// //
18b43626 94// The default number of events per file, which is usually set in the //
95// config file, can be changed for individual detectors and data types //
96// by calling //
97// //
98// sim.SetEventsPerFile("PHOS", "Reconstructed Points", 3); //
99// //
100// The first argument is the detector, the second one the data type and the //
101// last one the number of events per file. Valid data types are "Hits", //
102// "Summable Digits", "Digits", "Reconstructed Points" and "Tracks". //
103// The number of events per file has to be set before the simulation of //
104// hits. Otherwise it has no effect. //
105// //
85a5290f 106///////////////////////////////////////////////////////////////////////////////
107
e8d02863 108#include <TGeoManager.h>
85a5290f 109#include <TObjString.h>
fd46e2d2 110#include <TStopwatch.h>
af7ba10c 111#include <TSystem.h>
bfbd5665 112#include <TFile.h>
85a5290f 113
e8d02863 114#include "AliCDBStorage.h"
115#include "AliCDBEntry.h"
116#include "AliCDBManager.h"
117#include "AliAlignObj.h"
118#include "AliCentralTrigger.h"
362c9d61 119#include "AliDAQ.h"
af7ba10c 120#include "AliDigitizer.h"
121#include "AliGenerator.h"
e8d02863 122#include "AliLog.h"
af7ba10c 123#include "AliModule.h"
124#include "AliRun.h"
125#include "AliRunDigitizer.h"
126#include "AliRunLoader.h"
127#include "AliSimulation.h"
128#include "AliVertexGenFile.h"
a5a091ce 129#include "AliCentralTrigger.h"
007155cc 130#include "AliCTPRawData.h"
85a5290f 131
132ClassImp(AliSimulation)
133
134
135//_____________________________________________________________________________
024cf675 136AliSimulation::AliSimulation(const char* configFileName, const char* cdbUri,
95601147 137 const char* name, const char* title) :
138 TNamed(name, title),
139
140 fRunGeneration(kTRUE),
141 fRunSimulation(kTRUE),
024cf675 142 fLoadAlignFromCDB(kTRUE),
143 fLoadAlignData("ALL"),
95601147 144 fMakeSDigits("ALL"),
145 fMakeDigits("ALL"),
a5a091ce 146 fMakeTrigger(""),
95601147 147 fMakeDigitsFromHits(""),
0421c3d1 148 fWriteRawData(""),
06c7e392 149 fRawDataFileName(""),
150 fDeleteIntermediateFiles(kFALSE),
95601147 151 fStopOnError(kFALSE),
152
153 fNEvents(1),
154 fConfigFileName(configFileName),
155 fGAliceFileName("galice.root"),
18b43626 156 fEventsPerFile(),
95601147 157 fBkgrdFileNames(NULL),
995ad051 158 fAlignObjArray(NULL),
04bae0a0 159 fUseBkgrdVertex(kTRUE),
024cf675 160 fRegionOfInterest(kFALSE),
ec92bee0 161 fCDBUri(cdbUri),
162 fSpecCDBUri()
85a5290f 163{
164// create simulation object with default parameters
165
0421c3d1 166 SetGAliceFile("galice.root");
85a5290f 167}
168
169//_____________________________________________________________________________
170AliSimulation::AliSimulation(const AliSimulation& sim) :
95601147 171 TNamed(sim),
172
173 fRunGeneration(sim.fRunGeneration),
174 fRunSimulation(sim.fRunSimulation),
024cf675 175 fLoadAlignFromCDB(sim.fLoadAlignFromCDB),
176 fLoadAlignData(sim.fLoadAlignData),
95601147 177 fMakeSDigits(sim.fMakeSDigits),
178 fMakeDigits(sim.fMakeDigits),
a5a091ce 179 fMakeTrigger(sim.fMakeTrigger),
95601147 180 fMakeDigitsFromHits(sim.fMakeDigitsFromHits),
0421c3d1 181 fWriteRawData(sim.fWriteRawData),
06c7e392 182 fRawDataFileName(""),
183 fDeleteIntermediateFiles(kFALSE),
95601147 184 fStopOnError(sim.fStopOnError),
185
186 fNEvents(sim.fNEvents),
187 fConfigFileName(sim.fConfigFileName),
188 fGAliceFileName(sim.fGAliceFileName),
18b43626 189 fEventsPerFile(),
95601147 190 fBkgrdFileNames(NULL),
995ad051 191 fAlignObjArray(NULL),
04bae0a0 192 fUseBkgrdVertex(sim.fUseBkgrdVertex),
024cf675 193 fRegionOfInterest(sim.fRegionOfInterest),
ec92bee0 194 fCDBUri(sim.fCDBUri),
195 fSpecCDBUri()
85a5290f 196{
197// copy constructor
198
18b43626 199 for (Int_t i = 0; i < sim.fEventsPerFile.GetEntriesFast(); i++) {
200 if (!sim.fEventsPerFile[i]) continue;
201 fEventsPerFile.Add(sim.fEventsPerFile[i]->Clone());
202 }
203
85a5290f 204 fBkgrdFileNames = new TObjArray;
205 for (Int_t i = 0; i < sim.fBkgrdFileNames->GetEntriesFast(); i++) {
206 if (!sim.fBkgrdFileNames->At(i)) continue;
207 fBkgrdFileNames->Add(sim.fBkgrdFileNames->At(i)->Clone());
208 }
ec92bee0 209
210 for (Int_t i = 0; i < sim.fSpecCDBUri.GetEntriesFast(); i++) {
211 if (sim.fSpecCDBUri[i]) fSpecCDBUri.Add(sim.fSpecCDBUri[i]->Clone());
212 }
213
85a5290f 214}
215
216//_____________________________________________________________________________
217AliSimulation& AliSimulation::operator = (const AliSimulation& sim)
218{
219// assignment operator
220
221 this->~AliSimulation();
222 new(this) AliSimulation(sim);
223 return *this;
224}
225
226//_____________________________________________________________________________
227AliSimulation::~AliSimulation()
228{
229// clean up
230
18b43626 231 fEventsPerFile.Delete();
024cf675 232// if(fAlignObjArray) fAlignObjArray->Delete(); // fAlignObjArray->RemoveAll() ???
233// delete fAlignObjArray; fAlignObjArray=0;
18b43626 234
95601147 235 if (fBkgrdFileNames) {
236 fBkgrdFileNames->Delete();
237 delete fBkgrdFileNames;
238 }
ec92bee0 239
240 fSpecCDBUri.Delete();
85a5290f 241}
242
243
244//_____________________________________________________________________________
245void AliSimulation::SetNumberOfEvents(Int_t nEvents)
246{
247// set the number of events for one run
248
249 fNEvents = nEvents;
250}
251
024cf675 252//_____________________________________________________________________________
ec92bee0 253void AliSimulation::InitCDBStorage()
024cf675 254{
255// activate a default CDB storage
256// First check if we have any CDB storage set, because it is used
257// to retrieve the calibration and alignment constants
258
259 AliCDBManager* man = AliCDBManager::Instance();
ec92bee0 260 if (man->IsDefaultStorageSet())
024cf675 261 {
ec92bee0 262 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
263 AliWarning("Default CDB storage has been already set !");
264 AliWarning(Form("Ignoring the default storage declared in AliReconstruction: %s",fCDBUri.Data()));
265 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
266 fCDBUri = "";
267 }
268 else {
269 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
270 AliWarning(Form("Default CDB storage is set to: %s",fCDBUri.Data()));
271 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
272 man->SetDefaultStorage(fCDBUri);
273 }
274
275 // Now activate the detector specific CDB storage locations
276 for (Int_t i = 0; i < fSpecCDBUri.GetEntriesFast(); i++) {
277 TObject* obj = fSpecCDBUri[i];
278 if (!obj) continue;
279 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
280 AliWarning(Form("Specific CDB storage for %s is set to: %s",obj->GetName(),obj->GetTitle()));
281 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
282 man->SetSpecificStorage(obj->GetName(),obj->GetTitle());
283 }
024cf675 284
285}
286
287//_____________________________________________________________________________
288void AliSimulation::SetDefaultStorage(const char* uri) {
ec92bee0 289// Store the desired default CDB storage location
290// Activate it later within the Run() method
024cf675 291
ec92bee0 292 fCDBUri = uri;
024cf675 293
294}
295
296//_____________________________________________________________________________
297void AliSimulation::SetSpecificStorage(const char* detName, const char* uri) {
ec92bee0 298// Store a detector-specific CDB storage location
299// Activate it later within the Run() method
024cf675 300
ec92bee0 301 TObject* obj = fSpecCDBUri.FindObject(detName);
302 if (obj) fSpecCDBUri.Remove(obj);
303 fSpecCDBUri.Add(new TNamed(detName, uri));
024cf675 304
305}
306
85a5290f 307//_____________________________________________________________________________
308void AliSimulation::SetConfigFile(const char* fileName)
309{
310// set the name of the config file
311
312 fConfigFileName = fileName;
313}
314
0421c3d1 315//_____________________________________________________________________________
316void AliSimulation::SetGAliceFile(const char* fileName)
317{
318// set the name of the galice file
319// the path is converted to an absolute one if it is relative
320
321 fGAliceFileName = fileName;
322 if (!gSystem->IsAbsoluteFileName(fGAliceFileName)) {
323 char* absFileName = gSystem->ConcatFileName(gSystem->WorkingDirectory(),
324 fGAliceFileName);
325 fGAliceFileName = absFileName;
326 delete[] absFileName;
327 }
ac1671c0 328
329 AliDebug(2, Form("galice file name set to %s", fileName));
0421c3d1 330}
331
18b43626 332//_____________________________________________________________________________
333void AliSimulation::SetEventsPerFile(const char* detector, const char* type,
334 Int_t nEvents)
335{
336// set the number of events per file for the given detector and data type
337// ("Hits", "Summable Digits", "Digits", "Reconstructed Points" or "Tracks")
338
339 TNamed* obj = new TNamed(detector, type);
340 obj->SetUniqueID(nEvents);
341 fEventsPerFile.Add(obj);
342}
343
6bae477a 344//_____________________________________________________________________________
345Bool_t AliSimulation::ApplyAlignObjsToGeom(TObjArray* alObjArray)
346{
347 // Read collection of alignment objects (AliAlignObj derived) saved
348 // in the TClonesArray ClArrayName and apply them to the geometry
349 // manager singleton.
350 //
351 alObjArray->Sort();
352 Int_t nvols = alObjArray->GetEntriesFast();
353
b5874ce3 354 Bool_t flag = kTRUE;
355
6bae477a 356 for(Int_t j=0; j<nvols; j++)
357 {
358 AliAlignObj* alobj = (AliAlignObj*) alObjArray->UncheckedAt(j);
b5874ce3 359 if (alobj->ApplyToGeometry() == kFALSE) flag = kFALSE;
6bae477a 360 }
361
362 if (AliDebugLevelClass() >= 1) {
e0625db4 363 gGeoManager->GetTopNode()->CheckOverlaps(20);
6bae477a 364 TObjArray* ovexlist = gGeoManager->GetListOfOverlaps();
365 if(ovexlist->GetEntriesFast()){
366 AliError("The application of alignment objects to the geometry caused huge overlaps/extrusions!");
367 }
368 }
369
b5874ce3 370 return flag;
6bae477a 371
372}
373
b9bc0f82 374//_____________________________________________________________________________
178a256c 375Bool_t AliSimulation::ApplyAlignObjsToGeom(const char* fileName, const char* clArrayName)
b9bc0f82 376{
377 // read collection of alignment objects (AliAlignObj derived) saved
378 // in the TClonesArray ClArrayName in the file fileName and apply
379 // them to the TGeo geometry passed as argument
380 //
381
382 TFile* inFile = TFile::Open(fileName,"READ");
383 if (!inFile || !inFile->IsOpen()) {
384 AliErrorClass(Form("Could not open file %s !",fileName));
385 return kFALSE;
386 }
387
995ad051 388 TClonesArray* alObjArray = ((TClonesArray*) inFile->Get(clArrayName));
8c30cb69 389 inFile->Close();
995ad051 390 if (!alObjArray) {
391 AliErrorClass(Form("Could not get array (%s) from file (%s) !",clArrayName,fileName));
392 return kFALSE;
393 }
8c30cb69 394
6bae477a 395 return ApplyAlignObjsToGeom(alObjArray);
b9bc0f82 396
397}
398
399//_____________________________________________________________________________
178a256c 400Bool_t AliSimulation::ApplyAlignObjsToGeom(AliCDBParam* param, AliCDBId& Id)
b9bc0f82 401{
402 // read collection of alignment objects (AliAlignObj derived) saved
403 // in the TClonesArray ClArrayName in the AliCDBEntry identified by
404 // param (to get the AliCDBStorage) and Id; apply the alignment objects
405 // to the TGeo geometry passed as argument
406 //
407
408 AliCDBStorage* storage = AliCDBManager::Instance()->GetStorage(param);
409 AliCDBEntry* entry = storage->Get(Id);
410 TClonesArray* AlObjArray = ((TClonesArray*) entry->GetObject());
411
6bae477a 412 return ApplyAlignObjsToGeom(AlObjArray);
b9bc0f82 413
414}
415
416//_____________________________________________________________________________
178a256c 417Bool_t AliSimulation::ApplyAlignObjsToGeom(const char* uri, const char* path, Int_t runnum, Int_t version, Int_t sversion)
b9bc0f82 418{
419 // read collection of alignment objects (AliAlignObj derived) saved
420 // in the TClonesArray ClArrayName in the AliCDBEntry identified by
421 // param (to get the AliCDBStorage) and Id; apply the alignment objects
422 // to the TGeo geometry passed as argument
423 //
424
425 AliCDBParam* param = AliCDBManager::Instance()->CreateParameter(uri);
426 AliCDBId id(path, runnum, runnum, version, sversion);
b9bc0f82 427
178a256c 428 return ApplyAlignObjsToGeom(param, id);
b9bc0f82 429
430}
431
024cf675 432//_____________________________________________________________________________
433Bool_t AliSimulation::ApplyAlignObjsToGeom(const char* detName, Int_t runnum, Int_t version, Int_t sversion)
434{
435 // read collection of alignment objects (AliAlignObj derived) saved
436 // in the TClonesArray ClArrayName in the AliCDBEntry identified by
437 // param (to get the AliCDBStorage) and Id; apply the alignment objects
438 // to the TGeo geometry passed as argument
439 //
440
024cf675 441 AliCDBPath path(detName,"Align","Data");
442 AliCDBEntry* entry = AliCDBManager::Instance()->Get(path.GetPath(),runnum,version,sversion);
443
444 if(!entry) return kFALSE;
445 TClonesArray* AlObjArray = ((TClonesArray*) entry->GetObject());
446
6bae477a 447 return ApplyAlignObjsToGeom(AlObjArray);
024cf675 448}
449
450//_____________________________________________________________________________
451Bool_t AliSimulation::SetAlignObjArraySingleDet(const char* detName)
452{
453 // Fills array of single detector's alignable objects from CDB
454
455 AliDebug(2, Form("Loading alignment data for detector: %s",detName));
456
457 AliCDBEntry *entry;
458
459 AliCDBPath path(detName,"Align","Data");
460
461 entry=AliCDBManager::Instance()->Get(path.GetPath());
462 if(!entry){
463 AliDebug(2,Form("Couldn't load alignment data for detector %s",detName));
464 return kFALSE;
465 }
466 entry->SetOwner(1);
467 TClonesArray *alignArray = (TClonesArray*) entry->GetObject();
468 alignArray->SetOwner(0);
469 AliDebug(2,Form("Found %d alignment objects for %s",
470 alignArray->GetEntries(),detName));
471
472 AliAlignObj *alignObj=0;
473 TIter iter(alignArray);
474
475 // loop over align objects in detector
476 while( ( alignObj=(AliAlignObj *) iter.Next() ) ){
477 fAlignObjArray->Add(alignObj);
478 }
479 // delete entry --- Don't delete, it is cached!
480
481 AliDebug(2, Form("fAlignObjArray entries: %d",fAlignObjArray->GetEntries() ));
482 return kTRUE;
483
484}
485
6bae477a 486//_____________________________________________________________________________
487Bool_t AliSimulation::MisalignGeometry(AliRunLoader *runLoader)
488{
489 // Read the alignment objects from CDB.
490 // Each detector is supposed to have the
491 // alignment objects in DET/Align/Data CDB path.
492 // All the detector objects are then collected,
493 // sorted by geometry level (starting from ALIC) and
494 // then applied to the TGeo geometry.
495 // Finally an overlaps check is performed.
496
497 Bool_t delRunLoader = kFALSE;
498 if (!runLoader) {
499 runLoader = LoadRun("READ");
500 if (!runLoader) return kFALSE;
501 delRunLoader = kTRUE;
502 }
503
504 // Load alignment data from CDB and fill fAlignObjArray
505 if(fLoadAlignFromCDB){
506 if(!fAlignObjArray) fAlignObjArray = new TObjArray();
507
508 //fAlignObjArray->RemoveAll();
509 fAlignObjArray->Clear();
510 fAlignObjArray->SetOwner(0);
511
512 TString detStr = fLoadAlignData;
513 TString dataNotLoaded="";
514 TString dataLoaded="";
515
516 TObjArray* detArray = runLoader->GetAliRun()->Detectors();
517 for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
518 AliModule* det = (AliModule*) detArray->At(iDet);
519 if (!det || !det->IsActive()) continue;
520 if (IsSelected(det->GetName(), detStr)) {
521 if(!SetAlignObjArraySingleDet(det->GetName())){
522 dataNotLoaded += det->GetName();
523 dataNotLoaded += " ";
524 } else {
525 dataLoaded += det->GetName();
526 dataLoaded += " ";
527 }
528 }
529 } // end loop over detectors
530
531 if ((detStr.CompareTo("ALL") == 0)) detStr = "";
532 dataNotLoaded += detStr;
533 AliInfo(Form("Alignment data loaded for: %s",
534 dataLoaded.Data()));
535 AliInfo(Form("Didn't/couldn't load alignment data for: %s",
536 dataNotLoaded.Data()));
537 } // fLoadAlignFromCDB flag
538
539 // Check if the array with alignment objects was
540 // provided by the user. If yes, apply the objects
541 // to the present TGeo geometry
542 if (fAlignObjArray) {
543 if (gGeoManager && gGeoManager->IsClosed()) {
544 if (ApplyAlignObjsToGeom(fAlignObjArray) == kFALSE) {
b5874ce3 545 AliError("The misalignment of one or more volumes failed!"
546 "Compare the list of simulated detectors and the list of detector alignment data!");
547 if (delRunLoader) delete runLoader;
6bae477a 548 return kFALSE;
549 }
550 }
551 else {
552 AliError("Can't apply the misalignment! gGeoManager doesn't exist or it is still opened!");
b5874ce3 553 if (delRunLoader) delete runLoader;
6bae477a 554 return kFALSE;
555 }
556 }
557
558 if (delRunLoader) delete runLoader;
559
560 return kTRUE;
561}
562
563
564//_____________________________________________________________________________
565Bool_t AliSimulation::SetRunNumber()
566{
567 // Set the CDB manager run number
568 // The run number is retrieved from gAlice
569
570 if(AliCDBManager::Instance()->GetRun() < 0) {
571 AliRunLoader* runLoader = LoadRun("READ");
572 if (!runLoader) return kFALSE;
573 else {
574 AliCDBManager::Instance()->SetRun(runLoader->GetAliRun()->GetRunNumber());
575 AliInfo(Form("Run number: %d",AliCDBManager::Instance()->GetRun()));
576 delete runLoader;
577 }
578 }
579 return kTRUE;
580}
581
85a5290f 582//_____________________________________________________________________________
583void AliSimulation::MergeWith(const char* fileName, Int_t nSignalPerBkgrd)
584{
585// add a file with background events for merging
586
587 TObjString* fileNameStr = new TObjString(fileName);
588 fileNameStr->SetUniqueID(nSignalPerBkgrd);
95601147 589 if (!fBkgrdFileNames) fBkgrdFileNames = new TObjArray;
85a5290f 590 fBkgrdFileNames->Add(fileNameStr);
591}
592
593
594//_____________________________________________________________________________
595Bool_t AliSimulation::Run(Int_t nEvents)
596{
597// run the generation, simulation and digitization
598
ec92bee0 599 InitCDBStorage();
b0314964 600
85a5290f 601 if (nEvents > 0) fNEvents = nEvents;
602
603 // generation and simulation -> hits
604 if (fRunGeneration) {
85a5290f 605 if (!RunSimulation()) if (fStopOnError) return kFALSE;
606 }
607
6bae477a 608 // Set run number in CDBManager (if it is not already set in RunSimulation)
609 if (!SetRunNumber()) if (fStopOnError) return kFALSE;
610
611 // Load and misalign the geometry
612 if (!gGeoManager) {
613 TGeoManager::Import("geometry.root");
614 if (!gGeoManager) if (fStopOnError) return kFALSE;
615 if (!MisalignGeometry()) if (fStopOnError) return kFALSE;
616 }
617
85a5290f 618 // hits -> summable digits
619 if (!fMakeSDigits.IsNull()) {
620 if (!RunSDigitization(fMakeSDigits)) if (fStopOnError) return kFALSE;
621 }
622
623 // summable digits -> digits
624 if (!fMakeDigits.IsNull()) {
625 if (!RunDigitization(fMakeDigits, fMakeDigitsFromHits)) {
626 if (fStopOnError) return kFALSE;
627 }
628 }
629
630 // hits -> digits
631 if (!fMakeDigitsFromHits.IsNull()) {
95601147 632 if (fBkgrdFileNames && (fBkgrdFileNames->GetEntriesFast() > 0)) {
ac1671c0 633 AliWarning(Form("Merging and direct creation of digits from hits "
634 "was selected for some detectors. "
635 "No merging will be done for the following detectors: %s",
636 fMakeDigitsFromHits.Data()));
85a5290f 637 }
638 if (!RunHitsDigitization(fMakeDigitsFromHits)) {
639 if (fStopOnError) return kFALSE;
640 }
641 }
642
a5a091ce 643 // digits -> trigger
087426de 644 if (!RunTrigger(fMakeTrigger)) {
645 if (fStopOnError) return kFALSE;
a5a091ce 646 }
647
0421c3d1 648 // digits -> raw data
649 if (!fWriteRawData.IsNull()) {
06c7e392 650 if (!WriteRawData(fWriteRawData, fRawDataFileName,
651 fDeleteIntermediateFiles)) {
0421c3d1 652 if (fStopOnError) return kFALSE;
653 }
654 }
655
85a5290f 656 return kTRUE;
657}
658
a5a091ce 659//_____________________________________________________________________________
660Bool_t AliSimulation::RunTrigger(const char* descriptors)
661{
662 // run the trigger
663
bacbe0fd 664 TStopwatch stopwatch;
665 stopwatch.Start();
a5a091ce 666
bacbe0fd 667 AliRunLoader* runLoader = LoadRun("READ");
668 if (!runLoader) return kFALSE;
669 TString des = descriptors;
670
087426de 671 if (des.IsNull()) {
672 if (gAlice->GetTriggerDescriptor() != "") {
673 des = gAlice->GetTriggerDescriptor();
674 }
675 else {
676 AliWarning("No trigger descriptor is specified. Skipping the trigger simulation...");
677 return kTRUE;
678 }
679 }
680
bacbe0fd 681 runLoader->MakeTree( "CT" );
682 AliCentralTrigger* aCTP = runLoader->GetTrigger();
a5a091ce 683 // Load Descriptors
bacbe0fd 684 aCTP->LoadDescriptor( des );
a5a091ce 685
686 // digits -> trigger
bacbe0fd 687 if( !aCTP->RunTrigger( runLoader ) ) {
a5a091ce 688 if (fStopOnError) {
bacbe0fd 689 // delete aCTP;
a5a091ce 690 return kFALSE;
691 }
bacbe0fd 692 }
a5a091ce 693
bacbe0fd 694 AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
695 stopwatch.RealTime(),stopwatch.CpuTime()));
a5a091ce 696
bacbe0fd 697 delete runLoader;
a5a091ce 698
bacbe0fd 699 return kTRUE;
a5a091ce 700}
701
007155cc 702//_____________________________________________________________________________
703Bool_t AliSimulation::WriteTriggerRawData()
704{
705 // Writes the CTP (trigger) DDL raw data
706 // Details of the format are given in the
707 // trigger TDR - pages 134 and 135.
708 AliCTPRawData writer;
709 writer.RawData();
a5a091ce 710
007155cc 711 return kTRUE;
712}
a5a091ce 713
85a5290f 714//_____________________________________________________________________________
95601147 715Bool_t AliSimulation::RunSimulation(Int_t nEvents)
85a5290f 716{
717// run the generation and simulation
718
4df28b43 719 TStopwatch stopwatch;
720 stopwatch.Start();
721
95601147 722 if (!gAlice) {
ac1671c0 723 AliError("no gAlice object. Restart aliroot and try again.");
95601147 724 return kFALSE;
725 }
726 if (gAlice->Modules()->GetEntries() > 0) {
ac1671c0 727 AliError("gAlice was already run. Restart aliroot and try again.");
95601147 728 return kFALSE;
729 }
730
ac1671c0 731 AliInfo(Form("initializing gAlice with config file %s",
732 fConfigFileName.Data()));
815c2b38 733 StdoutToAliInfo(StderrToAliError(
734 gAlice->Init(fConfigFileName.Data());
735 ););
995ad051 736
087426de 737 // Get the trigger descriptor string
738 // Either from AliSimulation or from
739 // gAlice
740 if (fMakeTrigger.IsNull()) {
741 if (gAlice->GetTriggerDescriptor() != "")
742 fMakeTrigger = gAlice->GetTriggerDescriptor();
743 }
744 else
745 gAlice->SetTriggerDescriptor(fMakeTrigger.Data());
746
6bae477a 747 // Set run number in CDBManager
024cf675 748 AliCDBManager::Instance()->SetRun(gAlice->GetRunNumber());
749 AliInfo(Form("Run number: %d",AliCDBManager::Instance()->GetRun()));
750
751 AliRunLoader* runLoader = gAlice->GetRunLoader();
752 if (!runLoader) {
753 AliError(Form("gAlice has no run loader object. "
754 "Check your config file: %s", fConfigFileName.Data()));
755 return kFALSE;
756 }
757 SetGAliceFile(runLoader->GetFileName());
758
bea6ea16 759 // Export ideal geometry
760 if (gGeoManager) gGeoManager->Export("geometry.root");
024cf675 761
6bae477a 762 // Misalign geometry
e198134a 763// if (!MisalignGeometry(runLoader)) {
764// delete runLoader;
765// return kFALSE;
766// }
767 MisalignGeometry(runLoader);
995ad051 768
99673dbf 769 // Temporary fix by A.Gheata
770 // Could be removed with the next Root version (>5.11)
771 if (gGeoManager) {
772 TIter next(gGeoManager->GetListOfVolumes());
773 TGeoVolume *vol;
774 while ((vol = (TGeoVolume *)next())) {
775 if (vol->GetVoxels()) {
776 if (vol->GetVoxels()->NeedRebuild()) {
777 vol->GetVoxels()->Voxelize();
778 vol->FindOverlaps();
779 }
780 }
781 }
782 }
783
6bae477a 784 // Export (mis)aligned geometry
785 if (gGeoManager) gGeoManager->Export("misaligned_geometry.root");
786
024cf675 787// AliRunLoader* runLoader = gAlice->GetRunLoader();
788// if (!runLoader) {
789// AliError(Form("gAlice has no run loader object. "
790// "Check your config file: %s", fConfigFileName.Data()));
791// return kFALSE;
792// }
793// SetGAliceFile(runLoader->GetFileName());
85a5290f 794
04bae0a0 795 if (!gAlice->Generator()) {
ac1671c0 796 AliError(Form("gAlice has no generator object. "
797 "Check your config file: %s", fConfigFileName.Data()));
04bae0a0 798 return kFALSE;
799 }
05526d44 800 if (nEvents <= 0) nEvents = fNEvents;
04bae0a0 801
802 // get vertex from background file in case of merging
803 if (fUseBkgrdVertex &&
804 fBkgrdFileNames && (fBkgrdFileNames->GetEntriesFast() > 0)) {
05526d44 805 Int_t signalPerBkgrd = GetNSignalPerBkgrd(nEvents);
04bae0a0 806 const char* fileName = ((TObjString*)
807 (fBkgrdFileNames->At(0)))->GetName();
ac1671c0 808 AliInfo(Form("The vertex will be taken from the background "
809 "file %s with nSignalPerBackground = %d",
810 fileName, signalPerBkgrd));
04bae0a0 811 AliVertexGenFile* vtxGen = new AliVertexGenFile(fileName, signalPerBkgrd);
812 gAlice->Generator()->SetVertexGenerator(vtxGen);
813 }
814
85a5290f 815 if (!fRunSimulation) {
85a5290f 816 gAlice->Generator()->SetTrackingFlag(0);
817 }
818
18b43626 819 // set the number of events per file for given detectors and data types
820 for (Int_t i = 0; i < fEventsPerFile.GetEntriesFast(); i++) {
821 if (!fEventsPerFile[i]) continue;
822 const char* detName = fEventsPerFile[i]->GetName();
823 const char* typeName = fEventsPerFile[i]->GetTitle();
824 TString loaderName(detName);
825 loaderName += "Loader";
826 AliLoader* loader = runLoader->GetLoader(loaderName);
827 if (!loader) {
ac1671c0 828 AliError(Form("RunSimulation", "no loader for %s found\n"
829 "Number of events per file not set for %s %s",
830 detName, typeName, detName));
18b43626 831 continue;
832 }
833 AliDataLoader* dataLoader =
834 loader->GetDataLoader(typeName);
835 if (!dataLoader) {
ac1671c0 836 AliError(Form("no data loader for %s found\n"
837 "Number of events per file not set for %s %s",
838 typeName, detName, typeName));
18b43626 839 continue;
840 }
841 dataLoader->SetNumberOfEventsPerFile(fEventsPerFile[i]->GetUniqueID());
ac1671c0 842 AliDebug(1, Form("number of events per file set to %d for %s %s",
843 fEventsPerFile[i]->GetUniqueID(), detName, typeName));
18b43626 844 }
845
ac1671c0 846 AliInfo("running gAlice");
815c2b38 847 StdoutToAliInfo(StderrToAliError(
848 gAlice->Run(nEvents);
849 ););
95601147 850
851 delete runLoader;
85a5290f 852
5f8272e1 853 AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
854 stopwatch.RealTime(),stopwatch.CpuTime()));
4df28b43 855
85a5290f 856 return kTRUE;
857}
858
859//_____________________________________________________________________________
95601147 860Bool_t AliSimulation::RunSDigitization(const char* detectors)
85a5290f 861{
862// run the digitization and produce summable digits
863
4df28b43 864 TStopwatch stopwatch;
865 stopwatch.Start();
866
95601147 867 AliRunLoader* runLoader = LoadRun();
868 if (!runLoader) return kFALSE;
869
85a5290f 870 TString detStr = detectors;
95601147 871 TObjArray* detArray = runLoader->GetAliRun()->Detectors();
85a5290f 872 for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
873 AliModule* det = (AliModule*) detArray->At(iDet);
874 if (!det || !det->IsActive()) continue;
875 if (IsSelected(det->GetName(), detStr)) {
ac1671c0 876 AliInfo(Form("creating summable digits for %s", det->GetName()));
de76655b 877 TStopwatch stopwatchDet;
878 stopwatchDet.Start();
85a5290f 879 det->Hits2SDigits();
5f8272e1 880 AliInfo(Form("Execution time for %s: R:%.2fs C:%.2fs",
881 det->GetName(),stopwatchDet.RealTime(),stopwatchDet.CpuTime()));
85a5290f 882 }
883 }
884
885 if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
ac1671c0 886 AliError(Form("the following detectors were not found: %s",
887 detStr.Data()));
85a5290f 888 if (fStopOnError) return kFALSE;
889 }
890
95601147 891 delete runLoader;
892
5f8272e1 893 AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
894 stopwatch.RealTime(),stopwatch.CpuTime()));
4df28b43 895
85a5290f 896 return kTRUE;
897}
898
899
900//_____________________________________________________________________________
95601147 901Bool_t AliSimulation::RunDigitization(const char* detectors,
902 const char* excludeDetectors)
85a5290f 903{
904// run the digitization and produce digits from sdigits
905
4df28b43 906 TStopwatch stopwatch;
907 stopwatch.Start();
908
95601147 909 while (AliRunLoader::GetRunLoader()) delete AliRunLoader::GetRunLoader();
910 if (gAlice) delete gAlice;
911 gAlice = NULL;
912
913 Int_t nStreams = 1;
05526d44 914 if (fBkgrdFileNames) nStreams = fBkgrdFileNames->GetEntriesFast() + 1;
915 Int_t signalPerBkgrd = GetNSignalPerBkgrd();
85a5290f 916 AliRunDigitizer* manager = new AliRunDigitizer(nStreams, signalPerBkgrd);
917 manager->SetInputStream(0, fGAliceFileName.Data());
918 for (Int_t iStream = 1; iStream < nStreams; iStream++) {
919 const char* fileName = ((TObjString*)
920 (fBkgrdFileNames->At(iStream-1)))->GetName();
921 manager->SetInputStream(iStream, fileName);
922 }
923
924 TString detStr = detectors;
925 TString detExcl = excludeDetectors;
95601147 926 manager->GetInputStream(0)->ImportgAlice();
927 AliRunLoader* runLoader =
928 AliRunLoader::GetRunLoader(manager->GetInputStream(0)->GetFolderName());
929 TObjArray* detArray = runLoader->GetAliRun()->Detectors();
85a5290f 930 for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
931 AliModule* det = (AliModule*) detArray->At(iDet);
932 if (!det || !det->IsActive()) continue;
933 if (IsSelected(det->GetName(), detStr) &&
934 !IsSelected(det->GetName(), detExcl)) {
4df28b43 935 AliDigitizer* digitizer = det->CreateDigitizer(manager);
936 if (!digitizer) {
ac1671c0 937 AliError(Form("no digitizer for %s", det->GetName()));
85a5290f 938 if (fStopOnError) return kFALSE;
4df28b43 939 } else {
940 digitizer->SetRegionOfInterest(fRegionOfInterest);
85a5290f 941 }
942 }
943 }
944
945 if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
ac1671c0 946 AliError(Form("the following detectors were not found: %s",
947 detStr.Data()));
85a5290f 948 if (fStopOnError) return kFALSE;
949 }
950
951 if (!manager->GetListOfTasks()->IsEmpty()) {
ac1671c0 952 AliInfo("executing digitization");
85a5290f 953 manager->Exec("");
954 }
95601147 955
85a5290f 956 delete manager;
957
5f8272e1 958 AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
959 stopwatch.RealTime(),stopwatch.CpuTime()));
960
85a5290f 961 return kTRUE;
962}
963
964//_____________________________________________________________________________
95601147 965Bool_t AliSimulation::RunHitsDigitization(const char* detectors)
85a5290f 966{
967// run the digitization and produce digits from hits
968
4df28b43 969 TStopwatch stopwatch;
970 stopwatch.Start();
971
591f51bb 972 AliRunLoader* runLoader = LoadRun("READ");
95601147 973 if (!runLoader) return kFALSE;
974
85a5290f 975 TString detStr = detectors;
95601147 976 TObjArray* detArray = runLoader->GetAliRun()->Detectors();
85a5290f 977 for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
978 AliModule* det = (AliModule*) detArray->At(iDet);
979 if (!det || !det->IsActive()) continue;
980 if (IsSelected(det->GetName(), detStr)) {
ac1671c0 981 AliInfo(Form("creating digits from hits for %s", det->GetName()));
85a5290f 982 det->Hits2Digits();
983 }
984 }
985
986 if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
ac1671c0 987 AliError(Form("the following detectors were not found: %s",
988 detStr.Data()));
85a5290f 989 if (fStopOnError) return kFALSE;
990 }
991
b17cac0f 992 delete runLoader;
591f51bb 993 //PH Temporary fix to avoid interference with the PHOS loder/getter
994 //PH The problem has to be solved in more general way 09/06/05
95601147 995
5f8272e1 996 AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
997 stopwatch.RealTime(),stopwatch.CpuTime()));
4df28b43 998
85a5290f 999 return kTRUE;
1000}
1001
0421c3d1 1002//_____________________________________________________________________________
06c7e392 1003Bool_t AliSimulation::WriteRawData(const char* detectors,
1004 const char* fileName,
1005 Bool_t deleteIntermediateFiles)
0421c3d1 1006{
1007// convert the digits to raw data
06c7e392 1008// First DDL raw data files for the given detectors are created.
1009// If a file name is given, the DDL files are then converted to a DATE file.
1010// If deleteIntermediateFiles is true, the DDL raw files are deleted
1011// afterwards.
1012// If the file name has the extension ".root", the DATE file is converted
1013// to a root file.
1014// If deleteIntermediateFiles is true, the DATE file is deleted afterwards.
0421c3d1 1015
1016 TStopwatch stopwatch;
1017 stopwatch.Start();
1018
06c7e392 1019 if (!WriteRawFiles(detectors)) {
1020 if (fStopOnError) return kFALSE;
1021 }
1022
1023 TString dateFileName(fileName);
1024 if (!dateFileName.IsNull()) {
1025 Bool_t rootOutput = dateFileName.EndsWith(".root");
1026 if (rootOutput) dateFileName += ".date";
1027 if (!ConvertRawFilesToDate(dateFileName)) {
1028 if (fStopOnError) return kFALSE;
1029 }
1030 if (deleteIntermediateFiles) {
1031 AliRunLoader* runLoader = LoadRun("READ");
1032 if (runLoader) for (Int_t iEvent = 0;
1033 iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
1034 char command[256];
1035 sprintf(command, "rm -r raw%d", iEvent);
1036 gSystem->Exec(command);
1037 }
1038 }
1039
1040 if (rootOutput) {
1041 if (!ConvertDateToRoot(dateFileName, fileName)) {
1042 if (fStopOnError) return kFALSE;
1043 }
1044 if (deleteIntermediateFiles) {
1045 gSystem->Unlink(dateFileName);
1046 }
1047 }
1048 }
1049
5f8272e1 1050 AliInfo(Form("Execution time: R:%.2fs C:%.2fs",
1051 stopwatch.RealTime(),stopwatch.CpuTime()));
06c7e392 1052
1053 return kTRUE;
1054}
1055
1056//_____________________________________________________________________________
1057Bool_t AliSimulation::WriteRawFiles(const char* detectors)
1058{
1059// convert the digits to raw data DDL files
1060
1061 AliRunLoader* runLoader = LoadRun("READ");
0421c3d1 1062 if (!runLoader) return kFALSE;
1063
06c7e392 1064 // write raw data to DDL files
0421c3d1 1065 for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
ac1671c0 1066 AliInfo(Form("processing event %d", iEvent));
0421c3d1 1067 runLoader->GetEvent(iEvent);
1068 TString baseDir = gSystem->WorkingDirectory();
1069 char dirName[256];
1070 sprintf(dirName, "raw%d", iEvent);
1071 gSystem->MakeDirectory(dirName);
1072 if (!gSystem->ChangeDirectory(dirName)) {
ac1671c0 1073 AliError(Form("couldn't change to directory %s", dirName));
0421c3d1 1074 if (fStopOnError) return kFALSE; else continue;
1075 }
1076
1077 TString detStr = detectors;
1078 TObjArray* detArray = runLoader->GetAliRun()->Detectors();
1079 for (Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++) {
1080 AliModule* det = (AliModule*) detArray->At(iDet);
1081 if (!det || !det->IsActive()) continue;
1082 if (IsSelected(det->GetName(), detStr)) {
ac1671c0 1083 AliInfo(Form("creating raw data from digits for %s", det->GetName()));
0421c3d1 1084 det->Digits2Raw();
1085 }
1086 }
1087
007155cc 1088 if (!WriteTriggerRawData())
1089 if (fStopOnError) return kFALSE;
1090
0421c3d1 1091 gSystem->ChangeDirectory(baseDir);
1092 if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
ac1671c0 1093 AliError(Form("the following detectors were not found: %s",
1094 detStr.Data()));
0421c3d1 1095 if (fStopOnError) return kFALSE;
1096 }
1097 }
1098
1099 delete runLoader;
06c7e392 1100 return kTRUE;
1101}
0421c3d1 1102
06c7e392 1103//_____________________________________________________________________________
1104Bool_t AliSimulation::ConvertRawFilesToDate(const char* dateFileName)
1105{
1106// convert raw data DDL files to a DATE file with the program "dateStream"
1107
06c7e392 1108 char* path = gSystem->Which(gSystem->Getenv("PATH"), "dateStream");
1109 if (!path) {
ac1671c0 1110 AliError("the program dateStream was not found");
06c7e392 1111 if (fStopOnError) return kFALSE;
1112 } else {
1113 delete[] path;
1114 }
0421c3d1 1115
06c7e392 1116 AliRunLoader* runLoader = LoadRun("READ");
1117 if (!runLoader) return kFALSE;
1118
ac1671c0 1119 AliInfo(Form("converting raw data DDL files to DATE file %s", dateFileName));
06c7e392 1120 char command[256];
7e621354 1121 sprintf(command, "dateStream -D -o %s -# %d -C",
06c7e392 1122 dateFileName, runLoader->GetNumberOfEvents());
1123 FILE* pipe = gSystem->OpenPipe(command, "w");
1124
1125 for (Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++) {
1126 fprintf(pipe, "GDC\n");
1127 Float_t ldc = 0;
1128 Int_t prevLDC = -1;
1129
1130 // loop over detectors and DDLs
362c9d61 1131 for (Int_t iDet = 0; iDet < AliDAQ::kNDetectors; iDet++) {
1132 for (Int_t iDDL = 0; iDDL < AliDAQ::NumberOfDdls(iDet); iDDL++) {
06c7e392 1133
362c9d61 1134 Int_t ddlID = AliDAQ::DdlID(iDet,iDDL);
06c7e392 1135 Int_t ldcID = Int_t(ldc + 0.0001);
362c9d61 1136 ldc += AliDAQ::NumberOfLdcs(iDet) / AliDAQ::NumberOfDdls(iDet);
06c7e392 1137
1138 char rawFileName[256];
362c9d61 1139 sprintf(rawFileName, "raw%d/%s",
1140 iEvent, AliDAQ::DdlFileName(iDet,iDDL));
06c7e392 1141
1142 // check existence and size of raw data file
1143 FILE* file = fopen(rawFileName, "rb");
1144 if (!file) continue;
1145 fseek(file, 0, SEEK_END);
1146 unsigned long size = ftell(file);
aa04f81b 1147 fclose(file);
1148 if (!size) continue;
06c7e392 1149
1150 if (ldcID != prevLDC) {
1151 fprintf(pipe, " LDC Id %d\n", ldcID);
1152 prevLDC = ldcID;
1153 }
1154 fprintf(pipe, " Equipment Id %d Payload %s\n", ddlID, rawFileName);
1155 }
1156 }
1157 }
1158
1159 Int_t result = gSystem->ClosePipe(pipe);
1160
1161 delete runLoader;
1162 return (result == 0);
1163}
1164
1165//_____________________________________________________________________________
1166Bool_t AliSimulation::ConvertDateToRoot(const char* dateFileName,
1167 const char* rootFileName)
1168{
1169// convert a DATE file to a root file with the program "alimdc"
1170
1171 // ALIMDC setup
1172 const Int_t kDBSize = 1000000000;
c83eb8f4 1173 const Int_t kTagDBSize = 1000000000;
06c7e392 1174 const Bool_t kFilter = kFALSE;
1175 const Int_t kCompression = 1;
1176
1177 char* path = gSystem->Which(gSystem->Getenv("PATH"), "alimdc");
1178 if (!path) {
ac1671c0 1179 AliError("the program alimdc was not found");
06c7e392 1180 if (fStopOnError) return kFALSE;
1181 } else {
1182 delete[] path;
1183 }
1184
ac1671c0 1185 AliInfo(Form("converting DATE file %s to root file %s",
1186 dateFileName, rootFileName));
06c7e392 1187
f732bd7d 1188 const char* rawDBFS[2] = { "/tmp/mdc1", "/tmp/mdc2" };
1189 const char* tagDBFS = "/tmp/mdc1/tags";
1190 const char* runDBFS = "/tmp/mdc1/meta";
1191
1192 // User defined file system locations
1193 if (gSystem->Getenv("ALIMDC_RAWDB1"))
1194 rawDBFS[0] = gSystem->Getenv("ALIMDC_RAWDB1");
1195 if (gSystem->Getenv("ALIMDC_RAWDB2"))
1196 rawDBFS[1] = gSystem->Getenv("ALIMDC_RAWDB2");
1197 if (gSystem->Getenv("ALIMDC_TAGDB"))
1198 tagDBFS = gSystem->Getenv("ALIMDC_TAGDB");
1199 if (gSystem->Getenv("ALIMDC_RUNDB"))
1200 runDBFS = gSystem->Getenv("ALIMDC_RUNDB");
1201
1202 gSystem->Exec(Form("rm -rf %s",rawDBFS[0]));
1203 gSystem->Exec(Form("rm -rf %s",rawDBFS[1]));
1204 gSystem->Exec(Form("rm -rf %s",tagDBFS));
1205 gSystem->Exec(Form("rm -rf %s",runDBFS));
1206
1207 gSystem->Exec(Form("mkdir %s",rawDBFS[0]));
1208 gSystem->Exec(Form("mkdir %s",rawDBFS[1]));
1209 gSystem->Exec(Form("mkdir %s",tagDBFS));
1210 gSystem->Exec(Form("mkdir %s",runDBFS));
1211
1212 Int_t result = gSystem->Exec(Form("alimdc %d %d %d %d %s",
1213 kDBSize, kTagDBSize, kFilter, kCompression, dateFileName));
1214 gSystem->Exec(Form("mv %s/*.root %s", rawDBFS[0], rootFileName));
1215
1216 gSystem->Exec(Form("rm -rf %s",rawDBFS[0]));
1217 gSystem->Exec(Form("rm -rf %s",rawDBFS[1]));
1218 gSystem->Exec(Form("rm -rf %s",tagDBFS));
1219 gSystem->Exec(Form("rm -rf %s",runDBFS));
06c7e392 1220
1221 return (result == 0);
0421c3d1 1222}
1223
85a5290f 1224
95601147 1225//_____________________________________________________________________________
06c7e392 1226AliRunLoader* AliSimulation::LoadRun(const char* mode) const
95601147 1227{
1228// delete existing run loaders, open a new one and load gAlice
1229
1230 while (AliRunLoader::GetRunLoader()) delete AliRunLoader::GetRunLoader();
1231 AliRunLoader* runLoader =
1232 AliRunLoader::Open(fGAliceFileName.Data(),
06c7e392 1233 AliConfig::GetDefaultEventFolderName(), mode);
95601147 1234 if (!runLoader) {
ac1671c0 1235 AliError(Form("no run loader found in file %s", fGAliceFileName.Data()));
95601147 1236 return NULL;
1237 }
1238 runLoader->LoadgAlice();
1239 gAlice = runLoader->GetAliRun();
1240 if (!gAlice) {
ac1671c0 1241 AliError(Form("no gAlice object found in file %s",
1242 fGAliceFileName.Data()));
95601147 1243 return NULL;
1244 }
1245 return runLoader;
1246}
1247
85a5290f 1248//_____________________________________________________________________________
05526d44 1249Int_t AliSimulation::GetNSignalPerBkgrd(Int_t nEvents) const
1250{
1251// get or calculate the number of signal events per background event
1252
1253 if (!fBkgrdFileNames) return 1;
1254 Int_t nBkgrdFiles = fBkgrdFileNames->GetEntriesFast();
1255 if (nBkgrdFiles == 0) return 1;
1256
1257 // get the number of signal events
1258 if (nEvents <= 0) {
1259 AliRunLoader* runLoader =
1260 AliRunLoader::Open(fGAliceFileName.Data(), "SIGNAL");
1261 if (!runLoader) return 1;
1262 nEvents = runLoader->GetNumberOfEvents();
1263 delete runLoader;
1264 }
1265
1266 Int_t result = 0;
1267 for (Int_t iBkgrdFile = 0; iBkgrdFile < nBkgrdFiles; iBkgrdFile++) {
1268 // get the number of background events
1269 const char* fileName = ((TObjString*)
1270 (fBkgrdFileNames->At(iBkgrdFile)))->GetName();
1271 AliRunLoader* runLoader =
1272 AliRunLoader::Open(fileName, "BKGRD");
1273 if (!runLoader) continue;
1274 Int_t nBkgrdEvents = runLoader->GetNumberOfEvents();
1275 delete runLoader;
1276
1277 // get or calculate the number of signal per background events
1278 Int_t nSignalPerBkgrd = fBkgrdFileNames->At(iBkgrdFile)->GetUniqueID();
1279 if (nSignalPerBkgrd <= 0) {
1280 nSignalPerBkgrd = (nEvents-1) / nBkgrdEvents + 1;
1281 } else if (result && (result != nSignalPerBkgrd)) {
ac1671c0 1282 AliInfo(Form("the number of signal events per background event "
1283 "will be changed from %d to %d for stream %d",
1284 nSignalPerBkgrd, result, iBkgrdFile+1));
05526d44 1285 nSignalPerBkgrd = result;
1286 }
1287
1288 if (!result) result = nSignalPerBkgrd;
1289 if (nSignalPerBkgrd * nBkgrdEvents < nEvents) {
ac1671c0 1290 AliWarning(Form("not enough background events (%d) for %d signal events "
1291 "using %d signal per background events for stream %d",
1292 nBkgrdEvents, nEvents, nSignalPerBkgrd, iBkgrdFile+1));
05526d44 1293 }
1294 }
1295
1296 return result;
1297}
1298
1299//_____________________________________________________________________________
85a5290f 1300Bool_t AliSimulation::IsSelected(TString detName, TString& detectors) const
1301{
1302// check whether detName is contained in detectors
1303// if yes, it is removed from detectors
1304
1305 // check if all detectors are selected
1306 if ((detectors.CompareTo("ALL") == 0) ||
1307 detectors.BeginsWith("ALL ") ||
1308 detectors.EndsWith(" ALL") ||
1309 detectors.Contains(" ALL ")) {
1310 detectors = "ALL";
1311 return kTRUE;
1312 }
1313
1314 // search for the given detector
1315 Bool_t result = kFALSE;
1316 if ((detectors.CompareTo(detName) == 0) ||
1317 detectors.BeginsWith(detName+" ") ||
1318 detectors.EndsWith(" "+detName) ||
1319 detectors.Contains(" "+detName+" ")) {
1320 detectors.ReplaceAll(detName, "");
1321 result = kTRUE;
1322 }
1323
1324 // clean up the detectors string
1325 while (detectors.Contains(" ")) detectors.ReplaceAll(" ", " ");
1326 while (detectors.BeginsWith(" ")) detectors.Remove(0, 1);
1327 while (detectors.EndsWith(" ")) detectors.Remove(detectors.Length()-1, 1);
1328
1329 return result;
1330}