Modifications in SDD correction maps: new method to plot the corrections vs. drift...

diff --git a/ITS/AliITSCorrMapSDD.cxx b/ITS/AliITSCorrMapSDD.cxx
index a648187340e23bebb2d8289b63349c14fe78a206..75f62a1231aba06e3f4489f0778dc13b201713b9 100644 (file)
--- a/ITS/AliITSCorrMapSDD.cxx
+++ b/ITS/AliITSCorrMapSDD.cxx
@@ -56,7 +56,8 @@ Float_t AliITSCorrMapSDD::GetCorrection(Float_t z, Float_t x, AliITSsegmentation
   if(bina>nAnodes)  AliError("Wrong anode anumber!");
   if(bina>=nAnodesHybrid) bina-=nAnodesHybrid;
   Float_t stept = seg->Dx()*kMicronTocm/(Float_t)fNDriftPts;
-  Int_t bint = TMath::Abs((Int_t)(x/stept));
+  Float_t drLen= seg->Dx()*kMicronTocm-TMath::Abs(x);
+  Int_t bint = TMath::Abs((Int_t)(drLen/stept));
   if(bint==fNDriftPts) bint-=1;
   if(bint>=fNDriftPts) AliError("Wrong bin number along drift direction!");
   return kMicronTocm*GetCellContent(bina,bint);
@@ -75,6 +76,22 @@ TH2F* AliITSCorrMapSDD::GetMapHisto() const{
   return hmap;
 }
 //______________________________________________________________________
+TH1F* AliITSCorrMapSDD::GetMapProfile() const{
+  // Returns a TH1F with the projection of the map along drift coordinate
+  Char_t hname[50];
+  sprintf(hname,"p%s",GetName());
+  TH1F* hprof=new TH1F(hname,"",fNDriftPts,0.,35.);
+  for(Int_t iDr=0;iDr<fNDriftPts; iDr++){
+    Float_t meanval=0.;
+    for(Int_t iAn=0;iAn<fNAnodePts; iAn++){
+      meanval+=GetCellContent(iAn,iDr);
+    }
+    hprof->SetBinContent(iDr+1,meanval/fNAnodePts);
+  }
+  return hprof;
+  
+}
+//______________________________________________________________________
 TH1F* AliITSCorrMapSDD::GetResidualDistr(Float_t dmin, Float_t dmax) const{
   // Returns a TH1F histogram with distribution of residual
   Char_t hname[50];

diff --git a/ITS/AliITSCorrMapSDD.h b/ITS/AliITSCorrMapSDD.h
index f00d13592d1717acb69af385508f6eab7f2b6fbb..2f149c6c33e0c1881b7ebb24e45cdbbf0f9f15f8 100644 (file)
--- a/ITS/AliITSCorrMapSDD.h
+++ b/ITS/AliITSCorrMapSDD.h
@@ -72,6 +72,7 @@ class AliITSCorrMapSDD : public TNamed {
 
   Float_t GetCorrection(Float_t z, Float_t x, AliITSsegmentationSDD *seg);
   TH2F* GetMapHisto() const;
+  TH1F* GetMapProfile() const;
   TH1F* GetResidualDistr(Float_t dmin=-300., Float_t dmax=300.) const;