Move to new baseline subtraction scheme in the clusterizer

diff --git a/TRD/AliTRDarrayADC.cxx b/TRD/AliTRDarrayADC.cxx
index 54ecdad804ecd23d2ca0c301fd4b222a087ded45..5aaf9e285049d4e0dead1d5920bdb2b14049f72a 100644 (file)
--- a/TRD/AliTRDarrayADC.cxx
+++ b/TRD/AliTRDarrayADC.cxx
@@ -622,3 +622,19 @@ void AliTRDarrayADC::CreateLut()
        }\r
     }\r
 }\r
+//________________________________________________________________________________\r
+void AliTRDarrayADC::SubtractBaseline(Short_t baseline)\r
+{\r
+  //\r
+  // Subtracts the baseline to all the adc array\r
+  // if the value is equal or greater than baseline\r
+  //\r
+\r
+  for(Int_t b=0; b<fNAdim; b++)\r
+    {\r
+      //if(fADC[b]>=baseline)\r
+      //{\r
+         fADC[b]=fADC[b]-baseline;\r
+      //}\r
+    }\r
+}\r

diff --git a/TRD/AliTRDarrayADC.h b/TRD/AliTRDarrayADC.h
index b7249bc39158c8c8dd3d4940caa6c8f8f1ab57c7..e7c6e7cbf9662d34de6c9f77931da2eb007ba3b5 100644 (file)
--- a/TRD/AliTRDarrayADC.h
+++ b/TRD/AliTRDarrayADC.h
@@ -47,6 +47,7 @@ class AliTRDarrayADC: public TObject
   Short_t GetData(Int_t nrow, Int_t ncol, Int_t ntime) const;\r
   void    SetData(Int_t nrow, Int_t ncol, Int_t ntime, Short_t value);\r
   static  void    CreateLut(); \r
+  void    SubtractBaseline(Short_t baseline);\r
 \r
  protected:\r
 \r

diff --git a/TRD/AliTRDclusterizer.cxx b/TRD/AliTRDclusterizer.cxx
index 3abb0a910fe23c13be8446d6dbc631c61dfb6f21..defdcbc461deef77edb54d9786ff2c15bae89bbf 100644 (file)
--- a/TRD/AliTRDclusterizer.cxx
+++ b/TRD/AliTRDclusterizer.cxx
@@ -35,6 +35,7 @@
 #include "AliTRDarrayDictionary.h"
 #include "AliTRDarrayADC.h"
 #include "AliTRDdigitsManager.h"
+#include "AliTRDdigitsParam.h"
 #include "AliTRDrawData.h"
 #include "AliTRDcalibDB.h"
 #include "AliTRDtransform.h"
@@ -561,6 +562,9 @@ Bool_t AliTRDclusterizer::MakeClusters()
     // This is to take care of switched off super modules
     if (!digitsIn->HasData()) continue;
     digitsIn->Expand();
+    // Subtract the ADC baseline
+    printf("----> baseline = %d\n",fDigitsManager->GetDigitsParam()->GetADCbaseline());
+    digitsIn->SubtractBaseline(fDigitsManager->GetDigitsParam()->GetADCbaseline());
     digitsIn->DeleteNegatives();  // Restore digits array to >=0 values
     AliTRDSignalIndex* indexes = fDigitsManager->GetIndexes(i);
     if (indexes->IsAllocated() == kFALSE){
@@ -735,6 +739,10 @@ Bool_t AliTRDclusterizer::MakeClusters(Int_t det)
   // This is to take care of switched off super modules
   if (!fDigits->HasData()) return kFALSE;
 
+  // Subtract the ADC baseline
+  printf("----> baseline = %d\n",fDigitsManager->GetDigitsParam()->GetADCbaseline());
+  fDigits->SubtractBaseline(fDigitsManager->GetDigitsParam()->GetADCbaseline());
+
   fIndexes = fDigitsManager->GetIndexes(det);
   if (fIndexes->IsAllocated() == kFALSE) {
     AliError("Indexes do not exist!");