PndEmcHighLowPSA Class Reference

Wrapping class to manage two independent PSAs for high, resp. low gain. More...

#include <PndEmcHighLowPSA.h>

Inheritance diagram for PndEmcHighLowPSA:

Public Member Functions
	PndEmcHighLowPSA (Int_t verbose=0)
virtual void	Init (PndEmcPSAFPGASampleAnalyser *highgainPSA, PndEmcPSAFPGASampleAnalyser *lowgainPSA, Double_t overflowThreshold, Int_t highgainWfIndex=0, Int_t lowgainWfIndex=1)
virtual	~PndEmcHighLowPSA ()
virtual void	SetVerbose (Int_t verbose=0)
virtual Int_t	Process (const PndEmcWaveform *waveform)
	Find Hits in Waveform. More...
virtual void	Reset ()
	reset found hits More...
virtual void	GetHit (Int_t i, Double_t &energy, Double_t &time)
	Get energy and time of hit. More...
virtual Int_t	GetWaveformIdx (Int_t i)
virtual void	Process (const PndEmcWaveform *waveform, Double_t &amplitude, Double_t &time)
	Deprecated. Do not use. More...
	ClassDef (PndEmcAbsPSA, 1)

Private Attributes
PndEmcPSAFPGASampleAnalyser *	fHighgainPSA
PndEmcPSAFPGASampleAnalyser *	fLowgainPSA
Int_t	fIdx_high
Int_t	fIdx_low
std::vector< std::pair < PndEmcPSAFPGASampleAnalyser *, int > >	fHitsInFE
Double_t	fOverflowThreshold
Int_t	fVerbose

Detailed Description

Wrapping class to manage two independent PSAs for high, resp. low gain.

Author

Ph. Mahlberg mahlb.nosp@m.erg@.nosp@m.hiskp.nosp@m..uni.nosp@m.-bonn.nosp@m..de

Definition at line 24 of file PndEmcHighLowPSA.h.

Constructor & Destructor Documentation

PndEmcHighLowPSA::PndEmcHighLowPSA

(

Int_t

verbose = 0

)

Definition at line 9 of file PndEmcHighLowPSA.cxx.

                                                :  fHighgainPSA(NULL), fLowgainPSA(NULL),fVerbose(verbose) {
};

virtual PndEmcHighLowPSA::~PndEmcHighLowPSA ( )

inlinevirtual

Definition at line 29 of file PndEmcHighLowPSA.h.

{}

Member Function Documentation

PndEmcAbsPSA::ClassDef ( PndEmcAbsPSA  , 1    )

inherited

void PndEmcHighLowPSA::GetHit ( Int_t  i, Double_t &  energy, Double_t &  time  )

virtual

Get energy and time of hit.

Parameters

Int_t	Index of Hit
Double_t	Energy of Hit
Double_t	Time of Hit

Returns

Implements PndEmcAbsPSA.

Definition at line 145 of file PndEmcHighLowPSA.cxx.

References fHitsInFE, PndEmcPSAFPGASampleAnalyser::GetHit(), and i.

Referenced by PndEmcFWEndcapDigi::Exec().

                                                                       {

        PndEmcPSAFPGASampleAnalyser* psa = fHitsInFE[i].first;
        if(psa!=NULL) {
                psa->GetHit(fHitsInFE[i].second, energy, time);
        } else {
                //std::cerr << "-E PndEmcHighLowPSA::GetHit: No suitable PSA available" << std::cout;
                energy = 0;
                time = 0;
        }

}

Int_t PndEmcHighLowPSA::GetWaveformIdx ( Int_t  i )

virtual

Definition at line 158 of file PndEmcHighLowPSA.cxx.

References fHighgainPSA, fHitsInFE, fIdx_high, fIdx_low, fLowgainPSA, and i.

Referenced by PndEmcFWEndcapDigi::Exec().

                                              {

        if(i<(int)fHitsInFE.size() && i>=0) {
                PndEmcPSAFPGASampleAnalyser* psa = fHitsInFE[i].first;
                if(psa==fHighgainPSA) return fIdx_high;
                else if(psa==fLowgainPSA) return fIdx_low;
        }

        return -1;              //error
}

void PndEmcHighLowPSA::Init ( PndEmcPSAFPGASampleAnalyser *  highgainPSA, PndEmcPSAFPGASampleAnalyser *  lowgainPSA, Double_t  overflowThreshold, Int_t  highgainWfIndex = 0, Int_t  lowgainWfIndex = 1  )

virtual

Definition at line 12 of file PndEmcHighLowPSA.cxx.

References fHighgainPSA, fIdx_high, fIdx_low, fLowgainPSA, and fOverflowThreshold.

Referenced by PndEmcFWEndcapDigi::Init().

                                                                                                                                                                                      {
        fHighgainPSA = highgainPSA;
        fLowgainPSA = lowgainPSA;
        fIdx_high = highgainWfIndex;
        fIdx_low = lowgainWfIndex;
        fOverflowThreshold = overflowThreshold;
}

Int_t PndEmcHighLowPSA::Process ( const PndEmcWaveform *  waveform )

virtual

Find Hits in Waveform.

Parameters

PndEmcWaveform

Returns

Return number of hits, negative if an error occurred

Implements PndEmcAbsPSA.

Definition at line 28 of file PndEmcHighLowPSA.cxx.

References fHighgainPSA, fHitsInFE, fIdx_high, fIdx_low, fLowgainPSA, fOverflowThreshold, fVerbose, PndEmcMultiWaveform::GetActiveWaveform(), PndEmcWaveform::GetSignal(), PndEmcPSAFPGASampleAnalyser::GetStatus(), PndEmcPSAFPGASampleAnalyser::kPileupFinished, PndEmcPSAFPGASampleAnalyser::kPulseDetected, PndEmcPSAFPGASampleAnalyser::kPulseFinished, PndEmcPSAFPGASampleAnalyser::put(), Reset(), and PndEmcMultiWaveform::SetActiveWaveform().

Referenced by PndEmcFWEndcapDigi::Exec().

                                                              {

        Reset();

        PndEmcMultiWaveform* multiWf = const_cast<PndEmcMultiWaveform*>(dynamic_cast<const PndEmcMultiWaveform*>(waveform));

        if(!multiWf) {
                std::cerr << "passed waveform is not of type PndEmcMultiwaveform" << std::endl;
                return -1;
        }

        //active method

        Int_t activeWf = multiWf->GetActiveWaveform();

        multiWf->SetActiveWaveform(fIdx_low);
        std::vector<double> signal_low = waveform->GetSignal();
        std::vector<double>::iterator it_low=signal_low.begin();

        multiWf->SetActiveWaveform(fIdx_high);
        std::vector<double> signal_high = waveform->GetSignal();
        std::vector<double>::iterator it_high=signal_high.begin();

        multiWf->SetActiveWaveform(activeWf);

        bool active_high = false;
        bool active_low = false;

        Int_t counter_low=0;
        Int_t counter_high=0;

        PndEmcPSAFPGASampleAnalyser::status_t status_low;
        PndEmcPSAFPGASampleAnalyser::status_t status_high;

        std::set<Int_t> hits_lowgain;
        bool overflow = false;

        while(it_low!=signal_low.end() && it_high!=signal_high.end()) {

                fHighgainPSA->put(*it_high);
                fLowgainPSA->put(*it_low);
                status_high = fHighgainPSA->GetStatus();
                status_low = fLowgainPSA->GetStatus();

                if(status_high==PndEmcPSAFPGASampleAnalyser::kPulseFinished) {
                        active_high = false;
                        if(!overflow) {
                                fHitsInFE.push_back(std::make_pair(fHighgainPSA, counter_high));
                                if(fVerbose>=2) {
                                        std::cout << "I- PndEmcHighLowPSA: adding highgain hit: #:"<< counter_high << std::endl;
                                }
                        }
                        counter_high++;
                } else if(status_high==PndEmcPSAFPGASampleAnalyser::kPileupFinished) {
                        active_high = false;
                        if(!overflow) {
                                fHitsInFE.push_back(std::make_pair(fHighgainPSA, counter_high));
                                fHitsInFE.push_back(std::make_pair(fHighgainPSA, counter_high+1));
                                if(fVerbose>=2) {
                                        //std::cout << "I- PndEmcHighLowPSA: adding highgain hits #:" << counter_high << "," << counter_high++ << std::endl; //FIXME Implicit operation is bad style!
                                        std::cout << "I- PndEmcHighLowPSA: adding highgain hits #:" << counter_high << ",";
          counter_high++; // [R.K. 9/2018] made implicit operation explicit
          std::cout << counter_high << std::endl;
                                }
                        }
                        counter_high+=2;
                } else if(status_high==PndEmcPSAFPGASampleAnalyser::kPulseDetected) {
                        active_high = true;
                }

                if(active_high && ((*it_high) > fOverflowThreshold) && !overflow) {
                        if(fVerbose>=2) {
                                std::cout << "I- PndEmcHighLowPSA: overflow detected" << std::endl;
                        }
                        overflow = true;
                }

                if(status_low==PndEmcPSAFPGASampleAnalyser::kPulseFinished) {
                        hits_lowgain.insert(counter_low);
                        counter_low++;
                        active_low = false;
                } else if(status_low==PndEmcPSAFPGASampleAnalyser::kPileupFinished) {
                        if(fVerbose>=3) {
                                std::cout << "I- PndEmcHighLowPSA: pileup in lowgain candidate: #" << counter_low << ", #" << counter_low+1 << std::endl;
                }
                        hits_lowgain.insert(counter_low);
                        hits_lowgain.insert(counter_low+1);
                        counter_low+=2;
                        active_low = false;
                } else if(status_low==PndEmcPSAFPGASampleAnalyser::kPulseDetected) {
                        active_low = true;
                }

                if((!active_high) && (!active_low)) {
                        if(overflow) {
                                if(fVerbose>=2) {
                                        std::cout << "I- PndEmcHighLowPSA: adding lowgain hit(s) #:";
                                        for(std::set<Int_t>::iterator it=hits_lowgain.begin(); it!=hits_lowgain.end(); ++it) {
                                                std::cout << *it << ",";
                                        }
                                        std::cout << "\b\n" << std::flush;
                                }
                                for(std::set<Int_t>::iterator it=hits_lowgain.begin(); it!=hits_lowgain.end(); ++it) {
                                        fHitsInFE.push_back(std::make_pair(fLowgainPSA, *it));
                                }
                        }
                        hits_lowgain.clear();
                        overflow = false;
                }
                it_low++;
                it_high++;
        }

        return fHitsInFE.size();
}

void PndEmcAbsPSA::Process ( const PndEmcWaveform *  waveform, Double_t &  amplitude, Double_t &  time  )

virtualinherited

Deprecated. Do not use.

Parameters

PndEmcWaveform
Double_t
Double_t

Returns

Reimplemented in PndEmcPSAMatchedDigiFilter, PndEmcPSAParabolic, PndEmcPSAParabolicBaseline, and PndEmcPSATrapDigiFilter.

Definition at line 2 of file PndEmcAbsPSA.cxx.

References PndEmcAbsPSA::GetHit(), and PndEmcAbsPSA::Process().

                                                                                             {
                if(Process(waveform)>0){
                        GetHit(0,amplitude,time);
                }else{
                        amplitude = 0.;
                        time = 0.;
                }
        }; //deprecated do not use

void PndEmcHighLowPSA::Reset ( )

virtual

reset found hits

Returns

Reimplemented from PndEmcAbsPSA.

Definition at line 20 of file PndEmcHighLowPSA.cxx.

References fHighgainPSA, fHitsInFE, fLowgainPSA, and PndEmcPSAFPGASampleAnalyser::Reset().

Referenced by Process().

                             {
        //fActivePSA=0;
        fHitsInFE.clear();
        fHighgainPSA->Reset();
        fLowgainPSA->Reset();
}

virtual void PndEmcHighLowPSA::SetVerbose ( Int_t  verbose = 0 )

inlinevirtual

Definition at line 31 of file PndEmcHighLowPSA.h.

References fVerbose, and verbose.

{ fVerbose = verbose; };

Member Data Documentation

PndEmcPSAFPGASampleAnalyser* PndEmcHighLowPSA::fHighgainPSA

private

Definition at line 38 of file PndEmcHighLowPSA.h.

Referenced by GetWaveformIdx(), Init(), Process(), and Reset().

std::vector< std::pair<PndEmcPSAFPGASampleAnalyser*, int> > PndEmcHighLowPSA::fHitsInFE

private

Definition at line 44 of file PndEmcHighLowPSA.h.

Referenced by GetHit(), GetWaveformIdx(), Process(), and Reset().

Int_t PndEmcHighLowPSA::fIdx_high

private

Definition at line 41 of file PndEmcHighLowPSA.h.

Referenced by GetWaveformIdx(), Init(), and Process().

Int_t PndEmcHighLowPSA::fIdx_low

private

Definition at line 42 of file PndEmcHighLowPSA.h.

Referenced by GetWaveformIdx(), Init(), and Process().

PndEmcPSAFPGASampleAnalyser* PndEmcHighLowPSA::fLowgainPSA

private

Definition at line 39 of file PndEmcHighLowPSA.h.

Referenced by GetWaveformIdx(), Init(), Process(), and Reset().

Double_t PndEmcHighLowPSA::fOverflowThreshold

private

Definition at line 45 of file PndEmcHighLowPSA.h.

Referenced by Init(), and Process().

Int_t PndEmcHighLowPSA::fVerbose

private

Definition at line 47 of file PndEmcHighLowPSA.h.

Referenced by Process(), and SetVerbose().

---

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

• PndEmcHighLowPSA.h
• PndEmcHighLowPSA.cxx