#include <MuonCaloCleanerByDistance.h>

Inheritance diagram for MuonCaloCleanerByDistance:

Public Member Functions
	MuonCaloCleanerByDistance (const edm::ParameterSet &)

	~MuonCaloCleanerByDistance ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef std::map< uint32_t, float >	detIdToFloatMap

Private Member Functions
void	fillEnergyDepositMap (const reco::Candidate &muon, const detIdToFloatMap &, detIdToFloatMap &, double &, double &)

virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
DetNaming	detNaming_

std::map< std::string, double >	energyDepositCorrection_

std::string	moduleLabel_

edm::InputTag	srcDistanceMapMuMinus_

edm::InputTag	srcDistanceMapMuPlus_

edm::InputTag	srcSelectedMuons_

int	verbosity_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Produce collections of calorimeter recHits from which energy deposits of muons are subtracted.

This module attributes to the muon energy proportional to the distance traversed by the muon through each calorimeter cell.

Author: Tomasz Maciej Frueboes; Christian Veelken, LLR

Version

Revision:: 1.4

Id:: MuonCaloCleanerByDistance.h,v 1.4 2013/02/05 20:01:19 veelken Exp

Definition at line 36 of file MuonCaloCleanerByDistance.h.

Member Typedef Documentation

typedef std::map<uint32_t, float> MuonCaloCleanerByDistance::detIdToFloatMap

private

Definition at line 45 of file MuonCaloCleanerByDistance.h.

Constructor & Destructor Documentation

MuonCaloCleanerByDistance::MuonCaloCleanerByDistance ( const edm::ParameterSet & cfg )

explicit

Definition at line 11 of file MuonCaloCleanerByDistance.cc.

References energyDepositCorrection_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNamesForType(), and verbosity_.

   : moduleLabel_(cfg.getParameter<std::string>("@module_label")), 
     srcSelectedMuons_(cfg.getParameter<edm::InputTag>("muons")),
     srcDistanceMapMuPlus_(cfg.getParameter<edm::InputTag>("distanceMapMuPlus")),
     srcDistanceMapMuMinus_(cfg.getParameter<edm::InputTag>("distanceMapMuMinus"))
 {
   edm::ParameterSet cfgEnergyDepositCorrection = cfg.getParameter<edm::ParameterSet>("energyDepositCorrection");
   typedef std::vector<std::string> vstring;
   vstring detNames = cfgEnergyDepositCorrection.getParameterNamesForType<double>();
   for ( vstring::const_iterator detName = detNames.begin();
     detName != detNames.end(); ++detName ) {
     energyDepositCorrection_[*detName] = cfgEnergyDepositCorrection.getParameter<double>(*detName);
   }
 
   verbosity_ = ( cfg.exists("verbosity") ) ?
     cfg.getParameter<int>("verbosity") : 0;
 
   // maps of detId to expected energy deposits of muon
   produces<detIdToFloatMap>("energyDepositsMuPlus");
   produces<double>("totalDistanceMuPlus");
   produces<double>("totalEnergyDepositMuPlus");
   produces<detIdToFloatMap>("energyDepositsMuMinus");
   produces<double>("totalDistanceMuMinus");
   produces<double>("totalEnergyDepositMuMinus");
 }

MuonCaloCleanerByDistance::~MuonCaloCleanerByDistance ( )

Definition at line 37 of file MuonCaloCleanerByDistance.cc.

 {
 // nothing to be done yet...
 }

Member Function Documentation

void MuonCaloCleanerByDistance::fillEnergyDepositMap	(	const reco::Candidate &	muon,
		const detIdToFloatMap &	distanceMap,
		detIdToFloatMap &	energyDepositMap,
		double &	totalDistance,
		double &	totalEnergyDeposit
	)

private

Definition at line 77 of file MuonCaloCleanerByDistance.cc.

References DENSITY_BRASS, DENSITY_IRON, DENSITY_PBWO4, DetId::det(), detNaming_, HLT_FULL_cff::distance, DetId::Ecal, energyDepositCorrection_, Exception, getDeDxForPbWO4(), DetNaming::getKey(), DetId::Hcal, HcalOuter, relval_steps::key, reco::Candidate::p(), DetId::rawId(), rho, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), reco::Candidate::theta(), and theta().

Referenced by produce().

 {
   for ( detIdToFloatMap::const_iterator rawDetId_and_distance = distanceMap.begin();
     rawDetId_and_distance != distanceMap.end(); ++rawDetId_and_distance ) {
     DetId detId(rawDetId_and_distance->first);
 
     std::string key = detNaming_.getKey(detId);
 
     if ( energyDepositCorrection_.find(key) == energyDepositCorrection_.end() )
       throw cms::Exception("MuonCaloCleanerByDistance") 
     << "No energy deposit correction defined for detId = " << detId.rawId() << " (key = " << key << ") !!\n";
 
     double distance = rawDetId_and_distance->second;
     double energyDepositCorrection_value = energyDepositCorrection_[key];
 
     double dEdx = 0.;
     double rho = 0.;
     switch ( detId.det() ) {
     case DetId::Ecal:
       dEdx = getDeDxForPbWO4(muon.p());
       rho = DENSITY_PBWO4;
       break;
     case DetId::Hcal:
       // AB: We don't have a dedx curve for the HCAL. Use the PbWO4 one as an approximation,
       // the correction factors should be determined with respect to the PbWO4 curve.
       dEdx = getDeDxForPbWO4(muon.p());
       if ( detId.subdetId() == HcalOuter ) {
         rho = DENSITY_IRON; // iron coil and return yoke
         // HO uses magnet coil as additional absorber, add to flight distance:
         const double theta = muon.theta();
         distance += 31.2 / sin(theta); // 31.2cm is dr of cold mass of the magnet coil
       } else {
         rho = DENSITY_BRASS; // brass absorber
       }
       break;
     default:
       throw cms::Exception("MuonCaloCleanerByDistance") 
     << "Unknown detector type: " << key << ", detId = " << static_cast<unsigned int>(detId);
     }
 
     double energyDeposit = distance*dEdx*rho*energyDepositCorrection_value;
     energyDepositMap[rawDetId_and_distance->first] += energyDeposit;
     totalDistance += distance;
     totalEnergyDeposit += energyDeposit;
   }
 }

void MuonCaloCleanerByDistance::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 42 of file MuonCaloCleanerByDistance.cc.

References gather_cfg::cout, fillEnergyDepositMap(), edm::Event::getByLabel(), getSelMuons(), getTheMuMinus(), getTheMuPlus(), edm::RefToBase< T >::isNonnull(), moduleLabel_, edm::Event::put(), srcDistanceMapMuMinus_, srcDistanceMapMuPlus_, srcSelectedMuons_, and verbosity_.

 {
   std::auto_ptr<detIdToFloatMap> energyDepositsMuPlus(new detIdToFloatMap());
   std::auto_ptr<detIdToFloatMap> energyDepositsMuMinus(new detIdToFloatMap());
 
   edm::Handle<detIdToFloatMap> distanceMapMuPlus;
   evt.getByLabel(srcDistanceMapMuPlus_, distanceMapMuPlus);
   edm::Handle<detIdToFloatMap> distanceMapMuMinus;
   evt.getByLabel(srcDistanceMapMuMinus_, distanceMapMuMinus);
 
   std::vector<reco::CandidateBaseRef> selMuons = getSelMuons(evt, srcSelectedMuons_);
   const reco::CandidateBaseRef muPlus  = getTheMuPlus(selMuons);
   const reco::CandidateBaseRef muMinus = getTheMuMinus(selMuons);
 
   std::auto_ptr<double> totalDistanceMuPlus(new double(0.));
   std::auto_ptr<double> totalEnergyDepositMuPlus(new double(0.));
   if ( muPlus.isNonnull() ) fillEnergyDepositMap(*muPlus, *distanceMapMuPlus, *energyDepositsMuPlus, *totalDistanceMuPlus, *totalEnergyDepositMuPlus);
   std::auto_ptr<double> totalDistanceMuMinus(new double(0.));
   std::auto_ptr<double> totalEnergyDepositMuMinus(new double(0.));
   if ( muMinus.isNonnull() ) fillEnergyDepositMap(*muMinus, *distanceMapMuMinus, *energyDepositsMuMinus, *totalDistanceMuMinus, *totalEnergyDepositMuMinus);
 
   if ( verbosity_ ) {
     std::cout << "<MuonCaloCleanerByDistance::produce (" << moduleLabel_ << ")>:" << std::endl;
     std::cout << " mu+: distance = " << (*totalDistanceMuPlus) << ", expected(EnergyDeposits) = " << (*totalEnergyDepositMuPlus) << std::endl;
     std::cout << " mu-: distance = " << (*totalDistanceMuMinus) << ", expected(EnergyDeposits) = " << (*totalEnergyDepositMuMinus) << std::endl;
   }
 
   evt.put(energyDepositsMuPlus, "energyDepositsMuPlus");
   evt.put(totalDistanceMuPlus, "totalDistanceMuPlus");
   evt.put(totalEnergyDepositMuPlus, "totalEnergyDepositMuPlus");
   evt.put(energyDepositsMuMinus, "energyDepositsMuMinus");
   evt.put(totalDistanceMuMinus, "totalDistanceMuMinus");
   evt.put(totalEnergyDepositMuMinus, "totalEnergyDepositMuMinus");
 }