reco::Mustache Class Reference

#include <Mustache.h>

Public Member Functions
void	FillMustacheVar (const std::vector< CaloCluster > &clusters)
int	InsideMust ()
float	LowestMustClust ()
	Mustache (const EcalMustacheSCParameters *mustache_params)
void	MustacheClust (const std::vector< CaloCluster > &clusters, std::vector< unsigned int > &insideMust, std::vector< unsigned int > &outsideMust)
float	MustacheE ()
float	MustacheEOut ()
float	MustacheEtOut ()
void	MustacheID (const CaloClusterPtrVector &clusters, int &nclusters, float &EoutsideMustache)
void	MustacheID (const std::vector< const CaloCluster *> &, int &nclusers, float &EoutsideMustache)
void	MustacheID (const reco::SuperCluster &sc, int &nclusters, float &EoutsideMustache)
int	OutsideMust ()

Private Member Functions
template<class RandomAccessPtrIterator >
void	MustacheID (const RandomAccessPtrIterator &, const RandomAccessPtrIterator &, int &nclusters, float &EoutsideMustache)

Private Attributes
float	Energy_In_Mustache_
float	Energy_Outside_Mustache_
float	Et_Outside_Mustache_
int	excluded_
int	included_
float	LowestClusterEInMustache_
const EcalMustacheSCParameters *	mustache_params_

Detailed Description

Definition at line 28 of file Mustache.h.

Constructor & Destructor Documentation

Mustache()

reco::Mustache::Mustache

(

const EcalMustacheSCParameters *

mustache_params

)

Definition at line 89 of file Mustache.cc.

: mustache_params_(mustache_params) {}

Member Function Documentation

FillMustacheVar()

void reco::Mustache::FillMustacheVar

(

const std::vector< CaloCluster > &

clusters

)

Definition at line 178 of file Mustache.cc.

References bsc_activity_cfg::clusters, hcalRecHitTable_cff::energy, Energy_In_Mustache_, Energy_Outside_Mustache_, Et_Outside_Mustache_, excluded_, mps_fire::i, included_, LowestClusterEInMustache_, MustacheClust(), position, funct::sin(), and theta().

                                                                       {
    Energy_In_Mustache_ = 0;
    Energy_Outside_Mustache_ = 0;
    LowestClusterEInMustache_ = 0;
    excluded_ = 0;
    included_ = 0;
    std::multimap<float, unsigned int> OrderedClust;
    std::vector<unsigned int> insideMust;
    std::vector<unsigned int> outsideMust;
    MustacheClust(clusters, insideMust, outsideMust);
    included_ = insideMust.size();
    excluded_ = outsideMust.size();
    for (unsigned int i = 0; i < insideMust.size(); ++i) {
      unsigned int index = insideMust[i];
      Energy_In_Mustache_ = clusters[index].energy() + Energy_In_Mustache_;
      OrderedClust.insert(make_pair(clusters[index].energy(), index));
    }
    for (unsigned int i = 0; i < outsideMust.size(); ++i) {
      unsigned int index = outsideMust[i];
      Energy_Outside_Mustache_ = clusters[index].energy() + Energy_Outside_Mustache_;
      Et_Outside_Mustache_ = clusters[index].energy() * sin(clusters[index].position().theta()) + Et_Outside_Mustache_;
    }
    std::multimap<float, unsigned int>::iterator it;
    it = OrderedClust.begin();
    unsigned int lowEindex = (*it).second;
    LowestClusterEInMustache_ = clusters[lowEindex].energy();
  }

InsideMust()

int reco::Mustache::InsideMust ( )

inline

Definition at line 46 of file Mustache.h.

References included_.

{ return included_; }

LowestMustClust()

float reco::Mustache::LowestMustClust ( )

inline

Definition at line 45 of file Mustache.h.

References LowestClusterEInMustache_.

{ return LowestClusterEInMustache_; }

MustacheClust()

void reco::Mustache::MustacheClust	(	const std::vector< CaloCluster > &	clusters,
		std::vector< unsigned int > &	insideMust,
		std::vector< unsigned int > &	outsideMust
	)

Definition at line 143 of file Mustache.cc.

References bsc_activity_cfg::clusters, MillePedeFileConverter_cfg::e, hcalRecHitTable_cff::energy, PVValHelper::eta, mps_fire::i, reco::MustacheKernel::inMustache(), dqmdumpme::k, and mustache_params_.

Referenced by FillMustacheVar().

                                                                     {
    unsigned int ncl = clusters.size();
    if (!ncl)
      return;

    //loop over all clusters to find the one with highest energy
    float emax = 0;
    int imax = -1;
    for (unsigned int i = 0; i < ncl; ++i) {
      float e = (clusters[i]).energy();
      if (e > emax) {
        emax = e;
        imax = i;
      }
    }

    if (imax < 0)
      return;
    float eta0 = (clusters[imax]).eta();
    float phi0 = (clusters[imax]).phi();

    for (unsigned int k = 0; k < ncl; k++) {
      bool inMust = MustacheKernel::inMustache(
          mustache_params_, eta0, phi0, (clusters[k]).energy(), (clusters[k]).eta(), (clusters[k]).phi());
      //return indices of Clusters outside the Mustache
      if (!(inMust)) {
        outsideMust.push_back(k);
      } else {  //return indices of Clusters inside the Mustache
        insideMust.push_back(k);
      }
    }
  }

MustacheE()

float reco::Mustache::MustacheE ( )

inline

Definition at line 42 of file Mustache.h.

References Energy_In_Mustache_.

{ return Energy_In_Mustache_; }

MustacheEOut()

float reco::Mustache::MustacheEOut ( )

inline

Definition at line 43 of file Mustache.h.

References Energy_Outside_Mustache_.

{ return Energy_Outside_Mustache_; }

MustacheEtOut()

float reco::Mustache::MustacheEtOut ( )

inline

Definition at line 44 of file Mustache.h.

References Et_Outside_Mustache_.

{ return Et_Outside_Mustache_; }

MustacheID() [1/4]

void reco::Mustache::MustacheID	(	const CaloClusterPtrVector &	clusters,
		int &	nclusters,
		float &	EoutsideMustache
	)

Definition at line 95 of file Mustache.cc.

References bsc_activity_cfg::clusters.

Referenced by PFElectronTranslator::createGsfElectrons(), PFPhotonTranslator::createPhotons(), and MustacheID().

                                                                                                         {
    MustacheID(clusters.begin(), clusters.end(), nclusters, EoutsideMustache);
  }

MustacheID() [2/4]

void reco::Mustache::MustacheID	(	const std::vector< const CaloCluster *> &	clusters,
		int &	nclusers,
		float &	EoutsideMustache
	)

Definition at line 99 of file Mustache.cc.

References bsc_activity_cfg::clusters, and MustacheID().

                                                                                                                  {
    MustacheID(clusters.cbegin(), clusters.cend(), nclusters, EoutsideMustache);
  }

MustacheID() [3/4]

void reco::Mustache::MustacheID	(	const reco::SuperCluster &	sc,
		int &	nclusters,
		float &	EoutsideMustache
	)

Definition at line 91 of file Mustache.cc.

References reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), and MustacheID().

                                                                                               {
    MustacheID(sc.clustersBegin(), sc.clustersEnd(), nclusters, EoutsideMustache);
  }

MustacheID() [4/4]

template<class RandomAccessPtrIterator >

void reco::Mustache::MustacheID ( const RandomAccessPtrIterator &  begin, const RandomAccessPtrIterator &  end, int &  nclusters, float &  EoutsideMustache  )

private

Definition at line 104 of file Mustache.cc.

References MillePedeFileConverter_cfg::e, reco::MustacheKernel::inMustache(), createfilelist::int, and mustache_params_.

                                                     {
    nclusters = 0;
    EoutsideMustache = 0;

    unsigned int ncl = end - begin;
    if (!ncl)
      return;

    //loop over all clusters to find the one with highest energy
    RandomAccessPtrIterator icl = begin;
    RandomAccessPtrIterator clmax = end;
    float emax = 0;
    for (; icl != end; ++icl) {
      const float e = (*icl)->energy();
      if (e > emax) {
        emax = e;
        clmax = icl;
      }
    }

    if (end == clmax)
      return;

    float eta0 = (*clmax)->eta();
    float phi0 = (*clmax)->phi();

    bool inMust = false;
    icl = begin;
    for (; icl != end; ++icl) {
      inMust = MustacheKernel::inMustache(mustache_params_, eta0, phi0, (*icl)->energy(), (*icl)->eta(), (*icl)->phi());

      nclusters += (int)!inMust;
      EoutsideMustache += (!inMust) * ((*icl)->energy());
    }
  }

OutsideMust()

int reco::Mustache::OutsideMust ( )

inline

Definition at line 47 of file Mustache.h.

References excluded_.

{ return excluded_; }

Member Data Documentation

Energy_In_Mustache_

float reco::Mustache::Energy_In_Mustache_

private

Definition at line 56 of file Mustache.h.

Referenced by FillMustacheVar(), and MustacheE().

Energy_Outside_Mustache_

float reco::Mustache::Energy_Outside_Mustache_

private

Definition at line 57 of file Mustache.h.

Referenced by FillMustacheVar(), and MustacheEOut().

Et_Outside_Mustache_

float reco::Mustache::Et_Outside_Mustache_

private

Definition at line 58 of file Mustache.h.

Referenced by FillMustacheVar(), and MustacheEtOut().

excluded_

int reco::Mustache::excluded_

private

Definition at line 60 of file Mustache.h.

Referenced by FillMustacheVar(), and OutsideMust().

included_

int reco::Mustache::included_

private

Definition at line 61 of file Mustache.h.

Referenced by FillMustacheVar(), and InsideMust().

LowestClusterEInMustache_

float reco::Mustache::LowestClusterEInMustache_

private

Definition at line 59 of file Mustache.h.

Referenced by FillMustacheVar(), and LowestMustClust().

mustache_params_

const EcalMustacheSCParameters* reco::Mustache::mustache_params_

private

Definition at line 62 of file Mustache.h.

Referenced by MustacheClust(), and MustacheID().