[u/mrichter/AliRoot.git] / MUON / AliMUONClusterInput.h

#ifndef ALIMUONCLUSTERINPUT_H
#define ALIMUONCLUSTERINPUT_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id $*/
// Revision of includes 07/04/2006

/// \ingroup rec
/// \class AliMUONClusterInput
/// \brief Global data service for hit reconstruction

#include <TObject.h>
#include <TClonesArray.h> // needed for inline function Digit

class AliMUONDigit;
class AliMUONRawCluster;
class AliMUONMathieson;
class AliMUONGeometryTransformer;
class AliMUONSegmentation;
class AliMUONGeometrySegmentation;

class TMinuit;

class AliMUONClusterInput : public TObject 
{
 public:
    virtual ~AliMUONClusterInput();
    static AliMUONClusterInput* Instance();
//  Configuration
    void SetDigits(Int_t chamber, Int_t idDE, TClonesArray* dig1, TClonesArray* dig2);
    void SetDigits(Int_t chamber, Int_t idDE, TClonesArray* dig);
    void SetCluster(AliMUONRawCluster* cluster);
// Access functions
    /// Return the current chamber number
    Int_t Chamber() const {return fChamber;}
    /// Return i-th digit for given cath
    AliMUONDigit* Digit(Int_t cath, Int_t i) const {return (AliMUONDigit*) (fDigits[cath]->UncheckedAt(i));}
    /// Return the array of digits for given cathod 
    TClonesArray* Digits(Int_t cath) const {return fDigits[cath];}
    /// Return number of digits for given cathod 
    Int_t NDigits(Int_t cath) const {return fNDigits[cath];}
    /// Return geometry segmentation for given cathod 
    AliMUONGeometrySegmentation* Segmentation2(Int_t cath) const {return fSegmentation2[cath];}

    /// Return Mathieson    
    AliMUONMathieson* Mathieson() const {return fgMathieson;}    
    /// Return charge correlation
    Float_t ChargeCorrel() const {return fChargeCorrel;}    
    /// Return detection elt id
    Int_t DetElemId() const {return fDetElemId;}

// Fitting    
    /// Return the fitter
    TMinuit*      Fitter() const {return fgMinuit;}
// Current cluster information    
    /// Return the total charge for given cathod 
    Float_t       TotalCharge(Int_t cath) const {return fChargeTot[cath];}
    /// Return the charge for the given cluster and cathod
    Float_t       Charge(Int_t dig, Int_t cath) const {return fCharge[dig][cath];}
    /// Return the x-position for the given cluster and cathod
    Int_t         Ix(Int_t dig, Int_t cath) const {return fix[dig][cath];}
    /// Return the y-position for the given cluster and cathod
    Int_t         Iy(Int_t dig, Int_t cath) const {return fiy[dig][cath];}
    /// Return the cluster multiplicity for given cathod 
    Int_t         Nmul(Int_t cath) const {return fNmul[cath];}

//  Helpers for Fit     
    Float_t DiscrChargeS1(Int_t i,Double_t *par);
    Float_t DiscrChargeCombiS1(Int_t i,Double_t *par, Int_t cath);
    Float_t DiscrChargeS2(Int_t i,Double_t *par);
    Float_t DiscrChargeCombiS2(Int_t i,Double_t *par, Int_t cath);
// 
 private:
    /// Not implemented
    AliMUONClusterInput();
    /// Not implemented
    AliMUONClusterInput(const AliMUONClusterInput& clusterInput);
    /// Not implemented
    AliMUONClusterInput & operator = (const AliMUONClusterInput& rhs);

    static AliMUONClusterInput*  fgClusterInput;  //!< singleton instance
    static AliMUONMathieson*     fgMathieson;     //!< Mathieson
    static AliMUONGeometryTransformer*  fgTransformer;  //!< Geometry transformer
    static AliMUONSegmentation*  fgSegmentation;  //!< Segmentation

    // Digits
    TClonesArray*        fDigits[2];       //!< Array of pointers to digits
    Int_t                fNDigits[2];      //!< Number of digits
    AliMUONGeometrySegmentation*     fSegmentation2[2]; //!< Geometry Segmentation per cathode

    Int_t                fNseg;            //!< number of cathode planes
    Int_t                fChamber;         //!< Current chamber number
    
    // Current cluster
    AliMUONRawCluster*   fCluster;         //!< current cluster
    Int_t                fNmul[2];         //!< current cluster multiplicity
    // Digits contribuing to current cluster
    Int_t                fix[500][2];      //!< List of x-positions for current cluster
    Int_t                fiy[500][2];      //!< List of y-positions for current cluster
    Float_t              fCharge[500][2];  //!< List of charges for current cluster
    Int_t                fChargeTot[2];    //!< Total charge
    Float_t              fQtot[2];         //!< Total charge
    Float_t              fZ;               //!< Current z-position
    Float_t              fChargeCorrel;    //!< charge correlation 
    Int_t                fDetElemId;       //!< detection elt id    

    // Fitter
    static TMinuit*      fgMinuit;          //!< Fitter
    ClassDef(AliMUONClusterInput, 0)        // Global data service for hit reconstruction
};
#endif
Commit	Line	Data
	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/04/2006
	8
	9	/// \ingroup rec
	10	/// \class AliMUONClusterInput
	11	/// \brief Global data service for hit reconstruction
	12
	13	#include <TObject.h>
	14	#include <TClonesArray.h> // needed for inline function Digit
	15
	16	class AliMUONDigit;
	17	class AliMUONRawCluster;
	18	class AliMUONMathieson;
	19	class AliMUONGeometryTransformer;
	20	class AliMUONSegmentation;
	21	class AliMUONGeometrySegmentation;
	22
	23	class TMinuit;
	24
	25	class AliMUONClusterInput : public TObject
	26	{
	27	public:
	28	virtual ~AliMUONClusterInput();
	29	static AliMUONClusterInput* Instance();
	30	// Configuration
	31	void SetDigits(Int_t chamber, Int_t idDE, TClonesArray* dig1, TClonesArray* dig2);
	32	void SetDigits(Int_t chamber, Int_t idDE, TClonesArray* dig);
	33	void SetCluster(AliMUONRawCluster* cluster);
	34	// Access functions
	35	/// Return the current chamber number
	36	Int_t Chamber() const {return fChamber;}
	37	/// Return i-th digit for given cath
	38	AliMUONDigit* Digit(Int_t cath, Int_t i) const {return (AliMUONDigit*) (fDigits[cath]->UncheckedAt(i));}
	39	/// Return the array of digits for given cathod
	40	TClonesArray* Digits(Int_t cath) const {return fDigits[cath];}
	41	/// Return number of digits for given cathod
	42	Int_t NDigits(Int_t cath) const {return fNDigits[cath];}
	43	/// Return geometry segmentation for given cathod
	44	AliMUONGeometrySegmentation* Segmentation2(Int_t cath) const {return fSegmentation2[cath];}
	45
	46	/// Return Mathieson
	47	AliMUONMathieson* Mathieson() const {return fgMathieson;}
	48	/// Return charge correlation
	49	Float_t ChargeCorrel() const {return fChargeCorrel;}
	50	/// Return detection elt id
	51	Int_t DetElemId() const {return fDetElemId;}
	52
	53	// Fitting
	54	/// Return the fitter
	55	TMinuit* Fitter() const {return fgMinuit;}
	56	// Current cluster information
	57	/// Return the total charge for given cathod
	58	Float_t TotalCharge(Int_t cath) const {return fChargeTot[cath];}
	59	/// Return the charge for the given cluster and cathod
	60	Float_t Charge(Int_t dig, Int_t cath) const {return fCharge[dig][cath];}
	61	/// Return the x-position for the given cluster and cathod
	62	Int_t Ix(Int_t dig, Int_t cath) const {return fix[dig][cath];}
	63	/// Return the y-position for the given cluster and cathod
	64	Int_t Iy(Int_t dig, Int_t cath) const {return fiy[dig][cath];}
	65	/// Return the cluster multiplicity for given cathod
	66	Int_t Nmul(Int_t cath) const {return fNmul[cath];}
	67
	68	// Helpers for Fit
	69	Float_t DiscrChargeS1(Int_t i,Double_t *par);
	70	Float_t DiscrChargeCombiS1(Int_t i,Double_t *par, Int_t cath);
	71	Float_t DiscrChargeS2(Int_t i,Double_t *par);
	72	Float_t DiscrChargeCombiS2(Int_t i,Double_t *par, Int_t cath);
	73	//
	74	private:
	75	/// Not implemented
	76	AliMUONClusterInput();
	77	/// Not implemented
	78	AliMUONClusterInput(const AliMUONClusterInput& clusterInput);
	79	/// Not implemented
	80	AliMUONClusterInput & operator = (const AliMUONClusterInput& rhs);
	81
	82	static AliMUONClusterInput* fgClusterInput; //!< singleton instance
	83	static AliMUONMathieson* fgMathieson; //!< Mathieson
	84	static AliMUONGeometryTransformer* fgTransformer; //!< Geometry transformer
	85	static AliMUONSegmentation* fgSegmentation; //!< Segmentation
	86
	87	// Digits
	88	TClonesArray* fDigits[2]; //!< Array of pointers to digits
	89	Int_t fNDigits[2]; //!< Number of digits
	90	AliMUONGeometrySegmentation* fSegmentation2[2]; //!< Geometry Segmentation per cathode
	91
	92	Int_t fNseg; //!< number of cathode planes
	93	Int_t fChamber; //!< Current chamber number
	94
	95	// Current cluster
	96	AliMUONRawCluster* fCluster; //!< current cluster
	97	Int_t fNmul[2]; //!< current cluster multiplicity
	98	// Digits contribuing to current cluster
	99	Int_t fix[500][2]; //!< List of x-positions for current cluster
	100	Int_t fiy[500][2]; //!< List of y-positions for current cluster
	101	Float_t fCharge[500][2]; //!< List of charges for current cluster
	102	Int_t fChargeTot[2]; //!< Total charge
	103	Float_t fQtot[2]; //!< Total charge
	104	Float_t fZ; //!< Current z-position
	105	Float_t fChargeCorrel; //!< charge correlation
	106	Int_t fDetElemId; //!< detection elt id
	107
	108	// Fitter
	109	static TMinuit* fgMinuit; //!< Fitter
	110	ClassDef(AliMUONClusterInput, 0) // Global data service for hit reconstruction
	111	};
	112	#endif
	113