]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - MUON/AliMUONDataInterface.cxx
Optional geometry without CPV
[u/mrichter/AliRoot.git] / MUON / AliMUONDataInterface.cxx
index 2bfa83bbcedd4ea8e1c39a9f6ba257230c8e8a5c..b132191e1ad5b959faebb2898b4bc0b1ef2be3fa 100644 (file)
@@ -1,7 +1,21 @@
-
-// Author: Artur Szostak
-//  email: artur@alice.phy.uct.ac.za
-
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+/* $Id$ */
+
+#include <cassert>
 #include <TError.h>
 #include <TParticle.h>
 
 
 #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;
+
+///
+/// \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...
+///
+/// \author Artur Szostak
+/// email: artur@alice.phy.uct.ac.za
+
+
+/// \cond CLASSIMP
 ClassImp(AliMUONDataInterface)
-
+/// \endcond
 
 AliMUONDataInterface::AliMUONDataInterface()
-       : TObject(), fData(NULL, "MUON", "MUON")
+       : TObject(), 
+         fCreatedRunLoader(kFALSE),
+         fHitAddressSet(kFALSE),
+         fSDigitAddressSet(kFALSE),
+         fDigitAddressSet(kFALSE),
+         fClusterAddressSet(kFALSE),
+         fTriggerAddressSet(kFALSE),
+         fRecTracksAddressSet(kFALSE),
+         fRunloader(0x0),
+         fMuonloader(0x0),
+         fData(0x0, "MUON", "MUON"),
+         fFilename(),
+         fFoldername(),
+         fEventnumber(-1),
+         fTrack(-1),
+         fSCathode(-1),
+         fCathode(-1)
 {
-// Set all internal pointers to NULL and indices to -1.
+/// Set all internal pointers to 0x0 and indices to -1.
 
        Reset();
-};
-
-AliMUONDataInterface::AliMUONDataInterface(const AliMUONDataInterface& rhs)
-  : TObject(rhs)
-{
-// Protected copy constructor
-
-  AliFatal("Not implemented.");
 }
 
 AliMUONDataInterface::~AliMUONDataInterface()
 {
-// Delete the runloader when done.
-// If the runloader is not to be deleted then call Reset just before 
-// the destructor is called.
+/// 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)
-{
-// Protected assignement operator
-
-  if (this == &rhs) return *this;
-
-  AliFatal("Not implemented.");
-    
-  return *this;  
-}    
-          
+}
 
 void AliMUONDataInterface::Reset()
 {
-// Sets all internal pointers to NULL and indices to -1.
-// Note: No resources are released!
-// Specificaly AliRunLoader is not deleted.
+/// Sets all internal pointers to NULL and indices to -1.
+/// Note: No resources are released!
+/// Specificaly AliRunLoader is not deleted.
 
+       fCreatedRunLoader = kFALSE;
        fRunloader = NULL;
        fMuonloader = NULL;
        fEventnumber = -1;
@@ -75,59 +119,107 @@ void AliMUONDataInterface::Reset()
        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)
 {
-// Fetch the run loader and muon loader objects from memory if they already exist,
-// or from memory if they do not. 
-// If the currently loaded run loader (if any) is not refering to the file and folder
-// we are interested in then it is deleted and reopened with the required file and
-// folder names.
+/// Fetch the run loader and muon loader objects from memory if they already exist,
+/// or from memory if they do not. 
+/// If the currently loaded run loader (if any) is not refering to the file and folder
+/// we are interested in then it is deleted and reopened with the required file and
+/// folder names.
 
        if (fRunloader == NULL)
        {
                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.
@@ -135,15 +227,17 @@ Bool_t AliMUONDataInterface::FetchLoaders(TString filename, TString foldername)
        {
                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.
+/// 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)
        {
@@ -156,7 +250,8 @@ Bool_t AliMUONDataInterface::FetchEvent(Int_t event)
                fDigitAddressSet = kFALSE;
                fClusterAddressSet = kFALSE;
                fTriggerAddressSet = kFALSE;
-       };
+               fRecTracksAddressSet = kFALSE;
+       }
        if ( event != fEventnumber )
        {
                if ( fRunloader->GetEvent(event) < 0 ) return kFALSE;
@@ -169,14 +264,15 @@ Bool_t AliMUONDataInterface::FetchEvent(Int_t event)
                fDigitAddressSet = kFALSE;
                fClusterAddressSet = kFALSE;
                fTriggerAddressSet = kFALSE;
-       };
+               fRecTracksAddressSet = kFALSE;
+       }
        return kTRUE;
-};
+}
 
 
 Bool_t AliMUONDataInterface::FetchTreeK()
 {
-// Fetch the Kine tree from the current run loader.
+/// Fetch the Kine tree from the current run loader.
 
        if (fRunloader->TreeK() == NULL)
        {
@@ -185,16 +281,16 @@ Bool_t AliMUONDataInterface::FetchTreeK()
                {
                        AliError("Could not load TreeK.");
                        return kFALSE;
-               };
-       };
+               }
+       }
        return kTRUE;
-};
+}
 
 
 Bool_t AliMUONDataInterface::FetchTreeH()
 {
-// Fetch the Hits tree from the current muon loader.
-// Set all the required addresses etc...
+/// Fetch the Hits tree from the current muon loader.
+/// Set all the required addresses etc...
 
        if (fMuonloader->TreeH() == NULL)
        {
@@ -203,7 +299,7 @@ Bool_t AliMUONDataInterface::FetchTreeH()
                {
                        AliError("Could not load TreeH.");
                        return kFALSE;
-               };
+               }
                fData.SetTreeAddress("H");
                fHitAddressSet = kTRUE;
        }
@@ -211,15 +307,15 @@ Bool_t AliMUONDataInterface::FetchTreeH()
        {
                fData.SetTreeAddress("H");
                fHitAddressSet = kTRUE;
-       };
+       }
        return kTRUE;
-};
+}
 
 
 Bool_t AliMUONDataInterface::FetchTreeS()
 {
-// Fetch the S-Digits tree from the current muon loader.
-// Set all the required addresses etc...
+/// Fetch the S-Digits tree from the current muon loader.
+/// Set all the required addresses etc...
 
        if (fMuonloader->TreeS() == NULL)
        {
@@ -228,7 +324,7 @@ Bool_t AliMUONDataInterface::FetchTreeS()
                {
                        AliError("Could not load TreeS.");
                        return kFALSE;
-               };
+               }
                fData.SetTreeAddress("S");
                fSDigitAddressSet = kTRUE;
        }
@@ -236,15 +332,15 @@ Bool_t AliMUONDataInterface::FetchTreeS()
        {
                fData.SetTreeAddress("S");
                fSDigitAddressSet = kTRUE;
-       };
+       }
        return kTRUE;
-};
+}
 
 
 Bool_t AliMUONDataInterface::FetchTreeD()
 {
-// Fetch the digits tree from the current muon loader.
-// Set all the required addresses etc...
+/// Fetch the digits tree from the current muon loader.
+/// Set all the required addresses etc...
 
        if (fMuonloader->TreeD() == NULL)
        {
@@ -253,7 +349,7 @@ Bool_t AliMUONDataInterface::FetchTreeD()
                {
                        AliError("Could not load TreeD.");
                        return kFALSE;
-               };
+               }
                fData.SetTreeAddress("D");
                fDigitAddressSet = kTRUE;
        }
@@ -261,60 +357,79 @@ Bool_t AliMUONDataInterface::FetchTreeD()
        {
                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.
+/// 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)
 {
-// Returns the number of events in the specified file/folder, and -1 on error.
+/// Returns the number of events in the specified file/folder, and -1 on error.
 
        if ( ! FetchLoaders(filename, foldername) ) return -1;
        if ( ! FetchEvent(event) ) return -1;
        if ( ! FetchTreeK() ) return -1;
        return (Int_t) fRunloader->TreeK()->GetEntriesFast();
-};
+}
 
 
 TParticle* AliMUONDataInterface::Particle(
                TString filename, TString foldername, Int_t event, Int_t particle
        )
 {
-// Returns the specified particle in the given file, folder and event.
-// NULL is returned on error.
+/// Returns the specified particle in the given file, folder and event.
+/// NULL is returned on error.
 
        if ( ! FetchLoaders(filename, foldername) ) return NULL;
        if ( ! FetchEvent(event) ) return NULL;
@@ -325,27 +440,27 @@ TParticle* AliMUONDataInterface::Particle(
        treeK->GetBranch("Particles")->SetAddress(&p);
        treeK->GetEvent(particle);
        return p;
-};
+}
 
 
 Int_t AliMUONDataInterface::NumberOfTracks(TString filename, TString foldername, Int_t event)
 {
-// Returns the number of tracks in the specified file/folder and event.
-// -1 is returned on error.
+/// Returns the number of tracks in the specified file/folder and event.
+/// -1 is returned on error.
 
        if ( ! FetchLoaders(filename, foldername) ) return -1;
        if ( ! FetchEvent(event) ) return -1;
        if ( ! FetchTreeH() ) return -1;
        return fData.GetNtracks();
-};
+}
 
 
 Int_t AliMUONDataInterface::NumberOfHits(
                TString filename, TString foldername, Int_t event, Int_t track
        )
 {
-// Returns the number of hits in the specified file/folder, event and track.
-// -1 is returned on error.
+/// Returns the number of hits in the specified file/folder, event and track.
+/// -1 is returned on error.
 
        if ( ! FetchLoaders(filename, foldername) ) return -1;
        if ( ! FetchEvent(event) ) return -1;
@@ -356,9 +471,9 @@ Int_t AliMUONDataInterface::NumberOfHits(
                fData.ResetHits();
                fData.GetTrack(track);
                fTrack = track;
-       };
+       }
        return fData.Hits()->GetEntriesFast();
-};
+}
 
 
 AliMUONHit* AliMUONDataInterface::Hit(
@@ -366,8 +481,8 @@ AliMUONHit* AliMUONDataInterface::Hit(
                Int_t track, Int_t hit
        )
 {
-// Returns the specified hit in the given file, folder, event and track.
-// NULL is returned on error.
+/// Returns the specified hit in the given file, folder, event and track.
+/// NULL is returned on error.
 
        if ( ! FetchLoaders(filename, foldername) ) return NULL;
        if ( ! FetchEvent(event) ) return NULL;
@@ -378,9 +493,9 @@ AliMUONHit* AliMUONDataInterface::Hit(
                fData.ResetHits();
                fData.GetTrack(track);
                fTrack = track;
-       };
+       }
        return static_cast<AliMUONHit*>( fData.Hits()->At(hit) );
-};
+}
 
 
 Int_t AliMUONDataInterface::NumberOfSDigits(
@@ -388,11 +503,11 @@ Int_t AliMUONDataInterface::NumberOfSDigits(
                Int_t chamber, Int_t cathode
        )
 {
-// Returns the number of s-digits in the given file, folder, event,
-// chamber and cathode. -1 is returned on error.
+/// Returns the number of s-digits in the given file, folder, event,
+/// chamber and cathode. -1 is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
-       Assert( 0 <= cathode && cathode <= 1 );
+       assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= cathode && cathode <= 1 );
        
        if ( ! FetchLoaders(filename, foldername) ) return -1;
        if ( ! FetchEvent(event) ) return -1;
@@ -401,11 +516,11 @@ Int_t AliMUONDataInterface::NumberOfSDigits(
        if ( fSCathode != cathode )
        {
                fData.ResetSDigits();
-               fData.GetCathodeS(cathode);
+               fData.GetSDigits();
                fSCathode = cathode;
-       };
+       }
        return fData.SDigits(chamber)->GetEntriesFast();
-};
+}
 
 
 AliMUONDigit* AliMUONDataInterface::SDigit(
@@ -413,11 +528,11 @@ AliMUONDigit* AliMUONDataInterface::SDigit(
                Int_t chamber, Int_t cathode, Int_t sdigit
        )
 {
-// Returns the specified s-digit in the given file, folder, event,
-// chamber and cathode. NULL is returned on error.
+/// Returns the specified s-digit in the given file, folder, event,
+/// chamber and cathode. NULL is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
-       Assert( 0 <= cathode && cathode <= 1 );
+       assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= cathode && cathode <= 1 );
        
        if ( ! FetchLoaders(filename, foldername) ) return NULL;
        if ( ! FetchEvent(event) ) return NULL;
@@ -426,11 +541,11 @@ 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(
@@ -438,10 +553,10 @@ Int_t AliMUONDataInterface::NumberOfDigits(
                Int_t chamber, Int_t cathode
        )
 {
-// Returns the number of digits in the given file, folder, event,
-// chamber and cathode. -1 is returned on error.
-       Assert( 0 <= chamber && chamber <= 13 );
-       Assert( 0 <= cathode && cathode <= 1 );
+/// Returns the number of digits in the given file, folder, event,
+/// chamber and cathode. -1 is returned on error.
+       assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= cathode && cathode <= 1 );
        
        if ( ! FetchLoaders(filename, foldername) ) return -1;
        if ( ! FetchEvent(event) ) return -1;
@@ -450,11 +565,11 @@ Int_t AliMUONDataInterface::NumberOfDigits(
        if ( fCathode != cathode )
        {
                fData.ResetDigits();
-               fData.GetCathode(cathode);
+               fData.GetDigits();
                fCathode = cathode;
-       };
+       }
        return fData.Digits(chamber)->GetEntriesFast();
-};
+}
 
 
 AliMUONDigit* AliMUONDataInterface::Digit(
@@ -462,11 +577,11 @@ AliMUONDigit* AliMUONDataInterface::Digit(
                Int_t chamber, Int_t cathode, Int_t digit
        )
 {
-// Returns the specified digit in the given file, folder, event,
-// chamber and cathode. NULL is returned on error.
+/// Returns the specified digit in the given file, folder, event,
+/// chamber and cathode. NULL is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
-       Assert( 0 <= cathode && cathode <= 1 );
+       assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= cathode && cathode <= 1 );
        
        if ( ! FetchLoaders(filename, foldername) ) return NULL;
        if ( ! FetchEvent(event) ) return NULL;
@@ -475,21 +590,21 @@ 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(
                TString filename, TString foldername, Int_t event, Int_t chamber
        )
 {
-// Returns the number of raw clusters in the specified file, folder, event and chamber.
-// -1 is returned or error.
+/// Returns the number of raw clusters in the specified file, folder, event and chamber.
+/// -1 is returned or error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= chamber && chamber <= 13 );
        if ( ! FetchLoaders(filename, foldername) ) return -1;
        if ( ! FetchEvent(event) ) return -1;
        if ( ! FetchTreeR() ) return -1;
@@ -500,9 +615,9 @@ Int_t AliMUONDataInterface::NumberOfRawClusters(
                fData.ResetRawClusters();
                fData.GetRawClusters();
                fClusterAddressSet = kTRUE;
-       };
+       }
        return fData.RawClusters(chamber)->GetEntriesFast();
-};
+}
 
 
 AliMUONRawCluster* AliMUONDataInterface::RawCluster(
@@ -510,10 +625,10 @@ AliMUONRawCluster* AliMUONDataInterface::RawCluster(
                Int_t chamber, Int_t cluster
        )
 {
-// Fetch the specified raw cluster from the given file, folder, event and chamber number.
-// NULL is returned on error.
+/// Fetch the specified raw cluster from the given file, folder, event and chamber number.
+/// NULL is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= chamber && chamber <= 13 );
        if ( ! FetchLoaders(filename, foldername) ) return NULL;
        if ( ! FetchEvent(event) ) return NULL;
        if ( ! FetchTreeR() ) return NULL;
@@ -524,333 +639,338 @@ 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)
 {
-// Return the number of local trigger objects in the specified file, folder and
-// event number. -1 is returned on error.
+/// Return the number of local trigger objects in the specified file, folder and
+/// event number. -1 is returned on error.
 
        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(
                TString filename, TString foldername, Int_t event, Int_t trigger
        )
 {
-// Fetch the specified local trigger object from the given file, folder and event number.
-// NULL is returned on error.
+/// Fetch the specified local trigger object from the given file, folder and event number.
+/// NULL is returned on error.
 
        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)
 {
-// Set the current file and folder from which to fetch data.
-// kTRUE is returned if the run and muon loaders were found, else kFALSE. 
+/// Set the current file and folder from which to fetch data.
+/// 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.
+/// 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()
 {
-// Get the number of events in the currently selected file.
-// -1 is returned on error.
+/// Get the number of events in the currently selected file.
+/// -1 is returned on error.
 
        if (fRunloader == NULL)
        {
                AliError("File not set.");
                return -1;
-       };
+       }
        return fRunloader->GetNumberOfEvents();
-};
+}
 
 
 Int_t AliMUONDataInterface::NumberOfParticles()
 {
-// Get the number of particles in the current event.
-// -1 is returned on error.
+/// Get the number of particles in the current event.
+/// -1 is returned on error.
 
        if (fRunloader == NULL)
        {
                AliError("File not set.");
                return -1;
-       };
+       }
        if ( ! FetchTreeK() ) return -1;
        return (Int_t) fRunloader->TreeK()->GetEntriesFast();
-};
+}
 
 
 TParticle* AliMUONDataInterface::Particle(Int_t particle)
 {
-// Fetch the specified particle from the current event.
-// NULL is returned on error.
+/// Fetch the specified particle 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 ( ! FetchTreeK() ) return NULL;
        TTree* treeK = fRunloader->TreeK();
        TParticle* p = NULL;
        treeK->GetBranch("Particles")->SetAddress(&p);
        treeK->GetEvent(particle);
        return p;
-};
+}
 
 
 Int_t AliMUONDataInterface::NumberOfTracks()
 {
-// Get the number of tracks in the current event.
-// -1 is returned on error.
+/// Get the number of tracks 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 ( ! FetchTreeH() ) return -1;
        return fData.GetNtracks();
-};
+}
 
 
 Int_t AliMUONDataInterface::NumberOfHits(Int_t track)
 {
-// Get the number of hits for the given track in the current event.
-// -1 is returned on error.
+/// Get the number of hits for the given track 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 ( ! 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)
 {
-// Fetch the specified hit from the current event.
-// NULL is returned on error.
+/// Fetch the specified hit 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 ( ! 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)
 {
-// Get the number of s-digits on the chamber, cathode in the current event.
-// -1 is returned on error.
+/// Get the number of s-digits on the chamber, cathode in the current event.
+/// -1 is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
-       Assert( 0 <= cathode && cathode <= 1 );
+       assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= cathode && cathode <= 1 );
        
        if (fRunloader == NULL)
        {
                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)
 {
-// Fetch the specified s-digits on the chamber, cathode from the current event.
-// NULL is returned on error.
+/// Fetch the specified s-digits on the chamber, cathode from the current event.
+/// NULL is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
-       Assert( 0 <= cathode && cathode <= 1 );
+       assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= cathode && cathode <= 1 );
        
        if (fRunloader == NULL)
        {
                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)
 {
-// Get the number of digits on the chamber, cathode in the current event.
-// -1 is returned on error.
+/// Get the number of digits on the chamber, cathode in the current event.
+/// -1 is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
-       Assert( 0 <= cathode && cathode <= 1 );
+       assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= cathode && cathode <= 1 );
        
        if (fRunloader == NULL)
        {
                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)
 {
-// Fetch the specified digits on the chamber, cathode from the current event.
-// NULL is returned on error.
+/// Fetch the specified digits on the chamber, cathode from the current event.
+/// NULL is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
-       Assert( 0 <= cathode && cathode <= 1 );
+       assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= cathode && cathode <= 1 );
        
        if (fRunloader == NULL)
        {
                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)
 {
-// Get the number of raw clusters on the given chamber in the current event.
-// -1 is returned on error.
+/// Get the number of raw clusters on the given chamber in the current event.
+/// -1 is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= chamber && chamber <= 13 );
 
        if (fRunloader == NULL)
        {
                AliError("File not set.");
                return -1;
-       };
+       }
        if (fEventnumber < 0)
        {
                AliError("Event not chosen.");
                return -1;
-       };
+       }
 
        if ( ! FetchTreeR() ) return -1;
        if ( ! fClusterAddressSet )
@@ -859,28 +979,28 @@ Int_t AliMUONDataInterface::NumberOfRawClusters(Int_t chamber)
                fData.ResetRawClusters();
                fData.GetRawClusters();
                fClusterAddressSet = kTRUE;
-       };
+       }
        return fData.RawClusters(chamber)->GetEntriesFast();
-};
+}
 
 
 AliMUONRawCluster* AliMUONDataInterface::RawCluster(Int_t chamber, Int_t cluster)
 {
-// Fetch the specified raw cluster on the given chamber from the current event.
-// NULL is returned on error.
+/// Fetch the specified raw cluster on the given chamber from the current event.
+/// NULL is returned on error.
 
-       Assert( 0 <= chamber && chamber <= 13 );
+       assert( 0 <= chamber && chamber <= 13 );
 
        if (fRunloader == NULL)
        {
                AliError("File not set.");
                return NULL;
-       };
+       }
        if (fEventnumber < 0)
        {
                AliError("Event not chosen.");
                return NULL;
-       };
+       }
 
        if ( ! FetchTreeR() ) return NULL;
        if ( ! fClusterAddressSet )
@@ -889,62 +1009,171 @@ AliMUONRawCluster* AliMUONDataInterface::RawCluster(Int_t chamber, Int_t cluster
                fData.ResetRawClusters();
                fData.GetRawClusters();
                fClusterAddressSet = kTRUE;
-       };
+       }
        return static_cast<AliMUONRawCluster*>( fData.RawClusters(chamber)->At(cluster) );
-};
+}
 
 
 Int_t AliMUONDataInterface::NumberOfLocalTriggers()
 {
-// Get the number of local trigger objects in the current event.
-// -1 is returned on error.
+/// 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 ( ! 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)
 {
-// Fetch the specified local trigger object from the current event.
-// NULL is returned on error.
+/// 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 ( ! 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));
+}