Gstpar removed from geometry code
[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>
0a035be5 132#include <TUrl.h>
596a855f 133
f7a1cc68 134#include "AliAlignObj.h"
135#include "AliCDBEntry.h"
136#include "AliCDBManager.h"
137#include "AliCDBStorage.h"
138#include "AliCTPRawStream.h"
139#include "AliCascadeVertexer.h"
140#include "AliCentralTrigger.h"
87932dab 141#include "AliCodeTimer.h"
f7a1cc68 142#include "AliDAQ.h"
143#include "AliDetectorRecoParam.h"
144#include "AliESDCaloCells.h"
145#include "AliESDCaloCluster.h"
af885e0f 146#include "AliESDEvent.h"
faffd83e 147#include "AliESDMuonTrack.h"
f7a1cc68 148#include "AliESDPmdTrack.h"
149#include "AliESDTagCreator.h"
2257f27e 150#include "AliESDVertex.h"
faffd83e 151#include "AliESDcascade.h"
f7a1cc68 152#include "AliESDfriend.h"
faffd83e 153#include "AliESDkink.h"
596a855f 154#include "AliESDpid.h"
ff8bb5ae 155#include "AliESDtrack.h"
f7a1cc68 156#include "AliESDtrack.h"
157#include "AliEventInfo.h"
158#include "AliGRPObject.h"
159#include "AliGRPRecoParam.h"
160#include "AliGenEventHeader.h"
25be1e5c 161#include "AliGeomManager.h"
aa3c69a9 162#include "AliGlobalQADataMaker.h"
f7a1cc68 163#include "AliHeader.h"
164#include "AliLog.h"
165#include "AliMagF.h"
166#include "AliMultiplicity.h"
167#include "AliPID.h"
168#include "AliPlaneEff.h"
4e25ac79 169#include "AliQAv1.h"
f7a1cc68 170#include "AliQADataMakerRec.h"
b03591ab 171#include "AliQAManager.h"
33314186 172#include "AliRawVEvent.h"
f7a1cc68 173#include "AliRawEventHeaderBase.h"
cd0b062e 174#include "AliRawHLTManager.h"
f7a1cc68 175#include "AliRawReaderDate.h"
176#include "AliRawReaderFile.h"
177#include "AliRawReaderRoot.h"
178#include "AliReconstruction.h"
179#include "AliReconstructor.h"
180#include "AliRun.h"
7e88424f 181#include "AliRunInfo.h"
f7a1cc68 182#include "AliRunLoader.h"
183#include "AliSysInfo.h" // memory snapshots
184#include "AliTrackPointArray.h"
185#include "AliTracker.h"
186#include "AliTriggerClass.h"
187#include "AliTriggerCluster.h"
a6dd87ad 188#include "AliTriggerIR.h"
f7a1cc68 189#include "AliTriggerConfiguration.h"
190#include "AliV0vertexer.h"
191#include "AliVertexer.h"
192#include "AliVertexerTracks.h"
52dd4a8c 193#include "AliTriggerRunScalers.h"
194#include "AliCTPTimeParams.h"
8b12d288 195#include "AliESDHLTDecision.h"
6ef9caeb 196#include "AliTriggerInput.h"
596a855f 197ClassImp(AliReconstruction)
198
596a855f 199//_____________________________________________________________________________
ac4a7581 200const char* AliReconstruction::fgkDetectorName[AliReconstruction::kNDetectors] = {"ITS", "TPC", "TRD", "TOF", "PHOS", "HMPID", "EMCAL", "MUON", "FMD", "ZDC", "PMD", "T0", "VZERO", "ACORDE", "HLT"};
c757bafd 201
202//_____________________________________________________________________________
4b71572b 203AliReconstruction::AliReconstruction(const char* gAliceFilename) :
204 TSelector(),
2257f27e 205 fRunVertexFinder(kTRUE),
a84e2607 206 fRunVertexFinderTracks(kTRUE),
1f46a9ae 207 fRunHLTTracking(kFALSE),
e66fbafb 208 fRunMuonTracking(kFALSE),
d1683eef 209 fRunV0Finder(kTRUE),
210 fRunCascadeFinder(kTRUE),
c8025cc8 211 fStopOnError(kTRUE),
1d99986f 212 fWriteAlignmentData(kFALSE),
213 fWriteESDfriend(kFALSE),
b647652d 214 fFillTriggerESD(kTRUE),
1d99986f 215
7f68891d 216 fCleanESD(kTRUE),
a023d8d8 217 fV0DCAmax(3.),
218 fV0CsPmin(0.),
7f68891d 219 fDmax(50.),
220 fZmax(50.),
221
1d99986f 222 fRunLocalReconstruction("ALL"),
b8cd5251 223 fRunTracking("ALL"),
e583c30d 224 fFillESD("ALL"),
7d566c20 225 fLoadCDB(""),
48ce48d1 226 fUseTrackingErrorsForAlignment(""),
e583c30d 227 fGAliceFileName(gAliceFilename),
4b71572b 228 fRawInput(""),
975320a0 229 fESDOutput(""),
81d57268 230 fProofOutputFileName(""),
231 fProofOutputLocation(""),
232 fProofOutputDataset(kFALSE),
233 fProofOutputArchive(""),
35042093 234 fEquipIdMap(""),
b26c3770 235 fFirstEvent(0),
236 fLastEvent(-1),
9d705769 237 fNumberOfEventsPerFile((UInt_t)-1),
b8cd5251 238 fOptions(),
6bae477a 239 fLoadAlignFromCDB(kTRUE),
240 fLoadAlignData("ALL"),
cd0b062e 241 fUseHLTData(),
7e88424f 242 fRunInfo(NULL),
243 fEventInfo(),
522fdd91 244 fRunScalers(NULL),
52dd4a8c 245 fCTPTimeParams(NULL),
e583c30d 246
247 fRunLoader(NULL),
b649205a 248 fRawReader(NULL),
cd0b062e 249 fParentRawReader(NULL),
b8cd5251 250
7e88424f 251 fRecoParam(),
252
58e8dc31 253 fSPDTrackleter(NULL),
254
f2a195c1 255 fDiamondProfileSPD(NULL),
9178838a 256 fDiamondProfile(NULL),
43c9dae1 257 fDiamondProfileTPC(NULL),
77ba28ba 258 fListOfCosmicTriggers(NULL),
f2a195c1 259
6b6e4472 260 fGRPData(NULL),
444753c6 261
6bae477a 262 fAlignObjArray(NULL),
795e4a22 263 fCDBUri(),
f1c1204d 264 fQARefUri(),
759c1df1 265 fSpecCDBUri(),
795e4a22 266 fInitCDBCalled(kFALSE),
267 fSetRunNumberFromDataCalled(kFALSE),
ce43afbe 268 fQADetectors("ALL"),
269 fQATasks("ALL"),
7e963665 270 fRunQA(kTRUE),
6b150027 271 fRunGlobalQA(kTRUE),
e4a998ed 272 fSameQACycle(kFALSE),
f1c1204d 273 fInitQACalled(kFALSE),
75373542 274 fWriteQAExpertData(kTRUE),
21a3aa09 275 fRunPlaneEff(kFALSE),
276
277 fesd(NULL),
278 fhltesd(NULL),
279 fesdf(NULL),
280 ffile(NULL),
ee7c441c 281 ffileF(NULL),
21a3aa09 282 ftree(NULL),
ee7c441c 283 ftreeF(NULL),
21a3aa09 284 fhlttree(NULL),
21a3aa09 285 ftVertexer(NULL),
14dd053c 286 fIsNewRunLoader(kFALSE),
4b71572b 287 fRunAliEVE(kFALSE),
288 fChain(NULL)
596a855f 289{
290// create reconstruction object with default parameters
002c9d1b 291 gGeoManager = NULL;
b8cd5251 292
ac4a7581 293 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 294 fReconstructor[iDet] = NULL;
295 fLoader[iDet] = NULL;
296 fTracker[iDet] = NULL;
297 }
4e25ac79 298 for (Int_t iDet = 0; iDet < AliQAv1::kNDET; iDet++) {
252f8aa8 299 fQACycles[iDet] = 999999 ;
300 fQAWriteExpert[iDet] = kFALSE ;
301 }
302
e47c4c2e 303 AliPID pid;
596a855f 304}
305
306//_____________________________________________________________________________
307AliReconstruction::AliReconstruction(const AliReconstruction& rec) :
4b71572b 308 TSelector(),
2257f27e 309 fRunVertexFinder(rec.fRunVertexFinder),
a84e2607 310 fRunVertexFinderTracks(rec.fRunVertexFinderTracks),
1f46a9ae 311 fRunHLTTracking(rec.fRunHLTTracking),
e66fbafb 312 fRunMuonTracking(rec.fRunMuonTracking),
d1683eef 313 fRunV0Finder(rec.fRunV0Finder),
314 fRunCascadeFinder(rec.fRunCascadeFinder),
1d99986f 315 fStopOnError(rec.fStopOnError),
316 fWriteAlignmentData(rec.fWriteAlignmentData),
317 fWriteESDfriend(rec.fWriteESDfriend),
b647652d 318 fFillTriggerESD(rec.fFillTriggerESD),
1d99986f 319
7f68891d 320 fCleanESD(rec.fCleanESD),
a023d8d8 321 fV0DCAmax(rec.fV0DCAmax),
5e5c1aa9 322 fV0CsPmin(rec.fV0CsPmin),
7f68891d 323 fDmax(rec.fDmax),
324 fZmax(rec.fZmax),
325
1d99986f 326 fRunLocalReconstruction(rec.fRunLocalReconstruction),
e583c30d 327 fRunTracking(rec.fRunTracking),
328 fFillESD(rec.fFillESD),
7d566c20 329 fLoadCDB(rec.fLoadCDB),
48ce48d1 330 fUseTrackingErrorsForAlignment(rec.fUseTrackingErrorsForAlignment),
e583c30d 331 fGAliceFileName(rec.fGAliceFileName),
4b71572b 332 fRawInput(rec.fRawInput),
975320a0 333 fESDOutput(rec.fESDOutput),
81d57268 334 fProofOutputFileName(rec.fProofOutputFileName),
335 fProofOutputLocation(rec.fProofOutputLocation),
336 fProofOutputDataset(rec.fProofOutputDataset),
337 fProofOutputArchive(rec.fProofOutputArchive),
35042093 338 fEquipIdMap(rec.fEquipIdMap),
b26c3770 339 fFirstEvent(rec.fFirstEvent),
340 fLastEvent(rec.fLastEvent),
973388c2 341 fNumberOfEventsPerFile(rec.fNumberOfEventsPerFile),
b8cd5251 342 fOptions(),
6bae477a 343 fLoadAlignFromCDB(rec.fLoadAlignFromCDB),
344 fLoadAlignData(rec.fLoadAlignData),
2972d4eb 345 fUseHLTData(rec.fUseHLTData),
7e88424f 346 fRunInfo(NULL),
347 fEventInfo(),
522fdd91 348 fRunScalers(NULL),
52dd4a8c 349 fCTPTimeParams(NULL),
e583c30d 350
351 fRunLoader(NULL),
b649205a 352 fRawReader(NULL),
2972d4eb 353 fParentRawReader(NULL),
b8cd5251 354
4b71572b 355 fRecoParam(rec.fRecoParam),
7e88424f 356
58e8dc31 357 fSPDTrackleter(NULL),
358
f2a195c1 359 fDiamondProfileSPD(rec.fDiamondProfileSPD),
4b71572b 360 fDiamondProfile(rec.fDiamondProfile),
361 fDiamondProfileTPC(rec.fDiamondProfileTPC),
77ba28ba 362 fListOfCosmicTriggers(NULL),
f2a195c1 363
6b6e4472 364 fGRPData(NULL),
444753c6 365
6bae477a 366 fAlignObjArray(rec.fAlignObjArray),
ec92bee0 367 fCDBUri(rec.fCDBUri),
f1c1204d 368 fQARefUri(rec.fQARefUri),
7e963665 369 fSpecCDBUri(),
795e4a22 370 fInitCDBCalled(rec.fInitCDBCalled),
371 fSetRunNumberFromDataCalled(rec.fSetRunNumberFromDataCalled),
ce43afbe 372 fQADetectors(rec.fQADetectors),
373 fQATasks(rec.fQATasks),
aa3c69a9 374 fRunQA(rec.fRunQA),
375 fRunGlobalQA(rec.fRunGlobalQA),
e4a998ed 376 fSameQACycle(rec.fSameQACycle),
f1c1204d 377 fInitQACalled(rec.fInitQACalled),
75373542 378 fWriteQAExpertData(rec.fWriteQAExpertData),
21a3aa09 379 fRunPlaneEff(rec.fRunPlaneEff),
380
381 fesd(NULL),
382 fhltesd(NULL),
383 fesdf(NULL),
384 ffile(NULL),
ee7c441c 385 ffileF(NULL),
21a3aa09 386 ftree(NULL),
ee7c441c 387 ftreeF(NULL),
21a3aa09 388 fhlttree(NULL),
21a3aa09 389 ftVertexer(NULL),
14dd053c 390 fIsNewRunLoader(rec.fIsNewRunLoader),
4b71572b 391 fRunAliEVE(kFALSE),
392 fChain(NULL)
596a855f 393{
394// copy constructor
395
ec92bee0 396 for (Int_t i = 0; i < rec.fOptions.GetEntriesFast(); i++) {
efd2085e 397 if (rec.fOptions[i]) fOptions.Add(rec.fOptions[i]->Clone());
398 }
ac4a7581 399 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 400 fReconstructor[iDet] = NULL;
401 fLoader[iDet] = NULL;
402 fTracker[iDet] = NULL;
b1af1125 403 }
404
4e25ac79 405 for (Int_t iDet = 0; iDet < AliQAv1::kNDET; iDet++) {
b1af1125 406 fQACycles[iDet] = rec.fQACycles[iDet];
252f8aa8 407 fQAWriteExpert[iDet] = rec.fQAWriteExpert[iDet] ;
408 }
b1af1125 409
ec92bee0 410 for (Int_t i = 0; i < rec.fSpecCDBUri.GetEntriesFast(); i++) {
411 if (rec.fSpecCDBUri[i]) fSpecCDBUri.Add(rec.fSpecCDBUri[i]->Clone());
412 }
596a855f 413}
414
415//_____________________________________________________________________________
416AliReconstruction& AliReconstruction::operator = (const AliReconstruction& rec)
417{
418// assignment operator
4b71572b 419// Used in PROOF mode
420// Be very careful while modifing it!
421// Simple rules to follow:
422// for persistent data members - use their assignment operators
423// for non-persistent ones - do nothing or take the default values from constructor
424// TSelector members should not be touched
425 if(&rec == this) return *this;
426
4b71572b 427 fRunVertexFinder = rec.fRunVertexFinder;
428 fRunVertexFinderTracks = rec.fRunVertexFinderTracks;
429 fRunHLTTracking = rec.fRunHLTTracking;
430 fRunMuonTracking = rec.fRunMuonTracking;
431 fRunV0Finder = rec.fRunV0Finder;
432 fRunCascadeFinder = rec.fRunCascadeFinder;
433 fStopOnError = rec.fStopOnError;
434 fWriteAlignmentData = rec.fWriteAlignmentData;
435 fWriteESDfriend = rec.fWriteESDfriend;
436 fFillTriggerESD = rec.fFillTriggerESD;
437
438 fCleanESD = rec.fCleanESD;
439 fV0DCAmax = rec.fV0DCAmax;
440 fV0CsPmin = rec.fV0CsPmin;
441 fDmax = rec.fDmax;
442 fZmax = rec.fZmax;
443
444 fRunLocalReconstruction = rec.fRunLocalReconstruction;
445 fRunTracking = rec.fRunTracking;
446 fFillESD = rec.fFillESD;
7d566c20 447 fLoadCDB = rec.fLoadCDB;
4b71572b 448 fUseTrackingErrorsForAlignment = rec.fUseTrackingErrorsForAlignment;
449 fGAliceFileName = rec.fGAliceFileName;
450 fRawInput = rec.fRawInput;
975320a0 451 fESDOutput = rec.fESDOutput;
81d57268 452 fProofOutputFileName = rec.fProofOutputFileName;
453 fProofOutputLocation = rec.fProofOutputLocation;
454 fProofOutputDataset = rec.fProofOutputDataset;
455 fProofOutputArchive = rec.fProofOutputArchive;
4b71572b 456 fEquipIdMap = rec.fEquipIdMap;
457 fFirstEvent = rec.fFirstEvent;
458 fLastEvent = rec.fLastEvent;
459 fNumberOfEventsPerFile = rec.fNumberOfEventsPerFile;
460
461 for (Int_t i = 0; i < rec.fOptions.GetEntriesFast(); i++) {
462 if (rec.fOptions[i]) fOptions.Add(rec.fOptions[i]->Clone());
463 }
464
465 fLoadAlignFromCDB = rec.fLoadAlignFromCDB;
466 fLoadAlignData = rec.fLoadAlignData;
4b71572b 467 fUseHLTData = rec.fUseHLTData;
468
469 delete fRunInfo; fRunInfo = NULL;
470 if (rec.fRunInfo) fRunInfo = new AliRunInfo(*rec.fRunInfo);
471
472 fEventInfo = rec.fEventInfo;
473
522fdd91 474 delete fRunScalers; fRunScalers = NULL;
475 if (rec.fRunScalers) fRunScalers = new AliTriggerRunScalers(*rec.fRunScalers);
476
52dd4a8c 477 delete fCTPTimeParams; fCTPTimeParams = NULL;
478 if (rec.fCTPTimeParams) fCTPTimeParams = new AliCTPTimeParams(*rec.fCTPTimeParams);
479
4b71572b 480 fRunLoader = NULL;
481 fRawReader = NULL;
482 fParentRawReader = NULL;
483
484 fRecoParam = rec.fRecoParam;
485
ac4a7581 486 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
4b71572b 487 delete fReconstructor[iDet]; fReconstructor[iDet] = NULL;
488 delete fLoader[iDet]; fLoader[iDet] = NULL;
489 delete fTracker[iDet]; fTracker[iDet] = NULL;
4b71572b 490 }
b1af1125 491
4e25ac79 492 for (Int_t iDet = 0; iDet < AliQAv1::kNDET; iDet++) {
b1af1125 493 fQACycles[iDet] = rec.fQACycles[iDet];
252f8aa8 494 fQAWriteExpert[iDet] = rec.fQAWriteExpert[iDet] ;
495 }
58e8dc31 496
497 delete fSPDTrackleter; fSPDTrackleter = NULL;
252f8aa8 498
f2a195c1 499 delete fDiamondProfileSPD; fDiamondProfileSPD = NULL;
500 if (rec.fDiamondProfileSPD) fDiamondProfileSPD = new AliESDVertex(*rec.fDiamondProfileSPD);
4b71572b 501 delete fDiamondProfile; fDiamondProfile = NULL;
502 if (rec.fDiamondProfile) fDiamondProfile = new AliESDVertex(*rec.fDiamondProfile);
503 delete fDiamondProfileTPC; fDiamondProfileTPC = NULL;
504 if (rec.fDiamondProfileTPC) fDiamondProfileTPC = new AliESDVertex(*rec.fDiamondProfileTPC);
4b71572b 505
77ba28ba 506 delete fListOfCosmicTriggers; fListOfCosmicTriggers = NULL;
507 if (rec.fListOfCosmicTriggers) fListOfCosmicTriggers = (THashTable*)((rec.fListOfCosmicTriggers)->Clone());
508
4b71572b 509 delete fGRPData; fGRPData = NULL;
44e45fac 510 // if (rec.fGRPData) fGRPData = (TMap*)((rec.fGRPData)->Clone());
511 if (rec.fGRPData) fGRPData = (AliGRPObject*)((rec.fGRPData)->Clone());
4b71572b 512
513 delete fAlignObjArray; fAlignObjArray = NULL;
514
515 fCDBUri = "";
f1c1204d 516 fQARefUri = rec.fQARefUri;
4b71572b 517 fSpecCDBUri.Delete();
518 fInitCDBCalled = rec.fInitCDBCalled;
519 fSetRunNumberFromDataCalled = rec.fSetRunNumberFromDataCalled;
520 fQADetectors = rec.fQADetectors;
4b71572b 521 fQATasks = rec.fQATasks;
522 fRunQA = rec.fRunQA;
523 fRunGlobalQA = rec.fRunGlobalQA;
4b71572b 524 fSameQACycle = rec.fSameQACycle;
f1c1204d 525 fInitQACalled = rec.fInitQACalled;
75373542 526 fWriteQAExpertData = rec.fWriteQAExpertData;
4b71572b 527 fRunPlaneEff = rec.fRunPlaneEff;
528
529 fesd = NULL;
530 fhltesd = NULL;
531 fesdf = NULL;
532 ffile = NULL;
ee7c441c 533 ffileF = NULL;
4b71572b 534 ftree = NULL;
ee7c441c 535 ftreeF = NULL;
4b71572b 536 fhlttree = NULL;
537 ftVertexer = NULL;
538 fIsNewRunLoader = rec.fIsNewRunLoader;
539 fRunAliEVE = kFALSE;
540 fChain = NULL;
596a855f 541
596a855f 542 return *this;
543}
544
545//_____________________________________________________________________________
546AliReconstruction::~AliReconstruction()
547{
548// clean up
549
e583c30d 550 CleanUp();
77ba28ba 551 if (fListOfCosmicTriggers) {
552 fListOfCosmicTriggers->Delete();
553 delete fListOfCosmicTriggers;
554 }
4f3c479a 555 delete fGRPData;
522fdd91 556 delete fRunScalers;
52dd4a8c 557 delete fCTPTimeParams;
efd2085e 558 fOptions.Delete();
4b71572b 559 if (fAlignObjArray) {
560 fAlignObjArray->Delete();
561 delete fAlignObjArray;
562 }
ec92bee0 563 fSpecCDBUri.Delete();
e6d66370 564
87932dab 565 AliCodeTimer::Instance()->Print();
596a855f 566}
567
024cf675 568//_____________________________________________________________________________
f1c1204d 569void AliReconstruction::InitQA()
570{
571 //Initialize the QA and start of cycle
52dd4a8c 572 AliCodeTimerAuto("",0);
f1c1204d 573
574 if (fInitQACalled) return;
575 fInitQACalled = kTRUE;
576
634696f5 577 AliQAManager * qam = AliQAManager::QAManager(AliQAv1::kRECMODE) ;
75373542 578 if (fWriteQAExpertData)
eca4fa66 579 qam->SetWriteExpert() ;
75373542 580
eca4fa66 581 if (qam->IsDefaultStorageSet()) {
f1c1204d 582 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
583 AliWarning("Default QA reference storage has been already set !");
584 AliWarning(Form("Ignoring the default storage declared in AliReconstruction: %s",fQARefUri.Data()));
585 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
eca4fa66 586 fQARefUri = qam->GetDefaultStorage()->GetURI();
f1c1204d 587 } else {
588 if (fQARefUri.Length() > 0) {
589 AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
590 AliDebug(2, Form("Default QA reference storage is set to: %s", fQARefUri.Data()));
591 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
592 } else {
593 fQARefUri="local://$ALICE_ROOT/QAref";
594 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
595 AliWarning("Default QA refeference storage not yet set !!!!");
596 AliWarning(Form("Setting it now to: %s", fQARefUri.Data()));
597 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
598
599 }
eca4fa66 600 qam->SetDefaultStorage(fQARefUri);
f1c1204d 601 }
602
603 if (fRunQA) {
eca4fa66 604 qam->SetActiveDetectors(fQADetectors) ;
4e25ac79 605 for (Int_t det = 0 ; det < AliQAv1::kNDET ; det++) {
eca4fa66 606 qam->SetCycleLength(AliQAv1::DETECTORINDEX_t(det), fQACycles[det]) ;
607 qam->SetWriteExpert(AliQAv1::DETECTORINDEX_t(det)) ;
f1c1204d 608 }
514cb8c7 609 if (!fRawReader && !fInput && IsInTasks(AliQAv1::kRAWS))
4e25ac79 610 fQATasks.ReplaceAll(Form("%d",AliQAv1::kRAWS), "") ;
eca4fa66 611 qam->SetTasks(fQATasks) ;
612 qam->InitQADataMaker(AliCDBManager::Instance()->GetRun()) ;
f1c1204d 613 }
614 if (fRunGlobalQA) {
615 Bool_t sameCycle = kFALSE ;
eca4fa66 616 AliQADataMaker *qadm = qam->GetQADataMaker(AliQAv1::kGLOBAL);
f1c1204d 617 AliInfo(Form("Initializing the global QA data maker"));
514cb8c7 618 if (IsInTasks(AliQAv1::kRECPOINTS)) {
4e25ac79 619 qadm->StartOfCycle(AliQAv1::kRECPOINTS, AliCDBManager::Instance()->GetRun(), sameCycle) ;
620 TObjArray **arr=qadm->Init(AliQAv1::kRECPOINTS);
f1c1204d 621 AliTracker::SetResidualsArray(arr);
622 sameCycle = kTRUE ;
623 }
514cb8c7 624 if (IsInTasks(AliQAv1::kESDS)) {
4e25ac79 625 qadm->StartOfCycle(AliQAv1::kESDS, AliCDBManager::Instance()->GetRun(), sameCycle) ;
626 qadm->Init(AliQAv1::kESDS);
f1c1204d 627 }
f1c1204d 628 }
fc07289e 629 AliSysInfo::AddStamp("InitQA") ;
f1c1204d 630}
631
632//_____________________________________________________________________________
87da0921 633void AliReconstruction::MergeQA(const char *fileName)
fc07289e 634{
635 //Initialize the QA and start of cycle
52dd4a8c 636 AliCodeTimerAuto("",0) ;
eca4fa66 637 AliQAManager::QAManager()->Merge(AliCDBManager::Instance()->GetRun(),fileName) ;
fc07289e 638 AliSysInfo::AddStamp("MergeQA") ;
639}
640
641//_____________________________________________________________________________
795e4a22 642void AliReconstruction::InitCDB()
024cf675 643{
644// activate a default CDB storage
645// First check if we have any CDB storage set, because it is used
646// to retrieve the calibration and alignment constants
52dd4a8c 647 AliCodeTimerAuto("",0);
024cf675 648
795e4a22 649 if (fInitCDBCalled) return;
650 fInitCDBCalled = kTRUE;
651
024cf675 652 AliCDBManager* man = AliCDBManager::Instance();
ec92bee0 653 if (man->IsDefaultStorageSet())
024cf675 654 {
655 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
ec92bee0 656 AliWarning("Default CDB storage has been already set !");
657 AliWarning(Form("Ignoring the default storage declared in AliReconstruction: %s",fCDBUri.Data()));
024cf675 658 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
795e4a22 659 fCDBUri = man->GetDefaultStorage()->GetURI();
ec92bee0 660 }
661 else {
795e4a22 662 if (fCDBUri.Length() > 0)
663 {
664 AliDebug(2,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
665 AliDebug(2, Form("Default CDB storage is set to: %s", fCDBUri.Data()));
666 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
52dd4a8c 667 man->SetDefaultStorage(fCDBUri);
668 }
669 else if (!man->GetRaw()){
670 fCDBUri="local://$ALICE_ROOT/OCDB";
795e4a22 671 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
672 AliWarning("Default CDB storage not yet set !!!!");
673 AliWarning(Form("Setting it now to: %s", fCDBUri.Data()));
674 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
52dd4a8c 675 man->SetDefaultStorage(fCDBUri);
676 }
677 else {
678 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
679 AliWarning("Default storage will be set after setting the Run Number!!!");
680 AliWarning("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
795e4a22 681 }
ec92bee0 682 }
683
684 // Now activate the detector specific CDB storage locations
c3a7b59a 685 for (Int_t i = 0; i < fSpecCDBUri.GetEntriesFast(); i++) {
686 TObject* obj = fSpecCDBUri[i];
687 if (!obj) continue;
b8ec52f6 688 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
689 AliDebug(2, Form("Specific CDB storage for %s is set to: %s",obj->GetName(),obj->GetTitle()));
690 AliDebug(2, "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
c3a7b59a 691 man->SetSpecificStorage(obj->GetName(), obj->GetTitle());
ec92bee0 692 }
002c9d1b 693 AliSysInfo::AddStamp("InitCDB");
024cf675 694}
695
696//_____________________________________________________________________________
697void AliReconstruction::SetDefaultStorage(const char* uri) {
ec92bee0 698// Store the desired default CDB storage location
699// Activate it later within the Run() method
024cf675 700
ec92bee0 701 fCDBUri = uri;
024cf675 702
703}
704
705//_____________________________________________________________________________
f1c1204d 706void AliReconstruction::SetQARefDefaultStorage(const char* uri) {
707 // Store the desired default CDB storage location
708 // Activate it later within the Run() method
709
710 fQARefUri = uri;
4e25ac79 711 AliQAv1::SetQARefStorage(fQARefUri.Data()) ;
f1c1204d 712
713}
714//_____________________________________________________________________________
c3a7b59a 715void AliReconstruction::SetSpecificStorage(const char* calibType, const char* uri) {
ec92bee0 716// Store a detector-specific CDB storage location
717// Activate it later within the Run() method
024cf675 718
c3a7b59a 719 AliCDBPath aPath(calibType);
720 if(!aPath.IsValid()){
721 // if calibType is not wildcard but it is a valid detector, add "/*" to make it a valid path
ac4a7581 722 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
c3a7b59a 723 if(!strcmp(calibType, fgkDetectorName[iDet])) {
724 aPath.SetPath(Form("%s/*", calibType));
725 AliInfo(Form("Path for specific storage set to %s", aPath.GetPath().Data()));
726 break;
727 }
728 }
729 if(!aPath.IsValid()){
730 AliError(Form("Not a valid path or detector: %s", calibType));
731 return;
732 }
733 }
734
53dd3c3d 735// // check that calibType refers to a "valid" detector name
736// Bool_t isDetector = kFALSE;
ac4a7581 737// for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
53dd3c3d 738// TString detName = fgkDetectorName[iDet];
739// if(aPath.GetLevel0() == detName) {
740// isDetector = kTRUE;
741// break;
742// }
743// }
744//
745// if(!isDetector) {
746// AliError(Form("Not a valid detector: %s", aPath.GetLevel0().Data()));
747// return;
748// }
c3a7b59a 749
750 TObject* obj = fSpecCDBUri.FindObject(aPath.GetPath().Data());
ec92bee0 751 if (obj) fSpecCDBUri.Remove(obj);
c3a7b59a 752 fSpecCDBUri.Add(new TNamed(aPath.GetPath().Data(), uri));
024cf675 753
754}
755
6bae477a 756//_____________________________________________________________________________
795e4a22 757Bool_t AliReconstruction::SetRunNumberFromData()
6bae477a 758{
759 // The method is called in Run() in order
760 // to set a correct run number.
761 // In case of raw data reconstruction the
762 // run number is taken from the raw data header
763
795e4a22 764 if (fSetRunNumberFromDataCalled) return kTRUE;
765 fSetRunNumberFromDataCalled = kTRUE;
766
767 AliCDBManager* man = AliCDBManager::Instance();
4b71572b 768
769 if(fRawReader) {
770 if(fRawReader->NextEvent()) {
771 if(man->GetRun() > 0) {
21a3aa09 772 AliWarning("Run number is taken from raw-event header! Ignoring settings in AliCDBManager!");
4b71572b 773 }
774 man->SetRun(fRawReader->GetRunNumber());
775 fRawReader->RewindEvents();
776 }
777 else {
778 if(man->GetRun() > 0) {
779 AliWarning("No raw-data events are found ! Using settings in AliCDBManager !");
780 }
781 else {
782 AliWarning("Neither raw events nor settings in AliCDBManager are found !");
783 return kFALSE;
784 }
785 }
786 }
787 else {
788 AliRunLoader *rl = AliRunLoader::Open(fGAliceFileName.Data());
789 if (!rl) {
790 AliError(Form("No run loader found in file %s", fGAliceFileName.Data()));
6bae477a 791 return kFALSE;
792 }
ec92bee0 793 else {
4b71572b 794 rl->LoadHeader();
795 // read run number from gAlice
796 if(rl->GetHeader()) {
797 man->SetRun(rl->GetHeader()->GetRun());
798 rl->UnloadHeader();
799 delete rl;
ec92bee0 800 }
801 else {
4b71572b 802 AliError("Neither run-loader header nor RawReader objects are found !");
803 delete rl;
ec92bee0 804 return kFALSE;
805 }
4b71572b 806 }
6bae477a 807 }
795e4a22 808
809 man->Print();
810
6bae477a 811 return kTRUE;
812}
813
814//_____________________________________________________________________________
795e4a22 815void AliReconstruction::SetCDBLock() {
816 // Set CDB lock: from now on it is forbidden to reset the run number
817 // or the default storage or to activate any further storage!
818
819 AliCDBManager::Instance()->SetLock(1);
820}
821
822//_____________________________________________________________________________
6bae477a 823Bool_t AliReconstruction::MisalignGeometry(const TString& detectors)
824{
825 // Read the alignment objects from CDB.
826 // Each detector is supposed to have the
827 // alignment objects in DET/Align/Data CDB path.
828 // All the detector objects are then collected,
829 // sorted by geometry level (starting from ALIC) and
830 // then applied to the TGeo geometry.
831 // Finally an overlaps check is performed.
832
833 // Load alignment data from CDB and fill fAlignObjArray
834 if(fLoadAlignFromCDB){
6bae477a 835
25be1e5c 836 TString detStr = detectors;
837 TString loadAlObjsListOfDets = "";
838
ac4a7581 839 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
8cb26cdf 840 if(!IsSelected(fgkDetectorName[iDet], detStr)) continue;
f12d42ce 841 if(!strcmp(fgkDetectorName[iDet],"HLT")) continue;
4fbb8e9d 842
843 if(AliGeomManager::GetNalignable(fgkDetectorName[iDet]) != 0)
8cb26cdf 844 {
845 loadAlObjsListOfDets += fgkDetectorName[iDet];
846 loadAlObjsListOfDets += " ";
847 }
25be1e5c 848 } // end loop over detectors
4fbb8e9d 849
850 if(AliGeomManager::GetNalignable("GRP") != 0)
8cb26cdf 851 loadAlObjsListOfDets.Prepend("GRP "); //add alignment objects for non-sensitive modules
98e303d9 852 AliGeomManager::ApplyAlignObjsFromCDB(loadAlObjsListOfDets.Data());
e30a9b4d 853 AliCDBManager::Instance()->UnloadFromCache("*/Align/*");
25be1e5c 854 }else{
855 // Check if the array with alignment objects was
856 // provided by the user. If yes, apply the objects
857 // to the present TGeo geometry
858 if (fAlignObjArray) {
859 if (gGeoManager && gGeoManager->IsClosed()) {
98e303d9 860 if (AliGeomManager::ApplyAlignObjsToGeom(*fAlignObjArray) == kFALSE) {
25be1e5c 861 AliError("The misalignment of one or more volumes failed!"
862 "Compare the list of simulated detectors and the list of detector alignment data!");
863 return kFALSE;
864 }
865 }
866 else {
867 AliError("Can't apply the misalignment! gGeoManager doesn't exist or it is still opened!");
6bae477a 868 return kFALSE;
869 }
870 }
6bae477a 871 }
25be1e5c 872
4b71572b 873 if (fAlignObjArray) {
874 fAlignObjArray->Delete();
875 delete fAlignObjArray; fAlignObjArray=NULL;
876 }
a03b0371 877
6bae477a 878 return kTRUE;
879}
596a855f 880
881//_____________________________________________________________________________
882void AliReconstruction::SetGAliceFile(const char* fileName)
883{
884// set the name of the galice file
885
886 fGAliceFileName = fileName;
887}
888
efd2085e 889//_____________________________________________________________________________
b58710ef 890void AliReconstruction::SetInput(const char* input)
21a3aa09 891{
b58710ef 892 // In case the input string starts with 'mem://', we run in an online mode
893 // and AliRawReaderDateOnline object is created. In all other cases a raw-data
894 // file is assumed. One can give as an input:
895 // mem://: - events taken from DAQ monitoring libs online
896 // or
897 // mem://<filename> - emulation of the above mode (via DATE monitoring libs)
4b71572b 898 if (input) fRawInput = input;
21a3aa09 899}
900
901//_____________________________________________________________________________
975320a0 902void AliReconstruction::SetOutput(const char* output)
903{
904 // Set the output ESD filename
905 // 'output' is a normalt ROOT url
906 // The method is used in case of raw-data reco with PROOF
81d57268 907 if (output) fESDOutput = output;
975320a0 908}
909
910//_____________________________________________________________________________
efd2085e 911void AliReconstruction::SetOption(const char* detector, const char* option)
912{
913// set options for the reconstruction of a detector
914
915 TObject* obj = fOptions.FindObject(detector);
916 if (obj) fOptions.Remove(obj);
917 fOptions.Add(new TNamed(detector, option));
918}
919
596a855f 920//_____________________________________________________________________________
7e88424f 921void AliReconstruction::SetRecoParam(const char* detector, AliDetectorRecoParam *par)
922{
923 // Set custom reconstruction parameters for a given detector
924 // Single set of parameters for all the events
a00021a7 925
926 // First check if the reco-params are global
927 if(!strcmp(detector, "GRP")) {
928 par->SetAsDefault();
ac4a7581 929 fRecoParam.AddDetRecoParam(kNDetectors,par);
a00021a7 930 return;
931 }
932
ac4a7581 933 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
7e88424f 934 if(!strcmp(detector, fgkDetectorName[iDet])) {
935 par->SetAsDefault();
936 fRecoParam.AddDetRecoParam(iDet,par);
937 break;
938 }
939 }
940
941}
942
943//_____________________________________________________________________________
61807e09 944Bool_t AliReconstruction::InitGRP() {
945 //------------------------------------
946 // Initialization of the GRP entry
947 //------------------------------------
948 AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/GRP/Data");
949
4b71572b 950 if (entry) {
44e45fac 951
952 TMap* m = dynamic_cast<TMap*>(entry->GetObject()); // old GRP entry
953
954 if (m) {
955 AliInfo("Found a TMap in GRP/GRP/Data, converting it into an AliGRPObject");
956 m->Print();
957 fGRPData = new AliGRPObject();
958 fGRPData->ReadValuesFromMap(m);
959 }
960
961 else {
962 AliInfo("Found an AliGRPObject in GRP/GRP/Data, reading it");
963 fGRPData = dynamic_cast<AliGRPObject*>(entry->GetObject()); // new GRP entry
964 entry->SetOwner(0);
965 }
966
f168abba 967 // FIX ME: The unloading of GRP entry is temporarily disabled
968 // because ZDC and VZERO are using it in order to initialize
969 // their reconstructor objects. In the future one has to think
970 // of propagating AliRunInfo to the reconstructors.
971 // AliCDBManager::Instance()->UnloadFromCache("GRP/GRP/Data");
4b71572b 972 }
61807e09 973
974 if (!fGRPData) {
975 AliError("No GRP entry found in OCDB!");
976 return kFALSE;
977 }
978
44e45fac 979 TString lhcState = fGRPData->GetLHCState();
980 if (lhcState==AliGRPObject::GetInvalidString()) {
7e88424f 981 AliError("GRP/GRP/Data entry: missing value for the LHC state ! Using UNKNOWN");
44e45fac 982 lhcState = "UNKNOWN";
7e88424f 983 }
984
44e45fac 985 TString beamType = fGRPData->GetBeamType();
986 if (beamType==AliGRPObject::GetInvalidString()) {
7e88424f 987 AliError("GRP/GRP/Data entry: missing value for the beam type ! Using UNKNOWN");
44e45fac 988 beamType = "UNKNOWN";
7e88424f 989 }
990
44e45fac 991 Float_t beamEnergy = fGRPData->GetBeamEnergy();
992 if (beamEnergy==AliGRPObject::GetInvalidFloat()) {
7e88424f 993 AliError("GRP/GRP/Data entry: missing value for the beam energy ! Using 0");
44e45fac 994 beamEnergy = 0;
7e88424f 995 }
996
44e45fac 997 TString runType = fGRPData->GetRunType();
998 if (runType==AliGRPObject::GetInvalidString()) {
7e88424f 999 AliError("GRP/GRP/Data entry: missing value for the run type ! Using UNKNOWN");
44e45fac 1000 runType = "UNKNOWN";
7e88424f 1001 }
1002
44e45fac 1003 Int_t activeDetectors = fGRPData->GetDetectorMask();
a5dede5c 1004 if (activeDetectors==AliGRPObject::GetInvalidUInt()) {
7e88424f 1005 AliError("GRP/GRP/Data entry: missing value for the detector mask ! Using 1074790399");
44e45fac 1006 activeDetectors = 1074790399;
7e88424f 1007 }
1008
44e45fac 1009 fRunInfo = new AliRunInfo(lhcState, beamType, beamEnergy, runType, activeDetectors);
44e45fac 1010 fRunInfo->Dump();
1011
7e88424f 1012
1013 // Process the list of active detectors
44e45fac 1014 if (activeDetectors) {
1015 UInt_t detMask = activeDetectors;
7a317bef 1016 fRunLocalReconstruction = MatchDetectorList(fRunLocalReconstruction,detMask);
1017 fRunTracking = MatchDetectorList(fRunTracking,detMask);
1018 fFillESD = MatchDetectorList(fFillESD,detMask);
1019 fQADetectors = MatchDetectorList(fQADetectors,detMask);
7d566c20 1020 fLoadCDB.Form("%s %s %s %s",
1021 fRunLocalReconstruction.Data(),
1022 fRunTracking.Data(),
1023 fFillESD.Data(),
1024 fQADetectors.Data());
7d566c20 1025 fLoadCDB = MatchDetectorList(fLoadCDB,detMask);
d0f84c42 1026 if (!((detMask >> AliDAQ::DetectorID("ITSSPD")) & 0x1) &&
1027 !((detMask >> AliDAQ::DetectorID("ITSSDD")) & 0x1) &&
1028 !((detMask >> AliDAQ::DetectorID("ITSSSD")) & 0x1) ) {
91e2f025 1029 // switch off the vertexer
d0f84c42 1030 AliInfo("SPD,SDD,SSD is not in the list of active detectors. Vertexer switched off.");
91e2f025 1031 fRunVertexFinder = kFALSE;
1032 }
304d7f0b 1033 if (!((detMask >> AliDAQ::DetectorID("TRG")) & 0x1)) {
1034 // switch off the reading of CTP raw-data payload
1035 if (fFillTriggerESD) {
1036 AliInfo("CTP is not in the list of active detectors. CTP data reading switched off.");
1037 fFillTriggerESD = kFALSE;
1038 }
1039 }
7e88424f 1040 }
1041
1042 AliInfo("===================================================================================");
1043 AliInfo(Form("Running local reconstruction for detectors: %s",fRunLocalReconstruction.Data()));
1044 AliInfo(Form("Running tracking for detectors: %s",fRunTracking.Data()));
1045 AliInfo(Form("Filling ESD for detectors: %s",fFillESD.Data()));
ca13fb87 1046 AliInfo(Form("Quality assurance is active for detectors: %s",fQADetectors.Data()));
7d566c20 1047 AliInfo(Form("CDB and reconstruction parameters are loaded for detectors: %s",fLoadCDB.Data()));
7e88424f 1048 AliInfo("===================================================================================");
61807e09 1049
1050 //*** Dealing with the magnetic field map
5a004fb4 1051 if ( TGeoGlobalMagField::Instance()->IsLocked() ) {
1052 if (TGeoGlobalMagField::Instance()->GetField()->TestBit(AliMagF::kOverrideGRP)) {
1053 AliInfo("ExpertMode!!! GRP information will be ignored !");
1054 AliInfo("ExpertMode!!! Running with the externally locked B field !");
1055 }
1056 else {
1057 AliInfo("Destroying existing B field instance!");
1058 delete TGeoGlobalMagField::Instance();
1059 }
1060 }
1061 if ( !TGeoGlobalMagField::Instance()->IsLocked() ) {
61807e09 1062 // Construct the field map out of the information retrieved from GRP.
856024f0 1063 Bool_t ok = kTRUE;
61807e09 1064 // L3
44e45fac 1065 Float_t l3Current = fGRPData->GetL3Current((AliGRPObject::Stats)0);
1066 if (l3Current == AliGRPObject::GetInvalidFloat()) {
1067 AliError("GRP/GRP/Data entry: missing value for the L3 current !");
1068 ok = kFALSE;
1069 }
f7a1cc68 1070
44e45fac 1071 Char_t l3Polarity = fGRPData->GetL3Polarity();
1072 if (l3Polarity == AliGRPObject::GetInvalidChar()) {
1073 AliError("GRP/GRP/Data entry: missing value for the L3 polarity !");
1074 ok = kFALSE;
1075 }
1076
1077 // Dipole
1078 Float_t diCurrent = fGRPData->GetDipoleCurrent((AliGRPObject::Stats)0);
1079 if (diCurrent == AliGRPObject::GetInvalidFloat()) {
1080 AliError("GRP/GRP/Data entry: missing value for the dipole current !");
1081 ok = kFALSE;
1082 }
1083
1084 Char_t diPolarity = fGRPData->GetDipolePolarity();
1085 if (diPolarity == AliGRPObject::GetInvalidChar()) {
1086 AliError("GRP/GRP/Data entry: missing value for the dipole polarity !");
1087 ok = kFALSE;
1088 }
1089
99c7d495 1090 // read special bits for the polarity convention and map type
1091 Int_t polConvention = fGRPData->IsPolarityConventionLHC() ? AliMagF::kConvLHC : AliMagF::kConvDCS2008;
1092 Bool_t uniformB = fGRPData->IsUniformBMap();
61807e09 1093
856024f0 1094 if (ok) {
33fe5eb1 1095 AliMagF* fld = AliMagF::CreateFieldMap(TMath::Abs(l3Current) * (l3Polarity ? -1:1),
1096 TMath::Abs(diCurrent) * (diPolarity ? -1:1),
1097 polConvention,uniformB,beamEnergy, beamType.Data());
1098 if (fld) {
1099 TGeoGlobalMagField::Instance()->SetField( fld );
1100 TGeoGlobalMagField::Instance()->Lock();
1101 AliInfo("Running with the B field constructed out of GRP !");
1102 }
1103 else AliFatal("Failed to create a B field map !");
856024f0 1104 }
f7a1cc68 1105 else AliFatal("B field is neither set nor constructed from GRP ! Exitig...");
61807e09 1106 }
f7a1cc68 1107
f2a195c1 1108 //*** Get the diamond profiles from OCDB
1109 entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertexSPD");
1110 if (entry) {
1111 fDiamondProfileSPD = dynamic_cast<AliESDVertex*> (entry->GetObject());
1112 } else {
1113 AliError("No SPD diamond profile found in OCDB!");
1114 }
1115
61807e09 1116 entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertex");
1117 if (entry) {
f2a195c1 1118 fDiamondProfile = dynamic_cast<AliESDVertex*> (entry->GetObject());
61807e09 1119 } else {
1120 AliError("No diamond profile found in OCDB!");
1121 }
1122
1123 entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertexTPC");
1124 if (entry) {
f2a195c1 1125 fDiamondProfileTPC = dynamic_cast<AliESDVertex*> (entry->GetObject());
61807e09 1126 } else {
f2a195c1 1127 AliError("No TPC diamond profile found in OCDB!");
61807e09 1128 }
1129
77ba28ba 1130 entry = AliCDBManager::Instance()->Get("GRP/Calib/CosmicTriggers");
1131 if (entry) {
1132 fListOfCosmicTriggers = dynamic_cast<THashTable*>(entry->GetObject());
1133 entry->SetOwner(0);
1134 AliCDBManager::Instance()->UnloadFromCache("GRP/Calib/CosmicTriggers");
1135 }
1136
1137 if (!fListOfCosmicTriggers) {
1138 AliWarning("Can not get list of cosmic triggers from OCDB! Cosmic event specie will be effectively disabled!");
1139 }
1140
61807e09 1141 return kTRUE;
1142}
1143
1144//_____________________________________________________________________________
4b71572b 1145Bool_t AliReconstruction::LoadCDB()
1146{
0a035be5 1147 // Load CDB entries for all active detectors.
1148 // By default we load all the entries in <det>/Calib
1149 // folder.
1150
52dd4a8c 1151 AliCodeTimerAuto("",0);
4b71572b 1152
1153 AliCDBManager::Instance()->Get("GRP/CTP/Config");
1154
7d566c20 1155 TString detStr = fLoadCDB;
ac4a7581 1156 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
4b71572b 1157 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
1158 AliCDBManager::Instance()->GetAll(Form("%s/Calib/*",fgkDetectorName[iDet]));
1159 }
e122592f 1160
1161 // Temporary fix - one has to define the correct policy in order
1162 // to load the trigger OCDB entries only for the detectors that
1163 // in the trigger or that are needed in order to put correct
1164 // information in ESD
1165 AliCDBManager::Instance()->GetAll("TRIGGER/*/*");
1166
4b71572b 1167 return kTRUE;
1168}
522fdd91 1169//_____________________________________________________________________________
1170Bool_t AliReconstruction::LoadTriggerScalersCDB()
1171{
0a035be5 1172 // Load CTP scalers from OCDB.
1173 // The scalers are checked for consistency.
1174
52dd4a8c 1175 AliCodeTimerAuto("",0);
522fdd91 1176
1177 AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/CTP/Scalers");
4b71572b 1178
522fdd91 1179 if (entry) {
1180
1181 AliInfo("Found an AliTriggerRunScalers in GRP/CTP/Scalers, reading it");
1182 fRunScalers = dynamic_cast<AliTriggerRunScalers*> (entry->GetObject());
1183 entry->SetOwner(0);
82ebedd6 1184 if (fRunScalers->CorrectScalersOverflow() == 0) AliInfo("32bit Trigger counters corrected for overflow");
1185
522fdd91 1186 }
1187 return kTRUE;
1188}
4b71572b 1189//_____________________________________________________________________________
52dd4a8c 1190Bool_t AliReconstruction::LoadCTPTimeParamsCDB()
1191{
0a035be5 1192 // Load CTP timing information (alignment)
1193 // from OCDB.
1194
52dd4a8c 1195 AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/CTP/CTPtiming");
1196
1197 if (entry) {
1198
1199 AliInfo("Found an AliCTPTimeParams in GRP/CTP/CTPtiming, reading it");
1200 fCTPTimeParams = dynamic_cast<AliCTPTimeParams*> (entry->GetObject());
1201 entry->SetOwner(0);
1202 return kTRUE;
1203 }
1204
1205 return kFALSE;
1206}
1207//_____________________________________________________________________________
21a3aa09 1208Bool_t AliReconstruction::Run(const char* input)
596a855f 1209{
21a3aa09 1210 // Run Run Run
52dd4a8c 1211 AliCodeTimerAuto("",0);
596a855f 1212
4b71572b 1213 InitRun(input);
1214 if (GetAbort() != TSelector::kContinue) return kFALSE;
1215
1216 TChain *chain = NULL;
1217 if (fRawReader && (chain = fRawReader->GetChain())) {
2a284cef 1218 Long64_t nEntries = (fLastEvent < 0) ? (TChain::kBigNumber) : (fLastEvent - fFirstEvent + 1);
4b71572b 1219 // Proof mode
1220 if (gProof) {
52dd4a8c 1221 // Temporary fix for long raw-data runs (until socket timeout handling in PROOF is revised)
1222 gProof->Exec("gEnv->SetValue(\"Proof.SocketActivityTimeout\",-1)", kTRUE);
6c6f2624 1223
1224 if (gGrid)
1225 gProof->Exec("TGrid::Connect(\"alien://\")",kTRUE);
1226
1227 TMessage::EnableSchemaEvolutionForAll(kTRUE);
1228 gProof->Exec("TMessage::EnableSchemaEvolutionForAll(kTRUE)",kTRUE);
1229
4b71572b 1230 gProof->AddInput(this);
81d57268 1231
1232 if (!ParseOutput()) return kFALSE;
1233
975320a0 1234 gProof->SetParameter("PROOF_MaxSlavesPerNode", 9999);
4b71572b 1235 chain->SetProof();
2a284cef 1236 chain->Process("AliReconstruction","",nEntries,fFirstEvent);
4b71572b 1237 }
1238 else {
2a284cef 1239 chain->Process(this,"",nEntries,fFirstEvent);
4b71572b 1240 }
1241 }
1242 else {
1243 Begin(NULL);
1244 if (GetAbort() != TSelector::kContinue) return kFALSE;
1245 SlaveBegin(NULL);
1246 if (GetAbort() != TSelector::kContinue) return kFALSE;
1247 //******* The loop over events
44e45fac 1248 AliInfo("Starting looping over events");
4b71572b 1249 Int_t iEvent = 0;
1250 while ((iEvent < fRunLoader->GetNumberOfEvents()) ||
1251 (fRawReader && fRawReader->NextEvent())) {
1252 if (!ProcessEvent(iEvent)) {
930e6e3e 1253 Abort("ProcessEvent",TSelector::kAbortFile);
1254 return kFALSE;
4b71572b 1255 }
1256 iEvent++;
1257 }
1258 SlaveTerminate();
1259 if (GetAbort() != TSelector::kContinue) return kFALSE;
1260 Terminate();
1261 if (GetAbort() != TSelector::kContinue) return kFALSE;
21a3aa09 1262 }
21a3aa09 1263
1264 return kTRUE;
1265}
1266
1267//_____________________________________________________________________________
4b71572b 1268void AliReconstruction::InitRawReader(const char* input)
21a3aa09 1269{
4b71572b 1270 // Init raw-reader and
21a3aa09 1271 // set the input in case of raw data
0a035be5 1272
1273 AliCodeTimerAuto("",0);
1274
4b71572b 1275 if (input) fRawInput = input;
1276 fRawReader = AliRawReader::Create(fRawInput.Data());
c35f192e 1277 if (!fRawReader) {
1278 if (fRawInput.IsNull()) {
1279 AliInfo("Reconstruction will run over digits");
1280 }
1281 else {
1282 AliFatal("Can not create raw-data reader ! Exiting...");
1283 }
1284 }
e0027792 1285
35042093 1286 if (!fEquipIdMap.IsNull() && fRawReader)
1287 fRawReader->LoadEquipmentIdsMap(fEquipIdMap);
1288
cd0b062e 1289 if (!fUseHLTData.IsNull()) {
1290 // create the RawReaderHLT which performs redirection of HLT input data for
1291 // the specified detectors
1292 AliRawReader* pRawReader=AliRawHLTManager::CreateRawReaderHLT(fRawReader, fUseHLTData.Data());
1293 if (pRawReader) {
1294 fParentRawReader=fRawReader;
1295 fRawReader=pRawReader;
1296 } else {
1297 AliError(Form("can not create Raw Reader for HLT input %s", fUseHLTData.Data()));
1298 }
1299 }
4b71572b 1300 AliSysInfo::AddStamp("CreateRawReader");
1301}
cd0b062e 1302
4b71572b 1303//_____________________________________________________________________________
1304void AliReconstruction::InitRun(const char* input)
1305{
1306 // Initialization of raw-reader,
1307 // run number, CDB etc.
52dd4a8c 1308 AliCodeTimerAuto("",0);
4b71572b 1309 AliSysInfo::AddStamp("Start");
1310
002c9d1b 1311 // Initialize raw-reader if any
4b71572b 1312 InitRawReader(input);
596a855f 1313
ec92bee0 1314 // Initialize the CDB storage
795e4a22 1315 InitCDB();
ec92bee0 1316
6bae477a 1317 // Set run number in CDBManager (if it is not already set by the user)
4b71572b 1318 if (!SetRunNumberFromData()) {
1319 Abort("SetRunNumberFromData", TSelector::kAbortProcess);
1320 return;
1321 }
002c9d1b 1322
795e4a22 1323 // Set CDB lock: from now on it is forbidden to reset the run number
1324 // or the default storage or to activate any further storage!
1325 SetCDBLock();
1326
002c9d1b 1327}
1328
1329//_____________________________________________________________________________
1330void AliReconstruction::Begin(TTree *)
1331{
1332 // Initialize AlReconstruction before
1333 // going into the event loop
1334 // Should follow the TSelector convention
1335 // i.e. initialize only the object on the client side
52dd4a8c 1336 AliCodeTimerAuto("",0);
db4aeca1 1337
b4fef0d2 1338 AliReconstruction *reco = NULL;
002c9d1b 1339 if (fInput) {
8fb1f1dc 1340 if ((reco = (AliReconstruction*)fInput->FindObject("AliReconstruction"))) {
002c9d1b 1341 *this = *reco;
002c9d1b 1342 }
1343 AliSysInfo::AddStamp("ReadInputInBegin");
1344 }
1345
6bae477a 1346 // Import ideal TGeo geometry and apply misalignment
1347 if (!gGeoManager) {
1348 TString geom(gSystem->DirName(fGAliceFileName));
1349 geom += "/geometry.root";
98e303d9 1350 AliGeomManager::LoadGeometry(geom.Data());
4b71572b 1351 if (!gGeoManager) {
1352 Abort("LoadGeometry", TSelector::kAbortProcess);
1353 return;
1354 }
281270d3 1355 AliSysInfo::AddStamp("LoadGeom");
ff5970a3 1356 TString detsToCheck=fRunLocalReconstruction;
4b71572b 1357 if(!AliGeomManager::CheckSymNamesLUT(detsToCheck.Data())) {
1358 Abort("CheckSymNamesLUT", TSelector::kAbortProcess);
1359 return;
1360 }
281270d3 1361 AliSysInfo::AddStamp("CheckGeom");
6bae477a 1362 }
8e245d15 1363
4b71572b 1364 if (!MisalignGeometry(fLoadAlignData)) {
1365 Abort("MisalignGeometry", TSelector::kAbortProcess);
1366 return;
1367 }
1368 AliCDBManager::Instance()->UnloadFromCache("GRP/Geometry/Data");
281270d3 1369 AliSysInfo::AddStamp("MisalignGeom");
6bae477a 1370
8cb26cdf 1371 if (!InitGRP()) {
1372 Abort("InitGRP", TSelector::kAbortProcess);
1373 return;
1374 }
1375 AliSysInfo::AddStamp("InitGRP");
1376
4b71572b 1377 if (!LoadCDB()) {
1378 Abort("LoadCDB", TSelector::kAbortProcess);
1379 return;
1380 }
1381 AliSysInfo::AddStamp("LoadCDB");
325aa76f 1382
a6dd87ad 1383 if (!LoadTriggerScalersCDB()) {
1384 Abort("LoadTriggerScalersCDB", TSelector::kAbortProcess);
1385 return;
1386 }
1387 AliSysInfo::AddStamp("LoadTriggerScalersCDB");
1388
52dd4a8c 1389 if (!LoadCTPTimeParamsCDB()) {
1390 Abort("LoadCTPTimeParamsCDB", TSelector::kAbortProcess);
1391 return;
1392 }
1393 AliSysInfo::AddStamp("LoadCTPTimeParamsCDB");
a6dd87ad 1394
7e88424f 1395 // Read the reconstruction parameters from OCDB
1396 if (!InitRecoParams()) {
002c9d1b 1397 AliWarning("Not all detectors have correct RecoParam objects initialized");
4b71572b 1398 }
1399 AliSysInfo::AddStamp("InitRecoParams");
1400
fd2b273c 1401 if (fInput && gProof) {
b4fef0d2 1402 if (reco) *reco = *this;
fd2b273c 1403
45a6dee6 1404 gGeoManager->SetName("Geometry");
fd2b273c 1405 gProof->AddInputData(gGeoManager,kTRUE);
002c9d1b 1406 gGeoManager = NULL;
fd2b273c 1407 gProof->AddInputData(const_cast<TMap*>(AliCDBManager::Instance()->GetEntryCache()),kTRUE);
002c9d1b 1408 fInput->Add(new TParameter<Int_t>("RunNumber",AliCDBManager::Instance()->GetRun()));
f7a1cc68 1409 AliMagF *magFieldMap = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
ead7b391 1410 magFieldMap->SetName("MagneticFieldMap");
fd2b273c 1411 gProof->AddInputData(magFieldMap,kTRUE);
002c9d1b 1412 }
4b71572b 1413
1414}
1415
1416//_____________________________________________________________________________
1417void AliReconstruction::SlaveBegin(TTree*)
1418{
1419 // Initialization related to run-loader,
1420 // vertexer, trackers, recontructors
1421 // In proof mode it is executed on the slave
52dd4a8c 1422 AliCodeTimerAuto("",0);
4b71572b 1423
db4aeca1 1424 TProofOutputFile *outProofFile = NULL;
81d57268 1425 if (fInput) {
1426 if (AliDebugLevel() > 0) fInput->Print();
ef5e52f6 1427 if (AliDebugLevel() > 10) fInput->Dump();
db4aeca1 1428 if (AliReconstruction *reco = (AliReconstruction*)fInput->FindObject("AliReconstruction")) {
1429 *this = *reco;
1430 }
4b71572b 1431 if (TGeoManager *tgeo = (TGeoManager*)fInput->FindObject("Geometry")) {
1432 gGeoManager = tgeo;
1433 AliGeomManager::SetGeometry(tgeo);
7e88424f 1434 }
281270d3 1435 if (TMap *entryCache = (TMap*)fInput->FindObject("CDBEntryCache")) {
1436 Int_t runNumber = -1;
1437 if (TProof::GetParameter(fInput,"RunNumber",runNumber) == 0) {
1438 AliCDBManager *man = AliCDBManager::Instance(entryCache,runNumber);
1439 man->SetCacheFlag(kTRUE);
1440 man->SetLock(kTRUE);
1441 man->Print();
1442 }
4b71572b 1443 }
ead7b391 1444 if (AliMagF *map = (AliMagF*)fInput->FindObject("MagneticFieldMap")) {
ef5e52f6 1445 AliMagF *newMap = new AliMagF(*map);
1446 if (!newMap->LoadParameterization()) {
1447 Abort("AliMagF::LoadParameterization", TSelector::kAbortProcess);
1448 return;
1449 }
1450 TGeoGlobalMagField::Instance()->SetField(newMap);
1451 TGeoGlobalMagField::Instance()->Lock();
4b71572b 1452 }
81d57268 1453 if (TNamed *outputFileName = (TNamed*)fInput->FindObject("PROOF_OUTPUTFILE"))
1454 fProofOutputFileName = outputFileName->GetTitle();
1455 if (TNamed *outputLocation = (TNamed*)fInput->FindObject("PROOF_OUTPUTFILE_LOCATION"))
1456 fProofOutputLocation = outputLocation->GetTitle();
1457 if (fInput->FindObject("PROOF_OUTPUTFILE_DATASET"))
1458 fProofOutputDataset = kTRUE;
1459 if (TNamed *archiveList = (TNamed*)fInput->FindObject("PROOF_OUTPUTFILE_ARCHIVE"))
1460 fProofOutputArchive = archiveList->GetTitle();
1461 if (!fProofOutputFileName.IsNull() &&
1462 !fProofOutputLocation.IsNull() &&
1463 fProofOutputArchive.IsNull()) {
1464 if (!fProofOutputDataset) {
1465 outProofFile = new TProofOutputFile(fProofOutputFileName.Data(),"M");
1466 outProofFile->SetOutputFileName(Form("%s%s",fProofOutputLocation.Data(),fProofOutputFileName.Data()));
1467 }
1468 else {
1469 outProofFile = new TProofOutputFile(fProofOutputFileName.Data(),"DROV",fProofOutputLocation.Data());
1470 }
1471 if (AliDebugLevel() > 0) outProofFile->Dump();
db4aeca1 1472 fOutput->Add(outProofFile);
4b71572b 1473 }
002c9d1b 1474 AliSysInfo::AddStamp("ReadInputInSlaveBegin");
7e88424f 1475 }
325aa76f 1476
4b71572b 1477 // get the run loader
1478 if (!InitRunLoader()) {
1479 Abort("InitRunLoader", TSelector::kAbortProcess);
1480 return;
1481 }
1482 AliSysInfo::AddStamp("LoadLoader");
1483
325aa76f 1484 ftVertexer = new AliVertexerTracks(AliTracker::GetBz());
2257f27e 1485
f08fc9f5 1486 // get trackers
b8cd5251 1487 if (!fRunTracking.IsNull() && !CreateTrackers(fRunTracking)) {
4b71572b 1488 Abort("CreateTrackers", TSelector::kAbortProcess);
1489 return;
1490 }
1491 AliSysInfo::AddStamp("CreateTrackers");
24f7a148 1492
36711aa4 1493 // create the ESD output file and tree
db4aeca1 1494 if (!outProofFile) {
1495 ffile = TFile::Open("AliESDs.root", "RECREATE");
1496 ffile->SetCompressionLevel(2);
1497 if (!ffile->IsOpen()) {
1498 Abort("OpenESDFile", TSelector::kAbortProcess);
1499 return;
1500 }
1501 }
1502 else {
81d57268 1503 AliInfo(Form("Opening output PROOF file: %s/%s",
1504 outProofFile->GetDir(), outProofFile->GetFileName()));
db4aeca1 1505 if (!(ffile = outProofFile->OpenFile("RECREATE"))) {
1506 Abort(Form("Problems opening output PROOF file: %s/%s",
1507 outProofFile->GetDir(), outProofFile->GetFileName()),
1508 TSelector::kAbortProcess);
1509 return;
1510 }
596a855f 1511 }
46698ae4 1512
21a3aa09 1513 ftree = new TTree("esdTree", "Tree with ESD objects");
1514 fesd = new AliESDEvent();
1515 fesd->CreateStdContent();
8b12d288 1516 // add a so far non-std object to the ESD, this will
1517 // become part of the std content
1518 fesd->AddObject(new AliESDHLTDecision);
1519
32ba9c61 1520 fesd->WriteToTree(ftree);
1d99986f 1521 if (fWriteESDfriend) {
ee7c441c 1522 ffileF = TFile::Open("AliESDfriends.root", "RECREATE");
03a8000a 1523 ftreeF = new TTree("esdFriendTree", "Tree with ESD Friend objects");
ee7c441c 1524 fesdf = new AliESDfriend();
03a8000a 1525 ftreeF->Branch("ESDfriend.","AliESDfriend", &fesdf);
21a3aa09 1526 fesd->AddObject(fesdf);
ee7c441c 1527 ffile->cd();
1d99986f 1528 }
f9604a22 1529 ftree->GetUserInfo()->Add(fesd);
1530
1531 fhlttree = new TTree("HLTesdTree", "Tree with HLT ESD objects");
1532 fhltesd = new AliESDEvent();
1533 fhltesd->CreateStdContent();
7491aea6 1534 // read the ESD template from CDB
1535 // HLT is allowed to put non-std content to its ESD, the non-std
1536 // objects need to be created before invocation of WriteToTree in
1537 // order to create all branches. Initialization is done from an
1538 // ESD layout template in CDB
1539 AliCDBManager* man = AliCDBManager::Instance();
1540 AliCDBPath hltESDConfigPath("HLT/ConfigHLT/esdLayout");
1541 AliCDBEntry* hltESDConfig=NULL;
1542 if (man->GetId(hltESDConfigPath)!=NULL &&
1543 (hltESDConfig=man->Get(hltESDConfigPath))!=NULL) {
1544 AliESDEvent* pESDLayout=dynamic_cast<AliESDEvent*>(hltESDConfig->GetObject());
1545 if (pESDLayout) {
1546 // init all internal variables from the list of objects
1547 pESDLayout->GetStdContent();
1548
1549 // copy content and create non-std objects
1550 *fhltesd=*pESDLayout;
1551 fhltesd->Reset();
1552 } else {
1553 AliError(Form("error setting hltEsd layout from %s: invalid object type",
1554 hltESDConfigPath.GetPath().Data()));
1555 }
1556 }
1557
f9604a22 1558 fhltesd->WriteToTree(fhlttree);
1559 fhlttree->GetUserInfo()->Add(fhltesd);
5728d3d5 1560
ac4a7581 1561 ProcInfo_t procInfo;
1562 gSystem->GetProcInfo(&procInfo);
1563 AliInfo(Form("Current memory usage %d %d", procInfo.fMemResident, procInfo.fMemVirtual));
8661738e 1564
325aa76f 1565 //QA
930e6e3e 1566 //Initialize the QA and start of cycle
f1c1204d 1567 if (fRunQA || fRunGlobalQA)
1568 InitQA() ;
aa3c69a9 1569
7167ae53 1570 //Initialize the Plane Efficiency framework
1571 if (fRunPlaneEff && !InitPlaneEff()) {
4b71572b 1572 Abort("InitPlaneEff", TSelector::kAbortProcess);
1573 return;
7167ae53 1574 }
aa3c69a9 1575
14dd053c 1576 if (strcmp(gProgName,"alieve") == 0)
1577 fRunAliEVE = InitAliEVE();
1578
4b71572b 1579 return;
21a3aa09 1580}
1581
1582//_____________________________________________________________________________
4b71572b 1583Bool_t AliReconstruction::Process(Long64_t entry)
1584{
1585 // run the reconstruction over a single entry
1586 // from the chain with raw data
52dd4a8c 1587 AliCodeTimerAuto("",0);
4b71572b 1588
1589 TTree *currTree = fChain->GetTree();
33314186 1590 AliRawVEvent *event = NULL;
4b71572b 1591 currTree->SetBranchAddress("rawevent",&event);
1592 currTree->GetEntry(entry);
1593 fRawReader = new AliRawReaderRoot(event);
1594 fStatus = ProcessEvent(fRunLoader->GetNumberOfEvents());
1595 delete fRawReader;
1596 fRawReader = NULL;
1597 delete event;
1598
1599 return fStatus;
1600}
1601
1602//_____________________________________________________________________________
1603void AliReconstruction::Init(TTree *tree)
1604{
0a035be5 1605 // Implementation of TSelector::Init()
1606 // method
4b71572b 1607 if (tree == 0) {
1608 AliError("The input tree is not found!");
1609 return;
1610 }
1611 fChain = tree;
1612}
1613
1614//_____________________________________________________________________________
1615Bool_t AliReconstruction::ProcessEvent(Int_t iEvent)
21a3aa09 1616{
1617 // run the reconstruction over a single event
1618 // The event loop is steered in Run method
1619
10d100d4 1620
470f88b0 1621 static Long_t oldMres=0;
1622 static Long_t oldMvir=0;
1623 static Float_t oldCPU=0;
d97f1dbe 1624 static Long_t aveDMres=0;
1625 static Long_t aveDMvir=0;
1626 static Float_t aveDCPU=0;
470f88b0 1627
52dd4a8c 1628 AliCodeTimerAuto("",0);
21a3aa09 1629
0a035be5 1630 AliESDpid pid;
10d100d4 1631
21a3aa09 1632 if (iEvent >= fRunLoader->GetNumberOfEvents()) {
1633 fRunLoader->SetEventNumber(iEvent);
1634 fRunLoader->GetHeader()->Reset(fRawReader->GetRunNumber(),
1635 iEvent, iEvent);
21a3aa09 1636 fRunLoader->TreeE()->Fill();
4b71572b 1637 if (fRawReader && fRawReader->UseAutoSaveESD())
1638 fRunLoader->TreeE()->AutoSave("SaveSelf");
21a3aa09 1639 }
1640
1641 if ((iEvent < fFirstEvent) || ((fLastEvent >= 0) && (iEvent > fLastEvent))) {
21a3aa09 1642 return kTRUE;
1643 }
1644
aa3c69a9 1645
16447f18 1646 fRunLoader->GetEvent(iEvent);
1647
7e88424f 1648 // Fill Event-info object
1649 GetEventInfo();
77ba28ba 1650 fRecoParam.SetEventSpecie(fRunInfo,fEventInfo,fListOfCosmicTriggers);
d97f1dbe 1651
1652 ProcInfo_t procInfo;
1653 if(iEvent==fFirstEvent) {
1654 gSystem->GetProcInfo(&procInfo);
1655 oldMres=procInfo.fMemResident;
1656 oldMvir=procInfo.fMemVirtual;
1657 oldCPU=procInfo.fCpuUser+procInfo.fCpuSys;
1658 }
470f88b0 1659 AliInfo(Form("================================= Processing event %d of type %-10s ==================================", iEvent,fRecoParam.PrintEventSpecie()));
7e88424f 1660
a00021a7 1661 // Set the reco-params
1662 {
1663 TString detStr = fLoadCDB;
ac4a7581 1664 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
a00021a7 1665 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
1666 AliReconstructor *reconstructor = GetReconstructor(iDet);
1667 if (reconstructor && fRecoParam.GetDetRecoParamArray(iDet)) {
57acd2d2 1668 const AliDetectorRecoParam *par = fRecoParam.GetDetRecoParam(iDet);
1669 reconstructor->SetRecoParam(par);
0a035be5 1670 reconstructor->GetPidSettings(&pid);
10d100d4 1671 reconstructor->SetEventInfo(&fEventInfo);
57acd2d2 1672 if (fRunQA) {
eca4fa66 1673 AliQAManager::QAManager()->SetRecoParam(iDet, par) ;
619aafe6 1674 if (par) AliQAManager::QAManager()->SetEventSpecie(AliRecoParam::Convert(par->GetEventSpecie())) ;
57acd2d2 1675 }
a00021a7 1676 }
1677 }
9dad1020 1678 if (fRunQA) {
1679 const AliDetectorRecoParam *grppar = fRecoParam.GetDetRecoParam(kNDetectors);
1680 AliQAManager::QAManager()->SetRecoParam(AliQAv1::kGLOBAL, grppar) ;
1681 AliQAManager::QAManager()->SetEventSpecie(AliRecoParam::Convert(grppar->GetEventSpecie())) ;
1682 }
a00021a7 1683 }
1684
ca13fb87 1685 // QA on single raw
514cb8c7 1686 if (fRunQA && IsInTasks(AliQAv1::kRAWS)) {
eca4fa66 1687 AliQAManager::QAManager()->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
1688 AliQAManager::QAManager()->RunOneEvent(fRawReader) ;
57acd2d2 1689 }
d506c543 1690 // local single event reconstruction
b26c3770 1691 if (!fRunLocalReconstruction.IsNull()) {
d506c543 1692 TString detectors=fRunLocalReconstruction;
a441bf51 1693 // run HLT event reconstruction first
d506c543 1694 // ;-( IsSelected changes the string
1695 if (IsSelected("HLT", detectors) &&
1696 !RunLocalEventReconstruction("HLT")) {
4b71572b 1697 if (fStopOnError) {CleanUp(); return kFALSE;}
a441bf51 1698 }
1699 detectors=fRunLocalReconstruction;
1700 detectors.ReplaceAll("HLT", "");
1701 if (!RunLocalEventReconstruction(detectors)) {
13effe3f 1702 if (fStopOnError) {
1703 CleanUp();
1704 return kFALSE;
1705 }
b26c3770 1706 }
1707 }
1708
6b3a883f 1709
1710 // fill Event header information from the RawEventHeader
1711 if (fRawReader){FillRawEventHeaderESD(fesd);}
1712 if (fRawReader){FillRawEventHeaderESD(fhltesd);}
1713
21a3aa09 1714 fesd->SetRunNumber(fRunLoader->GetHeader()->GetRun());
1715 fhltesd->SetRunNumber(fRunLoader->GetHeader()->GetRun());
1716 fesd->SetEventNumberInFile(fRunLoader->GetHeader()->GetEventNrInRun());
1717 fhltesd->SetEventNumberInFile(fRunLoader->GetHeader()->GetEventNrInRun());
fd2e2210 1718
1719 fesd->SetEventSpecie(fRecoParam.GetEventSpecie());
1720 fhltesd->SetEventSpecie(fRecoParam.GetEventSpecie());
46698ae4 1721
d6ee376f 1722 // Set magnetic field from the tracker
21a3aa09 1723 fesd->SetMagneticField(AliTracker::GetBz());
1724 fhltesd->SetMagneticField(AliTracker::GetBz());
33fe5eb1 1725 //
c587bcde 1726 ((AliESDRun*)fesd->GetESDRun())->SetBeamEnergyIsSqrtSHalfGeV();
1727 ((AliESDRun*)fhltesd->GetESDRun())->SetBeamEnergyIsSqrtSHalfGeV();
1728 //
33fe5eb1 1729 AliMagF* fld = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
1730 if (fld) { // set info needed for field initialization
1731 fesd->SetCurrentL3(fld->GetCurrentSol());
1732 fesd->SetCurrentDip(fld->GetCurrentDip());
1733 fesd->SetBeamEnergy(fld->GetBeamEnergy());
1734 fesd->SetBeamType(fld->GetBeamTypeText());
1735 fesd->SetUniformBMap(fld->IsUniform());
1736 fesd->SetBInfoStored();
1737 //
1738 fhltesd->SetCurrentL3(fld->GetCurrentSol());
1739 fhltesd->SetCurrentDip(fld->GetCurrentDip());
1740 fhltesd->SetBeamEnergy(fld->GetBeamEnergy());
1741 fhltesd->SetBeamType(fld->GetBeamTypeText());
1742 fhltesd->SetUniformBMap(fld->IsUniform());
1743 fhltesd->SetBInfoStored();
1744 }
1745 //
71f6cda4 1746 // Set most probable pt, for B=0 tracking
9257a1bd 1747 // Get the global reco-params. They are atposition 16 inside the array of detectors in fRecoParam
ac4a7581 1748 const AliGRPRecoParam *grpRecoParam = dynamic_cast<const AliGRPRecoParam*>(fRecoParam.GetDetRecoParam(kNDetectors));
71f6cda4 1749 if (grpRecoParam) AliExternalTrackParam::SetMostProbablePt(grpRecoParam->GetMostProbablePt());
46698ae4 1750
2e3550da 1751 // Fill raw-data error log into the ESD
21a3aa09 1752 if (fRawReader) FillRawDataErrorLog(iEvent,fesd);
2e3550da 1753
2257f27e 1754 // vertex finder
1755 if (fRunVertexFinder) {
4b71572b 1756 if (!RunVertexFinder(fesd)) {
1757 if (fStopOnError) {CleanUp(); return kFALSE;}
2257f27e 1758 }
1759 }
1760
58e8dc31 1761 // For Plane Efficiency: run the SPD trackleter
1762 if (fRunPlaneEff && fSPDTrackleter) {
1763 if (!RunSPDTrackleting(fesd)) {
1764 if (fStopOnError) {CleanUp(); return kFALSE;}
1765 }
1766 }
1767
e66fbafb 1768 // Muon tracking
b8cd5251 1769 if (!fRunTracking.IsNull()) {
e66fbafb 1770 if (fRunMuonTracking) {
21a3aa09 1771 if (!RunMuonTracking(fesd)) {
4b71572b 1772 if (fStopOnError) {CleanUp(); return kFALSE;}
24f7a148 1773 }
596a855f 1774 }
1775 }
1776
e66fbafb 1777 // barrel tracking
1778 if (!fRunTracking.IsNull()) {
0a035be5 1779 if (!RunTracking(fesd,pid)) {
4b71572b 1780 if (fStopOnError) {CleanUp(); return kFALSE;}
e66fbafb 1781 }
1782 }
21c573b7 1783
596a855f 1784 // fill ESD
1785 if (!fFillESD.IsNull()) {
d506c543 1786 TString detectors=fFillESD;
f6806ad3 1787 // run HLT first and on hltesd
d506c543 1788 // ;-( IsSelected changes the string
1789 if (IsSelected("HLT", detectors) &&
1790 !FillESD(fhltesd, "HLT")) {
4b71572b 1791 if (fStopOnError) {CleanUp(); return kFALSE;}
f6806ad3 1792 }
1793 detectors=fFillESD;
d5105682 1794 // Temporary fix to avoid problems with HLT that overwrites the offline ESDs
1795 if (detectors.Contains("ALL")) {
1796 detectors="";
ac4a7581 1797 for (Int_t idet=0; idet<kNDetectors; ++idet){
d5105682 1798 detectors += fgkDetectorName[idet];
1799 detectors += " ";
1800 }
1801 }
f6806ad3 1802 detectors.ReplaceAll("HLT", "");
1803 if (!FillESD(fesd, detectors)) {
4b71572b 1804 if (fStopOnError) {CleanUp(); return kFALSE;}
596a855f 1805 }
1806 }
6b3a883f 1807
596a855f 1808 // combined PID
0a035be5 1809 pid.MakePID(fesd);
596a855f 1810
b647652d 1811 if (fFillTriggerESD) {
4b71572b 1812 if (!FillTriggerESD(fesd)) {
1813 if (fStopOnError) {CleanUp(); return kFALSE;}
b647652d 1814 }
1815 }
f7812afc 1816 // Always fill scalers
1817 if (!FillTriggerScalers(fesd)) {
1818 if (fStopOnError) {CleanUp(); return kFALSE;}
1819 }
1820
b647652d 1821
21a3aa09 1822 ffile->cd();
a6ee503a 1823
3c3709c4 1824 //
67be5c77 1825 // Propagate track to the beam pipe (if not already done by ITS)
3c3709c4 1826 //
21a3aa09 1827 const Int_t ntracks = fesd->GetNumberOfTracks();
3c3709c4 1828 const Double_t kRadius = 2.8; //something less than the beam pipe radius
1829
1830 TObjArray trkArray;
1831 UShort_t *selectedIdx=new UShort_t[ntracks];
1832
1833 for (Int_t itrack=0; itrack<ntracks; itrack++){
98a50ff3 1834 const Double_t kMaxStep = 1; //max step over the material
3c3709c4 1835 Bool_t ok;
1836
21a3aa09 1837 AliESDtrack *track = fesd->GetTrack(itrack);
3c3709c4 1838 if (!track) continue;
1839
1840 AliExternalTrackParam *tpcTrack =
1841 (AliExternalTrackParam *)track->GetTPCInnerParam();
bcabd6af 1842 ok = kFALSE;
1843 if (tpcTrack)
1844 ok = AliTracker::
52dd4a8c 1845 PropagateTrackToBxByBz(tpcTrack,kRadius,track->GetMass(),kMaxStep,kFALSE);
43c9dae1 1846
3c3709c4 1847 if (ok) {
1848 Int_t n=trkArray.GetEntriesFast();
1849 selectedIdx[n]=track->GetID();
1850 trkArray.AddLast(tpcTrack);
1851 }
1852
3d65e645 1853 //Tracks refitted by ITS should already be at the SPD vertex
1854 if (track->IsOn(AliESDtrack::kITSrefit)) continue;
1855
1856 AliTracker::
52dd4a8c 1857 PropagateTrackToBxByBz(track,kRadius,track->GetMass(),kMaxStep,kFALSE);
1858 Double_t x[3]; track->GetXYZ(x);
1859 Double_t b[3]; AliTracker::GetBxByBz(x,b);
1860 track->RelateToVertexBxByBz(fesd->GetPrimaryVertexSPD(), b, kVeryBig);
3c3709c4 1861
3c3709c4 1862 }
1863
1864 //
1865 // Improve the reconstructed primary vertex position using the tracks
1866 //
59224b2b 1867 Bool_t runVertexFinderTracks = fRunVertexFinderTracks;
1868 if(fesd->GetPrimaryVertexSPD()) {
1869 TString vtitle = fesd->GetPrimaryVertexSPD()->GetTitle();
1870 if(vtitle.Contains("cosmics")) {
1871 runVertexFinderTracks=kFALSE;
1872 }
c060d7fe 1873 }
a00021a7 1874
1875 if (runVertexFinderTracks) {
3c3709c4 1876 // TPC + ITS primary vertex
f09c879d 1877 ftVertexer->SetITSMode();
f2a195c1 1878 ftVertexer->SetConstraintOff();
a00021a7 1879 // get cuts for vertexer from AliGRPRecoParam
1880 if (grpRecoParam) {
1881 Int_t nCutsVertexer = grpRecoParam->GetVertexerTracksNCuts();
1882 Double_t *cutsVertexer = new Double_t[nCutsVertexer];
1883 grpRecoParam->GetVertexerTracksCutsITS(cutsVertexer);
1884 ftVertexer->SetCuts(cutsVertexer);
1885 delete [] cutsVertexer; cutsVertexer = NULL;
dd15203b 1886 if(fDiamondProfile && grpRecoParam->GetVertexerTracksConstraintITS()) {
1887 if(fDiamondProfile->GetXRes()<kRadius) ftVertexer->SetVtxStart(fDiamondProfile); // apply constraint only if sigmax is smaller than the beam pipe radius
1888 }
43c9dae1 1889 }
21a3aa09 1890 AliESDVertex *pvtx=ftVertexer->FindPrimaryVertex(fesd);
3c3709c4 1891 if (pvtx) {
1892 if (pvtx->GetStatus()) {
c264b61b 1893 fesd->SetPrimaryVertexTracks(pvtx);
3c3709c4 1894 for (Int_t i=0; i<ntracks; i++) {
21a3aa09 1895 AliESDtrack *t = fesd->GetTrack(i);
52dd4a8c 1896 Double_t x[3]; t->GetXYZ(x);
1897 Double_t b[3]; AliTracker::GetBxByBz(x,b);
1898 t->RelateToVertexBxByBz(pvtx, b, kVeryBig);
3c3709c4 1899 }
1900 }
9bcaa1d7 1901 delete pvtx; pvtx=NULL;
3c3709c4 1902 }
1903
1904 // TPC-only primary vertex
f09c879d 1905 ftVertexer->SetTPCMode();
f2a195c1 1906 ftVertexer->SetConstraintOff();
a00021a7 1907 // get cuts for vertexer from AliGRPRecoParam
1908 if (grpRecoParam) {
1909 Int_t nCutsVertexer = grpRecoParam->GetVertexerTracksNCuts();
1910 Double_t *cutsVertexer = new Double_t[nCutsVertexer];
1911 grpRecoParam->GetVertexerTracksCutsTPC(cutsVertexer);
1912 ftVertexer->SetCuts(cutsVertexer);
1913 delete [] cutsVertexer; cutsVertexer = NULL;
dd15203b 1914 if(fDiamondProfileTPC && grpRecoParam->GetVertexerTracksConstraintTPC()) {
1915 if(fDiamondProfileTPC->GetXRes()<kRadius) ftVertexer->SetVtxStart(fDiamondProfileTPC); // apply constraint only if sigmax is smaller than the beam pipe radius
1916 }
43c9dae1 1917 }
21a3aa09 1918 pvtx=ftVertexer->FindPrimaryVertex(&trkArray,selectedIdx);
3c3709c4 1919 if (pvtx) {
1920 if (pvtx->GetStatus()) {
21a3aa09 1921 fesd->SetPrimaryVertexTPC(pvtx);
3d65e645 1922 for (Int_t i=0; i<ntracks; i++) {
1923 AliESDtrack *t = fesd->GetTrack(i);
52dd4a8c 1924 Double_t x[3]; t->GetXYZ(x);
1925 Double_t b[3]; AliTracker::GetBxByBz(x,b);
1926 t->RelateToVertexTPCBxByBz(pvtx, b, kVeryBig);
3c3709c4 1927 }
1928 }
9bcaa1d7 1929 delete pvtx; pvtx=NULL;
3c3709c4 1930 }
1931
1932 }
1933 delete[] selectedIdx;
1934
21a3aa09 1935 if(fDiamondProfile) fesd->SetDiamond(fDiamondProfile);
17c86e90 1936
c5e3e5d1 1937
d1683eef 1938 if (fRunV0Finder) {
1939 // V0 finding
1940 AliV0vertexer vtxer;
21a3aa09 1941 vtxer.Tracks2V0vertices(fesd);
5e4ff34d 1942
d1683eef 1943 if (fRunCascadeFinder) {
1944 // Cascade finding
1945 AliCascadeVertexer cvtxer;
21a3aa09 1946 cvtxer.V0sTracks2CascadeVertices(fesd);
d1683eef 1947 }
5e4ff34d 1948 }
1949
596a855f 1950 // write ESD
21a3aa09 1951 if (fCleanESD) CleanESD(fesd);
854c6476 1952
514cb8c7 1953 if (fRunQA && IsInTasks(AliQAv1::kESDS)) {
eca4fa66 1954 AliQAManager::QAManager()->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
aeb8fc30 1955 AliQAManager::QAManager()->RunOneEvent(fesd, fhltesd) ;
57acd2d2 1956 }
bea94759 1957 if (fRunGlobalQA) {
eca4fa66 1958 AliQADataMaker *qadm = AliQAManager::QAManager()->GetQADataMaker(AliQAv1::kGLOBAL);
57acd2d2 1959 qadm->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
514cb8c7 1960 if (qadm && IsInTasks(AliQAv1::kESDS))
4e25ac79 1961 qadm->Exec(AliQAv1::kESDS, fesd);
b03591ab 1962 }
854c6476 1963
8b12d288 1964 // copy HLT decision from HLTesd to esd
1965 // the most relevant information is stored in a reduced container in the esd,
1966 // while the full information can be found in the HLTesd
1967 TObject* pHLTSrc=fhltesd->FindListObject(AliESDHLTDecision::Name());
1968 TObject* pHLTTgt=fesd->FindListObject(AliESDHLTDecision::Name());
1969 if (pHLTSrc && pHLTTgt) {
1970 pHLTSrc->Copy(*pHLTTgt);
1971 }
1972
ee7c441c 1973 if (fWriteESDfriend)
21a3aa09 1974 fesd->GetESDfriend(fesdf);
ee7c441c 1975
21a3aa09 1976 ftree->Fill();
ee7c441c 1977 if (fWriteESDfriend) ftreeF->Fill();
500d54ab 1978
2d91a353 1979 // Auto-save the ESD tree in case of prompt reco @P2
be103ac8 1980 if (fRawReader && fRawReader->UseAutoSaveESD()) {
2d91a353 1981 ftree->AutoSave("SaveSelf");
ee7c441c 1982 if (fWriteESDfriend) ftreeF->AutoSave("SaveSelf");
be103ac8 1983 TFile *friendfile = (TFile *)(gROOT->GetListOfFiles()->FindObject("AliESDfriends.root"));
1984 if (friendfile) friendfile->Save();
1985 }
2d91a353 1986
500d54ab 1987 // write HLT ESD
21a3aa09 1988 fhlttree->Fill();
1d99986f 1989
14dd053c 1990 // call AliEVE
1991 if (fRunAliEVE) RunAliEVE();
1992
21a3aa09 1993 fesd->Reset();
1994 fhltesd->Reset();
5728d3d5 1995 if (fWriteESDfriend) {
21a3aa09 1996 fesdf->~AliESDfriend();
1997 new (fesdf) AliESDfriend(); // Reset...
5728d3d5 1998 }
a5fa6165 1999
ac4a7581 2000 gSystem->GetProcInfo(&procInfo);
d97f1dbe 2001 Long_t dMres=(procInfo.fMemResident-oldMres)/1024;
2002 Long_t dMvir=(procInfo.fMemVirtual-oldMvir)/1024;
2003 Float_t dCPU=procInfo.fCpuUser+procInfo.fCpuSys-oldCPU;
2004 aveDMres+=(dMres-aveDMres)/(iEvent-fFirstEvent+1);
2005 aveDMvir+=(dMvir-aveDMvir)/(iEvent-fFirstEvent+1);
2006 aveDCPU+=(dCPU-aveDCPU)/(iEvent-fFirstEvent+1);
2007 AliInfo(Form("======================= End Event %d: Res %d(%3d <%3d>) Vir %d(%3d <%3d>) CPU %5.2f <%5.2f> ===================",
2008 iEvent, procInfo.fMemResident/1024, dMres, aveDMres, procInfo.fMemVirtual/1024, dMvir, aveDMvir, dCPU, aveDCPU));
470f88b0 2009 oldMres=procInfo.fMemResident;
2010 oldMvir=procInfo.fMemVirtual;
2011 oldCPU=procInfo.fCpuUser+procInfo.fCpuSys;
a5fa6165 2012
ca13fb87 2013 fEventInfo.Reset();
ac4a7581 2014 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
42ec5d3d 2015 if (fReconstructor[iDet]) {
a00021a7 2016 fReconstructor[iDet]->SetRecoParam(NULL);
42ec5d3d 2017 fReconstructor[iDet]->SetEventInfo(NULL);
2018 }
1f26f3e7 2019 if (fTracker[iDet]) fTracker[iDet]->SetEventInfo(NULL);
ca13fb87 2020 }
2021
53c8f690 2022 if (fRunQA || fRunGlobalQA)
eca4fa66 2023 AliQAManager::QAManager()->Increment() ;
53c8f690 2024
ca13fb87 2025 return kTRUE;
21a3aa09 2026}
2027
2028//_____________________________________________________________________________
4b71572b 2029void AliReconstruction::SlaveTerminate()
21a3aa09 2030{
4b71572b 2031 // Finalize the run on the slave side
21a3aa09 2032 // Called after the exit
2033 // from the event loop
52dd4a8c 2034 AliCodeTimerAuto("",0);
21a3aa09 2035
2036 if (fIsNewRunLoader) { // galice.root didn't exist
2037 fRunLoader->WriteHeader("OVERWRITE");
2038 fRunLoader->CdGAFile();
2039 fRunLoader->Write(0, TObject::kOverwrite);
2040 }
2041
f747912b 2042 const TMap *cdbMap = AliCDBManager::Instance()->GetStorageMap();
2043 const TList *cdbList = AliCDBManager::Instance()->GetRetrievedIds();
2044
2045 TMap *cdbMapCopy = new TMap(cdbMap->GetEntries());
2046 cdbMapCopy->SetOwner(1);
2047 cdbMapCopy->SetName("cdbMap");
2048 TIter iter(cdbMap->GetTable());
2049
2050 TPair* pair = 0;
2051 while((pair = dynamic_cast<TPair*> (iter.Next()))){
2052 TObjString* keyStr = dynamic_cast<TObjString*> (pair->Key());
2053 TObjString* valStr = dynamic_cast<TObjString*> (pair->Value());
2054 cdbMapCopy->Add(new TObjString(keyStr->GetName()), new TObjString(valStr->GetName()));
2055 }
2056
2057 TList *cdbListCopy = new TList();
2058 cdbListCopy->SetOwner(1);
2059 cdbListCopy->SetName("cdbList");
2060
2061 TIter iter2(cdbList);
2062
b940cb9b 2063 AliCDBId* id=0;
e84c88f5 2064 while((id = dynamic_cast<AliCDBId*> (iter2.Next()))){
a4970db9 2065 cdbListCopy->Add(new TObjString(id->ToString().Data()));
f747912b 2066 }
2067
21a3aa09 2068 ftree->GetUserInfo()->Add(cdbMapCopy);
2069 ftree->GetUserInfo()->Add(cdbListCopy);
abe0c04e 2070
46698ae4 2071
21a3aa09 2072 ffile->cd();
aa3c69a9 2073
562dd0b4 2074 // we want to have only one tree version number
21a3aa09 2075 ftree->Write(ftree->GetName(),TObject::kOverwrite);
63314086 2076 fhlttree->Write(fhlttree->GetName(),TObject::kOverwrite);
f3a97c86 2077
ee7c441c 2078 if (fWriteESDfriend) {
2079 ffileF->cd();
2080 ftreeF->Write(ftreeF->GetName(),TObject::kOverwrite);
2081 }
2082
a7a1e1c7 2083// Finish with Plane Efficiency evaluation: before of CleanUp !!!
2084 if (fRunPlaneEff && !FinishPlaneEff()) {
2085 AliWarning("Finish PlaneEff evaluation failed");
2086 }
2087
930e6e3e 2088 // End of cycle for the in-loop
87da0921 2089
2090 if (fRunQA || fRunGlobalQA) {
5cbed243 2091 AliQAManager::QAManager()->EndOfCycle() ;
81d57268 2092 if (fInput &&
2093 !fProofOutputLocation.IsNull() &&
2094 fProofOutputArchive.IsNull() &&
2095 !fProofOutputDataset) {
2096 TString qaOutputFile(Form("%sMerged.%s.Data.root",
2097 fProofOutputLocation.Data(),
2098 AliQAv1::GetQADataFileName()));
2099 TProofOutputFile *qaProofFile = new TProofOutputFile(Form("Merged.%s.Data.root",
2100 AliQAv1::GetQADataFileName()));
2101 qaProofFile->SetOutputFileName(qaOutputFile.Data());
2102 if (AliDebugLevel() > 0) qaProofFile->Dump();
2103 fOutput->Add(qaProofFile);
2104 MergeQA(qaProofFile->GetFileName());
87da0921 2105 }
2106 else {
2107 MergeQA();
2108 }
2109 }
2110
4b71572b 2111 gROOT->cd();
2112 CleanUp();
81d57268 2113
2114 if (fInput) {
2115 if (!fProofOutputFileName.IsNull() &&
2116 !fProofOutputLocation.IsNull() &&
2117 fProofOutputDataset &&
2118 !fProofOutputArchive.IsNull()) {
2119 TProofOutputFile *zipProofFile = new TProofOutputFile(fProofOutputFileName.Data(),
2120 "DROV",
2121 fProofOutputLocation.Data());
2122 if (AliDebugLevel() > 0) zipProofFile->Dump();
2123 fOutput->Add(zipProofFile);
2124 TString fileList(fProofOutputArchive.Data());
2125 fileList.ReplaceAll(","," ");
38c18bf1 2126 TString command;
2127#if ROOT_SVN_REVISION >= 30174
2128 command.Form("zip -n root %s/%s %s",zipProofFile->GetDir(kTRUE),zipProofFile->GetFileName(),fileList.Data());
2129#else
2130 command.Form("zip -n root %s/%s %s",zipProofFile->GetDir(),zipProofFile->GetFileName(),fileList.Data());
2131#endif
2132 AliInfo(Form("Executing: %s",command.Data()));
2133 gSystem->Exec(command.Data());
81d57268 2134 }
2135 }
4b71572b 2136}
2137
2138//_____________________________________________________________________________
2139void AliReconstruction::Terminate()
2140{
f3a97c86 2141 // Create tags for the events in the ESD tree (the ESD tree is always present)
2142 // In case of empty events the tags will contain dummy values
52dd4a8c 2143 AliCodeTimerAuto("",0);
4b71572b 2144
e6d66370 2145 // Do not call the ESD tag creator in case of PROOF-based reconstruction
2146 if (!fInput) {
2147 AliESDTagCreator *esdtagCreator = new AliESDTagCreator();
2148 esdtagCreator->CreateESDTags(fFirstEvent,fLastEvent,fGRPData, AliQAv1::Instance()->GetQA(), AliQAv1::Instance()->GetEventSpecies(), AliQAv1::kNDET, AliRecoParam::kNSpecies);
9bcaa1d7 2149 delete esdtagCreator;
e6d66370 2150 }
e84c88f5 2151
795e4a22 2152 // Cleanup of CDB manager: cache and active storages!
2153 AliCDBManager::Instance()->ClearCache();
596a855f 2154}
2155
2156//_____________________________________________________________________________
b26c3770 2157Bool_t AliReconstruction::RunLocalEventReconstruction(const TString& detectors)
2158{
2159// run the local reconstruction
aa3c69a9 2160
0f88822a 2161 static Int_t eventNr=0;
52dd4a8c 2162 AliCodeTimerAuto("",0)
b26c3770 2163
2164 TString detStr = detectors;
ac4a7581 2165 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b26c3770 2166 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
2167 AliReconstructor* reconstructor = GetReconstructor(iDet);
2168 if (!reconstructor) continue;
2169 AliLoader* loader = fLoader[iDet];
f6806ad3 2170 // Matthias April 2008: temporary fix to run HLT reconstruction
2171 // although the HLT loader is missing
2172 if (strcmp(fgkDetectorName[iDet], "HLT")==0) {
2173 if (fRawReader) {
44ed7a66 2174 reconstructor->Reconstruct(fRawReader, NULL);
f6806ad3 2175 } else {
44ed7a66 2176 TTree* dummy=NULL;
2177 reconstructor->Reconstruct(dummy, NULL);
f6806ad3 2178 }
2179 continue;
2180 }
d76c31f4 2181 if (!loader) {
2182 AliWarning(Form("No loader is defined for %s!",fgkDetectorName[iDet]));
2183 continue;
2184 }
b26c3770 2185 // conversion of digits
2186 if (fRawReader && reconstructor->HasDigitConversion()) {
2187 AliInfo(Form("converting raw data digits into root objects for %s",
2188 fgkDetectorName[iDet]));
30bbd491 2189// AliCodeTimerAuto(Form("converting raw data digits into root objects for %s",
52dd4a8c 2190// fgkDetectorName[iDet]),0);
b26c3770 2191 loader->LoadDigits("update");
2192 loader->CleanDigits();
2193 loader->MakeDigitsContainer();
2194 TTree* digitsTree = loader->TreeD();
2195 reconstructor->ConvertDigits(fRawReader, digitsTree);
2196 loader->WriteDigits("OVERWRITE");
2197 loader->UnloadDigits();
b26c3770 2198 }
b26c3770 2199 // local reconstruction
b26c3770 2200 AliInfo(Form("running reconstruction for %s", fgkDetectorName[iDet]));
52dd4a8c 2201 //AliCodeTimerAuto(Form("running reconstruction for %s", fgkDetectorName[iDet]),0);
b26c3770 2202 loader->LoadRecPoints("update");
2203 loader->CleanRecPoints();
2204 loader->MakeRecPointsContainer();
2205 TTree* clustersTree = loader->TreeR();
2206 if (fRawReader && !reconstructor->HasDigitConversion()) {
2207 reconstructor->Reconstruct(fRawReader, clustersTree);
2208 } else {
2209 loader->LoadDigits("read");
2210 TTree* digitsTree = loader->TreeD();
2211 if (!digitsTree) {
44ed7a66 2212 AliError(Form("Can't get the %s digits tree", fgkDetectorName[iDet]));
13effe3f 2213 if (fStopOnError)
2214 return kFALSE;
b26c3770 2215 } else {
44ed7a66 2216 reconstructor->Reconstruct(digitsTree, clustersTree);
514cb8c7 2217 if (fRunQA && IsInTasks(AliQAv1::kDIGITSR)) {
eca4fa66 2218 AliQAManager::QAManager()->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
2219 AliQAManager::QAManager()->RunOneEventInOneDetector(iDet, digitsTree) ;
44ed7a66 2220 }
b26c3770 2221 }
2222 loader->UnloadDigits();
2223 }
514cb8c7 2224 if (fRunQA && IsInTasks(AliQAv1::kRECPOINTS)) {
eca4fa66 2225 AliQAManager::QAManager()->SetEventSpecie(fRecoParam.GetEventSpecie()) ;
2226 AliQAManager::QAManager()->RunOneEventInOneDetector(iDet, clustersTree) ;
57acd2d2 2227 }
eca4fa66 2228 loader->WriteRecPoints("OVERWRITE");
2229 loader->UnloadRecPoints();
2230 AliSysInfo::AddStamp(Form("LRec%s_%d",fgkDetectorName[iDet],eventNr), iDet,1,eventNr);
a00021a7 2231 }
13effe3f 2232 IsSelected("CTP", detStr);
a00021a7 2233 if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
b26c3770 2234 AliError(Form("the following detectors were not found: %s",
2235 detStr.Data()));
13effe3f 2236 if (fStopOnError)
2237 return kFALSE;
b26c3770 2238 }
0f88822a 2239 eventNr++;
b26c3770 2240 return kTRUE;
2241}
58e8dc31 2242//_____________________________________________________________________________
2243Bool_t AliReconstruction::RunSPDTrackleting(AliESDEvent*& esd)
2244{
2245// run the SPD trackleting (for SPD efficiency purpouses)
2246
52dd4a8c 2247 AliCodeTimerAuto("",0)
58e8dc31 2248
2249 Double_t vtxPos[3] = {0, 0, 0};
2250 Double_t vtxErr[3] = {0.0, 0.0, 0.0};
2251/*
2252 TArrayF mcVertex(3);
2253 // if(MC)
2254 if (fRunLoader->GetHeader() && fRunLoader->GetHeader()->GenEventHeader()) {
2255 fRunLoader->GetHeader()->GenEventHeader()->PrimaryVertex(mcVertex);
2256 for (Int_t i = 0; i < 3; i++) vtxPos[i] = mcVertex[i];
2257 }
2258*/
2259 const AliESDVertex *vertex = esd->GetVertex();
2260 if(!vertex){
2261 AliWarning("Vertex not found");
2262 return kFALSE;
2263 }
2264 vertex->GetXYZ(vtxPos);
2265 vertex->GetSigmaXYZ(vtxErr);
2266 if (fSPDTrackleter) {
2267 AliInfo("running the SPD Trackleter for Plane Efficiency Evaluation");
2268
2269 // load clusters
2270 fLoader[0]->LoadRecPoints("read");
2271 TTree* tree = fLoader[0]->TreeR();
2272 if (!tree) {
2273 AliError("Can't get the ITS cluster tree");
2274 return kFALSE;
2275 }
2276 fSPDTrackleter->LoadClusters(tree);
2277 fSPDTrackleter->SetVertex(vtxPos, vtxErr);
2278 // run trackleting
2279 if (fSPDTrackleter->Clusters2Tracks(esd) != 0) {
84290fcc 2280 AliWarning("AliITSTrackleterSPDEff Clusters2Tracks failed");
58e8dc31 2281 // fLoader[0]->UnloadRecPoints();
2282 return kFALSE;
2283 }
2284//fSPDTrackleter->UnloadRecPoints();
2285 } else {
2286 AliWarning("SPDTrackleter not available");
2287 return kFALSE;
2288 }
2289 return kTRUE;
2290}
b26c3770 2291
2292//_____________________________________________________________________________
af885e0f 2293Bool_t AliReconstruction::RunVertexFinder(AliESDEvent*& esd)
596a855f 2294{
2295// run the barrel tracking
2296
52dd4a8c 2297 AliCodeTimerAuto("",0)
030b532d 2298
92bffc4d 2299 AliVertexer *vertexer = CreateVertexer();
2300 if (!vertexer) return kFALSE;
2301
2302 AliInfo("running the ITS vertex finder");
2257f27e 2303 AliESDVertex* vertex = NULL;
92bffc4d 2304 if (fLoader[0]) {
2305 fLoader[0]->LoadRecPoints();
2306 TTree* cltree = fLoader[0]->TreeR();
2307 if (cltree) {
2308 if(fDiamondProfileSPD) vertexer->SetVtxStart(fDiamondProfileSPD);
2309 vertex = vertexer->FindVertexForCurrentEvent(cltree);
308c2f7c 2310 }
2311 else {
92bffc4d 2312 AliError("Can't get the ITS cluster tree");
308c2f7c 2313 }
92bffc4d 2314 fLoader[0]->UnloadRecPoints();
2257f27e 2315 }
92bffc4d 2316 else {
2317 AliError("Can't get the ITS loader");
2318 }
2319 if(!vertex){
2320 AliWarning("Vertex not found");
2321 vertex = new AliESDVertex();
2322 vertex->SetName("default");
2323 }
2324 else {
2325 vertex->SetName("reconstructed");
2257f27e 2326 }
92bffc4d 2327
2328 Double_t vtxPos[3];
2329 Double_t vtxErr[3];
2330 vertex->GetXYZ(vtxPos);
2331 vertex->GetSigmaXYZ(vtxErr);
2332
06cc9d95 2333 esd->SetPrimaryVertexSPD(vertex);
73c51de2 2334 AliESDVertex *vpileup = NULL;
2335 Int_t novertices = 0;
2336 vpileup = vertexer->GetAllVertices(novertices);
2337 if(novertices>1){
2338 for (Int_t kk=1; kk<novertices; kk++)esd->AddPileupVertexSPD(&vpileup[kk]);
2339 }
32e449be 2340 // if SPD multiplicity has been determined, it is stored in the ESD
92bffc4d 2341 AliMultiplicity *mult = vertexer->GetMultiplicity();
32e449be 2342 if(mult)esd->SetMultiplicity(mult);
2343
ac4a7581 2344 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 2345 if (fTracker[iDet]) fTracker[iDet]->SetVertex(vtxPos, vtxErr);
2346 }
2257f27e 2347 delete vertex;
2348
92bffc4d 2349 delete vertexer;
2350
2257f27e 2351 return kTRUE;
2352}
2353
2354//_____________________________________________________________________________
af885e0f 2355Bool_t AliReconstruction::RunHLTTracking(AliESDEvent*& esd)
1f46a9ae 2356{
2357// run the HLT barrel tracking
2358
52dd4a8c 2359 AliCodeTimerAuto("",0)
1f46a9ae 2360
2361 if (!fRunLoader) {
2362 AliError("Missing runLoader!");
2363 return kFALSE;
2364 }
2365
2366 AliInfo("running HLT tracking");
2367
2368 // Get a pointer to the HLT reconstructor
ac4a7581 2369 AliReconstructor *reconstructor = GetReconstructor(kNDetectors-1);
1f46a9ae 2370 if (!reconstructor) return kFALSE;
2371
2372 // TPC + ITS
2373 for (Int_t iDet = 1; iDet >= 0; iDet--) {
2374 TString detName = fgkDetectorName[iDet];
2375 AliDebug(1, Form("%s HLT tracking", detName.Data()));
2376 reconstructor->SetOption(detName.Data());
d76c31f4 2377 AliTracker *tracker = reconstructor->CreateTracker();
1f46a9ae 2378 if (!tracker) {
2379 AliWarning(Form("couldn't create a HLT tracker for %s", detName.Data()));
2380 if (fStopOnError) return kFALSE;
9dcc06e1 2381 continue;
1f46a9ae 2382 }
2383 Double_t vtxPos[3];
2384 Double_t vtxErr[3]={0.005,0.005,0.010};
2385 const AliESDVertex *vertex = esd->GetVertex();
2386 vertex->GetXYZ(vtxPos);
2387 tracker->SetVertex(vtxPos,vtxErr);
2388 if(iDet != 1) {
2389 fLoader[iDet]->LoadRecPoints("read");
2390 TTree* tree = fLoader[iDet]->TreeR();
2391 if (!tree) {
2392 AliError(Form("Can't get the %s cluster tree", detName.Data()));
2393 return kFALSE;
2394 }
2395 tracker->LoadClusters(tree);
2396 }
2397 if (tracker->Clusters2Tracks(esd) != 0) {
2398 AliError(Form("HLT %s Clusters2Tracks failed", fgkDetectorName[iDet]));
2399 return kFALSE;
2400 }
2401 if(iDet != 1) {
2402 tracker->UnloadClusters();
2403 }
2404 delete tracker;
2405 }
2406
1f46a9ae 2407 return kTRUE;
2408}
2409
2410//_____________________________________________________________________________
af885e0f 2411Bool_t AliReconstruction::RunMuonTracking(AliESDEvent*& esd)
e66fbafb 2412{
2413// run the muon spectrometer tracking
2414
52dd4a8c 2415 AliCodeTimerAuto("",0)
e66fbafb 2416
2417 if (!fRunLoader) {
2418 AliError("Missing runLoader!");
2419 return kFALSE;
2420 }
2421 Int_t iDet = 7; // for MUON
2422
2423 AliInfo("is running...");
2424
2425 // Get a pointer to the MUON reconstructor
2426 AliReconstructor *reconstructor = GetReconstructor(iDet);
2427 if (!reconstructor) return kFALSE;
2428
2429
2430 TString detName = fgkDetectorName[iDet];
2431 AliDebug(1, Form("%s tracking", detName.Data()));
d76c31f4 2432 AliTracker *tracker = reconstructor->CreateTracker();
e66fbafb 2433 if (!tracker) {
2434 AliWarning(Form("couldn't create a tracker for %s", detName.Data()));
2435 return kFALSE;
2436 }
2437
e66fbafb 2438 // read RecPoints
761350a6 2439 fLoader[iDet]->LoadRecPoints("read");
c1954ee5 2440
761350a6 2441 tracker->LoadClusters(fLoader[iDet]->TreeR());
2442
2443 Int_t rv = tracker->Clusters2Tracks(esd);
2444
e66fbafb 2445 fLoader[iDet]->UnloadRecPoints();
2446
c1954ee5 2447 tracker->UnloadClusters();
2448
e66fbafb 2449 delete tracker;
2450
cb23c6ca 2451 if ( rv )
2452 {
2453 AliError(Form("%s Clusters2Tracks failed", fgkDetectorName[iDet]));
2454 return kFALSE;
2455 }
2456
e66fbafb 2457 return kTRUE;
2458}
2459
2460
2461//_____________________________________________________________________________
10d100d4 2462Bool_t AliReconstruction::RunTracking(AliESDEvent*& esd,AliESDpid &PID)
2257f27e 2463{
2464// run the barrel tracking
0f88822a 2465 static Int_t eventNr=0;
52dd4a8c 2466 AliCodeTimerAuto("",0)
24f7a148 2467
815c2b38 2468 AliInfo("running tracking");
596a855f 2469
1f26f3e7 2470 // Set the event info which is used
2471 // by the trackers in order to obtain
2472 // information about read-out detectors,
2473 // trigger etc.
2474 AliDebug(1, "Setting event info");
2475 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
2476 if (!fTracker[iDet]) continue;
2477 fTracker[iDet]->SetEventInfo(&fEventInfo);
2478 }
2479
91b876d1 2480 //Fill the ESD with the T0 info (will be used by the TOF)
d76c31f4 2481 if (fReconstructor[11] && fLoader[11]) {
2482 fLoader[11]->LoadRecPoints("READ");
2483 TTree *treeR = fLoader[11]->TreeR();
89916438 2484 if (treeR) {
2485 GetReconstructor(11)->FillESD((TTree *)NULL,treeR,esd);
2486 }
d76c31f4 2487 }
91b876d1 2488
b8cd5251 2489 // pass 1: TPC + ITS inwards
2490 for (Int_t iDet = 1; iDet >= 0; iDet--) {
2491 if (!fTracker[iDet]) continue;
2492 AliDebug(1, Form("%s tracking", fgkDetectorName[iDet]));
24f7a148 2493
b8cd5251 2494 // load clusters
2495 fLoader[iDet]->LoadRecPoints("read");
6efecea1 2496 AliSysInfo::AddStamp(Form("RLoadCluster%s_%d",fgkDetectorName[iDet],eventNr),iDet,1, eventNr);
b8cd5251 2497 TTree* tree = fLoader[iDet]->TreeR();
2498 if (!tree) {
2499 AliError(Form("Can't get the %s cluster tree", fgkDetectorName[iDet]));
24f7a148 2500 return kFALSE;
2501 }
b8cd5251 2502 fTracker[iDet]->LoadClusters(tree);
6efecea1 2503 AliSysInfo::AddStamp(Form("TLoadCluster%s_%d",fgkDetectorName[iDet],eventNr), iDet,2, eventNr);
b8cd5251 2504 // run tracking
2505 if (fTracker[iDet]->Clusters2Tracks(esd) != 0) {
2506 AliError(Form("%s Clusters2Tracks failed", fgkDetectorName[iDet]));
24f7a148 2507 return kFALSE;
2508 }
878e1fe1 2509 // preliminary PID in TPC needed by the ITS tracker
2510 if (iDet == 1) {
b26c3770 2511 GetReconstructor(1)->FillESD((TTree*)NULL, (TTree*)NULL, esd);
10d100d4 2512 PID.MakePID(esd,kTRUE);
0f88822a 2513 }
6efecea1 2514 AliSysInfo::AddStamp(Form("Tracking0%s_%d",fgkDetectorName[iDet],eventNr), iDet,3,eventNr);
b8cd5251 2515 }
596a855f 2516
b8cd5251 2517 // pass 2: ALL backwards
aa3c69a9 2518
ac4a7581 2519 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b8cd5251 2520 if (!fTracker[iDet]) continue;
2521 AliDebug(1, Form("%s back propagation", fgkDetectorName[iDet]));
2522
2523 // load clusters
2524 if (iDet > 1) { // all except ITS, TPC
2525 TTree* tree = NULL;
7b61cd9c 2526 fLoader[iDet]->LoadRecPoints("read");
6efecea1 2527 AliSysInfo::AddStamp(Form("RLoadCluster0%s_%d",fgkDetectorName[iDet],eventNr), iDet,1, eventNr);
7b61cd9c 2528 tree = fLoader[iDet]->TreeR();
b8cd5251 2529 if (!tree) {
eca4fa66 2530 AliError(Form("Can't get the %s cluster tree", fgkDetectorName[iDet]));
2531 return kFALSE;
24f7a148 2532 }
0f88822a 2533 fTracker[iDet]->LoadClusters(tree);
6efecea1 2534 AliSysInfo::AddStamp(Form("TLoadCluster0%s_%d",fgkDetectorName[iDet],eventNr), iDet,2, eventNr);
b8cd5251 2535 }
24f7a148 2536
b8cd5251 2537 // run tracking
283f39c6 2538 if (iDet>1) // start filling residuals for the "outer" detectors
eca4fa66 2539 if (fRunGlobalQA) {
2540 AliTracker::SetFillResiduals(fRecoParam.GetEventSpecie(), kTRUE);
2541 TObjArray ** arr = AliTracker::GetResidualsArray() ;
c9526f68 2542 if (arr) {
0a349581 2543 AliRecoParam::EventSpecie_t es=fRecoParam.GetEventSpecie();
2544 TObjArray * elem = arr[AliRecoParam::AConvert(es)];
c9526f68 2545 if ( elem && (! elem->At(0)) ) {
2546 AliQADataMaker *qadm = AliQAManager::QAManager()->GetQADataMaker(AliQAv1::kGLOBAL);
2547 if (qadm) qadm->InitRecPointsForTracker() ;
2548 }
2549 }
eca4fa66 2550 }
b8cd5251 2551 if (fTracker[iDet]->PropagateBack(esd) != 0) {
2552 AliError(Form("%s backward propagation failed", fgkDetectorName[iDet]));
49dfd67a 2553 // return kFALSE;
b8cd5251 2554 }
24f7a148 2555
b8cd5251 2556 // unload clusters
6e65648b 2557 if (iDet > 3) { // all except ITS, TPC, TRD and TOF
b8cd5251 2558 fTracker[iDet]->UnloadClusters();
7b61cd9c 2559 fLoader[iDet]->UnloadRecPoints();
b8cd5251 2560 }
8f37df88 2561 // updated PID in TPC needed by the ITS tracker -MI
2562 if (iDet == 1) {
10d100d4 2563 //GetReconstructor(1)->FillESD((TTree*)NULL, (TTree*)NULL, esd);
2564 //AliESDpid::MakePID(esd);
2565 PID.MakePID(esd,kTRUE);
8f37df88 2566 }
6efecea1 2567 AliSysInfo::AddStamp(Form("Tracking1%s_%d",fgkDetectorName[iDet],eventNr), iDet,3, eventNr);
b8cd5251 2568 }
283f39c6 2569 //stop filling residuals for the "outer" detectors
57acd2d2 2570 if (fRunGlobalQA) AliTracker::SetFillResiduals(fRecoParam.GetEventSpecie(), kFALSE);
596a855f 2571
b8cd5251 2572 // pass 3: TRD + TPC + ITS refit inwards
aa3c69a9 2573
b8cd5251 2574 for (Int_t iDet = 2; iDet >= 0; iDet--) {
2575 if (!fTracker[iDet]) continue;
2576 AliDebug(1, Form("%s inward refit", fgkDetectorName[iDet]));
596a855f 2577
b8cd5251 2578 // run tracking
283f39c6 2579 if (iDet<2) // start filling residuals for TPC and ITS
eca4fa66 2580 if (fRunGlobalQA) {
2581 AliTracker::SetFillResiduals(fRecoParam.GetEventSpecie(), kTRUE);
2582 TObjArray ** arr = AliTracker::GetResidualsArray() ;
c9526f68 2583 if (arr) {
0a349581 2584 AliRecoParam::EventSpecie_t es=fRecoParam.GetEventSpecie();
2585 TObjArray * elem = arr[AliRecoParam::AConvert(es)];
c9526f68 2586 if ( elem && (! elem->At(0)) ) {
2587 AliQADataMaker *qadm = AliQAManager::QAManager()->GetQADataMaker(AliQAv1::kGLOBAL);
2588 if (qadm) qadm->InitRecPointsForTracker() ;
2589 }
2590 }
eca4fa66 2591 }
2592
b8cd5251 2593 if (fTracker[iDet]->RefitInward(esd) != 0) {
2594 AliError(Form("%s inward refit failed", fgkDetectorName[iDet]));
49dfd67a 2595 // return kFALSE;
b8cd5251 2596 }
db2368d0 2597 // run postprocessing
2598 if (fTracker[iDet]->PostProcess(esd) != 0) {
2599 AliError(Form("%s postprocessing failed", fgkDetectorName[iDet]));
2600 // return kFALSE;
2601 }
6efecea1 2602 AliSysInfo::AddStamp(Form("Tracking2%s_%d",fgkDetectorName[iDet],eventNr), iDet,3, eventNr);
6e65648b 2603 }
2604
2605 // write space-points to the ESD in case alignment data output
2606 // is switched on
2607 if (fWriteAlignmentData)
2608 WriteAlignmentData(esd);
2609
2610 for (Int_t iDet = 3; iDet >= 0; iDet--) {
2611 if (!fTracker[iDet]) continue;
b8cd5251 2612 // unload clusters
2613 fTracker[iDet]->UnloadClusters();
6efecea1 2614 AliSysInfo::AddStamp(Form("TUnloadCluster%s_%d",fgkDetectorName[iDet],eventNr), iDet,4, eventNr);
b8cd5251 2615 fLoader[iDet]->UnloadRecPoints();
6efecea1 2616 AliSysInfo::AddStamp(Form("RUnloadCluster%s_%d",fgkDetectorName[iDet],eventNr), iDet,5, eventNr);
b8cd5251 2617 }
283f39c6 2618 // stop filling residuals for TPC and ITS
57acd2d2 2619 if (fRunGlobalQA) AliTracker::SetFillResiduals(fRecoParam.GetEventSpecie(), kFALSE);
854c6476 2620
0f88822a 2621 eventNr++;
596a855f 2622 return kTRUE;
2623}
2624
2625//_____________________________________________________________________________
d64bd07d 2626Bool_t AliReconstruction::CleanESD(AliESDEvent *esd){
2627 //
2628 // Remove the data which are not needed for the physics analysis.
2629 //
2630
d64bd07d 2631 Int_t nTracks=esd->GetNumberOfTracks();
18571674 2632 Int_t nV0s=esd->GetNumberOfV0s();
cf37fd88 2633 AliInfo
2634 (Form("Number of ESD tracks and V0s before cleaning: %d %d",nTracks,nV0s));
d64bd07d 2635
18571674 2636 Float_t cleanPars[]={fV0DCAmax,fV0CsPmin,fDmax,fZmax};
7f68891d 2637 Bool_t rc=esd->Clean(cleanPars);
d64bd07d 2638
7f68891d 2639 nTracks=esd->GetNumberOfTracks();
18571674 2640 nV0s=esd->GetNumberOfV0s();
cf37fd88 2641 AliInfo
ae5d5566 2642 (Form("Number of ESD tracks and V0s after cleaning %d %d",nTracks,nV0s));
d64bd07d 2643
7f68891d 2644 return rc;
d64bd07d 2645}
2646
2647//_____________________________________________________________________________
af885e0f 2648Bool_t AliReconstruction::FillESD(AliESDEvent*& esd, const TString& detectors)
596a855f 2649{
2650// fill the event summary data
2651
52dd4a8c 2652 AliCodeTimerAuto("",0)
0f88822a 2653 static Int_t eventNr=0;
596a855f 2654 TString detStr = detectors;
abe0c04e 2655
f1640d23 2656 AliSysInfo::AddStamp(Form("FillESDb%d",eventNr), -19,-19, eventNr);
ac4a7581 2657 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
abe0c04e 2658 if (!IsSelected(fgkDetectorName[iDet], detStr)) continue;
b8cd5251 2659 AliReconstructor* reconstructor = GetReconstructor(iDet);
2660 if (!reconstructor) continue;
4b71572b 2661 AliDebug(1, Form("filling ESD for %s", fgkDetectorName[iDet]));
2662 TTree* clustersTree = NULL;
2663 if (fLoader[iDet]) {
2664 fLoader[iDet]->LoadRecPoints("read");
2665 clustersTree = fLoader[iDet]->TreeR();
2666 if (!clustersTree) {
2667 AliError(Form("Can't get the %s clusters tree",
2668 fgkDetectorName[iDet]));
2669 if (fStopOnError) return kFALSE;
2670 }
2671 }
2672 if (fRawReader && !reconstructor->HasDigitConversion()) {
2673 reconstructor->FillESD(fRawReader, clustersTree, esd);
2674 } else {
2675 TTree* digitsTree = NULL;
d76c31f4 2676 if (fLoader[iDet]) {
4b71572b 2677 fLoader[iDet]->LoadDigits("read");
2678 digitsTree = fLoader[iDet]->TreeD();
2679 if (!digitsTree) {
2680 AliError(Form("Can't get the %s digits tree",
b26c3770 2681 fgkDetectorName[iDet]));
2682 if (fStopOnError) return kFALSE;
2683 }
2684 }
4b71572b 2685 reconstructor->FillESD(digitsTree, clustersTree, esd);
2686 if (fLoader[iDet]) fLoader[iDet]->UnloadDigits();
2687 }
2688 if (fLoader[iDet]) {
2689 fLoader[iDet]->UnloadRecPoints();
596a855f 2690 }
2691 }
13effe3f 2692
2693 IsSelected("CTP", detStr);
596a855f 2694 if ((detStr.CompareTo("ALL") != 0) && !detStr.IsNull()) {
815c2b38 2695 AliError(Form("the following detectors were not found: %s",
2696 detStr.Data()));
596a855f 2697 if (fStopOnError) return kFALSE;
2698 }
f1640d23 2699 AliSysInfo::AddStamp(Form("FillESDe%d",eventNr), -20,-20, eventNr);
0f88822a 2700 eventNr++;
596a855f 2701 return kTRUE;
2702}
2703
b647652d 2704//_____________________________________________________________________________
af885e0f 2705Bool_t AliReconstruction::FillTriggerESD(AliESDEvent*& esd)
b647652d 2706{
2707 // Reads the trigger decision which is
2708 // stored in Trigger.root file and fills
2709 // the corresponding esd entries
2710
52dd4a8c 2711 AliCodeTimerAuto("",0)
87932dab 2712
b647652d 2713 AliInfo("Filling trigger information into the ESD");
2714
2715 if (fRawReader) {
2716 AliCTPRawStream input(fRawReader);
2717 if (!input.Next()) {
7e88424f 2718 AliWarning("No valid CTP (trigger) DDL raw data is found ! The trigger info is taken from the event header!");
b647652d 2719 }
2720 else {
7e88424f 2721 if (esd->GetTriggerMask() != input.GetClassMask())
2722 AliError(Form("Invalid trigger pattern found in CTP raw-data: %llx %llx",
2723 input.GetClassMask(),esd->GetTriggerMask()));
2724 if (esd->GetOrbitNumber() != input.GetOrbitID())
2725 AliError(Form("Invalid orbit id found in CTP raw-data: %x %x",
2726 input.GetOrbitID(),esd->GetOrbitNumber()));
2727 if (esd->GetBunchCrossNumber() != input.GetBCID())
2728 AliError(Form("Invalid bunch-crossing id found in CTP raw-data: %x %x",
2729 input.GetBCID(),esd->GetBunchCrossNumber()));
e61ed4b1 2730 AliESDHeader* esdheader = esd->GetHeader();
2731 esdheader->SetL0TriggerInputs(input.GetL0Inputs());
2732 esdheader->SetL1TriggerInputs(input.GetL1Inputs());
2733 esdheader->SetL2TriggerInputs(input.GetL2Inputs());
a6dd87ad 2734 // IR
2735 UInt_t orbit=input.GetOrbitID();
2736 for(Int_t i=0 ; i<input.GetNIRs() ; i++ )
2737 if(TMath::Abs(Int_t(orbit-(input.GetIR(i))->GetOrbit()))<=1){
2738 esdheader->AddTriggerIR(input.GetIR(i));
2739 }
b647652d 2740 }
b024fd7f 2741 }
f7812afc 2742 return kTRUE;
2743}
2744//_____________________________________________________________________________
2745Bool_t AliReconstruction::FillTriggerScalers(AliESDEvent*& esd)
2746{
522fdd91 2747 //Scalers
82ebedd6 2748 //fRunScalers->Print();
2749 if(fRunScalers && fRunScalers->CheckRunScalers()){
a6dd87ad 2750 AliTimeStamp* timestamp = new AliTimeStamp(esd->GetOrbitNumber(), esd->GetPeriodNumber(), esd->GetBunchCrossNumber());
82ebedd6 2751 //AliTimeStamp* timestamp = new AliTimeStamp(10308000, 0, (ULong64_t)486238);
522fdd91 2752 AliESDHeader* esdheader = fesd->GetHeader();
2753 for(Int_t i=0;i<50;i++){
1e78ae8c 2754 if((1ull<<i) & esd->GetTriggerMask()){
6863d231 2755 AliTriggerScalersESD* scalesd = fRunScalers->GetScalersForEventClass( timestamp, i+1);
82ebedd6 2756 if(scalesd)esdheader->SetTriggerScalersRecord(scalesd);
522fdd91 2757 }
2758 }
2759 }
b647652d 2760 return kTRUE;
2761}
001397cd 2762//_____________________________________________________________________________
af885e0f 2763Bool_t AliReconstruction::FillRawEventHeaderESD(AliESDEvent*& esd)
001397cd 2764{
2765 //
2766 // Filling information from RawReader Header
2767 //
2768
151bea4e 2769 if (!fRawReader) return kFALSE;
2770
001397cd 2771 AliInfo("Filling information from RawReader Header");
151bea4e 2772
2773 esd->SetBunchCrossNumber(fRawReader->GetBCID());
2774 esd->SetOrbitNumber(fRawReader->GetOrbitID());
2775 esd->SetPeriodNumber(fRawReader->GetPeriod());
2776
2777 esd->SetTimeStamp(fRawReader->GetTimestamp());
2778 esd->SetEventType(fRawReader->GetType());
001397cd 2779
2780 return kTRUE;
2781}
2782
2783
596a855f 2784//_____________________________________________________________________________
2785Bool_t AliReconstruction::IsSelected(TString detName, TString& detectors) const
2786{
2787// check whether detName is contained in detectors
2788// if yes, it is removed from detectors
2789
2790 // check if all detectors are selected
2791 if ((detectors.CompareTo("ALL") == 0) ||
2792 detectors.BeginsWith("ALL ") ||
2793 detectors.EndsWith(" ALL") ||
2794 detectors.Contains(" ALL ")) {
2795 detectors = "ALL";
2796 return kTRUE;
2797 }
2798
2799 // search for the given detector
2800 Bool_t result = kFALSE;
2801 if ((detectors.CompareTo(detName) == 0) ||
2802 detectors.BeginsWith(detName+" ") ||
2803 detectors.EndsWith(" "+detName) ||
2804 detectors.Contains(" "+detName+" ")) {
2805 detectors.ReplaceAll(detName, "");
2806 result = kTRUE;
2807 }
2808
2809 // clean up the detectors string
2810 while (detectors.Contains(" ")) detectors.ReplaceAll(" ", " ");
2811 while (detectors.BeginsWith(" ")) detectors.Remove(0, 1);
2812 while (detectors.EndsWith(" ")) detectors.Remove(detectors.Length()-1, 1);
2813
2814 return result;
2815}
e583c30d 2816
2817//_____________________________________________________________________________
f08fc9f5 2818Bool_t AliReconstruction::InitRunLoader()
2819{
2820// get or create the run loader
2821
2822 if (gAlice) delete gAlice;
2823 gAlice = NULL;
2824
52dd4a8c 2825 TFile *gafile = TFile::Open(fGAliceFileName.Data());
2826 // if (!gSystem->AccessPathName(fGAliceFileName.Data())) { // galice.root exists
2827 if (gafile) { // galice.root exists
2828 gafile->Close();
2829 delete gafile;
2830
b26c3770 2831 // load all base libraries to get the loader classes
2832 TString libs = gSystem->GetLibraries();
ac4a7581 2833 for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
b26c3770 2834 TString detName = fgkDetectorName[iDet];
2835 if (detName == "HLT") continue;
2836 if (libs.Contains("lib" + detName + "base.so")) continue;
2837 gSystem->Load("lib" + detName + "base.so");
2838 }
f08fc9f5 2839 fRunLoader = AliRunLoader::Open(fGAliceFileName.Data());
2840 if (!fRunLoader) {
2841 AliError(Form("no run loader found in file %s", fGAliceFileName.Data()));
2842 CleanUp();
2843 return kFALSE;
2844 }
325aa76f 2845
b26c3770 2846 fRunLoader->CdGAFile();
325aa76f 2847 fRunLoader->LoadgAlice();
f08fc9f5 2848
6cae184e 2849 //PH This is a temporary fix to give access to the kinematics
2850 //PH that is needed for the labels of ITS clusters
f2ee4290 2851 fRunLoader->LoadHeader();
6cae184e 2852 fRunLoader->LoadKinematics();
2853
f08fc9f5 2854 } else { // galice.root does not exist
2855 if (!fRawReader) {
2856 AliError(Form("the file %s does not exist", fGAliceFileName.Data()));
f08fc9f5 2857 }
2858 fRunLoader = AliRunLoader::Open(fGAliceFileName.Data(),
2859 AliConfig::GetDefaultEventFolderName(),
2860 "recreate");
2861 if (!fRunLoader) {
2862 AliError(Form("could not create run loader in file %s",
2863 fGAliceFileName.Data()));
2864 CleanUp();
2865 return kFALSE;
2866 }
21a3aa09 2867 fIsNewRunLoader = kTRUE;
f08fc9f5 2868 fRunLoader->MakeTree("E");
21a3aa09 2869
973388c2 2870 if (fNumberOfEventsPerFile > 0)
2871 fRunLoader->SetNumberOfEventsPerFile(fNumberOfEventsPerFile);
2872 else
21a3aa09 2873 fRunLoader->SetNumberOfEventsPerFile((UInt_t)-1);
f08fc9f5 2874 }
2875
2876 return kTRUE;
2877}
2878
2879//_____________________________________________________________________________
b8cd5251 2880AliReconstructor* AliReconstruction::GetReconstructor(Int_t iDet)
c757bafd 2881{
f08fc9f5 2882// get the reconstructor object and the loader for a detector
c757bafd 2883
7e88424f 2884 if (fReconstructor[iDet]) {
2885 if (fRecoParam.GetDetRecoParamArray(iDet) && !AliReconstructor::GetRecoParam(iDet)) {
2886 const AliDetectorRecoParam *par = fRecoParam.GetDetRecoParam(iDet);
2887 fReconstructor[iDet]->SetRecoParam(par);
42ec5d3d 2888 fReconstructor[iDet]->SetRunInfo(fRunInfo);
7e88424f 2889 }
2890 return fReconstructor[iDet];
2891 }
b8cd5251 2892
2893 // load the reconstructor object
2894 TPluginManager* pluginManager = gROOT->GetPluginManager();
2895 TString detName = fgkDetectorName[iDet];
2896 TString recName = "Ali" + detName + "Reconstructor";
f0999a9a 2897
2898 if (!fIsNewRunLoader && !fRunLoader->GetLoader(detName+"Loader") && (detName != "HLT")) return NULL;
b8cd5251 2899
b8cd5251 2900 AliReconstructor* reconstructor = NULL;
2901 // first check if a plugin is defined for the reconstructor
2902 TPluginHandler* pluginHandler =
2903 pluginManager->FindHandler("AliReconstructor", detName);
f08fc9f5 2904 // if not, add a plugin for it
2905 if (!pluginHandler) {
b8cd5251 2906 AliDebug(1, Form("defining plugin for %s", recName.Data()));
b26c3770 2907 TString libs = gSystem->GetLibraries();
2908 if (libs.Contains("lib" + detName + "base.so") ||
2909 (gSystem->Load("lib" + detName + "base.so") >= 0)) {
b8cd5251 2910 pluginManager->AddHandler("AliReconstructor", detName,
2911 recName, detName + "rec", recName + "()");
2912 } else {
2913 pluginManager->AddHandler("AliReconstructor", detName,
2914 recName, detName, recName + "()");
c757bafd 2915 }
b8cd5251 2916 pluginHandler = pluginManager->FindHandler("AliReconstructor", detName);
2917 }
2918 if (pluginHandler && (pluginHandler->LoadPlugin() == 0)) {
2919 reconstructor = (AliReconstructor*) pluginHandler->ExecPlugin(0);
c757bafd 2920 }
b8cd5251 2921 if (reconstructor) {
2922 TObject* obj = fOptions.FindObject(detName.Data());
2923 if (obj) reconstructor->SetOption(obj->GetTitle());
1e500f25 2924 reconstructor->SetRunInfo(fRunInfo);
d76c31f4 2925 reconstructor->Init();
b8cd5251 2926 fReconstructor[iDet] = reconstructor;
2927 }
2928
f08fc9f5 2929 // get or create the loader
2930 if (detName != "HLT") {
2931 fLoader[iDet] = fRunLoader->GetLoader(detName + "Loader");
2932 if (!fLoader[iDet]) {
2933 AliConfig::Instance()
2934 ->CreateDetectorFolders(fRunLoader->GetEventFolder(),
2935 detName, detName);
2936 // first check if a plugin is defined for the loader
bb0901a4 2937 pluginHandler =
f08fc9f5 2938 pluginManager->FindHandler("AliLoader", detName);
2939 // if not, add a plugin for it
2940 if (!pluginHandler) {
2941 TString loaderName = "Ali" + detName + "Loader";
2942 AliDebug(1, Form("defining plugin for %s", loaderName.Data()));
2943 pluginManager->AddHandler("AliLoader", detName