#include <GaussianReceiver.hh>

Public Member Functions
	GaussianReceiver (OptSignal WDMSignal2, RanNumGen RNG2, int DemuxChannArrayIndex2, string InDir, string Job2)
	~GaussianReceiver ()
string	WriteTypeOutput (void)
string	WriteTypeGaussianNoise (void)
string	WriteTypePMDonSignal (void)
string	WriteTypeClockRecovery (void)
void	ReceiveSignal (string FileID)
void	SetFirstOrderPMDParameters (double StokesParametersPlusPlus, cplx StokesParametersPlusMinus, cplx c_plus2, cplx c_minus2)
void	ComputeFirstOrderPMDParameters (OptFiber *PMF)
void	ReceiveSignalDGD_NoisePartiallyPolarized (double DGD2, double ClockRecoveryTime2, double UnitStokesPlusPlus2, double UnitAverageStokesSignal, double DOP_Signal, double TotalNoiseSpectralDensity, double DOP_Noise, double *UnitStokesVectorPolarizedPartNoise, string OutFileID)
Private Member Functions
void	RecoverClockAndSetTimeSlices (void)
void	SetTimeSlices (void)
void	ReceiveSignalUnpolarizedNoise (string FileID)
void	ReceiveSignalNoPMD_NoisePartiallyPolarized (string FileID)
void	ReceiveSignalDGD_NoisePartiallyPolarized (string FileID)
void	ComputeMinimumBERAndDecisionThreshold (string FileID)
void	WriteFilePairVectors (double x, double y, int dim, string OutFileName)
void	WriteFilePDFTimeSlices (string OutFileName)
void	ComputeAveragePDFs (string FileID)
Private Attributes
typeOutput	TypeOutput
typeGaussianNoise	TypeGaussianNoise
typePMDonSignal	TypePMDonSignal
typeClockRecoveryGaussianReceiver	TypeClockRecovery
int	DebugLevel
double	ClockRecoveryTime
OptSignal *	WDMSignal
OptSignal *	SingleChannelSignal
OptSignal *	WDMSignalTmp
OptFilter *	TunableOptFilter
OptDemuxer *	oOptDemuxer
Photodetector *	oPhotodetector
ElecFilter *	oElecFilter
ElecSignalStat *	oElecSignalStat
int	DemuxChannArrayIndex
string	Job
int	ObjectInitialized
int	NumTimeSlices
int	NumSlicesPerBit
double	TimeBetweenSlices
int	NumBitsToPlot
double *	TimeSlices
int *	TimeSliceIndices
double *	VarianceCurrent
double *	MeanCurrent
double	MeanNoise
double	MinCurrent
double	MaxCurrent
double	DeltaCurrent
int	NumPointsCurrent
double *	CurrentVector
double *	BER
double *	AveragePDFZeros
double *	AveragePDFOnes
double **	PDFTimeSlices
double *	PDFTimeSlicesStorageVec
int *	BitString
bool	FirstOrderPMDParametersSet
double	DGD
double *	StokesPlusPlus
cplx *	StokesPlusMinus
cplx	c_plus
cplx	c_minus

Detailed Description

Given a noise free signal, appropriate noise statistics and receive filters this class computes the Gaussian approximation of the PDF of the current in each bit and the corresponding bit-error ratio

Constructor & Destructor Documentation

GaussianReceiver::GaussianReceiver	(	OptSignal *	WDMSignal2,
		RanNumGen *	RNG2,
		int	DemuxChannArrayIndex2,
		string	InDir,
		string	Job2
	)

References AveragePDFOnes, AveragePDFZeros, BER, ClockRecoveryTime, CurrentVector, DebugLevel, DeltaCurrent, DemuxChannArrayIndex, EYE, FIRST_ORDER_PMD, FirstOrderPMDParametersSet, OptSignal::GetBitLength(), OptSignal::GetQtPoints(), OptSignal::GetStringLength(), OptSignal::GetTimeWindow(), Job, LogFile, LogFileSeparator(), LOWER_AND_UPPER, LOWER_ONLY, MANDATORY, MANUALLY_SET, MaxCurrent, MeanCurrent, MinCurrent, NewDouble(), NewInt(), NO_BOUNDS, NumBitsToPlot, NumPointsCurrent, NumSlicesPerBit, NumTimeSlices, ObjectInitialized, oElecFilter, oElecSignalStat, oOptDemuxer, oPhotodetector, OPTIONAL_NO_WARNING, PARTIALLY_POLARIZED, PDFS, PDFTimeSlices, PDFTimeSlicesStorageVec, OptDemuxer::PrepareToDemuxChannel(), ReadDouble(), ReadInt(), ElecSignalStat::SetOptElecFiltersForGaussianStatistics(), SingleChannelSignal, StokesPlusMinus, StokesPlusPlus, TimeBetweenSlices, TimeSliceIndices, TimeSlices, TunableOptFilter, TypeClockRecovery, TypeGaussianNoise, TypeOutput, TypePMDonSignal, VarianceCurrent, WDMSignal, WDMSignalTmp, WriteTypeClockRecovery(), WriteTypeGaussianNoise(), WriteTypeOutput(), and WriteTypePMDonSignal().

Member Function Documentation

void GaussianReceiver::ComputeAveragePDFs ( string FileID ) [private]

void GaussianReceiver::ComputeFirstOrderPMDParameters ( OptFiber * PMF )

The OptFiber PMF must be associated to the OptSignal WDMSignal. You must have called PMF->GetFiberSampleCompleteScatt() before calling this method. The implementation of the theory given in this method assumes that WDMSignal has only been generated using a signal constructor or OptSignal::RegenerateSignalString() and is to be propagated through the PMF. In particular, the signal has not been sent through a PolarizTransformer. Instead, you can use OptSignal::SetInJones() and OptSignal::RegenerateSignalString()

References c_minus, c_plus, DGD, OptSignal::GetCenterFreq(), OptFiber::GetDGD(), OptSignal::GetInJones(), OptFiber::GetOutPD_Vector(), JonesPair2MixedStokes(), SetFirstOrderPMDParameters(), Stokes2Jones(), OptFiber::TransformJonesVectorThroughFiber(), and WDMSignal.

void GaussianReceiver::ComputeMinimumBERAndDecisionThreshold ( string FileID ) [private]

References AveragePDFOnes, AveragePDFZeros, BER, CurrentVector, DebugLevel, Job, max(), NumPointsCurrent, Quadrature(), and WriteFilePairVectors().

Referenced by ReceiveSignal(), and ReceiveSignalDGD_NoisePartiallyPolarized().

void GaussianReceiver::ReceiveSignal ( string FileID )

References ComputeAveragePDFs(), ComputeMinimumBERAndDecisionThreshold(), FIRST_ORDER_PMD, NO_PMD, PARTIALLY_POLARIZED, PDFS, ReceiveSignalDGD_NoisePartiallyPolarized(), ReceiveSignalNoPMD_NoisePartiallyPolarized(), ReceiveSignalUnpolarizedNoise(), TypeGaussianNoise, TypeOutput, TypePMDonSignal, UNPOLARIZED, WriteTypeGaussianNoise(), and WriteTypePMDonSignal().

void GaussianReceiver::ReceiveSignalDGD_NoisePartiallyPolarized	(	double	DGD2,
		double	ClockRecoveryTime2,
		double *	UnitStokesPlusPlus2,
		double *	UnitAverageStokesSignal,
		double	DOP_Signal,
		double	TotalNoiseSpectralDensity,
		double	DOP_Noise,
		double *	UnitStokesVectorPolarizedPartNoise,
		string	OutFileID
	)

Receives a signal that is the sum of (1) a noise-free signal that has gone through first-order PMD fiber (ie PMF) and (2) Partially polarized noise. To be used when know input parameters, say from an experiment. See Jan 5th 2006 "How to validate receiver formula" by J. Zweck

UnitStokesPlusPlus2, UnitAverageStokesSignal and UnitStokesVectorPolarizedPartNoise are 3x1 and have length one.

References ClockRecoveryTime, ComputeAveragePDFs(), ComputeMinimumBERAndDecisionThreshold(), DemuxChannArrayIndex, OptDemuxer::DemuxChannel(), Photodetector::DetectOptSignal(), DGD, ElecFilter::FilterElecSignal(), OptSignal::GetAutocorrelationFunctionOptSignalX(), Photodetector::GetConvFactOptElecPhotod(), ElecSignalStat::GetIntegralForSignalNoiseVariance(), ElecSignalStat::GetMeanNoiseGaussianStatistics(), ElecSignalStat::GetNoiseNoiseVarianceGaussianStatistics(), OptSignal::GetNumChannels(), InnerProduct(), jc, Job, JonesPair2MixedStokes(), MeanCurrent, MeanNoise, NumTimeSlices, oElecFilter, oElecSignalStat, oOptDemuxer, oPhotodetector, PDFS, OptDemuxer::PrepareToDemuxChannel(), OptSignal::RegenerateSignalString(), ElecSignalStat::ResetOptElecFiltersForGaussianStatistics(), OptSignal::SetInJones(), SetTimeSlices(), Photodetector::sfftPM, SingleChannelSignal, sq(), sqrt(), Stokes2Jones(), StokesPlusPlus, Photodetector::TimeShiftSignal(), TimeSliceIndices, TimeSlices, fftPowerMatrices::tPowerSignal, TypeOutput, VarianceCurrent, WDMSignal, and WDMSignalTmp.

Referenced by ReceiveSignal().

void GaussianReceiver::ReceiveSignalDGD_NoisePartiallyPolarized ( string FileID ) [private]

Receives a signal that is the sum of (1) a noise-free signal that has gone through first-order PMD fiber (ie PMF) and (2) Partially polarized noise.

For each different fiber realization and input polarization state, you must call GaussianReceiver::ComputeFirstOrderPMDParameters. You must also propagate the signal through an OptFiber object before calling this method. The clock will be recovered based on the signal after the PMF unless you use TypeClockRecovery = MANUALLY_SET.

References c_minus, c_plus, ClockRecoveryTime, DemuxChannArrayIndex, OptDemuxer::DemuxChannel(), Photodetector::DetectOptSignal(), DGD, ElecFilter::FilterElecSignal(), Photodetector::GetConvFactOptElecPhotod(), ElecSignalStat::GetIntegralForSignalNoiseVariance(), ElecSignalStat::GetMeanNoiseGaussianStatistics(), ElecSignalStat::GetNoiseNoiseVarianceGaussianStatistics(), OptSignal::GetNumChannels(), ReduOptSignal::GetStokesNoise(), ReduOptSignal::GetTotalNoiseSpectralDensity(), OptSignal::GetTypeSimulation(), InnerProduct(), Job, MeanCurrent, MeanNoise, NumTimeSlices, oElecFilter, oElecSignalStat, oOptDemuxer, oPhotodetector, OptSignal::oReduOptSignal, OptDemuxer::PrepareToDemuxChannel(), RecoverClockAndSetTimeSlices(), OptSignal::RegenerateSignalString(), ElecSignalStat::ResetOptElecFiltersForGaussianStatistics(), OptSignal::SetInJones(), Photodetector::sfftPM, SingleChannelSignal, sq(), sqrt(), StokesPlusMinus, StokesPlusPlus, TimeSliceIndices, TimeSlices, fftPowerMatrices::tPowerSignal, VarianceCurrent, WDMSignal, WDMSignalTmp, and Photodetector::WriteFileTimePower().

void GaussianReceiver::ReceiveSignalNoPMD_NoisePartiallyPolarized ( string FileID ) [private]

References sq(), and sqrt().

Referenced by ReceiveSignal().

void GaussianReceiver::ReceiveSignalUnpolarizedNoise ( string FileID ) [private]

void GaussianReceiver::RecoverClockAndSetTimeSlices ( void ) [private]

References AUTOMATIC, ClockRecoveryTime, ElecSignalStat::ComputePerformanceStatistics(), DemuxChannArrayIndex, OptDemuxer::DemuxChannel(), Photodetector::DetectOptSignal(), ElecFilter::FilterElecSignal(), ElecSignalStat::GetClockRecoveryTimeBeforeTimeShift(), oElecFilter, oElecSignalStat, oOptDemuxer, oPhotodetector, SetTimeSlices(), and TypeClockRecovery.

Referenced by ReceiveSignalDGD_NoisePartiallyPolarized(), and ReceiveSignalUnpolarizedNoise().

void GaussianReceiver::SetFirstOrderPMDParameters	(	double *	StokesParametersPlusPlus,
		cplx *	StokesParametersPlusMinus,
		cplx	c_plus2,
		cplx	c_minus2
	)

StokesParametersPlusPlus and StokesParametersPlusMinus are 4x1 vectors for the first-order PMD fiber. c_plus2 and c_minus2 are the cplx splitting factors.

References c_minus, c_plus, FirstOrderPMDParametersSet, StokesPlusMinus, and StokesPlusPlus.

Referenced by ComputeFirstOrderPMDParameters().

void GaussianReceiver::SetTimeSlices ( void ) [private]

References ClockRecoveryTime, DebugLevel, EYE, OptSignal::GetBitLength(), OptSignal::GetDeltaTime(), OptSignal::GetTimeWindow(), NumTimeSlices, SingleChannelSignal, TimeBetweenSlices, TimeSliceIndices, TimeSlices, and TypeOutput.