#include "AliMUONDataInterface.h"
#include "AliMUONLocalTrigger.h"
+#include "AliMUONGlobalTrigger.h"
#include "AliMUONHit.h"
#include "AliMUONDigit.h"
#include "AliMUONRawCluster.h"
+#include "AliMUONTrack.h"
#include "AliLog.h"
+#include <iostream>
+using std::endl;
+using std::cout;
+
ClassImp(AliMUONDataInterface)
+///
+/// \class AliMUONDataInterface
+///
+/// An easy to use interface to the MUON module data stored in
+/// TreeK, TreeH, TreeS, TreeD and TreeR
+/// One can fetch any of the data objects with all the calls to runloader,
+/// muon loader and AliMUONData done behind the scenes and automatically.
+///
+/// This interface in not necessarily the fastest way to fetch the data but
+/// it is the easiest.
+/// Note: If independant calls to the run loader, muon loader or
+/// AliMUONData objects are interspersed with calls to the
+/// AliMUONDataInterface to fetch data, one might need to call the Reset
+/// method between these method calls at some point to prevent
+/// AliMUONDataInterface from getting confused.
+/// This is necessary since this object assumes the state of runloader,
+/// muon loader nor AliMUONData has not changed between calls.
+/// If the state has changes then one must call Reset so that
+/// AliMUONDataInterface refreshes what it knows about the state
+/// of the loader and AliMUONData objects.
+///
+/// \deprecated We have to revisit all this AliMUONData stuff anyway,
+/// and probably make a real AliMUONLoader instead...
+///
AliMUONDataInterface::AliMUONDataInterface()
: TObject(), fData(NULL, "MUON", "MUON")
// Set all internal pointers to NULL and indices to -1.
Reset();
-};
+}
AliMUONDataInterface::AliMUONDataInterface(const AliMUONDataInterface& rhs)
: TObject(rhs)
AliMUONDataInterface::~AliMUONDataInterface()
{
-// Delete the runloader when done.
+// Delete the runloader if we created it.
// If the runloader is not to be deleted then call Reset just before
// the destructor is called.
- if (fRunloader != NULL)
+ if (fRunloader != NULL && fCreatedRunLoader)
delete fRunloader;
-};
+}
AliMUONDataInterface&
AliMUONDataInterface::operator=(const AliMUONDataInterface& rhs)
// Note: No resources are released!
// Specificaly AliRunLoader is not deleted.
+ fCreatedRunLoader = kFALSE;
fRunloader = NULL;
fMuonloader = NULL;
fEventnumber = -1;
fDigitAddressSet = kFALSE;
fClusterAddressSet = kFALSE;
fTriggerAddressSet = kFALSE;
-};
+ fRecTracksAddressSet = kFALSE;
+}
-Bool_t AliMUONDataInterface::LoadLoaders(TString filename, TString foldername)
+Bool_t AliMUONDataInterface::UseCurrentRunLoader()
{
-// Load the run and muon loaders from the specified file and folder.
-// kTRUE is returned on success and kFALSE on failure.
+// Tries to fetch the current runloader with AliRunLoader::GetRunLoader. If nothing is
+// currently loaded then kFALSE is returned and AliMUONDataInterface is reset.
- fRunloader = AliRunLoader::Open(filename, foldername, "READ");
- if (fRunloader == NULL)
+ Reset();
+ fRunloader = AliRunLoader::GetRunLoader();
+ if (fRunloader == NULL) return kFALSE;
+ // Fetch the current file name, folder name and event number.
+ fFilename = fRunloader->GetFileName();
+ fFoldername = fRunloader->GetEventFolder()->GetName();
+ fEventnumber = fRunloader->GetEventNumber();
+
+ if ( ! FetchMuonLoader(fFilename.Data(), fFoldername.Data()) )
{
- AliError(Form("Could not find or load the run loader for the file: %s and folder: %s",
- (const char*)filename, (const char*)foldername));
+ Reset();
return kFALSE;
- };
+ }
+
+ return kTRUE;
+}
+
+
+Bool_t AliMUONDataInterface::FetchMuonLoader(TString filename, TString foldername)
+{
+ // fetches the muon loader for the given filename/foldername
+
fMuonloader = fRunloader->GetLoader("MUONLoader");
if (fMuonloader == NULL)
{
AliError(Form("Could not find the MUON loader in file: %s and folder: %s",
(const char*)filename, (const char*)foldername));
- fRunloader = NULL;
return kFALSE;
- };
+ }
// Need to connect the muon loader to the AliMUONData object,
// else class to fData will return NULL.
fData.SetLoader(fMuonloader);
+ return kTRUE;
+}
+
+
+Bool_t AliMUONDataInterface::LoadLoaders(TString filename, TString foldername)
+{
+// Load the run and muon loaders from the specified file and folder.
+// kTRUE is returned on success and kFALSE on failure.
+
+ fRunloader = AliRunLoader::Open(filename, foldername, "READ");
+ if (fRunloader == NULL)
+ {
+ AliError(Form("Could not find or load the run loader for the file: %s and folder: %s",
+ (const char*)filename, (const char*)foldername));
+ return kFALSE;
+ }
+ fCreatedRunLoader = kTRUE;
+ if ( ! FetchMuonLoader(filename, foldername) )
+ {
+ fRunloader = NULL;
+ return kFALSE;
+ }
fFilename = filename;
fFoldername = foldername;
fEventnumber = -1; // Reset the event number to force the event to be loaded.
return kTRUE;
-};
+}
Bool_t AliMUONDataInterface::FetchLoaders(TString filename, TString foldername)
fRunloader = AliRunLoader::GetRunLoader();
if (fRunloader == NULL)
return LoadLoaders(filename, foldername);
+ else
+ {
+ if (fMuonloader == NULL)
+ {
+ if ( ! FetchMuonLoader(filename, foldername) )
+ {
+ fRunloader = NULL;
+ return kFALSE;
+ }
+ }
+ }
// Fetch the current file and folder names.
fFilename = fRunloader->GetFileName();
fFoldername = fRunloader->GetEventFolder()->GetName();
- };
+ }
// If filename or foldername are not the same as the ones currently selected then
// reopen the file.
{
delete fRunloader;
return LoadLoaders(filename, foldername);
- };
+ }
return kTRUE;
-};
+}
Bool_t AliMUONDataInterface::FetchEvent(Int_t event)
{
// Fetch the specified event from the runloader and reset all the track, cathode
// and address flags to force them to be reloaded.
+// If a negative event number is specified then the current runloader event
+// number is used.
if (fEventnumber < 0)
{
fDigitAddressSet = kFALSE;
fClusterAddressSet = kFALSE;
fTriggerAddressSet = kFALSE;
- };
+ fRecTracksAddressSet = kFALSE;
+ }
if ( event != fEventnumber )
{
if ( fRunloader->GetEvent(event) < 0 ) return kFALSE;
fDigitAddressSet = kFALSE;
fClusterAddressSet = kFALSE;
fTriggerAddressSet = kFALSE;
- };
+ fRecTracksAddressSet = kFALSE;
+ }
return kTRUE;
-};
+}
Bool_t AliMUONDataInterface::FetchTreeK()
{
AliError("Could not load TreeK.");
return kFALSE;
- };
- };
+ }
+ }
return kTRUE;
-};
+}
Bool_t AliMUONDataInterface::FetchTreeH()
{
AliError("Could not load TreeH.");
return kFALSE;
- };
+ }
fData.SetTreeAddress("H");
fHitAddressSet = kTRUE;
}
{
fData.SetTreeAddress("H");
fHitAddressSet = kTRUE;
- };
+ }
return kTRUE;
-};
+}
Bool_t AliMUONDataInterface::FetchTreeS()
{
AliError("Could not load TreeS.");
return kFALSE;
- };
+ }
fData.SetTreeAddress("S");
fSDigitAddressSet = kTRUE;
}
{
fData.SetTreeAddress("S");
fSDigitAddressSet = kTRUE;
- };
+ }
return kTRUE;
-};
+}
Bool_t AliMUONDataInterface::FetchTreeD()
{
AliError("Could not load TreeD.");
return kFALSE;
- };
+ }
fData.SetTreeAddress("D");
fDigitAddressSet = kTRUE;
}
{
fData.SetTreeAddress("D");
fDigitAddressSet = kTRUE;
- };
+ }
return kTRUE;
-};
+}
Bool_t AliMUONDataInterface::FetchTreeR()
{
-// Fetch the reconstructed objects tree from the current muon loader.
-// Nore: The addresses must still be set.
-
- if (fMuonloader->TreeR() == NULL)
+ // Fetch the reconstructed objects tree from the current muon loader.
+ // Note: The addresses must still be set.
+
+ if (fMuonloader->TreeR() == NULL)
+ {
+ fMuonloader->LoadRecPoints("READ");
+ if (fMuonloader->TreeR() == NULL)
{
- fMuonloader->LoadRecPoints("READ");
- if (fMuonloader->TreeR() == NULL)
- {
- AliError("Could not load TreeR.");
- return kFALSE;
- };
-
- // Need to reset these flags so that the cluster and trigger address
- // gets reset after this method.
- fClusterAddressSet = kFALSE;
- fTriggerAddressSet = kFALSE;
- };
- return kTRUE;
-};
-
+ AliError("Could not load TreeR.");
+ return kFALSE;
+ }
+
+ // Need to reset these flags so that the cluster and trigger address
+ // gets reset after this method.
+ fClusterAddressSet = kFALSE;
+ fTriggerAddressSet = kFALSE;
+ }
+ return kTRUE;
+}
+Bool_t AliMUONDataInterface::FetchTreeT()
+{
+ // fetch the reconstructed tracks tree from the current muon loader
+ // note : the addresses must still be set.
+ if (fMuonloader->TreeT() == NULL)
+ {
+ fMuonloader->LoadTracks("READ");
+ if (fMuonloader->TreeT() == NULL)
+ {
+ AliError("Could not load TreeT.");
+ return kFALSE;
+ }
+
+ // Need to reset these flags so that the rec tracks address
+ // gets reset after this method.
+ fRecTracksAddressSet = kFALSE;
+ }
+ return kTRUE;
+}
+
Int_t AliMUONDataInterface::NumberOfEvents(TString filename, TString foldername)
{
// Returns the number of events in the specified file/folder, and -1 on error.
if ( ! FetchLoaders(filename, foldername) ) return -1;
return fRunloader->GetNumberOfEvents();
-};
+}
Int_t AliMUONDataInterface::NumberOfParticles(TString filename, TString foldername, Int_t event)
if ( ! FetchEvent(event) ) return -1;
if ( ! FetchTreeK() ) return -1;
return (Int_t) fRunloader->TreeK()->GetEntriesFast();
-};
+}
TParticle* AliMUONDataInterface::Particle(
treeK->GetBranch("Particles")->SetAddress(&p);
treeK->GetEvent(particle);
return p;
-};
+}
Int_t AliMUONDataInterface::NumberOfTracks(TString filename, TString foldername, Int_t event)
if ( ! FetchEvent(event) ) return -1;
if ( ! FetchTreeH() ) return -1;
return fData.GetNtracks();
-};
+}
Int_t AliMUONDataInterface::NumberOfHits(
fData.ResetHits();
fData.GetTrack(track);
fTrack = track;
- };
+ }
return fData.Hits()->GetEntriesFast();
-};
+}
AliMUONHit* AliMUONDataInterface::Hit(
fData.ResetHits();
fData.GetTrack(track);
fTrack = track;
- };
+ }
return static_cast<AliMUONHit*>( fData.Hits()->At(hit) );
-};
+}
Int_t AliMUONDataInterface::NumberOfSDigits(
if ( fSCathode != cathode )
{
fData.ResetSDigits();
- fData.GetCathodeS(cathode);
+ fData.GetSDigits();
fSCathode = cathode;
- };
+ }
return fData.SDigits(chamber)->GetEntriesFast();
-};
+}
AliMUONDigit* AliMUONDataInterface::SDigit(
if ( fSCathode != cathode )
{
fData.ResetSDigits();
- fData.GetCathodeS(cathode);
+ fData.GetSDigits();
fSCathode = cathode;
- };
+ }
return static_cast<AliMUONDigit*>( fData.SDigits(chamber)->At(sdigit) );
-};
+}
Int_t AliMUONDataInterface::NumberOfDigits(
if ( fCathode != cathode )
{
fData.ResetDigits();
- fData.GetCathode(cathode);
+ fData.GetDigits();
fCathode = cathode;
- };
+ }
return fData.Digits(chamber)->GetEntriesFast();
-};
+}
AliMUONDigit* AliMUONDataInterface::Digit(
if ( fCathode != cathode )
{
fData.ResetDigits();
- fData.GetCathode(cathode);
+ fData.GetDigits();
fCathode = cathode;
- };
+ }
return static_cast<AliMUONDigit*>( fData.Digits(chamber)->At(digit) );
-};
+}
Int_t AliMUONDataInterface::NumberOfRawClusters(
fData.ResetRawClusters();
fData.GetRawClusters();
fClusterAddressSet = kTRUE;
- };
+ }
return fData.RawClusters(chamber)->GetEntriesFast();
-};
+}
AliMUONRawCluster* AliMUONDataInterface::RawCluster(
fData.ResetRawClusters();
fData.GetRawClusters();
fClusterAddressSet = kTRUE;
- };
+ }
return static_cast<AliMUONRawCluster*>( fData.RawClusters(chamber)->At(cluster) );
-};
+}
Int_t AliMUONDataInterface::NumberOfLocalTriggers(TString filename, TString foldername, Int_t event)
if ( ! FetchLoaders(filename, foldername) ) return -1;
if ( ! FetchEvent(event) ) return -1;
- if ( ! FetchTreeR() ) return -1;
+ if ( ! FetchTreeD() ) return -1;
if ( ! fTriggerAddressSet )
{
// If the local trigger address in TreeR is not set yet then set it now.
fData.SetTreeAddress("GLT");
fData.ResetTrigger();
- fData.GetTrigger();
+ fData.GetTriggerD();
fTriggerAddressSet = kTRUE;
- };
+ }
return fData.LocalTrigger()->GetEntriesFast();
-};
+}
AliMUONLocalTrigger* AliMUONDataInterface::LocalTrigger(
if ( ! FetchLoaders(filename, foldername) ) return NULL;
if ( ! FetchEvent(event) ) return NULL;
- if ( ! FetchTreeR() ) return NULL;
+ if ( ! FetchTreeD() ) return NULL;
if ( ! fTriggerAddressSet )
{
// If the local trigger address in TreeR is not set yet then set it now.
fData.SetTreeAddress("GLT");
fData.ResetTrigger();
- fData.GetTrigger();
+ fData.GetTriggerD();
fTriggerAddressSet = kTRUE;
- };
+ }
return static_cast<AliMUONLocalTrigger*>( fData.LocalTrigger()->At(trigger) );
-};
-
+}
Bool_t AliMUONDataInterface::SetFile(TString filename, TString foldername)
{
// kTRUE is returned if the run and muon loaders were found, else kFALSE.
return FetchLoaders(filename, foldername);
-};
+}
Bool_t AliMUONDataInterface::GetEvent(Int_t event)
// Select the current event from which to fetch data.
// kTRUE is returned if the event was found, else kFALSE is returned.
- return FetchEvent(event);
-};
+ if (fRunloader == NULL)
+ {
+ AliError("File not set.");
+ return kFALSE;
+ }
+ else
+ return FetchEvent(event);
+}
Int_t AliMUONDataInterface::NumberOfEvents()
{
AliError("File not set.");
return -1;
- };
+ }
return fRunloader->GetNumberOfEvents();
-};
+}
Int_t AliMUONDataInterface::NumberOfParticles()
{
AliError("File not set.");
return -1;
- };
+ }
if ( ! FetchTreeK() ) return -1;
return (Int_t) fRunloader->TreeK()->GetEntriesFast();
-};
+}
TParticle* AliMUONDataInterface::Particle(Int_t particle)
{
AliError("File not set.");
return NULL;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return NULL;
- };
+ }
if ( ! FetchTreeK() ) return NULL;
TTree* treeK = fRunloader->TreeK();
TParticle* p = NULL;
treeK->GetBranch("Particles")->SetAddress(&p);
treeK->GetEvent(particle);
return p;
-};
+}
Int_t AliMUONDataInterface::NumberOfTracks()
{
AliError("File not set.");
return -1;
- };
+ }
if (fEventnumber < 0)
{
AliError( "Event not chosen.");
return -1;
- };
+ }
if ( ! FetchTreeH() ) return -1;
return fData.GetNtracks();
-};
+}
Int_t AliMUONDataInterface::NumberOfHits(Int_t track)
{
AliError("File not set.");
return -1;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return -1;
- };
+ }
if ( ! FetchTreeH() ) return -1;
if (fTrack < 0 || fTrack != track)
{
fData.ResetHits();
fData.GetTrack(track);
fTrack = track;
- };
+ }
return fData.Hits()->GetEntriesFast();
-};
+}
AliMUONHit* AliMUONDataInterface::Hit(Int_t track, Int_t hit)
{
AliError("File not set.");
return NULL;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return NULL;
- };
+ }
if ( ! FetchTreeH() ) return NULL;
if (fTrack < 0 || fTrack != track)
{
fData.ResetHits();
fData.GetTrack(track);
fTrack = track;
- };
+ }
return static_cast<AliMUONHit*>( fData.Hits()->At(hit) );
-};
+}
Int_t AliMUONDataInterface::NumberOfSDigits(Int_t chamber, Int_t cathode)
{
AliError("File not set.");
return -1;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return -1;
- };
+ }
if ( ! FetchTreeS() ) return -1;
if ( fSCathode != cathode )
{
fData.ResetSDigits();
- fData.GetCathodeS(cathode);
+ fData.GetSDigits();
fSCathode = cathode;
- };
+ }
return fData.SDigits(chamber)->GetEntriesFast();
-};
+}
AliMUONDigit* AliMUONDataInterface::SDigit(Int_t chamber, Int_t cathode, Int_t sdigit)
{
AliError("File not set.");
return NULL;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return NULL;
- };
+ }
if ( ! FetchTreeS() ) return NULL;
if ( fSCathode != cathode )
{
fData.ResetSDigits();
- fData.GetCathodeS(cathode);
+ fData.GetSDigits();
fSCathode = cathode;
- };
+ }
return static_cast<AliMUONDigit*>( fData.SDigits(chamber)->At(sdigit) );
-};
+}
Int_t AliMUONDataInterface::NumberOfDigits(Int_t chamber, Int_t cathode)
{
AliError("File not set.");
return -1;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return -1;
- };
+ }
if ( ! FetchTreeD() ) return -1;
if ( fCathode != cathode )
{
fData.ResetDigits();
- fData.GetCathode(cathode);
+ fData.GetDigits();
fCathode = cathode;
- };
+ }
return fData.Digits(chamber)->GetEntriesFast();
-};
+}
AliMUONDigit* AliMUONDataInterface::Digit(Int_t chamber, Int_t cathode, Int_t digit)
{
AliError("File not set.");
return NULL;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return NULL;
- };
+ }
if ( ! FetchTreeD() ) return NULL;
if ( fCathode != cathode )
{
fData.ResetDigits();
- fData.GetCathode(cathode);
+ fData.GetDigits();
fCathode = cathode;
- };
+ }
return static_cast<AliMUONDigit*>( fData.Digits(chamber)->At(digit) );
-};
+}
Int_t AliMUONDataInterface::NumberOfRawClusters(Int_t chamber)
{
AliError("File not set.");
return -1;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return -1;
- };
+ }
if ( ! FetchTreeR() ) return -1;
if ( ! fClusterAddressSet )
fData.ResetRawClusters();
fData.GetRawClusters();
fClusterAddressSet = kTRUE;
- };
+ }
return fData.RawClusters(chamber)->GetEntriesFast();
-};
+}
AliMUONRawCluster* AliMUONDataInterface::RawCluster(Int_t chamber, Int_t cluster)
{
AliError("File not set.");
return NULL;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return NULL;
- };
+ }
if ( ! FetchTreeR() ) return NULL;
if ( ! fClusterAddressSet )
fData.ResetRawClusters();
fData.GetRawClusters();
fClusterAddressSet = kTRUE;
- };
+ }
return static_cast<AliMUONRawCluster*>( fData.RawClusters(chamber)->At(cluster) );
-};
+}
Int_t AliMUONDataInterface::NumberOfLocalTriggers()
{
AliError("File not set.");
return -1;
- };
+ }
if (fEventnumber < 0)
{
AliError("Event not chosen.");
return -1;
- };
+ }
- if ( ! FetchTreeR() ) return -1;
+ if ( ! FetchTreeD() ) return -1;
if ( ! fTriggerAddressSet )
{
fData.SetTreeAddress("GLT");
fData.ResetTrigger();
- fData.GetTrigger();
+ fData.GetTriggerD();
fTriggerAddressSet = kTRUE;
- };
+ }
return fData.LocalTrigger()->GetEntriesFast();
-};
+}
AliMUONLocalTrigger* AliMUONDataInterface::LocalTrigger(Int_t trigger)
{
AliError("File not set.");
return NULL;
- };
+ }
if (fEventnumber < 0)
{
AliError( "Event not chosen.");
return NULL;
- };
+ }
- if ( ! FetchTreeR() ) return NULL;
+ if ( ! FetchTreeD() ) return NULL;
if ( ! fTriggerAddressSet )
{
fData.SetTreeAddress("GLT");
fData.ResetTrigger();
- fData.GetTrigger();
+ fData.GetTriggerD();
fTriggerAddressSet = kTRUE;
- };
+ }
return static_cast<AliMUONLocalTrigger*>( fData.LocalTrigger()->At(trigger) );
-};
+}
+
+Int_t AliMUONDataInterface::NumberOfGlobalTriggers()
+{
+
+ // Get the number of local trigger objects in the current event.
+ // -1 is returned on error.
+
+ if (fRunloader == NULL)
+ {
+ AliError("File not set.");
+ return -1;
+ }
+ if (fEventnumber < 0)
+ {
+ AliError("Event not chosen.");
+ return -1;
+ }
+
+ if ( ! FetchTreeD() ) return -1;
+ if ( ! fTriggerAddressSet )
+ {
+ fData.SetTreeAddress("GLT");
+ fData.ResetTrigger();
+ fData.GetTriggerD();
+ fTriggerAddressSet = kTRUE;
+ }
+ return fData.GlobalTrigger()->GetEntriesFast();
+}
+
+AliMUONGlobalTrigger* AliMUONDataInterface::GlobalTrigger(Int_t trigger)
+{
+ // Fetch the specified local trigger object from the current event.
+ // NULL is returned on error.
+
+ if (fRunloader == NULL)
+ {
+ AliError("File not set.");
+ return NULL;
+ }
+ if (fEventnumber < 0)
+ {
+ AliError( "Event not chosen.");
+ return NULL;
+ }
+
+ if ( ! FetchTreeD() ) return NULL;
+ if ( ! fTriggerAddressSet )
+ {
+ fData.SetTreeAddress("GLT");
+ fData.ResetTrigger();
+ fData.GetTriggerD();
+ fTriggerAddressSet = kTRUE;
+ }
+ return static_cast<AliMUONGlobalTrigger*>( fData.GlobalTrigger()->At(trigger) );
+}
+
+Int_t AliMUONDataInterface::NumberOfRecTracks()
+{
+ // Fetch the number of reconstructed tracks from the current event.
+ // NULL is returned on error.
+
+ if (fRunloader == NULL)
+ {
+ AliError("File not set.");
+ return -1;
+ }
+ if (fEventnumber < 0)
+ {
+ AliError( "Event not chosen.");
+ return -1;
+ }
+
+ if ( ! FetchTreeT() ) return -1;
+ if ( ! fRecTracksAddressSet )
+ {
+ fData.SetTreeAddress("RT");
+ fData.ResetRecTracks();
+ fData.GetRecTracks();
+ fRecTracksAddressSet = kTRUE;
+ }
+ return fData.RecTracks()->GetEntriesFast();
+}
+
+AliMUONTrack* AliMUONDataInterface::RecTrack(Int_t rectrack)
+{
+ // Fetch the specified reconstructed track object from the current event.
+ // NULL is returned on error.
+
+ if (fRunloader == NULL)
+ {
+ AliError("File not set.");
+ return NULL;
+ }
+ if (fEventnumber < 0)
+ {
+ AliError( "Event not chosen.");
+ return NULL;
+ }
+
+ if ( ! FetchTreeT() ) return NULL;
+ if ( ! fRecTracksAddressSet )
+ {
+ fData.SetTreeAddress("RT");
+ fData.ResetRecTracks();
+ fData.GetRecTracks();
+ fRecTracksAddressSet = kTRUE;
+ }
+ return static_cast<AliMUONTrack*>( fData.RecTracks()->At(rectrack) );
+ // return (AliMUONTrack*)(fData.RecTracks()->At(rectrack));
+}