//______________________________________________________________________________
AliESD::AliESD():
- fEventNumber(0),
+ fEventNumberInFile(0),
+ fBunchCrossNumber(0),
+ fOrbitNumber(0),
+ fPeriodNumber(0),
fRunNumber(0),
+ fTimeStamp(0),
+ fEventType(0),
fTriggerMask(0),
fTriggerCluster(0),
fRecoVersion(0),
fZDCEMEnergy(0),
fZDCParticipants(0),
fT0zVertex(0),
- fSPDVertex(0),
- fPrimaryVertex(0),
+ fSPDVertex(),
+ fPrimaryVertex(),
+ fSPDMult(),
fT0timeStart(0),
fTracks("AliESDtrack",15000),
fHLTConfMapTracks("AliESDHLTtrack",25000),
fV0s("AliESDv0",200),
fCascades("AliESDcascade",20),
fKinks("AliESDkink",4000),
- fV0MIs("AliESDV0MI",4000),
fCaloClusters("AliESDCaloCluster",10000),
fEMCALClusters(0),
fFirstEMCALCluster(-1),
+ fEMCALTriggerPosition(0x0),
+ fEMCALTriggerAmplitudes(0x0),
fPHOSClusters(0),
fFirstPHOSCluster(-1),
- fESDFMD(0x0)
+ fPHOSTriggerPosition(0x0),
+ fPHOSTriggerAmplitudes(0x0),
+ fESDFMD(0x0),
+ fESDVZERO(0x0),
+ fErrorLogs("AliRawDataErrorLog",5)
+
{
for (Int_t i=0; i<24; i++) {
fT0time[i] = 0;
fT0amplitude[i] = 0;
}
+ for (Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
+ for (Int_t i=0; i<3; i++) fDiamondCovXY[i]=0.;
+}
+//______________________________________________________________________________
+AliESD::AliESD(const AliESD& esd):
+ TObject(esd),
+ fEventNumberInFile(esd.fEventNumberInFile),
+ fBunchCrossNumber(esd.fBunchCrossNumber),
+ fOrbitNumber(esd.fOrbitNumber),
+ fPeriodNumber(esd.fPeriodNumber),
+ fRunNumber(esd.fRunNumber),
+ fTimeStamp(esd.fTimeStamp),
+ fEventType(esd.fEventType),
+ fTriggerMask(esd.fTriggerMask),
+ fTriggerCluster(esd.fTriggerCluster),
+ fRecoVersion(esd.fRecoVersion),
+ fMagneticField(esd.fMagneticField),
+ fZDCN1Energy(esd.fZDCN1Energy),
+ fZDCP1Energy(esd.fZDCP1Energy),
+ fZDCN2Energy(esd.fZDCN2Energy),
+ fZDCP2Energy(esd.fZDCP2Energy),
+ fZDCEMEnergy(esd.fZDCEMEnergy),
+ fZDCParticipants(esd.fZDCParticipants),
+ fT0zVertex(esd.fT0zVertex),
+ fSPDVertex(esd.fSPDVertex),
+ fPrimaryVertex(esd.fPrimaryVertex),
+ fSPDMult(esd.fSPDMult),
+ fT0timeStart(esd.fT0timeStart),
+ fTracks(*((TClonesArray*)esd.fTracks.Clone())),
+ fHLTConfMapTracks(*((TClonesArray*)esd.fHLTConfMapTracks.Clone())),
+ fHLTHoughTracks(*((TClonesArray*)esd.fHLTHoughTracks.Clone())),
+ fMuonTracks(*((TClonesArray*)esd.fMuonTracks.Clone())),
+ fPmdTracks(*((TClonesArray*)esd.fPmdTracks.Clone())),
+ fTrdTracks(*((TClonesArray*)esd.fTrdTracks.Clone())),
+ fV0s(*((TClonesArray*)esd.fV0s.Clone())),
+ fCascades(*((TClonesArray*)esd.fCascades.Clone())),
+ fKinks(*((TClonesArray*)esd.fKinks.Clone())),
+ fCaloClusters(*((TClonesArray*)esd.fCaloClusters.Clone())),
+ fEMCALClusters(esd.fEMCALClusters),
+ fFirstEMCALCluster(esd.fFirstEMCALCluster),
+ fEMCALTriggerPosition(esd. fEMCALTriggerPosition),
+ fEMCALTriggerAmplitudes(esd.fEMCALTriggerAmplitudes),
+ fPHOSClusters(esd.fPHOSClusters),
+ fFirstPHOSCluster(esd.fFirstPHOSCluster),
+ fPHOSTriggerPosition(esd.fPHOSTriggerPosition),
+ fPHOSTriggerAmplitudes(esd.fPHOSTriggerAmplitudes),
+ fESDFMD(esd.fESDFMD),
+ fESDVZERO(esd.fESDVZERO),
+ fErrorLogs(*((TClonesArray*)esd.fErrorLogs.Clone()))
+{
+ for (Int_t i=0; i<24; i++) {
+ fT0time[i] = esd.fT0time[i];
+ fT0amplitude[i] = esd.fT0amplitude[i];
+ }
+ for (Int_t i=0; i<2; i++) fDiamondXY[i]=esd.fDiamondXY[i];
+ for (Int_t i=0; i<3; i++) fDiamondCovXY[i]=esd.fDiamondCovXY[i];
}
+//______________________________________________________________________________
+AliESD & AliESD::operator=(const AliESD& source) {
+
+ // Assignment operator
+
+ if(&source == this) return *this;
+
+ fEventNumberInFile = source.fEventNumberInFile;
+ fBunchCrossNumber = source.fBunchCrossNumber;
+ fOrbitNumber = source.fOrbitNumber;
+ fPeriodNumber = source.fPeriodNumber;
+ fRunNumber = source.fRunNumber;
+ fTimeStamp = source.fTimeStamp;
+ fEventType = source.fEventType;
+ fTriggerMask = source.fTriggerMask;
+ fTriggerCluster = source.fTriggerCluster;
+ fRecoVersion = source.fRecoVersion;
+ fMagneticField = source.fMagneticField;
+ fZDCN1Energy = source.fZDCN1Energy;
+ fZDCP1Energy = source.fZDCP1Energy;
+ fZDCN2Energy = source.fZDCN2Energy;
+ fZDCP2Energy = source.fZDCP2Energy;
+ fZDCEMEnergy = source.fZDCEMEnergy;
+ fZDCParticipants = source.fZDCParticipants;
+ fT0zVertex = source.fT0zVertex;
+ fSPDVertex = source.fSPDVertex;
+ fPrimaryVertex = source.fPrimaryVertex;
+ fSPDMult = source.fSPDMult;
+ fT0timeStart = source.fT0timeStart;
+ fTracks = *((TClonesArray*)source.fTracks.Clone());
+ fHLTConfMapTracks = *((TClonesArray*)source.fHLTConfMapTracks.Clone());
+ fHLTHoughTracks = *((TClonesArray*)source.fHLTHoughTracks.Clone());
+ fMuonTracks = *((TClonesArray*)source.fMuonTracks.Clone());
+ fPmdTracks = *((TClonesArray*)source.fPmdTracks.Clone());
+ fTrdTracks = *((TClonesArray*)source.fTrdTracks.Clone());
+ fV0s = *((TClonesArray*)source.fV0s.Clone());
+ fCascades = *((TClonesArray*)source.fCascades.Clone());
+ fKinks = *((TClonesArray*)source.fKinks.Clone());
+ fCaloClusters = *((TClonesArray*)source.fCaloClusters.Clone());
+ fEMCALClusters = source.fEMCALClusters;
+ fFirstEMCALCluster = source.fFirstEMCALCluster;
+ fPHOSClusters = source.fPHOSClusters;
+ fFirstPHOSCluster = source.fFirstPHOSCluster;
+ fESDFMD = source.fESDFMD;
+ fESDVZERO = source.fESDVZERO;
+ fEMCALTriggerPosition=source. fEMCALTriggerPosition;
+ fEMCALTriggerAmplitudes=source.fEMCALTriggerAmplitudes;
+ fPHOSTriggerPosition=source.fPHOSTriggerPosition;
+ fPHOSTriggerAmplitudes=source.fPHOSTriggerAmplitudes;
+ fErrorLogs = *((TClonesArray*)source.fErrorLogs.Clone());
+
+ for (Int_t i=0; i<24; i++) {
+ fT0time[i] = source.fT0time[i];
+ fT0amplitude[i] = source.fT0amplitude[i];
+ }
+ for (Int_t i=0; i<2; i++) fDiamondXY[i]=source.fDiamondXY[i];
+ for (Int_t i=0; i<3; i++) fDiamondCovXY[i]=source.fDiamondCovXY[i];
+
+ return *this;
+
+}
+
+
//______________________________________________________________________________
AliESD::~AliESD()
{
//
// Standard destructor
//
- delete fSPDVertex;
- delete fPrimaryVertex;
fTracks.Delete();
fHLTConfMapTracks.Delete();
fHLTHoughTracks.Delete();
fV0s.Delete();
fCascades.Delete();
fKinks.Delete();
- fV0MIs.Delete();
fCaloClusters.Delete();
delete fESDFMD;
-}
+ delete fESDVZERO;
+// fEMCALTriggerPosition->Delete();
+// fEMCALTriggerAmplitudes->Delete();
+// fPHOSTriggerPosition->Delete();
+// fPHOSTriggerAmplitudes->Delete();
+// delete fEMCALTriggerPosition;
+// delete fEMCALTriggerAmplitudes;
+// delete fPHOSTriggerPosition;
+// delete fPHOSTriggerAmplitudes;
+ fErrorLogs.Delete();
-void AliESD::UpdateV0PIDs()
-{
- //
- //
- //
- Int_t nV0 = GetNumberOfV0MIs();
- for (Int_t i=0;i<nV0;i++){
- AliESDV0MI * v0 = GetV0MI(i);
- AliESDtrack* tp = GetTrack(v0->GetIndex(0));
- AliESDtrack* tm = GetTrack(v0->GetIndex(1));
- if (!tm || !tp){
- printf("BBBUUUUUUUGGGG\n");
- }
- Double_t pp[5],pm[5];
- tp->GetESDpid(pp);
- tm->GetESDpid(pm);
- v0->UpdatePID(pp,pm);
- }
}
//______________________________________________________________________________
void AliESD::Reset()
{
- fEventNumber=0;
+ fEventNumberInFile=0;
+ fBunchCrossNumber=0;
+ fOrbitNumber=0;
+ fPeriodNumber=0;
fRunNumber=0;
+ fTimeStamp = 0;
+ fEventType = 0;
fTriggerMask=0;
fTriggerCluster=0;
fRecoVersion=0;
fZDCParticipants=0;
fT0zVertex=0;
fT0timeStart = 0;
- delete fSPDVertex; fSPDVertex=0;
- delete fPrimaryVertex; fPrimaryVertex=0;
+ new (&fSPDVertex) AliESDVertex();
+ new (&fPrimaryVertex) AliESDVertex();
+ new (&fSPDMult) AliMultiplicity();
fTracks.Clear();
fHLTConfMapTracks.Clear();
fHLTHoughTracks.Clear();
fPHOSClusters=0;
fFirstPHOSCluster=-1;
if (fESDFMD) fESDFMD->Clear();
+// fEMCALTriggerPosition->Clear();
+// fEMCALTriggerAmplitudes->Clear();
+// fPHOSTriggerPosition->Clear();
+// fPHOSTriggerAmplitudes->Clear();
+ fErrorLogs.Clear();
}
+Int_t AliESD::AddV0(const AliESDv0 *v) {
+ //
+ // Add V0
+ //
+ Int_t idx=fV0s.GetEntriesFast();
+ new(fV0s[idx]) AliESDv0(*v);
+ return idx;
+}
+
//______________________________________________________________________________
void AliESD::Print(Option_t *) const
{
// Print header information of the event
//
printf("ESD run information\n");
- printf("Event # %d Run # %d Trigger %lld Magnetic field %f \n",
- GetEventNumber(),
+ printf("Event # in file %d Bunch crossing # %d Orbit # %d Period # %d Run # %d Trigger %lld Magnetic field %f \n",
+ GetEventNumberInFile(),
+ GetBunchCrossNumber(),
+ GetOrbitNumber(),
+ GetPeriodNumber(),
GetRunNumber(),
GetTriggerMask(),
GetMagneticField() );
- if (fPrimaryVertex)
printf("Vertex: (%.4f +- %.4f, %.4f +- %.4f, %.4f +- %.4f) cm\n",
- fPrimaryVertex->GetXv(), fPrimaryVertex->GetXRes(),
- fPrimaryVertex->GetYv(), fPrimaryVertex->GetYRes(),
- fPrimaryVertex->GetZv(), fPrimaryVertex->GetZRes());
+ fPrimaryVertex.GetXv(), fPrimaryVertex.GetXRes(),
+ fPrimaryVertex.GetYv(), fPrimaryVertex.GetYRes(),
+ fPrimaryVertex.GetZv(), fPrimaryVertex.GetZRes());
+ printf("Mean vertex in RUN: X=%.4f Y=%.4f cm\n",
+ GetDiamondX(),GetDiamondY());
+ printf("SPD Multiplicity. Number of tracklets %d \n",
+ fSPDMult.GetNumberOfTracklets());
printf("Event from reconstruction version %d \n",fRecoVersion);
printf("Number of tracks: \n");
printf(" charged %d\n", GetNumberOfTracks());
printf(" v0 %d\n", GetNumberOfV0s());
printf(" cascades %d\n", GetNumberOfCascades());
printf(" kinks %d\n", GetNumberOfKinks());
- printf(" V0MIs %d\n", GetNumberOfV0MIs());
printf(" CaloClusters %d\n", GetNumberOfCaloClusters());
printf(" phos %d\n", GetNumberOfPHOSClusters());
printf(" emcal %d\n", GetNumberOfEMCALClusters());
printf(" FMD %s\n", (fESDFMD ? "yes" : "no"));
+ printf(" VZERO %s\n", (fESDVZERO ? "yes" : "no"));
}
void AliESD::SetESDfriend(const AliESDfriend *ev) {