FMD/AliFMDParameters.h

   1 #ifndef ALIFMDPARAMETERS_H
   2 #define ALIFMDPARAMETERS_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights
   4  * reserved.
   5  *
   6  * Latest changes by Christian Holm Christensen <cholm@nbi.dk>
   7  *
   8  * See cxx source for full Copyright notice
   9  */
  10 //____________________________________________________________________
  11 //
  12 //  Singleton class to handle various parameters (not geometry) of the
  13 //  FMD
  14 //  Should get ata fromm Conditions DB.
  15 //
  16 /** @file    AliFMDParameters.h
  17     @author  Christian Holm Christensen <cholm@nbi.dk>
  18     @date    Mon Mar 27 12:44:43 2006
  19     @brief   Manager of FMD parameters
  20 */
  21 #ifndef ROOT_TNamed
  22 # include <TNamed.h>
  23 #endif
  24 #ifndef ROOT_TArrayI
  25 # include <TArrayI.h>
  26 #endif
  27 #include "AliFMDCalibFwd.h"
  28 class AliCDBEntry;
  29 class AliFMDPreprocessor;
  30
  31 //____________________________________________________________________
  32 //
  33 //  Singleton class to handle various parameters (not geometry) of the
  34 //  FMD
  35 //  Should get ata fromm Conditions DB.
  36 //
  37
  38 /** @brief This class is a singleton that handles various parameters
  39     of the FMD detectors.
  40     This class reads from the Conditions DB to get the various
  41     parameters, which code can then request from here. In that way,
  42     all code uses the same data, and the interface is consistent.
  43
  44     Some of the parameter managed are
  45     - @c fPedestal, @c fPedestalWidth
  46       Mean and width of the pedestal.  The pedestal is simulated
  47       by a Guassian, but derived classes my override MakePedestal
  48       to simulate it differently (or pick it up from a database).
  49     - @c fVA1MipRange
  50       The dymamic MIP range of the VA1_ALICE pre-amplifier chip
  51     - @c fAltroChannelSize
  52       The largest number plus one that can be stored in one
  53       channel in one time step in the ALTRO ADC chip.
  54     - @c fSampleRate
  55       How many times the ALTRO ADC chip samples the VA1_ALICE
  56       pre-amplifier signal.   The VA1_ALICE chip is read-out at
  57       10MHz, while it's possible to drive the ALTRO chip at
  58       25MHz.  That means, that the ALTRO chip can have time to
  59       sample each VA1_ALICE signal up to 2 times.  Although it's
  60       not certain this feature will be used in the production,
  61       we'd like have the option, and so it should be reflected in
  62       the code.
  63
  64     @ingroup FMD_base
  65 */
  66 class AliFMDParameters : public TNamed
  67 {
  68 public:
  69   /** Enumeration of things to initialize */
  70   enum What {
  71     /** Pulser gain */
  72     kPulseGain = 0x1, // Pulser gain
  73     /** Pedestals and noise */
  74     kPedestal = 0x2, // Pedestal and noise
  75     /** Dead channel map */
  76     kDeadMap = 0x4, // Dead channel map
  77     /**  Over sampling rate */
  78     kSampleRate = 0x8, // Over sampling rate
  79     /** Zero suppression parameters */
  80     kZeroSuppression = 0x10, // Zero suppression parameters
  81     /** ALTRO data map */
  82     kAltroMap = 0x20, // Altro channel map
  83     /** Strip Range */
  84     kStripRange = 0x40 //Strip range,
  85   };
  86   enum {
  87     kAll = (kPulseGain|kPedestal|kDeadMap|kSampleRate|
  88             kZeroSuppression|kAltroMap|kStripRange)
  89   };
  90
  91   /**
  92    * Singleton access
  93    *
  94    *
  95    * @return singleton
  96    */
  97   static AliFMDParameters* Instance();
  98
  99   /**
 100    * Initialize the manager.  This tries to read the parameters from
 101    * CDB.  If that fails, the class uses the hard-coded parameters.
 102    *
 103    * @param forceReInit Force (re-)initalize flag
 104    * @param what        What to initialize
 105    */
 106   void Init(Bool_t forceReInit=kFALSE, UInt_t what=kAll );
 107   /**
 108    * Initialize the manager.  This tries to read the parameters from
 109    * CDB.  If that fails, the class uses the hard-coded parameters.
 110    *
 111    * @param pp          Preprocessor
 112    * @param forceReInit Force (re-)initalize flag
 113    * @param what        What to initialize
 114    */
 115   void Init(AliFMDPreprocessor* pp,
 116             Bool_t              forceReInit=kFALSE,
 117             UInt_t              what=kAll);
 118   /**
 119    * Initialize the manager.  This will try to read some calibrations
 120    * (sample rate, strip range, gains, pedestals) from local comma
 121    * separated value (CSV) files in the directory pointed at by @a
 122    * path.  If they are not found, then they will be retrieved from
 123    * OCDB as appropriately.   Other calibrations are always read from
 124    * OCDB.
 125    *
 126    * The CSV files should be named as
 127    *
 128    * - Pedestals: <tt>peds</tt><i>det_number</i><tt>.csv</tt>
 129    * - Gains: <tt>gains</tt><i>det_number</i><tt>.csv</tt>
 130    * - Sample Rate: <tt>conditions</tt><i>det_number</i><tt>.csv</tt>
 131    * - Strip Range: <tt>conditions</tt><i>det_number</i><tt>.csv</tt>
 132    *
 133    * where <i>det_number</i> is the detector number (1, 2, or 3).
 134    *
 135    * @param path        Where to look for the CSV files
 136    * @param forceReInit Always reinitialise
 137    * @param what        What calibrations to load.
 138    */
 139   void Init(const char* path,
 140             Bool_t      forceReInit=kFALSE,
 141             UInt_t      what=kAll);
 142
 143   /**
 144    * Automatically generate a dead map from the pedestals and gains.
 145    * A channel is marked as dead of the noise is too high (currently
 146    * more than 10 ADC counts), or the gain is unreasonable (currently
 147    * larger than 10, or smaller than 0.1).
 148    *
 149    * The procedure does not overwrite channels previously marked as
 150    * dead - e.g., channels marked as dead in the calibration loaded
 151    * from OCDB will continue to be marked as dead.  That is, this
 152    * procedure will never make a channel un-dead.
 153    *
 154    * @param maxNoise  Maximum noise value before a channel is marked
 155    * as dead.
 156    * @param minGain   Minimum value of the calibrated gain before a
 157    * channel is considered dead.
 158    * @param maxGain   Maximum value of the calibrated gain before a
 159    * channel is considered dead.
 160    */
 161   void MakeDeadMap(Float_t maxNoise=10, Float_t minGain=.1, Float_t maxGain=10);
 162   /**
 163    * Print all parameters.
 164    *
 165    * If option contains an 'A' then everything is printed.
 166    *
 167    * If the option contains the string "FMD" the function will search
 168    * for detector, ring, sector, and strip numbers to print, in
 169    * format
 170    *
 171    * @verbatim
 172    *    FMD<detector><ring>[<sector>,<string>]
 173    * @endverbatim
 174    *
 175    * The wild card '*' means all of <detector>, <ring>, <sector>, or
 176    * <strip>.
 177    *
 178    * @param option Option string
 179    */
 180   void Print(Option_t* option="A") const;
 181   /**
 182    * Draw parameters.
 183    *
 184    * @param option What to draw. Should be one of
 185    * - dead       Dead channels
 186    * - threshold Threshold
 187    * - gain       Gain
 188    * - pedestal  Pedestal
 189    * - noise      Noise (or pedestal width)
 190    * - zero       Zero suppression
 191    * - rate       Sampling rate (VA1 clock / ALTRO clock)
 192    * - min        Minimum strip read out
 193    * - max        Maximum strip read out
 194    * - map        hardware address
 195    */
 196   void Draw(Option_t* option="pedestal");
 197
 198   /** @{ */
 199   /** @name Set various `Fixed' parameters */
 200   /**
 201    * @param r How many MIP signals we can fit in the VA1
 202    * pre-amps. (default and design is 20)
 203    */
 204   void SetVA1MipRange(UShort_t r=20)          { fVA1MipRange = r; }
 205   /**
 206    * @param s Maximum number of the ADC (ALTRO).  This is a 10 bit
 207    * ADC so, the maximum number is 1024
 208    */
 209   void SetAltroChannelSize(UShort_t s=1024)   { fAltroChannelSize = s;}
 210   /**
 211    * @param size The number of strips multiplexed into one ALTRO
 212    * channel. That is, how many strips is connected to one VA1
 213    * pre-amp.
 214    */
 215   void SetChannelsPerAltro(UShort_t size=128) { fChannelsPerAltro = size; }
 216   /**
 217    * @param f Factor to use for accepting a signal.
 218    */
 219   void SetPedestalFactor(Float_t f=3)         { fPedestalFactor = f; }
 220   /**
 221    * @param n Number of pre-samples to keep during zero-suppression -
 222    * only used in simulation.
 223    */
 224   void SetZSPreSamples(UShort_t n=1) { fZSPre = (n & 0x3); }
 225   /**
 226    * @param n Number of post-samples to keep during zero-suppression -
 227    * only used in simulation.
 228    */
 229   void SetZSPostSamples(UShort_t n=1) { fZSPost = (n & 0x3); }
 230   /**
 231    * @param use If true, do pedestal subtraction before zero
 232    * suppression - only used in simulation
 233    */
 234   void SetZSPedSubtract(Bool_t use=kTRUE) { fZSPedSubtract = use; }
 235   /** @} */
 236
 237   /** @{ */
 238   /** @name Set various variable parameter defaults */
 239   /**
 240    * @param s Zero suppression threshold in ADC counts
 241    */
 242   void SetZeroSuppression(UShort_t s=0)       { fFixedZeroSuppression = s; }
 243   /**
 244    * @param r How many times we oversample each strip.
 245    */
 246   void SetSampleRate(UShort_t r=1)            { fFixedSampleRate = r ;}//(r>2?2:r);}
 247   /**
 248    * @param r How many times we oversample each strip.
 249    */
 250   void SetSampleRate(AliFMDCalibSampleRate* r) { fSampleRate = r; }
 251   /**
 252    * @param p Pedestal value in ADC counts
 253    */
 254   void SetPedestal(Float_t p=10)              { fFixedPedestal = p; }
 255   /**
 256    * @param p Pedestal map
 257    */
 258   void SetPedestal(AliFMDCalibPedestal* p) { fPedestal = p; }
 259   /**
 260    * @param w Pedestal width in ADC counts
 261    */
 262   void SetPedestalWidth(Float_t w=1)          { fFixedPedestalWidth = w; }
 263   /**
 264    * @param t Threshold used for 1 MIP acceptance.
 265    */
 266   void SetThreshold(Float_t t=0)              { fFixedThreshold = t; }
 267   /**
 268    * Range of strips read out
 269    *
 270    * @param min Minimum strip number (0-127).
 271    * @param max Maximum strip number (0-127).
 272    */
 273   void SetStripRange(UShort_t min=0, UShort_t max=127);
 274   /**
 275    * set the strip range from object
 276    *
 277    * @param r Strip range object
 278    */
 279   void SetStripRange(AliFMDCalibStripRange* r) { fStripRange = r; }
 280   /**
 281    * Whether raw data has full common data header (8 32bit words) or
 282    * the older invalid format (7 32bit words with bogus entries)
 283    *
 284    * @param yes if true the raw data has complete data header
 285    */
 286   void UseCompleteHeader(Bool_t yes=kTRUE) { fHasCompleteHeader = yes; }
 287   /**
 288    * @param g Gain map
 289    */
 290   void SetGain(AliFMDCalibGain* g) { fPulseGain = g; }
 291   /** @} */
 292
 293   /** @{ */
 294   /** @name Get `Fixed' various parameters */
 295   /**
 296    * @return Number of MIP signals that fit inside a VA1 channel
 297    */
 298   UShort_t GetVA1MipRange()          const { return fVA1MipRange; }
 299   /**
 300    * @return The maximum count in the ADC
 301    */
 302   UShort_t GetAltroChannelSize()     const { return fAltroChannelSize; }
 303   /**
 304    * @return Number of strips muliplexed into one ADC channel
 305    */
 306   UShort_t GetChannelsPerAltro()     const { return fChannelsPerAltro; }
 307   /**
 308    * @return The average energy deposited by one MIP
 309    */
 310   Float_t  GetEdepMip()              const;
 311   /**
 312    * This is the conversion from Digital-to-Analog-Converter setting
 313    * to the number of MIPs. The number was measured in the NBI lab during
 314    * August 2008.
 315    *
 316    * @return The conversion factor from DAC to ADC
 317    */
 318   Float_t  GetDACPerMIP()              const;
 319   /**
 320    * @return The factor used of signal acceptance
 321    */
 322   Float_t  GetPedestalFactor()       const { return fPedestalFactor; }
 323   /**
 324    * @param n Number of pre-samples to keep during zero-suppression -
 325    * only used in simulation.
 326    */
 327   UShort_t GetZSPreSamples() const { return fZSPre; }
 328   /**
 329    * @param n Number of post-samples to keep during zero-suppression -
 330    * only used in simulation.
 331    */
 332   UShort_t GetZSPostSamples() const { return fZSPost; }
 333   /**
 334    * @param use If true, do pedestal subtraction before zero
 335    * suppression - only used in simulation
 336    */
 337   Bool_t IsZSPedSubtract() const { return fZSPedSubtract; }
 338   /** @} */
 339
 340   /** @{ */
 341   /** @name Various varible conditions */
 342   /**
 343    * Whether the strip is considered dead
 344    *
 345    * @param detector Detector # (1-3)
 346    * @param ring     Ring ID ('I' or 'O')
 347    * @param sector   Sector number (0-39)
 348    * @param strip    Strip number (0-511)
 349    *
 350    * @return @c true if the strip is considered dead, @c false if it's
 351    * OK.
 352    */
 353   Bool_t   IsDead(UShort_t detector,
 354                   Char_t ring,
 355                   UShort_t sector,
 356                   UShort_t strip) const;
 357   /**
 358    * Get the threshold in the pulser gain
 359    *
 360    *
 361    * @return Threshold from pulser
 362    */
 363   Float_t  GetThreshold() const;
 364   /**
 365    * Gain of pre-amp. for strip, sector, ring, detector
 366    *
 367    * For simulations this is normally set to
 368    *
 369    * @f[
 370    *  \frac{\mbox{VA1_MIP_Range}{\mbox{ALTRO_channel_size}}\mbox{MIP_Energy_Loss}
 371    * @f]
 372    *
 373    *
 374    * @param detector Detector # (1-3)
 375    * @param ring     Ring ID ('I' or 'O')
 376    * @param sector   Sector number (0-39)
 377    * @param strip    Strip number (0-511)
 378    *
 379    * @return Gain of pre-amp.
 380    */
 381   Float_t  GetPulseGain(UShort_t detector,
 382                         Char_t ring,
 383                         UShort_t sector,
 384                         UShort_t strip) const;
 385   /**
 386    * Get mean of pedestal
 387    *
 388    * @param detector Detector # (1-3)
 389    * @param ring     Ring ID ('I' or 'O')
 390    * @param sector   Sector number (0-39)
 391    * @param strip    Strip number (0-511)
 392    *
 393    * @return Mean of pedestal
 394    */
 395   Float_t  GetPedestal(UShort_t detector,
 396                        Char_t ring,
 397                        UShort_t sector,
 398                        UShort_t strip) const;
 399   /**
 400    * Width of pedestal
 401    *
 402    * @param detector Detector # (1-3)
 403    * @param ring     Ring ID ('I' or 'O')
 404    * @param sector   Sector number (0-39)
 405    * @param strip    Strip number (0-511)
 406    *
 407    * @return Width of pedestal
 408    */
 409   Float_t  GetPedestalWidth(UShort_t detector,
 410                             Char_t ring,
 411                             UShort_t sector,
 412                             UShort_t strip) const;
 413   /**
 414    * zero suppression threshold (in ADC counts)
 415    *
 416    * @param detector Detector # (1-3)
 417    * @param ring     Ring ID ('I' or 'O')
 418    * @param sector   Sector number (0-39)
 419    * @param strip    Strip number (0-511)
 420    *
 421    * @return zero suppression threshold (in ADC counts)
 422    */
 423   UShort_t GetZeroSuppression(UShort_t detector,
 424                               Char_t ring,
 425                               UShort_t sector,
 426                               UShort_t strip) const;
 427   /**
 428    * Get the sampling rate
 429    *
 430    * @param detector Detector # (1-3)
 431    * @param ring     Ring ID ('I' or 'O')
 432    * @param sector   Sector number (0-39)
 433    * @param strip    Strip number (0-511)
 434    *
 435    * @return The sampling rate
 436    */
 437   UShort_t GetSampleRate(UShort_t detector,
 438                          Char_t ring,
 439                          UShort_t sector,
 440                          UShort_t strip) const;
 441   /**
 442    * Get the minimum strip in the read-out range
 443    *
 444    * @param detector Detector # (1-3)
 445    * @param ring     Ring ID ('I' or 'O')
 446    * @param sector   Sector number (0-39)
 447    * @param strip    Strip number (0-511)
 448    *
 449    * @return Minimum strip
 450    */
 451   UShort_t GetMinStrip(UShort_t detector,
 452                        Char_t ring,
 453                        UShort_t sector,
 454                        UShort_t strip) const;
 455   /**
 456    * Get the maximum strip in the read-out range
 457    *
 458    * @param detector Detector # (1-3)
 459    * @param ring     Ring ID ('I' or 'O')
 460    * @param sector   Sector number (0-39)
 461    * @param strip    Strip number (0-511)
 462    *
 463    * @return Maximum strip
 464    */
 465   UShort_t GetMaxStrip(UShort_t detector,
 466                        Char_t ring,
 467                        UShort_t sector,
 468                        UShort_t strip) const;
 469   /**
 470    * Get the number of pre-samples in ALTRO channels
 471    *
 472    * @param detector Detector # (1-3)
 473    * @param ring     Ring ID ('I' or 'O')
 474    * @param sector   Sector number (0-39)
 475    * @param strip    Strip number (0-511)
 476    *
 477    * @return Maximum strip
 478    */
 479   UShort_t GetPreSamples(UShort_t,
 480                          Char_t,
 481                          UShort_t,
 482                          UShort_t) const { return 14+5; }
 483   /* @}*/
 484
 485   /**
 486    * @{
 487    * @name Hardware to detector translation (and inverse)
 488    */
 489   /**
 490    * Map a hardware address into a detector index.
 491    *
 492    * @param ddl        Hardware DDL number
 493    * @param board      FEC number
 494    * @param altro      ALTRO number
 495    * @param channel    Channel number
 496    * @param timebin    Timebin
 497    * @param det        On return, the detector #
 498    * @param ring       On return, the ring ID
 499    * @param sec        On return, the sector #
 500    * @param str        On return, the base of strip #
 501    * @param sam        On return, the sample number for this strip
 502    *
 503    * @return @c true on success, false otherwise
 504    */
 505   Bool_t Hardware2Detector(UShort_t    ddl,        UShort_t    board,
 506                            UShort_t    altro,      UShort_t    chan,
 507                            UShort_t  timebin,
 508                            UShort_t& det,        Char_t&   ring,
 509                            UShort_t& sec,        Short_t& str,
 510                            UShort_t& sam) const;
 511   /**
 512    * Map a hardware address into a detector index.
 513    *
 514    * @param ddl        Hardware DDL number
 515    * @param hwaddr     Hardware address.
 516    * @param timebin    Timebin
 517    * @param det        On return, the detector #
 518    * @param ring       On return, the ring ID
 519    * @param sec        On return, the sector #
 520    * @param str        On return, the base of strip #
 521    * @param sam        On return, the sample number for this strip
 522    *
 523    * @return @c true on success, false otherwise
 524    */
 525   Bool_t Hardware2Detector(UShort_t    ddl,        UShort_t    hwaddr,
 526                            UShort_t  timebin,
 527                            UShort_t& det,        Char_t&   ring,
 528                            UShort_t& sec,        Short_t& str,
 529                            UShort_t& sam) const;
 530
 531   /**
 532    * Map a detector index into a hardware address.
 533    *
 534    * @param det         The detector #
 535    * @param ring        The ring ID
 536    * @param sec         The sector #
 537    * @param str         The strip #
 538    * @param sam         The sample number
 539    * @param ddl         On return, hardware DDL number
 540    * @param board       On return, the FEC board address (local to DDL)
 541    * @param altro       On return, the ALTRO number (local to FEC)
 542    * @param channel     On return, the channel number (local to ALTRO)
 543    * @param timebin     On return, the timebin number (local to ALTRO)
 544    *
 545    * @return @c true on success, false otherwise
 546    */
 547   Bool_t Detector2Hardware(UShort_t  det,        Char_t    ring,
 548                            UShort_t  sec,        UShort_t  str,
 549                            UShort_t  sam,
 550                            UShort_t&   ddl,        UShort_t&   board,
 551                            UShort_t&   altro,      UShort_t&   channel,
 552                            UShort_t& timebin) const;
 553   /**
 554    * Map a detector index into a hardware address.
 555    *
 556    * @param det         The detector #
 557    * @param ring        The ring ID
 558    * @param sec         The sector #
 559    * @param str         The strip #
 560    * @param sam         The sample number
 561    * @param ddl         On return, hardware DDL number
 562    * @param hwaddr      On return, hardware address.
 563    * @param timebin     On return, the timebin number (local to ALTRO)
 564    *
 565    * @return @c true on success, false otherwise
 566    */
 567   Bool_t Detector2Hardware(UShort_t  det,        Char_t    ring,
 568                            UShort_t  sec,        UShort_t  str,
 569                            UShort_t  sam,
 570                            UShort_t&   ddl,        UShort_t&   hwaddr,
 571                            UShort_t& timebin) const;
 572   /**
 573    * Get the map that translates hardware to detector coordinates
 574    *
 575    * @return Get the map that translates hardware to detector
 576    * coordinates
 577    */
 578   AliFMDAltroMapping* GetAltroMap() const;
 579   /**
 580    * Whether raw data has full common data header (8 32bit words) or
 581    * the older invalid format (7 32bit words with bogus entries)
 582    *
 583    * @return false if the raw data has incomplete data header
 584    */
 585   Bool_t HasCompleteHeader() const { return fHasCompleteHeader; }
 586
 587   /** @} */
 588
 589   static const char* PulseGainPath()       { return fgkPulseGain; }
 590   static const char* PedestalPath()        { return fgkPedestal; }
 591   static const char* DeadPath()            { return fgkDead; }
 592   static const char* SampleRatePath()      { return fgkSampleRate; }
 593   static const char* AltroMapPath()        { return fgkAltroMap; }
 594   static const char* ZeroSuppressionPath() { return fgkZeroSuppression; }
 595   static const char* StripRangePath()      { return fgkStripRange; }
 596   static const char* GetPedestalShuttleID()   {return fkPedestalShuttleID;}
 597   static const char* GetGainShuttleID()       {return fkGainShuttleID;}
 598   static const char* GetConditionsShuttleID()   {return fkConditionsShuttleID;}
 599
 600 protected:
 601   /**
 602    * CTOR
 603    */
 604   AliFMDParameters();
 605   /**
 606    * CTOR
 607    */
 608   AliFMDParameters(const AliFMDParameters& o)
 609     : TNamed(o),
 610       fIsInit(o.fIsInit),
 611       fkSiDeDxMip(o.fkSiDeDxMip),
 612       fVA1MipRange(o.fVA1MipRange),
 613       fAltroChannelSize(o.fAltroChannelSize),
 614       fChannelsPerAltro(o.fChannelsPerAltro),
 615       fPedestalFactor(o.fPedestalFactor),
 616       fZSPre(o.fZSPre),
 617       fZSPost(o.fZSPost),
 618       fZSPedSubtract(o.fZSPedSubtract),
 619       fFixedPedestal(o.fFixedPedestal),
 620       fFixedPedestalWidth(o.fFixedPedestalWidth),
 621       fFixedZeroSuppression(o.fFixedZeroSuppression),
 622       fFixedSampleRate(o.fFixedSampleRate),
 623       fFixedThreshold(o.fFixedThreshold),
 624       fFixedMinStrip(o.fFixedMinStrip),
 625       fFixedMaxStrip(o.fFixedMaxStrip),
 626       fFixedPulseGain(o.fFixedPulseGain),
 627       fEdepMip(o.fEdepMip),
 628       fHasCompleteHeader(o.fHasCompleteHeader),
 629       fZeroSuppression(o.fZeroSuppression),
 630       fSampleRate(o.fSampleRate),
 631       fPedestal(o.fPedestal),
 632       fPulseGain(o.fPulseGain),
 633       fDeadMap(o.fDeadMap),
 634       fAltroMap(o.fAltroMap),
 635       fStripRange(o.fStripRange)
 636   {}
 637   /**
 638    * Assignement operator
 639    *
 640    * @return Reference to this
 641    */
 642   AliFMDParameters& operator=(const AliFMDParameters&) { return *this; }
 643   /**
 644    * DTOR
 645    */
 646   virtual ~AliFMDParameters() {}
 647   /**
 648    * Singleton instance
 649    */
 650   static AliFMDParameters* fgInstance;   // Static singleton instance
 651   /**
 652    * Check if the file <i>prefix</i><i>number</i> exists in @a path,
 653    * and write the full path to @a f.
 654    *
 655    * @param prefix  File prefix (cond, peds, gains, ...)
 656    * @param path    Path to files
 657    * @param number  Detector number (1, 2, or 3)
 658    * @param f       On return full path to file (if found)
 659    *
 660    * @return @c true if file exists and is readable, @c false otherwise
 661    */
 662   Bool_t CheckFile(const char* prefix, const char* path,
 663                    int         number, TString&    f) const;
 664   /**
 665    * Get an entry from either global AliCDBManager or passed
 666    * AliFMDPreprocessor.
 667    *
 668    * @param path  Path to CDB object.
 669    * @param pp    AliFMDPreprocessor
 670    * @param fatal If true, raise a fatal flag if we didn't get the entry.
 671    * @return AliCDBEntry if found
 672    */
 673   AliCDBEntry* GetEntry(const char* path, AliFMDPreprocessor* pp,
 674                         Bool_t fatal=kTRUE) const;
 675   /**
 676    * Initialize gains.  Try to get them from CDB
 677    *
 678    * @param pp Pre-processor if called from shuttle
 679    */
 680   void InitPulseGain(AliFMDPreprocessor* pp=0);
 681   /**
 682    * Initialize pedestals.  Try to get them from CDB
 683    *
 684    * @param pp Pre-processor if called from shuttle
 685    */
 686   void InitPedestal(AliFMDPreprocessor* pp=0);
 687   /**
 688    * Initialize dead map.  Try to get it from CDB
 689    *
 690    * @param pp Pre-processor if called from shuttle
 691    */
 692   void InitDeadMap(AliFMDPreprocessor* pp=0);
 693   /**
 694    * Initialize sample rates.  Try to get them from CDB
 695    *
 696    * @param pp Pre-processor if called from shuttle
 697    */
 698   void InitSampleRate(AliFMDPreprocessor* pp=0);
 699   /**
 700    * Initialize zero suppression thresholds.  Try to get them from CDB
 701    *
 702    * @param pp Pre-processor if called from shuttle
 703    */
 704   void InitZeroSuppression(AliFMDPreprocessor* pp=0);
 705   /**
 706    * Initialize hardware map.  Try to get it from CDB
 707    *
 708    * @param pp Pre-processor if called from shuttle
 709    */
 710   void InitAltroMap(AliFMDPreprocessor* pp=0);
 711   /**
 712    * Initialize strip range.  Try to get it from CDB
 713    *
 714    * @param pp Pre-processor if called from shuttle
 715    */
 716   void InitStripRange(AliFMDPreprocessor* pp=0);
 717
 718   Bool_t          fIsInit;                   // Whether we've been initialised
 719
 720   static const char* fgkPulseGain;           // Path to PulseGain calib object
 721   static const char* fgkPedestal;            // Path to Pedestal calib object
 722   static const char* fgkDead;                // Path to Dead calib object
 723   static const char* fgkSampleRate;          // Path to SampleRate calib object
 724   static const char* fgkAltroMap;            // Path to AltroMap calib object
 725   static const char* fgkZeroSuppression;     // Path to ZeroSuppression cal object
 726   static const char* fgkStripRange;          // Path to strip range cal object
 727   const Float_t   fkSiDeDxMip;               // MIP dE/dx in Silicon
 728   UShort_t        fVA1MipRange;              // # MIPs the pre-amp can do
 729   UShort_t        fAltroChannelSize;         // Largest # to store in 1 ADC ch.
 730   UShort_t        fChannelsPerAltro;         // Number of pre-amp. chan/adc chan.
 731   Float_t         fPedestalFactor;           // Number of pedestal widths
 732   UShort_t        fZSPre;                    // Number of pre-samples in ZS
 733   UShort_t        fZSPost;                   // Number of post-samples in ZS
 734   Bool_t          fZSPedSubtract;            // Pedestal subtraction before ZS
 735
 736   Float_t         fFixedPedestal;            // Pedestal to subtract
 737   Float_t         fFixedPedestalWidth;       // Width of pedestal
 738   UShort_t        fFixedZeroSuppression;     // Threshold for zero-suppression
 739   UShort_t        fFixedSampleRate;          // Times the ALTRO samples pre-amp.
 740   Float_t         fFixedThreshold;           // Threshold in ADC counts
 741   UShort_t        fFixedMinStrip;            // Minimum strip read-out
 742   UShort_t        fFixedMaxStrip;            // Maximum strip read-out
 743   mutable Float_t fFixedPulseGain;           //! Gain (cached)
 744   mutable Float_t fEdepMip;                  //! Cache of energy loss for a MIP
 745   Bool_t          fHasCompleteHeader;        // raw data has incomplete data header
 746
 747   static const char* fkPedestalShuttleID;    // Shuttle/preprocessor ID for pedestals
 748   static const char* fkGainShuttleID;        // Shuttle/preprocessor ID for gains
 749   static const char* fkConditionsShuttleID;  // Shuttle/preprocessor ID for conditions
 750
 751   AliFMDCalibZeroSuppression* fZeroSuppression; // Zero suppression from CDB
 752   AliFMDCalibSampleRate*      fSampleRate;      // Sample rate from CDB
 753   AliFMDCalibPedestal*        fPedestal;        // Pedestals
 754   AliFMDCalibGain*            fPulseGain;       // Pulser gain
 755   AliFMDCalibDeadMap*         fDeadMap;         // Pulser gain
 756   AliFMDAltroMapping*         fAltroMap;        // Map of hardware
 757   AliFMDCalibStripRange*      fStripRange;      // Strip range
 758
 759   ClassDef(AliFMDParameters,6) // Manager of parameters
 760 };
 761
 762 //__________________________________________________________________
 763 inline void
 764 AliFMDParameters::SetStripRange(UShort_t min, UShort_t max)
 765 {
 766   // Set fixed strip range
 767   fFixedMinStrip = min;
 768   fFixedMaxStrip = max;
 769 }
 770
 771 #endif
 772 //____________________________________________________________________
 773 //
 774 // Local Variables:
 775 //   mode: C++
 776 // End:
 777 //
 778 // EOF
 779 //
 780