- fix number of slice allocation in ESD track
authorabercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 2 Jun 2009 09:24:47 +0000 (09:24 +0000)
committerabercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 2 Jun 2009 09:24:47 +0000 (09:24 +0000)
- fix warnings reported by Federico

TRD/AliTRDtrackV1.cxx
TRD/qaRec/AliTRDclusterResolution.cxx

index b9a92c3..ac1edff 100644 (file)
@@ -803,7 +803,7 @@ void AliTRDtrackV1::UpdateESDtrack(AliESDtrack *track)
   // pack the two numbers together and store them in the ESD
   track->SetTRDntracklets(nPID | (nTrk<<3));
   // allocate space to store raw PID signals dEdx & momentum
-  track->SetNumberOfTRDslices((nslices+2)*nPID);
+  track->SetNumberOfTRDslices((nslices+2)*AliTRDgeometry::kNlayer);
   // store raw signals
   Float_t p, sp; Double_t spd;
   for (Int_t ip = 0; ip < kNplane; ip++) {
index 16220a1..0ff75f0 100644 (file)
@@ -429,8 +429,10 @@ void AliTRDclusterResolution::Exec(Option_t *)
     x = (t+.5)*fgkTimeBinLength; // conservative approach !!
 
     // resolution as a function of y displacement from pad center
-    // only for phi equal exB
-    if(TMath::Abs(dydx-fExB) < kDtgPhi){
+    // only for phi equal exB and clusters close to cathode wire plane
+    // for ideal simulations time bins 4,5 and 6.
+    if(TMath::Abs(dydx-fExB) < kDtgPhi &&
+       TMath::Abs(x-0.675)<0.225){
       h3 = (TH3S*)arr0->At(AliTRDgeometry::GetLayer(det));
       h3->Fill(x, cli->GetYDisplacement(), dy);
     }
@@ -739,9 +741,9 @@ void AliTRDclusterResolution::ProcessSigma()
         gs.AddPoint(&tgg2, s2, s2e);
 
         if(!ggs) continue;
-        Int_t ip = ggs->GetN();
-        ggs->SetPoint(ip, tgg2, s2);
-        ggs->SetPointError(ip, 0., s2e);
+        Int_t jp = ggs->GetN();
+        ggs->SetPoint(jp, tgg2, s2);
+        ggs->SetPointError(jp, 0., s2e);
       }
       if(gs.Eval()) continue;
 
@@ -829,10 +831,9 @@ void AliTRDclusterResolution::ProcessMean()
         h1->Fit(&f, "QN");
 
         // Fill <Dy> = f(dydx - h*dzdx)
-        Int_t ip = gm->GetN();
-        gm->SetPoint(ip, tgl, f.GetParameter(1));
-        gm->SetPointError(ip, 0., f.GetParError(1));
-        ip++;
+        Int_t jp = gm->GetN();
+        gm->SetPoint(jp, tgl, f.GetParameter(1));
+        gm->SetPointError(jp, 0., f.GetParError(1));
       }
       if(gm->GetN()<4) continue;