1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 //_________________________________________________________________________
17 // Macros performing the full reconstruction chain starting from Digits
19 // root> .L Reconstruction.C++
20 // root> rec("RE", "PHOS EMCAL") --> does the reconstruction for
21 // PHOS and EMCAL and stores the
22 // reconstructed particles in
24 // author : Yves Schutz (CERN/SUBATECH)
26 //_________________________________________________________________________
28 #include "AliReconstruction.h"
30 #include "Riostream.h"
31 #include "AliPHOSGetter.h"
32 #include "AliEMCALGetter.h"
34 void reco(TString opt="TRE", TString name="all", Bool_t debug="kFALSE")
36 AliReconstruction rec ;
37 if ( !opt.Contains("T") ) {
38 rec.SetRunTracking(kFALSE) ;
39 rec.SetRunVertexFinder(kFALSE) ;
42 if ( opt.Contains("R") )
43 rec.SetRunLocalReconstruction(name.Data()) ;
45 rec.SetRunLocalReconstruction("") ;
47 if ( !opt.Contains("E") )
50 rec.SetFillESD(name.Data()) ;
54 if ( name.Contains("PHOS") ) {
55 cout << ">>>>>>>>>>>> PHOS " << endl ;
56 AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;
58 for (event = 0; event < gime->MaxEvent(); event++) {
59 cout << "event # " << event << endl ;
60 gime->Event(event, "RPT") ;
61 cout << " EMC RecPoints # " << gime->EmcRecPoints()->GetEntries() << endl ;
62 cout << " CPV RecPoints # " << gime->CpvRecPoints()->GetEntries() << endl ;
63 cout << " Track Segments # " << gime->TrackSegments()->GetEntries() << endl ;
64 cout << " Rec Particles # " << gime->RecParticles()->GetEntries() << endl ;
67 if ( name.Contains("EMCAL") ) {
68 cout << ">>>>>>>>>>>> EMCAL " << endl ;
69 AliEMCALGetter * gime = AliEMCALGetter::Instance("galice.root") ;
71 for (event = 0; event < gime->MaxEvent(); event++) {
72 cout << "event # " << event << endl ;
73 gime->Event(event, "RP") ;
74 cout << " RecPoints # " << gime->ECARecPoints()->GetEntries() << endl ;
75 cout << " Rec Particles # " << gime->RecParticles()->GetEntries() << endl ;