]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PHOS/AliPHOSTrackSegmentMakerv1.h
DP:Calculation of cluster properties dep. on vertex added
[u/mrichter/AliRoot.git] / PHOS / AliPHOSTrackSegmentMakerv1.h
index 8826c35afca97f8b19977518948864084e43dce2..2eb9253373b4e3b46087599b817b319166f5f1c5 100644 (file)
-#ifndef ALIPHOSSUBTRACKERV1_H
-#define ALIPHOSSUBTRACKERV1_H
+#ifndef ALIPHOSTRACKSEGMENTMAKERV1_H
+#define ALIPHOSTRACKSEGMENTMAKERV1_H
 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
  * See cxx source for full Copyright notice                               */
 
-///////////////////////////////////////////////////
-//  Track Segment Maker class for PHOS           //
-//  Version SUBATECH                             //
-//  Author Dmitri Peressounko RRC Ki             //
-//     comment: finds pairs of clusters EMC+PPSD //  
-//              performs unfolding.              //
-///////////////////////////////////////////////////
+/* $Id$ */
+/* History of cvs commits:
+ *
+ * $Log$
+ * Revision 1.49  2007/02/01 13:59:11  hristov
+ * Forward declaration
+ *
+ * Revision 1.48  2006/08/28 10:01:56  kharlov
+ * Effective C++ warnings fixed (Timur Pocheptsov)
+ *
+ * Revision 1.47  2005/11/17 12:35:27  hristov
+ * Use references instead of objects. Avoid to create objects when they are not really needed
+ *
+ * Revision 1.46  2005/05/28 14:19:05  schutz
+ * Compilation warnings fixed by T.P.
+ *
+ */
+
+//_________________________________________________________________________
+// Implementation version 1 of algorithm class to construct PHOS track segments
+// Associates EMC and CPV lusters
+// Unfolds the EMC cluster   
+//                  
+//*-- Author: Dmitri Peressounko (RRC Ki & SUBATECH)
 
 // --- ROOT system ---
+#include <TVector3.h>
 
 // --- Standard library ---
 
 // --- AliRoot header files ---
-
-#include "TObjArray.h"
-#include "AliPHOSClusterizer.h"
-#include "AliPHOSEmcRecPoint.h"
-#include "AliPHOSPpsdRecPoint.h"
 #include "AliPHOSTrackSegmentMaker.h"
 
+class AliPHOSEmcRecPoint ;
+class AliPHOSCpvRecPoint ;
+class TClonesArray;
 
 class  AliPHOSTrackSegmentMakerv1 : public AliPHOSTrackSegmentMaker {
 
 public:
 
   AliPHOSTrackSegmentMakerv1() ;                     
-  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(DigitsList * Dl, RecPointsList * emcIn, TObjArray * emcOut, RecPointsList * ppsdIn, 
-                       TObjArray * ppsdOutUp, TObjArray * ppsdOutLow, Int_t &PHOSModule, Int_t & emcStopedAt, 
-                       Int_t & ppsdStopedAt) ; // Unfolds clusters and fills temporary arrais   
-
-  Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSPpsdRecPoint * Ppsd , Bool_t & TooFar ) ; // see R0
-
-  void    MakeLinks(TObjArray * EmcRecPoints, TObjArray * PpsdRecPointsUp, TObjArray * PpsdRecPointsLow, 
-                   TClonesArray * LinkLowArray, TClonesArray *LinkUpArray) ; //Evaluates distances(links) between EMC and PPSD
-
-  void    MakePairs(TObjArray * EmcRecPoints, TObjArray * PpsdRecPointsUp, TObjArray * PpsdRecPointsLow, 
-                   TClonesArray * LinkLowArray, TClonesArray * LinkUpArray, TrackSegmentsList * trsl) ; 
-                    //Finds pairs(triplets) with smallest link
-
-  void    MakeTrackSegments(DigitsList * DL, RecPointsList * emcl, RecPointsList * ppsdl, TrackSegmentsList * trsl ) ; // does the job
-
-  void    SetMaxEmcPpsdDistance(Float_t r){ fR0 = r ;} //Radius within which we look for ppsd cluster
-
- Double_t ShowerShape(Double_t r) ; //Shape of shower used in unfolding
-
-  void    UnfoldClusters(DigitsList * DL, RecPointsList * emcIn, AliPHOSEmcRecPoint * iniEmc, Int_t Nmax, 
-                        int * maxAt, Float_t * maxAtEnergy, TObjArray * emclist) ; //Unfolds overlaping clusters using TMinuit packadge
-
-  void static UnfoldingChiSquare(Int_t &NPar, Double_t *Grad, Double_t & fret, Double_t *x, Int_t iflag); //used in TMinuit
+  AliPHOSTrackSegmentMakerv1(const TString & alirunFileNameFile, const TString & eventFolderName = AliConfig::GetDefaultEventFolderName());
+  AliPHOSTrackSegmentMakerv1(const AliPHOSTrackSegmentMakerv1 & tsm);
 
+  virtual ~ AliPHOSTrackSegmentMakerv1() ; // dtor
+  
+  //  virtual char*  GetRecPointsBranch    (void)const{return (char*)fRecPointsBranchTitle.Data() ;}
+  //  virtual char*  GetTrackSegmentsBranch(void)const{return (char*)fTrackSegmentsBranchTitle.Data() ;}
+  virtual Int_t GetTrackSegmentsInRun()const {return fTrackSegmentsInRun ;}  
+
+  virtual void   Exec(Option_t *option); // Does the job
+          void   FillOneModule() ;       // Finds range in which RecPoints belonging current PHOS module are
+
+          void   MakeLinks() const;      //Evaluates distances(links) between EMC and CPV
+          void   MakePairs() ;           //Finds pairs(triplets) with smallest link
+  virtual void   Print(const Option_t * = "") const ;
+  //Switch to "on flyght" mode, without writing to TreeR and file  
+  void SetWriting(Bool_t toWrite = kFALSE){fWrite = toWrite;} 
+  virtual void   SetMaxEmcCPVDistance(Float_t r){ fRcpv = r ;} //Maximal distance (in PHOS plane) 
+                                                               //between EMCrp and CPVrp
+  virtual void   SetMaxCPVTPCDistance(Float_t r){ fRtpc = r ;} //Maximal distance 
+                                                               //between EMCrp and extrapolation of TPC track
+  //  virtual void   SetRecPointsBranch(const char * title) { fRecPointsBranchTitle = title ;} 
+  //  virtual void   SetTrackSegmentsBranch(const char * title){ fTrackSegmentsBranchTitle = title ; }
+  virtual const char * Version() const { return "tsm-v1" ; }  
+
+  AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & )  {
+    // assignement operator requested by coding convention but not needed
+    Fatal("operator =", "not implemented") ;
+    return *this ; 
+  }
+  void Unload() ;
 
 private:
 
-  Float_t fDelta ;  // parameter used for sorting
-  Float_t fR0  ;    // Maximal distance between EMC and PPSD clusters of one Track Segment in module plane
+  const TString BranchName() const ; 
+  void  GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * EmcClu , AliPHOSCpvRecPoint * Cpv, 
+                               Int_t & track, Float_t &dx, Float_t &dz ) const ; // see R0
+  void    Init() ;
+  void    InitParameters() ;
+  void    PrintTrackSegments(Option_t *option) ;
+  virtual void   WriteTrackSegments() ;
+  void    GetVertex(void) ;
+  void    EvalRecPoints(void) ;
+
+private:  
+
+  Bool_t  fDefaultInit;               //! Says if the task was created by defaut ctor (only parameters are initialized)
+  Bool_t  fWrite ;                   // Write Tracks to TreeT  
+  Int_t fNTrackSegments ; // number of track segments found 
+
+  Float_t fRcpv ;        // Maximum distance between a EMC RecPoint and a CPV RecPoint   
+  Float_t fRtpc ;        // Maximum distance between a EMC RecPoint and extrapolation of a TPC track   
+  
+  TVector3 fVtx ;        //! Vertex in current position
 
-public: 
+  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 fModule ;      //! number of module being processed
+  Int_t fTrackSegmentsInRun ; //! Total number of track segments in one run
 
-ClassDef( AliPHOSTrackSegmentMakerv1,1)  // track segment maker implementation , version 1
+  ClassDef( AliPHOSTrackSegmentMakerv1,3)  // Implementation version 1 of algorithm class to make PHOS track segments 
 
-};
+ };
 
-#endif // AliPHOSSUBTRACKERV1_H
+#endif // AliPHOSTRACKSEGMENTMAKERV1_H