[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 EMC+PPSD
  void    FillOneModule(AliPHOSRecPoint::RecPointsList * emcIn, 
			TArrayI * emcOut, 
			AliPHOSRecPoint::RecPointsList * cpvIn, 
			TArrayI * cpvOut, 
			Int_t & PHOSModule, 
			Int_t & emcStopedAt, 
			Int_t & cpvStopedAt) ;  // Fills temporary arrais with clusters from one module EMC+CPV
  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
  void    MakeTrackSegmentsCPV(DigitsList * DL, 
			    AliPHOSRecPoint::RecPointsList * emcl, 
			    AliPHOSRecPoint::RecPointsList * ppsdl ) ; // just unfold EMC and CPV clusters
  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 & rvalue)  {
    // assignement operator requested by coding convention
    // but not needed
    assert(0==1) ;
    return *this ; 
  }

private:

  Float_t fDelta ;     // parameter used for sorting
  TMinuit * fMinuit ;  // Minuit object needed by cluster unfolding
  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,
0336aa72	50	Int_t & ppsdStopedAt) ; // Fills temporary arrais with clusters from one module EMC+PPSD
	51	void FillOneModule(AliPHOSRecPoint::RecPointsList * emcIn,
	52	TArrayI * emcOut,
	53	AliPHOSRecPoint::RecPointsList * cpvIn,
	54	TArrayI * cpvOut,
	55	Int_t & PHOSModule,
	56	Int_t & emcStopedAt,
	57	Int_t & cpvStopedAt) ; // Fills temporary arrais with clusters from one module EMC+CPV
d15a28e7	58	Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0
d15a28e7	59
28c3a259	60	void MakeLinks(TArrayI * EmcRecPoints, TArrayI * PpsdRecPointsUp, TArrayI * PpsdRecPointsLow,
d15a28e7	61	TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
28c3a259	62	void MakePairs(TArrayI * EmcRecPoints,
	63	TArrayI * PpsdRecPointsUp,
	64	TArrayI * PpsdRecPointsLow,
88714635	65	TClonesArray * LinkLowArray,
	66	TClonesArray * LinkUpArray,
	67	AliPHOSTrackSegment::TrackSegmentsList * trsl) ; //Finds pairs(triplets) with smallest link
	68	void MakeTrackSegments(DigitsList * DL,
	69	AliPHOSRecPoint::RecPointsList * emcl,
	70	AliPHOSRecPoint::RecPointsList * ppsdl,
	71	AliPHOSTrackSegment::TrackSegmentsList * trsl ) ; // does the job
0336aa72	72	void MakeTrackSegmentsCPV(DigitsList * DL,
	73	AliPHOSRecPoint::RecPointsList * emcl,
	74	AliPHOSRecPoint::RecPointsList * ppsdl ) ; // just unfold EMC and CPV clusters
d912c4d7	75	virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
6727be7e	76	virtual void SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ;
9f616d61	77	static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function)
98cbd830	78	void UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ;
0336aa72	79	// Unfolds and sorts all EMC clusters
88714635	80	void UnfoldClusters(DigitsList * DL,
	81	AliPHOSRecPoint::RecPointsList * emcIn,
	82	AliPHOSEmcRecPoint * iniEmc,
	83	Int_t Nmax,
	84	int * maxAt,
98cbd830	85	Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
d912c4d7	86	virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; }
d15a28e7	87
839ffcb3	88	AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & rvalue) {
	89	// assignement operator requested by coding convention
	90	// but not needed
	91	assert(0==1) ;
	92	return *this ;
	93	}
	94
d15a28e7	95	private:
d15a28e7	96
b2a60966	97	Float_t fDelta ; // parameter used for sorting
	98	TMinuit * fMinuit ; // Minuit object needed by cluster unfolding
	99	Float_t fR0 ; // Maximum distance between a EMC RecPoint and a PPSD RecPoint
	100	Bool_t fUnfoldFlag ; // Directive to unfold or not the clusters in case of multiple maxima
d15a28e7	101
b2a60966	102	ClassDef( AliPHOSTrackSegmentMakerv1,1) // Implementation version 1 of algorithm class to make PHOS track segments
d15a28e7	103
	104	};
	105
6ad0bfa0	106	#endif // AliPHOSTRACKSEGMENTMAKERV1_H