]> git.uio.no Git - u/mrichter/AliRoot.git/blame - STEER/AliReconstruction.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
remove unwanted AliCodeTimerxxx to stay compatible with release
[u/mrichter/AliRoot.git] / STEER / AliReconstruction.cxx
CommitLineData
596a855f 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
fc07289e 3 * *
596a855f 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 the reconstruction //
21// //
22// Clusters and tracks are created for all detectors and all events by //
23// typing: //
24// //
25// AliReconstruction rec; //
26// rec.Run(); //
27// //
28// The Run method returns kTRUE in case of successful execution. //
c71de921 29// //
30// If the input to the reconstruction are not simulated digits but raw data, //
31// this can be specified by an argument of the Run method or by the method //
32// //
33// rec.SetInput("..."); //
34// //
35// The input formats and the corresponding argument are: //
36// - DDL raw data files: directory name, ends with "/" //
37// - raw data root file: root file name, extension ".root" //
38// - raw data DATE file: DATE file name, any other non-empty string //
39// - MC root files : empty string, default //
40// //
b26c3770 41// By default all events are reconstructed. The reconstruction can be //
42// limited to a range of events by giving the index of the first and the //
43// last event as an argument to the Run method or by calling //
44// //
45// rec.SetEventRange(..., ...); //
46// //
47// The index -1 (default) can be used for the last event to indicate no //
48// upper limit of the event range. //
49// //
973388c2 50// In case of raw-data reconstruction the user can modify the default //
51// number of events per digits/clusters/tracks file. In case the option //
52// is not used the number is set 1. In case the user provides 0, than //
53// the number of events is equal to the number of events inside the //
54// raw-data file (i.e. one digits/clusters/tracks file): //
55// //
56// rec.SetNumberOfEventsPerFile(...); //
57// //
58// //
596a855f 59// The name of the galice file can be changed from the default //
e583c30d 60// "galice.root" by passing it as argument to the AliReconstruction //
61// constructor or by //
596a855f 62// //
63// rec.SetGAliceFile("..."); //
64// //
59697224 65// The local reconstruction can be switched on or off for individual //
66// detectors by //
596a855f 67// //
59697224 68// rec.SetRunLocalReconstruction("..."); //
596a855f 69// //
70// The argument is a (case sensitive) string with the names of the //
71// detectors separated by a space. The special string "ALL" selects all //
72// available detectors. This is the default. //
73// //
c71de921 74// The reconstruction of the primary vertex position can be switched off by //
75// //
76// rec.SetRunVertexFinder(kFALSE); //
77// //
b8cd5251 78// The tracking and the creation of ESD tracks can be switched on for //
79// selected detectors by //
596a855f 80// //
b8cd5251 81// rec.SetRunTracking("..."); //
596a855f 82// //
c84a5e9e 83// Uniform/nonuniform field tracking switches (default: uniform field) //
84// //
1d99986f 85// rec.SetUniformFieldTracking(); ( rec.SetUniformFieldTracking(kFALSE); ) //
c84a5e9e 86// //
596a855f 87// The filling of additional ESD information can be steered by //
88// //
89// rec.SetFillESD("..."); //
90// //
b8cd5251 91// Again, for both methods the string specifies the list of detectors. //
92// The default is "ALL". //
93// //
94// The call of the shortcut method //
95// //
96// rec.SetRunReconstruction("..."); //
97// //
98// is equivalent to calling SetRunLocalReconstruction, SetRunTracking and //
99// SetFillESD with the same detector selecting string as argument. //
596a855f 100// //
c71de921 101// The reconstruction requires digits or raw data as input. For the creation //
102// of digits and raw data have a look at the class AliSimulation. //
596a855f 103// //
cd0b062e 104// The input data of a detector can be replaced by the corresponding HLT //
105// data by calling (usual detector string) //
106// SetUseHLTData("..."); //
107// //
24f7a148 108// //
596a855f 109///////////////////////////////////////////////////////////////////////////////
110
f7a1cc68 111#include <TArrayD.h>
024a7e64 112#include <TArrayF.h>
f7a1cc68 113#include <TArrayS.h>
114#include <TChain.h>
024a7e64 115#include <TFile.h>
f7a1cc68 116#include <TGeoGlobalMagField.h>
3103d196 117#include <TGeoManager.h>
f7a1cc68 118#include <TList.h>
2bdb9d38 119#include <TLorentzVector.h>
325aa76f 120#include <TMap.h>
f7a1cc68 121#include <TObjArray.h>
122#include <TPRegexp.h>
123#include <TParameter.h>
124#include <TPluginManager.h>
4b71572b 125#include <TProof.h>
db4aeca1 126#include <TProofOutputFile.h>
f7a1cc68 127#include <TROOT.h>
128#include <TSystem.h>
77ba28ba 129#include <THashTable.h>
6c6f2624 130#include <TGrid.h>
131#include <TMessage.h>
596a855f 132
f7a1cc68 133#include "AliAlignObj.h"
134#include "AliCDBEntry.h"
135#include "AliCDBManager.h"
136#include "AliCDBStorage.h"
137#include "AliCTPRawStream.h"
138#include "AliCascadeVertexer.h"
139#include "AliCentralTrigger.h"
87932dab 140#include "AliCodeTimer.h"
f7a1cc68 141#include "AliDAQ.h"
142#include "AliDetectorRecoParam.h"
143#include "AliESDCaloCells.h"
144#include "AliESDCaloCluster.h"
af885e0f 145#include "AliESDEvent.h"
faffd83e 146#include "AliESDMuonTrack.h"
f7a1cc68 147#include "AliESDPmdTrack.h"
148#include "AliESDTagCreator.h"
2257f27e 149#include "AliESDVertex.h"
faffd83e 150#include "AliESDcascade.h"
f7a1cc68 151#include "AliESDfriend.h"
faffd83e 152#include "AliESDkink.h"
596a855f 153#include "AliESDpid.h"
ff8bb5ae 154#include "AliESDtrack.h"
f7a1cc68 155#include "AliESDtrack.h"
156#include "AliEventInfo.h"
157#include "AliGRPObject.h"
158#include "AliGRPRecoParam.h"
159#include "AliGenEventHeader.h"
25be1e5c 160#include "AliGeomManager.h"
aa3c69a9 161#include "AliGlobalQADataMaker.h"
f7a1cc68 162#include "AliHeader.h"
163#include "AliLog.h"
164#include "AliMagF.h"
165#include "AliMultiplicity.h"
166#include "AliPID.h"
167#include "AliPlaneEff.h"
4e25ac79 168#include "AliQAv1.h"
f7a1cc68 169#include "AliQADataMakerRec.h"
b03591ab 170#include "AliQAManager.h"
33314186 171#include "AliRawVEvent.h"
f7a1cc68 172#include "AliRawEventHeaderBase.h"
cd0b062e 173#include "AliRawHLTManager.h"
f7a1cc68 174#include "AliRawReaderDate.h"
175#include "AliRawReaderFile.h"
176#include "AliRawReaderRoot.h"
177#include "AliReconstruction.h"
178#include "AliReconstructor.h"
179#include "AliRun.h"
7e88424f 180#include "AliRunInfo.h"
f7a1cc68 181#include "AliRunLoader.h"
182#include "AliSysInfo.h" // memory snapshots
183#include "AliTrackPointArray.h"
184#include "AliTracker.h"
185#include "AliTriggerClass.h"
186#include "AliTriggerCluster.h"
a6dd87ad 187#include "AliTriggerIR.h"
f7a1cc68 188#include "AliTriggerConfiguration.h"
189#include "AliV0vertexer.h"
190#include "AliVertexer.h"
191#include "AliVertexerTracks.h"
52dd4a8c 192#include "AliTriggerRunScalers.h"
193#include "AliCTPTimeParams.h"
44e45fac 194
596a855f 195ClassImp(AliReconstruction)
196
c757bafd 197//_____________________________________________________________________________
ac4a7581 198const char* AliReconstruction::fgkDetectorName[AliReconstruction::kNDetectors] = {"ITS", "TPC", "TRD", "TOF", "PHOS", "HMPID", "EMCAL", "MUON", "FMD", "ZDC", "PMD", "T0", "VZERO", "ACORDE", "HLT"};
c757bafd 199
596a855f 200//_____________________________________________________________________________
4b71572b 201AliReconstruction::AliReconstruction(const char* gAliceFilename) :
202 TSelector(),
2257f27e 203 fRunVertexFinder(kTRUE),
a84e2607 204 fRunVertexFinderTracks(kTRUE),
1f46a9ae 205 fRunHLTTracking(kFALSE),
e66fbafb 206 fRunMuonTracking(kFALSE),
d1683eef 207 fRunV0Finder(kTRUE),
208 fRunCascadeFinder(kTRUE),
1d99986f 209 fStopOnError(kFALSE),
210 fWriteAlignmentData(kFALSE),
211 fWriteESDfriend(kFALSE),
b647652d 212 fFillTriggerESD(kTRUE),
1d99986f 213
7f68891d 214 fCleanESD(kTRUE),
a023d8d8 215 fV0DCAmax(3.),
216 fV0CsPmin(0.),
7f68891d 217 fDmax(50.),
218 fZmax(50.),
219
1d99986f 220 fRunLocalReconstruction("ALL"),
b8cd5251 221 fRunTracking("ALL"),
e583c30d 222 fFillESD("ALL"),
7d566c20 223 fLoadCDB(""),
48ce48d1 224 fUseTrackingErrorsForAlignment(""),
e583c30d 225 fGAliceFileName(gAliceFilename),
4b71572b 226 fRawInput(""),
975320a0 227 fESDOutput(""),
81d57268 228 fProofOutputFileName(""),
229 fProofOutputLocation(""),
230 fProofOutputDataset(kFALSE),
231 fProofOutputArchive(""),
35042093 232 fEquipIdMap(""),
b26c3770 233 fFirstEvent(0),
234 fLastEvent(-1),
9d705769 235 fNumberOfEventsPerFile((UInt_t)-1),
b8cd5251 236 fOptions(),
6bae477a 237 fLoadAlignFromCDB(kTRUE),
238 fLoadAlignData("ALL"),
cd0b062e 239 fUseHLTData(),
7e88424f 240 fRunInfo(NULL),
241 fEventInfo(),
522fdd91 242 fRunScalers(NULL),
52dd4a8c 243 fCTPTimeParams(NULL),
e583c30d 244
245 fRunLoader(NULL),
b649205a 246 fRawReader(NULL),
cd0b062e 247 fParentRawReader(NULL),
b8cd5251 248
7e88424f 249 fRecoParam(),
250
58e8dc31 251 fSPDTrackleter(NULL),
252
f2a195c1 253 fDiamondProfileSPD(NULL),
9178838a 254 fDiamondProfile(NULL),
43c9dae1 255 fDiamondProfileTPC(NULL),
77ba28ba 256 fListOfCosmicTriggers(NULL),
f2a195c1 257
6b6e4472 258 fGRPData(NULL),
444753c6 259
6bae477a 260 fAlignObjArray(NULL),
795e4a22 261 fCDBUri(),
f1c1204d 262 fQARefUri(),
759c1df1 263 fSpecCDBUri(),
795e4a22 264 fInitCDBCalled(kFALSE),
265 fSetRunNumberFromDataCalled(kFALSE),
ce43afbe 266 fQADetectors("ALL"),
267 fQATasks("ALL"),
7e963665 268 fRunQA(kTRUE),
6b150027 269 fRunGlobalQA(kTRUE),
e4a998ed 270 fSameQACycle(kFALSE),
f1c1204d 271 fInitQACalled(kFALSE),
75373542 272 fWriteQAExpertData(kTRUE),
21a3aa09 273 fRunPlaneEff(kFALSE),
274
275 fesd(NULL),
276 fhltesd(NULL),
277 fesdf(NULL),
278 ffile(NULL),
279 ftree(NULL),
280 fhlttree(NULL),
21a3aa09 281 ftVertexer(NULL),
14dd053c 282 fIsNewRunLoader(kFALSE),
4b71572b 283 fRunAliEVE(kFALSE),
284 fChain(NULL)
596a855f 285{
286// create reconstruction object with default parameters
002c9d1b 287 gGeoManager = NULL;
b8cd5251 288
ac4a7581 289 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 290 fReconstructor[iDet] = NULL;
291 fLoader[iDet] = NULL;
292 fTracker[iDet] = NULL;
293 }
4e25ac79 294 for (Int_t iDet = 0; iDet < AliQAv1::kNDET; iDet++) {
252f8aa8 295 fQACycles[iDet] = 999999 ;
296 fQAWriteExpert[iDet] = kFALSE ;
297 }
298
e47c4c2e 299 AliPID pid;
596a855f 300}
301
302//_____________________________________________________________________________
303AliReconstruction::AliReconstruction(const AliReconstruction& rec) :
4b71572b 304 TSelector(),
2257f27e 305 fRunVertexFinder(rec.fRunVertexFinder),
a84e2607 306 fRunVertexFinderTracks(rec.fRunVertexFinderTracks),
1f46a9ae 307 fRunHLTTracking(rec.fRunHLTTracking),
e66fbafb 308 fRunMuonTracking(rec.fRunMuonTracking),
d1683eef 309 fRunV0Finder(rec.fRunV0Finder),
310 fRunCascadeFinder(rec.fRunCascadeFinder),
1d99986f 311 fStopOnError(rec.fStopOnError),
312 fWriteAlignmentData(rec.fWriteAlignmentData),
313 fWriteESDfriend(rec.fWriteESDfriend),
b647652d 314 fFillTriggerESD(rec.fFillTriggerESD),
1d99986f 315
7f68891d 316 fCleanESD(rec.fCleanESD),
a023d8d8 317 fV0DCAmax(rec.fV0DCAmax),
5e5c1aa9 318 fV0CsPmin(rec.fV0CsPmin),
7f68891d 319 fDmax(rec.fDmax),
320 fZmax(rec.fZmax),
321
1d99986f 322 fRunLocalReconstruction(rec.fRunLocalReconstruction),
e583c30d 323 fRunTracking(rec.fRunTracking),
324 fFillESD(rec.fFillESD),
7d566c20 325 fLoadCDB(rec.fLoadCDB),
48ce48d1 326 fUseTrackingErrorsForAlignment(rec.fUseTrackingErrorsForAlignment),
e583c30d 327 fGAliceFileName(rec.fGAliceFileName),
4b71572b 328 fRawInput(rec.fRawInput),
975320a0 329 fESDOutput(rec.fESDOutput),
81d57268 330 fProofOutputFileName(rec.fProofOutputFileName),
331 fProofOutputLocation(rec.fProofOutputLocation),
332 fProofOutputDataset(rec.fProofOutputDataset),
333 fProofOutputArchive(rec.fProofOutputArchive),
35042093 334 fEquipIdMap(rec.fEquipIdMap),
b26c3770 335 fFirstEvent(rec.fFirstEvent),
336 fLastEvent(rec.fLastEvent),
973388c2 337 fNumberOfEventsPerFile(rec.fNumberOfEventsPerFile),
b8cd5251 338 fOptions(),
6bae477a 339 fLoadAlignFromCDB(rec.fLoadAlignFromCDB),
340 fLoadAlignData(rec.fLoadAlignData),
2972d4eb 341 fUseHLTData(rec.fUseHLTData),
7e88424f 342 fRunInfo(NULL),
343 fEventInfo(),
522fdd91 344 fRunScalers(NULL),
52dd4a8c 345 fCTPTimeParams(NULL),
e583c30d 346
347 fRunLoader(NULL),
b649205a 348 fRawReader(NULL),
2972d4eb 349 fParentRawReader(NULL),
b8cd5251 350
4b71572b 351 fRecoParam(rec.fRecoParam),
7e88424f 352
58e8dc31 353 fSPDTrackleter(NULL),
354
f2a195c1 355 fDiamondProfileSPD(rec.fDiamondProfileSPD),
4b71572b 356 fDiamondProfile(rec.fDiamondProfile),
357 fDiamondProfileTPC(rec.fDiamondProfileTPC),
77ba28ba 358 fListOfCosmicTriggers(NULL),
f2a195c1 359
6b6e4472 360 fGRPData(NULL),
444753c6 361
6bae477a 362 fAlignObjArray(rec.fAlignObjArray),
ec92bee0 363 fCDBUri(rec.fCDBUri),
f1c1204d 364 fQARefUri(rec.fQARefUri),
7e963665 365 fSpecCDBUri(),
795e4a22 366 fInitCDBCalled(rec.fInitCDBCalled),
367 fSetRunNumberFromDataCalled(rec.fSetRunNumberFromDataCalled),
ce43afbe 368 fQADetectors(rec.fQADetectors),
369 fQATasks(rec.fQATasks),
aa3c69a9 370 fRunQA(rec.fRunQA),
371 fRunGlobalQA(rec.fRunGlobalQA),
e4a998ed 372 fSameQACycle(rec.fSameQACycle),
f1c1204d 373 fInitQACalled(rec.fInitQACalled),
75373542 374 fWriteQAExpertData(rec.fWriteQAExpertData),
21a3aa09 375 fRunPlaneEff(rec.fRunPlaneEff),
376
377 fesd(NULL),
378 fhltesd(NULL),
379 fesdf(NULL),
380 ffile(NULL),
381 ftree(NULL),
382 fhlttree(NULL),
21a3aa09 383 ftVertexer(NULL),
14dd053c 384 fIsNewRunLoader(rec.fIsNewRunLoader),
4b71572b 385 fRunAliEVE(kFALSE),
386 fChain(NULL)
596a855f 387{
388// copy constructor
389
ec92bee0 390 for (Int_t i = 0; i < rec.fOptions.GetEntriesFast(); i++) {
efd2085e 391 if (rec.fOptions[i]) fOptions.Add(rec.fOptions[i]->Clone());
392 }
ac4a7581 393 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 394 fReconstructor[iDet] = NULL;
395 fLoader[iDet] = NULL;
396 fTracker[iDet] = NULL;
b1af1125 397 }
398
4e25ac79 399 for (Int_t iDet = 0; iDet < AliQAv1::kNDET; iDet++) {
b1af1125 400 fQACycles[iDet] = rec.fQACycles[iDet];
252f8aa8 401 fQAWriteExpert[iDet] = rec.fQAWriteExpert[iDet] ;
402 }
b1af1125 403
ec92bee0 404 for (Int_t i = 0; i < rec.fSpecCDBUri.GetEntriesFast(); i++) {
405 if (rec.fSpecCDBUri[i]) fSpecCDBUri.Add(rec.fSpecCDBUri[i]->Clone());
406 }
596a855f 407}
408
409//_____________________________________________________________________________
410AliReconstruction& AliReconstruction::operator = (const AliReconstruction& rec)
411{
412// assignment operator
4b71572b 413// Used in PROOF mode
414// Be very careful while modifing it!
415// Simple rules to follow:
416// for persistent data members - use their assignment operators
417// for non-persistent ones - do nothing or take the default values from constructor
418// TSelector members should not be touched
419 if(&rec == this) return *this;
420
4b71572b 421 fRunVertexFinder = rec.fRunVertexFinder;
422 fRunVertexFinderTracks = rec.fRunVertexFinderTracks;
423 fRunHLTTracking = rec.fRunHLTTracking;
424 fRunMuonTracking = rec.fRunMuonTracking;
425 fRunV0Finder = rec.fRunV0Finder;
426 fRunCascadeFinder = rec.fRunCascadeFinder;
427 fStopOnError = rec.fStopOnError;
428 fWriteAlignmentData = rec.fWriteAlignmentData;
429 fWriteESDfriend = rec.fWriteESDfriend;
430 fFillTriggerESD = rec.fFillTriggerESD;
431
432 fCleanESD = rec.fCleanESD;
433 fV0DCAmax = rec.fV0DCAmax;
434 fV0CsPmin = rec.fV0CsPmin;
435 fDmax = rec.fDmax;
436 fZmax = rec.fZmax;
437
438 fRunLocalReconstruction = rec.fRunLocalReconstruction;
439 fRunTracking = rec.fRunTracking;
440 fFillESD = rec.fFillESD;
7d566c20 441 fLoadCDB = rec.fLoadCDB;
4b71572b 442 fUseTrackingErrorsForAlignment = rec.fUseTrackingErrorsForAlignment;
443 fGAliceFileName = rec.fGAliceFileName;
444 fRawInput = rec.fRawInput;
975320a0 445 fESDOutput = rec.fESDOutput;
81d57268 446 fProofOutputFileName = rec.fProofOutputFileName;
447 fProofOutputLocation = rec.fProofOutputLocation;
448 fProofOutputDataset = rec.fProofOutputDataset;
449 fProofOutputArchive = rec.fProofOutputArchive;
4b71572b 450 fEquipIdMap = rec.fEquipIdMap;
451 fFirstEvent = rec.fFirstEvent;
452 fLastEvent = rec.fLastEvent;
453 fNumberOfEventsPerFile = rec.fNumberOfEventsPerFile;
454
455 for (Int_t i = 0; i < rec.fOptions.GetEntriesFast(); i++) {
456 if (rec.fOptions[i]) fOptions.Add(rec.fOptions[i]->Clone());
457 }
458
459 fLoadAlignFromCDB = rec.fLoadAlignFromCDB;
460 fLoadAlignData = rec.fLoadAlignData;
4b71572b 461 fUseHLTData = rec.fUseHLTData;
462
463 delete fRunInfo; fRunInfo = NULL;
464 if (rec.fRunInfo) fRunInfo = new AliRunInfo(*rec.fRunInfo);
465
466 fEventInfo = rec.fEventInfo;
467
522fdd91 468 delete fRunScalers; fRunScalers = NULL;
469 if (rec.fRunScalers) fRunScalers = new AliTriggerRunScalers(*rec.fRunScalers);
470
52dd4a8c 471 delete fCTPTimeParams; fCTPTimeParams = NULL;
472 if (rec.fCTPTimeParams) fCTPTimeParams = new AliCTPTimeParams(*rec.fCTPTimeParams);
473
4b71572b 474 fRunLoader = NULL;
475 fRawReader = NULL;
476 fParentRawReader = NULL;
477
478 fRecoParam = rec.fRecoParam;
479
ac4a7581 480 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
4b71572b 481 delete fReconstructor[iDet]; fReconstructor[iDet] = NULL;
482 delete fLoader[iDet]; fLoader[iDet] = NULL;
483 delete fTracker[iDet]; fTracker[iDet] = NULL;
4b71572b 484 }
b1af1125 485
4e25ac79 486 for (Int_t iDet = 0; iDet < AliQAv1::kNDET; iDet++) {
b1af1125 487 fQACycles[iDet] = rec.fQACycles[iDet];
252f8aa8 488 fQAWriteExpert[iDet] = rec.fQAWriteExpert[iDet] ;
489 }
58e8dc31 490
491 delete fSPDTrackleter; fSPDTrackleter = NULL;
252f8aa8 492
f2a195c1 493 delete fDiamondProfileSPD; fDiamondProfileSPD = NULL;
494 if (rec.fDiamondProfileSPD) fDiamondProfileSPD = new AliESDVertex(*rec.fDiamondProfileSPD);
4b71572b 495 delete fDiamondProfile; fDiamondProfile = NULL;
496 if (rec.fDiamondProfile) fDiamondProfile = new AliESDVertex(*rec.fDiamondProfile);
497 delete fDiamondProfileTPC; fDiamondProfileTPC = NULL;
498 if (rec.fDiamondProfileTPC) fDiamondProfileTPC = new AliESDVertex(*rec.fDiamondProfileTPC);
4b71572b 499
77ba28ba 500 delete fListOfCosmicTriggers; fListOfCosmicTriggers = NULL;
501 if (rec.fListOfCosmicTriggers) fListOfCosmicTriggers = (THashTable*)((rec.fListOfCosmicTriggers)->Clone());
502
4b71572b 503 delete fGRPData; fGRPData = NULL;
44e45fac 504 // if (rec.fGRPData) fGRPData = (TMap*)((rec.fGRPData)->Clone());
505 if (rec.fGRPData) fGRPData = (AliGRPObject*)((rec.fGRPData)->Clone());
4b71572b 506
507 delete fAlignObjArray; fAlignObjArray = NULL;
508
509 fCDBUri = "";
f1c1204d 510 fQARefUri = rec.fQARefUri;
4b71572b 511 fSpecCDBUri.Delete();
512 fInitCDBCalled = rec.fInitCDBCalled;
513 fSetRunNumberFromDataCalled = rec.fSetRunNumberFromDataCalled;
514 fQADetectors = rec.fQADetectors;
4b71572b 515 fQATasks = rec.fQATasks;
516 fRunQA = rec.fRunQA;
517 fRunGlobalQA = rec.fRunGlobalQA;
4b71572b 518 fSameQACycle = rec.fSameQACycle;
f1c1204d 519 fInitQACalled = rec.fInitQACalled;
75373542 520 fWriteQAExpertData = rec.fWriteQAExpertData;
4b71572b 521 fRunPlaneEff = rec.fRunPlaneEff;
522
523 fesd = NULL;
524 fhltesd = NULL;
525 fesdf = NULL;
526 ffile = NULL;
527 ftree = NULL;
528 fhlttree = NULL;
529 ftVertexer = NULL;
530 fIsNewRunLoader = rec.fIsNewRunLoader;
531 fRunAliEVE = kFALSE;
532 fChain = NULL;
596a855f 533
596a855f 534 return *this;
535}
536
537//_____________________________________________________________________________
538AliReconstruction::~AliReconstruction()
539{
540// clean up
541
e583c30d 542 CleanUp();
77ba28ba 543 if (fListOfCosmicTriggers) {
544 fListOfCosmicTriggers->Delete();
545 delete fListOfCosmicTriggers;
546 }
4f3c479a 547 delete fGRPData;
522fdd91 548 delete fRunScalers;
52dd4a8c 549 delete fCTPTimeParams;
efd2085e 550 fOptions.Delete();
4b71572b 551 if (fAlignObjArray) {
552 fAlignObjArray->Delete();
553 delete fAlignObjArray;
554 }
ec92bee0 555 fSpecCDBUri.Delete();
e6d66370 556
87932dab 557 AliCodeTimer::Instance()->Print();
596a855f 558}
559
f1c1204d 560//_____________________________________________________________________________
561void AliReconstruction::InitQA()
562{
563 //Initialize the QA and start of cycle
52dd4a8c 564 AliCodeTimerAuto("",0);
f1c1204d 565
566 if (fInitQACalled) return;
567 fInitQACalled = kTRUE;
568
634696f5 569 AliQAManager * qam = AliQAManager::QAManager(AliQAv1::kRECMODE) ;
75373542 570 if (fWriteQAExpertData)
eca4fa66 571 qam->SetWriteExpert() ;
75373542 572
eca4fa66 573 if (qam->IsDefaultStorageSet()) {
f1c1204d 574 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
575 AliWarning("Default QA reference storage has been already set !");
576 AliWarning(Form("Ignoring the default storage declared in AliReconstruction: %s",fQARefUri.Data()));
577 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
eca4fa66 578 fQARefUri = qam->GetDefaultStorage()->GetURI();
f1c1204d 579 } else {
580 if (fQARefUri.Length() > 0) {
581 AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
582 AliDebug(2, Form("Default QA reference storage is set to: %s", fQARefUri.Data()));
583 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
584 } else {
585 fQARefUri="local://$ALICE_ROOT/QAref";
586 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
587 AliWarning("Default QA refeference storage not yet set !!!!");
588 AliWarning(Form("Setting it now to: %s", fQARefUri.Data()));
589 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
590
591 }
eca4fa66 592 qam->SetDefaultStorage(fQARefUri);
f1c1204d 593 }
594
595 if (fRunQA) {
eca4fa66 596 qam->SetActiveDetectors(fQADetectors) ;
4e25ac79 597 for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++) {
eca4fa66 598 qam->SetCycleLength(AliQAv1::DETECTORINDEX_t(det), fQACycles[det]) ;
599 qam->SetWriteExpert(AliQAv1::DETECTORINDEX_t(det)) ;
f1c1204d 600 }
514cb8c7 601 if (!fRawReader && !fInput && IsInTasks(AliQAv1::kRAWS))
4e25ac79 602 fQATasks.ReplaceAll(Form("%d",AliQAv1::kRAWS), "") ;
eca4fa66 603 qam->SetTasks(fQATasks) ;
604 qam->InitQADataMaker(AliCDBManager::Instance()->GetRun()) ;
f1c1204d 605 }
606 if (fRunGlobalQA) {
607 Bool_t sameCycle = kFALSE ;
eca4fa66 608 AliQADataMaker *qadm = qam->GetQADataMaker(AliQAv1::kGLOBAL);
f1c1204d 609 AliInfo(Form("Initializing the global QA data maker"));
514cb8c7 610 if (IsInTasks(AliQAv1::kRECPOINTS)) {
4e25ac79 611 qadm->StartOfCycle(AliQAv1::kRECPOINTS, AliCDBManager::Instance()->GetRun(), sameCycle) ;
612 TObjArray **arr=qadm->Init(AliQAv1::kRECPOINTS);
f1c1204d 613 AliTracker::SetResidualsArray(arr);
614 sameCycle = kTRUE ;
615 }
514cb8c7 616 if (IsInTasks(AliQAv1::kESDS)) {
4e25ac79 617 qadm->StartOfCycle(AliQAv1::kESDS, AliCDBManager::Instance()->GetRun(), sameCycle) ;
618 qadm->Init(AliQAv1::kESDS);
f1c1204d 619 }
f1c1204d 620 }
fc07289e 621 AliSysInfo::AddStamp("InitQA") ;
f1c1204d 622}
623
fc07289e 624//_____________________________________________________________________________
87da0921 625void AliReconstruction::MergeQA(const char *fileName)
fc07289e 626{
627 //Initialize the QA and start of cycle
52dd4a8c 628 AliCodeTimerAuto("",0) ;
eca4fa66 629 AliQAManager::QAManager()->Merge(AliCDBManager::Instance()->GetRun(),fileName) ;
fc07289e 630 AliSysInfo::AddStamp("MergeQA") ;
631}
632
024cf675 633//_____________________________________________________________________________
795e4a22 634void AliReconstruction::InitCDB()
024cf675 635{
636// activate a default CDB storage
637// First check if we have any CDB storage set, because it is used
638// to retrieve the calibration and alignment constants
52dd4a8c 639 AliCodeTimerAuto("",0);
024cf675 640
795e4a22 641 if (fInitCDBCalled) return;
642 fInitCDBCalled = kTRUE;
643
024cf675 644 AliCDBManager* man = AliCDBManager::Instance();
ec92bee0 645 if (man->IsDefaultStorageSet())
024cf675 646 {
647 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
ec92bee0 648 AliWarning("Default CDB storage has been already set !");
649 AliWarning(Form("Ignoring the default storage declared in AliReconstruction: %s",fCDBUri.Data()));
024cf675 650 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
795e4a22 651 fCDBUri = man->GetDefaultStorage()->GetURI();
ec92bee0 652 }
653 else {
795e4a22 654 if (fCDBUri.Length() > 0)
655 {
656 AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
657 AliDebug(2, Form("Default CDB storage is set to: %s", fCDBUri.Data()));
658 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
52dd4a8c 659 man->SetDefaultStorage(fCDBUri);
660 }
661 else if (!man->GetRaw()){
662 fCDBUri="local://$ALICE_ROOT/OCDB";
795e4a22 663 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
664 AliWarning("Default CDB storage not yet set !!!!");
665 AliWarning(Form("Setting it now to: %s", fCDBUri.Data()));
666 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
52dd4a8c 667 man->SetDefaultStorage(fCDBUri);
668 }
669 else {
670 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
671 AliWarning("Default storage will be set after setting the Run Number!!!");
672 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
795e4a22 673 }
ec92bee0 674 }
675
676 // Now activate the detector specific CDB storage locations
c3a7b59a 677 for (Int_t i = 0; i < fSpecCDBUri.GetEntriesFast(); i++) {
678 TObject* obj = fSpecCDBUri[i];
679 if (!obj) continue;
b8ec52f6 680 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
681 AliDebug(2, Form("Specific CDB storage for %s is set to: %s",obj->GetName(),obj->GetTitle()));
682 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
c3a7b59a 683 man->SetSpecificStorage(obj->GetName(), obj->GetTitle());
ec92bee0 684 }
002c9d1b 685 AliSysInfo::AddStamp("InitCDB");
024cf675 686}
687
688//_____________________________________________________________________________
689void AliReconstruction::SetDefaultStorage(const char* uri) {
ec92bee0 690// Store the desired default CDB storage location
691// Activate it later within the Run() method
024cf675 692
ec92bee0 693 fCDBUri = uri;
024cf675 694
695}
696
f1c1204d 697//_____________________________________________________________________________
698void AliReconstruction::SetQARefDefaultStorage(const char* uri) {
699 // Store the desired default CDB storage location
700 // Activate it later within the Run() method
701
702 fQARefUri = uri;
4e25ac79 703 AliQAv1::SetQARefStorage(fQARefUri.Data()) ;
f1c1204d 704
705}
00aa02d5 706//_____________________________________________________________________________
c3a7b59a 707void AliReconstruction::SetSpecificStorage(const char* calibType, const char* uri) {
ec92bee0 708// Store a detector-specific CDB storage location
709// Activate it later within the Run() method
024cf675 710
c3a7b59a 711 AliCDBPath aPath(calibType);
712 if(!aPath.IsValid()){
713 // if calibType is not wildcard but it is a valid detector, add "/*" to make it a valid path
ac4a7581 714 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
c3a7b59a 715 if(!strcmp(calibType, fgkDetectorName[iDet])) {
716 aPath.SetPath(Form("%s/*", calibType));
717 AliInfo(Form("Path for specific storage set to %s", aPath.GetPath().Data()));
718 break;
719 }
720 }
721 if(!aPath.IsValid()){
722 AliError(Form("Not a valid path or detector: %s", calibType));
723 return;
724 }
725 }
726
53dd3c3d 727// // check that calibType refers to a "valid" detector name
728// Bool_t isDetector = kFALSE;
ac4a7581 729// for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
53dd3c3d 730// TString detName = fgkDetectorName[iDet];
731// if(aPath.GetLevel0() == detName) {
732// isDetector = kTRUE;
733// break;
734// }
735// }
736//
737// if(!isDetector) {
738// AliError(Form("Not a valid detector: %s", aPath.GetLevel0().Data()));
739// return;
740// }
c3a7b59a 741
742 TObject* obj = fSpecCDBUri.FindObject(aPath.GetPath().Data());
ec92bee0 743 if (obj) fSpecCDBUri.Remove(obj);
c3a7b59a 744 fSpecCDBUri.Add(new TNamed(aPath.GetPath().Data(), uri));
024cf675 745
746}
747
6bae477a 748//_____________________________________________________________________________
795e4a22 749Bool_t AliReconstruction::SetRunNumberFromData()
6bae477a 750{
751 // The method is called in Run() in order
752 // to set a correct run number.
753 // In case of raw data reconstruction the
754 // run number is taken from the raw data header
755
795e4a22 756 if (fSetRunNumberFromDataCalled) return kTRUE;
757 fSetRunNumberFromDataCalled = kTRUE;
758
759 AliCDBManager* man = AliCDBManager::Instance();
4b71572b 760
761 if(fRawReader) {
762 if(fRawReader->NextEvent()) {
763 if(man->GetRun() > 0) {
21a3aa09 764 AliWarning("Run number is taken from raw-event header! Ignoring settings in AliCDBManager!");
4b71572b 765 }
766 man->SetRun(fRawReader->GetRunNumber());
767 fRawReader->RewindEvents();
768 }
769 else {
770 if(man->GetRun() > 0) {
771 AliWarning("No raw-data events are found ! Using settings in AliCDBManager !");
772 }
773 else {
774 AliWarning("Neither raw events nor settings in AliCDBManager are found !");
775 return kFALSE;
776 }
777 }
778 }
779 else {
780 AliRunLoader *rl = AliRunLoader::Open(fGAliceFileName.Data());
781 if (!rl) {
782 AliError(Form("No run loader found in file %s", fGAliceFileName.Data()));
6bae477a 783 return kFALSE;
784 }
ec92bee0 785 else {
4b71572b 786 rl->LoadHeader();
787 // read run number from gAlice
788 if(rl->GetHeader()) {
789 man->SetRun(rl->GetHeader()->GetRun());
790 rl->UnloadHeader();
791 delete rl;
ec92bee0 792 }
793 else {
4b71572b 794 AliError("Neither run-loader header nor RawReader objects are found !");
795 delete rl;
ec92bee0 796 return kFALSE;
797 }
4b71572b 798 }
6bae477a 799 }
795e4a22 800
801 man->Print();
802
6bae477a 803 return kTRUE;
804}
805
795e4a22 806//_____________________________________________________________________________
807void AliReconstruction::SetCDBLock() {
808 // Set CDB lock: from now on it is forbidden to reset the run number
809 // or the default storage or to activate any further storage!
810
811 AliCDBManager::Instance()->SetLock(1);
812}
813
6bae477a 814//_____________________________________________________________________________
815Bool_t AliReconstruction::MisalignGeometry(const TString& detectors)
816{
817 // Read the alignment objects from CDB.
818 // Each detector is supposed to have the
819 // alignment objects in DET/Align/Data CDB path.
820 // All the detector objects are then collected,
821 // sorted by geometry level (starting from ALIC) and
822 // then applied to the TGeo geometry.
823 // Finally an overlaps check is performed.
824
825 // Load alignment data from CDB and fill fAlignObjArray
826 if(fLoadAlignFromCDB){
6bae477a 827
25be1e5c 828 TString detStr = detectors;
829 TString loadAlObjsListOfDets = "";
830
ac4a7581 831 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
8cb26cdf 832 if(!IsSelected(fgkDetectorName[iDet], detStr)) continue;
f12d42ce 833 if(!strcmp(fgkDetectorName[iDet],"HLT")) continue;
4fbb8e9d 834
835 if(AliGeomManager::GetNalignable(fgkDetectorName[iDet]) != 0)
8cb26cdf 836 {
837 loadAlObjsListOfDets += fgkDetectorName[iDet];
838 loadAlObjsListOfDets += " ";
839 }
25be1e5c 840 } // end loop over detectors
4fbb8e9d 841
842 if(AliGeomManager::GetNalignable("GRP") != 0)
8cb26cdf 843 loadAlObjsListOfDets.Prepend("GRP "); //add alignment objects for non-sensitive modules
98e303d9 844 AliGeomManager::ApplyAlignObjsFromCDB(loadAlObjsListOfDets.Data());
e30a9b4d 845 AliCDBManager::Instance()->UnloadFromCache("*/Align/*");
25be1e5c 846 }else{
847 // Check if the array with alignment objects was
848 // provided by the user. If yes, apply the objects
849 // to the present TGeo geometry
850 if (fAlignObjArray) {
851 if (gGeoManager && gGeoManager->IsClosed()) {
98e303d9 852 if (AliGeomManager::ApplyAlignObjsToGeom(*fAlignObjArray) == kFALSE) {
25be1e5c 853 AliError("The misalignment of one or more volumes failed!"
854 "Compare the list of simulated detectors and the list of detector alignment data!");
855 return kFALSE;
856 }
857 }
858 else {
859 AliError("Can't apply the misalignment! gGeoManager doesn't exist or it is still opened!");
6bae477a 860 return kFALSE;
861 }
862 }
6bae477a 863 }
25be1e5c 864
4b71572b 865 if (fAlignObjArray) {
866 fAlignObjArray->Delete();
867 delete fAlignObjArray; fAlignObjArray=NULL;
868 }
a03b0371 869
6bae477a 870 return kTRUE;
871}
596a855f 872
873//_____________________________________________________________________________
874void AliReconstruction::SetGAliceFile(const char* fileName)
875{
876// set the name of the galice file
877
878 fGAliceFileName = fileName;
879}
880
21a3aa09 881//_____________________________________________________________________________
b58710ef 882void AliReconstruction::SetInput(const char* input)
21a3aa09 883{
b58710ef 884 // In case the input string starts with 'mem://', we run in an online mode
885 // and AliRawReaderDateOnline object is created. In all other cases a raw-data
886 // file is assumed. One can give as an input:
887 // mem://: - events taken from DAQ monitoring libs online
888 // or
889 // mem://<filename> - emulation of the above mode (via DATE monitoring libs)
4b71572b 890 if (input) fRawInput = input;
21a3aa09 891}
892
975320a0 893//_____________________________________________________________________________
894void AliReconstruction::SetOutput(const char* output)
895{
896 // Set the output ESD filename
897 // 'output' is a normalt ROOT url
898 // The method is used in case of raw-data reco with PROOF
81d57268 899 if (output) fESDOutput = output;
975320a0 900}
901
efd2085e 902//_____________________________________________________________________________
903void AliReconstruction::SetOption(const char* detector, const char* option)
904{
905// set options for the reconstruction of a detector
906
907 TObject* obj = fOptions.FindObject(detector);
908 if (obj) fOptions.Remove(obj);
909 fOptions.Add(new TNamed(detector, option));
910}
911
7e88424f 912//_____________________________________________________________________________
913void AliReconstruction::SetRecoParam(const char* detector, AliDetectorRecoParam *par)
914{
915 // Set custom reconstruction parameters for a given detector
916 // Single set of parameters for all the events
a00021a7 917
918 // First check if the reco-params are global
919 if(!strcmp(detector, "GRP")) {
920 par->SetAsDefault();
ac4a7581 921 fRecoParam.AddDetRecoParam(kNDetectors,par);
a00021a7 922 return;
923 }
924
ac4a7581 925 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
7e88424f 926 if(!strcmp(detector, fgkDetectorName[iDet])) {
927 par->SetAsDefault();
928 fRecoParam.AddDetRecoParam(iDet,par);
929 break;
930 }
931 }
932
933}
934
61807e09 935//_____________________________________________________________________________
61807e09 936Bool_t AliReconstruction::InitGRP() {
937 //------------------------------------
938 // Initialization of the GRP entry
939 //------------------------------------
940 AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/GRP/Data");
941
4b71572b 942 if (entry) {
44e45fac 943
944 TMap* m = dynamic_cast<TMap*>(entry->GetObject()); // old GRP entry
945
946 if (m) {
947 AliInfo("Found a TMap in GRP/GRP/Data, converting it into an AliGRPObject");
948 m->Print();
949 fGRPData = new AliGRPObject();
950 fGRPData->ReadValuesFromMap(m);
951 }
952
953 else {
954 AliInfo("Found an AliGRPObject in GRP/GRP/Data, reading it");
955 fGRPData = dynamic_cast<AliGRPObject*>(entry->GetObject()); // new GRP entry
956 entry->SetOwner(0);
957 }
958
f168abba 959 // FIX ME: The unloading of GRP entry is temporarily disabled
960 // because ZDC and VZERO are using it in order to initialize
961 // their reconstructor objects. In the future one has to think
962 // of propagating AliRunInfo to the reconstructors.
963 // AliCDBManager::Instance()->UnloadFromCache("GRP/GRP/Data");
4b71572b 964 }
61807e09 965
966 if (!fGRPData) {
967 AliError("No GRP entry found in OCDB!");
968 return kFALSE;
969 }
970
44e45fac 971 TString lhcState = fGRPData->GetLHCState();
972 if (lhcState==AliGRPObject::GetInvalidString()) {
7e88424f 973 AliError("GRP/GRP/Data entry: missing value for the LHC state ! Using UNKNOWN");
44e45fac 974 lhcState = "UNKNOWN";
7e88424f 975 }
976
44e45fac 977 TString beamType = fGRPData->GetBeamType();
978 if (beamType==AliGRPObject::GetInvalidString()) {
7e88424f 979 AliError("GRP/GRP/Data entry: missing value for the beam type ! Using UNKNOWN");
44e45fac 980 beamType = "UNKNOWN";
7e88424f 981 }
982
44e45fac 983 Float_t beamEnergy = fGRPData->GetBeamEnergy();
984 if (beamEnergy==AliGRPObject::GetInvalidFloat()) {
7e88424f 985 AliError("GRP/GRP/Data entry: missing value for the beam energy ! Using 0");
44e45fac 986 beamEnergy = 0;
7e88424f 987 }
1e500f25 988 // LHC: "multiply by 120 to get the energy in MeV"
989 beamEnergy *= 0.120;
7e88424f 990
44e45fac 991 TString runType = fGRPData->GetRunType();
992 if (runType==AliGRPObject::GetInvalidString()) {
7e88424f 993 AliError("GRP/GRP/Data entry: missing value for the run type ! Using UNKNOWN");
44e45fac 994 runType = "UNKNOWN";
7e88424f 995 }
996
44e45fac 997 Int_t activeDetectors = fGRPData->GetDetectorMask();
a5dede5c 998 if (activeDetectors==AliGRPObject::GetInvalidUInt()) {
7e88424f 999 AliError("GRP/GRP/Data entry: missing value for the detector mask ! Using 1074790399");
44e45fac 1000 activeDetectors = 1074790399;
7e88424f 1001 }
1002
44e45fac 1003 fRunInfo = new AliRunInfo(lhcState, beamType, beamEnergy, runType, activeDetectors);
44e45fac 1004 fRunInfo->Dump();
1005
7e88424f 1006
1007 // Process the list of active detectors
44e45fac 1008 if (activeDetectors) {
1009 UInt_t detMask = activeDetectors;
7a317bef 1010 fRunLocalReconstruction = MatchDetectorList(fRunLocalReconstruction,detMask);
1011 fRunTracking = MatchDetectorList(fRunTracking,detMask);
1012 fFillESD = MatchDetectorList(fFillESD,detMask);
1013 fQADetectors = MatchDetectorList(fQADetectors,detMask);
7d566c20 1014 fLoadCDB.Form("%s %s %s %s",
1015 fRunLocalReconstruction.Data(),
1016 fRunTracking.Data(),
1017 fFillESD.Data(),
1018 fQADetectors.Data());
7d566c20 1019 fLoadCDB = MatchDetectorList(fLoadCDB,detMask);
d0f84c42 1020 if (!((detMask >> AliDAQ::DetectorID("ITSSPD")) & 0x1) &&
1021 !((detMask >> AliDAQ::DetectorID("ITSSDD")) & 0x1) &&
1022 !((detMask >> AliDAQ::DetectorID("ITSSSD")) & 0x1) ) {
91e2f025 1023 // switch off the vertexer
d0f84c42 1024 AliInfo("SPD,SDD,SSD is not in the list of active detectors. Vertexer switched off.");
91e2f025 1025 fRunVertexFinder = kFALSE;
1026 }
304d7f0b 1027 if (!((detMask >> AliDAQ::DetectorID("TRG")) & 0x1)) {
1028 // switch off the reading of CTP raw-data payload
1029 if (fFillTriggerESD) {
1030 AliInfo("CTP is not in the list of active detectors. CTP data reading switched off.");
1031 fFillTriggerESD = kFALSE;
1032 }
1033 }
7e88424f 1034 }
1035
1036 AliInfo("===================================================================================");
1037 AliInfo(Form("Running local reconstruction for detectors: %s",fRunLocalReconstruction.Data()));
1038 AliInfo(Form("Running tracking for detectors: %s",fRunTracking.Data()));
1039 AliInfo(Form("Filling ESD for detectors: %s",fFillESD.Data()));
ca13fb87 1040 AliInfo(Form("Quality assurance is active for detectors: %s",fQADetectors.Data()));
7d566c20 1041 AliInfo(Form("CDB and reconstruction parameters are loaded for detectors: %s",fLoadCDB.Data()));
7e88424f 1042 AliInfo("===================================================================================");
61807e09 1043
1044 //*** Dealing with the magnetic field map
5a004fb4 1045 if ( TGeoGlobalMagField::Instance()->IsLocked() ) {
1046 if (TGeoGlobalMagField::Instance()->GetField()->TestBit(AliMagF::kOverrideGRP)) {
1047 AliInfo("ExpertMode!!! GRP information will be ignored !");
1048 AliInfo("ExpertMode!!! Running with the externally locked B field !");
1049 }
1050 else {
1051 AliInfo("Destroying existing B field instance!");
1052 delete TGeoGlobalMagField::Instance();
1053 }
1054 }
1055 if ( !TGeoGlobalMagField::Instance()->IsLocked() ) {
61807e09 1056 // Construct the field map out of the information retrieved from GRP.
856024f0 1057 Bool_t ok = kTRUE;
61807e09 1058 // L3
44e45fac 1059 Float_t l3Current = fGRPData->GetL3Current((AliGRPObject::Stats)0);
1060 if (l3Current == AliGRPObject::GetInvalidFloat()) {
1061 AliError("GRP/GRP/Data entry: missing value for the L3 current !");
1062 ok = kFALSE;
1063 }
f7a1cc68 1064
44e45fac 1065 Char_t l3Polarity = fGRPData->GetL3Polarity();
1066 if (l3Polarity == AliGRPObject::GetInvalidChar()) {
1067 AliError("GRP/GRP/Data entry: missing value for the L3 polarity !");
1068 ok = kFALSE;
1069 }
1070
1071 // Dipole
1072 Float_t diCurrent = fGRPData->GetDipoleCurrent((AliGRPObject::Stats)0);
1073 if (diCurrent == AliGRPObject::GetInvalidFloat()) {
1074 AliError("GRP/GRP/Data entry: missing value for the dipole current !");
1075 ok = kFALSE;
1076 }
1077
1078 Char_t diPolarity = fGRPData->GetDipolePolarity();
1079 if (diPolarity == AliGRPObject::GetInvalidChar()) {
1080 AliError("GRP/GRP/Data entry: missing value for the dipole polarity !");
1081 ok = kFALSE;
1082 }
1083
99c7d495 1084 // read special bits for the polarity convention and map type
1085 Int_t polConvention = fGRPData->IsPolarityConventionLHC() ? AliMagF::kConvLHC : AliMagF::kConvDCS2008;
1086 Bool_t uniformB = fGRPData->IsUniformBMap();
61807e09 1087
856024f0 1088 if (ok) {
33fe5eb1 1089 AliMagF* fld = AliMagF::CreateFieldMap(TMath::Abs(l3Current) * (l3Polarity ? -1:1),
1090 TMath::Abs(diCurrent) * (diPolarity ? -1:1),
1091 polConvention,uniformB,beamEnergy, beamType.Data());
1092 if (fld) {
1093 TGeoGlobalMagField::Instance()->SetField( fld );
1094 TGeoGlobalMagField::Instance()->Lock();
1095 AliInfo("Running with the B field constructed out of GRP !");
1096 }
1097 else AliFatal("Failed to create a B field map !");
856024f0 1098 }
f7a1cc68 1099 else AliFatal("B field is neither set nor constructed from GRP ! Exitig...");
61807e09 1100 }
f7a1cc68 1101
f2a195c1 1102 //*** Get the diamond profiles from OCDB
1103 entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertexSPD");
1104 if (entry) {
1105 fDiamondProfileSPD = dynamic_cast<AliESDVertex*> (entry->GetObject());
1106 } else {
1107 AliError("No SPD diamond profile found in OCDB!");
1108 }
1109
61807e09 1110 entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertex");
1111 if (entry) {
f2a195c1 1112 fDiamondProfile = dynamic_cast<AliESDVertex*> (entry->GetObject());
61807e09 1113 } else {
1114 AliError("No diamond profile found in OCDB!");
1115 }
1116
1117 entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertexTPC");
1118 if (entry) {
f2a195c1 1119 fDiamondProfileTPC = dynamic_cast<AliESDVertex*> (entry->GetObject());
61807e09 1120 } else {
f2a195c1 1121 AliError("No TPC diamond profile found in OCDB!");
61807e09 1122 }
1123
77ba28ba 1124 entry = AliCDBManager::Instance()->Get("GRP/Calib/CosmicTriggers");
1125 if (entry) {
1126 fListOfCosmicTriggers = dynamic_cast<THashTable*>(entry->GetObject());
1127 entry->SetOwner(0);
1128 AliCDBManager::Instance()->UnloadFromCache("GRP/Calib/CosmicTriggers");
1129 }
1130
1131 if (!fListOfCosmicTriggers) {
1132 AliWarning("Can not get list of cosmic triggers from OCDB! Cosmic event specie will be effectively disabled!");
1133 }
1134
61807e09 1135 return kTRUE;
1136}
1137
4b71572b 1138//_____________________________________________________________________________
1139Bool_t AliReconstruction::LoadCDB()
1140{
52dd4a8c 1141 AliCodeTimerAuto("",0);
4b71572b 1142
1143 AliCDBManager::Instance()->Get("GRP/CTP/Config");
1144
7d566c20 1145 TString detStr = fLoadCDB;
ac4a7581 1146 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
4b71572b 1147 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
1148 AliCDBManager::Instance()->GetAll(Form("%s/Calib/*",fgkDetectorName[iDet]));
1149 }
e122592f 1150
1151 // Temporary fix - one has to define the correct policy in order
1152 // to load the trigger OCDB entries only for the detectors that
1153 // in the trigger or that are needed in order to put correct
1154 // information in ESD
1155 AliCDBManager::Instance()->GetAll("TRIGGER/*/*");
1156
4b71572b 1157 return kTRUE;
1158}
522fdd91 1159//_____________________________________________________________________________
1160Bool_t AliReconstruction::LoadTriggerScalersCDB()
1161{
52dd4a8c 1162 AliCodeTimerAuto("",0);
522fdd91 1163
1164 AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/CTP/Scalers");
4b71572b 1165
522fdd91 1166 if (entry) {
1167
1168 AliInfo("Found an AliTriggerRunScalers in GRP/CTP/Scalers, reading it");
1169 fRunScalers = dynamic_cast<AliTriggerRunScalers*> (entry->GetObject());
1170 entry->SetOwner(0);
82ebedd6 1171 if (fRunScalers->CorrectScalersOverflow() == 0) AliInfo("32bit Trigger counters corrected for overflow");
1172
522fdd91 1173 }
1174 return kTRUE;
1175}
596a855f 1176//_____________________________________________________________________________
52dd4a8c 1177Bool_t AliReconstruction::LoadCTPTimeParamsCDB()
1178{
1179 AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/CTP/CTPtiming");
1180
1181 if (entry) {
1182
1183 AliInfo("Found an AliCTPTimeParams in GRP/CTP/CTPtiming, reading it");
1184 fCTPTimeParams = dynamic_cast<AliCTPTimeParams*> (entry->GetObject());
1185 entry->SetOwner(0);
1186 return kTRUE;
1187 }
1188
1189 return kFALSE;
1190}
1191//_____________________________________________________________________________
21a3aa09 1192Bool_t AliReconstruction::Run(const char* input)
596a855f 1193{
21a3aa09 1194 // Run Run Run
52dd4a8c 1195 AliCodeTimerAuto("",0);
596a855f 1196
4b71572b 1197 InitRun(input);
1198 if (GetAbort() != TSelector::kContinue) return kFALSE;
1199
1200 TChain *chain = NULL;
1201 if (fRawReader && (chain = fRawReader->GetChain())) {
2a284cef 1202 Long64_t nEntries = (fLastEvent < 0) ? (TChain::kBigNumber) : (fLastEvent - fFirstEvent + 1);
4b71572b 1203 // Proof mode
1204 if (gProof) {
52dd4a8c 1205 // Temporary fix for long raw-data runs (until socket timeout handling in PROOF is revised)
1206 gProof->Exec("gEnv->SetValue(\"Proof.SocketActivityTimeout\",-1)", kTRUE);
6c6f2624 1207
1208 if (gGrid)
1209 gProof->Exec("TGrid::Connect(\"alien://\")",kTRUE);
1210
1211 TMessage::EnableSchemaEvolutionForAll(kTRUE);
1212 gProof->Exec("TMessage::EnableSchemaEvolutionForAll(kTRUE)",kTRUE);
1213
4b71572b 1214 gProof->AddInput(this);
81d57268 1215
1216 if (!ParseOutput()) return kFALSE;
1217
975320a0 1218 gProof->SetParameter("PROOF_MaxSlavesPerNode", 9999);
4b71572b 1219 chain->SetProof();
2a284cef 1220 chain->Process("AliReconstruction","",nEntries,fFirstEvent);
4b71572b 1221 }
1222 else {
2a284cef 1223 chain->Process(this,"",nEntries,fFirstEvent);
4b71572b 1224 }
1225 }
1226 else {
1227 Begin(NULL);
1228 if (GetAbort() != TSelector::kContinue) return kFALSE;
1229 SlaveBegin(NULL);
1230 if (GetAbort() != TSelector::kContinue) return kFALSE;
1231 //******* The loop over events
44e45fac 1232 AliInfo("Starting looping over events");
4b71572b 1233 Int_t iEvent = 0;
1234 while ((iEvent < fRunLoader->GetNumberOfEvents()) ||
1235 (fRawReader && fRawReader->NextEvent())) {
1236 if (!ProcessEvent(iEvent)) {
930e6e3e 1237 Abort("ProcessEvent",TSelector::kAbortFile);
1238 return kFALSE;
4b71572b 1239 }
1240 iEvent++;
1241 }
1242 SlaveTerminate();
1243 if (GetAbort() != TSelector::kContinue) return kFALSE;
1244 Terminate();
1245 if (GetAbort() != TSelector::kContinue) return kFALSE;
21a3aa09 1246 }
21a3aa09 1247
1248 return kTRUE;
1249}
1250
1251//_____________________________________________________________________________
4b71572b 1252void AliReconstruction::InitRawReader(const char* input)
21a3aa09 1253{
52dd4a8c 1254 AliCodeTimerAuto("",0);
21a3aa09 1255
4b71572b 1256 // Init raw-reader and
21a3aa09 1257 // set the input in case of raw data
4b71572b 1258 if (input) fRawInput = input;
1259 fRawReader = AliRawReader::Create(fRawInput.Data());
c35f192e 1260 if (!fRawReader) {
1261 if (fRawInput.IsNull()) {
1262 AliInfo("Reconstruction will run over digits");
1263 }
1264 else {
1265 AliFatal("Can not create raw-data reader ! Exiting...");
1266 }
1267 }
e0027792 1268
35042093 1269 if (!fEquipIdMap.IsNull() && fRawReader)
1270 fRawReader->LoadEquipmentIdsMap(fEquipIdMap);
1271
cd0b062e 1272 if (!fUseHLTData.IsNull()) {
1273 // create the RawReaderHLT which performs redirection of HLT input data for
1274 // the specified detectors
1275 AliRawReader* pRawReader=AliRawHLTManager::CreateRawReaderHLT(fRawReader, fUseHLTData.Data());
1276 if (pRawReader) {
1277 fParentRawReader=fRawReader;
1278 fRawReader=pRawReader;
1279 } else {
1280 AliError(Form("can not create Raw Reader for HLT input %s", fUseHLTData.Data()));
1281 }
1282 }
4b71572b 1283 AliSysInfo::AddStamp("CreateRawReader");
1284}
cd0b062e 1285
4b71572b 1286//_____________________________________________________________________________
1287void AliReconstruction::InitRun(const char* input)
1288{
1289 // Initialization of raw-reader,
1290 // run number, CDB etc.
52dd4a8c 1291 AliCodeTimerAuto("",0);
4b71572b 1292 AliSysInfo::AddStamp("Start");
1293
002c9d1b 1294 // Initialize raw-reader if any
4b71572b 1295 InitRawReader(input);
596a855f 1296
ec92bee0 1297 // Initialize the CDB storage
795e4a22 1298 InitCDB();
ec92bee0 1299
6bae477a 1300 // Set run number in CDBManager (if it is not already set by the user)
4b71572b 1301 if (!SetRunNumberFromData()) {
1302 Abort("SetRunNumberFromData", TSelector::kAbortProcess);
1303 return;
1304 }
002c9d1b 1305
795e4a22 1306 // Set CDB lock: from now on it is forbidden to reset the run number
1307 // or the default storage or to activate any further storage!
1308 SetCDBLock();
1309
002c9d1b 1310}
1311
1312//_____________________________________________________________________________
1313void AliReconstruction::Begin(TTree *)
1314{
1315 // Initialize AlReconstruction before
1316 // going into the event loop
1317 // Should follow the TSelector convention
1318 // i.e. initialize only the object on the client side
52dd4a8c 1319 AliCodeTimerAuto("",0);
db4aeca1 1320
b4fef0d2 1321 AliReconstruction *reco = NULL;
002c9d1b 1322 if (fInput) {
8fb1f1dc 1323 if ((reco = (AliReconstruction*)fInput->FindObject("AliReconstruction"))) {
002c9d1b 1324 *this = *reco;
002c9d1b 1325 }
1326 AliSysInfo::AddStamp("ReadInputInBegin");
1327 }
1328
6bae477a 1329 // Import ideal TGeo geometry and apply misalignment
1330 if (!gGeoManager) {
1331 TString geom(gSystem->DirName(fGAliceFileName));
1332 geom += "/geometry.root";
98e303d9 1333 AliGeomManager::LoadGeometry(geom.Data());
4b71572b 1334 if (!gGeoManager) {
1335 Abort("LoadGeometry", TSelector::kAbortProcess);
1336 return;
1337 }
281270d3 1338 AliSysInfo::AddStamp("LoadGeom");
ff5970a3 1339 TString detsToCheck=fRunLocalReconstruction;
4b71572b 1340 if(!AliGeomManager::CheckSymNamesLUT(detsToCheck.Data())) {
1341 Abort("CheckSymNamesLUT", TSelector::kAbortProcess);
1342 return;
1343 }
281270d3 1344 AliSysInfo::AddStamp("CheckGeom");
6bae477a 1345 }
8e245d15 1346
4b71572b 1347 if (!MisalignGeometry(fLoadAlignData)) {
1348 Abort("MisalignGeometry", TSelector::kAbortProcess);
1349 return;
1350 }
1351 AliCDBManager::Instance()->UnloadFromCache("GRP/Geometry/Data");
281270d3 1352 AliSysInfo::AddStamp("MisalignGeom");
6bae477a 1353
8cb26cdf 1354 if (!InitGRP()) {
1355 Abort("InitGRP", TSelector::kAbortProcess);
1356 return;
1357 }
1358 AliSysInfo::AddStamp("InitGRP");
1359
4b71572b 1360 if (!LoadCDB()) {
1361 Abort("LoadCDB", TSelector::kAbortProcess);
1362 return;
1363 }
1364 AliSysInfo::AddStamp("LoadCDB");
325aa76f 1365
a6dd87ad 1366 if (!LoadTriggerScalersCDB()) {
1367 Abort("LoadTriggerScalersCDB", TSelector::kAbortProcess);
1368 return;
1369 }
1370 AliSysInfo::AddStamp("LoadTriggerScalersCDB");
1371
52dd4a8c 1372 if (!LoadCTPTimeParamsCDB()) {
1373 Abort("LoadCTPTimeParamsCDB", TSelector::kAbortProcess);
1374 return;
1375 }
1376 AliSysInfo::AddStamp("LoadCTPTimeParamsCDB");
a6dd87ad 1377
7e88424f 1378 // Read the reconstruction parameters from OCDB
1379 if (!InitRecoParams()) {
002c9d1b 1380 AliWarning("Not all detectors have correct RecoParam objects initialized");
4b71572b 1381 }
1382 AliSysInfo::AddStamp("InitRecoParams");
1383
fd2b273c 1384 if (fInput && gProof) {
b4fef0d2 1385 if (reco) *reco = *this;
fd2b273c 1386
45a6dee6 1387 gGeoManager->SetName("Geometry");
fd2b273c 1388 gProof->AddInputData(gGeoManager,kTRUE);
002c9d1b 1389 gGeoManager = NULL;
fd2b273c 1390 gProof->AddInputData(const_cast<TMap*>(AliCDBManager::Instance()->GetEntryCache()),kTRUE);
002c9d1b 1391 fInput->Add(new TParameter<Int_t>("RunNumber",AliCDBManager::Instance()->GetRun()));
f7a1cc68 1392 AliMagF *magFieldMap = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
ead7b391 1393 magFieldMap->SetName("MagneticFieldMap");
fd2b273c 1394 gProof->AddInputData(magFieldMap,kTRUE);
002c9d1b 1395 }
4b71572b 1396
1397}
1398
1399//_____________________________________________________________________________
1400void AliReconstruction::SlaveBegin(TTree*)
1401{
1402 // Initialization related to run-loader,
1403 // vertexer, trackers, recontructors
1404 // In proof mode it is executed on the slave
52dd4a8c 1405 AliCodeTimerAuto("",0);
4b71572b 1406
db4aeca1 1407 TProofOutputFile *outProofFile = NULL;
81d57268 1408 if (fInput) {
1409 if (AliDebugLevel() > 0) fInput->Print();
ef5e52f6 1410 if (AliDebugLevel() > 10) fInput->Dump();
db4aeca1 1411 if (AliReconstruction *reco = (AliReconstruction*)fInput->FindObject("AliReconstruction")) {
1412 *this = *reco;
1413 }
4b71572b 1414 if (TGeoManager *tgeo = (TGeoManager*)fInput->FindObject("Geometry")) {
1415 gGeoManager = tgeo;
1416 AliGeomManager::SetGeometry(tgeo);
7e88424f 1417 }
281270d3 1418 if (TMap *entryCache = (TMap*)fInput->FindObject("CDBEntryCache")) {
1419 Int_t runNumber = -1;
1420 if (TProof::GetParameter(fInput,"RunNumber",runNumber) == 0) {
1421 AliCDBManager *man = AliCDBManager::Instance(entryCache,runNumber);
1422 man->SetCacheFlag(kTRUE);
1423 man->SetLock(kTRUE);
1424 man->Print();
1425 }
4b71572b 1426 }
ead7b391 1427 if (AliMagF *map = (AliMagF*)fInput->FindObject("MagneticFieldMap")) {
ef5e52f6 1428 AliMagF *newMap = new AliMagF(*map);
1429 if (!newMap->LoadParameterization()) {
1430 Abort("AliMagF::LoadParameterization", TSelector::kAbortProcess);
1431 return;
1432 }
1433 TGeoGlobalMagField::Instance()->SetField(newMap);
1434 TGeoGlobalMagField::Instance()->Lock();
4b71572b 1435 }
81d57268 1436 if (TNamed *outputFileName = (TNamed*)fInput->FindObject("PROOF_OUTPUTFILE"))
1437 fProofOutputFileName = outputFileName->GetTitle();
1438 if (TNamed *outputLocation = (TNamed*)fInput->FindObject("PROOF_OUTPUTFILE_LOCATION"))
1439 fProofOutputLocation = outputLocation->GetTitle();
1440 if (fInput->FindObject("PROOF_OUTPUTFILE_DATASET"))
1441 fProofOutputDataset = kTRUE;
1442 if (TNamed *archiveList = (TNamed*)fInput->FindObject("PROOF_OUTPUTFILE_ARCHIVE"))
1443 fProofOutputArchive = archiveList->GetTitle();
1444 if (!fProofOutputFileName.IsNull() &&
1445 !fProofOutputLocation.IsNull() &&
1446 fProofOutputArchive.IsNull()) {
1447 if (!fProofOutputDataset) {
1448 outProofFile = new TProofOutputFile(fProofOutputFileName.Data(),"M");
1449 outProofFile->SetOutputFileName(Form("%s%s",fProofOutputLocation.Data(),fProofOutputFileName.Data()));
1450 }
1451 else {
1452 outProofFile = new TProofOutputFile(fProofOutputFileName.Data(),"DROV",fProofOutputLocation.Data());
1453 }
1454 if (AliDebugLevel() > 0) outProofFile->Dump();
db4aeca1 1455 fOutput->Add(outProofFile);
4b71572b 1456 }
002c9d1b 1457 AliSysInfo::AddStamp("ReadInputInSlaveBegin");
7e88424f 1458 }
325aa76f 1459
4b71572b 1460 // get the run loader
1461 if (!InitRunLoader()) {
1462 Abort("InitRunLoader", TSelector::kAbortProcess);
1463 return;
1464 }
1465 AliSysInfo::AddStamp("LoadLoader");
1466
325aa76f 1467 ftVertexer = new AliVertexerTracks(AliTracker::GetBz());
2257f27e 1468
f08fc9f5 1469 // get trackers
b8cd5251 1470 if (!fRunTracking.IsNull() && !CreateTrackers(fRunTracking)) {
4b71572b 1471 Abort("CreateTrackers", TSelector::kAbortProcess);
1472 return;
1473 }
1474 AliSysInfo::AddStamp("CreateTrackers");
24f7a148 1475
36711aa4 1476 // create the ESD output file and tree
db4aeca1 1477 if (!outProofFile) {
1478 ffile = TFile::Open("AliESDs.root", "RECREATE");
1479 ffile->SetCompressionLevel(2);
1480 if (!ffile->IsOpen()) {
1481 Abort("OpenESDFile", TSelector::kAbortProcess);
1482 return;
1483 }
1484 }
1485 else {
81d57268 1486 AliInfo(Form("Opening output PROOF file: %s/%s",
1487 outProofFile->GetDir(), outProofFile->GetFileName()));
db4aeca1 1488 if (!(ffile = outProofFile->OpenFile("RECREATE"))) {
1489 Abort(Form("Problems opening output PROOF file: %s/%s",
1490 outProofFile->GetDir(), outProofFile->GetFileName()),
1491 TSelector::kAbortProcess);
1492 return;
1493 }
596a855f 1494 }
46698ae4 1495
21a3aa09 1496 ftree = new TTree("esdTree", "Tree with ESD objects");
1497 fesd = new AliESDEvent();
1498 fesd->CreateStdContent();
32ba9c61 1499
1500 fesd->WriteToTree(ftree);
1d99986f 1501 if (fWriteESDfriend) {
32ba9c61 1502 // careful:
1503 // Since we add the branch manually we must
1504 // book and add it after WriteToTree
1505 // otherwise it is created twice,
1506 // once via writetotree and once here.
1507 // The case for AliESDfriend is now
1508 // caught also in AlIESDEvent::WriteToTree but
1509 // be careful when changing the name (AliESDfriend is not
1510 // a TNamed so we had to hardwire it)
21a3aa09 1511 fesdf = new AliESDfriend();
1512 TBranch *br=ftree->Branch("ESDfriend.","AliESDfriend", &fesdf);
46698ae4 1513 br->SetFile("AliESDfriends.root");
21a3aa09 1514 fesd->AddObject(fesdf);
1d99986f 1515 }
f9604a22 1516 ftree->GetUserInfo()->Add(fesd);
1517
1518 fhlttree = new TTree("HLTesdTree", "Tree with HLT ESD objects");
1519 fhltesd = new AliESDEvent();
1520 fhltesd->CreateStdContent();
7491aea6 1521
1522 // read the ESD template from CDB
1523 // HLT is allowed to put non-std content to its ESD, the non-std
1524 // objects need to be created before invocation of WriteToTree in
1525 // order to create all branches. Initialization is done from an
1526 // ESD layout template in CDB
1527 AliCDBManager* man = AliCDBManager::Instance();
1528 AliCDBPath hltESDConfigPath("HLT/ConfigHLT/esdLayout");
1529 AliCDBEntry* hltESDConfig=NULL;
1530 if (man->GetId(hltESDConfigPath)!=NULL &&
1531 (hltESDConfig=man->Get(hltESDConfigPath))!=NULL) {
1532 AliESDEvent* pESDLayout=dynamic_cast<AliESDEvent*>(hltESDConfig->GetObject());
1533 if (pESDLayout) {
1534 // init all internal variables from the list of objects
1535 pESDLayout->GetStdContent();
1536
1537 // copy content and create non-std objects
1538 *fhltesd=*pESDLayout;
1539 fhltesd->Reset();
1540 } else {
1541 AliError(Form("error setting hltEsd layout from %s: invalid object type",
1542 hltESDConfigPath.GetPath().Data()));
1543 }
1544 }
1545
f9604a22 1546 fhltesd->WriteToTree(fhlttree);
1547 fhlttree->GetUserInfo()->Add(fhltesd);
5728d3d5 1548
ac4a7581 1549 ProcInfo_t procInfo;
1550 gSystem->GetProcInfo(&procInfo);
1551 AliInfo(Form("Current memory usage %d %d", procInfo.fMemResident, procInfo.fMemVirtual));
8661738e 1552
325aa76f 1553 //QA
930e6e3e 1554 //Initialize the QA and start of cycle
f1c1204d 1555 if (fRunQA || fRunGlobalQA)
1556 InitQA() ;
aa3c69a9 1557
7167ae53 1558 //Initialize the Plane Efficiency framework
1559 if (fRunPlaneEff && !InitPlaneEff()) {
4b71572b 1560 Abort("InitPlaneEff", TSelector::kAbortProcess);
1561 return;
7167ae53 1562 }
aa3c69a9 1563
14dd053c 1564 if (strcmp(gProgName,"alieve") == 0)
1565 fRunAliEVE = InitAliEVE();
1566
4b71572b 1567 return;
21a3aa09 1568}
1569
1570//_____________________________________________________________________________
4b71572b 1571Bool_t AliReconstruction::Process(Long64_t entry)
1572{
1573 // run the reconstruction over a single entry
1574 // from the chain with raw data
52dd4a8c 1575 AliCodeTimerAuto("",0);
4b71572b 1576
1577 TTree *currTree = fChain->GetTree();
33314186 1578 AliRawVEvent *event = NULL;
4b71572b 1579 currTree->SetBranchAddress("rawevent",&event);
1580 currTree->GetEntry(entry);
1581 fRawReader = new AliRawReaderRoot(event);
1582 fStatus = ProcessEvent(fRunLoader->GetNumberOfEvents());
1583 delete fRawReader;
1584 fRawReader = NULL;
1585 delete event;
1586
1587 return fStatus;
1588}
1589
1590//_____________________________________________________________________________
1591void AliReconstruction::Init(TTree *tree)
1592{
1593 if (tree == 0) {
1594 AliError("The input tree is not found!");
1595 return;
1596 }
1597 fChain = tree;
1598}
1599
1600//_____________________________________________________________________________
1601Bool_t AliReconstruction::ProcessEvent(Int_t iEvent)
21a3aa09 1602{
1603 // run the reconstruction over a single event
1604 // The event loop is steered in Run method
1605
52dd4a8c 1606 AliCodeTimerAuto("",0);
21a3aa09 1607
1608 if (iEvent >= fRunLoader->GetNumberOfEvents()) {
1609 fRunLoader->SetEventNumber(iEvent);
1610 fRunLoader->GetHeader()->Reset(fRawReader->GetRunNumber(),
1611 iEvent, iEvent);
21a3aa09 1612 fRunLoader->TreeE()->Fill();
4b71572b 1613 if (fRawReader && fRawReader->UseAutoSaveESD())
1614 fRunLoader->TreeE()->AutoSave("SaveSelf");
21a3aa09 1615 }
1616
1617 if ((iEvent < fFirstEvent) || ((fLastEvent >= 0) && (iEvent > fLastEvent))) {
21a3aa09 1618 return kTRUE;
1619 }
1620
1621 AliInfo(Form("processing event %d", iEvent));
aa3c69a9 1622
16447f18 1623 fRunLoader->GetEvent(iEvent);
1624
7e88424f 1625 // Fill Event-info object
1626 GetEventInfo();
77ba28ba 1627 fRecoParam.SetEventSpecie(fRunInfo,fEventInfo,fListOfCosmicTriggers);
48f5e52d 1628 AliInfo(Form("Current event specie: %s",fRecoParam.PrintEventSpecie()));
7e88424f 1629
a00021a7 1630 // Set the reco-params
1631 {
1632 TString detStr = fLoadCDB;
ac4a7581 1633 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
a00021a7 1634 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
1635 AliReconstructor *reconstructor = GetReconstructor(iDet);
1636 if (reconstructor && fRecoParam.GetDetRecoParamArray(iDet)) {
57acd2d2 1637 const AliDetectorRecoParam *par = fRecoParam.GetDetRecoParam(iDet);
1638 reconstructor->SetRecoParam(par);
52dd4a8c 1639 reconstructor->SetEventInfo(&fEventInfo);
57acd2d2 1640 if (fRunQA) {
eca4fa66 1641 AliQAManager::QAManager()->SetRecoParam(iDet, par) ;
6252ceeb 1642 AliQAManager::QAManager()->SetEventSpecie(AliRecoParam::Convert(par->GetEventSpecie())) ;
57acd2d2 1643 }
a00021a7 1644 }
1645 }
1646 }
1647
ca13fb87 1648 // QA on single raw
514cb8c7 1649 if (fRunQA && IsInTasks(AliQAv1::kRAWS)) {
eca4fa66 1650 AliQAManager::QAManager()->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
1651 AliQAManager::QAManager()->RunOneEvent(fRawReader) ;
57acd2d2 1652 }
d506c543 1653 // local single event reconstruction
b26c3770 1654 if (!fRunLocalReconstruction.IsNull()) {
d506c543 1655 TString detectors=fRunLocalReconstruction;
a441bf51 1656 // run HLT event reconstruction first
d506c543 1657 // ;-( IsSelected changes the string
1658 if (IsSelected("HLT", detectors) &&
1659 !RunLocalEventReconstruction("HLT")) {
4b71572b 1660 if (fStopOnError) {CleanUp(); return kFALSE;}
a441bf51 1661 }
1662 detectors=fRunLocalReconstruction;
1663 detectors.ReplaceAll("HLT", "");
1664 if (!RunLocalEventReconstruction(detectors)) {
4b71572b 1665 if (fStopOnError) {CleanUp(); return kFALSE;}
b26c3770 1666 }
1667 }
1668
21a3aa09 1669 fesd->SetRunNumber(fRunLoader->GetHeader()->GetRun());
1670 fhltesd->SetRunNumber(fRunLoader->GetHeader()->GetRun());
1671 fesd->SetEventNumberInFile(fRunLoader->GetHeader()->GetEventNrInRun());
1672 fhltesd->SetEventNumberInFile(fRunLoader->GetHeader()->GetEventNrInRun());
46698ae4 1673
d6ee376f 1674 // Set magnetic field from the tracker
21a3aa09 1675 fesd->SetMagneticField(AliTracker::GetBz());
1676 fhltesd->SetMagneticField(AliTracker::GetBz());
33fe5eb1 1677 //
1678 AliMagF* fld = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
1679 if (fld) { // set info needed for field initialization
1680 fesd->SetCurrentL3(fld->GetCurrentSol());
1681 fesd->SetCurrentDip(fld->GetCurrentDip());
1682 fesd->SetBeamEnergy(fld->GetBeamEnergy());
1683 fesd->SetBeamType(fld->GetBeamTypeText());
1684 fesd->SetUniformBMap(fld->IsUniform());
1685 fesd->SetBInfoStored();
1686 //
1687 fhltesd->SetCurrentL3(fld->GetCurrentSol());
1688 fhltesd->SetCurrentDip(fld->GetCurrentDip());
1689 fhltesd->SetBeamEnergy(fld->GetBeamEnergy());
1690 fhltesd->SetBeamType(fld->GetBeamTypeText());
1691 fhltesd->SetUniformBMap(fld->IsUniform());
1692 fhltesd->SetBInfoStored();
1693 }
1694 //
71f6cda4 1695 // Set most probable pt, for B=0 tracking
9257a1bd 1696 // Get the global reco-params. They are atposition 16 inside the array of detectors in fRecoParam
ac4a7581 1697 const AliGRPRecoParam *grpRecoParam = dynamic_cast<const AliGRPRecoParam*>(fRecoParam.GetDetRecoParam(kNDetectors));
71f6cda4 1698 if (grpRecoParam) AliExternalTrackParam::SetMostProbablePt(grpRecoParam->GetMostProbablePt());
46698ae4 1699
2e3550da 1700 // Fill raw-data error log into the ESD
21a3aa09 1701 if (fRawReader) FillRawDataErrorLog(iEvent,fesd);
2e3550da 1702
2257f27e 1703 // vertex finder
1704 if (fRunVertexFinder) {
4b71572b 1705 if (!RunVertexFinder(fesd)) {
1706 if (fStopOnError) {CleanUp(); return kFALSE;}
2257f27e 1707 }
1708 }
1709
58e8dc31 1710 // For Plane Efficiency: run the SPD trackleter
1711 if (fRunPlaneEff && fSPDTrackleter) {
1712 if (!RunSPDTrackleting(fesd)) {
1713 if (fStopOnError) {CleanUp(); return kFALSE;}
1714 }
1715 }
1716
e66fbafb 1717 // Muon tracking
b8cd5251 1718 if (!fRunTracking.IsNull()) {
e66fbafb 1719 if (fRunMuonTracking) {
21a3aa09 1720 if (!RunMuonTracking(fesd)) {
4b71572b 1721 if (fStopOnError) {CleanUp(); return kFALSE;}
24f7a148 1722 }
596a855f 1723 }
1724 }
1725
e66fbafb 1726 // barrel tracking
1727 if (!fRunTracking.IsNull()) {
4b71572b 1728 if (!RunTracking(fesd)) {
1729 if (fStopOnError) {CleanUp(); return kFALSE;}
e66fbafb 1730 }
1731 }
21c573b7 1732
596a855f 1733 // fill ESD
1734 if (!fFillESD.IsNull()) {
d506c543 1735 TString detectors=fFillESD;
f6806ad3 1736 // run HLT first and on hltesd
d506c543 1737 // ;-( IsSelected changes the string
1738 if (IsSelected("HLT", detectors) &&
1739 !FillESD(fhltesd, "HLT")) {
4b71572b 1740 if (fStopOnError) {CleanUp(); return kFALSE;}
f6806ad3 1741 }
1742 detectors=fFillESD;
d5105682 1743 // Temporary fix to avoid problems with HLT that overwrites the offline ESDs
1744 if (detectors.Contains("ALL")) {
1745 detectors="";
ac4a7581 1746 for (Int_t idet=0; idet<kNDetectors; ++idet){
d5105682 1747 detectors += fgkDetectorName[idet];
1748 detectors += " ";
1749 }
1750 }
f6806ad3 1751 detectors.ReplaceAll("HLT", "");
1752 if (!FillESD(fesd, detectors)) {
4b71572b 1753 if (fStopOnError) {CleanUp(); return kFALSE;}
596a855f 1754 }
1755 }
a5fa6165 1756
001397cd 1757 // fill Event header information from the RawEventHeader
21a3aa09 1758 if (fRawReader){FillRawEventHeaderESD(fesd);}
63314086 1759 if (fRawReader){FillRawEventHeaderESD(fhltesd);}
596a855f 1760
1761 // combined PID
21a3aa09 1762 AliESDpid::MakePID(fesd);
596a855f 1763
b647652d 1764 if (fFillTriggerESD) {
4b71572b 1765 if (!FillTriggerESD(fesd)) {
1766 if (fStopOnError) {CleanUp(); return kFALSE;}
b647652d 1767 }
1768 }
f7812afc 1769 // Always fill scalers
1770 if (!FillTriggerScalers(fesd)) {
1771 if (fStopOnError) {CleanUp(); return kFALSE;}
1772 }
1773
b647652d 1774
21a3aa09 1775 ffile->cd();
a6ee503a 1776
3c3709c4 1777 //
67be5c77 1778 // Propagate track to the beam pipe (if not already done by ITS)
3c3709c4 1779 //
21a3aa09 1780 const Int_t ntracks = fesd->GetNumberOfTracks();
3c3709c4 1781 const Double_t kRadius = 2.8; //something less than the beam pipe radius
1782
1783 TObjArray trkArray;
1784 UShort_t *selectedIdx=new UShort_t[ntracks];
1785
1786 for (Int_t itrack=0; itrack<ntracks; itrack++){
98a50ff3 1787 const Double_t kMaxStep = 1; //max step over the material
3c3709c4 1788 Bool_t ok;
1789
21a3aa09 1790 AliESDtrack *track = fesd->GetTrack(itrack);
3c3709c4 1791 if (!track) continue;
1792
1793 AliExternalTrackParam *tpcTrack =
1794 (AliExternalTrackParam *)track->GetTPCInnerParam();
bcabd6af 1795 ok = kFALSE;
1796 if (tpcTrack)
1797 ok = AliTracker::
52dd4a8c 1798 PropagateTrackToBxByBz(tpcTrack,kRadius,track->GetMass(),kMaxStep,kFALSE);
43c9dae1 1799
3c3709c4 1800 if (ok) {
1801 Int_t n=trkArray.GetEntriesFast();
1802 selectedIdx[n]=track->GetID();
1803 trkArray.AddLast(tpcTrack);
1804 }
1805
3d65e645 1806 //Tracks refitted by ITS should already be at the SPD vertex
1807 if (track->IsOn(AliESDtrack::kITSrefit)) continue;
1808
1809 AliTracker::
52dd4a8c 1810 PropagateTrackToBxByBz(track,kRadius,track->GetMass(),kMaxStep,kFALSE);
1811 Double_t x[3]; track->GetXYZ(x);
1812 Double_t b[3]; AliTracker::GetBxByBz(x,b);
1813 track->RelateToVertexBxByBz(fesd->GetPrimaryVertexSPD(), b, kVeryBig);
3c3709c4 1814
3c3709c4 1815 }
1816
1817 //
1818 // Improve the reconstructed primary vertex position using the tracks
1819 //
59224b2b 1820 Bool_t runVertexFinderTracks = fRunVertexFinderTracks;
1821 if(fesd->GetPrimaryVertexSPD()) {
1822 TString vtitle = fesd->GetPrimaryVertexSPD()->GetTitle();
1823 if(vtitle.Contains("cosmics")) {
1824 runVertexFinderTracks=kFALSE;
1825 }
c060d7fe 1826 }
a00021a7 1827
1828 if (runVertexFinderTracks) {
3c3709c4 1829 // TPC + ITS primary vertex
f09c879d 1830 ftVertexer->SetITSMode();
f2a195c1 1831 ftVertexer->SetConstraintOff();
a00021a7 1832 // get cuts for vertexer from AliGRPRecoParam
1833 if (grpRecoParam) {
1834 Int_t nCutsVertexer = grpRecoParam->GetVertexerTracksNCuts();
1835 Double_t *cutsVertexer = new Double_t[nCutsVertexer];
1836 grpRecoParam->GetVertexerTracksCutsITS(cutsVertexer);
1837 ftVertexer->SetCuts(cutsVertexer);
1838 delete [] cutsVertexer; cutsVertexer = NULL;
f2a195c1 1839 if(fDiamondProfile && grpRecoParam->GetVertexerTracksConstraintITS())
1840 ftVertexer->SetVtxStart(fDiamondProfile);
43c9dae1 1841 }
21a3aa09 1842 AliESDVertex *pvtx=ftVertexer->FindPrimaryVertex(fesd);
3c3709c4 1843 if (pvtx) {
1844 if (pvtx->GetStatus()) {
c264b61b 1845 fesd->SetPrimaryVertexTracks(pvtx);
3c3709c4 1846 for (Int_t i=0; i<ntracks; i++) {
21a3aa09 1847 AliESDtrack *t = fesd->GetTrack(i);
52dd4a8c 1848 Double_t x[3]; t->GetXYZ(x);
1849 Double_t b[3]; AliTracker::GetBxByBz(x,b);
1850 t->RelateToVertexBxByBz(pvtx, b, kVeryBig);
3c3709c4 1851 }
1852 }
1853 }
1854
1855 // TPC-only primary vertex
f09c879d 1856 ftVertexer->SetTPCMode();
f2a195c1 1857 ftVertexer->SetConstraintOff();
a00021a7 1858 // get cuts for vertexer from AliGRPRecoParam
1859 if (grpRecoParam) {
1860 Int_t nCutsVertexer = grpRecoParam->GetVertexerTracksNCuts();
1861 Double_t *cutsVertexer = new Double_t[nCutsVertexer];
1862 grpRecoParam->GetVertexerTracksCutsTPC(cutsVertexer);
1863 ftVertexer->SetCuts(cutsVertexer);
1864 delete [] cutsVertexer; cutsVertexer = NULL;
f2a195c1 1865 if(fDiamondProfileTPC && grpRecoParam->GetVertexerTracksConstraintTPC())
1866 ftVertexer->SetVtxStart(fDiamondProfileTPC);
43c9dae1 1867 }
21a3aa09 1868 pvtx=ftVertexer->FindPrimaryVertex(&trkArray,selectedIdx);
3c3709c4 1869 if (pvtx) {
1870 if (pvtx->GetStatus()) {
21a3aa09 1871 fesd->SetPrimaryVertexTPC(pvtx);
3d65e645 1872 for (Int_t i=0; i<ntracks; i++) {
1873 AliESDtrack *t = fesd->GetTrack(i);
52dd4a8c 1874 Double_t x[3]; t->GetXYZ(x);
1875 Double_t b[3]; AliTracker::GetBxByBz(x,b);
1876 t->RelateToVertexTPCBxByBz(pvtx, b, kVeryBig);
3c3709c4 1877 }
1878 }
1879 }
1880
1881 }
1882 delete[] selectedIdx;
1883
21a3aa09 1884 if(fDiamondProfile) fesd->SetDiamond(fDiamondProfile);
17c86e90 1885
c5e3e5d1 1886
d1683eef 1887 if (fRunV0Finder) {
1888 // V0 finding
1889 AliV0vertexer vtxer;
21a3aa09 1890 vtxer.Tracks2V0vertices(fesd);
5e4ff34d 1891
d1683eef 1892 if (fRunCascadeFinder) {
1893 // Cascade finding
1894 AliCascadeVertexer cvtxer;
21a3aa09 1895 cvtxer.V0sTracks2CascadeVertices(fesd);
d1683eef 1896 }
5e4ff34d 1897 }
1898
596a855f 1899 // write ESD
21a3aa09 1900 if (fCleanESD) CleanESD(fesd);
854c6476 1901
514cb8c7 1902 if (fRunQA && IsInTasks(AliQAv1::kESDS)) {
eca4fa66 1903 AliQAManager::QAManager()->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
1904 AliQAManager::QAManager()->RunOneEvent(fesd) ;
57acd2d2 1905 }
bea94759 1906 if (fRunGlobalQA) {
eca4fa66 1907 AliQADataMaker *qadm = AliQAManager::QAManager()->GetQADataMaker(AliQAv1::kGLOBAL);
57acd2d2 1908 qadm->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
514cb8c7 1909 if (qadm && IsInTasks(AliQAv1::kESDS))
4e25ac79 1910 qadm->Exec(AliQAv1::kESDS, fesd);
b03591ab 1911 }
854c6476 1912
1d99986f 1913 if (fWriteESDfriend) {
b090e6a3 1914 // fesdf->~AliESDfriend();
1915 // new (fesdf) AliESDfriend(); // Reset...
21a3aa09 1916 fesd->GetESDfriend(fesdf);
1d99986f 1917 }
21a3aa09 1918 ftree->Fill();
500d54ab 1919
2d91a353 1920 // Auto-save the ESD tree in case of prompt reco @P2
be103ac8 1921 if (fRawReader && fRawReader->UseAutoSaveESD()) {
2d91a353 1922 ftree->AutoSave("SaveSelf");
be103ac8 1923 TFile *friendfile = (TFile *)(gROOT->GetListOfFiles()->FindObject("AliESDfriends.root"));
1924 if (friendfile) friendfile->Save();
1925 }
2d91a353 1926
500d54ab 1927 // write HLT ESD
21a3aa09 1928 fhlttree->Fill();
1d99986f 1929
14dd053c 1930 // call AliEVE
1931 if (fRunAliEVE) RunAliEVE();
1932
21a3aa09 1933 fesd->Reset();
1934 fhltesd->Reset();
5728d3d5 1935 if (fWriteESDfriend) {
21a3aa09 1936 fesdf->~AliESDfriend();
1937 new (fesdf) AliESDfriend(); // Reset...
5728d3d5 1938 }
a5fa6165 1939
ac4a7581 1940 ProcInfo_t procInfo;
1941 gSystem->GetProcInfo(&procInfo);
1942 AliInfo(Form("Event %d -> Current memory usage %d %d",iEvent, procInfo.fMemResident, procInfo.fMemVirtual));
a5fa6165 1943
ca13fb87 1944 fEventInfo.Reset();
ac4a7581 1945 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
42ec5d3d 1946 if (fReconstructor[iDet]) {
a00021a7 1947 fReconstructor[iDet]->SetRecoParam(NULL);
42ec5d3d 1948 fReconstructor[iDet]->SetEventInfo(NULL);
1949 }
1f26f3e7 1950 if (fTracker[iDet]) fTracker[iDet]->SetEventInfo(NULL);
ca13fb87 1951 }
1952
53c8f690 1953 if (fRunQA || fRunGlobalQA)
eca4fa66 1954 AliQAManager::QAManager()->Increment() ;
53c8f690 1955
ca13fb87 1956 return kTRUE;
21a3aa09 1957}
1958
21a3aa09 1959//_____________________________________________________________________________
4b71572b 1960void AliReconstruction::SlaveTerminate()
21a3aa09 1961{
4b71572b 1962 // Finalize the run on the slave side
21a3aa09 1963 // Called after the exit
1964 // from the event loop
52dd4a8c 1965 AliCodeTimerAuto("",0);
21a3aa09 1966
1967 if (fIsNewRunLoader) { // galice.root didn't exist
1968 fRunLoader->WriteHeader("OVERWRITE");
1969 fRunLoader->CdGAFile();
1970 fRunLoader->Write(0, TObject::kOverwrite);
1971 }
1972
f747912b 1973 const TMap *cdbMap = AliCDBManager::Instance()->GetStorageMap();
1974 const TList *cdbList = AliCDBManager::Instance()->GetRetrievedIds();
1975
1976 TMap *cdbMapCopy = new TMap(cdbMap->GetEntries());
1977 cdbMapCopy->SetOwner(1);
1978 cdbMapCopy->SetName("cdbMap");
1979 TIter iter(cdbMap->GetTable());
1980
1981 TPair* pair = 0;
1982 while((pair = dynamic_cast<TPair*> (iter.Next()))){
1983 TObjString* keyStr = dynamic_cast<TObjString*> (pair->Key());
1984 TObjString* valStr = dynamic_cast<TObjString*> (pair->Value());
1985 cdbMapCopy->Add(new TObjString(keyStr->GetName()), new TObjString(valStr->GetName()));
1986 }
1987
1988 TList *cdbListCopy = new TList();
1989 cdbListCopy->SetOwner(1);
1990 cdbListCopy->SetName("cdbList");
1991
1992 TIter iter2(cdbList);
1993
b940cb9b 1994 AliCDBId* id=0;
e84c88f5 1995 while((id = dynamic_cast<AliCDBId*> (iter2.Next()))){
a4970db9 1996 cdbListCopy->Add(new TObjString(id->ToString().Data()));
f747912b 1997 }
1998
21a3aa09 1999 ftree->GetUserInfo()->Add(cdbMapCopy);
2000 ftree->GetUserInfo()->Add(cdbListCopy);
abe0c04e 2001
46698ae4 2002
21a3aa09 2003 ffile->cd();
aa3c69a9 2004
a9c0e6db 2005 if (fWriteESDfriend)
21a3aa09 2006 ftree->SetBranchStatus("ESDfriend*",0);
562dd0b4 2007 // we want to have only one tree version number
21a3aa09 2008 ftree->Write(ftree->GetName(),TObject::kOverwrite);
63314086 2009 fhlttree->Write(fhlttree->GetName(),TObject::kOverwrite);
f3a97c86 2010
a7a1e1c7 2011// Finish with Plane Efficiency evaluation: before of CleanUp !!!
2012 if (fRunPlaneEff && !FinishPlaneEff()) {
2013 AliWarning("Finish PlaneEff evaluation failed");
2014 }
2015
930e6e3e 2016 // End of cycle for the in-loop
57acd2d2 2017 if (fRunQA)
eca4fa66 2018 AliQAManager::QAManager()->EndOfCycle() ;
57acd2d2 2019
930e6e3e 2020 if (fRunGlobalQA) {
eca4fa66 2021 AliQADataMaker *qadm = AliQAManager::QAManager()->GetQADataMaker(AliQAv1::kGLOBAL);
4b71572b 2022 if (qadm) {
514cb8c7 2023 if (IsInTasks(AliQAv1::kRECPOINTS))
4e25ac79 2024 qadm->EndOfCycle(AliQAv1::kRECPOINTS);
514cb8c7 2025 if (IsInTasks(AliQAv1::kESDS))
4e25ac79 2026 qadm->EndOfCycle(AliQAv1::kESDS);
4b71572b 2027 qadm->Finish();
2028 }
a7807689 2029 }
87da0921 2030
2031 if (fRunQA || fRunGlobalQA) {
81d57268 2032 if (fInput &&
2033 !fProofOutputLocation.IsNull() &&
2034 fProofOutputArchive.IsNull() &&
2035 !fProofOutputDataset) {
2036 TString qaOutputFile(Form("%sMerged.%s.Data.root",
2037 fProofOutputLocation.Data(),
2038 AliQAv1::GetQADataFileName()));
2039 TProofOutputFile *qaProofFile = new TProofOutputFile(Form("Merged.%s.Data.root",
2040 AliQAv1::GetQADataFileName()));
2041 qaProofFile->SetOutputFileName(qaOutputFile.Data());
2042 if (AliDebugLevel() > 0) qaProofFile->Dump();
2043 fOutput->Add(qaProofFile);
2044 MergeQA(qaProofFile->GetFileName());
87da0921 2045 }
2046 else {
2047 MergeQA();
2048 }
2049 }
2050
4b71572b 2051 gROOT->cd();
2052 CleanUp();
81d57268 2053
2054 if (fInput) {
2055 if (!fProofOutputFileName.IsNull() &&
2056 !fProofOutputLocation.IsNull() &&
2057 fProofOutputDataset &&
2058 !fProofOutputArchive.IsNull()) {
2059 TProofOutputFile *zipProofFile = new TProofOutputFile(fProofOutputFileName.Data(),
2060 "DROV",
2061 fProofOutputLocation.Data());
2062 if (AliDebugLevel() > 0) zipProofFile->Dump();
2063 fOutput->Add(zipProofFile);
2064 TString fileList(fProofOutputArchive.Data());
2065 fileList.ReplaceAll(","," ");
38c18bf1 2066 TString command;
2067#if ROOT_SVN_REVISION >= 30174
2068 command.Form("zip -n root %s/%s %s",zipProofFile->GetDir(kTRUE),zipProofFile->GetFileName(),fileList.Data());
2069#else
2070 command.Form("zip -n root %s/%s %s",zipProofFile->GetDir(),zipProofFile->GetFileName(),fileList.Data());
2071#endif
2072 AliInfo(Form("Executing: %s",command.Data()));
2073 gSystem->Exec(command.Data());
81d57268 2074 }
2075 }
4b71572b 2076}
2077
2078//_____________________________________________________________________________
2079void AliReconstruction::Terminate()
2080{
f3a97c86 2081 // Create tags for the events in the ESD tree (the ESD tree is always present)
2082 // In case of empty events the tags will contain dummy values
52dd4a8c 2083 AliCodeTimerAuto("",0);
4b71572b 2084
e6d66370 2085 // Do not call the ESD tag creator in case of PROOF-based reconstruction
2086 if (!fInput) {
2087 AliESDTagCreator *esdtagCreator = new AliESDTagCreator();
2088 esdtagCreator->CreateESDTags(fFirstEvent,fLastEvent,fGRPData, AliQAv1::Instance()->GetQA(), AliQAv1::Instance()->GetEventSpecies(), AliQAv1::kNDET, AliRecoParam::kNSpecies);
2089 }
e84c88f5 2090
795e4a22 2091 // Cleanup of CDB manager: cache and active storages!
2092 AliCDBManager::Instance()->ClearCache();
596a855f 2093}
2094
b26c3770 2095//_____________________________________________________________________________
2096Bool_t AliReconstruction::RunLocalEventReconstruction(const TString& detectors)
2097{
2098// run the local reconstruction
aa3c69a9 2099
0f88822a 2100 static Int_t eventNr=0;
52dd4a8c 2101 AliCodeTimerAuto("",0)
b26c3770 2102
2103 TString detStr = detectors;
ac4a7581 2104 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b26c3770 2105 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
2106 AliReconstructor* reconstructor = GetReconstructor(iDet);
2107 if (!reconstructor) continue;
2108 AliLoader* loader = fLoader[iDet];
f6806ad3 2109 // Matthias April 2008: temporary fix to run HLT reconstruction
2110 // although the HLT loader is missing
2111 if (strcmp(fgkDetectorName[iDet], "HLT")==0) {
2112 if (fRawReader) {
44ed7a66 2113 reconstructor->Reconstruct(fRawReader, NULL);
f6806ad3 2114 } else {
44ed7a66 2115 TTree* dummy=NULL;
2116 reconstructor->Reconstruct(dummy, NULL);
f6806ad3 2117 }
2118 continue;
2119 }
d76c31f4 2120 if (!loader) {
2121 AliWarning(Form("No loader is defined for %s!",fgkDetectorName[iDet]));
2122 continue;
2123 }
b26c3770 2124 // conversion of digits
2125 if (fRawReader && reconstructor->HasDigitConversion()) {
2126 AliInfo(Form("converting raw data digits into root objects for %s",
2127 fgkDetectorName[iDet]));
30bbd491 2128// AliCodeTimerAuto(Form("converting raw data digits into root objects for %s",
52dd4a8c 2129// fgkDetectorName[iDet]),0);
b26c3770 2130 loader->LoadDigits("update");
2131 loader->CleanDigits();
2132 loader->MakeDigitsContainer();
2133 TTree* digitsTree = loader->TreeD();
2134 reconstructor->ConvertDigits(fRawReader, digitsTree);
2135 loader->WriteDigits("OVERWRITE");
2136 loader->UnloadDigits();
b26c3770 2137 }
b26c3770 2138 // local reconstruction
b26c3770 2139 AliInfo(Form("running reconstruction for %s", fgkDetectorName[iDet]));
52dd4a8c 2140 //AliCodeTimerAuto(Form("running reconstruction for %s", fgkDetectorName[iDet]),0);
b26c3770 2141 loader->LoadRecPoints("update");
2142 loader->CleanRecPoints();
2143 loader->MakeRecPointsContainer();
2144 TTree* clustersTree = loader->TreeR();
2145 if (fRawReader && !reconstructor->HasDigitConversion()) {
2146 reconstructor->Reconstruct(fRawReader, clustersTree);
2147 } else {
2148 loader->LoadDigits("read");
2149 TTree* digitsTree = loader->TreeD();
2150 if (!digitsTree) {
44ed7a66 2151 AliError(Form("Can't get the %s digits tree", fgkDetectorName[iDet]));
2152 if (fStopOnError) return kFALSE;
b26c3770 2153 } else {
44ed7a66 2154 reconstructor->Reconstruct(digitsTree, clustersTree);
514cb8c7 2155 if (fRunQA && IsInTasks(AliQAv1::kDIGITSR)) {
eca4fa66 2156 AliQAManager::QAManager()->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
2157 AliQAManager::QAManager()->RunOneEventInOneDetector(iDet, digitsTree) ;
44ed7a66 2158 }
b26c3770 2159 }
2160 loader->UnloadDigits();
2161 }
514cb8c7 2162 if (fRunQA && IsInTasks(AliQAv1::kRECPOINTS)) {
eca4fa66 2163 AliQAManager::QAManager()->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
2164 AliQAManager::QAManager()->RunOneEventInOneDetector(iDet, clustersTree) ;
57acd2d2 2165 }
eca4fa66 2166 loader->WriteRecPoints("OVERWRITE");
2167 loader->UnloadRecPoints();
2168 AliSysInfo::AddStamp(Form("LRec%s_%d",fgkDetectorName[iDet],eventNr), iDet,1,eventNr);
a00021a7 2169 }
2170 if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
b26c3770 2171 AliError(Form("the following detectors were not found: %s",
2172 detStr.Data()));
2173 if (fStopOnError) return kFALSE;
2174 }
0f88822a 2175 eventNr++;
b26c3770 2176 return kTRUE;
2177}
58e8dc31 2178//_____________________________________________________________________________
2179Bool_t AliReconstruction::RunSPDTrackleting(AliESDEvent*& esd)
2180{
2181// run the SPD trackleting (for SPD efficiency purpouses)
2182
52dd4a8c 2183 AliCodeTimerAuto("",0)
58e8dc31 2184
2185 Double_t vtxPos[3] = {0, 0, 0};
2186 Double_t vtxErr[3] = {0.0, 0.0, 0.0};
2187/*
2188 TArrayF mcVertex(3);
2189 // if(MC)
2190 if (fRunLoader->GetHeader() && fRunLoader->GetHeader()->GenEventHeader()) {
2191 fRunLoader->GetHeader()->GenEventHeader()->PrimaryVertex(mcVertex);
2192 for (Int_t i = 0; i < 3; i++) vtxPos[i] = mcVertex[i];
2193 }
2194*/
2195 const AliESDVertex *vertex = esd->GetVertex();
2196 if(!vertex){
2197 AliWarning("Vertex not found");
2198 return kFALSE;
2199 }
2200 vertex->GetXYZ(vtxPos);
2201 vertex->GetSigmaXYZ(vtxErr);
2202 if (fSPDTrackleter) {
2203 AliInfo("running the SPD Trackleter for Plane Efficiency Evaluation");
2204
2205 // load clusters
2206 fLoader[0]->LoadRecPoints("read");
2207 TTree* tree = fLoader[0]->TreeR();
2208 if (!tree) {
2209 AliError("Can't get the ITS cluster tree");
2210 return kFALSE;
2211 }
2212 fSPDTrackleter->LoadClusters(tree);
2213 fSPDTrackleter->SetVertex(vtxPos, vtxErr);
2214 // run trackleting
2215 if (fSPDTrackleter->Clusters2Tracks(esd) != 0) {
2216 AliError("AliITSTrackleterSPDEff Clusters2Tracks failed");
2217 // fLoader[0]->UnloadRecPoints();
2218 return kFALSE;
2219 }
2220//fSPDTrackleter->UnloadRecPoints();
2221 } else {
2222 AliWarning("SPDTrackleter not available");
2223 return kFALSE;
2224 }
2225 return kTRUE;
2226}
b26c3770 2227
596a855f 2228//_____________________________________________________________________________
af885e0f 2229Bool_t AliReconstruction::RunVertexFinder(AliESDEvent*& esd)
596a855f 2230{
2231// run the barrel tracking
2232
52dd4a8c 2233 AliCodeTimerAuto("",0)
030b532d 2234
92bffc4d 2235 AliVertexer *vertexer = CreateVertexer();
2236 if (!vertexer) return kFALSE;
2237
2238 AliInfo("running the ITS vertex finder");
2257f27e 2239 AliESDVertex* vertex = NULL;
92bffc4d 2240 if (fLoader[0]) {
2241 fLoader[0]->LoadRecPoints();
2242 TTree* cltree = fLoader[0]->TreeR();
2243 if (cltree) {
2244 if(fDiamondProfileSPD) vertexer->SetVtxStart(fDiamondProfileSPD);
2245 vertex = vertexer->FindVertexForCurrentEvent(cltree);
308c2f7c 2246 }
2247 else {
92bffc4d 2248 AliError("Can't get the ITS cluster tree");
308c2f7c 2249 }
92bffc4d 2250 fLoader[0]->UnloadRecPoints();
2257f27e 2251 }
92bffc4d 2252 else {
2253 AliError("Can't get the ITS loader");
2254 }
2255 if(!vertex){
2256 AliWarning("Vertex not found");
2257 vertex = new AliESDVertex();
2258 vertex->SetName("default");
2259 }
2260 else {
2261 vertex->SetName("reconstructed");
2257f27e 2262 }
92bffc4d 2263
2264 Double_t vtxPos[3];
2265 Double_t vtxErr[3];
2266 vertex->GetXYZ(vtxPos);
2267 vertex->GetSigmaXYZ(vtxErr);
2268
06cc9d95 2269 esd->SetPrimaryVertexSPD(vertex);
73c51de2 2270 AliESDVertex *vpileup = NULL;
2271 Int_t novertices = 0;
2272 vpileup = vertexer->GetAllVertices(novertices);
2273 if(novertices>1){
2274 for (Int_t kk=1; kk<novertices; kk++)esd->AddPileupVertexSPD(&vpileup[kk]);
2275 }
32e449be 2276 // if SPD multiplicity has been determined, it is stored in the ESD
92bffc4d 2277 AliMultiplicity *mult = vertexer->GetMultiplicity();
32e449be 2278 if(mult)esd->SetMultiplicity(mult);
2279
ac4a7581 2280 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 2281 if (fTracker[iDet]) fTracker[iDet]->SetVertex(vtxPos, vtxErr);
2282 }
2257f27e 2283 delete vertex;
2284
92bffc4d 2285 delete vertexer;
2286
2257f27e 2287 return kTRUE;
2288}
2289
1f46a9ae 2290//_____________________________________________________________________________
af885e0f 2291Bool_t AliReconstruction::RunHLTTracking(AliESDEvent*& esd)
1f46a9ae 2292{
2293// run the HLT barrel tracking
2294
52dd4a8c 2295 AliCodeTimerAuto("",0)
1f46a9ae 2296
2297 if (!fRunLoader) {
2298 AliError("Missing runLoader!");
2299 return kFALSE;
2300 }
2301
2302 AliInfo("running HLT tracking");
2303
2304 // Get a pointer to the HLT reconstructor
ac4a7581 2305 AliReconstructor *reconstructor = GetReconstructor(kNDetectors-1);
1f46a9ae 2306 if (!reconstructor) return kFALSE;
2307
2308 // TPC + ITS
2309 for (Int_t iDet = 1; iDet >= 0; iDet--) {
2310 TString detName = fgkDetectorName[iDet];
2311 AliDebug(1, Form("%s HLT tracking", detName.Data()));
2312 reconstructor->SetOption(detName.Data());
d76c31f4 2313 AliTracker *tracker = reconstructor->CreateTracker();
1f46a9ae 2314 if (!tracker) {
2315 AliWarning(Form("couldn't create a HLT tracker for %s", detName.Data()));
2316 if (fStopOnError) return kFALSE;
9dcc06e1 2317 continue;
1f46a9ae 2318 }
2319 Double_t vtxPos[3];
2320 Double_t vtxErr[3]={0.005,0.005,0.010};
2321 const AliESDVertex *vertex = esd->GetVertex();
2322 vertex->GetXYZ(vtxPos);
2323 tracker->SetVertex(vtxPos,vtxErr);
2324 if(iDet != 1) {
2325 fLoader[iDet]->LoadRecPoints("read");
2326 TTree* tree = fLoader[iDet]->TreeR();
2327 if (!tree) {
2328 AliError(Form("Can't get the %s cluster tree", detName.Data()));
2329 return kFALSE;
2330 }
2331 tracker->LoadClusters(tree);
2332 }
2333 if (tracker->Clusters2Tracks(esd) != 0) {
2334 AliError(Form("HLT %s Clusters2Tracks failed", fgkDetectorName[iDet]));
2335 return kFALSE;
2336 }
2337 if(iDet != 1) {
2338 tracker->UnloadClusters();
2339 }
2340 delete tracker;
2341 }
2342
1f46a9ae 2343 return kTRUE;
2344}
2345
e66fbafb 2346//_____________________________________________________________________________
af885e0f 2347Bool_t AliReconstruction::RunMuonTracking(AliESDEvent*& esd)
e66fbafb 2348{
2349// run the muon spectrometer tracking
2350
52dd4a8c 2351 AliCodeTimerAuto("",0)
e66fbafb 2352
2353 if (!fRunLoader) {
2354 AliError("Missing runLoader!");
2355 return kFALSE;
2356 }
2357 Int_t iDet = 7; // for MUON
2358
2359 AliInfo("is running...");
2360
2361 // Get a pointer to the MUON reconstructor
2362 AliReconstructor *reconstructor = GetReconstructor(iDet);
2363 if (!reconstructor) return kFALSE;
2364
2365
2366 TString detName = fgkDetectorName[iDet];
2367 AliDebug(1, Form("%s tracking", detName.Data()));
d76c31f4 2368 AliTracker *tracker = reconstructor->CreateTracker();
e66fbafb 2369 if (!tracker) {
2370 AliWarning(Form("couldn't create a tracker for %s", detName.Data()));
2371 return kFALSE;
2372 }
2373
e66fbafb 2374 // read RecPoints
761350a6 2375 fLoader[iDet]->LoadRecPoints("read");
c1954ee5 2376
761350a6 2377 tracker->LoadClusters(fLoader[iDet]->TreeR());
2378
2379 Int_t rv = tracker->Clusters2Tracks(esd);
2380
761350a6 2381 if ( rv )
2382 {
e66fbafb 2383 AliError(Form("%s Clusters2Tracks failed", fgkDetectorName[iDet]));
2384 return kFALSE;
2385 }
761350a6 2386
e66fbafb 2387 fLoader[iDet]->UnloadRecPoints();
2388
c1954ee5 2389 tracker->UnloadClusters();
2390
e66fbafb 2391 delete tracker;
2392
e66fbafb 2393 return kTRUE;
2394}
2395
2396
2257f27e 2397//_____________________________________________________________________________
af885e0f 2398Bool_t AliReconstruction::RunTracking(AliESDEvent*& esd)
2257f27e 2399{
2400// run the barrel tracking
0f88822a 2401 static Int_t eventNr=0;
52dd4a8c 2402 AliCodeTimerAuto("",0)
24f7a148 2403
815c2b38 2404 AliInfo("running tracking");
596a855f 2405
1f26f3e7 2406 // Set the event info which is used
2407 // by the trackers in order to obtain
2408 // information about read-out detectors,
2409 // trigger etc.
2410 AliDebug(1, "Setting event info");
2411 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
2412 if (!fTracker[iDet]) continue;
2413 fTracker[iDet]->SetEventInfo(&fEventInfo);
2414 }
2415
91b876d1 2416 //Fill the ESD with the T0 info (will be used by the TOF)
d76c31f4 2417 if (fReconstructor[11] && fLoader[11]) {
2418 fLoader[11]->LoadRecPoints("READ");
2419 TTree *treeR = fLoader[11]->TreeR();
89916438 2420 if (treeR) {
2421 GetReconstructor(11)->FillESD((TTree *)NULL,treeR,esd);
2422 }
d76c31f4 2423 }
91b876d1 2424
b8cd5251 2425 // pass 1: TPC + ITS inwards
2426 for (Int_t iDet = 1; iDet >= 0; iDet--) {
2427 if (!fTracker[iDet]) continue;
2428 AliDebug(1, Form("%s tracking", fgkDetectorName[iDet]));
24f7a148 2429
b8cd5251 2430 // load clusters
2431 fLoader[iDet]->LoadRecPoints("read");
6efecea1 2432 AliSysInfo::AddStamp(Form("RLoadCluster%s_%d",fgkDetectorName[iDet],eventNr),iDet,1, eventNr);
b8cd5251 2433 TTree* tree = fLoader[iDet]->TreeR();
2434 if (!tree) {
2435 AliError(Form("Can't get the %s cluster tree", fgkDetectorName[iDet]));
24f7a148 2436 return kFALSE;
2437 }
b8cd5251 2438 fTracker[iDet]->LoadClusters(tree);
6efecea1 2439 AliSysInfo::AddStamp(Form("TLoadCluster%s_%d",fgkDetectorName[iDet],eventNr), iDet,2, eventNr);
b8cd5251 2440 // run tracking
2441 if (fTracker[iDet]->Clusters2Tracks(esd) != 0) {
2442 AliError(Form("%s Clusters2Tracks failed", fgkDetectorName[iDet]));
24f7a148 2443 return kFALSE;
2444 }
878e1fe1 2445 // preliminary PID in TPC needed by the ITS tracker
2446 if (iDet == 1) {
b26c3770 2447 GetReconstructor(1)->FillESD((TTree*)NULL, (TTree*)NULL, esd);
878e1fe1 2448 AliESDpid::MakePID(esd);
0f88822a 2449 }
6efecea1 2450 AliSysInfo::AddStamp(Form("Tracking0%s_%d",fgkDetectorName[iDet],eventNr), iDet,3,eventNr);
b8cd5251 2451 }
596a855f 2452
b8cd5251 2453 // pass 2: ALL backwards
aa3c69a9 2454
ac4a7581 2455 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 2456 if (!fTracker[iDet]) continue;
2457 AliDebug(1, Form("%s back propagation", fgkDetectorName[iDet]));
2458
2459 // load clusters
2460 if (iDet > 1) { // all except ITS, TPC
2461 TTree* tree = NULL;
7b61cd9c 2462 fLoader[iDet]->LoadRecPoints("read");
6efecea1 2463 AliSysInfo::AddStamp(Form("RLoadCluster0%s_%d",fgkDetectorName[iDet],eventNr), iDet,1, eventNr);
7b61cd9c 2464 tree = fLoader[iDet]->TreeR();
b8cd5251 2465 if (!tree) {
eca4fa66 2466 AliError(Form("Can't get the %s cluster tree", fgkDetectorName[iDet]));
2467 return kFALSE;
24f7a148 2468 }
0f88822a 2469 fTracker[iDet]->LoadClusters(tree);
6efecea1 2470 AliSysInfo::AddStamp(Form("TLoadCluster0%s_%d",fgkDetectorName[iDet],eventNr), iDet,2, eventNr);
b8cd5251 2471 }
24f7a148 2472
b8cd5251 2473 // run tracking
283f39c6 2474 if (iDet>1) // start filling residuals for the "outer" detectors
eca4fa66 2475 if (fRunGlobalQA) {
2476 AliTracker::SetFillResiduals(fRecoParam.GetEventSpecie(), kTRUE);
2477 TObjArray ** arr = AliTracker::GetResidualsArray() ;
c9526f68 2478 if (arr) {
0a349581 2479 AliRecoParam::EventSpecie_t es=fRecoParam.GetEventSpecie();
2480 TObjArray * elem = arr[AliRecoParam::AConvert(es)];
c9526f68 2481 if ( elem && (! elem->At(0)) ) {
2482 AliQADataMaker *qadm = AliQAManager::QAManager()->GetQADataMaker(AliQAv1::kGLOBAL);
2483 if (qadm) qadm->InitRecPointsForTracker() ;
2484 }
2485 }
eca4fa66 2486 }
b8cd5251 2487 if (fTracker[iDet]->PropagateBack(esd) != 0) {
2488 AliError(Form("%s backward propagation failed", fgkDetectorName[iDet]));
49dfd67a 2489 // return kFALSE;
b8cd5251 2490 }
24f7a148 2491
b8cd5251 2492 // unload clusters
6e65648b 2493 if (iDet > 3) { // all except ITS, TPC, TRD and TOF
b8cd5251 2494 fTracker[iDet]->UnloadClusters();
7b61cd9c 2495 fLoader[iDet]->UnloadRecPoints();
b8cd5251 2496 }
8f37df88 2497 // updated PID in TPC needed by the ITS tracker -MI
2498 if (iDet == 1) {
8f37df88 2499 GetReconstructor(1)->FillESD((TTree*)NULL, (TTree*)NULL, esd);
2500 AliESDpid::MakePID(esd);
2501 }
6efecea1 2502 AliSysInfo::AddStamp(Form("Tracking1%s_%d",fgkDetectorName[iDet],eventNr), iDet,3, eventNr);
b8cd5251 2503 }
283f39c6 2504 //stop filling residuals for the "outer" detectors
57acd2d2 2505 if (fRunGlobalQA) AliTracker::SetFillResiduals(fRecoParam.GetEventSpecie(), kFALSE);
596a855f 2506
b8cd5251 2507 // pass 3: TRD + TPC + ITS refit inwards
aa3c69a9 2508
b8cd5251 2509 for (Int_t iDet = 2; iDet >= 0; iDet--) {
2510 if (!fTracker[iDet]) continue;
2511 AliDebug(1, Form("%s inward refit", fgkDetectorName[iDet]));
596a855f 2512
b8cd5251 2513 // run tracking
283f39c6 2514 if (iDet<2) // start filling residuals for TPC and ITS
eca4fa66 2515 if (fRunGlobalQA) {
2516 AliTracker::SetFillResiduals(fRecoParam.GetEventSpecie(), kTRUE);
2517 TObjArray ** arr = AliTracker::GetResidualsArray() ;
c9526f68 2518 if (arr) {
0a349581 2519 AliRecoParam::EventSpecie_t es=fRecoParam.GetEventSpecie();
2520 TObjArray * elem = arr[AliRecoParam::AConvert(es)];
c9526f68 2521 if ( elem && (! elem->At(0)) ) {
2522 AliQADataMaker *qadm = AliQAManager::QAManager()->GetQADataMaker(AliQAv1::kGLOBAL);
2523 if (qadm) qadm->InitRecPointsForTracker() ;
2524 }
2525 }
eca4fa66 2526 }
2527
b8cd5251 2528 if (fTracker[iDet]->RefitInward(esd) != 0) {
2529 AliError(Form("%s inward refit failed", fgkDetectorName[iDet]));
49dfd67a 2530 // return kFALSE;
b8cd5251 2531 }
db2368d0 2532 // run postprocessing
2533 if (fTracker[iDet]->PostProcess(esd) != 0) {
2534 AliError(Form("%s postprocessing failed", fgkDetectorName[iDet]));
2535 // return kFALSE;
2536 }
6efecea1 2537 AliSysInfo::AddStamp(Form("Tracking2%s_%d",fgkDetectorName[iDet],eventNr), iDet,3, eventNr);
6e65648b 2538 }
2539
2540 // write space-points to the ESD in case alignment data output
2541 // is switched on
2542 if (fWriteAlignmentData)
2543 WriteAlignmentData(esd);
2544
2545 for (Int_t iDet = 3; iDet >= 0; iDet--) {
2546 if (!fTracker[iDet]) continue;
b8cd5251 2547 // unload clusters
2548 fTracker[iDet]->UnloadClusters();
6efecea1 2549 AliSysInfo::AddStamp(Form("TUnloadCluster%s_%d",fgkDetectorName[iDet],eventNr), iDet,4, eventNr);
b8cd5251 2550 fLoader[iDet]->UnloadRecPoints();
6efecea1 2551 AliSysInfo::AddStamp(Form("RUnloadCluster%s_%d",fgkDetectorName[iDet],eventNr), iDet,5, eventNr);
b8cd5251 2552 }
283f39c6 2553 // stop filling residuals for TPC and ITS
57acd2d2 2554 if (fRunGlobalQA) AliTracker::SetFillResiduals(fRecoParam.GetEventSpecie(), kFALSE);
854c6476 2555
0f88822a 2556 eventNr++;
596a855f 2557 return kTRUE;
2558}
2559
d64bd07d 2560//_____________________________________________________________________________
2561Bool_t AliReconstruction::CleanESD(AliESDEvent *esd){
2562 //
2563 // Remove the data which are not needed for the physics analysis.
2564 //
2565
d64bd07d 2566 Int_t nTracks=esd->GetNumberOfTracks();
18571674 2567 Int_t nV0s=esd->GetNumberOfV0s();
cf37fd88 2568 AliInfo
2569 (Form("Number of ESD tracks and V0s before cleaning: %d %d",nTracks,nV0s));
d64bd07d 2570
18571674 2571 Float_t cleanPars[]={fV0DCAmax,fV0CsPmin,fDmax,fZmax};
7f68891d 2572 Bool_t rc=esd->Clean(cleanPars);
d64bd07d 2573
7f68891d 2574 nTracks=esd->GetNumberOfTracks();
18571674 2575 nV0s=esd->GetNumberOfV0s();
cf37fd88 2576 AliInfo
ae5d5566 2577 (Form("Number of ESD tracks and V0s after cleaning %d %d",nTracks,nV0s));
d64bd07d 2578
7f68891d 2579 return rc;
d64bd07d 2580}
2581
596a855f 2582//_____________________________________________________________________________
af885e0f 2583Bool_t AliReconstruction::FillESD(AliESDEvent*& esd, const TString& detectors)
596a855f 2584{
2585// fill the event summary data
2586
52dd4a8c 2587 AliCodeTimerAuto("",0)
0f88822a 2588 static Int_t eventNr=0;
596a855f 2589 TString detStr = detectors;
abe0c04e 2590
ac4a7581 2591 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
abe0c04e 2592 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
b8cd5251 2593 AliReconstructor* reconstructor = GetReconstructor(iDet);
2594 if (!reconstructor) continue;
4b71572b 2595 AliDebug(1, Form("filling ESD for %s", fgkDetectorName[iDet]));
2596 TTree* clustersTree = NULL;
2597 if (fLoader[iDet]) {
2598 fLoader[iDet]->LoadRecPoints("read");
2599 clustersTree = fLoader[iDet]->TreeR();
2600 if (!clustersTree) {
2601 AliError(Form("Can't get the %s clusters tree",
2602 fgkDetectorName[iDet]));
2603 if (fStopOnError) return kFALSE;
2604 }
2605 }
2606 if (fRawReader && !reconstructor->HasDigitConversion()) {
2607 reconstructor->FillESD(fRawReader, clustersTree, esd);
2608 } else {
2609 TTree* digitsTree = NULL;
d76c31f4 2610 if (fLoader[iDet]) {
4b71572b 2611 fLoader[iDet]->LoadDigits("read");
2612 digitsTree = fLoader[iDet]->TreeD();
2613 if (!digitsTree) {
2614 AliError(Form("Can't get the %s digits tree",
b26c3770 2615 fgkDetectorName[iDet]));
2616 if (fStopOnError) return kFALSE;
2617 }
2618 }
4b71572b 2619 reconstructor->FillESD(digitsTree, clustersTree, esd);
2620 if (fLoader[iDet]) fLoader[iDet]->UnloadDigits();
2621 }
2622 if (fLoader[iDet]) {
2623 fLoader[iDet]->UnloadRecPoints();
596a855f 2624 }
2625 }
2626
2627 if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
815c2b38 2628 AliError(Form("the following detectors were not found: %s",
2629 detStr.Data()));
596a855f 2630 if (fStopOnError) return kFALSE;
2631 }
6efecea1 2632 AliSysInfo::AddStamp(Form("FillESD%d",eventNr), 0,1, eventNr);
0f88822a 2633 eventNr++;
596a855f 2634 return kTRUE;
2635}
2636
b647652d 2637//_____________________________________________________________________________
af885e0f 2638Bool_t AliReconstruction::FillTriggerESD(AliESDEvent*& esd)
b647652d 2639{
2640 // Reads the trigger decision which is
2641 // stored in Trigger.root file and fills
2642 // the corresponding esd entries
2643
52dd4a8c 2644 AliCodeTimerAuto("",0)
87932dab 2645
b647652d 2646 AliInfo("Filling trigger information into the ESD");
2647
2648 if (fRawReader) {
2649 AliCTPRawStream input(fRawReader);
2650 if (!input.Next()) {
7e88424f 2651 AliWarning("No valid CTP (trigger) DDL raw data is found ! The trigger info is taken from the event header!");
b647652d 2652 }
2653 else {
7e88424f 2654 if (esd->GetTriggerMask() != input.GetClassMask())
2655 AliError(Form("Invalid trigger pattern found in CTP raw-data: %llx %llx",
2656 input.GetClassMask(),esd->GetTriggerMask()));
2657 if (esd->GetOrbitNumber() != input.GetOrbitID())
2658 AliError(Form("Invalid orbit id found in CTP raw-data: %x %x",
2659 input.GetOrbitID(),esd->GetOrbitNumber()));
2660 if (esd->GetBunchCrossNumber() != input.GetBCID())
2661 AliError(Form("Invalid bunch-crossing id found in CTP raw-data: %x %x",
2662 input.GetBCID(),esd->GetBunchCrossNumber()));
e61ed4b1 2663 AliESDHeader* esdheader = esd->GetHeader();
2664 esdheader->SetL0TriggerInputs(input.GetL0Inputs());
2665 esdheader->SetL1TriggerInputs(input.GetL1Inputs());
2666 esdheader->SetL2TriggerInputs(input.GetL2Inputs());
a6dd87ad 2667 // IR
2668 UInt_t orbit=input.GetOrbitID();
2669 for(Int_t i=0 ; i<input.GetNIRs() ; i++ )
2670 if(TMath::Abs(Int_t(orbit-(input.GetIR(i))->GetOrbit()))<=1){
2671 esdheader->AddTriggerIR(input.GetIR(i));
2672 }
b647652d 2673 }
b024fd7f 2674 }
f7812afc 2675 return kTRUE;
2676}
2677//_____________________________________________________________________________
2678Bool_t AliReconstruction::FillTriggerScalers(AliESDEvent*& esd)
2679{
522fdd91 2680 //Scalers
82ebedd6 2681 //fRunScalers->Print();
2682 if(fRunScalers && fRunScalers->CheckRunScalers()){
a6dd87ad 2683 AliTimeStamp* timestamp = new AliTimeStamp(esd->GetOrbitNumber(), esd->GetPeriodNumber(), esd->GetBunchCrossNumber());
82ebedd6 2684 //AliTimeStamp* timestamp = new AliTimeStamp(10308000, 0, (ULong64_t)486238);
522fdd91 2685 AliESDHeader* esdheader = fesd->GetHeader();
2686 for(Int_t i=0;i<50;i++){
1e78ae8c 2687 if((1ull<<i) & esd->GetTriggerMask()){
6863d231 2688 AliTriggerScalersESD* scalesd = fRunScalers->GetScalersForEventClass( timestamp, i+1);
82ebedd6 2689 if(scalesd)esdheader->SetTriggerScalersRecord(scalesd);
522fdd91 2690 }
2691 }
2692 }
b647652d 2693 return kTRUE;
2694}
001397cd 2695//_____________________________________________________________________________
af885e0f 2696Bool_t AliReconstruction::FillRawEventHeaderESD(AliESDEvent*& esd)
001397cd 2697{
2698 //
2699 // Filling information from RawReader Header
2700 //
2701
151bea4e 2702 if (!fRawReader) return kFALSE;
2703
001397cd 2704 AliInfo("Filling information from RawReader Header");
151bea4e 2705
2706 esd->SetBunchCrossNumber(fRawReader->GetBCID());
2707 esd->SetOrbitNumber(fRawReader->GetOrbitID());
2708 esd->SetPeriodNumber(fRawReader->GetPeriod());
2709
2710 esd->SetTimeStamp(fRawReader->GetTimestamp());
2711 esd->SetEventType(fRawReader->GetType());
001397cd 2712
2713 return kTRUE;
2714}
2715
2716
596a855f 2717//_____________________________________________________________________________
2718Bool_t AliReconstruction::IsSelected(TString detName, TString& detectors) const
2719{
2720// check whether detName is contained in detectors
2721// if yes, it is removed from detectors
2722
2723 // check if all detectors are selected
2724 if ((detectors.CompareTo("ALL") == 0) ||
2725 detectors.BeginsWith("ALL ") ||
2726 detectors.EndsWith(" ALL") ||
2727 detectors.Contains(" ALL ")) {
2728 detectors = "ALL";
2729 return kTRUE;
2730 }
2731
2732 // search for the given detector
2733 Bool_t result = kFALSE;
2734 if ((detectors.CompareTo(detName) == 0) ||
2735 detectors.BeginsWith(detName+" ") ||
2736 detectors.EndsWith(" "+detName) ||
2737 detectors.Contains(" "+detName+" ")) {
2738 detectors.ReplaceAll(detName, "");
2739 result = kTRUE;
2740 }
2741
2742 // clean up the detectors string
2743 while (detectors.Contains(" ")) detectors.ReplaceAll(" ", " ");
2744 while (detectors.BeginsWith(" ")) detectors.Remove(0, 1);
2745 while (detectors.EndsWith(" ")) detectors.Remove(detectors.Length()-1, 1);
2746
2747 return result;
2748}
e583c30d 2749
f08fc9f5 2750//_____________________________________________________________________________
2751Bool_t AliReconstruction::InitRunLoader()
2752{
2753// get or create the run loader
2754
2755 if (gAlice) delete gAlice;
2756 gAlice = NULL;
2757
52dd4a8c 2758 TFile *gafile = TFile::Open(fGAliceFileName.Data());
2759 // if (!gSystem->AccessPathName(fGAliceFileName.Data())) { // galice.root exists
2760 if (gafile) { // galice.root exists
2761 gafile->Close();
2762 delete gafile;
2763
b26c3770 2764 // load all base libraries to get the loader classes
2765 TString libs = gSystem->GetLibraries();
ac4a7581 2766 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b26c3770 2767 TString detName = fgkDetectorName[iDet];
2768 if (detName == "HLT") continue;
2769 if (libs.Contains("lib" + detName + "base.so")) continue;
2770 gSystem->Load("lib" + detName + "base.so");
2771 }
f08fc9f5 2772 fRunLoader = AliRunLoader::Open(fGAliceFileName.Data());
2773 if (!fRunLoader) {
2774 AliError(Form("no run loader found in file %s", fGAliceFileName.Data()));
2775 CleanUp();
2776 return kFALSE;
2777 }
325aa76f 2778
b26c3770 2779 fRunLoader->CdGAFile();
325aa76f 2780 fRunLoader->LoadgAlice();
f08fc9f5 2781
6cae184e 2782 //PH This is a temporary fix to give access to the kinematics
2783 //PH that is needed for the labels of ITS clusters
f2ee4290 2784 fRunLoader->LoadHeader();
6cae184e 2785 fRunLoader->LoadKinematics();
2786
f08fc9f5 2787 } else { // galice.root does not exist
2788 if (!fRawReader) {
2789 AliError(Form("the file %s does not exist", fGAliceFileName.Data()));
f08fc9f5 2790 }
2791 fRunLoader = AliRunLoader::Open(fGAliceFileName.Data(),
2792 AliConfig::GetDefaultEventFolderName(),
2793 "recreate");
2794 if (!fRunLoader) {
2795 AliError(Form("could not create run loader in file %s",
2796 fGAliceFileName.Data()));
2797 CleanUp();
2798 return kFALSE;
2799 }
21a3aa09 2800 fIsNewRunLoader = kTRUE;
f08fc9f5 2801 fRunLoader->MakeTree("E");
21a3aa09 2802
973388c2 2803 if (fNumberOfEventsPerFile > 0)
2804 fRunLoader->SetNumberOfEventsPerFile(fNumberOfEventsPerFile);
2805 else
21a3aa09 2806 fRunLoader->SetNumberOfEventsPerFile((UInt_t)-1);
f08fc9f5 2807 }
2808
2809 return kTRUE;
2810}
2811
c757bafd 2812//_____________________________________________________________________________
b8cd5251 2813AliReconstructor* AliReconstruction::GetReconstructor(Int_t iDet)
c757bafd 2814{
f08fc9f5 2815// get the reconstructor object and the loader for a detector
c757bafd 2816
7e88424f 2817 if (fReconstructor[iDet]) {
2818 if (fRecoParam.GetDetRecoParamArray(iDet) && !AliReconstructor::GetRecoParam(iDet)) {
2819 const AliDetectorRecoParam *par = fRecoParam.GetDetRecoParam(iDet);
2820 fReconstructor[iDet]->SetRecoParam(par);
42ec5d3d 2821 fReconstructor[iDet]->SetRunInfo(fRunInfo);
7e88424f 2822 }
2823 return fReconstructor[iDet];
2824 }
b8cd5251 2825
2826 // load the reconstructor object
2827 TPluginManager* pluginManager = gROOT->GetPluginManager();
2828 TString detName = fgkDetectorName[iDet];
2829 TString recName = "Ali" + detName + "Reconstructor";
f0999a9a 2830
2831 if (!fIsNewRunLoader && !fRunLoader->GetLoader(detName+"Loader") && (detName != "HLT")) return NULL;
b8cd5251 2832
b8cd5251 2833 AliReconstructor* reconstructor = NULL;
2834 // first check if a plugin is defined for the reconstructor
2835 TPluginHandler* pluginHandler =
2836 pluginManager->FindHandler("AliReconstructor", detName);
f08fc9f5 2837 // if not, add a plugin for it
2838 if (!pluginHandler) {
b8cd5251 2839 AliDebug(1, Form("defining plugin for %s", recName.Data()));
b26c3770 2840 TString libs = gSystem->GetLibraries();
2841 if (libs.Contains("lib" + detName + "base.so") ||
2842 (gSystem->Load("lib" + detName + "base.so") >= 0)) {
b8cd5251 2843 pluginManager->AddHandler("AliReconstructor", detName,
2844 recName, detName + "rec", recName + "()");
2845 } else {
2846 pluginManager->AddHandler("AliReconstructor", detName,
2847 recName, detName, recName + "()");
c757bafd 2848 }
b8cd5251 2849 pluginHandler = pluginManager->FindHandler("AliReconstructor", detName);
2850 }
2851 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
2852 reconstructor = (AliReconstructor*) pluginHandler->ExecPlugin(0);
c757bafd 2853 }
b8cd5251 2854 if (reconstructor) {
2855 TObject* obj = fOptions.FindObject(detName.Data());
2856 if (obj) reconstructor->SetOption(obj->GetTitle());
1e500f25 2857 reconstructor->SetRunInfo(fRunInfo);
d76c31f4 2858 reconstructor->Init();
b8cd5251 2859 fReconstructor[iDet] = reconstructor;
2860 }
2861
f08fc9f5 2862 // get or create the loader
2863 if (detName != "HLT") {
2864 fLoader[iDet] = fRunLoader->GetLoader(detName + "Loader");
2865 if (!fLoader[iDet]) {
2866 AliConfig::Instance()
2867 ->CreateDetectorFolders(fRunLoader->GetEventFolder(),
2868 detName, detName);
2869 // first check if a plugin is defined for the loader
bb0901a4 2870 pluginHandler =
f08fc9f5 2871 pluginManager->FindHandler("AliLoader", detName);
2872 // if not, add a plugin for it
2873 if (!pluginHandler) {
2874 TString loaderName = "Ali" + detName + "Loader";
2875 AliDebug(1, Form("defining plugin for %s", loaderName.Data()));
2876 pluginManager->AddHandler("AliLoader", detName,
2877 loaderName, detName + "base",
2878 loaderName + "(const char*, TFolder*)");
2879 pluginHandler = pluginManager->FindHandler("AliLoader", detName);
2880 }
2881 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
2882 fLoader[iDet] =
2883 (AliLoader*) pluginHandler->ExecPlugin(2, detName.Data(),
2884 fRunLoader->GetEventFolder());
2885 }
2886 if (!fLoader[iDet]) { // use default loader
2887 fLoader[iDet] = new AliLoader(detName, fRunLoader->GetEventFolder());
2888 }
2889 if (!fLoader[iDet]) {
2890 AliWarning(Form("couldn't get loader for %s", detName.Data()));
6667b602 2891 if (fStopOnError) return NULL;
f08fc9f5 2892 } else {
2893 fRunLoader->AddLoader(fLoader[iDet]);
2894 fRunLoader->CdGAFile();
2895 if (gFile && !gFile->IsWritable()) gFile->ReOpen("UPDATE");
2896 fRunLoader->Write(0, TObject::kOverwrite);
2897 }
2898 }
2899 }
2900
7e88424f 2901 if (fRecoParam.GetDetRecoParamArray(iDet) && !AliReconstructor::GetRecoParam(iDet)) {
2902 const AliDetectorRecoParam *par = fRecoParam.GetDetRecoParam(iDet);
2903 reconstructor->SetRecoParam(par);
42ec5d3d 2904 reconstructor->SetRunInfo(fRunInfo);
7e88424f 2905 }
b8cd5251 2906 return reconstructor;
c757bafd 2907}
2908
2257f27e 2909//_____________________________________________________________________________
92bffc4d 2910AliVertexer* AliReconstruction::CreateVertexer()
2257f27e 2911{
2912// create the vertexer
92bffc4d 2913// Please note that the caller is the owner of the
2914// vertexer
2257f27e 2915
92bffc4d 2916 AliVertexer* vertexer = NULL;
b8cd5251 2917 AliReconstructor* itsReconstructor = GetReconstructor(0);
5ec10f5d 2918 if (itsReconstructor && ((fRunLocalReconstruction.Contains("ITS")) || fRunTracking.Contains("ITS"))) {
92bffc4d 2919 vertexer = itsReconstructor->CreateVertexer();
2257f27e 2920 }
92bffc4d 2921 if (!vertexer) {
815c2b38 2922 AliWarning("couldn't create a vertexer for ITS");
2257f27e 2923 }
2924
92bffc4d 2925 return vertexer;
2257f27e 2926}
2927
24f7a148 2928//_____________________________________________________________________________
b8cd5251 2929Bool_t AliReconstruction::CreateTrackers(const TString& detectors)
24f7a148 2930{
f08fc9f5 2931// create the trackers
44e45fac 2932 AliInfo("Creating trackers");
24f7a148 2933
b8cd5251 2934 TString detStr = detectors;
ac4a7581 2935 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 2936 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
2937 AliReconstructor* reconstructor = GetReconstructor(iDet);
2938 if (!reconstructor) continue;
2939 TString detName = fgkDetectorName[iDet];
1f46a9ae 2940 if (detName == "HLT") {
2941 fRunHLTTracking = kTRUE;
2942 continue;
2943 }
e66fbafb 2944 if (detName == "MUON") {
2945 fRunMuonTracking = kTRUE;
2946 continue;
2947 }
2948
f08fc9f5 2949
d76c31f4 2950 fTracker[iDet] = reconstructor->CreateTracker();
f08fc9f5 2951 if (!fTracker[iDet] && (iDet < 7)) {
2952 AliWarning(Form("couldn't create a tracker for %s", detName.Data()));
8250d5f5 2953 if (fStopOnError) return kFALSE;
2954 }
6efecea1 2955 AliSysInfo::AddStamp(Form("LTracker%s",fgkDetectorName[iDet]), iDet,0);
8250d5f5 2956 }
2957
24f7a148 2958 return kTRUE;
2959}
2960
e583c30d 2961//_____________________________________________________________________________
4b71572b 2962void AliReconstruction::CleanUp()
e583c30d 2963{
2964// delete trackers and the run loader and close and delete the file
2965
ac4a7581 2966 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 2967 delete fReconstructor[iDet];
2968 fReconstructor[iDet] = NULL;
2969 fLoader[iDet] = NULL;
2970 delete fTracker[iDet];
2971 fTracker[iDet] = NULL;
2972 }
4b71572b 2973 delete fRunInfo;
7e88424f 2974 fRunInfo = NULL;
2975
58e8dc31 2976 delete fSPDTrackleter;
2977 fSPDTrackleter = NULL;
2978
4b71572b 2979 delete ftVertexer;
21a3aa09 2980 ftVertexer = NULL;
795e4a22 2981
e583c30d 2982 delete fRunLoader;
2983 fRunLoader = NULL;
b649205a 2984 delete fRawReader;
2985 fRawReader = NULL;
4b71572b 2986 delete fParentRawReader;
cd0b062e 2987 fParentRawReader=NULL;
e583c30d 2988
4b71572b 2989 if (ffile) {
2990 ffile->Close();
2991 delete ffile;
2992 ffile = NULL;
24f7a148 2993 }
87da0921 2994
bf76b847 2995 if (AliQAManager::QAManager())
2996 AliQAManager::QAManager()->ShowQA() ;
eca4fa66 2997 AliQAManager::Destroy() ;
2998
24f7a148 2999}
f3a97c86 3000
af885e0f 3001void AliReconstruction::WriteAlignmentData(AliESDEvent* esd)
98937d93 3002{
3003 // Write space-points which are then used in the alignment procedures
6e65648b 3004 // For the moment only ITS, TPC, TRD and TOF
98937d93 3005
98937d93 3006 Int_t ntracks = esd->GetNumberOfTracks();
3007 for (Int_t itrack = 0; itrack < ntracks; itrack++)
3008 {
3009 AliESDtrack *track = esd->GetTrack(itrack);
3010 Int_t nsp = 0;
ef7253ac 3011 Int_t idx[200];
76741755 3012 for (Int_t i=0; i<200; ++i) idx[i] = -1; //PH avoid uninitialized values
81aa7a0d 3013 for (Int_t iDet = 5; iDet >= 0; iDet--) {// TOF, TRD, TPC, ITS clusters
6e65648b 3014 nsp += track->GetNcls(iDet);
3015
3016 if (iDet==0) { // ITS "extra" clusters
3017 track->GetClusters(iDet,idx);
3018 for (Int_t i=6; i<12; i++) if(idx[i] >= 0) nsp++;
3019 }
3020 }
3021
98937d93 3022 if (nsp) {
3023 AliTrackPointArray *sp = new AliTrackPointArray(nsp);
3024 track->SetTrackPointArray(sp);
3025 Int_t isptrack = 0;
81aa7a0d 3026 for (Int_t iDet = 5; iDet >= 0; iDet--) {
98937d93 3027 AliTracker *tracker = fTracker[iDet];
3028 if (!tracker) continue;
6e65648b 3029 Int_t nspdet = track->GetClusters(iDet,idx);
3030
3031 if (iDet==0) // ITS "extra" clusters
3032 for (Int_t i=6; i<12; i++) if(idx[i] >= 0) nspdet++;
3033
98937d93 3034 if (nspdet <= 0) continue;
98937d93 3035 AliTrackPoint p;
3036 Int_t isp = 0;
3037 Int_t isp2 = 0;
4ed6fb1c 3038 while (isp2 < nspdet) {
f3c6e4c9 3039 Bool_t isvalid=kTRUE;
3040
3041 Int_t index=idx[isp++];
3042 if (index < 0) continue;
3043
c12b6e44 3044 TString dets = fgkDetectorName[iDet];
3045 if ((fUseTrackingErrorsForAlignment.CompareTo(dets) == 0) ||
3046 fUseTrackingErrorsForAlignment.BeginsWith(dets+" ") ||
3047 fUseTrackingErrorsForAlignment.EndsWith(" "+dets) ||
3048 fUseTrackingErrorsForAlignment.Contains(" "+dets+" ")) {
f3c6e4c9 3049 isvalid = tracker->GetTrackPointTrackingError(index,p,track);
48ce48d1 3050 } else {
f3c6e4c9 3051 isvalid = tracker->GetTrackPoint(index,p);
48ce48d1 3052 }
3053 isp2++;
98937d93 3054 if (!isvalid) continue;
eb2d90b9 3055 if (iDet==0 && (isp-1)>=6) p.SetExtra();
f3c6e4c9 3056 sp->AddPoint(isptrack,&p); isptrack++;
98937d93 3057 }
98937d93 3058 }
3059 }
3060 }
98937d93 3061}
2e3550da 3062
3063//_____________________________________________________________________________
af885e0f 3064void AliReconstruction::FillRawDataErrorLog(Int_t iEvent, AliESDEvent* esd)
2e3550da 3065{
3066 // The method reads the raw-data error log
3067 // accumulated within the rawReader.
3068 // It extracts the raw-data errors related to
3069 // the current event and stores them into
3070 // a TClonesArray inside the esd object.
3071
3072 if (!fRawReader) return;
3073
3074 for(Int_t i = 0; i < fRawReader->GetNumberOfErrorLogs(); i++) {
3075
3076 AliRawDataErrorLog *log = fRawReader->GetErrorLog(i);
3077 if (!log) continue;
3078 if (iEvent != log->GetEventNumber()) continue;
3079
3080 esd->AddRawDataErrorLog(log);
3081 }
3082
3083}
46698ae4 3084
8661738e 3085//_____________________________________________________________________________
3086void AliReconstruction::CheckQA()
3087{
3088// check the QA of SIM for this run and remove the detectors
3089// with status Fatal
3090
57acd2d2 3091// TString newRunLocalReconstruction ;
3092// TString newRunTracking ;
3093// TString newFillESD ;
3094//
4e25ac79 3095// for (Int_t iDet = 0; iDet < AliQAv1::kNDET; iDet++) {
3096// TString detName(AliQAv1::GetDetName(iDet)) ;
3097// AliQAv1 * qa = AliQAv1::Instance(AliQAv1::DETECTORINDEX_t(iDet)) ;
3098// if ( qa->IsSet(AliQAv1::DETECTORINDEX_t(iDet), AliQAv1::kSIM, specie, AliQAv1::kFATAL)) {
57acd2d2 3099// AliInfo(Form("QA status for %s %s in Hits and/or SDIGITS and/or Digits was Fatal; No reconstruction performed",
3100// detName.Data(), AliRecoParam::GetEventSpecieName(es))) ;
3101// } else {
4e25ac79 3102// if ( fRunLocalReconstruction.Contains(AliQAv1::GetDetName(iDet)) ||
57acd2d2 3103// fRunLocalReconstruction.Contains("ALL") ) {
3104// newRunLocalReconstruction += detName ;
3105// newRunLocalReconstruction += " " ;
3106// }
4e25ac79 3107// if ( fRunTracking.Contains(AliQAv1::GetDetName(iDet)) ||
57acd2d2 3108// fRunTracking.Contains("ALL") ) {
3109// newRunTracking += detName ;
3110// newRunTracking += " " ;
3111// }
4e25ac79 3112// if ( fFillESD.Contains(AliQAv1::GetDetName(iDet)) ||
57acd2d2 3113// fFillESD.Contains("ALL") ) {
3114// newFillESD += detName ;
3115// newFillESD += " " ;
3116// }
3117// }
3118// }
3119// fRunLocalReconstruction = newRunLocalReconstruction ;
3120// fRunTracking = newRunTracking ;
3121// fFillESD = newFillESD ;
a5fa6165 3122}
5b188f2f 3123
3124//_____________________________________________________________________________
3125Int_t AliReconstruction::GetDetIndex(const char* detector)
3126{
3127 // return the detector index corresponding to detector
3128 Int_t index = -1 ;
ac4a7581 3129 for (index = 0; index < kNDetectors ; index++) {
5b188f2f 3130 if ( strcmp(detector, fgkDetectorName[index]) == 0 )
3131 break ;
3132 }
3133 return index ;
3134}
7167ae53 3135//_____________________________________________________________________________
3136Bool_t AliReconstruction::FinishPlaneEff() {
3137 //
3138 // Here execute all the necessary operationis, at the end of the tracking phase,
d7f8fd68 3139 // in case that evaluation of PlaneEfficiencies was required for some detector.
3140 // E.g., write into a DataBase file the PlaneEfficiency which have been evaluated.
7167ae53 3141 //
3142 // This Preliminary version works only FOR ITS !!!!!
3143 // other detectors (TOF,TRD, etc. have to develop their specific codes)
3144 //
3145 // Input: none
d7f8fd68 3146 // Return: kTRUE if all operations have been done properly, kFALSE otherwise
7167ae53 3147 //
3148 Bool_t ret=kFALSE;
58e8dc31 3149 //for (Int_t iDet = 0; iDet < fgkNDetectors; iDet++) {
d7f8fd68 3150 for (Int_t iDet = 0; iDet < 1; iDet++) { // for the time being only ITS
7167ae53 3151 //if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
58e8dc31 3152 if(fTracker[iDet] && fTracker[iDet]->GetPlaneEff()) {
d7f8fd68 3153 AliPlaneEff *planeeff=fTracker[iDet]->GetPlaneEff();
3154 TString name=planeeff->GetName();
3155 name+=".root";
3156 TFile* pefile = TFile::Open(name, "RECREATE");
3157 ret=(Bool_t)planeeff->Write();
3158 pefile->Close();
5fbd4fd6 3159 if(planeeff->GetCreateHistos()) {
d7f8fd68 3160 TString hname=planeeff->GetName();
3161 hname+="Histo.root";
3162 ret*=planeeff->WriteHistosToFile(hname,"RECREATE");
5fbd4fd6 3163 }
3164 }
58e8dc31 3165 if(fSPDTrackleter) {
3166 AliPlaneEff *planeeff=fSPDTrackleter->GetPlaneEff();
3167 TString name="AliITSPlaneEffSPDtracklet.root";
3168 TFile* pefile = TFile::Open(name, "RECREATE");
3169 ret=(Bool_t)planeeff->Write();
3170 pefile->Close();
3171 AliESDEvent *dummy=NULL;
3172 ret=(Bool_t)fSPDTrackleter->PostProcess(dummy); // take care of writing other files
3173 }
7167ae53 3174 }
3175 return ret;
3176}
3177//_____________________________________________________________________________
3178Bool_t AliReconstruction::InitPlaneEff() {
3179//
3180 // Here execute all the necessary operations, before of the tracking phase,
3181 // for the evaluation of PlaneEfficiencies, in case required for some detectors.
58e8dc31 3182 // E.g., read from a DataBase file a first evaluation of the PlaneEfficiency
7167ae53 3183 // which should be updated/recalculated.
3184 //
3185 // This Preliminary version will work only FOR ITS !!!!!
3186 // other detectors (TOF,TRD, etc. have to develop their specific codes)
3187 //
3188 // Input: none
3189 // Return: kTRUE if all operations have been done properly, kFALSE otherwise
3190 //
58e8dc31 3191 AliWarning(Form("Implementation of this method not yet completed !! Method return kTRUE"));
3192
3193 fSPDTrackleter = NULL;
3194 AliReconstructor* itsReconstructor = GetReconstructor(0);
3195 if (itsReconstructor) {
3196 fSPDTrackleter = itsReconstructor->CreateTrackleter(); // this is NULL unless required in RecoParam
3197 }
3198 if (fSPDTrackleter) {
3199 AliInfo("Trackleter for SPD has been created");
3200 }
3201
7167ae53 3202 return kTRUE;
7520312d 3203}
14dd053c 3204
3205//_____________________________________________________________________________
3206Bool_t AliReconstruction::InitAliEVE()
3207{
3208 // This method should be called only in case
3209 // AliReconstruction is run
3210 // within the alieve environment.
3211 // It will initialize AliEVE in a way
3212 // so that it can visualize event processed
3213 // by AliReconstruction.
3214 // The return flag shows whenever the
3215 // AliEVE initialization was successful or not.
3216
3217 TString macroStr;
3218 macroStr.Form("%s/EVE/macros/alieve_online.C",gSystem->ExpandPathName("$ALICE_ROOT"));
3219 AliInfo(Form("Loading AliEVE macro: %s",macroStr.Data()));
3220 if (gROOT->LoadMacro(macroStr.Data()) != 0) return kFALSE;
3221
6a840f1e 3222 gROOT->ProcessLine("if (!AliEveEventManager::GetMaster()){new AliEveEventManager();AliEveEventManager::GetMaster()->AddNewEventCommand(\"alieve_online_on_new_event()\");gEve->AddEvent(AliEveEventManager::GetMaster());};");
de33999e 3223 gROOT->ProcessLine("alieve_online_init()");
14dd053c 3224
3225 return kTRUE;
3226}
3227
3228//_____________________________________________________________________________
3229void AliReconstruction::RunAliEVE()
3230{
3231 // Runs AliEVE visualisation of
3232 // the current event.
3233 // Should be executed only after
3234 // successful initialization of AliEVE.
3235
3236 AliInfo("Running AliEVE...");
519975fe 3237 gROOT->ProcessLine(Form("AliEveEventManager::GetMaster()->SetEvent((AliRunLoader*)0x%lx,(AliRawReader*)0x%lx,(AliESDEvent*)0x%lx,(AliESDfriend*)0x%lx);",fRunLoader,fRawReader,fesd,fesdf));
14dd053c 3238 gSystem->Run();
3239}
ce43afbe 3240
3241//_____________________________________________________________________________
3242Bool_t AliReconstruction::SetRunQA(TString detAndAction)
3243{
3244 // Allows to run QA for a selected set of detectors
44ed7a66 3245 // and a selected set of tasks among RAWS, DIGITSR, RECPOINTS and ESDS
ce43afbe 3246 // all selected detectors run the same selected tasks
3247
3248 if (!detAndAction.Contains(":")) {
3249 AliError( Form("%s is a wrong syntax, use \"DetectorList:ActionList\" \n", detAndAction.Data()) ) ;
3250 fRunQA = kFALSE ;
3251 return kFALSE ;
3252 }
3253 Int_t colon = detAndAction.Index(":") ;
3254 fQADetectors = detAndAction(0, colon) ;
e574cc31 3255 if (fQADetectors.Contains("ALL") ){
3256 TString tmp = fFillESD ;
3257 Int_t minus = fQADetectors.Last('-') ;
3258 TString toKeep = fFillESD ;
3259 TString toRemove("") ;
3260 while (minus >= 0) {
3261 toRemove = fQADetectors(minus+1, fQADetectors.Length()) ;
3262 toRemove = toRemove.Strip() ;
3263 toKeep.ReplaceAll(toRemove, "") ;
3264 fQADetectors.ReplaceAll(Form("-%s", toRemove.Data()), "") ;
3265 minus = fQADetectors.Last('-') ;
3266 }
3267 fQADetectors = toKeep ;
3268 }
3269 fQATasks = detAndAction(colon+1, detAndAction.Sizeof() ) ;
ce43afbe 3270 if (fQATasks.Contains("ALL") ) {
44ed7a66 3271 fQATasks = Form("%d %d %d %d", AliQAv1::kRAWS, AliQAv1::kDIGITSR, AliQAv1::kRECPOINTS, AliQAv1::kESDS) ;
ce43afbe 3272 } else {
3273 fQATasks.ToUpper() ;
3274 TString tempo("") ;
3275 if ( fQATasks.Contains("RAW") )
4e25ac79 3276 tempo = Form("%d ", AliQAv1::kRAWS) ;
44ed7a66 3277 if ( fQATasks.Contains("DIGIT") )
3278 tempo += Form("%d ", AliQAv1::kDIGITSR) ;
ce43afbe 3279 if ( fQATasks.Contains("RECPOINT") )
4e25ac79 3280 tempo += Form("%d ", AliQAv1::kRECPOINTS) ;
ce43afbe 3281 if ( fQATasks.Contains("ESD") )
4e25ac79 3282 tempo += Form("%d ", AliQAv1::kESDS) ;
ce43afbe 3283 fQATasks = tempo ;
3284 if (fQATasks.IsNull()) {
3285 AliInfo("No QA requested\n") ;
3286 fRunQA = kFALSE ;
3287 return kTRUE ;
3288 }
3289 }
3290 TString tempo(fQATasks) ;
4e25ac79 3291 tempo.ReplaceAll(Form("%d", AliQAv1::kRAWS), AliQAv1::GetTaskName(AliQAv1::kRAWS)) ;
44ed7a66 3292 tempo.ReplaceAll(Form("%d", AliQAv1::kDIGITSR), AliQAv1::GetTaskName(AliQAv1::kDIGITSR)) ;
4e25ac79 3293 tempo.ReplaceAll(Form("%d", AliQAv1::kRECPOINTS), AliQAv1::GetTaskName(AliQAv1::kRECPOINTS)) ;
3294 tempo.ReplaceAll(Form("%d", AliQAv1::kESDS), AliQAv1::GetTaskName(AliQAv1::kESDS)) ;
ce43afbe 3295 AliInfo( Form("QA will be done on \"%s\" for \"%s\"\n", fQADetectors.Data(), tempo.Data()) ) ;
3296 fRunQA = kTRUE ;
3297 return kTRUE;
3298}
3299
7e88424f 3300//_____________________________________________________________________________
3301Bool_t AliReconstruction::InitRecoParams()
3302{
3303 // The method accesses OCDB and retrieves all
3304 // the available reco-param objects from there.
3305
3306 Bool_t isOK = kTRUE;
3307
8b26452d 3308 if (fRecoParam.GetDetRecoParamArray(kNDetectors)) {
3309 AliInfo("Using custom GRP reconstruction parameters");
3310 }
3311 else {
3312 AliInfo("Loading GRP reconstruction parameter objects");
3313
3314 AliCDBPath path("GRP","Calib","RecoParam");
3315 AliCDBEntry *entry=AliCDBManager::Instance()->Get(path.GetPath());
3316 if(!entry){
3317 AliWarning("Couldn't find GRP RecoParam entry in OCDB");
3318 isOK = kFALSE;
3319 }
3320 else {
3321 TObject *recoParamObj = entry->GetObject();
3322 if (dynamic_cast<TObjArray*>(recoParamObj)) {
3323 // GRP has a normal TobjArray of AliDetectorRecoParam objects
3324 // Registering them in AliRecoParam
3325 fRecoParam.AddDetRecoParamArray(kNDetectors,dynamic_cast<TObjArray*>(recoParamObj));
3326 }
3327 else if (dynamic_cast<AliDetectorRecoParam*>(recoParamObj)) {
3328 // GRP has only onse set of reco parameters
3329 // Registering it in AliRecoParam
3330 AliInfo("Single set of GRP reconstruction parameters found");
3331 dynamic_cast<AliDetectorRecoParam*>(recoParamObj)->SetAsDefault();
3332 fRecoParam.AddDetRecoParam(kNDetectors,dynamic_cast<AliDetectorRecoParam*>(recoParamObj));
3333 }
3334 else {
3335 AliError("No valid GRP RecoParam object found in the OCDB");
3336 isOK = kFALSE;
3337 }
3338 entry->SetOwner(0);
3339 }
3340 }
3341
7d566c20 3342 TString detStr = fLoadCDB;
ac4a7581 3343 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
7e88424f 3344
7d566c20 3345 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
3346
7e88424f 3347 if (fRecoParam.GetDetRecoParamArray(iDet)) {
3348 AliInfo(Form("Using custom reconstruction parameters for detector %s",fgkDetectorName[iDet]));
3349 continue;
3350 }
3351
ac232c75 3352 AliInfo(Form("Loading reconstruction parameter objects for detector %s",fgkDetectorName[iDet]));
7e88424f 3353
3354 AliCDBPath path(fgkDetectorName[iDet],"Calib","RecoParam");
3355 AliCDBEntry *entry=AliCDBManager::Instance()->Get(path.GetPath());
3356 if(!entry){
3357 AliWarning(Form("Couldn't find RecoParam entry in OCDB for detector %s",fgkDetectorName[iDet]));
3358 isOK = kFALSE;
3359 }
3360 else {
3361 TObject *recoParamObj = entry->GetObject();
3362 if (dynamic_cast<TObjArray*>(recoParamObj)) {
3363 // The detector has a normal TobjArray of AliDetectorRecoParam objects
3364 // Registering them in AliRecoParam
3365 fRecoParam.AddDetRecoParamArray(iDet,dynamic_cast<TObjArray*>(recoParamObj));
3366 }
3367 else if (dynamic_cast<AliDetectorRecoParam*>(recoParamObj)) {
3368 // The detector has only onse set of reco parameters
3369 // Registering it in AliRecoParam
ac232c75 3370 AliInfo(Form("Single set of reconstruction parameters found for detector %s",fgkDetectorName[iDet]));
7e88424f 3371 dynamic_cast<AliDetectorRecoParam*>(recoParamObj)->SetAsDefault();
3372 fRecoParam.AddDetRecoParam(iDet,dynamic_cast<AliDetectorRecoParam*>(recoParamObj));
3373 }
3374 else {
3375 AliError(Form("No valid RecoParam object found in the OCDB for detector %s",fgkDetectorName[iDet]));
3376 isOK = kFALSE;
3377 }
3378 entry->SetOwner(0);
f168abba 3379 // FIX ME: We have to disable the unloading of reco-param CDB
3380 // entries because QA framework is using them. Has to be fix in
3381 // a way that the QA takes the objects already constructed in
3382 // this method.
3383 // AliCDBManager::Instance()->UnloadFromCache(path.GetPath());
7e88424f 3384 }
3385 }
3386
e30a9b4d 3387 if (AliDebugLevel() > 0) fRecoParam.Print();
ac232c75 3388
7e88424f 3389 return isOK;
3390}
3391
3392//_____________________________________________________________________________
3393Bool_t AliReconstruction::GetEventInfo()
3394{
3395 // Fill the event info object
3396 // ...
52dd4a8c 3397 AliCodeTimerAuto("",0)
7e88424f 3398
3399 AliCentralTrigger *aCTP = NULL;
3400 if (fRawReader) {
3401 fEventInfo.SetEventType(fRawReader->GetType());
3402
3403 ULong64_t mask = fRawReader->GetClassMask();
3404 fEventInfo.SetTriggerMask(mask);
3405 UInt_t clmask = fRawReader->GetDetectorPattern()[0];
3406 fEventInfo.SetTriggerCluster(AliDAQ::ListOfTriggeredDetectors(clmask));
3407
3408 aCTP = new AliCentralTrigger();
3409 TString configstr("");
3410 if (!aCTP->LoadConfiguration(configstr)) { // Load CTP config from OCDB
3411 AliError("No trigger configuration found in OCDB! The trigger configuration information will not be used!");
3412 delete aCTP;
3413 return kFALSE;
3414 }
3415 aCTP->SetClassMask(mask);
3416 aCTP->SetClusterMask(clmask);
3417 }
3418 else {
3419 fEventInfo.SetEventType(AliRawEventHeaderBase::kPhysicsEvent);
3420
3421 if (fRunLoader && (!fRunLoader->LoadTrigger())) {
3422 aCTP = fRunLoader->GetTrigger();
3423 fEventInfo.SetTriggerMask(aCTP->GetClassMask());
e61ed4b1 3424 // get inputs from actp - just get
3425 AliESDHeader* esdheader = fesd->GetHeader();
3426 esdheader->SetL0TriggerInputs(aCTP->GetL0TriggerInputs());
3427 esdheader->SetL1TriggerInputs(aCTP->GetL1TriggerInputs());
3428 esdheader->SetL2TriggerInputs(aCTP->GetL2TriggerInputs());
7e88424f 3429 fEventInfo.SetTriggerCluster(AliDAQ::ListOfTriggeredDetectors(aCTP->GetClusterMask()));
3430 }
3431 else {
3432 AliWarning("No trigger can be loaded! The trigger information will not be used!");
3433 return kFALSE;
3434 }
3435 }
3436
3437 AliTriggerConfiguration *config = aCTP->GetConfiguration();
3438 if (!config) {
3439 AliError("No trigger configuration has been found! The trigger configuration information will not be used!");
3440 if (fRawReader) delete aCTP;
3441 return kFALSE;
3442 }
3443
a0c2cf2d 3444 UChar_t clustmask = 0;
7e88424f 3445 TString trclasses;
3446 ULong64_t trmask = fEventInfo.GetTriggerMask();
3447 const TObjArray& classesArray = config->GetClasses();
3448 Int_t nclasses = classesArray.GetEntriesFast();
3449 for( Int_t iclass=0; iclass < nclasses; iclass++ ) {
3450 AliTriggerClass* trclass = (AliTriggerClass*)classesArray.At(iclass);
3451 if (trclass) {
a4b0683d 3452 Int_t trindex = TMath::Nint(TMath::Log2(trclass->GetMask()));
7e88424f 3453 fesd->SetTriggerClass(trclass->GetName(),trindex);
8a933107 3454 if (fRawReader) fRawReader->LoadTriggerClass(trclass->GetName(),trindex);
e61ed4b1 3455 if (trmask & (1ull << trindex)) {
7e88424f 3456 trclasses += " ";
3457 trclasses += trclass->GetName();
3458 trclasses += " ";
a0c2cf2d 3459 clustmask |= trclass->GetCluster()->GetClusterMask();
7e88424f 3460 }
3461 }
3462 }
3463 fEventInfo.SetTriggerClasses(trclasses);
3464
a0c2cf2d 3465 // Set the information in ESD
3466 fesd->SetTriggerMask(trmask);
3467 fesd->SetTriggerCluster(clustmask);
3468
7e88424f 3469 if (!aCTP->CheckTriggeredDetectors()) {
3470 if (fRawReader) delete aCTP;
3471 return kFALSE;
3472 }
3473
3474 if (fRawReader) delete aCTP;
3475
3476 // We have to fill also the HLT decision here!!
3477 // ...
3478
3479 return kTRUE;
3480}
3481
3482const char *AliReconstruction::MatchDetectorList(const char *detectorList, UInt_t detectorMask)
3483{
3484 // Match the detector list found in the rec.C or the default 'ALL'
3485 // to the list found in the GRP (stored there by the shuttle PP which
3486 // gets the information from ECS)
3487 static TString resultList;
3488 TString detList = detectorList;
3489
3490 resultList = "";
3491
3492 for(Int_t iDet = 0; iDet < (AliDAQ::kNDetectors-1); iDet++) {
3493 if ((detectorMask >> iDet) & 0x1) {
3494 TString det = AliDAQ::OfflineModuleName(iDet);
3495 if ((detList.CompareTo("ALL") == 0) ||
a101e1dd 3496 ((detList.BeginsWith("ALL ") ||
3497 detList.EndsWith(" ALL") ||
3498 detList.Contains(" ALL ")) &&
3499 !(detList.BeginsWith("-"+det+" ") ||
3500 detList.EndsWith(" -"+det) ||
3501 detList.Contains(" -"+det+" "))) ||
7e88424f 3502 (detList.CompareTo(det) == 0) ||
a101e1dd 3503 detList.BeginsWith(det+" ") ||
3504 detList.EndsWith(" "+det) ||
7e88424f 3505 detList.Contains( " "+det+" " )) {
3506 if (!resultList.EndsWith(det + " ")) {
3507 resultList += det;
3508 resultList += " ";
3509 }
3510 }
3511 }
3512 }
3513
3514 // HLT
3515 if ((detectorMask >> AliDAQ::kHLTId) & 0x1) {
3516 TString hltDet = AliDAQ::OfflineModuleName(AliDAQ::kNDetectors-1);
3517 if ((detList.CompareTo("ALL") == 0) ||
a101e1dd 3518 ((detList.BeginsWith("ALL ") ||
3519 detList.EndsWith(" ALL") ||
3520 detList.Contains(" ALL ")) &&
3521 !(detList.BeginsWith("-"+hltDet+" ") ||
3522 detList.EndsWith(" -"+hltDet) ||
3523 detList.Contains(" -"+hltDet+" "))) ||
7e88424f 3524 (detList.CompareTo(hltDet) == 0) ||
a101e1dd 3525 detList.BeginsWith(hltDet+" ") ||
3526 detList.EndsWith(" "+hltDet) ||
7e88424f 3527 detList.Contains( " "+hltDet+" " )) {
3528 resultList += hltDet;
3529 }
3530 }
3531
3532 return resultList.Data();
3533
3534}
4b71572b 3535
3536//______________________________________________________________________________
3537void AliReconstruction::Abort(const char *method, EAbort what)
3538{
3539 // Abort processing. If what = kAbortProcess, the Process() loop will be
3540 // aborted. If what = kAbortFile, the current file in a chain will be
3541 // aborted and the processing will continue with the next file, if there
3542 // is no next file then Process() will be aborted. Abort() can also be
3543 // called from Begin(), SlaveBegin(), Init() and Notify(). After abort
3544 // the SlaveTerminate() and Terminate() are always called. The abort flag
3545 // can be checked in these methods using GetAbort().
3546 //
3547 // The method is overwritten in AliReconstruction for better handling of
3548 // reco specific errors
3549
3550 if (!fStopOnError) return;
3551
3552 CleanUp();
3553
3554 TString whyMess = method;
3555 whyMess += " failed! Aborting...";
3556
3557 AliError(whyMess.Data());
3558
3559 fAbort = what;
3560 TString mess = "Abort";
3561 if (fAbort == kAbortProcess)
3562 mess = "AbortProcess";
3563 else if (fAbort == kAbortFile)
3564 mess = "AbortFile";
3565
3566 Info(mess, whyMess.Data());
3567}
3568
2f954aba 3569//______________________________________________________________________________
3570Bool_t AliReconstruction::ProcessEvent(void* event)
3571{
3572 // Method that is used in case the event loop
3573 // is steered from outside, for example by AMORE
3574 // 'event' is a pointer to the DATE event in the memory
3575
3576 if (fRawReader) delete fRawReader;
3577 fRawReader = new AliRawReaderDate(event);
3578 fStatus = ProcessEvent(fRunLoader->GetNumberOfEvents());
3579 delete fRawReader;
3580 fRawReader = NULL;
3581
3582 return fStatus;
3583}
81d57268 3584
3585//______________________________________________________________________________
3586Bool_t AliReconstruction::ParseOutput()
3587{
3588 // The method parses the output file
3589 // location string in order to steer
3590 // properly the selector
3591
3592 TPMERegexp re1("(\\w+\\.zip#\\w+\\.root):([,*\\w+\\.root,*]+)@dataset://(\\w++)");
3593 TPMERegexp re2("(\\w+\\.root)?@?dataset://(\\w++)");
3594
3595 if (re1.Match(fESDOutput) == 4) {
3596 // root archive with output files stored and regustered
3597 // in proof dataset
3598 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE",re1[1].Data()));
3599 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE_LOCATION",re1[3].Data()));
3600 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE_DATASET",""));
3601 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE_ARCHIVE",re1[2].Data()));
3602 AliInfo(Form("%s files will be stored within %s in dataset %s",
3603 re1[2].Data(),
3604 re1[1].Data(),
3605 re1[3].Data()));
3606 }
3607 else if (re2.Match(fESDOutput) == 3) {
3608 // output file stored and registered
3609 // in proof dataset
3610 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE",(re2[1].IsNull()) ? "AliESDs.root" : re2[1].Data()));
3611 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE_LOCATION",re2[2].Data()));
3612 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE_DATASET",""));
3613 AliInfo(Form("%s will be stored in dataset %s",
3614 (re2[1].IsNull()) ? "AliESDs.root" : re2[1].Data(),
3615 re2[2].Data()));
3616 }
3617 else {
3618 if (fESDOutput.IsNull()) {
3619 // Output location not given.
3620 // Assuming xrootd has been already started and
3621 // the output file has to be sent back
3622 // to the client machine
3623 TString esdUrl(Form("root://%s/%s/",
3624 TUrl(gSystem->HostName()).GetHostFQDN(),
3625 gSystem->pwd()));
3626 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE","AliESDs.root"));
3627 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE_LOCATION",esdUrl.Data()));
3628 AliInfo(Form("AliESDs.root will be stored in %s",
3629 esdUrl.Data()));
3630 }
3631 else {
3632 // User specified an output location.
3633 // Ones has just to parse it here
3634 TUrl outputUrl(fESDOutput.Data());
3635 TString outputFile(gSystem->BaseName(outputUrl.GetFile()));
3636 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE",outputFile.IsNull() ? "AliESDs.root" : outputFile.Data()));
3637 TString outputLocation(outputUrl.GetUrl());
3638 outputLocation.ReplaceAll(outputFile.Data(),"");
3639 gProof->AddInput(new TNamed("PROOF_OUTPUTFILE_LOCATION",outputLocation.Data()));
3640 AliInfo(Form("%s will be stored in %s",
3641 outputFile.IsNull() ? "AliESDs.root" : outputFile.Data(),
3642 outputLocation.Data()));
3643 }
3644 }
3645
3646 return kTRUE;
3647}