* provided "as is" without express or implied warranty. *
**************************************************************************/
-// $Id: AliEvent.cxx,v 1.5 2002/01/18 08:46:59 nick Exp $
+// $Id: AliEvent.cxx,v 1.22 2004/06/29 11:29:37 nick Exp $
///////////////////////////////////////////////////////////////////////////
// Class AliEvent
// Creation and investigation of an Alice physics event.
// An AliEvent can be constructed by adding AliTracks, Alivertices, AliJets
-// and/or AliCalorimeters.
+// and/or devices like AliCalorimeters or AliDevice (derived) objects.
+//
+// All objects which are derived from TObject can be regarded as a device.
+// However, AliDevice (or derived) objects profit from additional hit
+// handling facilities.
+// A "hit" is a generic name indicating an AliSignal (or derived) object.
+// Note that AliEvent does NOT own hits; it only provides references to hits
+// obtained from the various devices.
+// This implies that hits should be owned by the devices themselves.
//
// The basic functionality of AliEvent is identical to the one of AliVertex.
// So, an AliEvent may be used as the primary vertex with some additional
//
// To provide maximal flexibility to the user, the two modes of track/jet/vertex
// storage as described in AliJet and AliVertex can be used.
-// In addition an identical structure is provided for the storage of AliCalorimeter
-// objects, which can be selected by means of the memberfunction SetCalCopy().
-//
-// a) SetCalCopy(0) (which is the default).
-// Only the pointers of the 'added' calorimeters are stored.
-// This mode is typically used by making cal. studies based on a fixed set
-// of calorimeters which stays under user control or is kept on an external
+// In addition an identical structure is provided for the storage of devices like
+// AliCalorimeter objects, which can be selected by means of the memberfunction
+// SetDevCopy().
+//
+// a) SetDevCopy(0) (which is the default).
+// Only the pointers of the 'added' devices are stored.
+// This mode is typically used by making studies based on a fixed set
+// of devices which stays under user control or is kept on an external
// file/tree.
// In this way the AliEvent just represents a 'logical structure' for the
// physics analysis.
//
// Note :
-// Modifications made to the original calorimeters also affect the AliCalorimeter
+// Modifications made to the original devices also affect the device
// objects which are stored in the AliEvent.
//
-// b) SetCalCopy(1).
-// Of every 'added' calorimeter a private copy will be made of which the pointer
+// b) SetDevCopy(1).
+// Of every 'added' device a private copy will be made of which the pointer
// will be stored.
// In this way the AliEvent represents an entity on its own and modifications
// made to the original calorimeters do not affect the AliCalorimeter objects
// which are stored in the AliEvent.
-// This mode will allow 'adding' many different AliCalorimeters into an AliEvent by
-// creating only one AliCalorimeter instance in the main programme and using the
-// AliCalorimeter::Reset() and AliCalorimeter parameter setting memberfunctions.
+// This mode will allow 'adding' many different devices into an AliEvent by
+// creating only one device instance in the main programme and using the
+// Reset() and parameter setting memberfunctions of the object representing the device.
+//
+// Note :
+// The copy is made using the Clone() memberfunction.
+// All devices (i.e. classes derived from TObject) have the default TObject::Clone()
+// memberfunction.
+// However, devices generally contain an internal (signal) data structure
+// which may include pointers to other objects. Therefore it is recommended to provide
+// for all devices a specific copy constructor and override the default Clone()
+// memberfunction using this copy constructor.
+// Examples for this may be seen from AliCalorimeter, AliSignal and AliDevice.
+//
+// See also the documentation provided for the memberfunction SetOwner().
//
// Coding example to make an event consisting of a primary vertex,
// 2 secondary vertices and a calorimeter.
// Specify the event object as the repository of all objects
// for the event building and physics analysis.
//
-// evt.SetCalCopy(1);
+// evt.SetDevCopy(1);
// evt.SetTrackCopy(1);
//
// Fill the event structure with the basic objects
//
-// AliCalorimeter emcal;
+// AliCalorimeter emcal1;
+// AliCalorimeter emcal2;
// ...
-// ... // code to fill the calorimeter data
+// ... // code to fill the emcal1 and emcal2 calorimeter data
// ...
//
-// evt.AddCalorimeter(emcal);
+// evt.AddDevice(emcal1);
+// evt.AddDevice(emcal2);
+//
+// // Assume AliTOF has been derived from AliDevice
+// AliTOF tof1;
+// AliTOF tof2;
+// ...
+// ... // code to fill the tof1 and tof2 data
+// ...
+//
+// evt.AddDevice(tof1);
+// evt.AddDevice(tof2);
//
// AliTrack* tx=new AliTrack();
// for (Int_t i=0; i<10; i++)
// tx=0;
// }
//
+// Order and investigate all the hits of all the TOF devices
+//
+// TObjArray* hits=evt.GetHits("AliTOF");
+// TObjArray orderedtofs=evt.SortHits(hits);
+// Int_t nhits=orderedtofs.GetEntries();
+// for (Int_t i=0; i<nhits; i++)
+// {
+// AliSignal* sx=(AliSignal*)orderedtofs.At(i);
+// if (sx) sx->Data();
+// }
+//
+// Order and investigate all the hits of all the calorimeter devices
+//
+// TObjArray* hits=evt.GetHits("AliCalorimeter");
+// TObjArray orderedcals=evt.SortHits(hits);
+// Int_t nhits=orderedcals.GetEntries();
+// for (Int_t i=0; i<nhits; i++)
+// {
+// AliSignal* sx=(AliSignal*)orderedcals.At(i);
+// if (sx) sx->Data();
+// }
+//
// Build the event structure (vertices, jets, ...) for physics analysis
// based on the basic objects from the event repository.
//
// jx=v2.GetJet(2);
// evt.AddJet(jx,0);
//
-// evt.Info("sph");
+// evt.Data("sph");
// v1.ListAll();
// v2.List("cyl");
//
// Float_t loc[3];
// evt.GetPosition(loc,"sph");
// AliPosition r=v1.GetPosition();
-// r.Info();
+// r.Data();
// Int_t nt=v2.GetNtracks();
// AliTrack* tv=v2.GetTrack(1); // Access track number 1 of Vertex v2
//
// Note : All quantities are in GeV, GeV/c or GeV/c**2
//
//--- Author: Nick van Eijndhoven 27-may-2001 UU-SAP Utrecht
-//- Modified: NvE $Date: 2002/01/18 08:46:59 $ UU-SAP Utrecht
+//- Modified: NvE $Date: 2004/06/29 11:29:37 $ UU-SAP Utrecht
///////////////////////////////////////////////////////////////////////////
#include "AliEvent.h"
+#include "Riostream.h"
ClassImp(AliEvent) // Class implementation to enable ROOT I/O
-AliEvent::AliEvent()
+AliEvent::AliEvent() : AliVertex()
{
// Default constructor.
// All variables initialised to default values.
fAproj=0;
fZproj=0;
fPnucProj=0;
+ fIdProj=0;
fAtarg=0;
fZtarg=0;
fPnucTarg=0;
- fNcals=0;
- fCalorimeters=0;
- fCalCopy=0;
+ fIdTarg=0;
+ fDevices=0;
+ fDevCopy=0;
+ fHits=0;
}
///////////////////////////////////////////////////////////////////////////
-AliEvent::AliEvent(Int_t n): AliVertex(n)
+AliEvent::AliEvent(Int_t n) : AliVertex(n)
{
// Create an event to hold initially a maximum of n tracks
// All variables initialised to default values
+ if (n<=0)
+ {
+ cout << " *** This AliVertex initialisation was invoked via the AliEvent ctor." << endl;
+ }
fDaytime.Set();
fRun=0;
fEvent=0;
fAproj=0;
fZproj=0;
fPnucProj=0;
+ fIdProj=0;
fAtarg=0;
fZtarg=0;
fPnucTarg=0;
- fNcals=0;
- fCalorimeters=0;
- fCalCopy=0;
+ fIdTarg=0;
+ fDevices=0;
+ fDevCopy=0;
+ fHits=0;
}
///////////////////////////////////////////////////////////////////////////
AliEvent::~AliEvent()
{
// Default destructor
- if (fCalorimeters)
+ if (fDevices)
+ {
+ delete fDevices;
+ fDevices=0;
+ }
+ if (fHits)
+ {
+ delete fHits;
+ fHits=0;
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+AliEvent::AliEvent(const AliEvent& evt) : AliVertex(evt)
+{
+// Copy constructor.
+ fDaytime=evt.fDaytime;
+ fRun=evt.fRun;
+ fEvent=evt.fEvent;
+ fAproj=evt.fAproj;
+ fZproj=evt.fZproj;
+ fPnucProj=evt.fPnucProj;
+ fIdProj=evt.fIdProj;
+ fAtarg=evt.fAtarg;
+ fZtarg=evt.fZtarg;
+ fPnucTarg=evt.fPnucTarg;
+ fIdTarg=evt.fIdTarg;
+ fDevCopy=evt.fDevCopy;
+
+ fDevices=0;
+ Int_t ndevs=evt.GetNdevices();
+ if (ndevs)
+ {
+ fDevices=new TObjArray(ndevs);
+ if (fDevCopy) fDevices->SetOwner();
+ for (Int_t i=1; i<=ndevs; i++)
+ {
+ TObject* dev=evt.GetDevice(i);
+ if (dev)
+ {
+ if (fDevCopy)
+ {
+ fDevices->Add(dev->Clone());
+ }
+ else
+ {
+ fDevices->Add(dev);
+ }
+ }
+ }
+ }
+
+ fHits=0;
+ if (evt.fHits)
{
- delete fCalorimeters;
- fCalorimeters=0;
+ Int_t nhits=evt.fHits->GetEntries();
+ if (nhits)
+ {
+ fHits=new TObjArray(nhits);
+ for (Int_t ih=0; ih<nhits; ih++)
+ {
+ AliSignal* sx=(AliSignal*)evt.fHits->At(ih);
+ fHits->Add(sx);
+ }
+ }
}
}
///////////////////////////////////////////////////////////////////////////
// Reset all variables to default values
// The max. number of tracks is set to the initial value again
// The max. number of vertices is set to the default value again
+// Note : The CalCopy mode is maintained as it was set by the user before.
+
+ AliVertex::Reset();
+
fDaytime.Set();
fRun=0;
fEvent=0;
fAproj=0;
fZproj=0;
fPnucProj=0;
+ fIdProj=0;
fAtarg=0;
fZtarg=0;
fPnucTarg=0;
+ fIdTarg=0;
- fNcals=0;
- if (fCalorimeters)
+ if (fDevices)
{
- delete fCalorimeters;
- fCalorimeters=0;
+ delete fDevices;
+ fDevices=0;
}
+ if (fHits)
+ {
+ delete fHits;
+ fHits=0;
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+void AliEvent::SetOwner(Bool_t own)
+{
+// Set ownership of all added objects.
+// The default parameter is own=kTRUE.
+//
+// Invokation of this memberfunction also sets all the copy modes
+// (e.g. TrackCopy & co.) according to the value of own.
+//
+// This function (with own=kTRUE) is particularly useful when reading data
+// from a tree/file, since Reset() will then actually remove all the
+// added objects from memory irrespective of the copy mode settings
+// during the tree/file creation process. In this way it provides a nice way
+// of preventing possible memory leaks in the reading/analysis process.
+//
+// In addition this memberfunction can also be used as a shortcut to set all
+// copy modes in one go during a tree/file creation process.
+// However, in this case the user has to take care to only set/change the
+// ownership (and copy mode) for empty objects (e.g. newly created objects
+// or after invokation of the Reset() memberfunction) otherwise it will
+// very likely result in inconsistent destructor behaviour.
- AliVertex::Reset();
+ Int_t mode=1;
+ if (!own) mode=0;
+ if (fDevices) fDevices->SetOwner(own);
+ fDevCopy=mode;
+
+ AliVertex::SetOwner(own);
}
///////////////////////////////////////////////////////////////////////////
-void AliEvent::SetDayTime(TDatime& stamp)
+void AliEvent::SetDayTime(TTimeStamp& stamp)
{
-// Set the date and time stamp for this event
+// Set the date and time stamp for this event.
+// An exact copy of the entered date/time stamp will be saved with an
+// accuracy of 1 nanosecond.
fDaytime=stamp;
}
///////////////////////////////////////////////////////////////////////////
+void AliEvent::SetDayTime(TDatime& stamp)
+{
+// Set the date and time stamp for this event.
+// The entered date/time will be interpreted as being the local date/time
+// and the accuracy is 1 second.
+// This function with the TDatime argument is mainly kept for backward
+// compatibility reasons. It is recommended to use the corresponding
+// function with the TTimeStamp argument.
+
+ fDaytime.Set(stamp.GetDate(),stamp.GetTime(),0,kFALSE,0);
+}
+///////////////////////////////////////////////////////////////////////////
void AliEvent::SetRunNumber(Int_t run)
{
// Set the run number for this event
fEvent=evt;
}
///////////////////////////////////////////////////////////////////////////
-TDatime AliEvent::GetDayTime()
+TTimeStamp AliEvent::GetDayTime() const
{
// Provide the date and time stamp for this event
return fDaytime;
}
///////////////////////////////////////////////////////////////////////////
-Int_t AliEvent::GetRunNumber()
+Int_t AliEvent::GetRunNumber() const
{
// Provide the run number for this event
return fRun;
}
///////////////////////////////////////////////////////////////////////////
-Int_t AliEvent::GetEventNumber()
+Int_t AliEvent::GetEventNumber() const
{
// Provide the event number for this event
return fEvent;
}
///////////////////////////////////////////////////////////////////////////
-void AliEvent::SetProjectile(Int_t a,Int_t z,Double_t pnuc)
+void AliEvent::SetProjectile(Int_t a,Int_t z,Double_t pnuc,Int_t id)
{
-// Set the projectile A, Z and momentum value per nucleon.
+// Set the projectile A, Z, momentum per nucleon and user defined particle ID.
+// By default the particle ID is set to zero.
fAproj=a;
fZproj=z;
fPnucProj=pnuc;
+ fIdProj=id;
}
///////////////////////////////////////////////////////////////////////////
-Int_t AliEvent::GetProjectileA()
+Int_t AliEvent::GetProjectileA() const
{
// Provide the projectile A value.
return fAproj;
}
///////////////////////////////////////////////////////////////////////////
-Int_t AliEvent::GetProjectileZ()
+Int_t AliEvent::GetProjectileZ() const
{
// Provide the projectile Z value.
return fZproj;
}
///////////////////////////////////////////////////////////////////////////
-Double_t AliEvent::GetProjectilePnuc()
+Double_t AliEvent::GetProjectilePnuc() const
{
// Provide the projectile momentum value per nucleon.
return fPnucProj;
}
///////////////////////////////////////////////////////////////////////////
-void AliEvent::SetTarget(Int_t a,Int_t z,Double_t pnuc)
+Int_t AliEvent::GetProjectileId() const
+{
+// Provide the user defined particle ID of the projectile.
+ return fIdProj;
+}
+///////////////////////////////////////////////////////////////////////////
+void AliEvent::SetTarget(Int_t a,Int_t z,Double_t pnuc,Int_t id)
{
-// Set the target A, Z and momentum value per nucleon.
+// Set the target A, Z, momentum per nucleon and user defined particle ID.
+// By default the particle ID is set to zero.
fAtarg=a;
fZtarg=z;
fPnucTarg=pnuc;
+ fIdTarg=id;
}
///////////////////////////////////////////////////////////////////////////
-Int_t AliEvent::GetTargetA()
+Int_t AliEvent::GetTargetA() const
{
// Provide the target A value.
return fAtarg;
}
///////////////////////////////////////////////////////////////////////////
-Int_t AliEvent::GetTargetZ()
+Int_t AliEvent::GetTargetZ() const
{
// Provide the target Z value.
return fZtarg;
}
///////////////////////////////////////////////////////////////////////////
-Double_t AliEvent::GetTargetPnuc()
+Double_t AliEvent::GetTargetPnuc() const
{
// Provide the target momentum value per nucleon.
return fPnucTarg;
}
///////////////////////////////////////////////////////////////////////////
-void AliEvent::HeaderInfo()
+Int_t AliEvent::GetTargetId() const
+{
+// Provide the user defined particle ID of the target.
+ return fIdTarg;
+}
+///////////////////////////////////////////////////////////////////////////
+void AliEvent::HeaderData() const
{
// Provide event header information
- Int_t date=fDaytime.GetDate();
- Int_t time=fDaytime.GetTime();
+ const char* name=GetName();
+ const char* title=GetTitle();
+ Int_t ndevs=GetNdevices();
+ cout << " *" << ClassName() << "::Data*";
+ if (strlen(name)) cout << " Name : " << GetName();
+ if (strlen(title)) cout << " Title : " << GetTitle();
+ cout << endl;
+ cout << " " << fDaytime.AsString() << endl;
+ cout << " Run : " << fRun << " Event : " << fEvent
+ << " Number of devices : " << ndevs << endl;
- Int_t year=date/10000;
- Int_t month=(date%10000)/100;
- Int_t day=date%100;
- Int_t hh=time/10000;
- Int_t mm=(time%10000)/100;
- Int_t ss=time%100;
-
- char* c[12]={"jan","feb","mar","apr","may","jun",
- "jul","aug","sep","oct","nov","dec"};
-
- cout << " *AliEvent::Info* Run : " << fRun << " Event : " << fEvent;
- cout.fill('0');
- cout << " Date : " << setw(2) << day << "-" << c[month-1] << "-" << year
- << " Time : " << setw(2) << hh << ":" << setw(2) << mm << ":" << setw(2) << ss;
- cout.fill(' ');
- cout << " Ncalorimeters : " << fNcals << endl;
+ if (ndevs) ShowDevices();
}
///////////////////////////////////////////////////////////////////////////
-void AliEvent::Info(TString f)
+void AliEvent::Data(TString f)
{
// Provide event information within the coordinate frame f
- HeaderInfo();
- AliVertex::Info(f);
+ HeaderData();
+ AliVertex::Data(f);
}
///////////////////////////////////////////////////////////////////////////
-Int_t AliEvent::GetNcalorimeters()
+Int_t AliEvent::GetNdevices() const
{
-// Provide the number of stored calorimeter systems
- return fNcals;
+// Provide the number of stored devices
+ Int_t ndevs=0;
+ if (fDevices) ndevs=fDevices->GetEntries();
+ return ndevs;
}
///////////////////////////////////////////////////////////////////////////
-void AliEvent::AddCalorimeter(AliCalorimeter& c)
+void AliEvent::AddDevice(TObject& d)
{
-// Add a calorimeter system to the event
- if (!fCalorimeters)
+// Add a device to the event.
+//
+// Note :
+// In case a private copy is made, this is performed via the Clone() memberfunction.
+// All devices (i.e. classes derived from TObject) have the default TObject::Clone()
+// memberfunction.
+// However, devices generally contain an internal (signal) data structure
+// which may include pointers to other objects. Therefore it is recommended to provide
+// for all devices a specific copy constructor and override the default Clone()
+// memberfunction using this copy constructor.
+// An example for this may be seen from AliCalorimeter.
+
+ if (!fDevices)
{
- fCalorimeters=new TObjArray();
- if (fCalCopy) fCalorimeters->SetOwner();
+ fDevices=new TObjArray();
+ if (fDevCopy) fDevices->SetOwner();
}
- // Add the calorimeter system to this event
- fNcals++;
- if (fCalCopy)
+ // Add the device to this event
+ if (fDevCopy)
{
- fCalorimeters->AddLast((AliCalorimeter*)c.Clone());
+ fDevices->Add(d.Clone());
}
else
{
- fCalorimeters->AddLast(&c);
+ fDevices->Add(&d);
}
}
///////////////////////////////////////////////////////////////////////////
-void AliEvent::SetCalCopy(Int_t j)
+void AliEvent::SetDevCopy(Int_t j)
{
-// (De)activate the creation of private copies of the added calorimeters.
-// j=0 ==> No private copies are made; pointers of original cals. are stored.
-// j=1 ==> Private copies of the cals. are made and these pointers are stored.
-//
-// Note : Once the storage contains pointer(s) to AliCalorimeter(s) one cannot
-// change the CalCopy mode anymore.
-// To change the CalCopy mode for an existing AliEvent containing
-// calorimeters one first has to invoke Reset().
- if (!fCalorimeters)
+// (De)activate the creation of private copies of the added devices.
+// j=0 ==> No private copies are made; pointers of original devices are stored.
+// j=1 ==> Private copies of the devices are made and these pointers are stored.
+//
+//
+// Notes :
+// In case a private copy is made, this is performed via the Clone() memberfunction.
+// All devices (i.e. classes derived from TObject) have the default TObject::Clone()
+// memberfunction.
+// However, devices generally contain an internal (signal) data structure
+// which may include pointers to other objects. Therefore it is recommended to provide
+// for all devices a specific copy constructor and override the default Clone()
+// memberfunction using this copy constructor.
+// An example for this may be seen from AliCalorimeter.
+//
+// Once the storage contains pointer(s) to device(s) one cannot
+// change the DevCopy mode anymore.
+// To change the DevCopy mode for an existing AliEvent containing
+// devices one first has to invoke Reset().
+
+ if (!fDevices)
{
if (j==0 || j==1)
{
- fCalCopy=j;
+ fDevCopy=j;
}
else
{
- cout << "*AliEvent::SetCalCopy* Invalid argument : " << j << endl;
+ cout << " *" << ClassName() << "::SetDevCopy* Invalid argument : " << j << endl;
}
}
else
{
- cout << "*AliEvent::SetCalCopy* Storage already contained calorimeters."
- << " ==> CalCopy mode not changed." << endl;
+ cout << " *" << ClassName() << "::SetDevCopy* Storage already contained devices."
+ << " ==> DevCopy mode not changed." << endl;
}
}
///////////////////////////////////////////////////////////////////////////
-Int_t AliEvent::GetCalCopy()
+Int_t AliEvent::GetDevCopy() const
{
-// Provide value of the CalCopy mode.
-// 0 ==> No private copies are made; pointers of original cals. are stored.
-// 1 ==> Private copies of the cals. are made and these pointers are stored.
- return fCalCopy;
+// Provide value of the DevCopy mode.
+// 0 ==> No private copies are made; pointers of original devices are stored.
+// 1 ==> Private copies of the devices are made and these pointers are stored.
+//
+// Note :
+// In case a private copy is made, this is performed via the Clone() memberfunction.
+// All devices (i.e. classes derived from TObject) have the default TObject::Clone()
+// memberfunction.
+// However, devices generally contain an internal (signal) data structure
+// which may include pointers to other objects. Therefore it is recommended to provide
+// for all devices a specific copy constructor and override the default Clone()
+// memberfunction using this copy constructor.
+// An example for this may be seen from AliCalorimeter.
+
+ return fDevCopy;
}
///////////////////////////////////////////////////////////////////////////
-AliCalorimeter* AliEvent::GetCalorimeter(Int_t i)
+TObject* AliEvent::GetDevice(Int_t i) const
{
-// Return the i-th calorimeter of this event
- if (!fCalorimeters)
+// Return the i-th device of this event.
+// The first device corresponds to i=1.
+
+ if (!fDevices)
{
- cout << " *AliEvent::GetCalorimeter* No calorimeters present." << endl;
return 0;
}
else
{
- if (i<=0 || i>fNcals)
+ Int_t ndevs=GetNdevices();
+ if (i<=0 || i>ndevs)
{
- cout << " *AliEvent::GetCalorimeter* Invalid argument i : " << i
- << " Ncals = " << fNcals << endl;
+ cout << " *" << ClassName() << "::GetDevice* Invalid argument i : " << i
+ << " ndevs = " << ndevs << endl;
return 0;
}
else
{
- return (AliCalorimeter*)fCalorimeters->At(i-1);
+ return fDevices->At(i-1);
}
}
}
///////////////////////////////////////////////////////////////////////////
-AliCalorimeter* AliEvent::GetCalorimeter(TString name)
+TObject* AliEvent::GetDevice(TString name) const
{
-// Return the calorimeter with name tag "name"
- if (!fCalorimeters)
+// Return the device with name tag "name"
+ if (!fDevices)
{
- cout << " *AliEvent::GetCalorimeter* No calorimeters present." << endl;
return 0;
}
else
{
- AliCalorimeter* cx;
TString s;
- for (Int_t i=0; i<fNcals; i++)
+ Int_t ndevs=GetNdevices();
+ for (Int_t i=0; i<ndevs; i++)
{
- cx=(AliCalorimeter*)fCalorimeters->At(i);
- if (cx)
+ TObject* dev=fDevices->At(i);
+ if (dev)
{
- s=cx->GetName();
- if (s == name) return cx;
+ s=dev->GetName();
+ if (s == name) return dev;
}
}
}
}
///////////////////////////////////////////////////////////////////////////
+void AliEvent::ShowDevices() const
+{
+// Provide an overview of the available devices.
+ Int_t ndevs=GetNdevices();
+ if (ndevs)
+ {
+ cout << " The following " << ndevs << " devices are available :" << endl;
+ for (Int_t i=1; i<=ndevs; i++)
+ {
+ TObject* dev=GetDevice(i);
+ if (dev)
+ {
+ cout << " Device number : " << i
+ << " Class : " << dev->ClassName()
+ << " Name : " << dev->GetName() << endl;
+ }
+ }
+ }
+ else
+ {
+ cout << " No devices present for this event." << endl;
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+Int_t AliEvent::GetNhits(const char* classname)
+{
+// Provide the number of hits registered to the specified device class.
+// The specified device class has to be derived from AliDevice.
+// It is possible to indicate with the argument "classname" a specific
+// device instead of a whole class of devices. However, in such a case
+// it is more efficient to use the GetDevice() memberfunction directly.
+ LoadHits(classname);
+ Int_t nhits=0;
+ if (fHits) nhits=fHits->GetEntries();
+ return nhits;
+}
+///////////////////////////////////////////////////////////////////////////
+TObjArray* AliEvent::GetHits(const char* classname)
+{
+// Provide the references to all the hits registered to the specified
+// device class.
+// The specified device class has to be derived from AliDevice.
+// It is possible to indicate with the argument "classname" a specific
+// device instead of a whole class of devices. However, in such a case
+// it is more efficient to use the GetDevice() memberfunction directly.
+ LoadHits(classname);
+ return fHits;
+}
+///////////////////////////////////////////////////////////////////////////
+void AliEvent::LoadHits(const char* classname)
+{
+// Load the references to the various hits registered to the specified
+// device class.
+// The specified device class has to be derived from AliDevice.
+ if (fHits) fHits->Clear();
+
+ Int_t ndev=GetNdevices();
+ for (Int_t idev=1; idev<=ndev; idev++)
+ {
+ TObject* obj=GetDevice(idev);
+ if (!obj) continue;
+
+ if (obj->InheritsFrom(classname) && obj->InheritsFrom("AliDevice"))
+ {
+ AliDevice* dev=(AliDevice*)GetDevice(idev);
+ Int_t nhits=dev->GetNhits();
+ if (nhits)
+ {
+ if (!fHits) fHits=new TObjArray();
+ for (Int_t ih=1; ih<=nhits; ih++)
+ {
+ AliSignal* sx=dev->GetHit(ih);
+ if (sx) fHits->Add(sx);
+ }
+ }
+ }
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+TObjArray AliEvent::SortHits(TObjArray* hits,Int_t idx,Int_t mode) const
+{
+// Order the references to an array of hits by looping over the input array "hits"
+// and checking the signal value. The ordered array is returned as a TObjArray.
+// Note that the input array is not modified.
+// A "hit" represents an abstract object which is derived from AliSignal.
+// The user can specify the index of the signal slot to perform the sorting on.
+// By default the slotindex will be 1.
+// Via the "mode" argument the user can specify ordering in decreasing
+// order (mode=-1) or ordering in increasing order (mode=1).
+// The default is mode=-1.
+// Signals which were declared as "Dead" will be rejected.
+// The gain etc... corrected signals will be used in the ordering process.
+
+ if (idx<=0 || abs(mode)!=1 || !hits)
+ {
+ TObjArray ordered;
+ return ordered;
+ }
+ else
+ {
+ AliDevice dev;
+ TObjArray ordered=dev.SortHits(idx,mode,hits);
+ return ordered;
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+TObjArray AliEvent::SortHits(TObjArray* hits,TString name,Int_t mode) const
+{
+// Order the references to an array of hits by looping over the input array "hits"
+// and checking the signal value. The ordered array is returned as a TObjArray.
+// Note that the input array is not modified.
+// A "hit" represents an abstract object which is derived from AliSignal.
+// The user can specify the name of the signal slot to perform the sorting on.
+// In case no matching slotname is found, the signal will be skipped.
+// Via the "mode" argument the user can specify ordering in decreasing
+// order (mode=-1) or ordering in increasing order (mode=1).
+// The default is mode=-1.
+// Signals which were declared as "Dead" will be rejected.
+// The gain etc... corrected signals will be used in the ordering process.
+
+ if (abs(mode)!=1 || !hits)
+ {
+ TObjArray ordered;
+ return ordered;
+ }
+ else
+ {
+ AliDevice dev;
+ TObjArray ordered=dev.SortHits(name,mode,hits);
+ return ordered;
+ }
+}
+///////////////////////////////////////////////////////////////////////////
+TObject* AliEvent::Clone(const char* name) const
+{
+// Make a deep copy of the current object and provide the pointer to the copy.
+// This memberfunction enables automatic creation of new objects of the
+// correct type depending on the object type, a feature which may be very useful
+// for containers when adding objects in case the container owns the objects.
+// This feature allows to store either AliEvent objects or objects derived from
+// AliEvent via some generic AddEvent memberfunction, provided these derived
+// classes also have a proper Clone memberfunction.
+
+ AliEvent* evt=new AliEvent(*this);
+ if (name)
+ {
+ if (strlen(name)) evt->SetName(name);
+ }
+ return evt;
+}
+///////////////////////////////////////////////////////////////////////////