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

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

/* $Id$ */
/* $Id$ */


////////////////////////////////////////////////
//  AliDetector Class for MUON subsystem      //
////////////////////////////////////////////////
#include "AliDetector.h"

class AliMUONChamber;
class AliMUONLocalTrigger;
class AliMUONGlobalTrigger;
class AliMUONTriggerCircuit;
class AliMUONTriggerDecision;
class AliSegmentation;
class AliMUONResponse;
class AliMUONMerger;
class AliMUONHit;
class AliMUONPadHit;
class AliMUONRawCluster;
class AliMUONClusterFinderVS;
class AliMUONReconstHit;
class AliMUONMerger;

class TVector;
#include "TObjArray.h"
class TFile;
class TTree;


class AliMUON : public  AliDetector {
 public:
    AliMUON();
    AliMUON(const char *name, const char *title);
    AliMUON(const AliMUON& rMUON);
    virtual       ~AliMUON();
    virtual void   AddHit(Int_t track , Int_t *vol, Float_t *hits);
    virtual void   AddHit(Int_t fIshunt, Int_t track, Int_t iChamber, 
			  Int_t idpart, Float_t X, Float_t Y, Float_t Z, 
			  Float_t tof, Float_t momentum, Float_t theta, 
			  Float_t phi, Float_t length, Float_t destep);
    virtual void   AddPadHit(Int_t* clhits); // To be removed !
    virtual void   AddDigits(Int_t id, Int_t* tracks, Int_t* charges,
			     Int_t* digits);
    virtual void   AddRawCluster(Int_t id, const AliMUONRawCluster& clust);
    virtual void   BuildGeometry();
    void           AddGlobalTrigger(Int_t *singlePlus, Int_t *singleMinus,
				    Int_t *singleUndef, Int_t *pairUnlike, 
				    Int_t *pairLike);
    void           AddLocalTrigger(Int_t* ltrigger);
    Int_t          DistancetoPrimitive(Int_t px, Int_t py);
    virtual Int_t  IsVersion() const {return 0;}
    TClonesArray  *PadHits() {return fPadHits;}
    TClonesArray  *LocalTrigger() {return fLocalTrigger;}
    TClonesArray  *GlobalTrigger() {return fGlobalTrigger;}
    virtual void   MakeBranch(Option_t *opt=" ");
    virtual void   MakeBranchInTreeD(TTree *treeD, const char *file=0);
    void           SetTreeAddress();
    virtual void   ResetHits();
    virtual void   ResetDigits();
    virtual void   ResetTrigger();
    virtual void   ResetRawClusters();
    // Cluster Finding
    virtual void   Digits2Reco();
    virtual void   FindClusters();
    // Digitisation 
    virtual void   SDigits2Digits();      
// Configuration Methods (per station id)
//
// Set Chamber Segmentation Parameters
// id refers to the station and isec to the cathode plane
// Set Z values for all chambers
    virtual void SetChambersZ(const Float_t *Z);
    virtual void SetChambersZToDefault(void);
    virtual void SetPadSize(Int_t id, Int_t isec, Float_t p1, Float_t p2);
// Set Signal Generation Parameters
    virtual void   SetSigmaIntegration(Int_t id, Float_t p1);
    virtual void   SetChargeSlope(Int_t id, Float_t p1);
    virtual void   SetChargeSpread(Int_t id, Float_t p1, Float_t p2);
    virtual void   SetMaxAdc(Int_t id, Int_t p1);
// Set Segmentation and Response Model
    virtual void   SetSegmentationModel(Int_t id, Int_t isec,
					AliSegmentation *segmentation);
    virtual void   SetResponseModel(Int_t id, AliMUONResponse *response);
    virtual void   SetNsec(Int_t id, Int_t nsec);
// Set Reconstruction Model
    virtual void   SetReconstructionModel(Int_t id, AliMUONClusterFinderVS *reconstruction);
// Set Merger/Digitizer
    virtual void   SetMerger(AliMUONMerger* merger);
    virtual AliMUONMerger* Merger();
    
// Set Stepping Parameters
    virtual void   SetMaxStepGas(Float_t p1);
    virtual void   SetMaxStepAlu(Float_t p1);
    virtual void   SetMaxDestepGas(Float_t p1);
    virtual void   SetMaxDestepAlu(Float_t p1);
    virtual void   SetAcceptance(Bool_t acc=0, Float_t angmin=2, Float_t angmax=9);
// Response Simulation
    virtual void   MakePadHits(Float_t xhit,Float_t yhit, Float_t zhit,
			       Float_t eloss, Float_t tof, Int_t id);
// get Trigger answer
    void   Trigger(Int_t nev);
// Return reference to Chamber #id
    virtual AliMUONChamber& Chamber(Int_t id)
	{return *((AliMUONChamber *) (*fChambers)[id]);}
// Return reference to Circuit #id
    virtual AliMUONTriggerCircuit& TriggerCircuit(Int_t id)
      {return *((AliMUONTriggerCircuit *) (*fTriggerCircuits)[id]);}
// Retrieve pad hits for a given Hit
    virtual AliMUONPadHit* FirstPad(AliMUONHit *hit, TClonesArray *padHits);
    virtual AliMUONPadHit* NextPad(TClonesArray *padHits);
// Return pointers to digits
    TObjArray            *Dchambers() {return fDchambers;}
    Int_t                *Ndch() {return fNdch;}
    virtual TClonesArray *DigitsAddress(Int_t id)
	{return ((TClonesArray *) (*fDchambers)[id]);}
// Return pointers to reconstructed clusters
    TObjArray            *RawClusters() {return fRawClusters;}
    Int_t                *Nrawch() {return fNrawch;} 
    virtual TClonesArray *RawClustAddress(Int_t id) 
	{return ((TClonesArray *) (*fRawClusters)[id]);}

    AliMUONRawCluster    *RawCluster(Int_t ichamber, Int_t icathod,
				     Int_t icluster);
// Copy Operator
    AliMUON& operator = (const AliMUON& rhs);
    
	    
 protected:
    Int_t                 fNCh;                // Number of chambers   
    Int_t                 fNTrackingCh;        // Number of tracking chambers
    TObjArray            *fChambers;           //! List of Tracking Chambers
    TObjArray            *fTriggerCircuits;    //! List of Trigger Circuits
    Int_t                 fNPadHits;           // Number of pad hits
    TClonesArray         *fPadHits;            //! List of pad hits
    TObjArray            *fDchambers;          //! List of digits
    Int_t                *fNdch;               // [fNCh] Number of digits per chamber
    TObjArray            *fRawClusters;        //! List of raw clusters
    Int_t                *fNrawch;             // [fNTrackingCh] Number of raw clusters  per chamber
    Int_t                 fNLocalTrigger;      // Number of Local Trigger 
    TClonesArray         *fLocalTrigger;       //! List of Local Trigger      
    Int_t                 fNGlobalTrigger;     // Number of Global Trigger
    TClonesArray         *fGlobalTrigger;      //! List of Global Trigger  
    
//
    Bool_t   fAccCut;          //Transport acceptance cut
    Float_t  fAccMin;          //Minimum acceptance cut used during transport
    Float_t  fAccMax;          //Minimum acceptance cut used during transport
//  

//  Stepping Parameters
   Float_t fMaxStepGas;      // Maximum step size inside the chamber gas
   Float_t fMaxStepAlu;      // Maximum step size inside the chamber aluminum
   Float_t fMaxDestepGas;    // Maximum relative energy loss in gas
   Float_t fMaxDestepAlu;    // Maximum relative energy loss in aluminum


// Pad Iterator
   Int_t fMaxIterPad;        // Maximum pad index
   Int_t fCurIterPad;        // Current pad index
// Background eent for event mixing
   Text_t *fFileName;           // ! File with background hits
   AliMUONMerger *fMerger;   // ! pointer to merger

   ClassDef(AliMUON,3)  // MUON Detector base class
};
#endif
Commit	Line	Data
a9e2aefa	1	#ifndef ALIMUON_H
a9e2aefa	2	#define ALIMUON_H
3da30618	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
88cb7938	7	/* $Id$ */
88cb7938	8
3da30618	9
fe4da5cc	10	////////////////////////////////////////////////
5a35ac13	11	// AliDetector Class for MUON subsystem //
fe4da5cc	12	////////////////////////////////////////////////
fe4da5cc	13	#include "AliDetector.h"
a9e2aefa	14
a9e2aefa	15	class AliMUONChamber;
	16	class AliMUONLocalTrigger;
	17	class AliMUONGlobalTrigger;
	18	class AliMUONTriggerCircuit;
	19	class AliMUONTriggerDecision;
a30a000f	20	class AliSegmentation;
a9e2aefa	21	class AliMUONResponse;
375da359	22	class AliMUONMerger;
a9e2aefa	23	class AliMUONHit;
a9e2aefa	24	class AliMUONPadHit;
a897a37a	25	class AliMUONRawCluster;
30aaba74	26	class AliMUONClusterFinderVS;
a9e2aefa	27	class AliMUONReconstHit;
1cae9436	28	class AliMUONMerger;
1cae9436	29
a9e2aefa	30	class TVector;
94de3818	31	#include "TObjArray.h"
a9e2aefa	32	class TFile;
a9e2aefa	33	class TTree;
fe4da5cc	34
fe4da5cc	35
fe4da5cc	36	class AliMUON : public AliDetector {
	37	public:
	38	AliMUON();
	39	AliMUON(const char name, const char title);
a9e2aefa	40	AliMUON(const AliMUON& rMUON);
fe4da5cc	41	virtual ~AliMUON();
a9e2aefa	42	virtual void AddHit(Int_t track , Int_t vol, Float_t hits);
5a35ac13	43	virtual void AddHit(Int_t fIshunt, Int_t track, Int_t iChamber,
	44	Int_t idpart, Float_t X, Float_t Y, Float_t Z,
	45	Float_t tof, Float_t momentum, Float_t theta,
	46	Float_t phi, Float_t length, Float_t destep);
	47	virtual void AddPadHit(Int_t* clhits); // To be removed !
a9e2aefa	48	virtual void AddDigits(Int_t id, Int_t* tracks, Int_t* charges,
	49	Int_t* digits);
	50	virtual void AddRawCluster(Int_t id, const AliMUONRawCluster& clust);
fe4da5cc	51	virtual void BuildGeometry();
a9e2aefa	52	void AddGlobalTrigger(Int_t singlePlus, Int_t singleMinus,
	53	Int_t singleUndef, Int_t pairUnlike,
	54	Int_t *pairLike);
	55	void AddLocalTrigger(Int_t* ltrigger);
fe4da5cc	56	Int_t DistancetoPrimitive(Int_t px, Int_t py);
a9e2aefa	57	virtual Int_t IsVersion() const {return 0;}
	58	TClonesArray *PadHits() {return fPadHits;}
	59	TClonesArray *LocalTrigger() {return fLocalTrigger;}
	60	TClonesArray *GlobalTrigger() {return fGlobalTrigger;}
88cb7938	61	virtual void MakeBranch(Option_t *opt=" ");
cd4df77b	62	virtual void MakeBranchInTreeD(TTree treeD, const char file=0);
fe4da5cc	63	void SetTreeAddress();
	64	virtual void ResetHits();
	65	virtual void ResetDigits();
a9e2aefa	66	virtual void ResetTrigger();
a897a37a	67	virtual void ResetRawClusters();
a9e2aefa	68	// Cluster Finding
fe311ef1	69	virtual void Digits2Reco();
fe311ef1	70	virtual void FindClusters();
a9e2aefa	71	// Digitisation
2ab0c725	72	virtual void SDigits2Digits();
fe4da5cc	73	// Configuration Methods (per station id)
	74	//
	75	// Set Chamber Segmentation Parameters
a9e2aefa	76	// id refers to the station and isec to the cathode plane
	77	// Set Z values for all chambers
	78	virtual void SetChambersZ(const Float_t *Z);
	79	virtual void SetChambersZToDefault(void);
	80	virtual void SetPadSize(Int_t id, Int_t isec, Float_t p1, Float_t p2);
fe4da5cc	81	// Set Signal Generation Parameters
a897a37a	82	virtual void SetSigmaIntegration(Int_t id, Float_t p1);
	83	virtual void SetChargeSlope(Int_t id, Float_t p1);
	84	virtual void SetChargeSpread(Int_t id, Float_t p1, Float_t p2);
3b3ee6f4	85	virtual void SetMaxAdc(Int_t id, Int_t p1);
fe4da5cc	86	// Set Segmentation and Response Model
a9e2aefa	87	virtual void SetSegmentationModel(Int_t id, Int_t isec,
a30a000f	88	AliSegmentation *segmentation);
a9e2aefa	89	virtual void SetResponseModel(Int_t id, AliMUONResponse *response);
fe4da5cc	90	virtual void SetNsec(Int_t id, Int_t nsec);
a897a37a	91	// Set Reconstruction Model
30aaba74	92	virtual void SetReconstructionModel(Int_t id, AliMUONClusterFinderVS *reconstruction);
1cae9436	93	// Set Merger/Digitizer
	94	virtual void SetMerger(AliMUONMerger* merger);
	95	virtual AliMUONMerger* Merger();
	96
fe4da5cc	97	// Set Stepping Parameters
a897a37a	98	virtual void SetMaxStepGas(Float_t p1);
	99	virtual void SetMaxStepAlu(Float_t p1);
	100	virtual void SetMaxDestepGas(Float_t p1);
	101	virtual void SetMaxDestepAlu(Float_t p1);
5c1f55c5	102	virtual void SetAcceptance(Bool_t acc=0, Float_t angmin=2, Float_t angmax=9);
fe4da5cc	103	// Response Simulation
802a864d	104	virtual void MakePadHits(Float_t xhit,Float_t yhit, Float_t zhit,
a9e2aefa	105	Float_t eloss, Float_t tof, Int_t id);
	106	// get Trigger answer
	107	void Trigger(Int_t nev);
fe4da5cc	108	// Return reference to Chamber #id
a9e2aefa	109	virtual AliMUONChamber& Chamber(Int_t id)
	110	{return ((AliMUONChamber ) (*fChambers)[id]);}
	111	// Return reference to Circuit #id
	112	virtual AliMUONTriggerCircuit& TriggerCircuit(Int_t id)
	113	{return ((AliMUONTriggerCircuit ) (*fTriggerCircuits)[id]);}
fe4da5cc	114	// Retrieve pad hits for a given Hit
a9e2aefa	115	virtual AliMUONPadHit* FirstPad(AliMUONHit hit, TClonesArray padHits);
	116	virtual AliMUONPadHit* NextPad(TClonesArray *padHits);
	117	// Return pointers to digits
fe4da5cc	118	TObjArray *Dchambers() {return fDchambers;}
fe4da5cc	119	Int_t *Ndch() {return fNdch;}
a9e2aefa	120	virtual TClonesArray *DigitsAddress(Int_t id)
a9e2aefa	121	{return ((TClonesArray ) (fDchambers)[id]);}
fe4da5cc	122	// Return pointers to reconstructed clusters
a897a37a	123	TObjArray *RawClusters() {return fRawClusters;}
a9e2aefa	124	Int_t *Nrawch() {return fNrawch;}
	125	virtual TClonesArray *RawClustAddress(Int_t id)
	126	{return ((TClonesArray ) (fRawClusters)[id]);}
	127
	128	AliMUONRawCluster *RawCluster(Int_t ichamber, Int_t icathod,
	129	Int_t icluster);
	130	// Copy Operator
	131	AliMUON& operator = (const AliMUON& rhs);
a897a37a	132
a9e2aefa	133
fe4da5cc	134	protected:
1bd26093	135	Int_t fNCh; // Number of chambers
1bd26093	136	Int_t fNTrackingCh; // Number of tracking chambers
88cb7938	137	TObjArray *fChambers; //! List of Tracking Chambers
88cb7938	138	TObjArray *fTriggerCircuits; //! List of Trigger Circuits
a9e2aefa	139	Int_t fNPadHits; // Number of pad hits
88cb7938	140	TClonesArray *fPadHits; //! List of pad hits
88cb7938	141	TObjArray *fDchambers; //! List of digits
1bd26093	142	Int_t *fNdch; // [fNCh] Number of digits per chamber
88cb7938	143	TObjArray *fRawClusters; //! List of raw clusters
1bd26093	144	Int_t *fNrawch; // [fNTrackingCh] Number of raw clusters per chamber
a9e2aefa	145	Int_t fNLocalTrigger; // Number of Local Trigger
88cb7938	146	TClonesArray *fLocalTrigger; //! List of Local Trigger
a9e2aefa	147	Int_t fNGlobalTrigger; // Number of Global Trigger
88cb7938	148	TClonesArray *fGlobalTrigger; //! List of Global Trigger
1cae9436	149
fe4da5cc	150	//
	151	Bool_t fAccCut; //Transport acceptance cut
	152	Float_t fAccMin; //Minimum acceptance cut used during transport
	153	Float_t fAccMax; //Minimum acceptance cut used during transport
	154	//
	155
	156	// Stepping Parameters
	157	Float_t fMaxStepGas; // Maximum step size inside the chamber gas
	158	Float_t fMaxStepAlu; // Maximum step size inside the chamber aluminum
	159	Float_t fMaxDestepGas; // Maximum relative energy loss in gas
	160	Float_t fMaxDestepAlu; // Maximum relative energy loss in aluminum
a897a37a	161
fe4da5cc	162
a897a37a	163
a9e2aefa	164	// Pad Iterator
	165	Int_t fMaxIterPad; // Maximum pad index
	166	Int_t fCurIterPad; // Current pad index
	167	// Background eent for event mixing
30f212b7	168	Text_t *fFileName; // ! File with background hits
375da359	169	AliMUONMerger *fMerger; // ! pointer to merger
a897a37a	170
88cb7938	171	ClassDef(AliMUON,3) // MUON Detector base class
a897a37a	172	};
fe4da5cc	173	#endif
fe4da5cc	174