[u/mrichter/AliRoot.git] / PHOS / AliPHOSTrackSegmentMakerv1.h

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

/* $Id$ */

//_________________________________________________________________________
// Implementation version 1 of algorithm class to construct PHOS track segments
// Associates EMC and PPSD clusters
// Unfolds the EMC cluster   
//                  
//*-- Author: Dmitri Peressounko (RRC Ki & SUBATECH)

// --- ROOT system ---

// --- Standard library ---

// --- AliRoot header files ---

#include "TObjArray.h"
#include "AliPHOSClusterizer.h"
#include "AliPHOSEmcRecPoint.h"
#include "AliPHOSPpsdRecPoint.h"
#include "AliPHOSTrackSegmentMaker.h"
#include "TMinuit.h" 

class  AliPHOSTrackSegmentMakerv1 : public AliPHOSTrackSegmentMaker {

public:

  AliPHOSTrackSegmentMakerv1() ;                     
  AliPHOSTrackSegmentMakerv1(const AliPHOSTrackSegmentMakerv1 & tsm) {
    // cpy ctor: no implementation yet
    // requested by the Coding Convention
    assert(0==1) ; 
  }
   
  virtual ~ AliPHOSTrackSegmentMakerv1() ; // dtor
  
  Bool_t  FindFit(AliPHOSEmcRecPoint * emcRP, int * MaxAt, Float_t * maxAtEnergy, 
		  Int_t NPar, Float_t * FitParametres) ; //Used in UnfoldClusters, calls TMinuit
  void    FillOneModule(AliPHOSRecPoint::RecPointsList * emcIn, 
			TArrayI * emcOut, 
			AliPHOSRecPoint::RecPointsList * ppsdIn, 
			TArrayI * ppsdOutUp, 
			TArrayI * ppsdOutLow, 
			Int_t &PHOSModule, 
			Int_t & emcStopedAt, 
			Int_t & ppsdStopedAt) ; // Fills temporary arrais with clusters from one module  
  Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0

  void    MakeLinks(TArrayI * EmcRecPoints, TArrayI * PpsdRecPointsUp, TArrayI * PpsdRecPointsLow, 
		    TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
  void    MakePairs(TArrayI * EmcRecPoints, 
		    TArrayI * PpsdRecPointsUp, 
		    TArrayI * PpsdRecPointsLow, 
		    TClonesArray * LinkLowArray, 
		    TClonesArray * LinkUpArray, 
		    AliPHOSTrackSegment::TrackSegmentsList * trsl) ; //Finds pairs(triplets) with smallest link
  void    MakeTrackSegments(DigitsList * DL, 
			    AliPHOSRecPoint::RecPointsList * emcl, 
			    AliPHOSRecPoint::RecPointsList * ppsdl, 
			    AliPHOSTrackSegment::TrackSegmentsList * trsl ) ; // does the job
  virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
  virtual void    SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ; 
  static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function)
  void    UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ; 
                                             // Unfolds and sorts all EMC clusters
  void  UnfoldClusters(DigitsList * DL, 
		       AliPHOSRecPoint::RecPointsList * emcIn, 
		       AliPHOSEmcRecPoint * iniEmc, 
		       Int_t Nmax, 
		       int * maxAt, 
		       Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
  virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; } 

  AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & )  {
    // assignement operator requested by coding convention
    // but not needed
    assert(0==1) ;
    return *this ; 
  }

private:

  Float_t fDelta ;     // parameter used for sorting
  Float_t fR0 ;        // Maximum distance between a EMC RecPoint and a PPSD RecPoint   
  Bool_t fUnfoldFlag ; // Directive to unfold or not the clusters in case of multiple maxima

  ClassDef( AliPHOSTrackSegmentMakerv1,1)  // Implementation version 1 of algorithm class to make PHOS track segments 

};

#endif // AliPHOSTRACKSEGMENTMAKERV1_H
Commit	Line	Data
6ad0bfa0	1	#ifndef ALIPHOSTRACKSEGMENTMAKERV1_H
6ad0bfa0	2	#define ALIPHOSTRACKSEGMENTMAKERV1_H
d15a28e7	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
6ad0bfa0	6	/* $Id$ */
6ad0bfa0	7
b2a60966	8	//_________________________________________________________________________
	9	// Implementation version 1 of algorithm class to construct PHOS track segments
	10	// Associates EMC and PPSD clusters
	11	// Unfolds the EMC cluster
	12	//
	13	//*-- Author: Dmitri Peressounko (RRC Ki & SUBATECH)
d15a28e7	14
	15	// --- ROOT system ---
	16
	17	// --- Standard library ---
	18
	19	// --- AliRoot header files ---
	20
	21	#include "TObjArray.h"
	22	#include "AliPHOSClusterizer.h"
	23	#include "AliPHOSEmcRecPoint.h"
	24	#include "AliPHOSPpsdRecPoint.h"
	25	#include "AliPHOSTrackSegmentMaker.h"
9f616d61	26	#include "TMinuit.h"
d15a28e7	27
	28	class AliPHOSTrackSegmentMakerv1 : public AliPHOSTrackSegmentMaker {
	29
	30	public:
	31
	32	AliPHOSTrackSegmentMakerv1() ;
839ffcb3	33	AliPHOSTrackSegmentMakerv1(const AliPHOSTrackSegmentMakerv1 & tsm) {
	34	// cpy ctor: no implementation yet
	35	// requested by the Coding Convention
	36	assert(0==1) ;
	37	}
	38
9f616d61	39	virtual ~ AliPHOSTrackSegmentMakerv1() ; // dtor
d15a28e7	40
	41	Bool_t FindFit(AliPHOSEmcRecPoint * emcRP, int * MaxAt, Float_t * maxAtEnergy,
	42	Int_t NPar, Float_t * FitParametres) ; //Used in UnfoldClusters, calls TMinuit
98cbd830	43	void FillOneModule(AliPHOSRecPoint::RecPointsList * emcIn,
28c3a259	44	TArrayI * emcOut,
88714635	45	AliPHOSRecPoint::RecPointsList * ppsdIn,
28c3a259	46	TArrayI * ppsdOutUp,
28c3a259	47	TArrayI * ppsdOutLow,
88714635	48	Int_t &PHOSModule,
88714635	49	Int_t & emcStopedAt,
69183710	50	Int_t & ppsdStopedAt) ; // Fills temporary arrais with clusters from one module
d15a28e7	51	Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0
d15a28e7	52
28c3a259	53	void MakeLinks(TArrayI * EmcRecPoints, TArrayI * PpsdRecPointsUp, TArrayI * PpsdRecPointsLow,
d15a28e7	54	TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
28c3a259	55	void MakePairs(TArrayI * EmcRecPoints,
	56	TArrayI * PpsdRecPointsUp,
	57	TArrayI * PpsdRecPointsLow,
88714635	58	TClonesArray * LinkLowArray,
	59	TClonesArray * LinkUpArray,
	60	AliPHOSTrackSegment::TrackSegmentsList * trsl) ; //Finds pairs(triplets) with smallest link
	61	void MakeTrackSegments(DigitsList * DL,
	62	AliPHOSRecPoint::RecPointsList * emcl,
	63	AliPHOSRecPoint::RecPointsList * ppsdl,
	64	AliPHOSTrackSegment::TrackSegmentsList * trsl ) ; // does the job
d912c4d7	65	virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
6727be7e	66	virtual void SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ;
9f616d61	67	static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function)
98cbd830	68	void UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ;
69183710	69	// Unfolds and sorts all EMC clusters
88714635	70	void UnfoldClusters(DigitsList * DL,
	71	AliPHOSRecPoint::RecPointsList * emcIn,
	72	AliPHOSEmcRecPoint * iniEmc,
	73	Int_t Nmax,
	74	int * maxAt,
98cbd830	75	Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
d912c4d7	76	virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; }
d15a28e7	77
2f04ed65	78	AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & ) {
839ffcb3	79	// assignement operator requested by coding convention
	80	// but not needed
	81	assert(0==1) ;
	82	return *this ;
	83	}
	84
d15a28e7	85	private:
d15a28e7	86
b2a60966	87	Float_t fDelta ; // parameter used for sorting
b2a60966	88	Float_t fR0 ; // Maximum distance between a EMC RecPoint and a PPSD RecPoint
b2a60966	89	Bool_t fUnfoldFlag ; // Directive to unfold or not the clusters in case of multiple maxima
d15a28e7	90
b2a60966	91	ClassDef( AliPHOSTrackSegmentMakerv1,1) // Implementation version 1 of algorithm class to make PHOS track segments
d15a28e7	92
	93	};
	94
6ad0bfa0	95	#endif // AliPHOSTRACKSEGMENTMAKERV1_H