@@ -88,10 +88,6 @@ class opdet::OpDeconvolutionAlgWiener : opdet::OpDeconvolutionAlg {
8888 short unsigned int fBaseSampleBins ;
8989 double fBaseVarCut ;
9090
91- bool fUseLowPassFilter ;
92- double fFrequencyCutoff ;
93-
94-
9591 // Declare member functions
9692 void ApplyExpoAvSmoothing (std::vector<double >& wf);
9793 void ApplyUnAvSmoothing (std::vector<double >& wf);
@@ -137,8 +133,6 @@ opdet::OpDeconvolutionAlgWiener::OpDeconvolutionAlgWiener(fhicl::ParameterSet co
137133 fElectronics = p.get < std::string >(" Electronics" );
138134 fDaphne_Freq = p.get < double >(" DaphneFreq" );
139135 fScaleHypoSignal = p.get < bool >(" ScaleHypoSignal" );
140- fUseLowPassFilter = p.get < bool >(" UseLowPassFilter" );
141- fFrequencyCutoff = p.get < double >(" FrequencyCutoff" );
142136 fUseParamFilter = p.get < bool >(" UseParamFilter" );
143137 fUseParamFilterInidividualChannel = p.get < bool >(" UseParamFilterInidividualChannel" );
144138 fCorrectBaselineOscillations = p.get < bool >(" CorrectBaselineOscillations" );
@@ -520,15 +514,6 @@ std::vector<TComplex> opdet::OpDeconvolutionAlgWiener::DeconvolutionKernel(size_
520514 kernel[k]= fFilterTF1 ->Eval (k*freq_step) / serfft[k] ;
521515 }
522516 }
523- else if (fUseLowPassFilter )
524- {
525- double freq_step=fSamplingFreq /size;
526- for (size_t k=0 ; k<size/2 ; k++){
527- double freq = k*freq_step;
528- if (freq<=fFrequencyCutoff ) kernel[k]= fFilterTF1 ->Eval (k*freq_step) / serfft[k];
529- else kernel[k]=0 ;
530- }
531- }
532517 else {
533518 // Build Wiener filter kernel: G = Conj(R) / ( |R|^2 + |N|^2/|L|^2)
534519 // R=Detector resopnse FFT
0 commit comments