Automatically set the lower pt cut for hough tracks according to the magnetic field...
authorcvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 8 Nov 2004 12:48:01 +0000 (12:48 +0000)
committercvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 8 Nov 2004 12:48:01 +0000 (12:48 +0000)
MONITOR/monitorGDC.cxx
RAW/AliMDC.cxx

index 586ff52..f6d8928 100644 (file)
@@ -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);
index 99a7ead..1206823 100644 (file)
@@ -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();