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