]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PHOS/AliPHOSTrackSegmentMakerv1.h
can now read from HPSS
[u/mrichter/AliRoot.git] / PHOS / AliPHOSTrackSegmentMakerv1.h
index ce203ded756901623a7fc541ac376ad64a407e28..2e58f7b54a2e8d983dc68b751b757d9ff1576b18 100644 (file)
 // --- Standard library ---
 
 // --- AliRoot header files ---
-
-#include "TObjArray.h"
-#include "AliPHOSClusterizer.h"
-#include "AliPHOSEmcRecPoint.h"
-#include "AliPHOSPpsdRecPoint.h"
 #include "AliPHOSTrackSegmentMaker.h"
-#include "TMinuit.h" 
+
+class AliPHOSEmcRecPoint ;
+class AliPHOSRecPoint ;
+
 
 class  AliPHOSTrackSegmentMakerv1 : public AliPHOSTrackSegmentMaker {
 
 public:
 
   AliPHOSTrackSegmentMakerv1() ;                     
+  AliPHOSTrackSegmentMakerv1(const char* headerFile,const char* branchTitle = 0) ;                     
   AliPHOSTrackSegmentMakerv1(const AliPHOSTrackSegmentMakerv1 & tsm) {
     // cpy ctor: no implementation yet
     // requested by the Coding Convention
-    assert(0==1) ; 
+    abort() ; 
   }
    
   virtual ~ AliPHOSTrackSegmentMakerv1() ; // dtor
   
-  Bool_t  FindFit(AliPHOSEmcRecPoint * emcRP, int * MaxAt, Float_t * maxAtEnergy, 
-                 Int_t NPar, Float_t * FitParametres) ; //Used in UnfoldClusters, calls TMinuit
-  void    FillOneModule(AliPHOSRecPoint::RecPointsList * emcIn, 
-                       TArrayI * emcOut, 
-                       AliPHOSRecPoint::RecPointsList * ppsdIn, 
-                       TArrayI * ppsdOutUp, 
-                       TArrayI * ppsdOutLow, 
-                       Int_t &PHOSModule, 
-                       Int_t & emcStopedAt, 
-                       Int_t & ppsdStopedAt) ; // Fills temporary arrais with clusters from one module  
-  Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0
-
-  void    MakeLinks(TArrayI * EmcRecPoints, TArrayI * PpsdRecPointsUp, TArrayI * PpsdRecPointsLow, 
-                   TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
-  void    MakePairs(TArrayI * EmcRecPoints, 
-                   TArrayI * PpsdRecPointsUp, 
-                   TArrayI * PpsdRecPointsLow, 
-                   TClonesArray * LinkLowArray, 
-                   TClonesArray * LinkUpArray, 
-                   AliPHOSTrackSegment::TrackSegmentsList * trsl) ; //Finds pairs(triplets) with smallest link
-  void    MakeTrackSegments(DigitsList * DL, 
-                           AliPHOSRecPoint::RecPointsList * emcl, 
-                           AliPHOSRecPoint::RecPointsList * ppsdl, 
-                           AliPHOSTrackSegment::TrackSegmentsList * trsl ) ; // does the job
-  virtual void MakeTrackSegmentsCPV(DigitsList * DL, 
-                                AliPHOSRecPoint::RecPointsList * emcl, 
-                                AliPHOSRecPoint::RecPointsList * ppsdl ); // just unfold EMC and CPV clusters
-  virtual void SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
-  virtual void    SetUnfoldFlag() { fUnfoldFlag = kTRUE ; } ; 
-  static Double_t ShowerShape(Double_t r) ; // Shape of shower used in unfolding; class member function (not object member function)
-  void    UnfoldAll(DigitsList * Dl, AliPHOSRecPoint::RecPointsList * emcIn) ; 
-                                             // Unfolds and sorts all EMC clusters
-  void  UnfoldClusters(DigitsList * DL, 
-                      AliPHOSRecPoint::RecPointsList * emcIn, 
-                      AliPHOSEmcRecPoint * iniEmc, 
-                      Int_t Nmax, 
-                      int * maxAt, 
-                      Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
-  virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; } 
-
-  AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & rvalue)  {
+  virtual char*  GetRecPointsBranch    (void)const{return (char*)fRecPointsBranchTitle.Data() ;}
+  virtual char*  GetTrackSegmentsBranch(void)const{return (char*)fTSBranchTitle.Data() ;}
+
+  virtual void   Exec(Option_t * option) ;
+          void   FillOneModule() ;       // Finds range in which RecPoints belonging current PHOS module are
+
+          void   MakeLinks() ;           //Evaluates distances(links) between EMC and PPSD
+          void   MakePairs() ;           //Finds pairs(triplets) with smallest link
+  virtual void   Print(Option_t * option) const ;
+  virtual Bool_t ReadRecPoints() ;
+  virtual void   SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;}
+  virtual void   SetRecPointsBranch(const char * title) { fRecPointsBranchTitle = title ;} 
+  virtual void   SetTrackSegmentsBranch(const char * title){ fTSBranchTitle = title ; }
+  virtual void   WriteTrackSegments() ;
+
+  AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & )  {
     // assignement operator requested by coding convention
     // but not needed
-    assert(0==1) ;
+    abort() ;
     return *this ; 
   }
 
 private:
+  Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0
+  void    Init() ;
+  void    PrintTrackSegments(Option_t *option) ;
+
+private:  
+
+  TString fHeaderFileName ;          // name of the file which contains gAlice, Tree headers etc.
+  TString fRecPointsBranchTitle ; // name of the file, where RecPoints branchs are stored
+  TString fTSBranchTitle ;        // name of the file, where TrackSegment branchs is stored
+  AliPHOSClusterizer * fClusterizer ; // !  
+  Int_t                  fNTrackSegments ; // number of track segments found 
+  AliPHOSGeometry      * fGeom ;           //! pointer to PHOS geometry  
+  Int_t          fEvent ;            // ! event being precessed
+  TObjArray    * fEmcRecPoints ;     // ! List of EMC Rec Points
+  TObjArray    * fCpvRecPoints ;     // ! List of CPV/PPSD recPoints
+  TClonesArray * fTrackSegments;     // ! list of final track segments
+
+
+  Bool_t  fIsInitialized ; //
 
-  Float_t fDelta ;     // parameter used for sorting
   Float_t fR0 ;        // Maximum distance between a EMC RecPoint and a PPSD RecPoint   
-  Bool_t fUnfoldFlag ; // Directive to unfold or not the clusters in case of multiple maxima
+
+  TClonesArray * fLinkLowArray ;  //!
+  TClonesArray * fLinkUpArray  ;  //!
+
+
+  Int_t fEmcFirst;     //! Index of first EMC RecPoint belonging to currect PHOS module
+  Int_t fEmcLast ;     //!
+  Int_t fCpvFirst;     //! Cpv upper layer     
+  Int_t fCpvLast;      //! 
+  Int_t fPpsdFirst;    //! Cpv low layer     
+  Int_t fPpsdLast;     //!
+  Int_t fModule ;      //! number of module being processed
 
   ClassDef( AliPHOSTrackSegmentMakerv1,1)  // Implementation version 1 of algorithm class to make PHOS track segments