PndEmcWaveformToDigi Class Reference

Takes list of PndEmcWaveform and creates PndEmcDigi. More...

#include <PndEmcWaveformToDigi.h>

Inheritance diagram for PndEmcWaveformToDigi:

Public Member Functions
	PndEmcWaveformToDigi (Int_t verbose=0, Bool_t storedigis=kTRUE)
virtual	~PndEmcWaveformToDigi ()
virtual InitStatus	Init ()
	Init Task. More...
virtual void	Exec (Option_t *opt)
	Runs the task. More...
void	SetStorageOfData (Bool_t val)
void	RunTimeBased ()
void	UseDigitizationVersion2 ()
void	SetFakeOnline (bool d=true)
virtual void	SetPSAAlgorithm (PndEmcAbsPSA *psa)
	Set PSA Algorithm to be used for Barrel and Backward Endcap. More...
virtual void	SetPSAAlgorithmFWD (PndEmcAbsPSA *psa)
	Set PSA Algorithm to be used for Forward Endcap. More...
virtual void	SetPSAAlgorithmPMT (PndEmcAbsPSA *psa)
	Set PSA Algorithm to be used for Forward Shashlik. More...
virtual void	SetCrystalCalibrator (PndEmcAbsCrystalCalibrator *Cal)
	Set The Crystal Calibrator to use. More...
virtual void	FinishTask ()
	Called at end of task. More...
void	SetPersistency (Bool_t val=kTRUE)
Bool_t	GetPersistency ()

Protected Member Functions
virtual void	SetParContainers ()

Private Member Functions
	PndEmcWaveformToDigi (const PndEmcWaveformToDigi &)
PndEmcWaveformToDigi &	operator= (const PndEmcWaveformToDigi &)
Float_t	SmearFakeOnline (Float_t energy)
	ClassDef (PndEmcWaveformToDigi, 1)

Private Attributes
TClonesArray *	fWaveformArray
TClonesArray *	fEvtHeaderArray
TClonesArray *	fDigiArray
Double_t	fSampleRate
Double_t	fSampleRate_PMT
Double_t	fSampleRate_FWD
Double_t	fEnergyDigiThreshold
Double_t	fASIC_Shaping_int_time
Double_t	fPMT_Shaping_int_time
Double_t	fPMT_Shaping_diff_time
Double_t	fCrystal_time_constant
Double_t	fShashlyk_time_constant
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
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
Int_t	fVerbose
Bool_t	fTimeOrderedDigi
	set to kTRUE to use the time ordering of the output data. More...
Double_t	fWfNormalisation
Double_t	fWfNormalisation_fwd
Double_t	fWfNormalisation_pmt
BinaryFunctor *	fFunctor
TVectorD	CoeffMod3
TVectorD	CoeffMod5
TVectorD	CoeffModo
Int_t	fEventNo
Double_t	fTimeWindow
Double_t	fTimeShift
Int_t	totDigisAboveThreshold
Int_t	totHits
Int_t	totExpHits
Int_t	totNumOfWave
Bool_t	fDigitizationVersion2
bool	fFakeOnline

Detailed Description

Takes list of PndEmcWaveform and creates PndEmcDigi.

Definition at line 50 of file PndEmcWaveformToDigi.h.

Constructor & Destructor Documentation

PndEmcWaveformToDigi::PndEmcWaveformToDigi	(	Int_t	verbose = 0,
		Bool_t	storedigis = kTRUE
	)

Definition at line 50 of file PndEmcWaveformToDigi.cxx.

References fDigiPosMethod, fEmcDigiRescaleFactor, and PndPersistencyTask::SetPersistency().

  : fWaveformArray(0)
  , fDigiArray(0)
  , fSampleRate(0)
  , fSampleRate_PMT(0)
  , fEnergyDigiThreshold(0)
  ,     fASIC_Shaping_int_time(0)
  , fPMT_Shaping_int_time(0)
  , fPMT_Shaping_diff_time(0)
  ,     fCrystal_time_constant(0)
  , fShashlyk_time_constant(0)
  , fNumber_of_samples_in_waveform(0)
  , fNumber_of_samples_in_waveform_pmt(0)
  , fDigiPosMethod(0)
  , fEmcDigiRescaleFactor(0)
  ,     fEmcDigiPositionDepthPWO(0)
  , fEmcDigiPositionDepthShashlyk(0)
  , fPulseshape(0)
  ,     fPulseshape_pmt(0)
  , fpsaAlgorithm(0)
  , fpsaAlgorithm_pmt(0)
  , fpsaAlgorithm_fwd(0)
  , fCalibrator()
  , fDigiPar(new PndEmcDigiPar())
  ,     fRecoPar(new PndEmcRecoPar())
  , fFpgaPar(new PndEmcFpgaPar())
  , fVerbose(verbose)
  , fTimeOrderedDigi(kFALSE)
  , fWfNormalisation(0)
  , fWfNormalisation_fwd(0)
  , fWfNormalisation_pmt(0)
  , fDigitizationVersion2(kFALSE)
  , fFakeOnline(false)
{
        fDigiPosMethod="depth";// "surface" or "depth"
        fEmcDigiRescaleFactor=1.08;
        SetPersistency(storedigis);
        //fPndEmcDigiPositionDepth=6.2;
}

PndEmcWaveformToDigi::~PndEmcWaveformToDigi ( )

virtual

Definition at line 93 of file PndEmcWaveformToDigi.cxx.

References fPulseshape, fPulseshape_fwd, and fPulseshape_pmt.

{
        delete fPulseshape;
        delete fPulseshape_pmt;
        delete fPulseshape_fwd;

}

PndEmcWaveformToDigi::PndEmcWaveformToDigi ( const PndEmcWaveformToDigi &  )

private

Member Function Documentation

PndEmcWaveformToDigi::ClassDef ( PndEmcWaveformToDigi  , 1    )

private

void PndEmcWaveformToDigi::Exec ( Option_t *  opt )

virtual

Runs the task.

The task loops over the waveforms and uses the pulse shape analyser (PndEmcAbsPSA) to extract signal height and timing. From this the PndEmcDigi are created.

Parameters

opt

unused

Returns

void

Definition at line 334 of file PndEmcWaveformToDigi.cxx.

References At, PndEmcAbsCrystalCalibrator::Calibrate(), ctime, Double_t, fCalibrator, fDigiArray, PndEmcDigi::fDigiArrayTBD, fEnergyDigiThreshold, fEventNo, PndEmcDigi::fEvtNo, fFakeOnline, fFunctor, fpsaAlgorithm, fpsaAlgorithm_fwd, fpsaAlgorithm_pmt, fSampleRate, fSampleRate_FWD, fSampleRate_PMT, fTimeOrderedDigi, fTimeShift, fVerbose, fWaveformArray, PndEmcWaveform::GetDetectorId(), PndEmcDigi::GetEnergy(), PndEmcWaveform::GetEvtList(), PndEmcAbsPSA::GetHit(), PndEmcWaveform::GetHitIndex(), PndEmcWaveform::GetModule(), PndEmcWaveform::GetTrackId(), i, PndEmcAbsCrystalCalibrator::kCALOK, nHits, PndEmcAbsPSA::Process(), rtime, PndEmcDigi::SetEnergy(), SmearFakeOnline(), timer, totDigisAboveThreshold, totExpHits, totHits, and totNumOfWave.

{
        TStopwatch timer;
        if (fVerbose>2){
                timer.Start();
        }

        fDigiArray->Delete();

        Double_t fevtTime = FairRootManager::Instance()->GetEventTime();

        if(fVerbose>0){
                cout<<endl;
                cout<<"=======================PndEmcWaveformToDigi======================="<<endl;
                std::cout<<"Event NO. #"<<fEventNo<<", event time # "<<fevtTime <<std::endl;
        }
        if(fTimeOrderedDigi){
                if(FairRunAna::Instance()->IsTimeStamp()){
                        fWaveformArray->Delete();
                        Double_t time_length = 40.;//99.98%
                        if(fVerbose >0)
                                cout<<"--I-- time-based simulation, read data to later #"<<time_length<<" ns"<<endl;
                        fWaveformArray  = FairRootManager::Instance()->GetData("EmcSortedWaveform"
                                        , fFunctor
                                        , fevtTime + time_length);
                }
                if(fVerbose>0)
                        std::cout<<"fDigiArrayTBD size #"<<PndEmcDigi::fDigiArrayTBD->GetEntriesFast()<<std::endl;
        }
        Int_t nWaveforms = fWaveformArray->GetEntriesFast();
        if(fVerbose>0){
                std::cout<<"fWaveformArray size #"<<nWaveforms<<std::endl;
        }
        //std::cout<<"fEvtHeaderArray size #"<<fEvtHeaderArray->GetEntriesFast()<<std::endl;
        //Int_t evtNo = FairRunAna::Instance()->GetEventHeader()->GetMCEntryNumber();
        //for(Int_t iheader=0;iheader<fEvtHeaderArray->GetEntriesFast();++iheader){
        //      FairEventHeader* evtHeader= (FairEventHeader*) fEvtHeaderArray->At(iheader);
        //      cout<<"===> Event NO. #"<<evtHeader->GetMCEntryNumber()<<", event time#"<<evtHeader->GetEventTime()<<endl;
        //}


        Double_t fdigiEnergy(0.), fdigiTime(0.), theSampleRate, minT(1e10), maxT(-1e10);
        Int_t hitIndex, nHits, detId, trackId, i_digi(0), fMod;
        PndEmcAbsPSA *thePSA(0);

        totNumOfWave += nWaveforms;
        for (Int_t iWaveform=0; iWaveform<nWaveforms; iWaveform++)
        {
                PndEmcWaveform* theWaveform = (PndEmcWaveform*) fWaveformArray->At(iWaveform);
                hitIndex=theWaveform->GetHitIndex();
                detId=theWaveform->GetDetectorId();
                trackId=theWaveform->GetTrackId();
                fMod=theWaveform->GetModule();

                if(theWaveform->GetTimeStamp() < minT ) minT = theWaveform->GetTimeStamp();
                if(theWaveform->GetTimeStamp() > maxT ) maxT = theWaveform->GetTimeStamp();
                //Double_t timeshift; // how maximum is shifted
                if(fMod== 5){
                        thePSA = fpsaAlgorithm_pmt;
                        theSampleRate = fSampleRate_PMT;
                } else if(fMod == 3){
                        thePSA = fpsaAlgorithm_fwd;
                        theSampleRate = fSampleRate_FWD;
                }else
                {
                        thePSA = fpsaAlgorithm;
                        theSampleRate = fSampleRate;
                }
                //cout<<endl;
                const std::vector<Int_t>& evtList = theWaveform->GetEvtList();
                if(fVerbose > 1){
                        cout<<endl;
                        std::cout<<"The wave#"<<iWaveform<<" includes "<<evtList.size()<<" events #";
                        for(size_t i=0; i< evtList.size(); ++i){
                                std::cout<<evtList[i]<<" ";
                        }
                        cout<<endl;
                }

                nHits = thePSA->Process(theWaveform);

                if(fVerbose > 1){
                        std::cout<<"detID#"<<theWaveform->GetDetectorId()<<", hits#"<<nHits<<std::endl;
                }
                totHits += nHits;
                totExpHits += evtList.size();

                for(Int_t i = 0 ; i< nHits; i++)
                {
                        thePSA->GetHit(i, fdigiEnergy, fTimeShift);
                        if(fCalibrator->Calibrate(fdigiEnergy,detId)!=PndEmcAbsCrystalCalibrator::kCALOK){
                                continue;
                        }
                        fdigiTime = theWaveform->GetTimeStamp();// + fTimeShift/theSampleRate*1.e9;//translate position to global time //TS: I do not unerstand this line waveform time is in ns and the tof of the particle but TimeShift is divided by samplingtime??

                        if(fVerbose>1){
                                //if(fdigiTime-fevtTime>200. || fdigiTime - fevtTime < 0.)
                                {
                                        cout<<"wave #"<<iWaveform<<", WaveformTime#" << theWaveform->GetTimeStamp() << ", mod#"<<fMod<<", theSampleRate#"<<theSampleRate<<", shift#"<<fTimeShift<<endl;
                                        cout<<"fdigiEnergy#"<<fdigiEnergy<<", fdigiTime#"<<fdigiTime<<", evtTime#"<<fevtTime<<", diffT#"<<(fdigiTime-fevtTime)<<endl;
                                }
                        }
      // [R.K.] Assume there is either time based or non-time based, so smearung runs only once
      if(fFakeOnline) fdigiEnergy = SmearFakeOnline(fdigiEnergy);

                        if (fdigiEnergy>fEnergyDigiThreshold)
                        {
                                ++totDigisAboveThreshold;

                                PndEmcDigi* myDigi(0);
                                //std::cout<<"PndEmcWaveformToDigi::trackId #"<<trackId<<endl;
                                myDigi = new((*fDigiArray)[i_digi++]) PndEmcDigi(trackId,detId, fdigiEnergy, fdigiTime, hitIndex);
                                myDigi->fEvtNo = evtList[i];//fEventNo;

                                myDigi->AddLink(FairLink("EmcWaveform", iWaveform));
                        }
                }
        }
        Int_t nOldSize = fDigiArray->GetEntriesFast();
        Int_t nDigi =  PndEmcDigi::fDigiArrayTBD->GetEntriesFast();
        if(fTimeOrderedDigi){
                for(Int_t iDigi = 0; iDigi < nDigi; ++iDigi)
                { // [R.K.] Assume there is either time based or non-time based, so smearung runs only once
                        PndEmcDigi* copy = (PndEmcDigi*)PndEmcDigi::fDigiArrayTBD->At(iDigi);
      if(fFakeOnline) copy->SetEnergy(SmearFakeOnline(copy->GetEnergy()));
                        new((*fDigiArray)[i_digi++]) PndEmcDigi(*copy);
                }
                //free buffer
                PndEmcDigi::fDigiArrayTBD->Delete();
        }
        if(fVerbose >0){
                Int_t nNewSize = fDigiArray->GetEntriesFast();
                std::cout<<nNewSize<<" digis, "<<nOldSize<<" <--waveform, "<<nDigi<<" <--previous event."<<std::endl;
                cout<<"=======================PndEmcWaveformToDigi======================="<<endl;
        }

        //if(fVerbose >1){
        //      cout<<"=======================PndEmcWaveformToDigi======================="<<endl;
        //      cout<<"time region #"<<minT<<", "<<maxT<<endl;
        //      cout<<"=================================================================="<<endl;
        //}

        if (fVerbose>2){
                timer.Stop();
                Double_t rtime = timer.RealTime();
                Double_t ctime = timer.CpuTime();
                cout << "PndEmcWaveformToDigi, Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
        }
        fEventNo ++;
}

void PndEmcWaveformToDigi::FinishTask ( )

virtual

Called at end of task.

Outputs some statistics.

Returns

void

Definition at line 517 of file PndEmcWaveformToDigi.cxx.

References fEnergyDigiThreshold, totDigisAboveThreshold, totExpHits, totHits, and totNumOfWave.

{

        //unmatchCount += fBufferofDigisToBeDetermined.size();
        std::cout<<"========================================================="<<std::endl;
        std::cout<<"PndEmcWaveformToDigi::FinishTask"<<std::endl;
        std::cout<<"========================================================="<<std::endl;
        std::cout<<"Total waveforms# "<<totNumOfWave<<std::endl;
        std::cout<<"Total expected hits# "<<totExpHits<<std::endl;
        std::cout<<"Total fpga hits# "<<totHits<<std::endl;
        std::cout<<"Total digis Above th. #"<<totDigisAboveThreshold<<", threshold# "<<fEnergyDigiThreshold<<std::endl;
        std::cout<<"========================================================="<<std::endl;
        //fpsaAlgorithm->Print();
        //fpsaAlgorithm_pmt->Print();
        //std::cout<<"pileup events #"<<PileupEventSet.size()<<std::endl;
        //std::cout<<"Correct Ana#"<<CorrectAna<<std::endl;
        //std::cout<<"UnCorrect Ana#"<<unCorrectAna<<std::endl;
        //std::set<Int_t>::iterator it = PileupEventSet.begin();
        //std::set<Int_t>::iterator end = PileupEventSet.end();
        //Int_t LineNo= 0;
        //for(; it != end; ++it){
        //      std::cout<<(*it)<<' ';
        //      if(++LineNo % 10 == 0)
        //              std::cout<<std::endl;
        //}
        //std::cout<<std::endl;

}

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

Referenced by PndLmdPixelHitProducerFast::GetPersistance(), PndMdtDigitization::Init(), PndMdtHitProducerIdeal::Init(), PndMdtClusterTask::Init(), PndFtsHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndRichHitProducer::Init(), PndSttHitProducerRealFast::Init(), PndSttHelixHitProducer::Init(), PndDiscTaskPID::Init(), PndIdealTrackFinder::Init(), PndSttMvdGemTracking::Init(), PndMdtTrkProducer::Init(), PndFtsHitProducerRealFull::Init(), PndLmdPixelClusterTask::Init(), PndSttHitProducerRealFull::Init(), PndLmdStripClusterTask::Init(), PndEmcApdHitProducer::Init(), PndMissingPzCleanerTask::Init(), PndEmcMakeRecoHit::Init(), PndEmcMakeClusterOnline::Init(), PndTrackSmearTask::Init(), PndEmcFWEndcapTimebasedWaveforms::Init(), PndSttHitProducerIdeal::Init(), PndEmcFWEndcapDigi::Init(), PndFtsHitProducerIdeal::Init(), PndEmcMakeCluster::Init(), PndMdtPointsToWaveform::Init(), PndDiscTaskDigitization::Init(), PndEmcMakeDigi::Init(), PndSdsTimeWalkCorrTask::Init(), PndLmdPixelHitProducerFast::Init(), PndDrcHitFinder::Init(), PndRichHitFinder::Init(), PndEmcMakeCorr::Init(), PndFtofHitProducerIdeal::Init(), PndEmcHitsToWaveform::Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), PndEmcHitProducer::Init(), PndRecoMultiKalmanTask2::Init(), PndDrcHitProducerReal::Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndDrcTimeDigiTask::Init(), PndRecoKalmanTask2::Init(), PndEmcExpClusterSplitter::Init(), PndSdsNoiseProducer::Init(), PndFtsHoughTrackerTask::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsHybridHitProducer::Init(), PndSdsIdealRecoTask::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndSdsStripHitProducer::Init(), PndGemDigitize::Init(), PndGemFindHits::Init(), PndSdsPixelClusterTask::Init(), PndSdsStripClusterTask::Init(), PndMvdGemTrackFinderOnHits::Init(), PndBarrelTrackFinder::Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

{ return fPersistency; }

InitStatus PndEmcWaveformToDigi::Init ( )

virtual

Init Task.

Prepares the TClonesArrays of PndEmcWaveform for reading and of PndEmcDigi for writing. Also reads the EMC parameters and prepares the pulseshapes (PndEmcAbsPulseshape) and pulse shape analyser (PndEmcAbsPSA).

Returns

InitStatus

Return values

kSUCCESS

success

Definition at line 112 of file PndEmcWaveformToDigi.cxx.

References PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::cf_delay, PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::cf_fit_offset, PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::cf_fitter_length, PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::cf_ratio, PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::clock_unit, PndEmcDigi::depth, Double_t, PndEmcPSAFPGADigitalFilterAnalyser::enable_mwd_filter(), fASIC_Shaping_int_time, fCalibrator, fCrystal_time_constant, fDigiArray, fDigiPar, fDigiPosMethod, fDigitizationVersion2, fEmcDigiPositionDepthPWO, fEmcDigiPositionDepthShashlyk, fEmcDigiRescaleFactor, fEnergyDigiThreshold, fEventNo, fFpgaPar, fFunctor, fFWD_Shaping_int_time, fFWD_time_constant, fNumber_of_samples_in_waveform, fNumber_of_samples_in_waveform_fwd, fNumber_of_samples_in_waveform_pmt, fPMT_Shaping_diff_time, fPMT_Shaping_int_time, fpsaAlgorithm, fpsaAlgorithm_fwd, fpsaAlgorithm_pmt, fPulseshape, fPulseshape_fwd, fPulseshape_pmt, fRecoPar, fSampleRate, fSampleRate_FWD, fSampleRate_PMT, fShashlyk_time_constant, fTimeOrderedDigi, fVerbose, fWaveformArray, fWfNormalisation, fWfNormalisation_fwd, fWfNormalisation_pmt, PndEmcDigiPar::GetASIC_Shaping_int_time(), PndEmcFpgaPar::GetBarrelCFDelayLength(), PndEmcFpgaPar::GetBarrelCFFilterLength(), PndEmcFpgaPar::GetBarrelCFFilterRatio(), PndEmcFpgaPar::GetBarrelCFFitterLength(), PndEmcFpgaPar::GetBarrelMAFilterLength(), PndEmcFpgaPar::GetBarrelMWDFilterLength(), PndEmcFpgaPar::GetBarrelMWDFilterLifeT(), PndEmcFpgaPar::GetBarrelMWDFilterUsed(), PndEmcFpgaPar::GetBarrelPulseThreshold(), PndEmcDigiPar::GetCrystal_time_constant(), PndEmcRecoPar::GetEmcDigiPositionDepthPWO(), PndEmcRecoPar::GetEmcDigiPositionDepthShashlyk(), PndEmcDigiPar::GetEnergyDigiThreshold(), PndEmcFpgaPar::GetForwardCFDelayLength(), PndEmcFpgaPar::GetForwardCFFilterLength(), PndEmcFpgaPar::GetForwardCFFilterRatio(), PndEmcFpgaPar::GetForwardCFFitterLength(), PndEmcFpgaPar::GetForwardMAFilterLength(), PndEmcFpgaPar::GetForwardMWDFilterLength(), PndEmcFpgaPar::GetForwardMWDFilterLifeT(), PndEmcFpgaPar::GetForwardMWDFilterUsed(), PndEmcFpgaPar::GetForwardPulseThreshold(), PndEmcDigiPar::GetFWD_Shaping_int_time(), PndEmcDigiPar::GetFWD_time_constant(), PndEmcAbsPSA::GetHit(), PndEmcDigiPar::GetNumber_of_samples_in_waveform(), PndEmcDigiPar::GetNumber_of_samples_in_waveform_fwd(), PndEmcDigiPar::GetNumber_of_samples_in_waveform_pmt(), PndPersistencyTask::GetPersistency(), PndEmcDigiPar::GetPMT_Shaping_diff_time(), PndEmcDigiPar::GetPMT_Shaping_int_time(), PndEmcDigiPar::GetSampleRate(), PndEmcDigiPar::GetSampleRate_FWD(), PndEmcDigiPar::GetSampleRate_PMT(), PndEmcDigiPar::GetShashlyk_time_constant(), PndEmcFpgaPar::GetShashylikCFDelayLength(), PndEmcFpgaPar::GetShashylikCFFilterLength(), PndEmcFpgaPar::GetShashylikCFFilterRatio(), PndEmcFpgaPar::GetShashylikCFFitterLength(), PndEmcFpgaPar::GetShashylikMAFilterLength(), PndEmcFpgaPar::GetShashylikMWDFilterLength(), PndEmcFpgaPar::GetShashylikMWDFilterLifeT(), PndEmcFpgaPar::GetShashylikMWDFilterUsed(), PndEmcFpgaPar::GetShashylikPulseThreshold(), PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::hit_threshold, PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::iafactor, PndEmcPSAFPGADigitalFilterAnalyser::init(), PndEmcAbsCrystalCalibrator::Init(), PndEmcDigi::InitDigiArrayTBD(), PndEmcStructure::Instance(), PndEmcWaveform::Integral(), PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::ma_trig_M, PndEmcWaveform::Max(), PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::mwd_length, PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams::mwd_tau, PndEmcAbsPSA::Process(), PndEmcDigi::selectDigiPositionMethod(), PndEmcPSAFPGADigitalFilterAnalyser::setBaselineInterval(), PndEmcSimCrystalCalibrator::SetCalibration(), PndEmcHit::SetEnergy(), PndEmcHit::SetTime(), PndEmcDigi::surface, totDigisAboveThreshold, totExpHits, totHits, totNumOfWave, and PndEmcWaveform::UpdateWaveform().

{
        PndEmcDigi::InitDigiArrayTBD();
        // Get RootManager
        FairRootManager* ioman = FairRootManager::Instance();
        if ( ! ioman )
        {
                cout << "-E- PndEmcWaveformToDigi::Init: "
                        << "RootManager not instantiated!" << endl;
                return kFATAL;
        }

        // Get input array
        if(fTimeOrderedDigi){
                fWaveformArray = dynamic_cast<TClonesArray *>( ioman->GetObject("EmcSortedWaveform"));
                if (!FairRunAna::Instance()->IsTimeStamp()) {
                        cout << "-W- PndEmcWaveformToDigi::Init: "
                                  << "Task running timebased, but run not running timebased" << endl;
                }
        }else{
                fWaveformArray =dynamic_cast<TClonesArray *> (ioman->GetObject("EmcWaveform"));
        }
//              fWaveformArrayOrg = (TClonesArray*) ioman->GetObject("EmcWaveform");
        //fEvtHeaderArray = (TClonesArray*) ioman->GetObject("EventHeader.");

        if ( ! fWaveformArray ) {
                cout << "-W- PndEmcWaveformToDigi::Init: "
                        << "No PndEmcWaveform array!" << endl;
                return kERROR;
        }

        // Create and register output array
        fDigiArray = new TClonesArray("PndEmcDigi");
        ioman->Register("EmcDigi","Emc",fDigiArray,GetPersistency());


        fSampleRate=fDigiPar->GetSampleRate();
        fSampleRate_PMT=fDigiPar->GetSampleRate_PMT();
        fSampleRate_FWD = fDigiPar->GetSampleRate_FWD();
        fASIC_Shaping_int_time=fDigiPar->GetASIC_Shaping_int_time();      //s
        fPMT_Shaping_int_time=fDigiPar->GetPMT_Shaping_int_time();      //s
        fPMT_Shaping_diff_time=fDigiPar->GetPMT_Shaping_diff_time();      //s
        fCrystal_time_constant=fDigiPar->GetCrystal_time_constant();  //s
        fShashlyk_time_constant=fDigiPar->GetShashlyk_time_constant();  //s
        fNumber_of_samples_in_waveform=fDigiPar->GetNumber_of_samples_in_waveform();
        fNumber_of_samples_in_waveform_pmt=fDigiPar->GetNumber_of_samples_in_waveform_pmt();
        fEnergyDigiThreshold=fDigiPar->GetEnergyDigiThreshold();
        fEmcDigiPositionDepthPWO=fRecoPar->GetEmcDigiPositionDepthPWO();
        fEmcDigiPositionDepthShashlyk=fRecoPar->GetEmcDigiPositionDepthShashlyk();
        fNumber_of_samples_in_waveform_fwd=fDigiPar->GetNumber_of_samples_in_waveform_fwd();
        fFWD_Shaping_int_time=fDigiPar->GetFWD_Shaping_int_time();
        fFWD_time_constant=fDigiPar->GetFWD_time_constant();

        if(fVerbose>2){
                cout<<"fEnergyDigiThreshold: "<<fEnergyDigiThreshold<<endl;
                cout<<"fEmcDigiPositionDepthPWO: "<<fEmcDigiPositionDepthPWO<<endl;
                cout<<"fEmcDigiPositionDepthShashlyk: "<<fEmcDigiPositionDepthShashlyk<<endl;
        }

        if (!fDigiPosMethod.CompareTo("surface"))
        {
                PndEmcDigi::selectDigiPositionMethod( PndEmcDigi::surface, 0., 0., 1.0 );
        }
        else if (!fDigiPosMethod.CompareTo("depth"))
        {
                PndEmcDigi::selectDigiPositionMethod( PndEmcDigi::depth,
                                fEmcDigiPositionDepthPWO, fEmcDigiPositionDepthShashlyk, fEmcDigiRescaleFactor);
        }
        else
        {
                cout << "-W- PndEmcWaveformToDigi::Init: "
                        << "Unknown digi position method!" << endl;
                return kERROR;
        }

        fPulseshape= new PndEmcAsicPulseshape(fASIC_Shaping_int_time, fCrystal_time_constant);
        fPulseshape_pmt= new PndEmcCRRCPulseshape(fPMT_Shaping_int_time,fPMT_Shaping_diff_time,fShashlyk_time_constant);
        fPulseshape_fwd =new PndEmcAsicPulseshape(fFWD_Shaping_int_time, fFWD_time_constant);//LNP raw signal, decay time constant 25 microseconds

        // Determine normalisation constant for PndEmcWaveform
        PndEmcWaveform *tmpwaveform1=new PndEmcWaveform(0,101010001, fSampleRate, fNumber_of_samples_in_waveform);
        PndEmcWaveform *tmpwaveform2=new PndEmcWaveform(0,101010001, fSampleRate_PMT, fNumber_of_samples_in_waveform_pmt);
        PndEmcWaveform *tmpwaveform3=new PndEmcWaveform(0,101010001, fSampleRate_FWD, fNumber_of_samples_in_waveform_fwd);//length 10

        PndEmcHit *gevHit=new PndEmcHit();
        gevHit->SetEnergy(1.0);
        gevHit->SetTime(0.);
        tmpwaveform1->UpdateWaveform(gevHit, 0, false, 1., 0., fSampleRate, fPulseshape);
        tmpwaveform2->UpdateWaveform(gevHit, 0, false, 1., 0., fSampleRate_PMT, fPulseshape_pmt);
        tmpwaveform3->UpdateWaveform(gevHit, 0, false, 1., 0., fSampleRate_FWD, fPulseshape_fwd);

        fWfNormalisation = tmpwaveform1->Max();
        fWfNormalisation_pmt = tmpwaveform2->Max();
        fWfNormalisation_fwd = tmpwaveform3->Max();

        fFpgaPar->printParams();
        if(fVerbose > 2){
                cout<<"fWfNormalisation#"<<fWfNormalisation<<", I/A ="<<(tmpwaveform1->Integral()/fWfNormalisation)<<endl;
                cout<<"fWfNormalisation#"<<fWfNormalisation<<", I ="<<tmpwaveform1->Integral()<<endl;
                cout<<"fWfNormalisation_pmt#"<<fWfNormalisation_pmt<<", I/A ="<<(tmpwaveform2->Integral()/fWfNormalisation_pmt)<<endl;
                cout<<"fWfNormalisation_fwd#"<<fWfNormalisation_fwd<<", I ="<<(tmpwaveform3->Integral())<<endl;
                cout<<"fWfNormalisation_fwd#"<<fWfNormalisation_fwd<<", I/A ="<<(tmpwaveform3->Integral()/fWfNormalisation_fwd)<<endl;
        }
        //setup barrel backward
        PndEmcPSAFPGADigitalFilterAnalyser* theAlg2 = new PndEmcPSAFPGADigitalFilterAnalyser("BarBckFilter");
        PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams newParams;
        newParams.ma_trig_M = fFpgaPar->GetBarrelMAFilterLength();
        newParams.cf_delay = fFpgaPar->GetBarrelCFFilterLength();
        newParams.cf_ratio = fFpgaPar->GetBarrelCFFilterRatio();
        newParams.cf_fitter_length = fFpgaPar->GetBarrelCFFitterLength();
        newParams.cf_fit_offset = fFpgaPar->GetBarrelCFDelayLength();
        newParams.mwd_length = fFpgaPar->GetBarrelMWDFilterLength();
        newParams.mwd_tau = fFpgaPar->GetBarrelMWDFilterLifeT();
        newParams.hit_threshold = fFpgaPar->GetBarrelPulseThreshold()*fWfNormalisation*0.001;
        newParams.clock_unit = 1./fSampleRate*1.e9;//in unit of ns
        newParams.iafactor = tmpwaveform1->Integral()/fWfNormalisation;//31.4166;
        theAlg2->enable_mwd_filter(fFpgaPar->GetBarrelMWDFilterUsed());
        theAlg2->init(newParams);
        theAlg2->setBaselineInterval(0,0);
        fpsaAlgorithm= theAlg2;
        //setup forward digital filter
        PndEmcPSAFPGADigitalFilterAnalyser* theAlg1 = new PndEmcPSAFPGADigitalFilterAnalyser("FwdFilter");
        //PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams newParams;
        newParams.ma_trig_M = fFpgaPar->GetForwardMAFilterLength();
        newParams.cf_delay = fFpgaPar->GetForwardCFFilterLength();
        newParams.cf_ratio = fFpgaPar->GetForwardCFFilterRatio();
        newParams.cf_fitter_length = fFpgaPar->GetForwardCFFitterLength();
        newParams.cf_fit_offset = fFpgaPar->GetForwardCFDelayLength();
        newParams.mwd_length = fFpgaPar->GetForwardMWDFilterLength();
        newParams.mwd_tau = fFpgaPar->GetForwardMWDFilterLifeT();
        newParams.hit_threshold = fFpgaPar->GetForwardPulseThreshold()*fWfNormalisation_fwd*0.001;//five times larger than noise,
        newParams.clock_unit = 1./fSampleRate_FWD*1.e9;//in unit of ns
        newParams.iafactor = tmpwaveform3->Integral()/fWfNormalisation_fwd;//151.408;
        theAlg1->enable_mwd_filter(fFpgaPar->GetForwardMWDFilterUsed());
        theAlg1->init(newParams);
        theAlg1->setBaselineInterval(0,0);
        fpsaAlgorithm_fwd = theAlg1;

        //setup shashylik
        PndEmcPSAFPGADigitalFilterAnalyser* theAlg3 = new PndEmcPSAFPGADigitalFilterAnalyser("ShashylikFilter");
        //PndEmcPSAFPGASampleAnalyser::SampleAnalyserParams newParams;
        newParams.ma_trig_M = fFpgaPar->GetShashylikMAFilterLength();
        newParams.cf_delay = fFpgaPar->GetShashylikCFFilterLength();
        newParams.cf_ratio = fFpgaPar->GetShashylikCFFilterRatio();
        newParams.cf_fitter_length = fFpgaPar->GetShashylikCFFitterLength();
        newParams.cf_fit_offset = fFpgaPar->GetShashylikCFDelayLength();
        newParams.mwd_length = fFpgaPar->GetShashylikMWDFilterLength();
        newParams.mwd_tau = fFpgaPar->GetShashylikMWDFilterLifeT();
        newParams.hit_threshold = fFpgaPar->GetShashylikPulseThreshold()*fWfNormalisation_pmt*0.001;
        newParams.clock_unit = 1./fSampleRate_PMT*1.e9;//in unit of ns
        newParams.iafactor = tmpwaveform2->Integral()/fWfNormalisation_pmt;//5.08306;//this factor depends on the pulseshape function
        theAlg3->enable_mwd_filter(fFpgaPar->GetShashylikMWDFilterUsed());
        theAlg3->init(newParams);
        theAlg3->setBaselineInterval(0,0);
        fpsaAlgorithm_pmt = theAlg3;

        //Int_t iMapperVersion;
        //PndEmcMapper::Init(iMapperVersion=1);
        PndEmcStructure::Instance();

        fFunctor = new StopTime();
        fEventNo = 0;

        if(fDigitizationVersion2){
                // Matched digital filter
                // Parameters of the filter are hardcoded at the moment
                // For different pulseshape in barrel and endcaps different filters should be implemented
                std::vector<Double_t> params;
                params.push_back(30); // width
                params.push_back(fSampleRate); // Sample rate
                fpsaAlgorithm = new PndEmcPSAMatchedDigiFilter(params,fPulseshape);
                fpsaAlgorithm_fwd = fpsaAlgorithm;
                // Pulse shape analysis algorithm.
                // Simple parabolic fit.
                fpsaAlgorithm_pmt = new PndEmcPSAParabolic();
        }

        if(fCalibrator == NULL){
                PndEmcSimCrystalCalibrator *SimCalibrator = new PndEmcSimCrystalCalibrator();
                // Determine normalisation constant for PndEmcWaveform
                Double_t tmpPeakPosition;
                Double_t WfNormalisation;
                fpsaAlgorithm->Process(tmpwaveform1);
                fpsaAlgorithm->GetHit(0, WfNormalisation, tmpPeakPosition);
                SimCalibrator->SetCalibration(1,WfNormalisation); // Barrel
                SimCalibrator->SetCalibration(2,WfNormalisation); // Barrel
                SimCalibrator->SetCalibration(4,WfNormalisation); // Backward endcap
                fpsaAlgorithm_pmt->Process(tmpwaveform2);
                fpsaAlgorithm_pmt->GetHit(0, WfNormalisation, tmpPeakPosition);
                SimCalibrator->SetCalibration(5,WfNormalisation); // Shashlyk
                fpsaAlgorithm_fwd->Process(tmpwaveform3);
                fpsaAlgorithm_fwd->GetHit(0, WfNormalisation, tmpPeakPosition);
                SimCalibrator->SetCalibration(3,WfNormalisation); // Forward endcap
                fCalibrator=SimCalibrator;
        }
        fCalibrator->Init();


        delete tmpwaveform1;
        delete tmpwaveform2;
        delete tmpwaveform3;
        //delete tmpwaveform4;
        delete gevHit;

        totDigisAboveThreshold = 0;
        totHits = 0;
        totExpHits = 0;
        totNumOfWave = 0;

        cout << "-I- PndEmcWaveformToDigi: Intialization successfull" << endl;
        return kSUCCESS;
}

PndEmcWaveformToDigi& PndEmcWaveformToDigi::operator= ( const PndEmcWaveformToDigi &  )

private

void PndEmcWaveformToDigi::RunTimeBased ( )

inline

Definition at line 62 of file PndEmcWaveformToDigi.h.

References fTimeOrderedDigi.

{fTimeOrderedDigi = kTRUE;}

virtual void PndEmcWaveformToDigi::SetCrystalCalibrator ( PndEmcAbsCrystalCalibrator *  Cal )

inlinevirtual

Set The Crystal Calibrator to use.

Parameters

Cal	The Calibrator to be used

Returns

Definition at line 96 of file PndEmcWaveformToDigi.h.

References fCalibrator.

{fCalibrator=Cal;}

void PndEmcWaveformToDigi::SetFakeOnline ( bool  d = true )

inline

Definition at line 65 of file PndEmcWaveformToDigi.h.

References d, and fFakeOnline.

Referenced by PndMasterDigiTask::PndMasterDigiTask().

{fFakeOnline=d;};

void PndEmcWaveformToDigi::SetParContainers ( )

protectedvirtual

Get parameter containers

Definition at line 485 of file PndEmcWaveformToDigi.cxx.

References fDigiPar, fFpgaPar, fRecoPar, and run.

                                            {

        // Get run and runtime database
        FairRun* run = FairRun::Instance();
        if ( ! run ) Fatal("SetParContainers", "No analysis run");

        FairRuntimeDb* db = run->GetRuntimeDb();
        if ( ! db ) Fatal("SetParContainers", "No runtime database");

        // Get Emc digitisation parameter container
        fDigiPar = (PndEmcDigiPar*) db->getContainer("PndEmcDigiPar");

        // Get Emc reconstruction parameter container
        fRecoPar = (PndEmcRecoPar*) db->getContainer("PndEmcRecoPar");

        // Get Emc FPGA parameter container
        fFpgaPar = (PndEmcFpgaPar*) db->getContainer("PndEmcFpgaPar");
}

void PndPersistencyTask::SetPersistency ( Bool_t  val = kTRUE )

inlineinherited

Definition at line 31 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency, and val.

Referenced by barrelTrackFinder(), digi_complete(), digi_complete_newSTT(), digiOnly_complete(), PndBarrelTrackFinder::PndBarrelTrackFinder(), PndCATracking::PndCATracking(), PndDrcHitFinder::PndDrcHitFinder(), PndEmc2DLocMaxFinder::PndEmc2DLocMaxFinder(), PndEmcExpClusterSplitter::PndEmcExpClusterSplitter(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcFWEndcapDigi::PndEmcFWEndcapDigi(), PndEmcFWEndcapTimebasedWaveforms::PndEmcFWEndcapTimebasedWaveforms(), PndEmcHitProducer::PndEmcHitProducer(), PndEmcHitsToWaveform::PndEmcHitsToWaveform(), PndEmcMakeBump::PndEmcMakeBump(), PndEmcMakeCluster::PndEmcMakeCluster(), PndEmcMakeClusterOnline::PndEmcMakeClusterOnline(), PndEmcMakeDigi::PndEmcMakeDigi(), PndEmcMakeRecoHit::PndEmcMakeRecoHit(), PndEmcMultiWaveformToCalibratedDigi::PndEmcMultiWaveformToCalibratedDigi(), PndEmcPhiBumpSplitter::PndEmcPhiBumpSplitter(), PndEmcTmpWaveformToDigi::PndEmcTmpWaveformToDigi(), PndEmcWaveformToCalibratedDigi::PndEmcWaveformToCalibratedDigi(), PndEmcWaveformToDigi(), PndFtofHitProducerIdeal::PndFtofHitProducerIdeal(), PndFtsCATracking::PndFtsCATracking(), PndFtsHitProducerIdeal::PndFtsHitProducerIdeal(), PndFtsHitProducerRealFast::PndFtsHitProducerRealFast(), PndFtsHitProducerRealFull::PndFtsHitProducerRealFull(), PndFtsHoughTrackerTask::PndFtsHoughTrackerTask(), PndGemDigitize::PndGemDigitize(), PndGemFindHits::PndGemFindHits(), PndIdealTrackFinder::PndIdealTrackFinder(), PndLmdPixelClusterTask::PndLmdPixelClusterTask(), PndLmdPixelHitProducerFast::PndLmdPixelHitProducerFast(), PndMdtClusterTask::PndMdtClusterTask(), PndMdtDigitization::PndMdtDigitization(), PndMdtHitProducerIdeal::PndMdtHitProducerIdeal(), PndMdtPointsToWaveform::PndMdtPointsToWaveform(), PndMdtTrkProducer::PndMdtTrkProducer(), PndMissingPzCleanerTask::PndMissingPzCleanerTask(), PndMvdGemTrackFinderOnHits::PndMvdGemTrackFinderOnHits(), PndMvdHitProducerIdeal::PndMvdHitProducerIdeal(), PndMvdPixelClusterTask::PndMvdPixelClusterTask(), PndMvdTimeWalkCorrTask::PndMvdTimeWalkCorrTask(), PndMvdToPix4ClusterTask::PndMvdToPix4ClusterTask(), PndRecoKalmanTask::PndRecoKalmanTask(), PndRecoKalmanTask2::PndRecoKalmanTask2(), PndRecoMultiKalmanTask::PndRecoMultiKalmanTask(), PndRecoMultiKalmanTask2::PndRecoMultiKalmanTask2(), PndRichHitFinder::PndRichHitFinder(), PndRichHitProducer::PndRichHitProducer(), PndSciTDigiTask::PndSciTDigiTask(), PndSciTHitProducerIdeal::PndSciTHitProducerIdeal(), PndSdsHitProducerIdeal::PndSdsHitProducerIdeal(), PndSdsHybridHitProducer::PndSdsHybridHitProducer(), PndSdsIdealClusterTask::PndSdsIdealClusterTask(), PndSdsIdealRecoTask::PndSdsIdealRecoTask(), PndSdsNoiseProducer::PndSdsNoiseProducer(), PndSdsPixelClusterTask::PndSdsPixelClusterTask(), PndSdsStripClusterTask::PndSdsStripClusterTask(), PndSdsStripHitProducer::PndSdsStripHitProducer(), PndSdsTimeWalkCorrTask::PndSdsTimeWalkCorrTask(), PndSttFindTracks::PndSttFindTracks(), PndSttHelixHitProducer::PndSttHelixHitProducer(), PndSttHitProducerIdeal::PndSttHitProducerIdeal(), PndSttHitProducerRealFast::PndSttHitProducerRealFast(), PndSttHitProducerRealFull::PndSttHitProducerRealFull(), PndSttMatchTracks::PndSttMatchTracks(), PndSttMvdGemTracking::PndSttMvdGemTracking(), PndTrackSmearTask::PndTrackSmearTask(), PndTrkTracking2::PndTrkTracking2(), reco(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), recoideal_complete(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndLmdPixelHitProducerFast::SetPersistance(), PndSdsHitProducerIdeal::SetPersistance(), PndSttMvdGemTracking::SetPersistenc(), PndMdtClusterTask::SetPersistence(), PndSttHelixHitProducer::SetPersistence(), PndMissingPzCleanerTask::SetPersistence(), PndFtsHitProducerRealFast::SetPersistence(), PndFtsHitProducerRealFull::SetPersistence(), PndSttHitProducerIdeal::SetPersistence(), PndSttHitProducerRealFull::SetPersistence(), PndSttHitProducerRealFast::SetPersistence(), PndFtsHitProducerIdeal::SetPersistence(), PndTrackSmearTask::SetPersistence(), PndSciTHitProducerIdeal::SetPersistence(), PndIdealTrackFinder::SetPersistence(), PndSttMatchTracks::SetPersistence(), PndSttFindTracks::SetPersistence(), PndFtsHoughTrackerTask::SetPersistence(), PndTrkTracking2::SetPersistence(), PndEmcMakeRecoHit::SetStorageOfData(), PndEmcMakeClusterOnline::SetStorageOfData(), PndEmcFWEndcapDigi::SetStorageOfData(), PndEmcFWEndcapTimebasedWaveforms::SetStorageOfData(), PndEmcMakeDigi::SetStorageOfData(), PndMdtPointsToWaveform::SetStorageOfData(), PndEmc2DLocMaxFinder::SetStorageOfData(), PndEmcMakeCluster::SetStorageOfData(), PndEmcHitsToWaveform::SetStorageOfData(), PndEmcMakeBump::SetStorageOfData(), PndEmcTmpWaveformToDigi::SetStorageOfData(), SetStorageOfData(), PndEmcWaveformToCalibratedDigi::SetStorageOfData(), PndEmcMultiWaveformToCalibratedDigi::SetStorageOfData(), PndEmcExpClusterSplitter::SetStorageOfData(), PndEmcPhiBumpSplitter::SetStorageOfData(), standard_tracking(), and PndEmcFullDigiTask::StoreDigi().

{ fPersistency = val; }

virtual void PndEmcWaveformToDigi::SetPSAAlgorithm ( PndEmcAbsPSA *  psa )

inlinevirtual

Set PSA Algorithm to be used for Barrel and Backward Endcap.

Parameters

psa	The psa to be used

Definition at line 73 of file PndEmcWaveformToDigi.h.

References fpsaAlgorithm.

{fpsaAlgorithm=psa;}

virtual void PndEmcWaveformToDigi::SetPSAAlgorithmFWD ( PndEmcAbsPSA *  psa )

inlinevirtual

Set PSA Algorithm to be used for Forward Endcap.

Parameters

psa	The psa to be used

Definition at line 80 of file PndEmcWaveformToDigi.h.

References fpsaAlgorithm_fwd.

{fpsaAlgorithm_fwd=psa;}

virtual void PndEmcWaveformToDigi::SetPSAAlgorithmPMT ( PndEmcAbsPSA *  psa )

inlinevirtual

Set PSA Algorithm to be used for Forward Shashlik.

Parameters

psa	The psa to be used

Returns

Definition at line 88 of file PndEmcWaveformToDigi.h.

References fpsaAlgorithm_pmt.

{fpsaAlgorithm_pmt=psa;}

void PndEmcWaveformToDigi::SetStorageOfData

(

Bool_t

val

)

Definition at line 504 of file PndEmcWaveformToDigi.cxx.

References PndPersistencyTask::SetPersistency().

Referenced by digi_sttcombi(), emc(), emc_complete(), emc_complete_corr(), emc_correction_data_production(), emc_correction_QA_data_production(), and runDigiReco().

{
        SetPersistency(val);
        return;
}

Float_t PndEmcWaveformToDigi::SmearFakeOnline ( Float_t  energy )

private

Definition at line 597 of file PndEmcWaveformToDigi.cxx.

References energy, fFakeOnline, and sqrt().

Referenced by Exec().

{ // Additional fake energy resolution to get about a factor 2 worse
  // dE/E =~ 1% + 1.63%/(E/GeV)
  if(!fFakeOnline) return energy;
  Float_t resolution=0.01*energy+0.0163*sqrt(energy);
  return gRandom->Gaus(energy,resolution);
}

void PndEmcWaveformToDigi::UseDigitizationVersion2 ( )

inline

Definition at line 64 of file PndEmcWaveformToDigi.h.

References fDigitizationVersion2.

{ fDigitizationVersion2 = kTRUE;}

Member Data Documentation

TVectorD PndEmcWaveformToDigi::CoeffMod3

private

Definition at line 167 of file PndEmcWaveformToDigi.h.

TVectorD PndEmcWaveformToDigi::CoeffMod5

private

Definition at line 168 of file PndEmcWaveformToDigi.h.

TVectorD PndEmcWaveformToDigi::CoeffModo

private

Definition at line 169 of file PndEmcWaveformToDigi.h.

Double_t PndEmcWaveformToDigi::fASIC_Shaping_int_time

private

Definition at line 124 of file PndEmcWaveformToDigi.h.

Referenced by Init().

PndEmcAbsCrystalCalibrator* PndEmcWaveformToDigi::fCalibrator

private

Definition at line 148 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), Init(), and SetCrystalCalibrator().

Double_t PndEmcWaveformToDigi::fCrystal_time_constant

private

Definition at line 127 of file PndEmcWaveformToDigi.h.

Referenced by Init().

TClonesArray* PndEmcWaveformToDigi::fDigiArray

private

output array of EmcDigis

Definition at line 118 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), and Init().

PndEmcDigiPar* PndEmcWaveformToDigi::fDigiPar

private

Definition at line 150 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and SetParContainers().

TString PndEmcWaveformToDigi::fDigiPosMethod

private

Definition at line 136 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and PndEmcWaveformToDigi().

Bool_t PndEmcWaveformToDigi::fDigitizationVersion2

private

Definition at line 181 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and UseDigitizationVersion2().

Double_t PndEmcWaveformToDigi::fEmcDigiPositionDepthPWO

private

Definition at line 138 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Double_t PndEmcWaveformToDigi::fEmcDigiPositionDepthShashlyk

private

Definition at line 139 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Double_t PndEmcWaveformToDigi::fEmcDigiRescaleFactor

private

Definition at line 137 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and PndEmcWaveformToDigi().

Double_t PndEmcWaveformToDigi::fEnergyDigiThreshold

private

Definition at line 123 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), FinishTask(), and Init().

Int_t PndEmcWaveformToDigi::fEventNo

private

Definition at line 171 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcWaveformToDigi::fEvtHeaderArray

private

Definition at line 115 of file PndEmcWaveformToDigi.h.

bool PndEmcWaveformToDigi::fFakeOnline

private

Definition at line 183 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), SetFakeOnline(), and SmearFakeOnline().

PndEmcFpgaPar* PndEmcWaveformToDigi::fFpgaPar

private

Reconstruction parameter container

Definition at line 152 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and SetParContainers().

BinaryFunctor* PndEmcWaveformToDigi::fFunctor

private

Definition at line 164 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), and Init().

Double_t PndEmcWaveformToDigi::fFWD_Shaping_int_time

private

Definition at line 132 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Double_t PndEmcWaveformToDigi::fFWD_time_constant

private

Definition at line 133 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Int_t PndEmcWaveformToDigi::fNumber_of_samples_in_waveform

private

Definition at line 129 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Int_t PndEmcWaveformToDigi::fNumber_of_samples_in_waveform_fwd

private

Definition at line 131 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Int_t PndEmcWaveformToDigi::fNumber_of_samples_in_waveform_pmt

private

Definition at line 130 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Double_t PndEmcWaveformToDigi::fPMT_Shaping_diff_time

private

Definition at line 126 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Double_t PndEmcWaveformToDigi::fPMT_Shaping_int_time

private

Definition at line 125 of file PndEmcWaveformToDigi.h.

Referenced by Init().

PndEmcAbsPSA* PndEmcWaveformToDigi::fpsaAlgorithm

private

Definition at line 144 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), Init(), and SetPSAAlgorithm().

PndEmcAbsPSA* PndEmcWaveformToDigi::fpsaAlgorithm_fwd

private

Definition at line 146 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), Init(), and SetPSAAlgorithmFWD().

PndEmcAbsPSA* PndEmcWaveformToDigi::fpsaAlgorithm_pmt

private

Definition at line 145 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), Init(), and SetPSAAlgorithmPMT().

PndEmcAbsPulseshape* PndEmcWaveformToDigi::fPulseshape

private

Definition at line 141 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and ~PndEmcWaveformToDigi().

PndEmcAbsPulseshape* PndEmcWaveformToDigi::fPulseshape_fwd

private

Definition at line 143 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and ~PndEmcWaveformToDigi().

PndEmcAbsPulseshape* PndEmcWaveformToDigi::fPulseshape_pmt

private

Definition at line 142 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and ~PndEmcWaveformToDigi().

PndEmcRecoPar* PndEmcWaveformToDigi::fRecoPar

private

Digitisation parameter container

Definition at line 151 of file PndEmcWaveformToDigi.h.

Referenced by Init(), and SetParContainers().

Double_t PndEmcWaveformToDigi::fSampleRate

private

Definition at line 120 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), and Init().

Double_t PndEmcWaveformToDigi::fSampleRate_FWD

private

Definition at line 122 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), and Init().

Double_t PndEmcWaveformToDigi::fSampleRate_PMT

private

Definition at line 121 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), and Init().

Double_t PndEmcWaveformToDigi::fShashlyk_time_constant

private

Definition at line 128 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Bool_t PndEmcWaveformToDigi::fTimeOrderedDigi

private

set to kTRUE to use the time ordering of the output data.

Definition at line 158 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), Init(), and RunTimeBased().

Double_t PndEmcWaveformToDigi::fTimeShift

private

Definition at line 173 of file PndEmcWaveformToDigi.h.

Referenced by Exec().

Double_t PndEmcWaveformToDigi::fTimeWindow

private

Definition at line 172 of file PndEmcWaveformToDigi.h.

Int_t PndEmcWaveformToDigi::fVerbose

private

FPGA parameter container Verbosity level

Definition at line 157 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcWaveformToDigi::fWaveformArray

private

Input array of PndEmcWaveforms

Definition at line 113 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), and Init().

Double_t PndEmcWaveformToDigi::fWfNormalisation

private

Definition at line 160 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Double_t PndEmcWaveformToDigi::fWfNormalisation_fwd

private

Definition at line 161 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Double_t PndEmcWaveformToDigi::fWfNormalisation_pmt

private

Definition at line 162 of file PndEmcWaveformToDigi.h.

Referenced by Init().

Int_t PndEmcWaveformToDigi::totDigisAboveThreshold

private

Definition at line 176 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), FinishTask(), and Init().

Int_t PndEmcWaveformToDigi::totExpHits

private

Definition at line 178 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), FinishTask(), and Init().

Int_t PndEmcWaveformToDigi::totHits

private

Definition at line 177 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), FinishTask(), and Init().

Int_t PndEmcWaveformToDigi::totNumOfWave

private

Definition at line 179 of file PndEmcWaveformToDigi.h.

Referenced by Exec(), FinishTask(), and Init().

---

The documentation for this class was generated from the following files:

• PndEmcWaveformToDigi.h
• PndEmcWaveformToDigi.cxx