- Bool_t isEmc = kFALSE;
- Int_t type = -1;
- if (esdMode) {
- AliESDtrack *esdtrack = static_cast<AliESDtrack*>(track);
- if (fESDtrackCuts && !fESDtrackCuts->AcceptTrack(esdtrack))
- continue;
- type = 0;
- if (esdtrack->TestBit(BIT(27))) {
- if (esdtrack->TestBit(BIT(28)))
- type = 2;
- else
- type =1;
- }
- if (!fIncludeNoITS && (type==2))
- continue;
- if (TMath::Abs(esdtrack->GetTrackEtaOnEMCal()) < 0.75 &&
- esdtrack->GetTrackPhiOnEMCal() > 70 * TMath::DegToRad() &&
- esdtrack->GetTrackPhiOnEMCal() < 190 * TMath::DegToRad())
- isEmc = kTRUE;
- } else {
- AliAODTrack *aodtrack = static_cast<AliAODTrack*>(track);
- if (fAODfilterBits[0] < 0) {
- if (aodtrack->IsHybridGlobalConstrainedGlobal())
- type = 3;
- else /*not a good track*/
- continue;
- } else {
- if (aodtrack->TestFilterBit(fAODfilterBits[0])) {
- type = 0;
- } else if (aodtrack->TestFilterBit(fAODfilterBits[1])) {
- if ((aodtrack->GetStatus()&AliESDtrack::kITSrefit)==0) {
- if (fIncludeNoITS)
- type = 2;
- else
- continue;
- } else {
- type = 1;
- }
- }
- else {/*not a good track*/
- continue;
- }
- }
- if (fCutMaxFractionSharedTPCClusters > 0) {
- Double_t frac = Double_t(aodtrack->GetTPCnclsS()) / Double_t(aodtrack->GetTPCncls());
- if (frac > fCutMaxFractionSharedTPCClusters)
- continue;
- }
- if (TMath::Abs(track->GetTrackEtaOnEMCal()) < 0.75 &&
- track->GetTrackPhiOnEMCal() > 70 * TMath::DegToRad() &&
- track->GetTrackPhiOnEMCal() < 190 * TMath::DegToRad())
- isEmc = kTRUE;
- }
-