fSPDVertex(),
fPrimaryVertex(),
fSPDMult(),
+ fT0clock(0),
fT0timeStart(0),
+ fT0trig(0),
fTracks("AliESDtrack",15000),
fHLTConfMapTracks("AliESDHLTtrack",25000),
fHLTHoughTracks("AliESDHLTtrack",15000),
fPHOSTriggerAmplitudes(0x0),
fESDFMD(0x0),
fESDVZERO(0x0),
+ fESDACORDE(0x0),
fErrorLogs("AliRawDataErrorLog",5)
{
//
// Standar constructor
//
- for (Int_t i=0; i<24; i++) {
+ for (Int_t i=0; i<3; i++) fT0TOF[i] = 0;
+ for (Int_t i=0; i<24; i++) {
fT0time[i] = 0;
fT0amplitude[i] = 0;
}
fSPDVertex(esd.fSPDVertex),
fPrimaryVertex(esd.fPrimaryVertex),
fSPDMult(esd.fSPDMult),
+ fT0clock(esd.fT0clock),
fT0timeStart(esd.fT0timeStart),
+ fT0trig(esd.fT0trig),
fTracks(*((TClonesArray*)esd.fTracks.Clone())),
fHLTConfMapTracks(*((TClonesArray*)esd.fHLTConfMapTracks.Clone())),
fHLTHoughTracks(*((TClonesArray*)esd.fHLTHoughTracks.Clone())),
fPHOSTriggerAmplitudes(esd.fPHOSTriggerAmplitudes),
fESDFMD(esd.fESDFMD),
fESDVZERO(esd.fESDVZERO),
+ fESDACORDE(esd.fESDACORDE),
fErrorLogs(*((TClonesArray*)esd.fErrorLogs.Clone()))
{
//
// copy constructor
//
+ for (Int_t i=0; i<3; i++)fT0TOF[i] = esd.fT0TOF[i];
for (Int_t i=0; i<24; i++) {
fT0time[i] = esd.fT0time[i];
fT0amplitude[i] = esd.fT0amplitude[i];
delete fEMCALTriggerAmplitudes;
delete fPHOSTriggerPosition;
delete fPHOSTriggerAmplitudes;
+ delete fESDACORDE;
fErrorLogs.Delete();
fZDCEMEnergy=0;
fZDCParticipants=0;
fT0zVertex=0;
- fSPDVertex.~AliESDVertex();
- new (&fSPDVertex) AliESDVertex();
- fPrimaryVertex.~AliESDVertex();
- new (&fPrimaryVertex) AliESDVertex();
+
for (Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
for (Int_t i=0; i<3; i++) fDiamondCovXY[i]=0.;
- fSPDMult.~AliMultiplicity();
- new (&fSPDMult) AliMultiplicity();
+
for (Int_t i=0; i<24; i++) {
fT0time[i] = 0;
fT0amplitude[i] = 0;
}
fT0timeStart = 0;
+ fT0clock = 0;
+ for (Int_t i=0; i<3; i++) fT0TOF[i] = 0;
+//
+ fSPDMult.~AliMultiplicity();
+ new (&fSPDMult) AliMultiplicity();
+ fSPDVertex.~AliESDVertex();
+ new (&fSPDVertex) AliESDVertex();
+ fPrimaryVertex.~AliESDVertex();
+ new (&fPrimaryVertex) AliESDVertex();
+//
fTracks.Delete();
fHLTConfMapTracks.Delete();
fHLTHoughTracks.Delete();
fTrdTracks.Delete();
fV0s.Delete();
fCascades.Delete();
+ fKinks.Delete();
fCaloClusters.Delete();
+//
fEMCALClusters=0;
fFirstEMCALCluster=-1;
fPHOSClusters=0;
fFirstPHOSCluster=-1;
+//
+ if (fEMCALTriggerPosition) fEMCALTriggerPosition ->Reset();
+ if (fEMCALTriggerAmplitudes) fEMCALTriggerAmplitudes->Reset();
+ if (fPHOSTriggerPosition) fPHOSTriggerPosition ->Reset();
+ if (fPHOSTriggerAmplitudes) fPHOSTriggerAmplitudes ->Reset();
+//
if (fESDFMD) fESDFMD->Clear();
+//
if (fESDVZERO){
- fESDVZERO->~AliESDVZERO();
- new (fESDVZERO) AliESDVZERO();
+ fESDVZERO->~AliESDVZERO();
+ new (fESDVZERO) AliESDVZERO();
}
-// fEMCALTriggerPosition->Clear();
-// fEMCALTriggerAmplitudes->Clear();
-// fPHOSTriggerPosition->Clear();
-// fPHOSTriggerAmplitudes->Clear();
+//
+ if (fESDACORDE){
+ fESDACORDE->~AliESDACORDE();
+ new (fESDACORDE) AliESDACORDE();
+ }
+//
fErrorLogs.Delete();
}
Int_t used=0;
+ // Check if this track comes from the reconstructed primary vertex
+ if (fPrimaryVertex.GetStatus()) {
+ UShort_t *primIdx=fPrimaryVertex.GetIndices();
+ Int_t n=fPrimaryVertex.GetNIndices();
+ while (n--) {
+ Int_t idx=Int_t(primIdx[n]);
+ if (rm==idx) return kFALSE;
+ if (idx==last) used++;
+ }
+ }
+
// Check if this track comes from a reconstructed decay
Int_t nv0=GetNumberOfV0s();
for (Int_t n=0; n<nv0; n++) {
new (a[rm]) AliESDtrack(*t);
delete a.RemoveAt(last);
+
if (!used) return kTRUE;
+ // Remap the indices of the tracks used for the primary vertex reconstruction
+ if (fPrimaryVertex.GetStatus()) {
+ UShort_t *primIdx=fPrimaryVertex.GetIndices();
+ Int_t n=fPrimaryVertex.GetNIndices();
+ while (n--) {
+ Int_t idx=Int_t(primIdx[n]);
+ if (idx==last) {
+ primIdx[n]=Short_t(rm);
+ used--;
+ if (!used) return kTRUE;
+ }
+ }
+ }
+
// Remap the indices of the daughters of reconstructed decays
for (Int_t n=0; n<nv0; n++) {
AliESDv0 *v0=GetV0(n);