primary track and eta cuts for TRD influence on pt resolution @ vertex
authorabercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 12 Mar 2010 10:46:45 +0000 (10:46 +0000)
committerabercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 12 Mar 2010 10:46:45 +0000 (10:46 +0000)
PWG1/TRD/AliTRDcheckESD.cxx

index 49b2af4..29ade2a 100644 (file)
@@ -332,10 +332,13 @@ void AliTRDcheckESD::UserExec(Option_t *){
       if(kBarrel && (status & AliESDtrack::kTRDpid)) h->Fill(ptTRD, kTRDpid);
       if(kBarrel && (status & AliESDtrack::kTRDrefit)) h->Fill(ptTRD, kTRDref);
     }
-    if(HasMC() && kBarrel && (status & AliESDtrack::kTRDout)) {
+    if(HasMC() && 
+      kBarrel && kPhysPrim &&
+      TMath::Abs(eta0) < 0.9 &&
+      (status & AliESDtrack::kTRDrefit)) {
       TH3 *h3 = (TH3S*)fHistos->At(kPtRes);
       Int_t sgn = mcParticle->Charge()>0?1:-1;
-      h3->Fill(pt0, 1.e2*pt/pt0-1.e2, sgn*Pdg2Idx(TMath::Abs(mcParticle->PdgCode())));
+      h3->Fill(pt0, 1.e2*(pt/pt0-1.), sgn*Pdg2Idx(TMath::Abs(mcParticle->PdgCode())));
     }
     if(ip){
       h = (TH2I*)fHistos->At(kTRDmom);