6ad0bfa0 |
1 | #ifndef ALIPHOSTRACKSEGMENTMAKERV1_H |
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$ */ |
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 |
98cbd830 |
43 | void FillOneModule(AliPHOSRecPoint::RecPointsList * emcIn, |
28c3a259 |
44 | TArrayI * emcOut, |
88714635 |
45 | AliPHOSRecPoint::RecPointsList * ppsdIn, |
28c3a259 |
46 | TArrayI * ppsdOutUp, |
47 | TArrayI * ppsdOutLow, |
88714635 |
48 | Int_t &PHOSModule, |
49 | Int_t & emcStopedAt, |
0336aa72 |
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 |
d15a28e7 |
58 | Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0 |
59 | |
28c3a259 |
60 | void MakeLinks(TArrayI * EmcRecPoints, TArrayI * PpsdRecPointsUp, TArrayI * PpsdRecPointsLow, |
d15a28e7 |
61 | TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD |
28c3a259 |
62 | void MakePairs(TArrayI * EmcRecPoints, |
63 | TArrayI * PpsdRecPointsUp, |
64 | TArrayI * PpsdRecPointsLow, |
88714635 |
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 |
0336aa72 |
72 | void MakeTrackSegmentsCPV(DigitsList * DL, |
73 | AliPHOSRecPoint::RecPointsList * emcl, |
74 | AliPHOSRecPoint::RecPointsList * ppsdl ) ; // just unfold EMC and CPV clusters |
d912c4d7 |
75 | virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;} |
6727be7e |
76 | virtual void SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ; |
9f616d61 |
77 | static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function) |
98cbd830 |
78 | void UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ; |
0336aa72 |
79 | // Unfolds and sorts all EMC clusters |
88714635 |
80 | void UnfoldClusters(DigitsList * DL, |
81 | AliPHOSRecPoint::RecPointsList * emcIn, |
82 | AliPHOSEmcRecPoint * iniEmc, |
83 | Int_t Nmax, |
84 | int * maxAt, |
98cbd830 |
85 | Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package |
d912c4d7 |
86 | virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; } |
d15a28e7 |
87 | |
839ffcb3 |
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 | |
d15a28e7 |
95 | private: |
96 | |
b2a60966 |
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 |
d15a28e7 |
101 | |
b2a60966 |
102 | ClassDef( AliPHOSTrackSegmentMakerv1,1) // Implementation version 1 of algorithm class to make PHOS track segments |
d15a28e7 |
103 | |
104 | }; |
105 | |
6ad0bfa0 |
106 | #endif // AliPHOSTRACKSEGMENTMAKERV1_H |