]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
New conditions on the kITSrefit and kTPCrefit status bits
authorbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 11 Sep 2007 13:00:39 +0000 (13:00 +0000)
committerbelikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 11 Sep 2007 13:00:39 +0000 (13:00 +0000)
STEER/AliV0vertexer.cxx
STEER/AliV0vertexer.h

index 758b582d46a68e3660f195cc4ea9673a2c80c0bf..96a63fa745ea561f593c2ef85c01921f91c68209 100644 (file)
@@ -38,7 +38,7 @@ Double_t AliV0vertexer::fgChi2max=33.; //max chi2
 Double_t AliV0vertexer::fgDNmin=0.05;  //min imp parameter for the 1st daughter
 Double_t AliV0vertexer::fgDPmin=0.05;  //min imp parameter for the 2nd daughter
 Double_t AliV0vertexer::fgDCAmax=0.5;  //max DCA between the daughter tracks
-Double_t AliV0vertexer::fgCPAmax=0.99; //max cosine of V0's pointing angle
+Double_t AliV0vertexer::fgCPAmin=0.99; //min cosine of V0's pointing angle
 Double_t AliV0vertexer::fgRmin=0.2;    //min radius of the fiducial volume
 Double_t AliV0vertexer::fgRmax=100.;   //max radius of the fiducial volume
 
@@ -62,12 +62,10 @@ Int_t AliV0vertexer::Tracks2V0vertices(AliESDEvent *event) {
    Int_t i;
    for (i=0; i<nentr; i++) {
      AliESDtrack *esdTrack=event->GetTrack(i);
-     UInt_t status=esdTrack->GetStatus();
-     UInt_t flags=AliESDtrack::kITSin|AliESDtrack::kTPCin|
-                  AliESDtrack::kTPCpid|AliESDtrack::kESDpid;
+     ULong_t status=esdTrack->GetStatus();
 
      if ((status&AliESDtrack::kITSrefit)==0)
-        if (flags!=status) continue;
+        if ((status&AliESDtrack::kTPCrefit)==0) continue;
 
      Double_t d=esdTrack->GetD(xv,yv,b);
      if (TMath::Abs(d)<fDPmin) continue;
@@ -117,7 +115,7 @@ Int_t AliV0vertexer::Tracks2V0vertices(AliESDEvent *event) {
          if (vertex.GetChi2V0() > fChi2max) continue;
         
         Float_t cpa=vertex.GetV0CosineOfPointingAngle(xv,yv,zv);
-        if (cpa < fCPAmax) continue;
+        if (cpa < fCPAmin) continue;
         vertex.SetDcaV0Daughters(dca);
          vertex.SetV0CosineOfPointingAngle(cpa);
          vertex.ChangeMassHypothesis(kK0Short);
index 2366b7e4ba0cf7b68d29d62b493471598456004f..0945e56c3aa654eb05f41cd93bcc3b8bccfa8084 100644 (file)
@@ -36,7 +36,7 @@ private:
   static
   Double_t fgDCAmax;       // maximal allowed DCA between the daughter tracks 
   static
-  Double_t fgCPAmax;       // maximal allowed cosine of V0's pointing angle
+  Double_t fgCPAmin;       // minimal allowed cosine of V0's pointing angle
   static
   Double_t fgRmin, fgRmax; // max & min radii of the fiducial volume
   
@@ -44,10 +44,10 @@ private:
   Double_t fDNmin;        // min allowed impact parameter for the 1st daughter
   Double_t fDPmin;        // min allowed impact parameter for the 2nd daughter
   Double_t fDCAmax;       // maximal allowed DCA between the daughter tracks 
-  Double_t fCPAmax;       // maximal allowed cosine of V0's pointing angle
+  Double_t fCPAmin;       // minimal allowed cosine of V0's pointing angle
   Double_t fRmin, fRmax;  // max & min radii of the fiducial volume
   
-  ClassDef(AliV0vertexer,2)  // V0 verterxer 
+  ClassDef(AliV0vertexer,3)  // V0 verterxer 
 };
 
 inline AliV0vertexer::AliV0vertexer() :
@@ -56,7 +56,7 @@ inline AliV0vertexer::AliV0vertexer() :
   fDNmin(fgDNmin),
   fDPmin(fgDPmin),
   fDCAmax(fgDCAmax),
-  fCPAmax(fgCPAmax), 
+  fCPAmin(fgCPAmin), 
   fRmin(fgRmin),
   fRmax(fgRmax) 
 {
@@ -65,28 +65,28 @@ inline AliV0vertexer::AliV0vertexer() :
 inline void AliV0vertexer::SetCuts(const Double_t cuts[7]) {
   fChi2max=cuts[0]; 
   fDNmin=cuts[1];   fDPmin=cuts[2];
-  fDCAmax=cuts[3];  fCPAmax=cuts[4];
+  fDCAmax=cuts[3];  fCPAmin=cuts[4];
   fRmin=cuts[5];    fRmax=cuts[6]; 
 }
 
 inline void AliV0vertexer::SetDefaultCuts(const Double_t cuts[7]) {
   fgChi2max=cuts[0]; 
   fgDNmin=cuts[1];   fgDPmin=cuts[2];
-  fgDCAmax=cuts[3];  fgCPAmax=cuts[4];
+  fgDCAmax=cuts[3];  fgCPAmin=cuts[4];
   fgRmin=cuts[5];    fgRmax=cuts[6]; 
 }
 
 inline void AliV0vertexer::GetCuts(Double_t cuts[7]) const {
   cuts[0]=fChi2max; 
   cuts[1]=fDNmin;   cuts[2]=fDPmin; 
-  cuts[3]=fDCAmax;  cuts[4]=fCPAmax;
+  cuts[3]=fDCAmax;  cuts[4]=fCPAmin;
   cuts[5]=fRmin;    cuts[6]=fRmax; 
 }
 
 inline void AliV0vertexer::GetDefaultCuts(Double_t cuts[7]) {
   cuts[0]=fgChi2max; 
   cuts[1]=fgDNmin;   cuts[2]=fgDPmin; 
-  cuts[3]=fgDCAmax;  cuts[4]=fgCPAmax;
+  cuts[3]=fgDCAmax;  cuts[4]=fgCPAmin;
   cuts[5]=fgRmin;    cuts[6]=fgRmax; 
 }