Automatically set the lower pt cut for hough tracks according to the magnetic field...

diff --git a/MONITOR/monitorGDC.cxx b/MONITOR/monitorGDC.cxx
index 586ff5275f1e1239a3fdb8d33d50c60dcb6e3f4c..f6d8928fc4813470846afca5a5ef5bcfc2edb533 100644 (file)
--- a/MONITOR/monitorGDC.cxx
+++ b/MONITOR/monitorGDC.cxx
@@ -135,12 +135,14 @@ int main(int argc, char** argv)
        AliL3Benchmark *fBenchmark = new AliL3Benchmark();
        fBenchmark->Start("Overall timing");
 
+       Float_t ptmin = 0.1*AliL3Transform::GetSolenoidField();
+
        // Run the Hough Transformer
        fBenchmark->Start("Init");
        AliL3Hough *hough1 = new AliL3Hough();
 
        hough1->SetThreshold(4);
-       hough1->CalcTransformerParams(0.4);
+       hough1->CalcTransformerParams(ptmin);
        hough1->SetPeakThreshold(70,-1);
        //      printf("Pointer is %x\n",ptr);
        hough1->Init("./", kFALSE, 100, kFALSE,4,0,(Char_t*)ptr,0.0);
@@ -151,7 +153,7 @@ int main(int argc, char** argv)
        AliL3Hough *hough2 = new AliL3Hough();
 
        hough2->SetThreshold(4);
-       hough2->CalcTransformerParams(0.4);
+       hough2->CalcTransformerParams(ptmin);
        hough2->SetPeakThreshold(70,-1);
        //      printf("Pointer is %x\n",ptr);
        hough2->Init("./", kFALSE, 100, kFALSE,4,0,(Char_t*)ptr,0.0);

diff --git a/RAW/AliMDC.cxx b/RAW/AliMDC.cxx
index 99a7ead319c76e091c14b26d395a31d7eb83150c..1206823e7dad343ae3ad908eeaec171b7d174034 100644 (file)
--- a/RAW/AliMDC.cxx
+++ b/RAW/AliMDC.cxx
@@ -806,10 +806,12 @@ Int_t AliMDC::Filter(
     TStopwatch timer;
     timer.Start();
 
+    Float_t ptmin = 0.1*AliL3Transform::GetSolenoidField();
+
     AliL3Hough *hough1 = new AliL3Hough();
     
     hough1->SetThreshold(4);
-    hough1->CalcTransformerParams(0.4);
+    hough1->CalcTransformerParams(ptmin);
     hough1->SetPeakThreshold(70,-1);
     // Attention Z of the vertex to be taken from the event head!
     // So far for debug purposes it is fixed by hand...
@@ -819,7 +821,7 @@ Int_t AliMDC::Filter(
     AliL3Hough *hough2 = new AliL3Hough();
 
     hough2->SetThreshold(4);
-    hough2->CalcTransformerParams(0.4);
+    hough2->CalcTransformerParams(ptmin);
     hough2->SetPeakThreshold(70,-1);
     hough2->Init(100,4,event,3.82147);
     hough2->SetAddHistograms();