Add some more parameters for stand alone cosmics reconstruction
authorcblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 2 Jul 2008 12:47:49 +0000 (12:47 +0000)
committercblume <cblume@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 2 Jul 2008 12:47:49 +0000 (12:47 +0000)
TRD/AliTRDrecoParam.cxx
TRD/AliTRDrecoParam.h
TRD/AliTRDtrackingChamber.cxx

index 5e261de3f38874d066df662d1bdf2bbf0ab544fd..0e558f693946f2877bf46c8232eb5b7f92846281 100644 (file)
@@ -54,6 +54,8 @@ AliTRDrecoParam::AliTRDrecoParam()
   ,fClusSigThresh(3.5)
   ,fTCnexp(1)
   ,fADCbaseline(0)
+  ,fNumberOfPresamples(0)
+  ,fNumberOfPostsamples(0)
 {
   //
   // Default constructor
@@ -95,6 +97,8 @@ AliTRDrecoParam::AliTRDrecoParam(const AliTRDrecoParam &ref)
   ,fClusSigThresh(ref.fClusSigThresh)
   ,fTCnexp(ref.fTCnexp)
   ,fADCbaseline(ref.fADCbaseline)
+  ,fNumberOfPresamples(ref.fNumberOfPresamples)
+  ,fNumberOfPostsamples(ref.fNumberOfPostsamples)
 {
   //
   // Copy constructor
index a12024a204a5c6d1517aacb9acba066387436c70..3ed298511783abd3e4bb4eeec3eb0baa9f21957a 100644 (file)
@@ -50,6 +50,8 @@ public:
   Double_t GetClusMaxThresh() const         { return fClusMaxThresh;   };
   Double_t GetClusSigThresh() const         { return fClusSigThresh;   };
   Int_t    GetTCnexp() const                { return fTCnexp;          };
+  Int_t     GetNumberOfPresamples()  const {return fNumberOfPresamples;}
+  Int_t    GetNumberOfPostsamples() const {return fNumberOfPostsamples;}
 
         
   static   AliTRDrecoParam *GetLowFluxParam();
@@ -63,6 +65,7 @@ public:
   Bool_t   IsVertexConstrained() const      { return TestBit(kVertexConstrained); }
 
 
+  void     SetFindableClusters(Double_t r) {fkFindable = r;}
   void     SetClusterSharing(Bool_t share = kTRUE)            { SetBit(kClusterSharing, share);  };
   void     SetPIDMethod(AliTRDpidMethod pid)                  { fkPIDMethod = pid; };
   void     SetSeeding(Bool_t so = kTRUE)             { SetBit(kSeeding, so); }
@@ -77,6 +80,8 @@ public:
   void     SetNexponential(Int_t nexp)                        { fTCnexp          = nexp;   };
   void     SetADCbaseline(Int_t base)                         { fADCbaseline     = base;   };
   inline void SetSysCovMatrix(Double_t *sys);
+  void     SetNumberOfPresamples(Int_t n) {fNumberOfPresamples = n;}
+  void     SetNumberOfPostsamples(Int_t n) {fNumberOfPostsamples = n;}
 
 private:
   enum{
@@ -123,6 +128,8 @@ private:
   
   // ADC parameter
   Int_t     fADCbaseline;            // ADC baseline to be subtracted
+  Int_t     fNumberOfPresamples;     // number of presamples 
+  Int_t     fNumberOfPostsamples;     // number of postsamples 
 
   ClassDef(AliTRDrecoParam, 4)       // Reconstruction parameters for TRD detector
 
index 1cc2349c27da947e0c1b4847daefad99a2262628..d06b2a7d9b481a726914529f477a3e0d231b2fed 100644 (file)
@@ -13,7 +13,7 @@
  * provided "as is" without express or implied warranty.                  *\r
  **************************************************************************/\r
 \r
-/* $Id: AliTRDtrackingDebug.cxx 23810 2008-02-08 09:00:27Z hristov $ */\r
+/* $Id: AliTRDtrackingChamber.cxx 23810 2008-02-08 09:00:27Z hristov $ */\r
 \r
 ////////////////////////////////////////////////////////////////////////////\r
 //                                                                        //\r
@@ -166,7 +166,9 @@ Bool_t AliTRDtrackingChamber::GetSeedingLayer(AliTRDchamberTimeBin *&fakeLayer,
        const Int_t kMaxRows = 16;\r
        const Int_t kMaxCols = 144;\r
        const Int_t kMaxPads = 2304;\r
-               \r
+       Int_t timeBinMin = AliTRDReconstructor::RecoParam()->GetNumberOfPresamples();\r
+       Int_t timeBinMax = AliTRDReconstructor::RecoParam()->GetNumberOfPostsamples();\r
+\r
        // Get the geometrical data of the chamber\r
        Int_t layer = geo->GetLayer(fDetector);\r
        Int_t stack = geo->GetStack(fDetector);\r
@@ -193,8 +195,8 @@ Bool_t AliTRDtrackingChamber::GetSeedingLayer(AliTRDchamberTimeBin *&fakeLayer,
                histogram[irs] = &hvals[irs*kMaxCols];\r
                sigmas[irs] = &svals[irs*kMaxCols];\r
        }\r
-       for(Int_t iTime = 0; iTime < kNTimeBins; iTime++){\r
-               if(!(nClusters = fTB[iTime].GetNClusters())) continue;\r
+       for(Int_t iTime = timeBinMin; iTime < kNTimeBins-timeBinMax; iTime++){\r
+    if(!(nClusters = fTB[iTime].GetNClusters())) continue;\r
                z0 = fTB[iTime].GetZ0();\r
                zl = fTB[iTime].GetDZ0();\r
                for(Int_t incl = 0; incl < nClusters; incl++){\r