PndEmcHitsToWaveform Class Reference

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

#include <PndEmcHitsToWaveform.h>

Inheritance diagram for PndEmcHitsToWaveform:

Public Member Functions
	PndEmcHitsToWaveform (Int_t verbose=0, Bool_t storewaves=kTRUE)
virtual	~PndEmcHitsToWaveform ()
virtual InitStatus	Init ()
	Init Task. More...
virtual void	Exec (Option_t *opt)
	Runs the task. More...
PndEmcWaveform *	AddWaveform (Int_t detId, Int_t hitIndex, Int_t numOfSamples, Double_t timeStamp, Double_t sampleRate, Int_t MCTrackID)
	Create a new PndEmcWaveform from the passed parameters. More...
void	SetStorageOfData (Bool_t val)
void	RunTimeBased ()
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
	PndEmcHitsToWaveform (const PndEmcHitsToWaveform &)
PndEmcHitsToWaveform &	operator= (const PndEmcHitsToWaveform &)
	ClassDef (PndEmcHitsToWaveform, 2)

Private Attributes
TClonesArray *	fHitArray
TClonesArray *	fWaveformArray
PndEmcWaveformWriteoutBuffer *	fDataBuffer
Bool_t	fTimeOrderedWaveform
Double_t	fOneBitResolution
Double_t	fOneBitResolutionBW
Double_t	fOneBitResolutionPMT
Double_t	fOneBitResolutionFWD
Int_t	fNBits
Double_t	fDetectedPhotonsPerMeV
Double_t	fDetectedPhotonsPerMeV_PMT
Double_t	fNPhotoElectronsPerMeVAPDBarrel
Double_t	fNPhotoElectronsPerMeVAPDBWD
Double_t	fNPhotoElectronsPerMeVVPT
Double_t	fNPhotoElectronsPerMeVPMT
Double_t	fSensitiveAreaAPD
Double_t	fSensitiveAreaVPT
Double_t	fQuantumEfficiencyAPD
Double_t	fQuantumEfficiencyVPT
Double_t	fQuantumEfficiencyPMT
Double_t	fExcessNoiseFactorAPD
Double_t	fExcessNoiseFactorVPT
Double_t	fExcessNoiseFactorPMT
Double_t	fIncoherent_elec_noise_width_GeV_APD
Double_t	fIncoherent_elec_noise_width_GeV_VPT
Double_t	fEnergyRange
Double_t	fEnergyRangeBW
Double_t	fFirstSamplePhase
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	fASIC_Shaping_int_time
Double_t	fPMT_Shaping_int_time
Double_t	fPMT_Shaping_diff_time
Double_t	fFWD_Shaping_int_time
Double_t	fFWD_time_constant
Double_t	fCrystal_time_constant
Double_t	fShashlyk_time_constant
Double_t	fShashlykSamplingFactor
Double_t	fSampleRate
Double_t	fSampleRate_PMT
Double_t	fSampleRate_FWD
Int_t	fUse_shaped_noise
Int_t	fUse_photon_statistic
Int_t	fNoiseAllChannels
Int_t	fMapVersion
Double_t	fFirstADCBinTime
Double_t	fGevPeakAnalogue
Double_t	fGevPeakAnalogue_PMT
Double_t	fGevPeakAnalogue_FWD
PndEmcDigiPar *	fDigiPar
PndEmcGeoPar *	fGeoPar
Int_t	fVerbose
Int_t	HowManyHit
Int_t	nWaveformProduced
Int_t	HowManyEventPileup
PndEmcAbsPulseshape *	pulseshape1
PndEmcAbsPulseshape *	pulseshape2
PndEmcAbsPulseshape *	pulseshape3

Detailed Description

Takes list of PndEmcHits and creates PndEmcWaveform.

Definition at line 41 of file PndEmcHitsToWaveform.h.

Constructor & Destructor Documentation

PndEmcHitsToWaveform::PndEmcHitsToWaveform	(	Int_t	verbose = 0,
		Bool_t	storewaves = kTRUE
	)

Definition at line 55 of file PndEmcHitsToWaveform.cxx.

References HowManyHit, and PndPersistencyTask::SetPersistency().

                                                                          :
        fHitArray(0), fWaveformArray(0), fDataBuffer(0), fTimeOrderedWaveform(kFALSE), fOneBitResolution(0), fOneBitResolutionBW(0), fOneBitResolutionPMT(0), fNBits(0), fDetectedPhotonsPerMeV(0), fDetectedPhotonsPerMeV_PMT(0), fNPhotoElectronsPerMeVAPDBarrel(0), fNPhotoElectronsPerMeVAPDBWD(0), fNPhotoElectronsPerMeVVPT(0), fNPhotoElectronsPerMeVPMT(0), fSensitiveAreaAPD(0), fSensitiveAreaVPT(0), fQuantumEfficiencyAPD(0), fQuantumEfficiencyVPT(0), fQuantumEfficiencyPMT(0), fExcessNoiseFactorAPD(0), fExcessNoiseFactorVPT(0), fExcessNoiseFactorPMT(0), fIncoherent_elec_noise_width_GeV_APD(0), fIncoherent_elec_noise_width_GeV_VPT(0), fEnergyRange(0), fEnergyRangeBW(0), fFirstSamplePhase(0), fNumber_of_samples_in_waveform(0), fNumber_of_samples_in_waveform_pmt(0), fASIC_Shaping_int_time(0), fPMT_Shaping_int_time(0), fPMT_Shaping_diff_time(0), fCrystal_time_constant(0), fShashlyk_time_constant(0), fShashlykSamplingFactor(0), fSampleRate(0), fSampleRate_PMT(0), fUse_shaped_noise(0), fUse_photon_statistic(0), fNoiseAllChannels(0), fMapVersion(0), fFirstADCBinTime(0), fGevPeakAnalogue(0), fGevPeakAnalogue_PMT(0), fDigiPar(new PndEmcDigiPar()), fGeoPar(new PndEmcGeoPar()), fVerbose(verbose)
{
        HowManyHit = 0;
        SetPersistency(storewaves);
}

PndEmcHitsToWaveform::~PndEmcHitsToWaveform ( )

virtual

Definition at line 66 of file PndEmcHitsToWaveform.cxx.

References pulseshape1, pulseshape2, and pulseshape3.

{
//      delete fw_endcap_noise;
        delete pulseshape1;
        delete pulseshape2;
        delete pulseshape3;
}

PndEmcHitsToWaveform::PndEmcHitsToWaveform ( const PndEmcHitsToWaveform &  )

private

Member Function Documentation

PndEmcWaveform * PndEmcHitsToWaveform::AddWaveform	(	Int_t	detID,
		Int_t	iHit,
		Int_t	numOfSamples,
		Double_t	timeStamp,
		Double_t	sampleRate,
		Int_t	MCTrackID
	)

Create a new PndEmcWaveform from the passed parameters.

Parameters

detID	Detector ID
iHit	Index of PndEmcHit in hit TClonesArray
numOfSamples	Number of samples in waveform
timeStamp	Timestamp from the hit
sampleRate	Sample rate (in Hz)
MCTrackID	MC track ID if available, else -1

Returns

PndEmcWaveform* The created waveform object

Definition at line 594 of file PndEmcHitsToWaveform.cxx.

References fTimeOrderedWaveform, fWaveformArray, and nWaveformProduced.

Referenced by Exec().

{
        nWaveformProduced ++;
        PndEmcWaveform* thisWave(0);
        if(fTimeOrderedWaveform){
                thisWave  = new PndEmcWaveform(MCTrackID,detID,sampleRate, numOfSamples, iHit, timeStamp);
        }else{
                TClonesArray& clref = *fWaveformArray;
                Int_t size = clref.GetEntriesFast();
                thisWave  = new(clref[size]) PndEmcWaveform(MCTrackID,detID,sampleRate, numOfSamples,iHit, timeStamp);
        }
        return thisWave;
}

PndEmcHitsToWaveform::ClassDef ( PndEmcHitsToWaveform  , 2    )

private

void PndEmcHitsToWaveform::Exec ( Option_t *  opt )

virtual

Runs the task.

The tasks consists of using the PndEmcHit to create PndEmcWaveform by using the pulseshape generators (PndEmcAbsPulseshape derived) for the different parts of the EMC. Also the noise is added to the pulseshape here.

Parameters

opt

unused

Returns

void

Definition at line 257 of file PndEmcHitsToWaveform.cxx.

References PndEmcWaveform::AddElecNoiseAndDigitise(), PndEmcWaveform::AddEvt(), PndEmcWaveform::AddShapedElecNoiseAndDigitise(), AddWaveform(), ctime, Double_t, fDataBuffer, fExcessNoiseFactorAPD, fExcessNoiseFactorPMT, fExcessNoiseFactorVPT, fFirstADCBinTime, fFirstSamplePhase, fGevPeakAnalogue, fGevPeakAnalogue_FWD, fGevPeakAnalogue_PMT, fHitArray, fIncoherent_elec_noise_width_GeV_APD, fIncoherent_elec_noise_width_GeV_VPT, fNoiseAllChannels, fNPhotoElectronsPerMeVAPDBarrel, fNPhotoElectronsPerMeVPMT, fNPhotoElectronsPerMeVVPT, fNumber_of_samples_in_waveform, fNumber_of_samples_in_waveform_fwd, fNumber_of_samples_in_waveform_pmt, fOneBitResolution, fOneBitResolutionBW, fOneBitResolutionFWD, fOneBitResolutionPMT, fSampleRate, fSampleRate_FWD, fSampleRate_PMT, fShashlykSamplingFactor, fTimeOrderedWaveform, fUse_photon_statistic, fUse_shaped_noise, fVerbose, fWaveformArray, PndEmcWaveform::GetActiveTime(), PndEmcHit::GetEnergy(), PndEmcHit::GetMcList(), PndEmcHit::GetModule(), PndEmcMapper::GetTciMap(), PndEmcHit::GetTime(), HowManyHit, PndEmcMapper::Instance(), nHits, pulseshape1, pulseshape2, pulseshape3, rtime, PndEmcHit::SetEnergy(), sqrt(), timer, and PndEmcWaveform::UpdateWaveform().

{
        TStopwatch timer;
        if (fVerbose>2){
                timer.Start();
        }
        // Reset output array
        if(fTimeOrderedWaveform) {
                if ( ! fDataBuffer ) Fatal("Exec", "No Waveform Data Buffer");
        }else{
                if ( ! fWaveformArray ) Fatal("Exec", "No Waveform Array");
                fWaveformArray->Delete();
        }

        Double_t EventTime = FairRootManager::Instance()->GetEventTime();//nano seconds
        Int_t nHits = fHitArray->GetEntriesFast();
        Int_t evtNo = FairRun::Instance()->GetEventHeader()->GetMCEntryNumber();
        if (fVerbose>1){
                cout<<"**************************************"<<endl;
                cout<<"Event No. #"<<evtNo<<", EvtTime #"<<EventTime<<std::endl;
                cout<<"PndEmcHitsToWaveform:: Hit array contains "<<nHits<< " hits"<<endl;
                cout<<"fDataBuffer size #"<< (fTimeOrderedWaveform ? fDataBuffer->GetNData() : fWaveformArray->GetEntriesFast()) <<std::endl;
                cout<<"**************************************"<<endl;
        }

        // Variable declaration
        PndEmcHit* theHit = NULL;
        PndEmcWaveform* theWaveform = NULL;
        PndEmcHit* tmpHit = NULL;
        std::set<Int_t> waveformInd;
        Int_t NumOfSamples;

        // Loop over PndEmcHits to add them to correspondent waveforms
        // <set> fWaveformInd contains indexes of detectors for which Waveforms are created
        //PndEmcAsicPulseshape *pulseshape= new PndEmcAsicPulseshape(fASIC_Shaping_int_time,fCrystal_time_constant);
        //PndEmcAbsPulseshape *pulseshape2=new PndEmcCRRCPulseshape(fPMT_Shaping_int_time,fPMT_Shaping_diff_time,fShashlyk_time_constant);
        //PndEmcAbsPulseshape *pulseshape3=new PndEmcAsicPulseshape(10.e-9,25e-6);//fCrystal_time_constant = decay time constant 25 micorseconds

        Double_t sampleRate;
        Double_t TimeMin = 99999.;
        Double_t TimeMax = 0.;
        Double_t WaveformTimeStamp(0.);
        Double_t TimeError(0.);
        Double_t EnergyError(0.);
        Int_t detId, module, MCTrackID;

        HowManyHit += nHits;
        for (Int_t iHit=0; iHit<nHits; iHit++) 
        {
                theHit = (PndEmcHit*) fHitArray->At(iHit);
                module = theHit->GetModule();

                if(module > 5) {
                                std::cout<<" UpdateWaveform: Unknown module number "<<module<<" in EMC digitization. Detector ID = "<<detId<<std::endl;
                                continue;
                }
                detId=theHit->GetDetectorID();
                //if(theHit->GetEnergy() < 0.001) continue;//1MeV
                TimeError = gRandom->Gaus(0, 0.55 + 5.5*TMath::Exp(-27.7*theHit->GetEnergy()));//ns
                EnergyError = sqrt(0.01*0.01+0.02*0.02/theHit->GetEnergy());
                const std::vector<Int_t>& mcTrack = theHit->GetMcList();
                MCTrackID = mcTrack.size() > 0 ? mcTrack[0] : -1;


                waveformInd.insert(detId);
                if(module == 5){
                        sampleRate = fSampleRate_PMT;
                        NumOfSamples = fNumber_of_samples_in_waveform_pmt;
                }else if(module == 3){
                        sampleRate = fSampleRate_FWD;//100 MHz
                        NumOfSamples = fNumber_of_samples_in_waveform_fwd;
                }else{
                        sampleRate = fSampleRate;
                        NumOfSamples = fNumber_of_samples_in_waveform;
                }
                WaveformTimeStamp = EventTime + theHit->GetTime()*1.0e9 ;//to nano seconds
                theWaveform = AddWaveform(detId,iHit,NumOfSamples, WaveformTimeStamp, sampleRate, MCTrackID);
                theWaveform->AddEvt(evtNo);//very important
                theWaveform->SetTimeStampError(TimeError);

                //WaveformMap[detId].push_back(Element(evtNo, theWaveform->GetTimeStamp()));
                if(theWaveform->GetTimeStamp() > TimeMax) TimeMax = theWaveform->GetTimeStamp();
                if(theWaveform->GetTimeStamp() < TimeMin) TimeMin = theWaveform->GetTimeStamp();
                //              Int_t module_wf = theWaveform->GetModule();

                switch (module){
                        case 1: // Barrel
                                theWaveform->UpdateWaveform(theHit, fNPhotoElectronsPerMeVAPDBarrel, fUse_photon_statistic, fExcessNoiseFactorAPD, fFirstADCBinTime, fSampleRate, pulseshape1, EnergyError);
                                if (fUse_shaped_noise==0)
                                {
                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution, EnergyError);
                                }
                                else {
                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution, pulseshape1, fFirstADCBinTime, fSampleRate, EnergyError);
                                }
                                break;
                        case 2: // Barrel
                                theWaveform->UpdateWaveform(theHit, fNPhotoElectronsPerMeVAPDBarrel, fUse_photon_statistic, fExcessNoiseFactorAPD, fFirstADCBinTime, fSampleRate, pulseshape1, EnergyError);
                                if (fUse_shaped_noise==0)
                                {
                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution, EnergyError);
                                }
                                else {
                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution, pulseshape1, fFirstADCBinTime, fSampleRate, EnergyError);
                                }
                                break;
                        case 3: // Fwd endcap
                                theWaveform->UpdateWaveform(theHit, fNPhotoElectronsPerMeVVPT, fUse_photon_statistic, fExcessNoiseFactorVPT, fFirstADCBinTime, fSampleRate_FWD, pulseshape3, EnergyError);//use a different shape
                                if (fUse_shaped_noise==0)
                                {
                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_VPT*fGevPeakAnalogue_FWD,fOneBitResolutionFWD, EnergyError);
                                }
                                else {
                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_VPT*fGevPeakAnalogue_FWD,fOneBitResolutionFWD, pulseshape3, fFirstADCBinTime, fSampleRate_FWD, EnergyError);
                                }
                                break;
                        case 4: // Bwd endcap
                                theWaveform->UpdateWaveform(theHit, fNPhotoElectronsPerMeVVPT, fUse_photon_statistic, fExcessNoiseFactorVPT, fFirstADCBinTime, fSampleRate, pulseshape1, EnergyError);
                                if (fUse_shaped_noise==0)
                                {
                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolutionBW, EnergyError);
                                }
                                else {
                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolutionBW, pulseshape1, fFirstADCBinTime, fSampleRate, EnergyError);
                                }
                                break;
                        case 5: // Shashlyk calorimetr
                                tmpHit = theHit;
                                tmpHit->SetEnergy(tmpHit->GetEnergy()*fShashlykSamplingFactor);
                                EnergyError = sqrt(0.01*0.01+0.02*0.02/tmpHit->GetEnergy());
                                theWaveform->UpdateWaveform(theHit, fNPhotoElectronsPerMeVPMT, fUse_photon_statistic,fExcessNoiseFactorPMT, fFirstADCBinTime, fSampleRate_PMT, pulseshape2, EnergyError);
                                if (fUse_shaped_noise==0)
                                {
                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue_PMT,fOneBitResolutionPMT, EnergyError);
                                }
                                else {
                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue_PMT,fOneBitResolutionPMT, pulseshape2, fFirstADCBinTime, fSampleRate_PMT, EnergyError);
                                }
                                break;
                        default:
                                std::cout<<" UpdateWaveform: Unknown module number "<<module<<" in EMC digitization. Detector ID = "<<detId<<std::endl;
                                //abort();
                }
                if(fTimeOrderedWaveform){
                        fDataBuffer->FillNewData(theWaveform, theWaveform->GetTimeStamp(),  theWaveform->GetActiveTime());
                }
        }
        // Produce waveforms in all the crystals, not only where hits took place 
        // Since it is time consuming, by default it is off
        if (fNoiseAllChannels)
        {
                Int_t detId_tmp, modId_tmp;
                std::map<Int_t,PndEmcTwoCoordIndex*>  intTwoCoordMap =  PndEmcMapper::Instance()->GetTciMap();
                for(std::map<Int_t,PndEmcTwoCoordIndex* >::iterator iter = intTwoCoordMap.begin();
                                iter != intTwoCoordMap.end(); ++iter){
                        detId_tmp=(*iter).first;
                        modId_tmp = detId_tmp/100000000;
                        if(modId_tmp == 5){
                                NumOfSamples = fNumber_of_samples_in_waveform_pmt;
                                sampleRate = fSampleRate_PMT;
                        }else if(modId_tmp == 3){
                                NumOfSamples = fNumber_of_samples_in_waveform_fwd;
                                sampleRate = fSampleRate_FWD;
                        }else{
                                NumOfSamples = fNumber_of_samples_in_waveform;
                                sampleRate = fSampleRate;
                        }

                        if (waveformInd.insert(detId_tmp).second){
                                WaveformTimeStamp = gRandom->Uniform(TimeMin, TimeMax);
                                theWaveform =  AddWaveform(detId_tmp,-1,NumOfSamples, WaveformTimeStamp, sampleRate, -1); // -1 correponds to Waveform produced not from EmcHit but from Noise
                                switch (modId_tmp){
                                        case 1: // Barrel
                                                if (fUse_shaped_noise==0)
                                                {
                                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution);
                                                }
                                                else {
                                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution, pulseshape1, fFirstSamplePhase, fSampleRate);
                                                }
                                                break;
                                        case 2: // Barrel
                                                if (fUse_shaped_noise==0)
                                                {
                                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution);
                                                }
                                                else {
                                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution, pulseshape1, fFirstSamplePhase, fSampleRate);
                                                }
                                                break;
                                        case 3: // Fwd endcap
                                                if (fUse_shaped_noise==0)
                                                {
                                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_VPT*fGevPeakAnalogue_FWD,fOneBitResolutionFWD);
                                                }
                                                else {
                                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_VPT*fGevPeakAnalogue_FWD,fOneBitResolutionFWD, pulseshape3, fFirstSamplePhase, fSampleRate_FWD);
                                                }
                                                break;
                                        case 4: // Bwd endcap
                                                if (fUse_shaped_noise==0)
                                                {
                                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolutionBW);
                                                }
                                                else {
                                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolutionBW, pulseshape1, fFirstSamplePhase, fSampleRate);
                                                }
                                                break;
                                        case 5: // Shashlyk calorimetr
                                                if (fUse_shaped_noise==0)
                                                {
                                                        theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue_PMT,fOneBitResolutionPMT);
                                                }
                                                else {
                                                        theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue_PMT,fOneBitResolutionPMT, pulseshape2, fFirstSamplePhase, fSampleRate_PMT);
                                                }
                                                break;
                                        default:
                                                std::cout<<" UpdateWaveform: Unknown module number "<<modId_tmp<<" in EMC digitization. Detector ID = "<<detId_tmp<<std::endl;
                                                //abort();
                                }
                                if(fTimeOrderedWaveform){
                                        fDataBuffer->FillNewData(theWaveform, theWaveform->GetTimeStamp(),  theWaveform->GetActiveTime());
                                }
                        }
                }
        }

        // Add electronic noise
        // There are two options how to add noise (before and after shaping) 

        /*Int_t nWf = fWaveformArray->GetEntriesFast();
                if (fVerbose>2){
                cout << "Number of waveforms processed= "<<nWf<<endl;
                }

                for (Int_t iWf=0; iWf<nWf; iWf++) {
                theWaveform = (PndEmcWaveform*) fWaveformArray->At(iWf);
                Int_t module = theWaveform->GetModule();
                switch (module){
                case 1: // Barrel 
                if (fUse_shaped_noise==0)
                {
                theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution);
                }
                else {
                theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution, pulseshape, fFirstSamplePhase, fSampleRate);
                }
                break;
                case 2: // Barrel 
                if (fUse_shaped_noise==0)
                {
                theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution);
                }
                else {
                theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolution, pulseshape, fFirstSamplePhase, fSampleRate);
                }
                break;
                case 3: // FWD Endcap 
                if (fUse_shaped_noise==0)
                {
                theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_VPT*fGevPeakAnalogue,fOneBitResolution);
                }
                else {
                theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_VPT*fGevPeakAnalogue,fOneBitResolution, pulseshape, fFirstSamplePhase, fSampleRate);
                }
                break;
                case 4: // BWD Endcap 
                if (fUse_shaped_noise==0)
                {
                theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolutionBW);
                }
                else {
                theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue,fOneBitResolutionBW, pulseshape, fFirstSamplePhase, fSampleRate);
                }
                break;
                case 5: // shashlyk calorimetr 
                if (fUse_shaped_noise==0)
                {
                theWaveform->AddElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue_PMT,fOneBitResolutionPMT);
                }
                else {
                theWaveform->AddShapedElecNoiseAndDigitise(fIncoherent_elec_noise_width_GeV_APD*fGevPeakAnalogue_PMT,fOneBitResolutionPMT, pulseshape2, fFirstSamplePhase, fSampleRate_PMT);
                }
                break;
                default:
                std::cout<<"Add Noise: Unknown module number in EMC digitization"<<std::endl;
        //abort();
        }
        }*/

        //std::cout<<"nWaveformProduced = "<<nWaveformProduced<<std::endl;
        //delete pulseshape;
        //delete pulseshape2;
        //delete pulseshape3;

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

void PndEmcHitsToWaveform::FinishTask

(

)

Called at end of task.

Outputs statistics and writes PndEmcWaveformWriteoutBuffer.

Returns

void

Definition at line 626 of file PndEmcHitsToWaveform.cxx.

References fDataBuffer, fTimeOrderedWaveform, HowManyHit, nWaveformProduced, and PndEmcWaveformWriteoutBuffer::Write().

{
        std::cout<<"==================================================="<<std::endl;
        std::cout<<"PndEmcHitsToWaveform::FinishTask"<<std::endl;
        std::cout<<"***************************************************"<<std::endl;
        std::cout<<"Read Hits# "<<HowManyHit<<std::endl;
        std::cout<<"Produc waveforms# "<<nWaveformProduced<<std::endl;
        std::cout<<"***************************************************"<<std::endl;

        if(fTimeOrderedWaveform){
                fDataBuffer->Write();
        }
}

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(), PndSttHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndRichHitProducer::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(), Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), PndRecoMultiKalmanTask2::Init(), PndEmcHitProducer::Init(), PndDrcHitProducerReal::Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndRecoKalmanTask2::Init(), PndDrcTimeDigiTask::Init(), PndEmcExpClusterSplitter::Init(), PndFtsHoughTrackerTask::Init(), PndSdsNoiseProducer::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsIdealRecoTask::Init(), PndSdsHybridHitProducer::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndGemDigitize::Init(), PndSdsStripHitProducer::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 PndEmcHitsToWaveform::Init ( )

virtual

Init Task.

Prepares the TClonesArrays of PndEmcHit for reading and of PndEmcWaveform for writing. Also reads all relevant EMC parameters and initializes the mapper (PndEmcMapper). Prepares the pulseshapes (PndEmcAbsPulseshape) for the different parts of the EMC and calculates their normalization.

Returns

InitStatus

Return values

kSUCCESS

success

Definition at line 86 of file PndEmcHitsToWaveform.cxx.

References fASIC_Shaping_int_time, fCrystal_time_constant, fDataBuffer, fDetectedPhotonsPerMeV, fDetectedPhotonsPerMeV_PMT, fDigiPar, fEnergyRange, fEnergyRangeBW, fExcessNoiseFactorAPD, fExcessNoiseFactorPMT, fExcessNoiseFactorVPT, fFirstADCBinTime, fFirstSamplePhase, fFWD_Shaping_int_time, fFWD_time_constant, fGeoPar, fGevPeakAnalogue, fGevPeakAnalogue_FWD, fGevPeakAnalogue_PMT, fHitArray, fIncoherent_elec_noise_width_GeV_APD, fIncoherent_elec_noise_width_GeV_VPT, fNBits, fNoiseAllChannels, fNPhotoElectronsPerMeVAPDBarrel, fNPhotoElectronsPerMeVAPDBWD, fNPhotoElectronsPerMeVPMT, fNPhotoElectronsPerMeVVPT, fNumber_of_samples_in_waveform, fNumber_of_samples_in_waveform_fwd, fNumber_of_samples_in_waveform_pmt, fOneBitResolution, fOneBitResolutionBW, fOneBitResolutionFWD, fOneBitResolutionPMT, fPMT_Shaping_diff_time, fPMT_Shaping_int_time, fQuantumEfficiencyAPD, fQuantumEfficiencyPMT, fQuantumEfficiencyVPT, fSampleRate, fSampleRate_FWD, fSampleRate_PMT, fSensitiveAreaAPD, fSensitiveAreaVPT, fShashlyk_time_constant, fShashlykSamplingFactor, fTimeOrderedWaveform, fUse_photon_statistic, fUse_shaped_noise, fVerbose, fWaveformArray, PndEmcDigiPar::GetASIC_Shaping_int_time(), PndEmcDigiPar::GetCrystal_time_constant(), PndEmcDigiPar::GetDetectedPhotonsPerMeV(), PndEmcDigiPar::GetDetectedPhotonsPerMeV_PMT(), PndEmcDigiPar::GetEnergyRange(), PndEmcDigiPar::GetEnergyRangeBW(), PndEmcDigiPar::GetExcessNoiseFactorAPD(), PndEmcDigiPar::GetExcessNoiseFactorPMT(), PndEmcDigiPar::GetExcessNoiseFactorVPT(), PndEmcDigiPar::GetFirstSamplePhase(), PndEmcDigiPar::GetFWD_Shaping_int_time(), PndEmcDigiPar::GetFWD_time_constant(), PndEmcDigiPar::GetIncoherent_elec_noise_width_GeV_APD(), PndEmcDigiPar::GetIncoherent_elec_noise_width_GeV_VPT(), PndEmcGeoPar::GetMapperVersion(), PndEmcDigiPar::GetNBits(), PndEmcDigiPar::GetNoiseAllChannels(), 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::GetQuantumEfficiencyAPD(), PndEmcDigiPar::GetQuantumEfficiencyPMT(), PndEmcDigiPar::GetQuantumEfficiencyVPT(), PndEmcDigiPar::GetSampleRate(), PndEmcDigiPar::GetSampleRate_FWD(), PndEmcDigiPar::GetSampleRate_PMT(), PndEmcWaveform::GetScale(), PndEmcDigiPar::GetSensitiveAreaAPD(), PndEmcDigiPar::GetSensitiveAreaVPT(), PndEmcDigiPar::GetShashlyk_time_constant(), PndEmcDigiPar::GetShashlykSamplingFactor(), PndEmcDigiPar::GetUse_photon_statistic(), PndEmcDigiPar::GetUse_shaped_noise(), PndEmcGeoPar::InitEmcMapper(), PndEmcStructure::Instance(), nWaveformProduced, pulseshape1, pulseshape2, pulseshape3, and PndEmcWaveformWriteoutBuffer::SaveToTree().

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

        // Get input array
        fHitArray = (TClonesArray*) ioman->GetObject("EmcHit");
        if ( ! fHitArray ) {
                cout << "-W- PndEmcHitsToWaveform::Init: "
                        << "No EmcHit array!" << endl;
                return kERROR;
        }
        if(fTimeOrderedWaveform){
                // Create and register output buffer
                fDataBuffer = new PndEmcWaveformWriteoutBuffer("EmcWaveform", "Emc", GetPersistency());
                fDataBuffer->ActivateBuffering(fTimeOrderedWaveform);
                fDataBuffer->SetVerbose(fVerbose);
                ioman->RegisterWriteoutBuffer("EmcWaveform", fDataBuffer);
                fDataBuffer->SaveToTree(kTRUE);
        }else{
                // Create and register output array
                fWaveformArray = new TClonesArray("PndEmcWaveform");
                ioman->Register("EmcWaveform","Emc",fWaveformArray,GetPersistency());
        }

        cout << "-I- PndEmcHitsToWaveform: Intialization successfull" << endl;

        fNBits=fDigiPar->GetNBits();
        fDetectedPhotonsPerMeV=fDigiPar->GetDetectedPhotonsPerMeV();
        fDetectedPhotonsPerMeV_PMT=fDigiPar->GetDetectedPhotonsPerMeV_PMT();
        fSensitiveAreaAPD=fDigiPar->GetSensitiveAreaAPD();
        fSensitiveAreaVPT=fDigiPar->GetSensitiveAreaVPT();
        fQuantumEfficiencyAPD=fDigiPar->GetQuantumEfficiencyAPD();
        fQuantumEfficiencyVPT=fDigiPar->GetQuantumEfficiencyVPT();
        fQuantumEfficiencyPMT=fDigiPar->GetQuantumEfficiencyPMT();
        fExcessNoiseFactorAPD=fDigiPar->GetExcessNoiseFactorAPD();
        fExcessNoiseFactorVPT=fDigiPar->GetExcessNoiseFactorVPT();
        fExcessNoiseFactorPMT = fDigiPar->GetExcessNoiseFactorPMT();
        fIncoherent_elec_noise_width_GeV_APD=fDigiPar->GetIncoherent_elec_noise_width_GeV_APD(); //GeV
        fIncoherent_elec_noise_width_GeV_VPT=fDigiPar->GetIncoherent_elec_noise_width_GeV_VPT(); //GeV
        fEnergyRange=fDigiPar->GetEnergyRange(); //GeV
        fEnergyRangeBW=fDigiPar->GetEnergyRangeBW(); //GeV      
        fFirstSamplePhase=fDigiPar->GetFirstSamplePhase();
        fNumber_of_samples_in_waveform=fDigiPar->GetNumber_of_samples_in_waveform();
        fNumber_of_samples_in_waveform_fwd=fDigiPar->GetNumber_of_samples_in_waveform_fwd();
        fNumber_of_samples_in_waveform_pmt=fDigiPar->GetNumber_of_samples_in_waveform_pmt();
        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
        fFWD_Shaping_int_time=fDigiPar->GetFWD_Shaping_int_time();      //s
        fFWD_time_constant=fDigiPar->GetFWD_time_constant();  //s
        fCrystal_time_constant=fDigiPar->GetCrystal_time_constant();  //s
        fShashlyk_time_constant=fDigiPar->GetShashlyk_time_constant();  //s
        fShashlykSamplingFactor=fDigiPar->GetShashlykSamplingFactor();
        fSampleRate=fDigiPar->GetSampleRate();
        fSampleRate_PMT=fDigiPar->GetSampleRate_PMT();
        fSampleRate_FWD=fDigiPar->GetSampleRate_FWD();
        fUse_shaped_noise=fDigiPar->GetUse_shaped_noise();
        fUse_photon_statistic=fDigiPar->GetUse_photon_statistic();
        fNoiseAllChannels=fDigiPar->GetNoiseAllChannels();

        fDigiPar->printParams();
        // Test how parameters were read from DB.
        cout<<"EMC digitisation parameters "<<endl;
        cout<<"  nBits "<<fNBits<<endl;
        cout<<"  detectedPhotonsPerMeV "<<fDetectedPhotonsPerMeV<<endl;
        cout<<"  detectedPhotonsPerMeV_PMT "<<fDetectedPhotonsPerMeV_PMT<<endl;
        cout<<"  excessNoiseFactor APD"<<fExcessNoiseFactorAPD<<endl;
        cout<<"  excessNoiseFactor VPT"<<fExcessNoiseFactorVPT<<endl;
        cout<<"  excessNoiseFactor PMT"<<fExcessNoiseFactorPMT<<endl;
        cout<<"  incoherent_elec_noise_width_GeV_APD "<<fIncoherent_elec_noise_width_GeV_APD<<endl;
        cout<<"  incoherent_elec_noise_width_GeV_VPT "<<fIncoherent_elec_noise_width_GeV_VPT<<endl;
        cout<<"  energyRange "<<fEnergyRange<<endl;
        cout<<"  energyRangeBW "<<fEnergyRangeBW<<endl; 
        cout<<"  firstSamplePhase "<<fFirstSamplePhase<<endl;
        cout<<"  number_of_samples_in_waveform "<<fNumber_of_samples_in_waveform<<endl;
        cout<<"  number_of_samples_in_waveform_pmt "<<fNumber_of_samples_in_waveform_pmt<<endl;
        cout<<"  number_of_samples_in_waveform_fwd "<<fNumber_of_samples_in_waveform_fwd<<endl;
        cout<<"  ASIC_Shaping_int_time "<<fASIC_Shaping_int_time<<endl;
        cout<<"  PMT_Shaping_int_time "<<fPMT_Shaping_int_time<<endl;
        cout<<"  PMT_Shaping_diff_time "<<fPMT_Shaping_diff_time<<endl;
        cout<<"  crystal_time_constant "<<fCrystal_time_constant<<endl;
        cout<<"  shashlyk_time_constant "<<fShashlyk_time_constant<<endl;
        cout<<"  ShashlykSamplingFactor "<<fShashlykSamplingFactor<<endl;
        cout<<"  sampleRate "<<fSampleRate<<endl;
        cout<<"  sampleRate_PMT "<<fSampleRate_PMT<<endl;
        cout<<"  use_shaped_noise "<<fUse_shaped_noise<<endl;
        cout<<"  use_photon_statistic "<<fUse_photon_statistic<<endl;
        cout<<"  EMC mapper "<<fGeoPar->GetMapperVersion()<<endl;

        fGeoPar->InitEmcMapper();
        PndEmcStructure::Instance();

        //fUse_shaped_noise = 0;
        //fNumber_of_samples_in_waveform = 64;
        //fNumber_of_samples_in_waveform_FWD = 64;
        // Calculate 1 bit resolution (in units of FADC amplitude)
        PndEmcWaveform *tmpwaveform1=new PndEmcWaveform(0,101010001, fNumber_of_samples_in_waveform);
        pulseshape1 =new PndEmcAsicPulseshape(fASIC_Shaping_int_time,fCrystal_time_constant);

        PndEmcWaveform *tmpwaveform2=new PndEmcWaveform(0,101010001, fNumber_of_samples_in_waveform_pmt);
        pulseshape2 =new PndEmcCRRCPulseshape(fPMT_Shaping_int_time,fPMT_Shaping_diff_time,fShashlyk_time_constant);

        PndEmcWaveform *tmpwaveform3=new PndEmcWaveform(0,101010001, fNumber_of_samples_in_waveform_fwd);//length 10
        pulseshape3 =new PndEmcAsicPulseshape(fFWD_Shaping_int_time, fFWD_time_constant);//LNP raw signal, decay time constant 25 microseconds


        //temp parameters;
        //fSampleRate_FWD = 100e6;//100 MHz

        fGevPeakAnalogue     = tmpwaveform1->GetScale(fSampleRate, pulseshape1);
        fGevPeakAnalogue_PMT = tmpwaveform2->GetScale(fSampleRate_PMT, pulseshape2);
        fGevPeakAnalogue_FWD = tmpwaveform3->GetScale(fSampleRate_FWD, pulseshape3);//forward endcap,100 MHz

        //      cout<<" -I- PndEmcHitsToWaveform::Init:=========== "<<endl;
        fOneBitResolution=fEnergyRange/((double) (1<<fNBits))*fGevPeakAnalogue;
        fOneBitResolutionBW=fEnergyRangeBW/((double) (1<<fNBits))*fGevPeakAnalogue;
        fOneBitResolutionPMT=fEnergyRange/((double) (1<<fNBits))*fGevPeakAnalogue_PMT;
        fOneBitResolutionFWD=fEnergyRange/((double) (1<<fNBits))*fGevPeakAnalogue_FWD;

        cout<<"  fGevPeakAnalogue= "<<fGevPeakAnalogue<<endl;
        cout<<"  fGevPeakAnalogue_PMT= "<<fGevPeakAnalogue_PMT<<endl;
        cout<<"  fGevPeakAnalogue_FWD= "<<fGevPeakAnalogue_FWD<<endl;
        cout<<"  fOneBitResolution= "<<fOneBitResolution<<endl;
        cout<<"  fOneBitResolutionBW= "<<fOneBitResolutionBW<<endl;
        cout<<"  fOneBitResolutionPMT= "<<fOneBitResolutionPMT<<endl;

        fFirstADCBinTime=fFirstSamplePhase/fSampleRate;
        //      cout<<"  fFirstADCBinTime= "<<fFirstADCBinTime<<endl;

        // Calculate number of photoelectrons for APD and VPT
        // The number fDetectedPhotonsPerMeV is the measured number of photoelectrons with PM covering the whole rear surface divided by quantum efficiency of PM (18%)
        // To estimate Number of photoelectrons in barrel the rare surface is taken equal for all the crystals 745 mm^2, which is average surface, hovewer it varies depending on the type of the crystal
        // For forward and backward endcap rear surface is equal 26x26=676 mm^2
        // Therefore the different number of photoelectrons are used with APD for barrel and backward endcap
        fNPhotoElectronsPerMeVAPDBarrel=fDetectedPhotonsPerMeV*fSensitiveAreaAPD/745.*fQuantumEfficiencyAPD;
        fNPhotoElectronsPerMeVAPDBWD=fDetectedPhotonsPerMeV*fSensitiveAreaAPD/676.*fQuantumEfficiencyAPD;
        fNPhotoElectronsPerMeVVPT=fDetectedPhotonsPerMeV*fSensitiveAreaVPT/676.*fQuantumEfficiencyVPT;
        fNPhotoElectronsPerMeVPMT=fDetectedPhotonsPerMeV_PMT*fQuantumEfficiencyPMT;
        cout<<"  fNPhotoElectronsPerMeVAPDBarrel= "<<fNPhotoElectronsPerMeVAPDBarrel<<endl;
        cout<<"  fNPhotoElectronsPerMeVAPDBWD= "<<fNPhotoElectronsPerMeVAPDBWD<<endl;
        cout<<"  fNPhotoElectronsPerMeVVPT= "<<fNPhotoElectronsPerMeVVPT<<endl;
        cout<<"  fNPhotoElectronsPerMeVPMT= "<<fNPhotoElectronsPerMeVPMT<<endl;
        //delete pulseshape1;
        delete tmpwaveform1;
        //delete pulseshape2;
        delete tmpwaveform2;
        //delete pulseshape3;
        delete tmpwaveform3;

        nWaveformProduced = 0;

        return kSUCCESS;
}

PndEmcHitsToWaveform& PndEmcHitsToWaveform::operator= ( const PndEmcHitsToWaveform &  )

private

void PndEmcHitsToWaveform::RunTimeBased ( )

inline

Definition at line 57 of file PndEmcHitsToWaveform.h.

References fTimeOrderedWaveform.

{ fTimeOrderedWaveform = kTRUE;}

void PndEmcHitsToWaveform::SetParContainers ( )

protectedvirtual

Get parameter containers

Definition at line 561 of file PndEmcHitsToWaveform.cxx.

References fDigiPar, fGeoPar, 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 geometry parameter container
        fGeoPar = (PndEmcGeoPar*) db->getContainer("PndEmcGeoPar");
        // Get Emc digitisation parameter container
        fDigiPar = (PndEmcDigiPar*) db->getContainer("PndEmcDigiPar");

}

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(), PndEmcMakeBump::PndEmcMakeBump(), PndEmcMakeCluster::PndEmcMakeCluster(), PndEmcMakeClusterOnline::PndEmcMakeClusterOnline(), PndEmcMakeDigi::PndEmcMakeDigi(), PndEmcMakeRecoHit::PndEmcMakeRecoHit(), PndEmcMultiWaveformToCalibratedDigi::PndEmcMultiWaveformToCalibratedDigi(), PndEmcPhiBumpSplitter::PndEmcPhiBumpSplitter(), PndEmcTmpWaveformToDigi::PndEmcTmpWaveformToDigi(), PndEmcWaveformToCalibratedDigi::PndEmcWaveformToCalibratedDigi(), PndEmcWaveformToDigi::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(), PndSttHitProducerRealFull::SetPersistence(), PndSttHitProducerIdeal::SetPersistence(), PndSttHitProducerRealFast::SetPersistence(), PndFtsHitProducerIdeal::SetPersistence(), PndTrackSmearTask::SetPersistence(), PndSciTHitProducerIdeal::SetPersistence(), PndIdealTrackFinder::SetPersistence(), PndSttMatchTracks::SetPersistence(), PndSttFindTracks::SetPersistence(), PndFtsHoughTrackerTask::SetPersistence(), PndTrkTracking2::SetPersistence(), PndEmcMakeRecoHit::SetStorageOfData(), PndEmcFWEndcapDigi::SetStorageOfData(), PndEmcMakeClusterOnline::SetStorageOfData(), PndEmcFWEndcapTimebasedWaveforms::SetStorageOfData(), PndEmcMakeDigi::SetStorageOfData(), PndMdtPointsToWaveform::SetStorageOfData(), PndEmc2DLocMaxFinder::SetStorageOfData(), PndEmcMakeCluster::SetStorageOfData(), SetStorageOfData(), PndEmcMakeBump::SetStorageOfData(), PndEmcTmpWaveformToDigi::SetStorageOfData(), PndEmcWaveformToDigi::SetStorageOfData(), PndEmcWaveformToCalibratedDigi::SetStorageOfData(), PndEmcMultiWaveformToCalibratedDigi::SetStorageOfData(), PndEmcExpClusterSplitter::SetStorageOfData(), PndEmcPhiBumpSplitter::SetStorageOfData(), standard_tracking(), and PndEmcFullDigiTask::StoreDigi().

{ fPersistency = val; }

void PndEmcHitsToWaveform::SetStorageOfData

(

Bool_t

val

)

Definition at line 613 of file PndEmcHitsToWaveform.cxx.

References PndPersistencyTask::SetPersistency().

Referenced by digi(), digi_complete(), digi_complete_runs(), digi_complete_tb(), digi_rich(), digi_sttcombi(), DigiComplete(), digihit(), emc(), emc_complete(), emc_complete_corr(), emc_correction_data_production(), emc_correction_QA_data_production(), error_matrix_data_production(), prod_dig(), run_digi_complete(), run_digi_sttcombi_timebased(), runDigiReco(), and tb_digi_complete().

{
        SetPersistency(val);
        return;
}

Member Data Documentation

Double_t PndEmcHitsToWaveform::fASIC_Shaping_int_time

private

Definition at line 105 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fCrystal_time_constant

private

Definition at line 110 of file PndEmcHitsToWaveform.h.

Referenced by Init().

PndEmcWaveformWriteoutBuffer* PndEmcHitsToWaveform::fDataBuffer

private

Definition at line 74 of file PndEmcHitsToWaveform.h.

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

Double_t PndEmcHitsToWaveform::fDetectedPhotonsPerMeV

private

Definition at line 83 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fDetectedPhotonsPerMeV_PMT

private

Definition at line 84 of file PndEmcHitsToWaveform.h.

Referenced by Init().

PndEmcDigiPar* PndEmcHitsToWaveform::fDigiPar

private

Definition at line 127 of file PndEmcHitsToWaveform.h.

Referenced by Init(), and SetParContainers().

Double_t PndEmcHitsToWaveform::fEnergyRange

private

Definition at line 99 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fEnergyRangeBW

private

Definition at line 100 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fExcessNoiseFactorAPD

private

Definition at line 94 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fExcessNoiseFactorPMT

private

Definition at line 96 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fExcessNoiseFactorVPT

private

Definition at line 95 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fFirstADCBinTime

private

Definition at line 121 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fFirstSamplePhase

private

Definition at line 101 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fFWD_Shaping_int_time

private

Definition at line 108 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fFWD_time_constant

private

Definition at line 109 of file PndEmcHitsToWaveform.h.

Referenced by Init().

PndEmcGeoPar* PndEmcHitsToWaveform::fGeoPar

private

Digitisation parameter container

Definition at line 128 of file PndEmcHitsToWaveform.h.

Referenced by Init(), and SetParContainers().

Double_t PndEmcHitsToWaveform::fGevPeakAnalogue

private

Definition at line 123 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fGevPeakAnalogue_FWD

private

Definition at line 125 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fGevPeakAnalogue_PMT

private

Definition at line 124 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcHitsToWaveform::fHitArray

private

Input array of PndEmcHits

Definition at line 70 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fIncoherent_elec_noise_width_GeV_APD

private

Definition at line 97 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fIncoherent_elec_noise_width_GeV_VPT

private

Definition at line 98 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Int_t PndEmcHitsToWaveform::fMapVersion

private

Definition at line 119 of file PndEmcHitsToWaveform.h.

Int_t PndEmcHitsToWaveform::fNBits

private

Definition at line 82 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Int_t PndEmcHitsToWaveform::fNoiseAllChannels

private

Definition at line 118 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fNPhotoElectronsPerMeVAPDBarrel

private

Definition at line 85 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fNPhotoElectronsPerMeVAPDBWD

private

Definition at line 86 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fNPhotoElectronsPerMeVPMT

private

Definition at line 88 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fNPhotoElectronsPerMeVVPT

private

Definition at line 87 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Int_t PndEmcHitsToWaveform::fNumber_of_samples_in_waveform

private

Definition at line 102 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Int_t PndEmcHitsToWaveform::fNumber_of_samples_in_waveform_fwd

private

Definition at line 104 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Int_t PndEmcHitsToWaveform::fNumber_of_samples_in_waveform_pmt

private

Definition at line 103 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fOneBitResolution

private

Definition at line 77 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fOneBitResolutionBW

private

Definition at line 78 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fOneBitResolutionFWD

private

Definition at line 80 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fOneBitResolutionPMT

private

Definition at line 79 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fPMT_Shaping_diff_time

private

Definition at line 107 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fPMT_Shaping_int_time

private

Definition at line 106 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fQuantumEfficiencyAPD

private

Definition at line 91 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fQuantumEfficiencyPMT

private

Definition at line 93 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fQuantumEfficiencyVPT

private

Definition at line 92 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fSampleRate

private

Definition at line 113 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fSampleRate_FWD

private

Definition at line 115 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fSampleRate_PMT

private

Definition at line 114 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Double_t PndEmcHitsToWaveform::fSensitiveAreaAPD

private

Definition at line 89 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fSensitiveAreaVPT

private

Definition at line 90 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fShashlyk_time_constant

private

Definition at line 111 of file PndEmcHitsToWaveform.h.

Referenced by Init().

Double_t PndEmcHitsToWaveform::fShashlykSamplingFactor

private

Definition at line 112 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Bool_t PndEmcHitsToWaveform::fTimeOrderedWaveform

private

Definition at line 75 of file PndEmcHitsToWaveform.h.

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

Int_t PndEmcHitsToWaveform::fUse_photon_statistic

private

Definition at line 117 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Int_t PndEmcHitsToWaveform::fUse_shaped_noise

private

Definition at line 116 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

Int_t PndEmcHitsToWaveform::fVerbose

private

Geometry parameter container Verbosity level

Definition at line 131 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcHitsToWaveform::fWaveformArray

private

Output array of PndEmcWaveforms

Definition at line 73 of file PndEmcHitsToWaveform.h.

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

Int_t PndEmcHitsToWaveform::HowManyEventPileup

private

Definition at line 136 of file PndEmcHitsToWaveform.h.

Int_t PndEmcHitsToWaveform::HowManyHit

private

Definition at line 134 of file PndEmcHitsToWaveform.h.

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

Int_t PndEmcHitsToWaveform::nWaveformProduced

private

Definition at line 135 of file PndEmcHitsToWaveform.h.

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

PndEmcAbsPulseshape* PndEmcHitsToWaveform::pulseshape1

private

Definition at line 139 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), Init(), and ~PndEmcHitsToWaveform().

PndEmcAbsPulseshape* PndEmcHitsToWaveform::pulseshape2

private

Definition at line 140 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), Init(), and ~PndEmcHitsToWaveform().

PndEmcAbsPulseshape* PndEmcHitsToWaveform::pulseshape3

private

Definition at line 141 of file PndEmcHitsToWaveform.h.

Referenced by Exec(), Init(), and ~PndEmcHitsToWaveform().

---

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

• PndEmcHitsToWaveform.h
• PndEmcHitsToWaveform.cxx