]>
Commit | Line | Data |
---|---|---|
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 |