Bool_t desc1 = (mask1 >> 3) & 0x1;
Bool_t desc2 = (mask2 >> 3) & 0x1;
if (desc1==0 || desc2==0) {
- AliError(Form("TPC not in DAQ/RECO: %u (%u)/%u (%u)",
- mask1, esdevent->GetESDRun()->GetDetectorsInReco(),
- mask2, esdevent->GetESDRun()->GetDetectorsInDAQ()));
+// AliError(Form("TPC not in DAQ/RECO: %u (%u)/%u (%u)",
+// mask1, esdevent->GetESDRun()->GetDetectorsInReco(),
+// mask2, esdevent->GetESDRun()->GetDetectorsInDAQ()));
fITSTrackSA=kTRUE;
}
}
if(!ExtrapolateTrackToEMCalSurface(&emcalParam, 430., fMass, fStepSurface, eta, phi))
{
if(aodevent && trackParam) delete trackParam;
+ if(fITSTrackSA && trackParam) delete trackParam;
continue;
}
if(TMath::Abs(eta)>0.75 || (phi) < 70*TMath::DegToRad() || (phi) > 190*TMath::DegToRad())
{
if(aodevent && trackParam) delete trackParam;
+ if(fITSTrackSA && trackParam) delete trackParam;
continue;
}
matched++;
}
if(aodevent && trackParam) delete trackParam;
+ if(fITSTrackSA && trackParam) delete trackParam;
}//track loop
if(clusterArray)
if(!trackParam) return index;
AliExternalTrackParam emcalParam(*trackParam);
Float_t eta, phi;
- if(!ExtrapolateTrackToEMCalSurface(&emcalParam, 430., fMass, fStepSurface, eta, phi)) return index;
- if(TMath::Abs(eta)>0.75 || (phi) < 70*TMath::DegToRad() || (phi) > 190*TMath::DegToRad()) return index;
+ if(!ExtrapolateTrackToEMCalSurface(&emcalParam, 430., fMass, fStepSurface, eta, phi)) {
+ if(fITSTrackSA) delete trackParam;
+ return index;
+ }
+ if(TMath::Abs(eta)>0.75 || (phi) < 70*TMath::DegToRad() || (phi) > 190*TMath::DegToRad()){
+ if(fITSTrackSA) delete trackParam;
+ return index;
+ }
+
TObjArray *clusterArr = new TObjArray(event->GetNumberOfCaloClusters());
for(Int_t icl=0; icl<event->GetNumberOfCaloClusters(); icl++)
index = FindMatchedClusterInClusterArr(&emcalParam, &emcalParam, clusterArr, dEta, dPhi);
clusterArr->Clear();
delete clusterArr;
-
+ if(fITSTrackSA) delete trackParam;
+
return index;
}
// Otherwise use the TPCInner point
dEta = -999, dPhi = -999;
+ Bool_t ITSTrackSA = 0;
AliExternalTrackParam *trkParam = 0;
trkParam = const_cast<AliExternalTrackParam*>(friendTrack->GetTPCOut());
else if(track->GetInnerParam())
trkParam = const_cast<AliExternalTrackParam*>(track->GetInnerParam());
- else
+ else{
trkParam = new AliExternalTrackParam(*track); //If there is ITSSa track
+ ITSTrackSA = 1;
+ }
if(!trkParam) return kFALSE;
AliExternalTrackParam trkParamTmp (*trkParam);
- if(!AliEMCALRecoUtils::ExtrapolateTrackToCluster(&trkParamTmp, cluster, track->GetMass(kTRUE), GetRecParam()->GetExtrapolateStep(), dEta, dPhi)) return kFALSE;
+ if(!AliEMCALRecoUtils::ExtrapolateTrackToCluster(&trkParamTmp, cluster, track->GetMass(kTRUE), GetRecParam()->GetExtrapolateStep(), dEta, dPhi)){
+ if(ITSTrackSA) delete trkParam;
+ return kFALSE;
+ }
+ if(ITSTrackSA) delete trkParam;
return kTRUE;
}
Bool_t desc1 = (mask1 >> 3) & 0x1;
Bool_t desc2 = (mask2 >> 3) & 0x1;
if (desc1==0 || desc2==0) {
- AliError(Form("TPC not in DAQ/RECO: %u (%u)/%u (%u)",
- mask1, esd->GetESDRun()->GetDetectorsInReco(),
- mask2, esd->GetESDRun()->GetDetectorsInDAQ()));
+// AliError(Form("TPC not in DAQ/RECO: %u (%u)/%u (%u)",
+// mask1, esd->GetESDRun()->GetDetectorsInReco(),
+// mask2, esd->GetESDRun()->GetDetectorsInDAQ()));
fITSTrackSA = kTRUE;
}
AliExternalTrackParam trkParamTmp(*trkParam);
Float_t eta, phi;
- if(!AliEMCALRecoUtils::ExtrapolateTrackToEMCalSurface(&trkParamTmp, 430., track->GetMass(kTRUE), fStep, eta, phi)) return index;
+ if(!AliEMCALRecoUtils::ExtrapolateTrackToEMCalSurface(&trkParamTmp, 430., track->GetMass(kTRUE), fStep, eta, phi)) {
+ if(fITSTrackSA) delete trkParam;
+ return index;
+ }
track->SetTrackPhiEtaOnEMCal(phi,eta);
- if(TMath::Abs(eta)>0.75 || (phi) < 70*TMath::DegToRad() || (phi) > 190*TMath::DegToRad()) return index;
+ if(TMath::Abs(eta)>0.75 || (phi) < 70*TMath::DegToRad() || (phi) > 190*TMath::DegToRad()){
+ if(fITSTrackSA) delete trkParam;
+ return index;
+ }
//Perform extrapolation
Double_t trkPos[3];
index=ic;
}
}
+
+ if(fITSTrackSA) delete trkParam;
return index;
}