Line data    Source code

       1             : /**************************************************************************
       2             :  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
       3             :  *                                                                        *
       4             :  * Author: The ALICE Off-line Project.                                    *
       5             :  * Contributors are mentioned in the code where appropriate.              *
       6             :  *                                                                        *
       7             :  * Permission to use, copy, modify and distribute this software and its   *
       8             :  * documentation strictly for non-commercial purposes is hereby granted   *
       9             :  * without fee, provided that the above copyright notice appears in all   *
      10             :  * copies and that both the copyright notice and this permission notice   *
      11             :  * appear in the supporting documentation. The authors make no claims     *
      12             :  * about the suitability of this software for any purpose. It is          *
      13             :  * provided "as is" without express or implied warranty.                  *
      14             :  **************************************************************************/
      15             : //
      16             : // Input: ESDevent
      17             : // Functionality: find ESDtrack pairs according to some criteria to form one cosmic ray; AliTPCCosmicTrackfit then performs the combined track fit for the pair
      18             : //
      19             : //  Xianguo Lu 
      20             : //  lu@physi.uni-heidelberg.de
      21             : //  Xianguo.Lu@cern.ch
      22             : //
      23             : 
      24             : #ifndef ALICOSMICTRACKER_H
      25             : #define ALICOSMICTRACKER_H
      26             : 
      27             : #include "AliESDtrackCuts.h"
      28             : #include "AliExternalTrackParam.h"
      29             : #include "TTreeStream.h"
      30             : 
      31             : #ifndef TVECTOR3_H
      32             : #include "TVector3.h"
      33             : #endif
      34             : 
      35             : class AliESDCosmicTrack;
      36             : class AliTPCCosmicTrackfit;
      37             : 
      38             : class AliESDEvent;
      39             : 
      40             : class AliCosmicTracker
      41             : {
      42             :  public:
      43             : 
      44             :   typedef Bool_t (*CutFunc)(const AliESDtrack *trk);
      45             : 
      46             :   AliCosmicTracker(const Int_t dlev=0, const TString tag="test"); 
      47             :   ~AliCosmicTracker();
      48             : 
      49             :   void SetESDEvent(AliESDEvent *esd);
      50             :   Int_t Process(const TString tag="", const Bool_t kprint=kFALSE);
      51          16 :   TClonesArray * GetTrackStack() const {return fTrackStack;}
      52             : 
      53           0 :   TTreeSRedirector * GetStreamer() const {return fStreamer;}
      54           0 :   Int_t GetDebugLevel() const {return fDebugLevel;}
      55             :   Int_t GetErrFlag() const;
      56             : 
      57           0 :   void SetCutPull(     const Int_t ii, const Double_t cut){ fCutPull[ii] = cut;}
      58           0 :   void SetCutDelta(    const Int_t ii, const Double_t cut){ fCutDelta[ii] = cut;}
      59             : 
      60           0 :   void SetCutdPhi(const Double_t cut){fCutdPhi = cut;}
      61           0 :   void SetCutdTheta(const Double_t cut){fCutdTheta = cut;}
      62             : 
      63           0 :   void SetUserESDtrackCut(CutFunc func){fUserCut = func;}
      64             : 
      65             :   static TClonesArray *FindCosmic(AliESDEvent *event, const Bool_t kadd);
      66             : 
      67             :  private:
      68             :   AliCosmicTracker(const AliCosmicTracker &p);
      69             :   AliCosmicTracker & operator=(const AliCosmicTracker &p);
      70             : 
      71        1392 :   static Double_t CutFindable(){return 0.5;}          //cut on findable ratio of TPC cluster; DCA-anormlay is caused by laser!! should check trigger!!
      72             : 
      73             :   Bool_t ESDtrackCut(AliESDtrack * trk, Double_t &findabler);
      74             : 
      75             :   Bool_t IsPair(AliESDtrack* trk0, AliESDtrack*trk1);
      76             :   void WriteStreamer(Int_t ntrk, AliESDCosmicTrack *costrk);
      77             : 
      78             :   CutFunc fUserCut;                               //user ESDtrack function
      79             :   TTreeSRedirector *fStreamer;                   //debug streamer
      80             :   Int_t fDebugLevel;                             //debug level
      81             : 
      82             :   AliESDEvent *fESDEvent;                        //esd
      83             :   AliTPCCosmicTrackfit *fCosmicTrackfit;          //cosmictrackfit
      84             :   TClonesArray *fTrackStack;                     //storing cosmic ray
      85             : 
      86             :   AliExternalTrackParam fTrack0;                //upper track param estimated at x=0 from IsPair
      87             :   AliExternalTrackParam fTrack1;                //lower track param estimated at x=0 from IsPair
      88             : 
      89             :   TVector3 fRawVtx;                             //raw vertex position calculated only from upper and lower inner most TPC cluster
      90             :   Double_t fRawDCA;                              //raw DCA (2d) calculated only from upper and lower inner most TPC cluster 
      91             :   Double_t fdPhi;                                //phi0-phi1-pi of the EDS tracks
      92             :   Double_t fCutdPhi;                           //cut
      93             : 
      94             :   Double_t fdTheta;                              //theta0+theta1-pi of the ESD tracks
      95             :   Double_t fCutdTheta;                         //cut
      96             : 
      97             :   Double_t fPull[5];                             //pull of the two ESD tracks at x=0
      98             :   Double_t fCutPull[5];                            //pull
      99             : 
     100             :   Double_t fDelta[5];                             //delta of the two ESD tracks at x=0
     101             :   Double_t fCutDelta[5];                          //delta cut
     102             : 
     103             :   Int_t fErrFlagESDtrackCut;                       //error status in ESDtrackCut()
     104             :   Int_t fErrFlagIsPair;                            //error status in IsPair()
     105             :   Int_t fErrFlagCosmicTrackfit;                    //error status in fCosmicTrackfit
     106             :   //
     107             :   Bool_t fBFieldOn;                                // flag field ON event
     108             : };
     109             : 
     110             : #endif
     111             : 
     112             :