the MIXT geometry (IHEP+GPS2) has been introduced
[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                             AliPHOSRecPoint::RecPointsList * cpvl, 
65                             AliPHOSTrackSegment::TrackSegmentsList * trsl ) ; // does the job
66   virtual void MakeTrackSegmentsCPV(DigitsList * DL, 
67                                  AliPHOSRecPoint::RecPointsList * emcl, 
68                                  AliPHOSRecPoint::RecPointsList * ppsdl ); // just unfold EMC and CPV clusters
69   virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
70   virtual void    SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ; 
71   static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function)
72   void    UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ; 
73                                              // Unfolds and sorts all EMC clusters
74   void  UnfoldClusters(DigitsList * DL, 
75                        AliPHOSRecPoint::RecPointsList * emcIn, 
76                        AliPHOSEmcRecPoint * iniEmc, 
77                        Int_t Nmax, 
78                        int * maxAt, 
79                        Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
80   virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; } 
81
82   AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & rvalue)  {
83     // assignement operator requested by coding convention
84     // but not needed
85     assert(0==1) ;
86     return *this ; 
87   }
88
89 private:
90
91   Float_t fDelta ;     // parameter used for sorting
92   Float_t fR0 ;        // Maximum distance between a EMC RecPoint and a PPSD RecPoint   
93   Bool_t fUnfoldFlag ; // Directive to unfold or not the clusters in case of multiple maxima
94
95   ClassDef( AliPHOSTrackSegmentMakerv1,1)  // Implementation version 1 of algorithm class to make PHOS track segments 
96
97 };
98
99 #endif // AliPHOSTRACKSEGMENTMAKERV1_H