PixelBarrelTemplateSmearerPlugin Class Reference

Inheritance diagram for PixelBarrelTemplateSmearerPlugin:

Public Member Functions
	PixelBarrelTemplateSmearerPlugin (const std::string &name, const edm::ParameterSet &config, edm::ConsumesCollector &consumesCollector)
	~PixelBarrelTemplateSmearerPlugin () override
Public Member Functions inherited from PixelTemplateSmearerBase
bool	hitsMerge (const PSimHit &simHit1, const PSimHit &simHit2) const
	PixelTemplateSmearerBase (const std::string &name, const edm::ParameterSet &config, edm::ConsumesCollector &consumesCollector)
TrackingRecHitProductPtr	process (TrackingRecHitProductPtr product) const override
TrackingRecHitProductPtr	processMergeGroups (std::vector< MergeGroup * > &mergeGroups, TrackingRecHitProductPtr product, const PixelGeomDetUnit *detUnit, const double boundX, const double boundY, RandomEngineAndDistribution const *random) const
TrackingRecHitProductPtr	processUnmergedHits (std::vector< TrackingRecHitProduct::SimHitIdPair > &unmergedHits, TrackingRecHitProductPtr product, const PixelGeomDetUnit *detUnit, const double boundX, const double boundY, RandomEngineAndDistribution const *random) const
FastSingleTrackerRecHit	smearHit (const PSimHit &simHit, const PixelGeomDetUnit *detUnit, const double boundX, const double boundY, RandomEngineAndDistribution const *) const
FastSingleTrackerRecHit	smearMergeGroup (MergeGroup *mg, const PixelGeomDetUnit *detUnit, const double boundX, const double boundY, const RandomEngineAndDistribution *random) const
	~PixelTemplateSmearerBase () override
Public Member Functions inherited from TrackingRecHitAlgorithm
virtual void	beginEvent (edm::Event &event, const edm::EventSetup &eventSetup)
virtual void	beginStream (const edm::StreamID &id)
virtual void	endEvent (edm::Event &event, const edm::EventSetup &eventSetup)
virtual void	endStream ()
const TrackerGeometry &	getMisalignedGeometry () const
const std::string &	getName () const
const RandomEngineAndDistribution &	getRandomEngine () const
const std::string &	getSelectionString () const
const TrackerGeometry &	getTrackerGeometry () const
const TrackerTopology &	getTrackerTopology () const
	TrackingRecHitAlgorithm (const std::string &name, const edm::ParameterSet &config, edm::ConsumesCollector &consumesCollector)
virtual	~TrackingRecHitAlgorithm ()

Private Member Functions
void	initializeBarrel ()

Additional Inherited Members
Protected Member Functions inherited from PixelTemplateSmearerBase
bool	isFlipped (const PixelGeomDetUnit *theDet) const
Protected Attributes inherited from PixelTemplateSmearerBase
bool	isForward
bool	mergeHitsOn
double	rescotAlpha_binMin
unsigned int	rescotAlpha_binN
double	rescotAlpha_binWidth
double	rescotBeta_binMin
unsigned int	rescotBeta_binN
double	rescotBeta_binWidth
int	resqbin_binMin
unsigned int	resqbin_binN
int	resqbin_binWidth
int	templateId
std::unique_ptr< TFile >	theBigPixelResolutionFile
std::string	theBigPixelResolutionFileName
std::unique_ptr< TFile >	theEdgePixelResolutionFile
std::string	theEdgePixelResolutionFileName
unsigned int	theLayer
std::unique_ptr< TFile >	theMergedPixelResolutionXFile
std::string	theMergedPixelResolutionXFileName
std::unique_ptr< TFile >	theMergedPixelResolutionYFile
std::string	theMergedPixelResolutionYFileName
std::unique_ptr< TFile >	theMergingProbabilityFile
std::string	theMergingProbabilityFileName
std::vector< SiPixelTemplateStore >	thePixelTemp_
std::unique_ptr< TFile >	theRegularPixelResolutionFile
std::string	theRegularPixelResolutionFileName
std::map< unsigned int, const SimpleHistogramGenerator * >	theXHistos
std::map< unsigned int, const SimpleHistogramGenerator * >	theYHistos

Detailed Description

Definition at line 29 of file PixelBarrelTemplateSmearerPlugin.cc.

Constructor & Destructor Documentation

PixelBarrelTemplateSmearerPlugin::PixelBarrelTemplateSmearerPlugin ( const std::string &  name, const edm::ParameterSet &  config, edm::ConsumesCollector &  consumesCollector  )

explicit

Definition at line 45 of file PixelBarrelTemplateSmearerPlugin.cc.

References Exception, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), initializeBarrel(), PixelTemplateSmearerBase::isForward, SiPixelTemplate::pushfile(), PixelTemplateSmearerBase::templateId, PixelTemplateSmearerBase::theBigPixelResolutionFile, PixelTemplateSmearerBase::theBigPixelResolutionFileName, PixelTemplateSmearerBase::theEdgePixelResolutionFile, PixelTemplateSmearerBase::theEdgePixelResolutionFileName, PixelTemplateSmearerBase::theMergedPixelResolutionXFile, PixelTemplateSmearerBase::theMergedPixelResolutionXFileName, PixelTemplateSmearerBase::theMergedPixelResolutionYFile, PixelTemplateSmearerBase::theMergedPixelResolutionYFileName, PixelTemplateSmearerBase::theMergingProbabilityFile, PixelTemplateSmearerBase::theMergingProbabilityFileName, PixelTemplateSmearerBase::thePixelTemp_, PixelTemplateSmearerBase::theRegularPixelResolutionFile, and PixelTemplateSmearerBase::theRegularPixelResolutionFileName.

 :
    PixelTemplateSmearerBase(name,config,consumesCollector)
{
    isForward = false;

    theBigPixelResolutionFileName = config.getParameter<string>( "BigPixelBarrelResolutionFile" );
    theBigPixelResolutionFile = std::make_unique<TFile>( edm::FileInPath( theBigPixelResolutionFileName ).fullPath().c_str()  ,"READ");
    theEdgePixelResolutionFileName =  config.getParameter<string>( "EdgePixelBarrelResolutionFile" );
    theEdgePixelResolutionFile = std::make_unique<TFile>( edm::FileInPath( theEdgePixelResolutionFileName ).fullPath().c_str()  ,"READ");
    theRegularPixelResolutionFileName = config.getParameter<string>( "RegularPixelBarrelResolutionFile" );
    theRegularPixelResolutionFile = std::make_unique<TFile>( edm::FileInPath( theRegularPixelResolutionFileName ).fullPath().c_str()  ,"READ");
    
    theMergingProbabilityFileName = config.getParameter<string>( "MergingProbabilityBarrelFile" );
    theMergingProbabilityFile = std::make_unique<TFile>( edm::FileInPath( theMergingProbabilityFileName ).fullPath().c_str()  ,"READ");
    theMergedPixelResolutionXFileName = config.getParameter<string>( "MergedPixelBarrelResolutionXFile" );
    theMergedPixelResolutionXFile = std::make_unique<TFile>( edm::FileInPath( theMergedPixelResolutionXFileName ).fullPath().c_str()  ,"READ");
    theMergedPixelResolutionYFileName = config.getParameter<string>( "MergedPixelBarrelResolutionYFile" );
    theMergedPixelResolutionYFile = std::make_unique<TFile>( edm::FileInPath( theMergedPixelResolutionYFileName ).fullPath().c_str()  ,"READ");

    initializeBarrel();
    
    if (!SiPixelTemplate::pushfile(templateId, thePixelTemp_))
    {
        throw cms::Exception("PixelBarrelTemplateSmearerPlugin:")<<"SiPixel Barrel Template Not Loaded Correctly!"<<endl;
    }

}

PixelBarrelTemplateSmearerPlugin::~PixelBarrelTemplateSmearerPlugin ( )

override

Definition at line 78 of file PixelBarrelTemplateSmearerPlugin.cc.

{
}

Member Function Documentation

void PixelBarrelTemplateSmearerPlugin::initializeBarrel ( )

private

Definition at line 83 of file PixelBarrelTemplateSmearerPlugin.cc.

References DEFINE_EDM_PLUGIN, PixelTemplateSmearerBase::rescotAlpha_binMin, PixelTemplateSmearerBase::rescotAlpha_binN, PixelTemplateSmearerBase::rescotAlpha_binWidth, PixelTemplateSmearerBase::rescotBeta_binMin, PixelTemplateSmearerBase::rescotBeta_binN, PixelTemplateSmearerBase::rescotBeta_binWidth, PixelTemplateSmearerBase::resqbin_binMin, PixelTemplateSmearerBase::resqbin_binN, PixelTemplateSmearerBase::resqbin_binWidth, PixelTemplateSmearerBase::theBigPixelResolutionFile, PixelTemplateSmearerBase::theEdgePixelResolutionFile, PixelTemplateSmearerBase::theRegularPixelResolutionFile, PixelTemplateSmearerBase::theXHistos, and PixelTemplateSmearerBase::theYHistos.

Referenced by PixelBarrelTemplateSmearerPlugin().

{
    rescotAlpha_binMin = -0.2;
    rescotAlpha_binWidth = 0.08 ;
    rescotAlpha_binN = 5;
    rescotBeta_binMin = -5.5;
    rescotBeta_binWidth = 1.0;
    rescotBeta_binN = 11;
    resqbin_binMin = 0;
    resqbin_binWidth = 1;
    resqbin_binN = 4;

    for (unsigned cotalphaHistBin=1; cotalphaHistBin<=rescotAlpha_binN; ++cotalphaHistBin)
    {
        for (unsigned cotbetaHistBin=1; cotbetaHistBin<=rescotBeta_binN; ++cotbetaHistBin)
        {
            unsigned int singleBigPixelHistN = 1*100000 + cotalphaHistBin*100 + cotbetaHistBin;
            theXHistos[singleBigPixelHistN] = new SimpleHistogramGenerator((TH1F*)theBigPixelResolutionFile->Get(Form("DQMData/clustBPIX/hx%u",singleBigPixelHistN)));
            theYHistos[singleBigPixelHistN] = new SimpleHistogramGenerator((TH1F*)theBigPixelResolutionFile->Get(Form("DQMData/clustBPIX/hy%u",singleBigPixelHistN)));

            unsigned int singlePixelHistN = 1*10000 + cotbetaHistBin*10 + cotalphaHistBin;
            theXHistos[singlePixelHistN] = new SimpleHistogramGenerator((TH1F*)theRegularPixelResolutionFile->Get(Form("hx%u",singlePixelHistN)));
            theYHistos[singlePixelHistN] = new SimpleHistogramGenerator((TH1F*)theRegularPixelResolutionFile->Get(Form("hy%u",singlePixelHistN)));
            
            for(unsigned qbinBin=1; qbinBin<=resqbin_binN; ++qbinBin )
            {
                unsigned int edgePixelHistN = cotalphaHistBin*1000 + cotbetaHistBin*10 + qbinBin;
                theXHistos[edgePixelHistN] = new SimpleHistogramGenerator((TH1F*)theEdgePixelResolutionFile->Get(Form("DQMData/clustBPIX/hx0%u",edgePixelHistN)));
                theYHistos[edgePixelHistN] = new SimpleHistogramGenerator((TH1F*)theEdgePixelResolutionFile->Get(Form("DQMData/clustBPIX/hy0%u",edgePixelHistN)));
                
                unsigned int multiPixelBigHistN = 1*1000000 + 1*100000 + cotalphaHistBin*1000 + cotbetaHistBin * 10 + qbinBin;
                theXHistos[multiPixelBigHistN] = new SimpleHistogramGenerator((TH1F*)theBigPixelResolutionFile->Get(Form("DQMData/clustBPIX/hx%u",multiPixelBigHistN)));
                theYHistos[multiPixelBigHistN] = new SimpleHistogramGenerator((TH1F*)theBigPixelResolutionFile->Get(Form("DQMData/clustBPIX/hy%u",multiPixelBigHistN)));
                
                unsigned int multiPixelHistN = 1*100000 + 1*10000 + cotbetaHistBin*100 + cotalphaHistBin*10 + qbinBin;
                theXHistos[multiPixelHistN] = new SimpleHistogramGenerator((TH1F*)theRegularPixelResolutionFile->Get(Form("hx%u",multiPixelHistN)));
                theYHistos[multiPixelHistN] = new SimpleHistogramGenerator((TH1F*)theRegularPixelResolutionFile->Get(Form("hy%u",multiPixelHistN)));
            }
        }
    }
}