-#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