bug is fixed

diff --git a/PMD/AliPMDClusteringV1.cxx b/PMD/AliPMDClusteringV1.cxx
index e5eeb128774d8e377b0ee7733e3d0ae48ab1c73a..b2e0a076588e24e241754ed0d79e7b7f70653a34 100644 (file)
--- a/PMD/AliPMDClusteringV1.cxx
+++ b/PMD/AliPMDClusteringV1.cxx
@@ -246,7 +246,7 @@ void AliPMDClusteringV1::DoClust(Int_t idet, Int_t ismn,
       clusdata[3]     = cluCELLS;
       clusdata[4]     = cluSIGX;
       clusdata[5]     = cluSIGY;
-      
+
       //
       // Cells associated with a cluster
       //
@@ -292,6 +292,7 @@ void AliPMDClusteringV1::DoClust(Int_t idet, Int_t ismn,
                  celldataTr[ihit]  = celltrack[irow][icol];
                  celldataPid[ihit] = cellpid[irow][icol];
                  celldataAdc[ihit] = (Float_t) celladc[irow][icol];
+
                }
              else
                {
@@ -302,7 +303,7 @@ void AliPMDClusteringV1::DoClust(Int_t idet, Int_t ismn,
            }
          
        }
-
+      
       pmdcl = new AliPMDcluster(idet, ismn, clusdata, celldataX, celldataY,
                                celldataTr, celldataPid, celldataAdc);
       pmdcont->Add(pmdcl);
@@ -527,6 +528,7 @@ void AliPMDClusteringV1::RefClust(Int_t incr, Double_t edepcell[])
            {
              clxy[icltr] = -1;
            }
+
          pmdcludata  = new AliPMDcludata(clusdata,clxy);
          fPMDclucont->Add(pmdcludata);
        }
@@ -847,10 +849,19 @@ void AliPMDClusteringV1::RefClust(Int_t incr, Double_t edepcell[])
                if(totaladc[kcl]>0.)xclust[kcl] = (ax[kcl])/ totaladc[kcl];
                if(totaladc[kcl]>0.)yclust[kcl] = (ay[kcl])/ totaladc[kcl];
                
-               if(totaladc[kcl]>0.)sigxclust[kcl] = (totaladc[kcl]/(pow(totaladc[kcl],2)-totaladc2[kcl]))*ax2[kcl];
-               if(totaladc[kcl]>0.)sigyclust[kcl] = (totaladc[kcl]/(pow(totaladc[kcl],2)-totaladc2[kcl]))*ay2[kcl];
-             } 
-           
+               //natasha
+               if(totaladc2[kcl] >= pow(totaladc[kcl],2))
+                 {
+                   sigxclust[kcl] = 0.25;
+                   sigyclust[kcl] = 0.25;
+                 }
+               else
+                 {
+                   sigxclust[kcl] = (totaladc[kcl]/(pow(totaladc[kcl],2)-totaladc2[kcl]))*ax2[kcl];
+                   sigyclust[kcl] = (totaladc[kcl]/(pow(totaladc[kcl],2)-totaladc2[kcl]))*ay2[kcl];
+                 }     
+             }
+             
              for(j = 0; j < cellCount[kcl]; j++) clno++; 
              
              if (clno >= 4608) 
@@ -862,6 +873,8 @@ void AliPMDClusteringV1::RefClust(Int_t incr, Double_t edepcell[])
              clusdata[1] = yclust[kcl];
              clusdata[2] = totaladc[kcl];
              clusdata[3] = ncell[kcl];
+
+
              if(sigxclust[kcl] > sigyclust[kcl]) 
                {
                  clusdata[4] = pow(sigxclust[kcl],0.5);
@@ -872,7 +885,7 @@ void AliPMDClusteringV1::RefClust(Int_t incr, Double_t edepcell[])
                  clusdata[4] = pow(sigyclust[kcl],0.5);
                  clusdata[5] = pow(sigxclust[kcl],0.5);
                }
-
+             
              clxy[0] = tc[kcl];
 
              Int_t Ncell=1;

diff --git a/PMD/AliPMDtracker.cxx b/PMD/AliPMDtracker.cxx
index 92eff59467ab2649aaeaf0d9755cda5d22548ef0..016f5fe2489ecaf31f24d22dd2c356feaf342d24 100644 (file)
--- a/PMD/AliPMDtracker.cxx
+++ b/PMD/AliPMDtracker.cxx
@@ -194,7 +194,7 @@ void AliPMDtracker::Clusters2Tracks(AliESDEvent *event)
          clusdata[4] = fPMDrecpoint->GetClusSigmaX();
          clusdata[5] = fPMDrecpoint->GetClusSigmaY();
          
-         if (clusdata[4] < -90. && clusdata[5] < -90.)
+         if (clusdata[4] >= 0. && clusdata[5] >= 0.)
            { 
              // extract the associated cell information
              branch1->GetEntry(ncrhit);