PixelForwardTemplateSmearerPlugin Class Reference

Inheritance diagram for PixelForwardTemplateSmearerPlugin:

Public Member Functions
	PixelForwardTemplateSmearerPlugin (const std::string &name, const edm::ParameterSet &config, edm::ConsumesCollector &consumesCollector)
virtual	~PixelForwardTemplateSmearerPlugin ()
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)
virtual TrackingRecHitProductPtr	process (TrackingRecHitProductPtr product) const
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
virtual	~PixelTemplateSmearerBase ()
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	initializeForward ()

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 45 of file PixelForwardTemplateSmearerPlugin.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 61 of file PixelForwardTemplateSmearerPlugin.cc.

References Exception, edm::FileInPath::fullPath(), edm::ParameterSet::getParameter(), initializeForward(), 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 = true;

    theEdgePixelResolutionFileName = config.getParameter<string>( "EdgePixelForwardResolutionFile" );
    theEdgePixelResolutionFile = std::make_unique<TFile>( edm::FileInPath( theEdgePixelResolutionFileName ).fullPath().c_str()  ,"READ");
    theBigPixelResolutionFileName = config.getParameter<string>( "BigPixelForwardResolutionFile" );
    theBigPixelResolutionFile = std::make_unique<TFile>( edm::FileInPath( theBigPixelResolutionFileName ).fullPath().c_str()  ,"READ");
    theRegularPixelResolutionFileName = config.getParameter<string>( "RegularPixelForwardResolutionFile" );
    theRegularPixelResolutionFile = std::make_unique<TFile>( edm::FileInPath( theRegularPixelResolutionFileName ).fullPath().c_str()  ,"READ");
    
    theMergingProbabilityFileName = config.getParameter<string>( "MergingProbabilityForwardFile" );
    theMergingProbabilityFile =std::make_unique<TFile>( edm::FileInPath( theMergingProbabilityFileName ).fullPath().c_str()  ,"READ");
    theMergedPixelResolutionXFileName = config.getParameter<string>( "MergedPixelForwardResolutionXFile" );
    theMergedPixelResolutionXFile = std::make_unique<TFile>( edm::FileInPath( theMergedPixelResolutionXFileName ).fullPath().c_str()  ,"READ");
    theMergedPixelResolutionYFileName = config.getParameter<string>( "MergedPixelForwardResolutionYFile" );
    theMergedPixelResolutionYFile = std::make_unique<TFile>( edm::FileInPath( theMergedPixelResolutionYFileName ).fullPath().c_str()  ,"READ");

    initializeForward();
    
    
    
    if( ! SiPixelTemplate::pushfile(templateId, thePixelTemp_) )
    {
        throw cms::Exception("PixelForwardTemplateSmearerPlugin:") <<"SiPixel Forward Template Not Loaded Correctly!";
    }
}

PixelForwardTemplateSmearerPlugin::~PixelForwardTemplateSmearerPlugin ( )

virtual

Definition at line 96 of file PixelForwardTemplateSmearerPlugin.cc.

{
}

Member Function Documentation

void PixelForwardTemplateSmearerPlugin::initializeForward ( )

private

Definition at line 101 of file PixelForwardTemplateSmearerPlugin.cc.

References 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 PixelForwardTemplateSmearerPlugin().

{
    rescotAlpha_binMin = 0.1;
    rescotAlpha_binWidth = 0.1;
    rescotAlpha_binN = 4;
    rescotBeta_binMin = 0.;
    rescotBeta_binWidth = 0.15;
    rescotBeta_binN = 4;
    resqbin_binMin = 0;
    resqbin_binWidth = 1;
    resqbin_binN = 4;

    // Initialize the forward histos once and for all, and prepare the random generation
    for ( unsigned cotalphaHistBin=1; cotalphaHistBin<=rescotAlpha_binN; ++cotalphaHistBin)
    {
        for ( unsigned cotbetaHistBin=1; cotbetaHistBin<=rescotBeta_binN; ++cotbetaHistBin)
        {
            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/clustFPIX/fhx0%u",edgePixelHistN)));
                theYHistos[edgePixelHistN] = new SimpleHistogramGenerator((TH1F*)theEdgePixelResolutionFile->Get(Form("DQMData/clustFPIX/fhy0%u",edgePixelHistN)));
                
                unsigned int PixelHistN = 10000 + cotbetaHistBin*100 +  cotalphaHistBin*10 +  qbinBin;
                theXHistos[PixelHistN] = new SimpleHistogramGenerator((TH1F*) theRegularPixelResolutionFile->Get(Form("hx0%u",PixelHistN)));
                theYHistos[PixelHistN] = new SimpleHistogramGenerator((TH1F*) theRegularPixelResolutionFile->Get(Form("hy0%u",PixelHistN)));
            }
        }
    }

    for ( unsigned cotalphaHistBin=1; cotalphaHistBin<=rescotAlpha_binN; ++cotalphaHistBin)
    {
        for ( unsigned cotbetaHistBin=1; cotbetaHistBin<=rescotBeta_binN; ++cotbetaHistBin)
        {
            unsigned int SingleBigPixelHistN = 100000 + cotalphaHistBin*100 + cotbetaHistBin;
            theXHistos[SingleBigPixelHistN] = new SimpleHistogramGenerator((TH1F*)theBigPixelResolutionFile->Get(Form("DQMData/clustFPIX/fhx%u",SingleBigPixelHistN)));
            theYHistos[SingleBigPixelHistN] = new SimpleHistogramGenerator((TH1F*)theBigPixelResolutionFile->Get(Form("DQMData/clustFPIX/fhy%u",SingleBigPixelHistN)));
            
            unsigned int SinglePixelHistN = cotbetaHistBin*10 + cotalphaHistBin;
            theXHistos[SinglePixelHistN]  = new SimpleHistogramGenerator((TH1F*)theRegularPixelResolutionFile->Get(Form("hx000%u",SinglePixelHistN)));
            theYHistos[SinglePixelHistN]  = new SimpleHistogramGenerator((TH1F*)theRegularPixelResolutionFile->Get(Form("hy000%u",SinglePixelHistN)));
        }
    }
}