CbmStsRealDigitize.h

Go to the documentation of this file.

//* $Id: */
 
// -------------------------------------------------------------------------
// -----                    CbmStsRealDigitize header file             -----
// -----                  Created 01/07/2008  by R. Karabowicz         -----
// -------------------------------------------------------------------------
 
/** CbmStsDigitise
 *@author Volker Friese <v.friese@gsi.de>
 *@since 30.08.06
 *@version 1.0
 **
 ** CBM task class for digitising STS
 ** Task level SIM
 ** Produces objects of type CbmStsDigi out of CbmStsPoint.
 **/
 
#ifndef CBMSTSREALDIGITIZE_H
#define CBMSTSREALDIGITIZE_H 1
 
#include "FairTask.h"
#include "TStopwatch.h"
 
#include <set>
 
using std::map;
using std::pair;
using std::set;
 
class TClonesArray;
class CbmGeostsPar;
class CbmStsDigiPar;
class CbmStsStation;
 
class CbmStsPoint;
 
class CbmStsRealDigitize : public FairTask
{
 
  public:
    /** Default constructor **/
    CbmStsRealDigitize();
 
    /** Standard constructor **/
    CbmStsRealDigitize(Int_t iVerbose);
 
    /** Constructor with name **/
    CbmStsRealDigitize(const char* name, Int_t iVerbose);
 
    /** Destructor **/
    virtual ~CbmStsRealDigitize();
 
    void SetStep(Double_t tStep) { fStep = tStep; }
 
    /** Execution **/
    virtual void Exec(Option_t* opt);
 
    /** Virtual method Finish **/
    virtual void Finish();
 
    virtual void SetRealisticResponse(Bool_t real = kTRUE) { fRealistic = real; }
 
    void FindFiredStrips(CbmStsPoint* pnt, Int_t& nofStr, Int_t*& strips, Double_t*& signals, Int_t side);
    void ProduceHitResponse(CbmStsSensor* sensor);
 
    void SetFrontThreshold(Double_t frontThr = 0.) { fFThreshold = frontThr; }
    void SetBackThreshold(Double_t backThr = 0.) { fBThreshold = backThr; }
    void SetFrontNoiseWidth(Double_t frontNoW = 0.) { fFNoiseWidth = frontNoW; }
    void SetBackNoiseWidth(Double_t backNoW = 0.) { fBNoiseWidth = backNoW; }
 
    void SetFrontNofBits(Int_t frontNB = 0) { fFNofBits = frontNB; }
    void SetBackNofBits(Int_t backNB = 0) { fBNofBits = backNB; }
    void SetFrontMinStep(Double_t frontMS = 0.) { fFMinStep = frontMS; }
    void SetBackMinStep(Double_t backMS = 0.) { fBMinStep = backMS; }
 
    void SetStripDeadTime(Double_t StripDeadTime = 0.) { fStripDeadTime = StripDeadTime; }
 
  private:
    CbmGeoStsPar* fGeoPar;         /** Geometry parameter container **/
    CbmStsDigiPar* fDigiPar;       /** Digitisation parameter container **/
    CbmStsDigiScheme* fDigiScheme; /** Digitisation scheme **/
    TClonesArray* fPoints;         /** Input array of CbmStsPoint **/
    TClonesArray* fDigis;          /** Output array of CbmStsDigi **/
    TClonesArray* fDigiMatches;    /** Output array of CbmStsDigiMatches**/
 
    // statistics
    Int_t fNDigis;
    Int_t fNMulti;
 
    Double_t fNEvents;
    Double_t fNPoints;
    Double_t fNOutside;
    Double_t fNDigisFront;
    Double_t fNDigisBack;
    Double_t fTime;
 
    // settings
    Double_t fStep;
 
    TStopwatch fTimer;
 
    Bool_t fRealistic;
    Double_t fEnergyLossToSignal;
 
    Double_t fFThreshold;
    Double_t fBThreshold;
    Double_t fFNoiseWidth;
    Double_t fBNoiseWidth;
 
    Double_t fStripDeadTime;
 
    Int_t fFNofBits;
    Int_t fBNofBits;
    Double_t fFMinStep;
    Double_t fBMinStep;
    Int_t fFNofSteps;
    Int_t fBNofSteps;
 
    Double_t* fStripSignalF;
    Double_t* fStripSignalB;
 
    // map of the point indices contributing to a fired strip
    map<Int_t, set<Int_t>> fFChannelPointsMap;
    map<Int_t, set<Int_t>> fBChannelPointsMap;
 
    map<CbmStsSensor*, set<Int_t>> fPointMap; /** sensor points **/
    Float_t occupancy[10][1000][20];
 
    /** Make sensorwise set for points **/
    void MakeSets();
 
    /** Get parameter containers **/
    virtual void SetParContainers();
 
    /** Intialisation **/
    virtual InitStatus Init();
 
    /** Reinitialisation **/
    virtual InitStatus ReInit();
 
    /** Reset eventwise counters **/
    void Reset();
 
    ClassDef(CbmStsRealDigitize, 1);
};
 
#endif