#include <AFP_NoisyPixelTool.h>

Inheritance diagram for AFP_NoisyPixelTool:

Collaboration diagram for AFP_NoisyPixelTool:

Public Member Functions
	AFP_NoisyPixelTool ()
int	Identify (std::shared_ptr< const TH2F > input, std::vector< TH2F > &output) const override
std::vector< std::vector< double > >	makeFits (std::shared_ptr< const TH2F > input) const
void	setThreshNoisy (double t)
double	getThreshNoisy () const
void	setThreshLeff (double t)
double	getThreshLeff () const
void	setSensitivity (double s)
double	getSensitivity () const
void	setMethods (const std::vector< std::string > &m)
const std::vector< std::string > &	getMethods () const

Private Member Functions
std::vector< std::pair< int, int > >	getLegitPixels (std::shared_ptr< const TH2F > input, const int col_ID, const int row_ID, const std::string &method) const
TH2F	countInactivePixelsAround (std::shared_ptr< const TH2F > input) const
double	getNeighbours (std::shared_ptr< const TH2F > input, int row_ID, int col_ID) const
std::tuple< TH2F, TH2F, TH2F, TH2F >	findLEffAndNoisyPixels (std::shared_ptr< const TH2F > input, const TH2F &inact, const std::string &method) const
void	filterLEffPixelsAroundNoisy (std::shared_ptr< const TH2F > input, const TH2F &inact, TH2F &fleff, TH2F &fnoisy, TH2F &eleff, TH2F &enoisy, const std::string &method) const

Private Attributes
double	m_threshNoisy
double	m_threshLEff
double	m_sensitivity
std::vector< std::string >	m_methods

Detailed Description

Definition at line 22 of file AFP_NoisyPixelTool.h.

Constructor & Destructor Documentation

◆ AFP_NoisyPixelTool()

AFP_NoisyPixelTool::AFP_NoisyPixelTool ( )

inline

Definition at line 24 of file AFP_NoisyPixelTool.h.

24: m_threshNoisy(1.4), m_threshLEff(0.6), m_sensitivity(0.0), m_methods({"2_ROW","2_COL"}) {}

AFP_NoisyPixelTool::m_sensitivity

double m_sensitivity

Definition AFP_NoisyPixelTool.h:46

AFP_NoisyPixelTool::m_threshLEff

double m_threshLEff

Definition AFP_NoisyPixelTool.h:45

AFP_NoisyPixelTool::m_methods

std::vector< std::string > m_methods

Definition AFP_NoisyPixelTool.h:48

AFP_NoisyPixelTool::m_threshNoisy

double m_threshNoisy

Definition AFP_NoisyPixelTool.h:44

Member Function Documentation

◆ countInactivePixelsAround()

TH2F AFP_NoisyPixelTool::countInactivePixelsAround ( std::shared_ptr< const TH2F > input ) const

private

Definition at line 183 of file AFP_NoisyPixelTool.cxx.

{
    TH2F tmp_inact_pix_around(*input);
    tmp_inact_pix_around.SetName("tmp_inact_pix_around");
    tmp_inact_pix_around.Reset();
            
    for(int col_ID=1; col_ID<=input->GetNbinsY(); col_ID++)
    {
        for(int row_ID=1; row_ID<=input->GetNbinsX(); row_ID++)
        {                   
            int inactive_pixels_around=0;
            auto legit_pixels=getLegitPixels(input, col_ID, row_ID);
            
            for(auto legpix : legit_pixels)
            {
                if(input->GetBinContent(legpix.first, legpix.second)<1) ++inactive_pixels_around;
            }
 
            tmp_inact_pix_around.SetBinContent(row_ID,col_ID, inactive_pixels_around+0.01);
        }
    }
    
    return tmp_inact_pix_around;
}

◆ filterLEffPixelsAroundNoisy()

void AFP_NoisyPixelTool::filterLEffPixelsAroundNoisy	(	std::shared_ptr< const TH2F >	input,
		const TH2F &	inact,
		TH2F &	fleff,
		TH2F &	fnoisy,
		TH2F &	eleff,
		TH2F &	enoisy,
		const std::string &	method ) const

private

Definition at line 289 of file AFP_NoisyPixelTool.cxx.

{   
    for(int col_ID=1; col_ID<=input->GetNbinsY(); col_ID++)
    {
        for(int row_ID=1; row_ID<=input->GetNbinsX(); row_ID++)
        {   
            if(input->GetBinContent(row_ID,col_ID)==0) continue;
                    
            int inactive_pixels_around=tmp_inact_pix_around.GetBinContent(row_ID, col_ID);
            int leff_pixels_around=0;
            auto legit_pixels=getLegitPixels(input, col_ID,row_ID);
            for(auto legpix : legit_pixels)
            {
                if(tmp_found_leff.GetBinContent(legpix.first,legpix.second)>0) leff_pixels_around++;
            }
            
            if(leff_pixels_around!=0  && inactive_pixels_around==0)
            {
                //if low eff. pixels around > 0, recalculate avarage
                double re_av=0.;
                int npix=0;
                
                auto legit_pixels=getLegitPixels(input, col_ID,row_ID, method);
                for(auto legpix : legit_pixels)
                {
                    if(tmp_found_leff.GetBinContent(legpix.first,legpix.second)==0) 
                    {
                        re_av+=input->GetBinContent(legpix.first,legpix.second);
                        npix++;
                    }
                }
            
                if(npix==0)
                {
                    re_av=0;
                    for(const auto & legpix : legit_pixels)
                    {
                        re_av+=input->GetBinContent(legpix.first,legpix.second);
                    }
                    
                    if(re_av/input->GetBinContent(row_ID,col_ID)<1.5 && re_av/input->GetBinContent(row_ID,col_ID)>0.5) re_av=1;
                }
                else
                {
                    re_av/=npix;
                }
                if (re_av == 0)[[unlikely]]{
                  std::cerr<< "AFP_NoisyPixelTool::findLEffAndNoisyPixels: re_av denominator is zero\n";
                  continue;
                } 
                double ratio_noisy = m_sensitivity+input->GetBinContent(row_ID,col_ID)/re_av;
                double ratio_leff  = input->GetBinContent(row_ID,col_ID)/re_av; 
                if(re_av==1)
                {
                    tmp_eff_leff.SetBinContent(row_ID,col_ID,1);
                    tmp_eff_noisy.SetBinContent(row_ID,col_ID,1);
                }
                        
                if(row_ID!=input->GetNbinsX() && row_ID!=1 && col_ID!=1 && col_ID!=input->GetNbinsY())
                {
                    if( ((ratio_noisy>m_threshNoisy && re_av!=1) || re_av==1) )
                    {
                        tmp_found_noisy.SetBinContent(row_ID,col_ID,1);
                        if(re_av!=1) tmp_eff_noisy.SetBinContent(row_ID,col_ID,100.*ratio_noisy);       
                    }
                    if( ((ratio_leff<m_threshLEff && re_av!=1) || re_av==1) )
                    {
                        tmp_found_leff.SetBinContent(row_ID,col_ID,1);      
                        if(re_av!=1) tmp_eff_leff.SetBinContent(row_ID,col_ID,100.*ratio_leff);
                    }
                }
            }
        }
    }
    
    return;
}

◆ findLEffAndNoisyPixels()

std::tuple< TH2F, TH2F, TH2F, TH2F > AFP_NoisyPixelTool::findLEffAndNoisyPixels	(	std::shared_ptr< const TH2F >	input,
		const TH2F &	inact,
		const std::string &	method ) const

private

Definition at line 217 of file AFP_NoisyPixelTool.cxx.

{
    TH2F tmp_found_leff(*input);
    tmp_found_leff.SetName("tmp_found_leff");
    tmp_found_leff.Reset();
    
    TH2F tmp_found_noisy(*input);
    tmp_found_noisy.SetName("tmp_found_noisy");
    tmp_found_noisy.Reset();
    
    TH2F tmp_eff_leff(*input);
    tmp_eff_leff.SetName("tmp_eff_leff");
    tmp_eff_leff.Reset();
    
    TH2F tmp_eff_noisy(*input);
    tmp_eff_noisy.SetName("tmp_eff_noisy");
    tmp_eff_noisy.Reset();
    
    for(int col_ID=1; col_ID<=input->GetNbinsY(); col_ID++)
    {
        for(int row_ID=1; row_ID<=input->GetNbinsX(); row_ID++)
        {
            if(input->GetBinContent(row_ID, col_ID)!=0 && tmp_inact_pix_around.GetBinContent(row_ID, col_ID)<0.1)
            {   
                double npixels=getNeighbours(input, row_ID, col_ID);
                
                double sum=0.;
                auto legit_pixels=getLegitPixels(input, col_ID, row_ID, method);
                for(auto legpix : legit_pixels)
                {
                    sum+=input->GetBinContent(legpix.first,legpix.second);
                }
                
                double av_noisy=0;
                double av_leff=0;
                
                // TODO: something about edges?
                if(method=="2_ROW" || method=="2_COL")
                {
                    av_noisy=sum/(std::abs(npixels-4)*0.5);
                    av_leff=sum/2.;
                }
                else if(method=="4_PIX") 
                {
                    av_noisy=sum/(std::ceil(npixels*0.5));
                    av_leff=sum/4.;
                }
                else if(method=="8_PIX")
                {
                    av_noisy=sum/(npixels);
                    av_leff=sum/8.;
                }
                
                if ((av_leff == 0) or (av_noisy == 0))[[unlikely]]{
                  std::cerr<< "AFP_NoisyPixelTool::findLEffAndNoisyPixels: denominator is zero\n";
                  continue;
                } 
                double ratio_leff  = input->GetBinContent(row_ID,col_ID)/av_leff;
                double ratio_noisy = input->GetBinContent(row_ID,col_ID)/av_noisy;
                tmp_eff_leff.SetBinContent(row_ID,col_ID,100.*ratio_leff);
                tmp_eff_noisy.SetBinContent(row_ID,col_ID,100.*ratio_noisy);
                
                if(ratio_leff  < m_threshLEff) tmp_found_leff.SetBinContent(row_ID,col_ID,1.); 
                if(ratio_noisy > m_threshNoisy) tmp_found_noisy.SetBinContent(row_ID,col_ID,1.);   
            }
        }
    }
    
    return {tmp_found_leff,tmp_found_noisy,tmp_eff_leff,tmp_eff_noisy};
}

◆ getLegitPixels()

std::vector< std::pair< int, int > > AFP_NoisyPixelTool::getLegitPixels	(	std::shared_ptr< const TH2F >	input,
		const int	col_ID,
		const int	row_ID,
		const std::string &	method = "8_PIX" ) const

private

Definition at line 158 of file AFP_NoisyPixelTool.cxx.

{
    // method="8_PIX" means to investigate all pixels around
    
    std::vector<std::pair<int,int>> legit_pixels={};
    
    for(int r=-1; r<2; r++)
    {
        for(int c=-1; c<2; c++)
        {
            if( (row_ID+r)>=1 && (col_ID+c)>=1 && (col_ID+c)<=input->GetNbinsY() && (row_ID+r)<=input->GetNbinsX() && (r!=0 || c!=0))
            {
                if( (method=="2_ROW" && c==0) || (method=="2_COL" && r==0)
                    || (method=="4_PIX" && (c==0 || r==0)) || (method=="8_PIX") )
                {
                    legit_pixels.push_back(std::pair<int,int>(row_ID+r,col_ID+c));
                }
            }
        }
    }
    
    return legit_pixels;
}

◆ getMethods()

const std::vector< std::string > & AFP_NoisyPixelTool::getMethods ( ) const

inline

Definition at line 40 of file AFP_NoisyPixelTool.h.

40{return m_methods;}

◆ getNeighbours()

double AFP_NoisyPixelTool::getNeighbours	(	std::shared_ptr< const TH2F >	input,
		int	row_ID,
		int	col_ID ) const

private

Definition at line 209 of file AFP_NoisyPixelTool.cxx.

{
    if( row_ID!=input->GetNbinsX() && row_ID!=1 && col_ID!=1 && col_ID!=input->GetNbinsY()) return 8.;
    else if( (row_ID==input->GetNbinsX() || row_ID==1) && (col_ID==1 || col_ID==input->GetNbinsY())) return 3.;
    else return 5;
}

◆ getSensitivity()

double AFP_NoisyPixelTool::getSensitivity ( ) const

inline

Definition at line 37 of file AFP_NoisyPixelTool.h.

37{return m_sensitivity;}

◆ getThreshLeff()

double AFP_NoisyPixelTool::getThreshLeff ( ) const

inline

Definition at line 34 of file AFP_NoisyPixelTool.h.

34{return m_threshLEff;}

◆ getThreshNoisy()

double AFP_NoisyPixelTool::getThreshNoisy ( ) const

inline

Definition at line 31 of file AFP_NoisyPixelTool.h.

31{return m_threshNoisy;}

◆ Identify()

int AFP_NoisyPixelTool::Identify	(	std::shared_ptr< const TH2F >	input,
		std::vector< TH2F > &	output ) const

overridevirtual

Implements IAFP_GenericPixelTool.

Definition at line 8 of file AFP_NoisyPixelTool.cxx.

{
    if(output.size()!=1)
    {
        return 0;
    }
    
    output[0].Reset();
    output.resize(4, output[0]);//container now has 4 copies of the 0'th element
    auto name = [&output](const char * n)->TString{
      return Form(n,output[0].GetName());
    };
    auto title = [&output](const char * t)->TString{
      return Form(t,output[0].GetTitle());
    };
    output[1].SetNameTitle(name("leffpixels_found_%s"), title("low efficiency pixels, found, %s"));
    output[2].SetNameTitle(name("noisypixels_eff_%s"), title("noisy pixels, efficiency, %s"));
    output[3].SetNameTitle(name("leffpixels_eff_%s"), title("low efficiency pixels, efficiency, %s"));
    output[0].SetNameTitle(name("noisypixels_found_%s"), title("noisy pixels, found, %s"));
    
  TH2F& noisypixels_found_output = output[0];
    TH2F& leffpixels_found_output = output[1];
    TH2F& noisypixels_eff_output = output[2];
    TH2F& leffpixels_eff_output = output[3];
        
    
    if(input->GetMaximum()<0.5) return 0;
 
    std::vector<TH2F> vec_found_leff{},vec_found_noisy{},vec_eff_leff{},vec_eff_noisy{};
 
    for(const auto& method : m_methods)
    {
        TH2F tmp_eff_noisy, tmp_eff_leff;
        TH2F tmp_found_leff, tmp_found_noisy;
        
        if(method!="FIT")
        {
            // count inactive pixels around each pixel
            TH2F tmp_inact_pix_around=countInactivePixelsAround(input);
            
            // find low eff. and noisy pixels
            std::tie(tmp_found_leff,tmp_found_noisy,tmp_eff_leff, tmp_eff_noisy) = findLEffAndNoisyPixels(input, tmp_inact_pix_around, method);
            
            // reconsider some pixels
            filterLEffPixelsAroundNoisy(input, tmp_inact_pix_around, tmp_found_leff,tmp_found_noisy,tmp_eff_leff,tmp_eff_noisy, method);
        }
        else
        {
            // method is "FIT"
            
            tmp_eff_noisy=TH2F(*input);
            tmp_eff_noisy.SetName("tmp_eff_noisy");
            tmp_eff_noisy.Reset();
            tmp_eff_leff=TH2F(*input);
            tmp_eff_leff.SetName("tmp_eff_leff");
            tmp_eff_leff.Reset();
            tmp_found_noisy=TH2F(*input);
            tmp_found_noisy.SetName("tmp_found_noisy");
            tmp_found_noisy.Reset();
            tmp_found_leff=TH2F(*input);
            tmp_found_leff.SetName("tmp_found_leff");
            tmp_found_leff.Reset();
            
            auto p=makeFits(input);
            
            const double Threshold_LEFF=5.;
            const double Threshold_NOISY=(m_threshNoisy-1)*100;
            
            for(int col_ID=1; col_ID<=input->GetNbinsY(); col_ID++)
            {
                for(int row_ID=1; row_ID<=input->GetNbinsX(); row_ID++)
                {
                    double ratio_leff=0, ratio_noisy=0;
                    double gauss_helper_12 = (col_ID-p[1][row_ID])/p[2][row_ID];
                    double gauss_helper_45 = (col_ID-p[4][row_ID])/p[5][row_ID];
                    double fit = p[0][row_ID]*exp(-0.5*gauss_helper_12*gauss_helper_12)
                               + p[3][row_ID]*exp(-0.5*gauss_helper_45*gauss_helper_45) + p[6][row_ID];
                    double sigma=std::sqrt(fit);
                    if(sigma<1) sigma=1;
                    
                    double bin_content=input->GetBinContent(row_ID,col_ID);
                    
                    if(fit>=bin_content) ratio_leff = std::abs(fit-bin_content)/sigma;
                    else                 ratio_noisy = std::abs(fit-bin_content)/sigma;
                    
                    if(bin_content!=0)
                    {
                        if( ratio_leff  > Threshold_LEFF) tmp_found_leff.SetBinContent(row_ID,col_ID,1);
                        if( ratio_noisy > Threshold_NOISY) tmp_found_noisy.SetBinContent(row_ID,col_ID,1);
                    }
                    tmp_eff_leff.SetBinContent(row_ID,col_ID,100.*ratio_leff);
                    tmp_eff_noisy.SetBinContent(row_ID,col_ID,100.*ratio_noisy);
                }
            }
        }
        
        // save output for this method
        vec_found_leff.push_back(tmp_found_leff);
        vec_found_noisy.push_back(tmp_found_noisy);
        vec_eff_leff.push_back(tmp_eff_leff);
        vec_eff_noisy.push_back(tmp_eff_noisy); 
    }
    
    // sum found pixels for all methods, set efficiency to the maximum of all methods
    for(unsigned int m=0; m<m_methods.size(); m++)
    {   
        for(int col_ID=1; col_ID<=input->GetNbinsY(); col_ID++)
        {
            for(int row_ID=1; row_ID<=input->GetNbinsX(); row_ID++)
            {   
                noisypixels_found_output.Fill(row_ID,col_ID, vec_found_noisy[m].GetBinContent(row_ID,col_ID));
                leffpixels_found_output.Fill(row_ID,col_ID, vec_found_leff[m].GetBinContent(row_ID,col_ID));
                
                if(m==0)
                {
                    noisypixels_eff_output.SetBinContent(row_ID,col_ID, vec_eff_noisy[m].GetBinContent(row_ID,col_ID));
                    leffpixels_eff_output.SetBinContent(row_ID,col_ID, vec_eff_leff[m].GetBinContent(row_ID,col_ID));
                }
                else
                { 
                    if(noisypixels_eff_output.GetBinContent(row_ID,col_ID)<vec_eff_noisy[m].GetBinContent(row_ID,col_ID))
                    {
                        noisypixels_eff_output.SetBinContent(row_ID,col_ID, vec_eff_noisy[m].GetBinContent(row_ID,col_ID));
                    }
                    if(leffpixels_eff_output.GetBinContent(row_ID,col_ID)>vec_eff_leff[m].GetBinContent(row_ID,col_ID))
                    {
                        leffpixels_eff_output.SetBinContent(row_ID,col_ID, vec_eff_leff[m].GetBinContent(row_ID,col_ID));
                    }
                }
            }
        }
    }
 
    // report only pixels that are identified by all methods
    for(int col_ID=1; col_ID<=input->GetNbinsY(); col_ID++)
    {
        for(int row_ID=1; row_ID<=input->GetNbinsX(); row_ID++)
        {   
            if(noisypixels_found_output.GetBinContent(row_ID,col_ID)==m_methods.size())    noisypixels_found_output.SetBinContent(row_ID,col_ID,1);
            else noisypixels_found_output.SetBinContent(row_ID,col_ID,0);
            
            if(leffpixels_found_output.GetBinContent(row_ID,col_ID)==m_methods.size()) leffpixels_found_output.SetBinContent(row_ID,col_ID,1);
            else leffpixels_found_output.SetBinContent(row_ID,col_ID,0);
        }
    }
 
    return 0;
}

◆ makeFits()

std::vector< std::vector< double > > AFP_NoisyPixelTool::makeFits ( std::shared_ptr< const TH2F > input ) const

Definition at line 368 of file AFP_NoisyPixelTool.cxx.

{
    const int nParams=7;
    std::vector<std::vector<double>> params(nParams, std::vector<double>(input->GetNbinsX()+1, 0.));
    
    for(int row_ID=0; row_ID<=input->GetNbinsX(); row_ID++)
    {
        std::string srow_ID = std::to_string(row_ID);
        std::unique_ptr<TH1F> hist(new TH1F(("ROW"+srow_ID).c_str(), ("ROW"+srow_ID).c_str(), input->GetNbinsY(),  0.5,  input->GetNbinsX()+0.5));
        
        for(int col_ID=1; col_ID<=input->GetNbinsY(); col_ID++)
        {  
            hist->Fill(col_ID, 1.0*input->GetBinContent(row_ID,col_ID)); 
        }
        
        //Set double gaussian fit with starting parameters
        std::shared_ptr<TF1> FIT_2(new TF1("gaus","gaus(0)+gaus(3)+[6]",0,input->GetNbinsY(),nParams));
        if(row_ID==1)
        {
            FIT_2->SetParameters(25000,hist->GetMean(),hist->GetStdDev()/std::sqrt(2),
                                  5000,hist->GetMean(),hist->GetStdDev()/std::sqrt(2),
                                 std::max(hist->GetBinContent(2),std::max(hist->GetBinContent(3),hist->GetBinContent(4))) );
        }
        else
        {
            for(int par=0;par<nParams;++par)
            {
                FIT_2->SetParameter(par,params.at(par).at(row_ID-1));
            }
        }
        
        //fit
        hist->Fit(FIT_2.get(),"Q");
        
        for(int par=0;par<nParams;++par)
        {
            params.at(par).at(row_ID-1)=FIT_2->GetParameter(par);
        }
    }
    
    return params;
}

◆ setMethods()

void AFP_NoisyPixelTool::setMethods ( const std::vector< std::string > & m )

inline

Definition at line 39 of file AFP_NoisyPixelTool.h.

39{m_methods=m;}

◆ setSensitivity()

void AFP_NoisyPixelTool::setSensitivity ( double s )

inline

Definition at line 36 of file AFP_NoisyPixelTool.h.

36{m_sensitivity=s;}

python.SystemOfUnits.s

float s

Definition SystemOfUnits.py:147

◆ setThreshLeff()

void AFP_NoisyPixelTool::setThreshLeff ( double t )

inline

Definition at line 33 of file AFP_NoisyPixelTool.h.

33{m_threshLEff=t;}

◆ setThreshNoisy()

void AFP_NoisyPixelTool::setThreshNoisy ( double t )

inline

Definition at line 30 of file AFP_NoisyPixelTool.h.

30{m_threshNoisy=t;}

Member Data Documentation

◆ m_methods

std::vector<std::string> AFP_NoisyPixelTool::m_methods

private

Definition at line 48 of file AFP_NoisyPixelTool.h.

◆ m_sensitivity

double AFP_NoisyPixelTool::m_sensitivity

private

Definition at line 46 of file AFP_NoisyPixelTool.h.

◆ m_threshLEff

double AFP_NoisyPixelTool::m_threshLEff

private

Definition at line 45 of file AFP_NoisyPixelTool.h.

◆ m_threshNoisy

double AFP_NoisyPixelTool::m_threshNoisy

private

Definition at line 44 of file AFP_NoisyPixelTool.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ AFP_NoisyPixelTool()

Member Function Documentation

◆ countInactivePixelsAround()

◆ filterLEffPixelsAroundNoisy()

◆ findLEffAndNoisyPixels()

◆ getLegitPixels()

◆ getMethods()

◆ getNeighbours()

◆ getSensitivity()

◆ getThreshLeff()

◆ getThreshNoisy()

◆ Identify()

◆ makeFits()

◆ setMethods()

◆ setSensitivity()

◆ setThreshLeff()

◆ setThreshNoisy()

Member Data Documentation

◆ m_methods

◆ m_sensitivity

◆ m_threshLEff

◆ m_threshNoisy