ERBeamDetDigitizer.h

/********************************************************************************
 *              Copyright (C) Joint Institute for Nuclear Research              *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *         GNU Lesser General Public Licence version 3 (LGPL) version 3,        *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/

#ifndef ERBeamDetDigitizer_H
#define ERBeamDetDigitizer_H

#include "TClonesArray.h"

#include "ERDigitizer.h"
#include "ERBeamDetMWPCPoint.h"
#include "ERBeamDetTOFPoint.h"
#include "ERDigi.h"
#include "ERSupport.h"

class ERBeamDetDigitizer : public ERDigitizer {
public:
  ERBeamDetDigitizer();
  ERBeamDetDigitizer(Int_t verbose);
  virtual ~ERBeamDetDigitizer() = default;
  void SetToFElossThreshold(Float_t treshold){fToFElossThreshold = treshold;}
  void SetMWPCElossThreshold(Float_t treshold){fMWPCElossThreshold = treshold;}
  void SetToFElossSigma(Float_t sigma) {fElossSigmaToF = sigma;}
  void SetMWPCElossSigma(Float_t sigma) {fElossSigmaMWPC = sigma;}
  void SetToFTimeSigma(Float_t sigma) {fTimeSigmaToF = sigma;}
  void SetMWPCTimeSigma(Float_t sigma) {fTimeSigmaMWPC = sigma;}
  void SetToFElossSigmaOverEloss(Float_t a) {
    fElossSigmaOverElossToF = a; 
    fSigmaEOverEToFIsSet = kTRUE;
  }
public:
  virtual InitStatus Init();
  virtual void Exec(Option_t* opt);
  virtual void Reset();
protected:
  ERDigi* AddMWPCDigi(float edep, float time, int mwpcNb, int planeNb, ERChannel wireNb);
  ERDigi*  AddToFDigi(float edep, float time, int tofNb);
protected:
  //Input arrays
  TClonesArray* fBeamDetToFPoints = nullptr;  
  TClonesArray* fBeamDetMWPCPoints = nullptr; 
  //Output arrays
  TClonesArray* fBeamDetToFDigi1 = nullptr; 
  TClonesArray* fBeamDetToFDigi2 = nullptr; 
  TClonesArray* fBeamDetMWPCDigiX1 = nullptr;
  TClonesArray* fBeamDetMWPCDigiX2 = nullptr;
  TClonesArray* fBeamDetMWPCDigiY1 = nullptr;
  TClonesArray* fBeamDetMWPCDigiY2 = nullptr; 
  Float_t fElossSigmaOverElossToF = 0.; 
  Float_t fElossSigmaToF = 0.; 
  Float_t fTimeSigmaToF = 0.; 
  Float_t fElossSigmaMWPC = 0.; 
  Float_t fTimeSigmaMWPC = 0.; 
  Float_t fToFElossThreshold = 0.; 
  Float_t fMWPCElossThreshold = 0.; 
  Bool_t  fSigmaEOverEToFIsSet = kFALSE; 
private:
  ClassDef(ERBeamDetDigitizer,1)
};

#endif