icarusutil::SignalShapingICARUSService Class Reference

#include <SignalShapingICARUSService_service.h>

Public Member Functions
	SignalShapingICARUSService (const fhicl::ParameterSet &pset, art::ActivityRegistry &reg)
void	reconfigure (const fhicl::ParameterSet &pset)
DoubleVec2	GetNoiseFactVec ()
double	GetASICGain (unsigned int const channel) const
double	GetShapingTime (unsigned int const planeIdx) const
double	GetRawNoise (unsigned int const channel) const
double	GetDeconNoise (unsigned int const channel) const
const icarus_tool::IResponse &	GetResponse (size_t channel) const
int	ResponseTOffset (unsigned int const channel) const
void	SetDecon (double samplingRate, size_t fftsize, size_t channel)
double	GetDeconNorm ()

Private Types
using	IResponsePtr = std::unique_ptr< icarus_tool::IResponse >
using	ResponseVec = std::vector< IResponsePtr >
using	PlaneToResponseMap = std::map< size_t, ResponseVec >

Private Member Functions
void	init () const
void	init ()

Private Attributes
bool	fInit
	Initialization flag. More...
size_t	fPlaneForNormalization
	Normalize responses to this plane. More...
double	fDeconNorm
	Final normalization to apply. More...
DoubleVec2	fNoiseFactVec
	RMS noise in ADCs for lowest gain setting. More...
size_t	fInitialFFTSize
	Size we initially initalize the responses. More...
bool	fPrintResponses
bool	fStoreHistograms
PlaneToResponseMap	fPlaneToResponseMap

Detailed Description

Definition at line 96 of file SignalShapingICARUSService_service.h.

Member Typedef Documentation

using icarusutil::SignalShapingICARUSService::IResponsePtr = std::unique_ptr<icarus_tool::IResponse>

private

Definition at line 127 of file SignalShapingICARUSService_service.h.

using icarusutil::SignalShapingICARUSService::PlaneToResponseMap = std::map<size_t, ResponseVec>

private

Definition at line 129 of file SignalShapingICARUSService_service.h.

using icarusutil::SignalShapingICARUSService::ResponseVec = std::vector<IResponsePtr>

private

Definition at line 128 of file SignalShapingICARUSService_service.h.

Constructor & Destructor Documentation

icarusutil::SignalShapingICARUSService::SignalShapingICARUSService	(	const fhicl::ParameterSet &	pset,
		art::ActivityRegistry &	reg
	)

Definition at line 43 of file SignalShapingICARUSService_service.cc.

: fInit(false)
{
    reconfigure(pset);
}

Member Function Documentation

double icarusutil::SignalShapingICARUSService::GetASICGain

(

unsigned int const

channel

)

const

Definition at line 177 of file SignalShapingICARUSService_service.cc.

{
    static const double fcToElectrons(6241.50975);

    art::ServiceHandle<geo::Geometry> geom;
    size_t planeIdx = geom->ChannelToWire(channel)[0].Plane;
    
    double gain = fPlaneToResponseMap.at(planeIdx).front()->getElectronicsResponse()->getFCperADCMicroS() * fcToElectrons;
    
    return gain;
}

double icarusutil::SignalShapingICARUSService::GetDeconNoise

(

unsigned int const

channel

)

const

Definition at line 225 of file SignalShapingICARUSService_service.cc.

{
    art::ServiceHandle<geo::Geometry> geom;
    size_t planeIdx = geom->ChannelToWire(channel)[0].Plane;
    
    double shaping_time = fPlaneToResponseMap.at(planeIdx).front()->getElectronicsResponse()->getASICShapingTime();
    int temp;
    
    if (std::abs(shaping_time - 0.6)<1e-6){
        temp = 0;
    }else if (std::abs(shaping_time - 1.3)<1e-6){
        temp = 1;
    }else if (std::abs(shaping_time - 2.0)<1.e-6){
        temp = 2;
    }else{
        temp = 3;
    }
    auto   tempNoise  = fNoiseFactVec.at(planeIdx);
    double deconNoise = tempNoise.at(temp);
    
    //deconNoise = deconNoise /4096.*2000./4.7 *6.241*1000/fDeconNorm; <== I don't know where these numbers come from...

    return deconNoise;
}

double icarusutil::SignalShapingICARUSService::GetDeconNorm ( )

inline

Definition at line 121 of file SignalShapingICARUSService_service.h.

{return fDeconNorm;};

DoubleVec2 icarusutil::SignalShapingICARUSService::GetNoiseFactVec ( )

inline

Definition at line 108 of file SignalShapingICARUSService_service.h.

{ return fNoiseFactVec; }

double icarusutil::SignalShapingICARUSService::GetRawNoise

(

unsigned int const

channel

)

const

Definition at line 197 of file SignalShapingICARUSService_service.cc.

{
    art::ServiceHandle<geo::Geometry> geom;
    size_t planeIdx = geom->ChannelToWire(channel)[0].Plane;

    double gain         = fPlaneToResponseMap.at(planeIdx).front()->getElectronicsResponse()->getFCperADCMicroS();
    double shaping_time = fPlaneToResponseMap.at(planeIdx).front()->getElectronicsResponse()->getASICShapingTime();
    int    temp;
    
    if (std::abs(shaping_time - 0.6)<1e-6){
        temp = 0;
    }else if (std::abs(shaping_time - 1.3)<1e-6){
        temp = 1;
    }else if (std::abs(shaping_time - 2.0)<1e-6){
        temp = 2;
    }else{
        temp = 3;
    }
    
    double rawNoise;
    
    auto tempNoise = fNoiseFactVec.at(planeIdx);
    rawNoise = tempNoise.at(temp);
    
    rawNoise *= gain/4.7;
    return rawNoise;
}

const icarus_tool::IResponse & icarusutil::SignalShapingICARUSService::GetResponse

(

size_t

channel

)

const

Definition at line 98 of file SignalShapingICARUSService_service.cc.

{
    if (!fInit) init();
    
    art::ServiceHandle<geo::Geometry> geom;
    
        //use channel number to set some useful numbers
    size_t planeIdx = geom->ChannelToWire(channel)[0].Plane;

    return *fPlaneToResponseMap.at(planeIdx).front();
}

double icarusutil::SignalShapingICARUSService::GetShapingTime

(

unsigned int const

planeIdx

)

const

Definition at line 190 of file SignalShapingICARUSService_service.cc.

{
    double shaping_time = fPlaneToResponseMap.at(planeIdx).front()->getElectronicsResponse()->getASICShapingTime();

    return shaping_time;
}

void icarusutil::SignalShapingICARUSService::init ( ) const

inlineprivate

Definition at line 132 of file SignalShapingICARUSService_service.h.

{const_cast<SignalShapingICARUSService*>(this)->init();}

void icarusutil::SignalShapingICARUSService::init ( )

private

Definition at line 115 of file SignalShapingICARUSService_service.cc.

{
    // First get a lock to make sure we are clear to run
    tbb::spin_mutex::scoped_lock lock(signalShapingSpinMutex);

    if(!fInit)
    {
        // Do ICARUS-specific configuration of SignalShaping by providing
        // ICARUS response and filter functions.
        art::ServiceHandle<geo::Geometry> geo;

        auto const samplingRate = sampling_rate(art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob());

        // Get the normalization from the field response for the collection plane
        double integral = fPlaneToResponseMap.at(fPlaneForNormalization).front().get()->getFieldResponse()->getIntegral();
        double weight   = 1. / integral;
        
        for(size_t planeIdx = 0; planeIdx < geo->Nplanes(); planeIdx++)
        {
          fPlaneToResponseMap[planeIdx].front().get()->setResponse(samplingRate, weight);
        }
        
        // Check to see if we want histogram output
        if (fStoreHistograms)
        {
            art::ServiceHandle<art::TFileService> tfs;
            
            // Make sure we are at the top level
            tfs->file().cd();
            
            // Make a directory for these histograms
            art::TFileDirectory dir = tfs->mkdir("SignalShaping");
            
            // Loop through response tools first
            for(const auto& response: fPlaneToResponseMap) response.second.front().get()->outputHistograms(samplingRate, dir);

        }

        fInit = true;
    }
    
    return;
}

void icarusutil::SignalShapingICARUSService::reconfigure

(

const fhicl::ParameterSet &

pset

)

Definition at line 52 of file SignalShapingICARUSService_service.cc.

{
    mf::LogInfo("fSignalShapingICARUSService") << " reconfiguring setup " ;
    
    // If called again, then we need to clear out the existing tools...
    fPlaneToResponseMap.clear();
    
    // Implement the tools for handling the responses
    const fhicl::ParameterSet& responseTools = pset.get<fhicl::ParameterSet>("ResponseTools");
    
    for(const std::string& responseTool : responseTools.get_pset_names())
    {
        const fhicl::ParameterSet& responseToolParamSet = responseTools.get<fhicl::ParameterSet>(responseTool);
        size_t                     planeIdx             = responseToolParamSet.get<size_t>("Plane");
        
        fPlaneToResponseMap[planeIdx].push_back(art::make_tool<icarus_tool::IResponse>(responseToolParamSet));
    }
    
    fPlaneForNormalization  = pset.get<size_t>(    "PlaneForNormalization");
    fPrintResponses         = pset.get<bool>(      "PrintResponses"       );
    fDeconNorm              = pset.get<double>(    "DeconNorm"            );
    fInitialFFTSize         = pset.get<size_t>(    "InitialFFTSize"       );
    fNoiseFactVec           = pset.get<DoubleVec2>("NoiseFactVec"         );
    fStoreHistograms        = pset.get<bool>(      "StoreHistograms"      );
    
    fInit = false;
    
    return;
}

int icarusutil::SignalShapingICARUSService::ResponseTOffset

(

unsigned int const

channel

)

const

Definition at line 250 of file SignalShapingICARUSService_service.cc.

{
    if (!fInit) init();
    
    art::ServiceHandle<geo::Geometry> geom;
    
    size_t planeIdx = geom->ChannelToWire(channel)[0].Plane;

    return fPlaneToResponseMap.at(planeIdx).front()->getTOffset();
}

void icarusutil::SignalShapingICARUSService::SetDecon	(	double	samplingRate,
		size_t	fftsize,
		size_t	channel
	)

Definition at line 159 of file SignalShapingICARUSService_service.cc.

{
    art::ServiceHandle<geo::Geometry> geo;
    
    if (!fInit) init();
    
    // Assume we need to reset the kernels
    double integral = fPlaneToResponseMap.at(fPlaneForNormalization).front().get()->getFieldResponse()->getIntegral();
    double weight   = 1. / integral;
    
    for(size_t planeIdx = 0; planeIdx < geo->Nplanes(); planeIdx++)
    {
        fPlaneToResponseMap.at(planeIdx).front()->setResponse(samplingRate, weight);
    }
}

Member Data Documentation

double icarusutil::SignalShapingICARUSService::fDeconNorm

private

Final normalization to apply.

Definition at line 140 of file SignalShapingICARUSService_service.h.

bool icarusutil::SignalShapingICARUSService::fInit

private

Initialization flag.

Definition at line 136 of file SignalShapingICARUSService_service.h.

size_t icarusutil::SignalShapingICARUSService::fInitialFFTSize

private

Size we initially initalize the responses.

Definition at line 142 of file SignalShapingICARUSService_service.h.

DoubleVec2 icarusutil::SignalShapingICARUSService::fNoiseFactVec

private

RMS noise in ADCs for lowest gain setting.

Definition at line 141 of file SignalShapingICARUSService_service.h.

size_t icarusutil::SignalShapingICARUSService::fPlaneForNormalization

private

Normalize responses to this plane.

Definition at line 139 of file SignalShapingICARUSService_service.h.

PlaneToResponseMap icarusutil::SignalShapingICARUSService::fPlaneToResponseMap

private

Definition at line 147 of file SignalShapingICARUSService_service.h.

bool icarusutil::SignalShapingICARUSService::fPrintResponses

private

Definition at line 143 of file SignalShapingICARUSService_service.h.

bool icarusutil::SignalShapingICARUSService::fStoreHistograms

private

Definition at line 144 of file SignalShapingICARUSService_service.h.

---

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

• SignalShapingICARUSService_service.h
• SignalShapingICARUSService_service.cc