Unfolding slightly reorganized (added finction UnfoldAll) to be consistent with Index...
[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                         TObjArray * emcOut, 
45                         AliPHOSRecPoint::RecPointsList * ppsdIn, 
46                         TObjArray * ppsdOutUp, 
47                         TObjArray * ppsdOutLow, 
48                         Int_t &PHOSModule, 
49                         Int_t & emcStopedAt, 
50                         Int_t & ppsdStopedAt) ; // Fills temporary arrais with clusters from one module  
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
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
65   virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
66   virtual void    SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ; 
67   static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function)
68   void    UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ; 
69                                              // Unfolds and sorts all EMC clusters
70   void  UnfoldClusters(DigitsList * DL, 
71                        AliPHOSRecPoint::RecPointsList * emcIn, 
72                        AliPHOSEmcRecPoint * iniEmc, 
73                        Int_t Nmax, 
74                        int * maxAt, 
75                        Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
76   virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; } 
77
78   AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & rvalue)  {
79     // assignement operator requested by coding convention
80     // but not needed
81     assert(0==1) ;
82     return *this ; 
83   }
84
85 private:
86
87   Float_t fDelta ;     // parameter used for sorting
88   TMinuit * fMinuit ;  // Minuit object needed by cluster unfolding
89   Float_t fR0 ;        // Maximum distance between a EMC RecPoint and a PPSD RecPoint   
90   Bool_t fUnfoldFlag ; // Directive to unfold or not the clusters in case of multiple maxima
91
92   ClassDef( AliPHOSTrackSegmentMakerv1,1)  // Implementation version 1 of algorithm class to make PHOS track segments 
93
94 };
95
96 #endif // AliPHOSTRACKSEGMENTMAKERV1_H