Dimitri just makes it work
[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  
51   Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0
52
53   void    MakeLinks(TArrayI * EmcRecPoints, TArrayI * PpsdRecPointsUp, TArrayI * PpsdRecPointsLow, 
54                     TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
55   void    MakePairs(TArrayI * EmcRecPoints, 
56                     TArrayI * PpsdRecPointsUp, 
57                     TArrayI * 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 MakeTrackSegmentsCPV(DigitsList * DL, 
66                                  AliPHOSRecPoint::RecPointsList * emcl, 
67                                  AliPHOSRecPoint::RecPointsList * ppsdl ); // just unfold EMC and CPV clusters
68   virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
69   virtual void    SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ; 
70   static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function)
71   void    UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ; 
72                                              // Unfolds and sorts all EMC clusters
73   void  UnfoldClusters(DigitsList * DL, 
74                        AliPHOSRecPoint::RecPointsList * emcIn, 
75                        AliPHOSEmcRecPoint * iniEmc, 
76                        Int_t Nmax, 
77                        int * maxAt, 
78                        Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
79   virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; } 
80
81   AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & rvalue)  {
82     // assignement operator requested by coding convention
83     // but not needed
84     assert(0==1) ;
85     return *this ; 
86   }
87
88 private:
89
90   Float_t fDelta ;     // parameter used for sorting
91   Float_t fR0 ;        // Maximum distance between a EMC RecPoint and a PPSD RecPoint   
92   Bool_t fUnfoldFlag ; // Directive to unfold or not the clusters in case of multiple maxima
93
94   ClassDef( AliPHOSTrackSegmentMakerv1,1)  // Implementation version 1 of algorithm class to make PHOS track segments 
95
96 };
97
98 #endif // AliPHOSTRACKSEGMENTMAKERV1_H