- Int_t nTrack = fESD->GetNumberOfTracks();
- for(Int_t i=0; i<nTrack; ++i){
- AliESDtrack* esdtrack = fESD->GetTrack(i); // pointer to reconstructed to track
- if(!esdtrack) {
- AliError(Form("ERROR: Could not retrieve esdtrack %d",i));
+ Int_t nCluster = fESD->GetNumberOfCaloClusters();
+ for(Int_t i=0; i<nCluster; ++i){
+ AliESDCaloCluster *c1 = fESD->GetCaloCluster(i);
+ if(!c1->IsEMCAL()) continue;
+ if(!IsGoodCluster(c1)) continue;
+ for(Int_t j=i+1; j<nCluster; ++j){
+ AliESDCaloCluster *c2 = fESD->GetCaloCluster(j);
+ if(!c2->IsEMCAL()) continue;
+ if(!IsGoodCluster(c2)) continue;
+ TLorentzVector p1;
+ GetMom(p1, c1, vertex);
+ TLorentzVector p2;
+ GetMom(p2, c2, vertex);
+ FillPion(p1, p2, fEPV0r, fEPV0A, fEPV0C, fEPTPC);
+ }
+ }
+
+ //for track analysis.
+ fTrackCuts->SetAcceptKinkDaughters(kFALSE);
+ fTrackCuts->SetRequireTPCRefit(kTRUE);
+ fTrackCuts->SetRequireITSRefit(kTRUE);
+ fTrackCuts->SetEtaRange(-0.7,0.7);
+ fTrackCuts->SetRequireSigmaToVertex(kTRUE);
+ fTrackCuts->SetMaxChi2PerClusterTPC(3.5);
+ fTrackCuts->SetMinNClustersTPC(100);
+
+ Int_t nTrack = fESD->GetNumberOfTracks();
+ for(Int_t i=0; i<nTrack; ++i){
+ AliESDtrack* esdtrack = fESD->GetTrack(i); // pointer to reconstructed to track
+ if(!fTrackCuts->AcceptTrack(esdtrack))
+ continue;
+ if(!esdtrack) {
+ AliError(Form("ERROR: Could not retrieve esdtrack %d",i));