Adding some track quality selections
authorbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 18 Oct 2007 07:23:52 +0000 (07:23 +0000)
committerbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 18 Oct 2007 07:23:52 +0000 (07:23 +0000)
STEER/AliESDComparison.C

index ddf97b914425755041722a82603d7a18a68ff9f5..522ff9b45cb56928220b3a31008a55d31f6f2942 100644 (file)
@@ -149,15 +149,19 @@ Int_t AliESDComparison(const Char_t *dir=".") {
       AliStack *stack = rl->Stack();
 
       while (ntrk--) {
-         AliESDtrack *t=event->GetTrack(ntrk);
-
-         UInt_t status=AliESDtrack::kESDpid;
-         status|=AliESDtrack::kITSpid; 
-         status|=AliESDtrack::kTPCpid; 
-         status|=AliESDtrack::kTRDpid; 
-         status|=AliESDtrack::kTOFpid; 
-
-        if ((t->GetStatus()&status) == status) {
+        AliESDtrack *t=event->GetTrack(ntrk);
+
+       //*** Some track quality cuts ****
+        if (t->GetITSclusters(0) < 6 ) continue;
+        if (t->GetTPCclusters(0) < 60) continue;
+        //if (t->GetTRDclusters(0) < 60) continue;
+
+        if (!t->IsOn(AliESDtrack::kESDpid)) continue;
+        if (!t->IsOn(AliESDtrack::kITSpid)) continue;
+        if (!t->IsOn(AliESDtrack::kTPCpid)) continue;
+        //if (!t->IsOn(AliESDtrack::kTRDpid)) continue;
+        if (!t->IsOn(AliESDtrack::kTOFpid)) continue;
+        {
            nsel++;
 
            Double_t p=t->GetP();
@@ -169,8 +173,10 @@ Int_t AliESDComparison(const Char_t *dir=".") {
            Int_t code=part->GetPdgCode();
 
            Double_t r[10]; t->GetESDpid(r);
-           //t->GetTRDpid(r);
+           //t->GetITSpid(r);
            //t->GetTPCpid(r);
+           //t->GetTRDpid(r);
+           //t->GetTOFpid(r);
 
            AliPID pid(r);