PndEmcWaveformToDigi.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

//---------------------------------------------------------------------- // File and Version Information:
//      $Id: Exp $
//
// Description:
//      Class PndEmcWaveformToDigi. Module to take the ADC waveforms and produces digi.
//
//   Software developed for the BaBar Detector at the SLAC B-Factory.
// Adapted for the PANDA experiment at GSI      
//      
// Author List:
//      Phil Strother                  Original Author
// Dima Melnichuk - adaption for PANDA              
//
// Copyright Information:
//      Copyright (C) 1996             Imperial College
//
//----------------------------------------------------------------------
//#pragma once
#ifndef PndEmcWaveformToDigi_H
#define PndEmcWaveformToDigi_H

#include <PndPersistencyTask.h>
#include "PndEmcDigiWriteoutBuffer.h"
#include <string>       
#include "FairTSBufferFunctional.h"
#include <map>
#include <list>
#include <set>
#include "TTree.h"
#include "TFile.h"
#include "TVectorD.h"

class PndEmcHit;
class PndEmcTwoCoordIndex;
class PndEmcWaveform;

class TClonesArray;
class PndEmcDigiPar;
class PndEmcRecoPar;
class PndEmcFpgaPar;
class PndEmcAbsPSA;
class PndEmcAbsPulseshape;
class PndEmcAbsCrystalCalibrator;


class PndEmcWaveformToDigi : public PndPersistencyTask
{
public:
    // Constructors
    PndEmcWaveformToDigi(Int_t verbose=0, Bool_t storedigis=kTRUE);
    // Destructor
    virtual ~PndEmcWaveformToDigi();

    virtual InitStatus Init();
    virtual void Exec(Option_t* opt);

    void SetStorageOfData(Bool_t val); // Method to specify whether digis are stored or not.
    void RunTimeBased(){fTimeOrderedDigi = kTRUE;}

    void UseDigitizationVersion2() { fDigitizationVersion2 = kTRUE;}
  void SetFakeOnline(bool d=true){fFakeOnline=d;};

    virtual void SetPSAAlgorithm(PndEmcAbsPSA *psa){fpsaAlgorithm=psa;}
    virtual void SetPSAAlgorithmFWD(PndEmcAbsPSA *psa){fpsaAlgorithm_fwd=psa;}
    virtual void SetPSAAlgorithmPMT(PndEmcAbsPSA *psa){fpsaAlgorithm_pmt=psa;}
    virtual void SetCrystalCalibrator(PndEmcAbsCrystalCalibrator *Cal){fCalibrator=Cal;}

    virtual void FinishTask();

protected:
    virtual void SetParContainers();
private:
    // don't allow copying (-Weffc++)
    PndEmcWaveformToDigi(const PndEmcWaveformToDigi&);  // no implementation
    PndEmcWaveformToDigi& operator= (const PndEmcWaveformToDigi&);  // no implementation

    //Double_t GetEventTimebyDigiTime(Double_t digiT, Double_t digiE, Int_t detID, bool PrintOut=false) const;
    //Int_t GetIdxByEnergy(Double_t energy) const;

private:
    TClonesArray* fWaveformArray;  
    //TClonesArray* fMcTrackArray;  
    TClonesArray* fEvtHeaderArray;  

    TClonesArray* fDigiArray;  

    Double_t fSampleRate;
    Double_t fSampleRate_PMT;
    Double_t fSampleRate_FWD;
    Double_t fEnergyDigiThreshold;
    Double_t fASIC_Shaping_int_time;      //s
    Double_t fPMT_Shaping_int_time;    //s
    Double_t fPMT_Shaping_diff_time;    //s
    Double_t fCrystal_time_constant;  //s
    Double_t fShashlyk_time_constant;  //s
    Int_t fNumber_of_samples_in_waveform;
    Int_t fNumber_of_samples_in_waveform_pmt;
    Int_t fNumber_of_samples_in_waveform_fwd;
    Double_t  fFWD_Shaping_int_time;
    Double_t    fFWD_time_constant;

    //std::string fDigiPosMethod;// "surface" or "depth"
    TString fDigiPosMethod;
    Double_t fEmcDigiRescaleFactor;
    Double_t fEmcDigiPositionDepthPWO;
    Double_t fEmcDigiPositionDepthShashlyk;

    PndEmcAbsPulseshape* fPulseshape;
    PndEmcAbsPulseshape* fPulseshape_pmt;
    PndEmcAbsPulseshape* fPulseshape_fwd;
    PndEmcAbsPSA *fpsaAlgorithm;
    PndEmcAbsPSA *fpsaAlgorithm_pmt;
    PndEmcAbsPSA *fpsaAlgorithm_fwd;

    PndEmcAbsCrystalCalibrator *fCalibrator;

    PndEmcDigiPar*    fDigiPar;      
    PndEmcRecoPar*    fRecoPar;      
    PndEmcFpgaPar*    fFpgaPar;      
    //PndEmcDigiWriteoutBuffer* fDataBuffer;

    Int_t fVerbose;
    Bool_t fStoreDigis;
    Bool_t fTimeOrderedDigi; 

    Double_t fWfNormalisation; // Waveform normalisation constant
    Double_t fWfNormalisation_fwd; // Waveform normalisation constant
    Double_t fWfNormalisation_pmt;

    BinaryFunctor* fFunctor;

    //coefficients
    TVectorD CoeffMod3;
    TVectorD CoeffMod5;
    TVectorD CoeffModo;

    Int_t fEventNo;
    Double_t fTimeWindow;
    Double_t fTimeShift;
    //buffer for undetermined event

    Int_t totDigisAboveThreshold;
    Int_t totHits ;
    Int_t totExpHits ;
    Int_t totNumOfWave ;

    Bool_t fDigitizationVersion2;
    //
  bool fFakeOnline; // Flag to add ~factor2 of worse resolution .
  Float_t SmearFakeOnline(Float_t energy);

    ClassDef(PndEmcWaveformToDigi,1);
};

#endif