]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliESD.h
Use default errors in case the vertexer didn't find any
[u/mrichter/AliRoot.git] / STEER / AliESD.h
index d0fbae66e4b987ae61e019e7e86565937657b162..44fccd1bf7e862a8177403b48a5a5c571da14a46 100644 (file)
+// -*- mode: C++ -*- 
 #ifndef ALIESD_H
 #define ALIESD_H
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
-/* $Id$ */
-
-///////////////////////////////////////////////////////////////////////////////
-//                                                                           //
-//   Alice ESD object                                                        //
-//                                                                           //
-///////////////////////////////////////////////////////////////////////////////
-
-#include "TObject.h"
-#include "TObjArray.h"
-#include "TArrayI.h"
-#include "TArrayF.h"
-#include "TArrayD.h"
-
-class AliESDTrack : public TObject
-{
-public:
-  AliESDTrack();
-  virtual ~AliESDTrack() {}
-  
-protected:
-  Int_t     fTrackID;        // Track number
 
-  // Track parameters at Vertex
-  TArrayD   fPVertex;        // (5) Track parameters
-  TArrayD   fPEVertex;       // (15) Track parameter errors
-  
-  // Track parameters at first measured point
-  TArrayD   fPFMeasPoint;        // (6) Track parameters
-  TArrayD   fPFMeasPointErr;     // (15) Track parameter error
-  
-  // Track parameters at last measured point
-  TArrayD   fPLMeasPoint;        // (6) Track parameters
-  TArrayD   fPLMeasPointErr;     // (15) Track parameter error
-
-  Float_t   fTrackLength;        // Track length
-  Float_t   fTrackLengthErr;     // Track length error
-  Int_t     fStopVertex;         // Index of stop vertex
-  
-  Int_t     fNPointsITS;         // Number of points in ITS
-  Int_t     fNPointsTPC;         // Number of points in TPC
-  Int_t     fNPointsTRD;         // Number of points in TRD
-  Float_t   fMeanResITS;         // Mean residual in ITS
-  Float_t   fMeanResTPC;         // Mean residual in TPC
-  Float_t   fMeanResTRD;         // Mean residual in TRD
-  Float_t   fGlobalChi2;         // Global chi square
-  Int_t     fParticleType;       // PDG code
-
-  Float_t   fPIDprobPi;          // PID probability for pi
-  Float_t   fPIDprobK;           // PID probability for K
-  Float_t   fPIDprobP;           // PID probability for p
-  Float_t   fPIDprobE;           // PID probability for e
-
-private:
-  AliESDTrack(const AliESDTrack &) {}
-  AliESDTrack & operator=(const AliESDTrack &) {return (*this);}
-  
-  ClassDef(AliESDTrack,1)  //ESDTrack 
-};
-
-
-class AliESDVertex : public TObject
-{
-public:
-  AliESDVertex();
-  virtual ~AliESDVertex() {}
-  
-protected:
-  Int_t        fNPrimary;               // Number of primary tracks
-  TArrayF      fCoordinates;            // (3) Vertex coordinates
-  TArrayF      fErrorMatrix;            // (6) Error Matrix
-  TObjArray    fPrimaryTracks;          // List of primary tracks
-  Float_t      fEffectiveMass;          // Effective Mass
-  Float_t      fEffectiveMassError;     // Effective Mass Error
-private:
-  AliESDVertex(const AliESDVertex &) {}
-  AliESDVertex & operator=(const AliESDVertex &) {return (*this);}
-  
-  ClassDef(AliESDVertex,1)  //ESDVertex 
-};
+/* $Id$ */
 
-class AliESD : public TObject
-{
+//-------------------------------------------------------------------------
+//                          Class AliESD
+//   This is the class to deal with during the physical analysis of data
+//      
+//         Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch 
+//-------------------------------------------------------------------------
+
+#include <TClonesArray.h>
+#include <TObject.h>
+
+#include "AliESDMuonTrack.h"
+#include "AliESDPmdTrack.h"
+#include "AliESDTrdTrack.h"
+#include "AliESDVertex.h"
+#include "AliESDcascade.h"
+#include "AliESDkink.h"
+#include "AliESDtrack.h"
+#include "AliESDHLTtrack.h"
+#include "AliESDCaloCluster.h"
+#include "AliESDv0.h"
+#include "AliESDV0MI.h"
+#include "AliESDFMD.h"
+
+class AliESDfriend;
+
+class AliESD : public TObject {
 public:
   AliESD();
-  virtual ~AliESD() {}
+  virtual ~AliESD()
 
-  Int_t EventNumber() const {return fEventNumber;}
-  Int_t RunNumber() const {return fRunNumber;}
-  Long_t Trigger() const {return fTrigger;}
-  
-  Int_t BitDDL() const {return fBitDDL;}
-  Int_t NSecVertex() const {return fNSecVertex;}
-  Float_t NParticipants() const {return fNParticipants;}
+  void SetESDfriend(const AliESDfriend *f);
+  void GetESDfriend(AliESDfriend *f) const;
+
+  void SetEventNumber(Int_t n) {fEventNumber=n;}
+  void SetRunNumber(Int_t n) {fRunNumber=n;}
+  void SetTriggerMask(ULong64_t n) {fTriggerMask=n;}
+  void SetTriggerCluster(UChar_t n) {fTriggerCluster = n;}
+  void SetMagneticField(Float_t mf){fMagneticField = mf;}
+  Float_t GetMagneticField() const {return fMagneticField;}
   
+  AliESDtrack *GetTrack(Int_t i) const {
+    return (AliESDtrack *)fTracks.UncheckedAt(i);
+  }
+  AliESDHLTtrack *GetHLTConfMapTrack(Int_t i) const {
+    return (AliESDHLTtrack *)fHLTConfMapTracks.UncheckedAt(i);
+  }
+  AliESDHLTtrack *GetHLTHoughTrack(Int_t i) const {
+    return (AliESDHLTtrack *)fHLTHoughTracks.UncheckedAt(i);
+  }
+  AliESDMuonTrack *GetMuonTrack(Int_t i) const {
+    return (AliESDMuonTrack *)fMuonTracks.UncheckedAt(i);
+  }
+  AliESDPmdTrack *GetPmdTrack(Int_t i) const {
+    return (AliESDPmdTrack *)fPmdTracks.UncheckedAt(i);
+  }
+  AliESDTrdTrack *GetTrdTrack(Int_t i) const {
+    return (AliESDTrdTrack *)fTrdTracks.UncheckedAt(i);
+  }
+
+  Int_t  AddTrack(const AliESDtrack *t) {
+    AliESDtrack * track = new(fTracks[fTracks.GetEntriesFast()]) AliESDtrack(*t);
+    track->SetID(fTracks.GetEntriesFast()-1);
+    return  track->GetID();
+    
+  }
+  void AddHLTConfMapTrack(const AliESDHLTtrack *t) {
+    new(fHLTConfMapTracks[fHLTConfMapTracks.GetEntriesFast()]) AliESDHLTtrack(*t);
+  }
+  void AddHLTHoughTrack(const AliESDHLTtrack *t) {
+    new(fHLTHoughTracks[fHLTHoughTracks.GetEntriesFast()]) AliESDHLTtrack(*t);
+  }
+  void AddMuonTrack(const AliESDMuonTrack *t) {
+    new(fMuonTracks[fMuonTracks.GetEntriesFast()]) AliESDMuonTrack(*t);
+  }
+  void AddPmdTrack(const AliESDPmdTrack *t) {
+    new(fPmdTracks[fPmdTracks.GetEntriesFast()]) AliESDPmdTrack(*t);
+  }
+  void AddTrdTrack(const AliESDTrdTrack *t) {
+    new(fTrdTracks[fTrdTracks.GetEntriesFast()]) AliESDTrdTrack(*t);
+  }
+
+  AliESDv0 *GetV0(Int_t i) const {
+    return (AliESDv0 *)fV0s.UncheckedAt(i);
+  }
+  void AddV0(const AliESDv0 *v) {
+    new(fV0s[fV0s.GetEntriesFast()]) AliESDv0(*v);
+  }
+  void UpdateV0PIDs();
+
+  AliESDcascade *GetCascade(Int_t i) const {
+    return (AliESDcascade *)fCascades.UncheckedAt(i);
+  }
+  void AddCascade(const AliESDcascade *c) {
+    new(fCascades[fCascades.GetEntriesFast()]) AliESDcascade(*c);
+  }
+
+  AliESDkink *GetKink(Int_t i) const {
+    return (AliESDkink *)fKinks.UncheckedAt(i);
+  }
+  Int_t AddKink(const AliESDkink *c) {
+    AliESDkink * kink = new(fKinks[fKinks.GetEntriesFast()]) AliESDkink(*c);
+    kink->SetID(fKinks.GetEntriesFast());
+    return fKinks.GetEntriesFast()-1;
+  }
+
+  AliESDV0MI *GetV0MI(Int_t i) const {
+    return (AliESDV0MI *)fV0MIs.UncheckedAt(i);
+  }
+  Int_t AddV0MI(const AliESDV0MI *c) {
+    AliESDV0MI * v0 = new(fV0MIs[fV0MIs.GetEntriesFast()]) AliESDV0MI(*c);
+    v0->SetID(fV0MIs.GetEntriesFast()-1);
+    return fV0MIs.GetEntriesFast()-1;
+  }
+
+  AliESDCaloCluster *GetCaloCluster(Int_t i) const {
+    return (AliESDCaloCluster *)fCaloClusters.UncheckedAt(i);
+  }
+  Int_t AddCaloCluster(const AliESDCaloCluster *c) {
+    AliESDCaloCluster *clus = new(fCaloClusters[fCaloClusters.GetEntriesFast()]) AliESDCaloCluster(*c);
+    clus->SetID(fCaloClusters.GetEntriesFast()-1);
+    return fCaloClusters.GetEntriesFast()-1;
+  }
+    
+  void SetVertex(const AliESDVertex* vertex) {
+     fSPDVertex=new AliESDVertex(*vertex);
+  }
+  const AliESDVertex* GetVertex() const {return fSPDVertex;};
+
+  void SetPrimaryVertex(const AliESDVertex* vertex) {
+     fPrimaryVertex=new AliESDVertex(*vertex);
+  }
+  const AliESDVertex* GetPrimaryVertex() const {return fPrimaryVertex;};
+
+  Int_t  GetEventNumber() const {return fEventNumber;}
+  Int_t  GetRunNumber() const {return fRunNumber;}
+  ULong64_t GetTriggerMask() const {return fTriggerMask;}
+  UChar_t  GetTriggerCluster() const {return fTriggerCluster;}
   
+  Int_t GetNumberOfTracks()     const {return fTracks.GetEntriesFast();}
+  Int_t GetNumberOfHLTConfMapTracks()     const {return fHLTConfMapTracks.GetEntriesFast();}
+  Int_t GetNumberOfHLTHoughTracks()     const {return fHLTHoughTracks.GetEntriesFast();}
+  Int_t GetNumberOfMuonTracks() const {return fMuonTracks.GetEntriesFast();}
+  Int_t GetNumberOfPmdTracks() const {return fPmdTracks.GetEntriesFast();}
+  Int_t GetNumberOfTrdTracks() const {return fTrdTracks.GetEntriesFast();}
+  Int_t GetNumberOfV0s()      const {return fV0s.GetEntriesFast();}
+  Int_t GetNumberOfCascades() const {return fCascades.GetEntriesFast();}
+  Int_t GetNumberOfKinks() const {return fKinks.GetEntriesFast();}
+  Int_t GetNumberOfV0MIs() const {return fV0MIs.GetEntriesFast();}
+  Int_t GetNumberOfCaloClusters() const {return fCaloClusters.GetEntriesFast();}
+
+  Int_t GetNumberOfEMCALClusters() const {return fEMCALClusters;}
+  void  SetNumberOfEMCALClusters(Int_t clus) {fEMCALClusters = clus;}
+  Int_t GetFirstEMCALCluster() const {return fFirstEMCALCluster;}
+  void  SetFirstEMCALCluster(Int_t index) {fFirstEMCALCluster = index;}
+
+  Int_t GetNumberOfPHOSClusters() const {return fPHOSClusters;}
+  void  SetNumberOfPHOSClusters(Int_t part) { fPHOSClusters = part ; }
+  void  SetFirstPHOSCluster(Int_t index) { fFirstPHOSCluster = index ; } 
+  Int_t GetFirstPHOSCluster() const  { return fFirstPHOSCluster ; }
+
+  Float_t GetT0zVertex() const {return fT0zVertex;}
+  void SetT0zVertex(Float_t z) {fT0zVertex=z;}
+  Float_t GetT0() const {return fT0timeStart;}
+  void SetT0(Float_t timeStart) {fT0timeStart = timeStart;}
+  const Float_t * GetT0time() const {return fT0time;}
+  void SetT0time(Float_t time[24]) {
+    for (Int_t i=0; i<24; i++) fT0time[i] = time[i];
+  }
+  const Float_t * GetT0amplitude() const {return fT0amplitude;}
+  void SetT0amplitude(Float_t amp[24]) {
+    for (Int_t i=0; i<24; i++) fT0amplitude[i] = amp[i];
+  }
+
+  Float_t GetZDCN1Energy() const {return fZDCN1Energy;}
+  Float_t GetZDCP1Energy() const {return fZDCP1Energy;}
+  Float_t GetZDCN2Energy() const {return fZDCN2Energy;}
+  Float_t GetZDCP2Energy() const {return fZDCP2Energy;}
+  Float_t GetZDCEMEnergy() const {return fZDCEMEnergy;}
+  Int_t   GetZDCParticipants() const {return fZDCParticipants;}
+  void    SetZDC(Float_t n1Energy, Float_t p1Energy, Float_t emEnergy,
+                 Float_t n2Energy, Float_t p2Energy, Int_t participants) 
+   {fZDCN1Energy=n1Energy; fZDCP1Energy=p1Energy; fZDCEMEnergy=emEnergy;
+    fZDCN2Energy=n2Energy; fZDCP2Energy=p2Energy; fZDCParticipants=participants;}
+
+  void ResetV0s() { fV0s.Clear(); }
+  void ResetCascades() { fCascades.Clear(); }
+  void Reset();
+
+  void  Print(Option_t *option="") const;
+
+  void SetFMDData(AliESDFMD * obj) {
+    fESDFMD = new AliESDFMD(*obj);
+  }
+
+  AliESDFMD * GetFMDData(){ return fESDFMD;}
+   
 protected:
 
   // Event Identification
-  Int_t        fEventNumber;            // Event Number
-  Int_t        fRunNumber;              // Run Number
-  Long_t       fTrigger;                // Trigger Type (cfg Transverse Energy&Max trans ch mom)
-  Int_t        fRecoVersion;            // Version of reconstruction 
-
-  // Summary Information
-  Int_t        fBitDDL;                 // Bitmap of active DDL
-  Int_t        fNSecVertex;             // Number of Secondary Vertexes
-  Float_t      fNParticipants;          // Estimated Number of participants
-  Float_t      fNPartError;             // N of participant error
-  Int_t        fNElectron;              // N of electrons
-  Int_t        fNMuons;                 // N of muons
-  Int_t        fNPions;                 // N of pions
-  Int_t        fNKaons;                 // N of kaons
-  Int_t        fNProtons;               // N of protons
-  Int_t        fNPHOSPhotons;           // N of photons in PHOS
-  Int_t        fNPHOSNeutrons;          // N of neutrons in PHOS
-  Int_t        fNPHOSCCluster;          // N of charged clusters in PHOS
-  Int_t        fNEMCALCluster;          // N of clusters in EMCAL
-  Int_t        fNPMDCluster;            // N of clusters in PMD
-  Float_t      fTMaxClusterEnergy;      // Transverse energy of biggest cluster
-  Float_t      fTMaxPCharged;           // Biggest transverse momentum of charged particles
-  TArrayI      fNCharged;               // Charged Multiplicity
-  Float_t      fTotTranEnergy;          // Total transverse energy
-
-  // Primary Vertex Object
-  AliESDVertex fESDVertex;              // Primary Vertex Object
-  TObjArray    fSecVertex;              // List secondary vertexes
-  TObjArray    fNonAssTrack;            // List of non assigned tracks
-  TObjArray    fPhoton;                 // List of photons
-  TObjArray    fNeutron;                // List of neutrons
-  TObjArray    fEMCALCluster;           // List of EMCAL clusters
-  TObjArray    fPMDCluster;             // List of PMD clusters
-
-private:
-  AliESD(const AliESD &) {}
-  AliESD & operator=(const AliESD &) {return (*this);}
-  
-  ClassDef(AliESD,1)  //ESD 
+  Int_t        fEventNumber;     // Event Number
+  Int_t        fRunNumber;       // Run Number
+  ULong64_t    fTriggerMask;     // Trigger Type (mask)
+  UChar_t      fTriggerCluster;  // Trigger cluster (mask)
+  Int_t        fRecoVersion;     // Version of reconstruction 
+  Float_t      fMagneticField;   // Solenoid Magnetic Field in kG : for compatibility with AliMagF
+
+  Float_t      fZDCN1Energy;      // reconstructed energy in the neutron ZDC
+  Float_t      fZDCP1Energy;      // reconstructed energy in the proton ZDC
+  Float_t      fZDCN2Energy;      // reconstructed energy in the neutron ZDC
+  Float_t      fZDCP2Energy;      // reconstructed energy in the proton ZDC
+  Float_t      fZDCEMEnergy;     // reconstructed energy in the electromagnetic ZDC
+  Int_t        fZDCParticipants; // number of participants estimated by the ZDC
+
+  Float_t      fT0zVertex;       // vertex z position estimated by the START
+  AliESDVertex *fSPDVertex;      // Primary vertex estimated by the SPD
+  AliESDVertex *fPrimaryVertex;  // Primary vertex estimated using ESD tracks
+
+  Float_t      fT0timeStart;     // interaction time estimated by the START
+  Float_t      fT0time[24];      // best TOF on each START PMT
+  Float_t      fT0amplitude[24]; // number of particles(MIPs) on each START PMT
+
+  TClonesArray fTracks;          // ESD tracks
+  TClonesArray fHLTConfMapTracks;// HLT ESD tracks from Conformal Mapper method
+  TClonesArray fHLTHoughTracks;  // HLT ESD tracks from Hough Transform method
+  TClonesArray fMuonTracks;      // MUON ESD tracks
+  TClonesArray fPmdTracks;       // PMD ESD tracks
+  TClonesArray fTrdTracks;       // TRD ESD tracks (triggered)
+  TClonesArray fV0s;             // V0 vertices
+  TClonesArray fCascades;        // Cascade vertices
+  TClonesArray fKinks;           // Kinks
+  TClonesArray fV0MIs;           // V0MI
+  TClonesArray fCaloClusters;    // Calorimeter clusters for PHOS/EMCAL
+  Int_t        fEMCALClusters;   // Number of EMCAL clusters (subset of caloclusters)
+  Int_t        fFirstEMCALCluster; // First EMCAL cluster in the fCaloClusters list 
+
+  Int_t        fPHOSClusters;     // Number of PHOS clusters (subset of caloclusters)
+  Int_t        fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list 
+  AliESDFMD *  fESDFMD; // FMD object containing rough multiplicity
+
+  ClassDef(AliESD,11)  //ESD class 
 };
-
 #endif