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")
36 AliReconstruction rec ;
37 if ( !opt.Contains("T") )
38 rec.SetRunTracking(kFALSE) ;
39 if ( opt.Contains("R") )
40 rec.SetRunReconstruction(name.Data()) ;
42 rec.SetRunReconstruction("") ;
43 if ( !opt.Contains("E") )
46 rec.SetFillESD(name.Data()) ;
49 if ( name.Contains("PHOS") ) {
50 cout << ">>>>>>>>>>>> PHOS " << endl ;
51 AliPHOSGetter * gime = AliPHOSGetter::Instance("galice.root") ;
53 for (event = 0; event < gime->MaxEvent(); event++) {
54 cout << "event # " << event << endl ;
55 gime->Event(event, "RPT") ;
56 cout << " EMC RecPoints # " << gime->EmcRecPoints()->GetEntries() << endl ;
57 cout << " CPV RecPoints # " << gime->CpvRecPoints()->GetEntries() << endl ;
58 cout << " Track Segments # " << gime->TrackSegments()->GetEntries() << endl ;
59 cout << " Rec Particles # " << gime->RecParticles()->GetEntries() << endl ;
62 if ( name.Contains("EMCAL") ) {
63 cout << ">>>>>>>>>>>> EMCAL " << endl ;
64 AliEMCALGetter * gime = AliEMCALGetter::Instance("galice.root") ;
66 for (event = 0; event < gime->MaxEvent(); event++) {
67 cout << "event # " << event << endl ;
68 gime->Event(event, "RP") ;
69 cout << " RecPoints # " << gime->ECARecPoints()->GetEntries() << endl ;
70 cout << " Rec Particles # " << gime->RecParticles()->GetEntries() << endl ;