fe37af8aa9a6f81a64ea919449192bef84090b1b
[u/mrichter/AliRoot.git] / PHOS / AliPHOSTrackSegmentMakerv1.h
1 #ifndef ALIPHOSTRACKSEGMENTMAKERV1_H
2 #define ALIPHOSTRACKSEGMENTMAKERV1_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
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)
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"
26 #include "TMinuit.h" 
27
28 class  AliPHOSTrackSegmentMakerv1 : public AliPHOSTrackSegmentMaker {
29
30 public:
31
32   AliPHOSTrackSegmentMakerv1() ;                     
33   AliPHOSTrackSegmentMakerv1(const AliPHOSTrackSegmentMakerv1 & tsm) {
34     // cpy ctor: no implementation yet
35     // requested by the Coding Convention
36     assert(0==1) ; 
37   }
38    
39   virtual ~ AliPHOSTrackSegmentMakerv1() ; // dtor
40   
41   Bool_t  FindFit(AliPHOSEmcRecPoint * emcRP, int * MaxAt, Float_t * maxAtEnergy, 
42                   Int_t NPar, Float_t * FitParametres) ; //Used in UnfoldClusters, calls TMinuit
43   void    FillOneModule(AliPHOSRecPoint::RecPointsList * emcIn, 
44                         TArrayI * emcOut, 
45                         AliPHOSRecPoint::RecPointsList * ppsdIn, 
46                         TArrayI * ppsdOutUp, 
47                         TArrayI * ppsdOutLow, 
48                         Int_t &PHOSModule, 
49                         Int_t & emcStopedAt, 
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
58   Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0
59
60   void    MakeLinks(TArrayI * EmcRecPoints, TArrayI * PpsdRecPointsUp, TArrayI * PpsdRecPointsLow, 
61                     TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
62   void    MakePairs(TArrayI * EmcRecPoints, 
63                     TArrayI * PpsdRecPointsUp, 
64                     TArrayI * PpsdRecPointsLow, 
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
72   void    MakeTrackSegmentsCPV(DigitsList * DL, 
73                             AliPHOSRecPoint::RecPointsList * emcl, 
74                             AliPHOSRecPoint::RecPointsList * ppsdl ) ; // just unfold EMC and CPV clusters
75   virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
76   virtual void    SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ; 
77   static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function)
78   void    UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ; 
79                                             // Unfolds and sorts all EMC clusters
80   void  UnfoldClusters(DigitsList * DL, 
81                        AliPHOSRecPoint::RecPointsList * emcIn, 
82                        AliPHOSEmcRecPoint * iniEmc, 
83                        Int_t Nmax, 
84                        int * maxAt, 
85                        Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
86   virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; } 
87
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
95 private:
96
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
101
102   ClassDef( AliPHOSTrackSegmentMakerv1,1)  // Implementation version 1 of algorithm class to make PHOS track segments 
103
104 };
105
106 #endif // AliPHOSTRACKSEGMENTMAKERV1_H