]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDfeeParam.h
Use debug stream only if requested
[u/mrichter/AliRoot.git] / TRD / AliTRDfeeParam.h
index 96b632bece057a19f2cc2d593a6561ffdb4319ab..282483a3ae2b2a89f863e056939735fa91903042 100644 (file)
@@ -5,16 +5,21 @@
 
 /* $Id$ */
 
-//////////////////////////////////////////////////
-//                                              //
-//  TRD front end electronics parameters class  //
-//  Contains all FEE (MCM, TRAP, PASA) related  //
-//  parameters, constants, and mapping.         //
-//                                              //
-//////////////////////////////////////////////////
+////////////////////////////////////////////////////////////////////////////
+//                                                                        //
+//  TRD front end electronics parameters class                            //
+//  Contains all FEE (MCM, TRAP, PASA) related                            //
+//  parameters, constants, and mapping.                                   //
+//                                                                        //
+//  Author:                                                               //
+//    Ken Oyama (oyama@physi.uni-heidelberg.de)                           //
+//                                                                        //
+////////////////////////////////////////////////////////////////////////////
 
 #include <TObject.h>
 
+class TRootIoCtor;
+
 class AliTRDCommonParam;
 class AliTRDpadPlane;
 class AliTRDgeometry;
@@ -25,6 +30,7 @@ class AliTRDfeeParam : public TObject
 
  public:
 
+  AliTRDfeeParam(TRootIoCtor *);
   AliTRDfeeParam(const AliTRDfeeParam &p);
   virtual           ~AliTRDfeeParam();
   AliTRDfeeParam    &operator=(const AliTRDfeeParam &p);
@@ -42,35 +48,88 @@ class AliTRDfeeParam : public TObject
   virtual Int_t    GetColSide(Int_t icol) const;
 
   static  Float_t  GetSamplingFrequency() { return (Float_t)fgkLHCfrequency / 4000000.0; }
-  static  Int_t    GetNmcmRob()           { return fgkNmcmRob; }
+  static  Int_t    GetNmcmRob()           { return fgkNmcmRob;      }
   static  Int_t    GetNmcmRobInRow()      { return fgkNmcmRobInRow; }
   static  Int_t    GetNmcmRobInCol()      { return fgkNmcmRobInCol; }
-  static  Int_t    GetNrobC0()            { return fgkNrobC0; }
-  static  Int_t    GetNrobC1()            { return fgkNrobC1; }
-  static  Int_t    GetNadcMcm()           { return fgkNadcMcm; }
-  static  Int_t    GetNtimebin()          { return fgkNtimebin; }
-  static  Int_t    GetNcol()              { return fgkNcol; }
-  static  Int_t    GetNcolMcm()           { return fgkNcolMcm; }
-  static  Int_t    GetNrowC0()            { return fgkNrowC0; }
-  static  Int_t    GetNrowC1()            { return fgkNrowC1; }
-
-  //          Float_t  GetClusThr()           { return fClusThr; };
+  static  Int_t    GetNrobC0()            { return fgkNrobC0;       }
+  static  Int_t    GetNrobC1()            { return fgkNrobC1;       }
+  static  Int_t    GetNadcMcm()           { return fgkNadcMcm;      }
+  // static  Int_t    GetNtimebin()       { return fgkNtimebin;     }
+  static  Int_t    GetNcol()              { return fgkNcol;         }
+  static  Int_t    GetNcolMcm()           { return fgkNcolMcm;      }
+  static  Int_t    GetNrowC0()            { return fgkNrowC0;       }
+  static  Int_t    GetNrowC1()            { return fgkNrowC1;       }
+
+  // static  Int_t    GetADCpedestal()    { return fgkADCpedestal;  }
+  // static  Int_t    GetADCnoise()       { return fgkADCnoise;     }
+  static  Int_t    GetADCDAC()            { return fgkADCDAC;       }
+
+  static  Bool_t   IsPFon()               { return fgkPFon;         }
+  static  Bool_t   IsGFon()               { return fgkGFon;         }
+  static  Bool_t   IsTFon()               { return fgkTFon;         }
+
+  static  Int_t    GetPFtimeConstant()    { return fgkPFtimeConstant;   }
+  static  Int_t    GetPFeffectPedestal()  { return fgkPFeffectPedestal; }
+
+  //new
+  static  Int_t    GetQacc0Start()        {  return fgkPREPqAcc0Start; }
+  static  Int_t    GetQacc0End()          {  return fgkPREPqAcc0End; }
+  static  Int_t    GetQacc1Start()        {  return fgkPREPqAcc1Start; }
+  static  Int_t    GetQacc1End()          {  return fgkPREPqAcc1End; }
+          Float_t  GetMinClusterCharge() const {  return fgkMinClusterCharge; }
+  static  Int_t    GetLinearFitStart()    {  return fgkPREPLinearFitStart; }
+  static  Int_t    GetLinearFitEnd()      {  return fgkPREPLinearFitEnd;  }
+
+  //        Float_t  GetClusThr()           { return fClusThr; };
   //        Float_t  GetPadThr() const { return fPadThr; };
   //        Int_t    GetTailCancelation() const { return fTCOn; };
   //        Int_t    GetNexponential() const { return fTCnexp; };
   //virtual void     GetFilterParam(Float_t &r1, Float_t &r2, Float_t &c1, Float_t &c2, Float_t &ped) const;
   //        Int_t    GetFilterType() const { return fFilterType; };
 
-  static  Float_t  GetTFattenuationParam() { return ((Float_t)fgkTFattenuationParameter1) / ((Float_t)fgkTFattenuationParameter2) ; }
-  static  Float_t  GetTFf0()               { return 1 + fgkTFon*(-1+GetTFattenuationParam()); }   // 1 if TC off
+  static  Int_t    GetTFtype()            { return fgkTFtype;       }
+  //static  Int_t    GetTFnExp()            { return fgkTFnExp;       }
+          Int_t    GetTFnExp()            { return fTFnExp;         }
+          Float_t  GetTFr1()        const { return fTFr1;           }
+          Float_t  GetTFr2()        const { return fTFr2;           }
+          Float_t  GetTFc1()        const { return fTFc1;           }
+          Float_t  GetTFc2()        const { return fTFc2;           }
+
+ // for tracklets
+         Bool_t   GetTracklet()         const { return fgkTracklet; } 
+          Int_t    GetMaxNrOfTracklets() const { return fgkMaxNrOfTracklets; } 
+         Bool_t    GetMCTrackletOutput() const { return fgkMCTrackletOutput; }
+
+  static  Float_t  GetTFattPar()          { return ((Float_t) fgkTFattPar1) / ((Float_t) fgkTFattPar2); }
+          Float_t  GetTFf0()        const { return 1.0 + fgkTFon*(-1.0+GetTFattPar()); }   // 1 if TC off
+
+          void     SetEBsglIndThr(Int_t val);  
+          Int_t    GetEBsglIndThr() const { return fEBsglIndThr;    }
+
+          void     SetEBsumIndThr(Int_t val);
+          Int_t    GetEBsumIndThr() const { return fEBsumIndThr;    }
+
+          void     SetEBindLUT(Int_t val);
+          Int_t    GetEBindLUT()    const { return fEBindLUT;       }
+
+          void     SetEBignoreNeighbour(Int_t val);
+          Int_t    GetEBignoreNeighbour() const             { return fEBignoreNeighbour; }
+
+  // Concerning raw data format
+          Int_t    GetRAWversion() const                    { return fRAWversion;        }
+          void     SetRAWversion( Int_t rawver );
+          Bool_t   GetRAWstoreRaw() const                   { return fRAWstoreRaw;       }
+          void     SetRAWstoreRaw( Bool_t storeraw )        { fRAWstoreRaw = storeraw;   }
+
+          void     SetXenon();
+          void     SetArgon();
 
  protected:
 
-  static AliTRDfeeParam *fgInstance;
-  static Bool_t          fgTerminated;       //  Defines if this class has already been terminated                                                        
+  static AliTRDfeeParam *fgInstance;         // Singleton instance
+  static Bool_t          fgTerminated;       // Defines if this class has already been terminated
 
-  //  AliTRDgeometry    *fGeo;     // TRD geometry class
-  AliTRDCommonParam *fCP;      // TRD common parameters class
+  AliTRDCommonParam     *fCP;                // TRD common parameters class
 
   // Remark: ISO C++ allows initialization of static const values only for integer.
 
@@ -82,15 +141,13 @@ class AliTRDfeeParam : public TObject
   static const Int_t    fgkNrobC0            = 6;         // Number of ROBs per C0 chamber  (old fgkROBmaxC0)
   static const Int_t    fgkNrobC1            = 8;         // Number of ROBs per C1 chamber  (old fgkROBmaxC1)
   static const Int_t    fgkNadcMcm           = 21;        // Number of ADC channels per MCM (old fgkADCmax)
-  static const Int_t    fgkNtimebin          = 30;        // Number of Time bins            (old fgkTBmax)
+  // static const Int_t    fgkNtimebin       = 24;        // Number of Time bins should come from calibDB
   static const Int_t    fgkNcol              = 144;       // Number of pads per padplane row(old fgkColmax)
   static const Int_t    fgkNcolMcm           = 18;        // Number of pads per MCM         (old fgkPadmax)
   static const Int_t    fgkNrowC0            = 12;        // Number of Rows per C0 chamber  (old fgkRowmaxC0)
   static const Int_t    fgkNrowC1            = 16;        // Number of Rows per C1 chamber  (old fgkRowmaxC1)
 
   // ADC intrinsic parameters
-  static const Int_t    fgkADCpedestal       = 100000;    // ADC baseline * 100 (old name fPedestal)
-  static const Int_t    fgkADCnoise          = 10;        // ADC noise    * 100 (not contained in the digitizer) [in ADC] 
   static const Int_t    fgkADCDAC            = 0;         // 5 bit ADC gain parameter
 
   // TRAP filter global setup
@@ -106,46 +163,68 @@ class AliTRDfeeParam : public TObject
   static const Int_t    fgkGFnoise           =  0;        // Noise level increased by gain filter x 100 [in ADC] (to be measured)
 
   // TF setup
-  static const Int_t    fgkTFtype                 = 2;    // TC type (0=analog, 1=digital, 2=TRAPsim) (old name fFilterType)
-  static const Int_t    fgkTFlongDecayWeight      = 270;  // 0 to 1024 corresponds to 0 to 0.5
-  static const Int_t    fgkTFlongDecayParameter   = 348;  // 0 to 511 corresponds to 0.75 to 1
-  static const Int_t    fgkTFshortDecayParameter  = 449;  // 0 to 511 correponds to 0.25 to 0.5
-  static const Int_t    fgkTFattenuationParameter1= 45;   // attenuationParameter = fgkTFattenuationParameter1/fgkTFattenuationParameter2
-  static const Int_t    fgkTFattenuationParameter2= 14;   //                      = -alphaL/ln(lambdaL)-(1-alphaL)/ln(lambdaS)
-
-  // OLD TF setup (calculated from above)  (valid only for fgkTFtype = 0 or 1)
-  static const Int_t    fgkTFaNExp                 = 1;    // Number of exponential
-               Float_t  fTFaR1;                        // Time constant [microseconds] long (old name fR1)
-               Float_t  fTFaR2;                        // Time constant [microseconds] short(old name fR2)
-               Float_t  fTFaC1;                        // Weight long  (old name fC1)
-               Float_t  fTFaC2;                        // Weight short (old name fC2)
-
-  // Zero suppression parameters
-  static const Int_t    fgkEBsingleIndicatorThreshold = 3;    // used in EBIS, in ADC units above the pedestal
-  static const Int_t    fgkEBsumIndicatorThreshold    = 4;    // used in EBIT, in ADC units above the pedestal
-  static const Int_t    fgkEBindicatorLookupTable     = 0xF0; // see the TRAP user manual, used in EBIL
-  static const Int_t    fgkEBmarkIgnoreNeighbour      = 1;    // used in EBIN
+  static const Int_t    fgkTFtype            = 1;         // TC type (0=analog, 1=digital, 2=MI, 3=close to electronics) (old name fFilterType)
+
+  // OLD TF setup (calculated from above)  (valid only for fgkTFsimType = 0 or 1)
+  //static const Int_t    fgkTFnExp          = 1;           // Number of exponential for simType 0 and 1
+               Int_t    fTFnExp;                            // Number of exponential for simType 0 and 1
+
+ // Tracklet  processing on/off 
+  static const Bool_t   fgkTracklet         = kFALSE; // default should be kFALSE
+
+  // max. nr of tracklet words for one mcm
+  static const Int_t    fgkMaxNrOfTracklets = 4; 
+
+  // additional tracklet folder structure output, containing all necessary Monte Carlo information; maybe this should go somewhere else;
+  static const Bool_t   fgkMCTrackletOutput = kTRUE; // default should be kTRUE
+
+  // following need Instance because initialized in constructor
+               Float_t  fTFr1;                            // Time constant [us] long (old name fR1)
+               Float_t  fTFr2;                            // Time constant [us] short(old name fR2)
+               Float_t  fTFc1;                            // Weight long  (old name fC1)
+               Float_t  fTFc2;                            // Weight short (old name fC2)
+
+  // here is for TRAP simulation (not yet used)
+  static const Int_t    fgkTFdecayWeightL     = 270;      // 0 to 1024 corresponds to 0 to 0.5
+  static const Int_t    fgkTFdecayParL        = 348;      // 0 to 511 corresponds to 0.75 to 1
+  static const Int_t    fgkTFdecayParS        = 449;      // 0 to 511 correponds to 0.25 to 0.5
+  static const Int_t    fgkTFattPar1          = 45;       // attenuationParameter = fgkTFattenuationParameter1/fgkTFattenuationParameter2
+  static const Int_t    fgkTFattPar2          = 14;       //                      = -alphaL/ln(lambdaL)-(1-alphaL)/ln(lambdaS)
+
+  // ZS parameters
+               Int_t    fEBsglIndThr;                     // EBIS in ADC units
+               Int_t    fEBsumIndThr;                     // EBIT in ADC units
+               Int_t    fEBindLUT;                        // EBIL lookup table
+               Int_t    fEBignoreNeighbour;               // EBIN 0:include neighbor
 
   // Charge accumulators
-  static const Int_t    fgkPREPqAcc0Start             =  0;   // Preprocessor Charge Accumulator 0 Start
-  static const Int_t    fgkPREPqAcc0End               = 10;   // Preprocessor Charge Accumulator 0 End
-  static const Int_t    fgkPREPqAcc1Start             = 11;   // Preprocessor Charge Accumulator 1 Start
-  static const Int_t    fgkPREPqAcc1End               = 20;   // Preprocessor Charge Accumulator 1 End
-  static const Int_t    fgkMinClusterCharge           = 20;   // Hit detection [in ADC units]
+  static const Int_t    fgkPREPqAcc0Start     =  5;       // Preprocessor Charge Accumulator 0 Start
+  static const Int_t    fgkPREPqAcc0End       = 10;       // Preprocessor Charge Accumulator 0 End
+  static const Int_t    fgkPREPqAcc1Start     = 11;       // Preprocessor Charge Accumulator 1 Start
+  static const Int_t    fgkPREPqAcc1End       = 20;       // Preprocessor Charge Accumulator 1 End
+  static const Int_t    fgkMinClusterCharge   = 20;       // Hit detection [in ADC units]
+
+  //new
+  //time constants for linear fit
+  static const Int_t    fgkPREPLinearFitStart = 5;
+  static const Int_t    fgkPREPLinearFitEnd   = 20;
+
+
 
   // OLD TRAP processing parameters calculated from above
-  //static const Float_t  fClusThr;                     // Cluster threshold
-  //static const Float_t  fPadThr;                      // Pad threshold
+  //static const Float_t  fClusThr;                       // Cluster threshold
+  //static const Float_t  fPadThr;                        // Pad threshold
 
   // For raw production
-  static const Int_t    fgkRAWversion            = 1;         // Raw data production version
-  static const Bool_t   fgkRAWstoreRaw           = kTRUE;     // Store unfiltered data for raw data stream
+               Int_t    fRAWversion;                      // Raw data production version
+  static const Int_t    fgkMaxRAWversion      = 3;        // Maximum raw version number supported
+               Bool_t   fRAWstoreRaw;                     // Store unfiltered data for raw data stream
 
  private:
 
   AliTRDfeeParam();
 
-  ClassDef(AliTRDfeeParam,1)  //
-};
+  ClassDef(AliTRDfeeParam,3)                              // The TRD front end electronics parameter
 
+};
 #endif