MUON/AliMUONClusterInput.h

   1 #ifndef ALIMUONCLUSTERINPUT_H
   2 #define ALIMUONCLUSTERINPUT_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id $*/
   7 // Revision of includes 07/05/2004
   8
   9 #include <TObject.h>
  10 #include <TClonesArray.h> // needed for inline function Digit
  11 #include "AliMUONGeometrySegmentation.h"
  12
  13 class TMinuit;
  14
  15 class AliMUONDigit;
  16 class AliMUONRawCluster;
  17 class AliSegmentation;
  18 class AliMUONMathieson;
  19
  20 class AliMUONClusterInput : public TObject
  21 {
  22  public:
  23     virtual ~AliMUONClusterInput();
  24     static AliMUONClusterInput* Instance();
  25 //  Configuration
  26     void SetDigits(Int_t chamber, TClonesArray* dig1, TClonesArray* dig2);
  27     void SetDigits(Int_t chamber, TClonesArray* dig);
  28     void SetDigits(Int_t chamber, Int_t idDE, TClonesArray* dig1, TClonesArray* dig2);
  29     void SetDigits(Int_t chamber, Int_t idDE, TClonesArray* dig);
  30     void SetCluster(AliMUONRawCluster* cluster);
  31 // Access functions
  32     Int_t Chamber() const {return fChamber;}
  33     AliMUONDigit* Digit(Int_t cath, Int_t i) const {return (AliMUONDigit*) (fDigits[cath]->UncheckedAt(i));}
  34     TClonesArray* Digits(Int_t cath) const {return fDigits[cath];}
  35     Int_t NDigits(Int_t cath) const {return fNDigits[cath];}
  36     AliSegmentation* Segmentation(Int_t cath) const {return fSegmentation[cath];}
  37     AliMUONGeometrySegmentation* Segmentation2(Int_t cath) const {return fSegmentation2[cath];}
  38
  39     AliMUONMathieson* Mathieson() const {return fgMathieson;}
  40     Float_t ChargeCorrel() const {return fChargeCorrel;}
  41     Int_t DetElemId() const {return fDetElemId;}
  42     Int_t GetSegmentationType() const {return fSegmentationType;}
  43
  44 // Fitting
  45     TMinuit*      Fitter() const {return fgMinuit;}
  46 // Current cluster information
  47     Float_t       TotalCharge(Int_t cath) const {return fChargeTot[cath];}
  48     Float_t       Charge(Int_t dig, Int_t cath) const {return fCharge[dig][cath];}
  49     Int_t         Ix(Int_t dig, Int_t cath) const {return fix[dig][cath];}
  50     Int_t         Iy(Int_t dig, Int_t cath) const {return fiy[dig][cath];}
  51     Int_t         Nmul(Int_t cath) const {return fNmul[cath];}
  52 //  Helpers for Fit
  53     Float_t DiscrChargeS1(Int_t i,Double_t *par);
  54     Float_t DiscrChargeCombiS1(Int_t i,Double_t *par, Int_t cath);
  55     Float_t DiscrChargeS2(Int_t i,Double_t *par);
  56     Float_t DiscrChargeCombiS2(Int_t i,Double_t *par, Int_t cath);
  57 //
  58  protected:
  59     AliMUONClusterInput();
  60     AliMUONClusterInput(const AliMUONClusterInput& clusterInput);
  61     AliMUONClusterInput & operator = (const AliMUONClusterInput& rhs);
  62  private:
  63     static AliMUONClusterInput* fgClusterInput; // ! singleton instance
  64     static AliMUONMathieson*    fgMathieson;     // ! Mathieson
  65
  66     // Digits
  67     TClonesArray*        fDigits[2];       // ! Array of pointers to digits
  68     Int_t                fNDigits[2];      // ! Number of digits
  69     AliSegmentation*     fSegmentation[2]; // ! Segmentation per cathode
  70     AliMUONGeometrySegmentation*     fSegmentation2[2]; // ! Geometry Segmentation per cathode
  71
  72     Int_t                fNseg;            // ! number of cathode planes
  73     Int_t                fChamber;         // ! Current chamber number
  74
  75     // Current cluster
  76     AliMUONRawCluster*   fCluster;         // ! current cluster
  77     Int_t                fNmul[2];         // ! current cluster multiplicity
  78     // Digits contribuing to current cluster
  79     Int_t                fix[500][2];      // ! List of x-positions for current cluster
  80     Int_t                fiy[500][2];      // ! List of y-positions for current cluster
  81     Float_t              fCharge[500][2];  // ! List of charges for current cluster
  82     Int_t                fChargeTot[2];    // ! Total charge
  83     Float_t              fQtot[2];         // ! Total charge
  84     Float_t              fZ;               // ! Current z-position
  85     Float_t              fChargeCorrel;    // ! charge correlation
  86     Int_t                fSegmentationType;// ! Segmentation type old or new
  87     Int_t                fDetElemId;       // ! detection elt id
  88
  89     // Fitter
  90     static TMinuit*      fgMinuit;          // ! Fitter
  91     ClassDef(AliMUONClusterInput, 0)        // Global data service for hit reconstruction
  92 };
  93 #endif
  94