X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=STEER%2FAliESD.h;h=c92aa446e98077b7bfd4d6d5bf5f2682bc0fe051;hb=80e93fef1261a29fb11dc3ab89698110ec74dbc9;hp=0498a26e5088be03f1f7cf51a73816074ef27c62;hpb=5d8718b8c08f9365868aedf62eb815945720f4da;p=u%2Fmrichter%2FAliRoot.git diff --git a/STEER/AliESD.h b/STEER/AliESD.h index 0498a26e508..c92aa446e98 100644 --- a/STEER/AliESD.h +++ b/STEER/AliESD.h @@ -1,147 +1,307 @@ +// -*- 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 // -// // -/////////////////////////////////////////////////////////////////////////////// +//------------------------------------------------------------------------- +// Class AliESD +// This is the class to deal with during the physical analysis of data +// +// Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch +//------------------------------------------------------------------------- + +#include +#include +#include + +#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 "AliESDFMD.h" +#include "AliESDVZERO.h" +#include "AliMultiplicity.h" +#include "AliRawDataErrorLog.h" +#include "AliESDACORDE.h" -#include "TObject.h" -#include "TObjArray.h" -#include "TArrayI.h" -#include "TArrayF.h" -#include "TArrayD.h" +class AliESDfriend; -class AliESDTrack : public TObject -{ +class AliESD : public TObject { public: - AliESDTrack(); - virtual ~AliESDTrack() {} - -protected: - Int_t fTrackID; // Track number + AliESD(); + AliESD(const AliESD&); + virtual ~AliESD(); - // 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 + void SetESDfriend(const AliESDfriend *f); + void GetESDfriend(AliESDfriend *f) const; - 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 & esdt); - AliESDTrack & operator=(const AliESDTrack &) {return (*this);} + void SetEventNumberInFile(Int_t n) {fEventNumberInFile=n;} + void SetBunchCrossNumber(UShort_t n) {fBunchCrossNumber=n;} + void SetOrbitNumber(UInt_t n) {fOrbitNumber=n;} + void SetPeriodNumber(UInt_t n) {fPeriodNumber=n;} + void SetRunNumber(Int_t n) {fRunNumber=n;} + void SetTimeStamp(UInt_t timeStamp){fTimeStamp = timeStamp;} + void SetEventType(UInt_t eventType){fEventType = eventType;} + 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;} - ClassDef(AliESDTrack,1) //ESDTrack -}; + 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); + } + Bool_t Clean(Float_t *cleanPars); + Bool_t RemoveKink(Int_t i); + Bool_t RemoveV0(Int_t i); + Bool_t RemoveTrack(Int_t i); -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 & esdv); - AliESDVertex & operator=(const AliESDVertex &) {return (*this);} - - ClassDef(AliESDVertex,1) //ESDVertex -}; + 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); + } -class AliESD : public TObject -{ -public: - AliESD(); - virtual ~AliESD() {} + AliESDv0 *GetV0(Int_t i) const { + return (AliESDv0 *)fV0s.UncheckedAt(i); + } + Int_t AddV0(const AliESDv0 *v); - 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;} - + 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; + } + + 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 AddPHOSTriggerPosition(TArrayF array) { fPHOSTriggerPosition = new TArrayF(array) ; } + void AddPHOSTriggerAmplitudes(TArrayF array) { fPHOSTriggerAmplitudes = new TArrayF(array) ; } + void AddEMCALTriggerPosition(TArrayF array) { fEMCALTriggerPosition = new TArrayF(array) ; } + void AddEMCALTriggerAmplitudes(TArrayF array){ fEMCALTriggerAmplitudes = new TArrayF(array) ; } + + void SetVertex(const AliESDVertex *vertex) { + new (&fSPDVertex) AliESDVertex(*vertex); + } + const AliESDVertex *GetVertex() const {return &fSPDVertex;} + + void SetMultiplicity(const AliMultiplicity *mul) { + new (&fSPDMult) AliMultiplicity(*mul); + } + const AliMultiplicity *GetMultiplicity() const {return &fSPDMult;} + + void SetPrimaryVertex(const AliESDVertex *vertex) { + new (&fPrimaryVertex) AliESDVertex(*vertex); + } + const AliESDVertex *GetPrimaryVertex() const {return &fPrimaryVertex;} + + void SetDiamond(const AliESDVertex *vertex); + Float_t GetDiamondX() const {return fDiamondXY[0];} + Float_t GetDiamondY() const {return fDiamondXY[1];} + Float_t GetSigma2DiamondX() const {return fDiamondCovXY[0];} + Float_t GetSigma2DiamondY() const {return fDiamondCovXY[2];} + void GetDiamondCovXY(Float_t cov[3]) const { + for(Int_t i=0;i<3;i++) cov[i]=fDiamondCovXY[i]; return; + } + + Int_t GetEventNumberInFile() const {return fEventNumberInFile;} + UShort_t GetBunchCrossNumber() const {return fBunchCrossNumber;} + UInt_t GetOrbitNumber() const {return fOrbitNumber;} + UInt_t GetPeriodNumber() const {return fPeriodNumber;} + Int_t GetRunNumber() const {return fRunNumber;} + UInt_t GetTimeStamp() const { return fTimeStamp;} + UInt_t GetEventType() const { return fEventType;} + 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 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;} + TArrayF *GetEMCALTriggerPosition() const {return fEMCALTriggerPosition;} + TArrayF *GetEMCALTriggerAmplitudes() const {return fEMCALTriggerAmplitudes;} + + 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 ; } + TArrayF *GetPHOSTriggerPosition() const {return fPHOSTriggerPosition;} + TArrayF *GetPHOSTriggerAmplitudes() const {return fPHOSTriggerAmplitudes;} + + + 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; } + + void SetVZEROData(AliESDVZERO * obj) { fESDVZERO = new AliESDVZERO(*obj); } + AliESDVZERO *GetVZEROData(){ return fESDVZERO; } + void SetACORDEData(AliESDACORDE * obj){ fESDACORDE = new AliESDACORDE(*obj); } + AliESDACORDE *GetACORDEDAta(){ return fESDACORDE; } + AliRawDataErrorLog *GetErrorLog(Int_t i) const { + return (AliRawDataErrorLog *)fErrorLogs.UncheckedAt(i); + } + void AddRawDataErrorLog(const AliRawDataErrorLog *log) { + new(fErrorLogs[fErrorLogs.GetEntriesFast()]) AliRawDataErrorLog(*log); + } + Int_t GetNumberOfErrorLogs() const {return fErrorLogs.GetEntriesFast();} + protected: + AliESD &operator=(const AliESD& source); // 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 &esd); - AliESD & operator=(const AliESD &) {return (*this);} - - ClassDef(AliESD,1) //ESD -}; + Int_t fEventNumberInFile;// running Event count in the file + UShort_t fBunchCrossNumber;// Bunch Crossing Number + UInt_t fOrbitNumber; // Orbit Number + UInt_t fPeriodNumber; // Period Number + Int_t fRunNumber; // Run Number + UInt_t fTimeStamp; // Time stamp + UInt_t fEventType; // Type of Event + 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 T0 + AliESDVertex fSPDVertex; // Primary vertex estimated by the SPD + AliESDVertex fPrimaryVertex; // Primary vertex estimated using ESD tracks + Float_t fDiamondXY[2]; // Interaction diamond (x,y) in RUN + Float_t fDiamondCovXY[3]; // Interaction diamond covariance (x,y) in RUN + AliMultiplicity fSPDMult; // SPD tracklet multiplicity + + Float_t fT0timeStart; // interaction time estimated by the T0 + Float_t fT0time[24]; // best TOF on each T0 PMT + Float_t fT0amplitude[24]; // number of particles(MIPs) on each T0 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 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 + TArrayF *fEMCALTriggerPosition; ///(x,y,z of 2x2 and x,y,z of nxn) not position of centroid but of patch corner + TArrayF *fEMCALTriggerAmplitudes; //(2x2 max ampl, 2x2 amp out of patch, nxn max ampl, nxn amp out of patch) + + Int_t fPHOSClusters; // Number of PHOS clusters (subset of caloclusters) + Int_t fFirstPHOSCluster; // First PHOS cluster in the fCaloClusters list + TArrayF *fPHOSTriggerPosition; //(x,y,z of 2x2 and x,y,z of nxn), not position of centroid but of patch corner + TArrayF *fPHOSTriggerAmplitudes; //(2x2 max ampl, 2x2 amp out of patch, nxn max ampl, nxn amp out of patch) + + AliESDFMD *fESDFMD; // FMD object containing rough multiplicity + AliESDVZERO *fESDVZERO; // VZERO object containing rough multiplicity + AliESDACORDE *fESDACORDE; // ACORDE ESD object containing bit pattern + + TClonesArray fErrorLogs; // Raw-data reading error messages + + ClassDef(AliESD,22) //ESD class +}; #endif