Correction for SDD charge dependence on drift distance (E. Biolcati)

diff --git a/ITS/AliITSClusterFinderV2SDD.cxx b/ITS/AliITSClusterFinderV2SDD.cxx
index 29eeb00a76c744645cf0abe85619ef24d8a3a3a9..d8f971b4e3c5c9d2cc55c2d66a77b720ac08a5cd 100644 (file)
--- a/ITS/AliITSClusterFinderV2SDD.cxx
+++ b/ITS/AliITSClusterFinderV2SDD.cxx
@@ -250,7 +250,8 @@ FindClustersSDD(AliBin* bins[2], Int_t nMaxBin, Int_t nzBins,
        y=trk[1];
        z=trk[2]; 
 
-       q/=rsdd->GetADC2keV();  //to have MPV 1 MIP = 86.4 KeV
+       q/=rsdd->GetADC2keV();
+       q+=(driftTime*rsdd->GetChargevsTime()); // correction for zero supp.
        if(cal-> IsAMAt20MHz()) q*=2.; // account for 1/2 sampling freq.
        if(q<repa->GetMinClusterChargeSDD()) continue; // remove noise clusters
 

diff --git a/ITS/AliITSresponseSDD.cxx b/ITS/AliITSresponseSDD.cxx
index 118dd5dc981af904507e42f2714bdde517e8a023..1ab76ba4d93f60471e0f5fc3aa957dbfb62c7f96 100644 (file)
--- a/ITS/AliITSresponseSDD.cxx
+++ b/ITS/AliITSresponseSDD.cxx
@@ -30,7 +30,8 @@
 #include <AliITSgeomTGeo.h>
 
 const Float_t AliITSresponseSDD::fgkTimeOffsetDefault = 54.30;
-const Float_t AliITSresponseSDD::fgkADC2keVDefault = 3.25;
+const Float_t AliITSresponseSDD::fgkADC2keVDefault = 3.34;
+const Float_t AliITSresponseSDD::fgkChargevsTimeDefault = 0.00355;
 const Float_t AliITSresponseSDD::fgkCarlosRXClockPeriod = 25.;
 ClassImp(AliITSresponseSDD)
 
@@ -38,7 +39,8 @@ ClassImp(AliITSresponseSDD)
 AliITSresponseSDD::AliITSresponseSDD():
 TObject(),
 fTimeOffset(fgkTimeOffsetDefault),
-fADC2keV(fgkADC2keVDefault){
+fADC2keV(fgkADC2keVDefault),
+fChargevsTime(fgkChargevsTimeDefault){
   // default constructor
   for(Int_t i=0; i<kNSDDmods;i++){
     fTimeZero[i]=fgkTimeOffsetDefault;

diff --git a/ITS/AliITSresponseSDD.h b/ITS/AliITSresponseSDD.h
index b9d2f4e3349da051f77d21401cff15d19ec5e9de..0aba552e455e130831d474ff145309446cd9520d 100644 (file)
--- a/ITS/AliITSresponseSDD.h
+++ b/ITS/AliITSresponseSDD.h
@@ -73,8 +73,12 @@ class AliITSresponseSDD : public TObject {
   static Float_t DefaultTimeOffset() {return fgkTimeOffsetDefault;}
 
   virtual void SetADC2keV(Float_t conv){fADC2keV=conv;}
+  virtual void SetChargevsTime(Float_t slope){fChargevsTime=slope;}
   virtual Float_t GetADC2keV()const {return fADC2keV;}
-  static Float_t DefaulttADC2keV() {return fgkADC2keVDefault;}
+  virtual Float_t GetChargevsTime()const {return fChargevsTime;}
+
+  static Float_t DefaultADC2keV() {return fgkADC2keVDefault;}
+  static Float_t DefaultChargevsTime() {return fgkChargevsTimeDefault;}
 
   static Float_t GetCarlosRXClockPeriod() {return fgkCarlosRXClockPeriod;}
  
@@ -89,6 +93,7 @@ class AliITSresponseSDD : public TObject {
 
   static const Float_t fgkTimeOffsetDefault;   // default for fTimeOffset
   static const Float_t fgkADC2keVDefault;      // default for fADC2keV
+  static const Float_t fgkChargevsTimeDefault; // default for fChargevsTime
   static const Float_t fgkCarlosRXClockPeriod; // clock period for CarlosRX
 
   Float_t  fTimeOffset;             // Time offset due to electronic delays 
@@ -96,13 +101,14 @@ class AliITSresponseSDD : public TObject {
   Float_t  fTimeZero[kNSDDmods];    // Time Zero for each module
   Float_t  fDeltaVDrift[kNSDDmods]; // Vdrift correction (um/ns) for each module
   Float_t  fADC2keV;                // Conversion factor from ADC to keV
-
+  Float_t  fChargevsTime;           // Correction for zero suppression effect
+  
  private:
 
   AliITSresponseSDD(const AliITSresponseSDD &ob); // copy constructor
   AliITSresponseSDD& operator=(const AliITSresponseSDD & /* source */); // ass. op.
 
-  ClassDef(AliITSresponseSDD,17) 
+  ClassDef(AliITSresponseSDD,18) 
      
     };
 #endif