[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(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn, 
			TObjArray * emcOut, 
			AliPHOSRecPoint::RecPointsList * ppsdIn, 
			TObjArray * ppsdOutUp, 
			TObjArray * ppsdOutLow, 
			Int_t &PHOSModule, 
			Int_t & emcStopedAt, 
			Int_t & ppsdStopedAt) ; // Unfolds clusters and fills temporary arrais   
  Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0

  void    MakeLinks(TObjArray * EmcRecPoints, TObjArray * PpsdRecPointsUp, TObjArray * PpsdRecPointsLow, 
		    TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
  void    MakePairs(TObjArray * EmcRecPoints, 
		    TObjArray * PpsdRecPointsUp, 
		    TObjArray * 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  UnfoldClusters(DigitsList * DL, 
		       AliPHOSRecPoint::RecPointsList * emcIn, 
		       AliPHOSEmcRecPoint * iniEmc, 
		       Int_t Nmax, 
		       int * maxAt, 
		       Float_t * maxAtEnergy, 
		       TObjArray * emclist) ; //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
88714635	43	void FillOneModule(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn,
	44	TObjArray * emcOut,
	45	AliPHOSRecPoint::RecPointsList * ppsdIn,
	46	TObjArray * ppsdOutUp,
	47	TObjArray * ppsdOutLow,
	48	Int_t &PHOSModule,
	49	Int_t & emcStopedAt,
d15a28e7	50	Int_t & ppsdStopedAt) ; // Unfolds clusters and fills temporary arrais
d15a28e7	51	Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0
	52
	53	void MakeLinks(TObjArray * EmcRecPoints, TObjArray * PpsdRecPointsUp, TObjArray * PpsdRecPointsLow,
	54	TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
88714635	55	void MakePairs(TObjArray * EmcRecPoints,
	56	TObjArray * PpsdRecPointsUp,
	57	TObjArray * PpsdRecPointsLow,
	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)
88714635	68	void UnfoldClusters(DigitsList * DL,
	69	AliPHOSRecPoint::RecPointsList * emcIn,
	70	AliPHOSEmcRecPoint * iniEmc,
	71	Int_t Nmax,
	72	int * maxAt,
	73	Float_t * maxAtEnergy,
	74	TObjArray * emclist) ; //Unfolds overlaping clusters using TMinuit package
d912c4d7	75	virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; }
d15a28e7	76
839ffcb3	77	AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & rvalue) {
	78	// assignement operator requested by coding convention
	79	// but not needed
	80	assert(0==1) ;
	81	return *this ;
	82	}
	83
d15a28e7	84	private:
d15a28e7	85
b2a60966	86	Float_t fDelta ; // parameter used for sorting
	87	TMinuit * fMinuit ; // Minuit object needed by cluster unfolding
	88	Float_t fR0 ; // Maximum distance between a EMC RecPoint and a PPSD RecPoint
	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