#include <L3MuonCombinedRelativeIsolationProducer.h>

Inheritance diagram for L3MuonCombinedRelativeIsolationProducer:

Public Member Functions
virtual void	beginJob ()
	initialisation More...

	L3MuonCombinedRelativeIsolationProducer (const edm::ParameterSet &)
	constructor with config More...

virtual void	produce (edm::Event &, const edm::EventSetup &)
	Produce isolation maps. More...

virtual	~L3MuonCombinedRelativeIsolationProducer ()
	destructor More...

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
reco::isodeposit::IsoDepositExtractor *	caloExtractor

bool	optOutputIsoDeposits

bool	printDebug

bool	theApplyCutsORmaxNTracks
	apply or not the maxN cut on top of the sumPt (or nominall eff) < cuts More...

edm::InputTag	theCaloDepsLabel

edm::ParameterSet	theConfig

muonisolation::Cuts	theCuts

int	theMaxNTracks

edm::InputTag	theMuonCollectionLabel

double	theTrackPt_Min

reco::isodeposit::IsoDepositExtractor *	trkExtractor

bool	useRhoCorrectedCaloDeps

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

typedef WorkerT< EDProducer >	WorkerType

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

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::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 18 of file L3MuonCombinedRelativeIsolationProducer.h.

Constructor & Destructor Documentation

L3MuonCombinedRelativeIsolationProducer::L3MuonCombinedRelativeIsolationProducer ( const edm::ParameterSet & par )

constructor with config

Definition at line 36 of file L3MuonCombinedRelativeIsolationProducer.cc.

References LogDebug, optOutputIsoDeposits, and useRhoCorrectedCaloDeps.

                                                                                                         :
   theConfig(par),
   theMuonCollectionLabel(par.getParameter<InputTag>("inputMuonCollection")),
   optOutputIsoDeposits(par.getParameter<bool>("OutputMuIsoDeposits")),
   useRhoCorrectedCaloDeps(par.existsAs<bool>("UseRhoCorrectedCaloDeposits") ? 
               par.getParameter<bool>("UseRhoCorrectedCaloDeposits") : false),
   theCaloDepsLabel(par.existsAs<InputTag>("CaloDepositsLabel") ? 
            par.getParameter<InputTag>("CaloDepositsLabel") :
            InputTag("hltL3CaloMuonCorrectedIsolations")),
   caloExtractor(0),
   trkExtractor(0),
   theTrackPt_Min(-1),
   printDebug (par.getParameter<bool>("printDebug"))
   {
   LogDebug("RecoMuon|L3MuonCombinedRelativeIsolationProducer")<<" L3MuonCombinedRelativeIsolationProducer CTOR";
 
   if (optOutputIsoDeposits) {
     produces<reco::IsoDepositMap>("trkIsoDeposits");
     if( useRhoCorrectedCaloDeps==false ) // otherwise, calo deposits have been previously computed
       produces<reco::IsoDepositMap>("caloIsoDeposits");
     //produces<std::vector<double> >("combinedRelativeIsoDeposits");
     produces<edm::ValueMap<double> >("combinedRelativeIsoDeposits");
   }
   produces<edm::ValueMap<bool> >();
 
 }

L3MuonCombinedRelativeIsolationProducer::~L3MuonCombinedRelativeIsolationProducer ( )

virtual

destructor

Definition at line 64 of file L3MuonCombinedRelativeIsolationProducer.cc.

References caloExtractor, LogDebug, and trkExtractor.

                                                                                  {
   LogDebug("RecoMuon|L3MuonCombinedRelativeIsolationProducer")<<" L3MuonCombinedRelativeIsolationProducer DTOR";
   if (caloExtractor) delete caloExtractor;
   if (trkExtractor) delete trkExtractor;
 }

Member Function Documentation

void L3MuonCombinedRelativeIsolationProducer::beginJob ( void )

virtual

initialisation

test cutsName only. The depositType is informational only (has not been used so far) [VK]

Reimplemented from edm::EDProducer.

Definition at line 70 of file L3MuonCombinedRelativeIsolationProducer.cc.

References caloExtractor, L3NominalEfficiencyConfigurator::cuts(), reco::get(), edm::ParameterSet::getParameter(), LogTrace, muonisolation::Cuts::print(), theApplyCutsORmaxNTracks, theConfig, theCuts, theMaxNTracks, theTrackPt_Min, trkExtractor, and useRhoCorrectedCaloDeps.

 {
 
   //
   // Extractor
   //
   // Calorimeters (ONLY if not previously computed)
   //
   if( useRhoCorrectedCaloDeps==false ) {
     edm::ParameterSet caloExtractorPSet = theConfig.getParameter<edm::ParameterSet>("CaloExtractorPSet");
   
     theTrackPt_Min = theConfig.getParameter<double>("TrackPt_Min");
     std::string caloExtractorName = caloExtractorPSet.getParameter<std::string>("ComponentName");
     caloExtractor = IsoDepositExtractorFactory::get()->create( caloExtractorName, caloExtractorPSet);
     //std::string caloDepositType = caloExtractorPSet.getUntrackedParameter<std::string>("DepositLabel"); // N.B. Not used in the following!
   }
 
   // Tracker
   //
   edm::ParameterSet trkExtractorPSet = theConfig.getParameter<edm::ParameterSet>("TrkExtractorPSet");
 
   std::string trkExtractorName = trkExtractorPSet.getParameter<std::string>("ComponentName");
   trkExtractor = IsoDepositExtractorFactory::get()->create( trkExtractorName, trkExtractorPSet);
   //std::string trkDepositType = trkExtractorPSet.getUntrackedParameter<std::string>("DepositLabel"); // N.B. Not used in the following!
 
 
   
   //
   // Cuts for track isolation
   //
   edm::ParameterSet cutsPSet = theConfig.getParameter<edm::ParameterSet>("CutsPSet");
   std::string cutsName = cutsPSet.getParameter<std::string>("ComponentName");
   if (cutsName == "SimpleCuts") {
     theCuts = Cuts(cutsPSet);
   } 
   else if (
        //        (cutsName== "L3NominalEfficiencyCuts_PXLS" && depositType=="PXLS")
        //     || (cutsName== "L3NominalEfficiencyCuts_TRKS" && depositType=="TRKS") 
        (cutsName== "L3NominalEfficiencyCuts_PXLS" )
        || (cutsName== "L3NominalEfficiencyCuts_TRKS") ) {
     theCuts = L3NominalEfficiencyConfigurator(cutsPSet).cuts();
   } 
   else {
     LogError("L3MuonCombinedRelativeIsolationProducer::beginJob")
       <<"cutsName: "<<cutsPSet<<" is not recognized:"
       <<" theCuts not set!";
   }
   LogTrace("")<< theCuts.print();
 
   // (kludge) additional cut on the number of tracks
   theMaxNTracks = cutsPSet.getParameter<int>("maxNTracks");
   theApplyCutsORmaxNTracks = cutsPSet.getParameter<bool>("applyCutsORmaxNTracks");
 }

void L3MuonCombinedRelativeIsolationProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

virtual

Produce isolation maps.

Implements edm::EDProducer.

Definition at line 125 of file L3MuonCombinedRelativeIsolationProducer.cc.

References caloExtractor, muonisolation::Cuts::CutSpec::conesize, prof2calltree::count, gather_cfg::cout, GOODCOLL_filter_cfg::cut, reco::isodeposit::IsoDepositExtractor::deposit(), reco::IsoDeposit::depositAndCountWithin(), reco::IsoDeposit::depositWithin(), reco::TrackBase::eta(), muonisolation::Cuts::CutSpec::etaRange, edm::helper::Filler< Map >::fill(), i, edm::helper::Filler< Map >::insert(), muonisolation::Range< T >::max(), metname, muonisolation::Range< T >::min(), RPCpg::mu, patZpeak::muons, optOutputIsoDeposits, reco::IsoDeposit::print(), printDebug, reco::TrackBase::pt(), query::result, theApplyCutsORmaxNTracks, theCaloDepsLabel, theCuts, theMaxNTracks, theMuonCollectionLabel, theTrackPt_Min, muonisolation::Cuts::CutSpec::threshold, trkExtractor, useRhoCorrectedCaloDeps, and reco::IsoDeposit::veto().

                                                                                                {
   std::string metname = "RecoMuon|L3MuonCombinedRelativeIsolationProducer";
   
   if (printDebug) std::cout  <<" L3 Muon Isolation producing..."
             <<" BEGINING OF EVENT " <<"================================" <<std::endl;
 
   // Take the SA container
   if (printDebug) std::cout  <<" Taking the muons: "<<theMuonCollectionLabel << std::endl;
   Handle<TrackCollection> muons;
   event.getByLabel(theMuonCollectionLabel,muons);
 
   // Take calo deposits with rho corrections (ONLY if previously computed)
   Handle< edm::ValueMap<float> > caloDepWithCorrMap;
   if( useRhoCorrectedCaloDeps )
     event.getByLabel(theCaloDepsLabel, caloDepWithCorrMap);
 
   std::auto_ptr<reco::IsoDepositMap> caloDepMap( new reco::IsoDepositMap());
   std::auto_ptr<reco::IsoDepositMap> trkDepMap( new reco::IsoDepositMap());
 
   std::auto_ptr<edm::ValueMap<bool> > comboIsoDepMap( new edm::ValueMap<bool> ());
 
   //std::auto_ptr<std::vector<double> > combinedRelativeDeps(new std::vector<double>());
   std::auto_ptr<edm::ValueMap<double> > combinedRelativeDepMap(new edm::ValueMap<double>());
 
   
   //
   // get Vetos and deposits
   //
   unsigned int nMuons = muons->size();
 
   IsoDeposit::Vetos trkVetos(nMuons);  
   std::vector<IsoDeposit> trkDeps(nMuons);
   
 
   // IsoDeposit::Vetos caloVetos(nMuons);  
   // std::vector<IsoDeposit> caloDeps(nMuons);
   // std::vector<float> caloCorrDeps(nMuons, 0.);  // if calo deposits with corrections available
 
   IsoDeposit::Vetos caloVetos;  
   std::vector<IsoDeposit> caloDeps;
   std::vector<float> caloCorrDeps;  // if calo deposits with corrections available
 
   if(useRhoCorrectedCaloDeps) {
     caloCorrDeps.resize(nMuons, 0.);
   }
   else {
     caloVetos.resize(nMuons);
     caloDeps.resize(nMuons);
   }
   
   std::vector<double> combinedRelativeDeps(nMuons, 0.);
   std::vector<bool> combinedRelativeIsos(nMuons, false);
 
   for (unsigned int i=0; i<nMuons; i++) {
     
     TrackRef mu(muons,i);
     
     trkDeps[i] = trkExtractor->deposit(event, eventSetup, *mu);
     trkVetos[i] = trkDeps[i].veto();
 
     if( useRhoCorrectedCaloDeps ) {
       caloCorrDeps[i] = (*caloDepWithCorrMap)[mu];
     }
     else {
       caloDeps[i] = caloExtractor->deposit(event, eventSetup, *mu);
       caloVetos[i] = caloDeps[i].veto();
     }
 
   }
 
   //
   // add here additional vetos
   //
   //.....
 
   //
   // actual cut step
   //
 
   if (printDebug) std::cout  << "Looping over deposits...." << std::endl;
   for(unsigned int iMu=0; iMu < nMuons; ++iMu){
 
     if (printDebug) std::cout  << "Muon number = " << iMu << std::endl;
     const reco::Track* mu = &(*muons)[iMu];
 
     // cuts 
     const Cuts::CutSpec & cut = theCuts( mu->eta());
 
 
     if (printDebug) std::cout << "CUTDEBUG: Muon eta = " << mu->eta() << std::endl
                               << "CUTDEBUG: Muon pt  = " <<  mu->pt() << std::endl
                               << "CUTDEBUG: minEta   = " << cut.etaRange.min() << std::endl
                               << "CUTDEBUG: maxEta   = " << cut.etaRange.max() << std::endl
                               << "CUTDEBUG: consize  = " << cut.conesize << std::endl
                               << "CUTDEBUG: thresho  = " << cut.threshold << std::endl;
     
     const IsoDeposit & trkDeposit = trkDeps[iMu];
     if (printDebug) std::cout  << trkDeposit.print();
     std::pair<double, int> trkIsoSumAndCount = trkDeposit.depositAndCountWithin(cut.conesize, trkVetos, theTrackPt_Min);
 
     double caloIsoSum = 0.;
     if( useRhoCorrectedCaloDeps ) {
       caloIsoSum = caloCorrDeps[iMu];
       if(caloIsoSum<0.) caloIsoSum = 0.;
       if(printDebug) std::cout << "Rho-corrected calo deposit (min. 0) = " << caloIsoSum << std::endl;
     }
     else {
       const IsoDeposit & caloDeposit = caloDeps[iMu];
       if (printDebug) std::cout  << caloDeposit.print();
       caloIsoSum = caloDeposit.depositWithin(cut.conesize, caloVetos);
     }
 
     double trkIsoSum = trkIsoSumAndCount.first;
     int count = trkIsoSumAndCount.second;
 
     double muPt = mu->pt();
     if( muPt<1. ) muPt = 1.;
     double combinedRelativeDeposit = ((trkIsoSum + caloIsoSum ) / muPt);
     bool result = ( combinedRelativeDeposit < cut.threshold); 
     if (theApplyCutsORmaxNTracks ) result |= count <= theMaxNTracks;
     if (printDebug) std::cout  <<"  trk dep in cone:  " << trkIsoSum << "  with count "<<count <<std::endl
               <<"  calo dep in cone: " << caloIsoSum << std::endl
               <<"  muPt: " << muPt << std::endl
               <<"  relIso:  " <<combinedRelativeDeposit  << std::endl
               <<"  is isolated: "<<result << std::endl;
 
     combinedRelativeIsos[iMu] = result;
     //combinedRelativeDeps->push_back(combinedRelativeDeposit);
     combinedRelativeDeps[iMu] = combinedRelativeDeposit;
   }
 
   //
   // store
   //
   if (optOutputIsoDeposits){
     
     reco::IsoDepositMap::Filler depFillerTrk(*trkDepMap);
     depFillerTrk.insert(muons, trkDeps.begin(), trkDeps.end());
     depFillerTrk.fill();
     event.put(trkDepMap, "trkIsoDeposits");
 
     if( useRhoCorrectedCaloDeps==false ) {
       reco::IsoDepositMap::Filler depFillerCalo(*caloDepMap);
       depFillerCalo.insert(muons, caloDeps.begin(), caloDeps.end());
       depFillerCalo.fill();
       event.put(caloDepMap, "caloIsoDeposits");
     }
 
     //event.put(combinedRelativeDeps, "combinedRelativeIsoDeposits");
     edm::ValueMap<double>::Filler depFillerCombRel(*combinedRelativeDepMap);
     depFillerCombRel.insert(muons, combinedRelativeDeps.begin(), combinedRelativeDeps.end());
     depFillerCombRel.fill();
     event.put(combinedRelativeDepMap, "combinedRelativeIsoDeposits");
 
   }
   edm::ValueMap<bool>::Filler isoFiller(*comboIsoDepMap);
   isoFiller.insert(muons, combinedRelativeIsos.begin(), combinedRelativeIsos.end());
   isoFiller.fill();
   event.put(comboIsoDepMap);
 
   if (printDebug) std::cout  <<" END OF EVENT " <<"================================";
 }