CbmStsClusterFinder.h

Go to the documentation of this file.

//* $Id: */
 
// -------------------------------------------------------------------------
// -----                     CbmStsClusterFinder header file       -----
// -----                  Created 26/06/2008  by R. Karabowicz         -----
// -------------------------------------------------------------------------
 
/** CbmStsClusterFinder
 *@author Volker Friese <v.friese@gsi.de>
 *@since 11.09.06
 *@version 1.0
 **
 ** CBM task class for finding hits in the STS
 ** Task level RECO
 ** Produces objects of type CbmStsHits out of CbmStsDigi.
 **/
 
#ifndef CBMSTSCLUSTERFINDER_H
#define CBMSTSCLUSTERFINDER_H 1
 
#include "FairTask.h"
#include "TStopwatch.h"
 
#include <map>
#include <set>
 
using std::map;
using std::pair;
using std::set;
 
class TClonesArray;
class CbmGeoStsPar;
class CbmStsDigiPar;
class CbmStsDigiScheme;
class CbmStsSector;
class CbmStsStation;
class CbmStsCluster;      // AZ
class CbmStsDigiScheme;   // GP
 
class CbmStsClusterFinder : public FairTask
{
  public:
    /** Default constructor **/
    CbmStsClusterFinder();
 
    /** Standard constructor **/
    CbmStsClusterFinder(Int_t iVerbose);
 
    /** Constructor with task name **/
    CbmStsClusterFinder(const char* name, Int_t iVerbose);
 
    /** Destructor **/
    virtual ~CbmStsClusterFinder();
 
    /** Execution **/
    virtual void Exec(Option_t* opt);
 
    /** Virtual method Finish **/
    virtual void Finish();
 
    Bool_t MC = false;
 
  private:
    CbmGeoStsPar* fGeoPar;                                       /** Geometry parameters **/
    CbmStsDigiPar* fDigiPar;                                     /** Digitisation parameters **/
    CbmStsDigiScheme* fDigiScheme; /** Digitisation scheme **/   // FIX GP
    TClonesArray* fDigis;                                        /** Input array of CbmStsDigi **/
    TClonesArray* fClustersCand;                                 /** Input array of CbmStsCluster **/
    TClonesArray* fClusters;                                     /** Input array of CbmStsCluster **/
    Int_t fNofDigis;
    map<CbmStsSector*, set<Int_t>> fDigiMapF; /** sector digis (front) **/
    map<CbmStsSector*, set<Int_t>> fDigiMapB; /** sector digis (back)  **/
    TStopwatch fTimer;
 
    Int_t fNofClustersCand;
    Int_t fNofClusters;
    Int_t fNofClustersGood;
    Int_t fNofClustersWP;
    Int_t fNofClustersWM;
    Int_t fLongestCluster;
    Int_t fLongestGoodCluster;
    Bool_t fUseLinks;
    /** Get parameter containers **/
    virtual void SetParContainers();
 
    /** Intialisation **/
    virtual InitStatus Init();
 
    /** Reinitialisation **/
    virtual InitStatus ReInit();
 
    /** Make sectorwise sets for sigis  **/
    void MakeSets();
 
    void AnalyzeClusters();
    void AnalyzeCluster(Int_t iclus);
    void EvalErrors(CbmStsCluster* clus, Int_t chan2);   // AZ
    void SplitCluster(Int_t iclus);                      // AZ
 
    /** Sort digis sectorwise  **/
    void SortDigis();
 
    /** Find hits in one sector **/
    /*   Int_t FindClusters(CbmStsStation* station, CbmStsSector* sector, */
    /*           set<Int_t>& fSet, set<Int_t>& bSet); */
 
    //  void FindClusters(Int_t stationNr, Int_t sectorNr, Int_t iSide, set<Int_t>& digiSet);
    Int_t FindClusters(Int_t stationNr, Int_t sectorNr, Int_t iSide, set<Int_t>& digiSet);
 
    CbmStsClusterFinder(const CbmStsClusterFinder&);
    CbmStsClusterFinder operator=(const CbmStsClusterFinder&);
 
    ClassDef(CbmStsClusterFinder, 1);
};
 
#endif