#include <MuonTCMETValueMapProducer.h>

Inheritance diagram for cms::MuonTCMETValueMapProducer:

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

	~MuonTCMETValueMapProducer ()

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

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	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
	EDConsumerBase ()

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

bool	isGoodCaloMuon (const reco::Muon *, const unsigned int)

bool	isGoodMuon (const reco::Muon *)

bool	isGoodTrack (const reco::Muon *)

bool	isValidVertex ()

int	nLayers (const reco::TrackRef)

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

class TVector3	propagateTrack (const reco::Muon *)

Private Attributes
edm::Handle< reco::BeamSpot >	beamSpot_h

edm::InputTag	beamSpotInputTag_

class MagneticField *	bField

double	d0cuta_

double	d0cutb_

bool	hasValidVertex

double	maxchi2_

double	maxchi2_tight_

double	maxd0cut_

double	maxeta_

double	maxpt_

double	maxpt_eta20_

double	maxpt_eta25_

double	maxPtErr_

double	maxPtErr_tight_

int	maxTrackAlgo_

double	minhits_

double	minhits_tight_

double	minpt_

double	muon_dptrel_

edm::Handle< reco::MuonCollection >	muon_h

double	muonchi2_

double	muond0_

double	muonDeltaR_

double	muoneta_

bool	muonGlobal_

double	muonhits_

edm::InputTag	muonInputTag_

int	muonMinValidStaHits_

double	muonpt_

bool	muonTracker_

int	nLayers_

int	nLayersTight_

class TH2D *	response_function

int	rfType_

TCMETAlgo *	tcmetAlgo_

std::vector< int >	trkAlgos_

std::vector< int >	trkQuality_

bool	useCaloMuons_

bool	usePvtxd0_

const reco::VertexCollection *	vertexColl

edm::Handle < reco::VertexCollection >	VertexHandle

edm::InputTag	vertexInputTag_

double	vertexMaxDZ_

int	vertexNdof_

double	vertexRho_

double	vertexZ_

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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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

Definition at line 38 of file MuonTCMETValueMapProducer.h.

Constructor & Destructor Documentation

MuonTCMETValueMapProducer::MuonTCMETValueMapProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 58 of file MuonTCMETValueMapProducer.cc.

References beamSpotInputTag_, d0cuta_, d0cutb_, edm::ParameterSet::getParameter(), maxchi2_, maxchi2_tight_, maxd0cut_, maxeta_, maxpt_, maxpt_eta20_, maxpt_eta25_, maxPtErr_, maxPtErr_tight_, maxTrackAlgo_, minhits_, minhits_tight_, minpt_, muon_dptrel_, muonchi2_, muond0_, muonDeltaR_, muoneta_, muonGlobal_, muonhits_, muonInputTag_, muonMinValidStaHits_, muonpt_, muonTracker_, nLayers_, nLayersTight_, response_function, rfType_, tcmetAlgo_, trkAlgos_, trkQuality_, useCaloMuons_, usePvtxd0_, vertexInputTag_, vertexMaxDZ_, vertexNdof_, vertexRho_, and vertexZ_.

                                                                                      {
   
     produces<edm::ValueMap<reco::MuonMETCorrectionData> > ("muCorrData");
 
     // get input collections
     muonInputTag_     = iConfig.getParameter<edm::InputTag>("muonInputTag"    );
     beamSpotInputTag_ = iConfig.getParameter<edm::InputTag>("beamSpotInputTag");
     vertexInputTag_   = iConfig.getParameter<edm::InputTag>("vertexInputTag");
 
     rfType_     = iConfig.getParameter<int>("rf_type");
 
     nLayers_                = iConfig.getParameter<int>      ("nLayers");
     nLayersTight_           = iConfig.getParameter<int>      ("nLayersTight");
     vertexNdof_             = iConfig.getParameter<int>      ("vertexNdof");
     vertexZ_                = iConfig.getParameter<double>   ("vertexZ");
     vertexRho_              = iConfig.getParameter<double>   ("vertexRho");
     vertexMaxDZ_            = iConfig.getParameter<double>   ("vertexMaxDZ");
     maxpt_eta20_            = iConfig.getParameter<double>   ("maxpt_eta20");
     maxpt_eta25_            = iConfig.getParameter<double>   ("maxpt_eta25");
 
     // get configuration parameters
     maxTrackAlgo_    = iConfig.getParameter<int>("trackAlgo_max");
     maxd0cut_        = iConfig.getParameter<double>("d0_max"       );
     minpt_           = iConfig.getParameter<double>("pt_min"       );
     maxpt_           = iConfig.getParameter<double>("pt_max"       );
     maxeta_          = iConfig.getParameter<double>("eta_max"      );
     maxchi2_         = iConfig.getParameter<double>("chi2_max"     );
     minhits_         = iConfig.getParameter<double>("nhits_min"    );
     maxPtErr_        = iConfig.getParameter<double>("ptErr_max"    );
 
     trkQuality_      = iConfig.getParameter<std::vector<int> >("track_quality");
     trkAlgos_        = iConfig.getParameter<std::vector<int> >("track_algos"  );
     maxchi2_tight_   = iConfig.getParameter<double>("chi2_max_tight");
     minhits_tight_   = iConfig.getParameter<double>("nhits_min_tight");
     maxPtErr_tight_  = iConfig.getParameter<double>("ptErr_max_tight");
     usePvtxd0_       = iConfig.getParameter<bool>("usePvtxd0");
     d0cuta_          = iConfig.getParameter<double>("d0cuta");
     d0cutb_          = iConfig.getParameter<double>("d0cutb");
 
     muon_dptrel_  = iConfig.getParameter<double>("muon_dptrel");
     muond0_     = iConfig.getParameter<double>("d0_muon"    );
     muonpt_     = iConfig.getParameter<double>("pt_muon"    );
     muoneta_    = iConfig.getParameter<double>("eta_muon"   );
     muonchi2_   = iConfig.getParameter<double>("chi2_muon"  );
     muonhits_   = iConfig.getParameter<double>("nhits_muon" );
     muonGlobal_   = iConfig.getParameter<bool>("global_muon");
     muonTracker_  = iConfig.getParameter<bool>("tracker_muon");
     muonDeltaR_ = iConfig.getParameter<double>("deltaR_muon");
     useCaloMuons_ = iConfig.getParameter<bool>("useCaloMuons");
     muonMinValidStaHits_ = iConfig.getParameter<int>("muonMinValidStaHits");
 
     response_function = 0;
     tcmetAlgo_=new TCMETAlgo();
   }

MuonTCMETValueMapProducer::~MuonTCMETValueMapProducer ( )

Definition at line 113 of file MuonTCMETValueMapProducer.cc.

References tcmetAlgo_.

   {
  
     // do anything here that needs to be done at desctruction time
     // (e.g. close files, deallocate resources etc.)
     delete tcmetAlgo_;
   }

Member Function Documentation

void MuonTCMETValueMapProducer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 220 of file MuonTCMETValueMapProducer.cc.

References TCMETAlgo::getResponseFunction_fit(), TCMETAlgo::getResponseFunction_mode(), response_function, rfType_, and tcmetAlgo_.

   {
 
     if( rfType_ == 1 )
          response_function = tcmetAlgo_->getResponseFunction_fit();
     else if( rfType_ == 2 )
          response_function = tcmetAlgo_->getResponseFunction_mode();
   }

void MuonTCMETValueMapProducer::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 230 of file MuonTCMETValueMapProducer.cc.

230 {

231 }

bool MuonTCMETValueMapProducer::isGoodCaloMuon	(	const reco::Muon *	muon,
		const unsigned int	index
	)

private

Definition at line 271 of file MuonTCMETValueMapProducer.cc.

References funct::cos(), HLTFastRecoForTau_cff::deltaEta, SiPixelRawToDigiRegional_cfi::deltaPhi, deltaR(), edm::Ref< C, T, F >::get(), reco::Muon::innerTrack(), isGoodTrack(), edm::Ref< C, T, F >::isNonnull(), RPCpg::mu, muon_h, muonDeltaR_, and reco::LeafCandidate::pt().

Referenced by produce().

                                                                                                  {
 
     if( muon->pt() < 10 ) return false;
 
     if( !isGoodTrack( muon ) ) return false;
 
     const reco::TrackRef inputSiliconTrack = muon->innerTrack();
     if( !inputSiliconTrack.isNonnull() ) return false;
 
     //check if it is in the vicinity of a global or tracker muon
     unsigned int nMuons = muon_h->size();
     for (unsigned int iMu = 0; iMu < nMuons; iMu++) {
 
       if( iMu == index ) continue;
 
       const reco::Muon* mu = &(*muon_h)[iMu];
 
       const reco::TrackRef testSiliconTrack = mu->innerTrack();
       if( !testSiliconTrack.isNonnull() ) continue;
 
       double deltaEta = inputSiliconTrack.get()->eta() - testSiliconTrack.get()->eta();
       double deltaPhi = acos( cos( inputSiliconTrack.get()->phi() - testSiliconTrack.get()->phi() ) );
       double deltaR   = TMath::Sqrt( deltaEta * deltaEta + deltaPhi * deltaPhi );
 
       if( deltaR < muonDeltaR_ ) return false;
     }
 
     return true;
   }

bool MuonTCMETValueMapProducer::isGoodMuon ( const reco::Muon * muon )

private

Definition at line 234 of file MuonTCMETValueMapProducer.cc.

References beamSpot_h, reco::LeafCandidate::eta(), reco::Muon::globalTrack(), reco::Muon::innerTrack(), reco::Muon::isGlobalMuon(), edm::Ref< C, T, F >::isNonnull(), reco::Muon::isTrackerMuon(), edm::HandleBase::isValid(), muon_dptrel_, muonchi2_, muond0_, muoneta_, muonGlobal_, muonhits_, muonMinValidStaHits_, muonpt_, muonTracker_, and reco::LeafCandidate::pt().

Referenced by produce().

                                                                    {
     double d0    = -999;
     double nhits = 0;
     double chi2  = 999;  
 
     // get d0 corrected for beam spot
     bool haveBeamSpot = true;
     if( !beamSpot_h.isValid() ) haveBeamSpot = false;
 
     if( muonGlobal_  && !muon->isGlobalMuon()  ) return false;
     if( muonTracker_ && !muon->isTrackerMuon() ) return false;
 
     const reco::TrackRef siTrack     = muon->innerTrack();
     const reco::TrackRef globalTrack = muon->globalTrack();
 
     Point bspot = haveBeamSpot ? beamSpot_h->position() : Point(0,0,0);
     if( siTrack.isNonnull() ) nhits = siTrack->numberOfValidHits();
     if( globalTrack.isNonnull() ) {
       d0   = -1 * globalTrack->dxy( bspot );
       chi2 = globalTrack->normalizedChi2();
     }
 
     if( fabs( d0 ) > muond0_ )                          return false;
     if( muon->pt() < muonpt_ )                          return false;
     if( fabs( muon->eta() ) > muoneta_ )                return false;
     if( nhits < muonhits_ )                             return false;
     if( chi2 > muonchi2_ )                              return false;
     if( globalTrack->hitPattern().numberOfValidMuonHits() < muonMinValidStaHits_ ) return false;
 
     //reject muons with tracker dpt/pt > X
     if( !siTrack.isNonnull() )                                return false;
     if( siTrack->ptError() / siTrack->pt() > muon_dptrel_ )   return false;
 
     else return true;
   }

bool MuonTCMETValueMapProducer::isGoodTrack ( const reco::Muon * muon )

private

Definition at line 302 of file MuonTCMETValueMapProducer.cc.

References beamSpot_h, GOODCOLL_filter_cfg::cut, d0cuta_, d0cutb_, hasValidVertex, i, reco::Muon::innerTrack(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), maxchi2_, maxchi2_tight_, maxd0cut_, maxeta_, maxpt_, maxpt_eta20_, maxpt_eta25_, maxPtErr_, maxPtErr_tight_, maxTrackAlgo_, minhits_, minhits_tight_, nLayers(), nLayers_, nLayersTight_, funct::pow(), mathSSE::sqrt(), trkAlgos_, trkQuality_, vertexColl, and vertexMaxDZ_.

Referenced by isGoodCaloMuon(), and produce().

                                                                     {
     double d0    = -999;
 
     const reco::TrackRef siTrack = muon->innerTrack();
     if (!siTrack.isNonnull())
       return false;
 
     if( hasValidVertex ){
       //get d0 corrected for primary vertex
             
       const Point pvtx = Point(vertexColl->begin()->x(),
                                vertexColl->begin()->y(), 
                                vertexColl->begin()->z());
             
       d0 = -1 * siTrack->dxy( pvtx );
             
       double dz = siTrack->dz( pvtx );
             
       if( fabs( dz ) < vertexMaxDZ_ ){
               
         //get d0 corrected for pvtx
         d0 = -1 * siTrack->dxy( pvtx );
               
       }else{
               
         // get d0 corrected for beam spot
         bool haveBeamSpot = true;
         if( !beamSpot_h.isValid() ) haveBeamSpot = false;
               
         Point bspot = haveBeamSpot ? beamSpot_h->position() : Point(0,0,0);
         d0 = -1 * siTrack->dxy( bspot );
               
       }
     }else{
        
       // get d0 corrected for beam spot
       bool haveBeamSpot = true;
       if( !beamSpot_h.isValid() ) haveBeamSpot = false;
        
       Point bspot = haveBeamSpot ? beamSpot_h->position() : Point(0,0,0);
       d0 = -1 * siTrack->dxy( bspot );
     }
 
     if( siTrack->algo() < maxTrackAlgo_ ){
       //1st 4 tracking iterations (pT-dependent d0 cut)
        
       float d0cut = sqrt(std::pow(d0cuta_,2) + std::pow(d0cutb_/siTrack->pt(),2)); 
       if(d0cut > maxd0cut_) d0cut = maxd0cut_;
        
       if( fabs( d0 ) > d0cut )            return false;    
       if( nLayers( siTrack ) < nLayers_ ) return false;
     }
     else{
       //last 2 tracking iterations (tighten chi2, nhits, pt error cuts)
      
       if( siTrack->normalizedChi2() > maxchi2_tight_ )               return false;
       if( siTrack->numberOfValidHits() < minhits_tight_ )            return false;
       if( (siTrack->ptError() / siTrack->pt()) > maxPtErr_tight_ )   return false;
       if( nLayers( siTrack ) < nLayersTight_ )                       return false;
     }
 
     if( siTrack->numberOfValidHits() < minhits_ )                         return false;
     if( siTrack->normalizedChi2() > maxchi2_ )                            return false;
     if( fabs( siTrack->eta() ) > maxeta_ )                                return false;
     if( siTrack->pt() > maxpt_ )                                          return false;
     if( (siTrack->ptError() / siTrack->pt()) > maxPtErr_ )                return false;
     if( fabs( siTrack->eta() ) > 2.5 && siTrack->pt() > maxpt_eta25_ )    return false;
     if( fabs( siTrack->eta() ) > 2.0 && siTrack->pt() > maxpt_eta20_ )    return false;
 
     int cut = 0;      
     for( unsigned int i = 0; i < trkQuality_.size(); i++ ) {
 
       cut |= (1 << trkQuality_.at(i));
     }
 
     if( !( (siTrack->qualityMask() & cut) == cut ) ) return false;
       
     bool isGoodAlgo = false;    
     if( trkAlgos_.size() == 0 ) isGoodAlgo = true;
     for( unsigned int i = 0; i < trkAlgos_.size(); i++ ) {
 
       if( siTrack->algo() == trkAlgos_.at(i) ) isGoodAlgo = true;
     }
 
     if( !isGoodAlgo ) return false;
       
     return true;
   }

bool MuonTCMETValueMapProducer::isValidVertex ( )

private

Definition at line 453 of file MuonTCMETValueMapProducer.cc.

References funct::pow(), mathSSE::sqrt(), vertexColl, vertexNdof_, vertexRho_, and vertexZ_.

Referenced by produce().

                                                {
     
     if( vertexColl->begin()->isFake()                ) return false;
     if( vertexColl->begin()->ndof() < vertexNdof_    ) return false;
     if( fabs( vertexColl->begin()->z() ) > vertexZ_  ) return false;
     if( sqrt( std::pow( vertexColl->begin()->x() , 2 ) + std::pow( vertexColl->begin()->y() , 2 ) ) > vertexRho_ ) return false;
     
     return true;
     
   }

int MuonTCMETValueMapProducer::nLayers ( const reco::TrackRef track )

private

Definition at line 446 of file MuonTCMETValueMapProducer.cc.

References AlCaHLTBitMon_ParallelJobs::p, and reco::HitPattern::trackerLayersWithMeasurement().

Referenced by isGoodTrack().

                                                                 {
     const reco::HitPattern& p = track->hitPattern();
     return p.trackerLayersWithMeasurement();
   }

void MuonTCMETValueMapProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 126 of file MuonTCMETValueMapProducer.cc.

References beamSpot_h, beamSpotInputTag_, bField, funct::cos(), edm::helper::Filler< Map >::fill(), edm::EventSetup::get(), edm::Event::getByLabel(), hasValidVertex, reco::Muon::innerTrack(), edm::helper::Filler< Map >::insert(), reco::Muon::isCaloMuon(), reco::Muon::isGlobalMuon(), isGoodCaloMuon(), isGoodMuon(), isGoodTrack(), reco::Muon::isTrackerMuon(), edm::HandleBase::isValid(), isValidVertex(), minpt_, RPCpg::mu, muon_h, reco::MuonMETCorrectionData::MuonCandidateValuesUsed, muonInputTag_, reco::MuonMETCorrectionData::NotUsed, edm::Handle< T >::product(), edm::ESHandle< class >::product(), propagateTrack(), edm::Event::put(), response_function, funct::sin(), reco::MuonMETCorrectionData::TreatedAsPion, useCaloMuons_, usePvtxd0_, vertexColl, VertexHandle, and vertexInputTag_.

                                                                                        {
   
     //get input collections
     iEvent.getByLabel(muonInputTag_    , muon_h    );
     iEvent.getByLabel(beamSpotInputTag_, beamSpot_h);
 
     //get vertex collection
     hasValidVertex = false;
     if( usePvtxd0_ ){
       iEvent.getByLabel( vertexInputTag_  , VertexHandle  );
        
       if( VertexHandle.isValid() ) {
         vertexColl = VertexHandle.product();
         hasValidVertex = isValidVertex();
       } 
     }
 
     //get the Bfield
     edm::ESHandle<MagneticField> theMagField;
     iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
     bField = theMagField.product();
 
     //make a ValueMap of ints => flags for 
     //met correction. The values and meanings of the flags are :
     // flag==0 --->    The muon is not used to correct the MET by default
     // flag==1 --->    The muon is used to correct the MET. The Global pt is used.
     // flag==2 --->    The muon is used to correct the MET. The tracker pt is used.
     // flag==3 --->    The muon is used to correct the MET. The standalone pt is used.
     // flag==4 --->    The muon is used to correct the MET as pion using the tcMET ZSP RF.
     // flag==5 --->    The muon is used to correct the MET.  The default fit is used; i.e. we get the pt from muon->pt().
     // In general, the flag should never be 3. You do not want to correct the MET using
     // the pt measurement from the standalone system (unless you really know what you're 
     // doing
 
     std::auto_ptr<edm::ValueMap<reco::MuonMETCorrectionData> > vm_muCorrData(new edm::ValueMap<reco::MuonMETCorrectionData>());
 
     std::vector<reco::MuonMETCorrectionData> v_muCorrData;
 
     unsigned int nMuons = muon_h->size();
 
     for (unsigned int iMu = 0; iMu < nMuons; iMu++) {
 
       const reco::Muon* mu = &(*muon_h)[iMu];
       double deltax = 0.0;
       double deltay = 0.0;
 
       reco::MuonMETCorrectionData muMETCorrData(reco::MuonMETCorrectionData::NotUsed, deltax, deltay);
     
       reco::TrackRef mu_track;
       if( mu->isGlobalMuon() || mu->isTrackerMuon() || mu->isCaloMuon() )
         mu_track = mu->innerTrack();
       else {
         v_muCorrData.push_back( muMETCorrData );
         continue;
       }
 
       // figure out depositions muons would make if they were treated as pions
       if( isGoodTrack( mu ) ) {
 
         if( mu_track->pt() < minpt_ ) 
           muMETCorrData = reco::MuonMETCorrectionData(reco::MuonMETCorrectionData::TreatedAsPion, deltax, deltay);
 
         else {
           int bin_index   = response_function->FindBin( mu_track->eta(), mu_track->pt() );
           double response = response_function->GetBinContent( bin_index );
 
           TVector3 outerTrkPosition = propagateTrack( mu );
 
           deltax = response * mu_track->p() * sin( outerTrkPosition.Theta() ) * cos( outerTrkPosition.Phi() );
           deltay = response * mu_track->p() * sin( outerTrkPosition.Theta() ) * sin( outerTrkPosition.Phi() );
 
           muMETCorrData = reco::MuonMETCorrectionData(reco::MuonMETCorrectionData::TreatedAsPion, deltax, deltay);
         }
       }
 
       // figure out muon flag
       if( isGoodMuon( mu ) )
         v_muCorrData.push_back( reco::MuonMETCorrectionData(reco::MuonMETCorrectionData::MuonCandidateValuesUsed, deltax, deltay) );
 
       else if( useCaloMuons_ && isGoodCaloMuon( mu, iMu ) )
         v_muCorrData.push_back( reco::MuonMETCorrectionData(reco::MuonMETCorrectionData::MuonCandidateValuesUsed, deltax, deltay) );
 
       else v_muCorrData.push_back( muMETCorrData );
     }
     
     edm::ValueMap<reco::MuonMETCorrectionData>::Filler dataFiller(*vm_muCorrData);
 
     dataFiller.insert( muon_h, v_muCorrData.begin(), v_muCorrData.end());
     dataFiller.fill();
     
     iEvent.put(vm_muCorrData, "muCorrData");    
   }

TVector3 MuonTCMETValueMapProducer::propagateTrack ( const reco::Muon * muon )

private

Definition at line 392 of file MuonTCMETValueMapProducer.cc.

References alongMomentum, AnalyticalPropagator_cfi::AnalyticalPropagator, Reference_intrackfit_cff::barrel, bField, Plane::build(), newFWLiteAna::build, edm::Ref< C, T, F >::get(), TrajectoryStateOnSurface::globalPosition(), reco::Muon::innerTrack(), edm::Ref< C, T, F >::isNonnull(), TrajectoryStateOnSurface::isValid(), M_PI, Pi, CosmicsPD_Skims::radius, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by produce().

                                                                           {
 
     TVector3 outerTrkPosition;
 
     outerTrkPosition.SetPtEtaPhi( 999., -10., 2 * TMath::Pi() );
 
     const reco::TrackRef track = muon->innerTrack();
 
     if( !track.isNonnull() ) {
       return outerTrkPosition;
     }
 
     GlobalPoint  tpVertex ( track->vx(), track->vy(), track->vz() );
     GlobalVector tpMomentum ( track.get()->px(), track.get()->py(), track.get()->pz() );
     int tpCharge ( track->charge() );
 
     FreeTrajectoryState fts ( tpVertex, tpMomentum, tpCharge, bField);
 
     const double zdist = 314.;
 
     const double radius = 130.;
 
     const double corner = 1.479;
 
     Plane::PlanePointer lendcap = Plane::build( Plane::PositionType (0, 0, -zdist), Plane::RotationType () );
     Plane::PlanePointer rendcap = Plane::build( Plane::PositionType (0, 0, zdist),  Plane::RotationType () );
 
     Cylinder::CylinderPointer barrel = Cylinder::build( Cylinder::PositionType (0, 0, 0), Cylinder::RotationType (), radius);
 
     AnalyticalPropagator myAP (bField, alongMomentum, 2*M_PI);
 
     TrajectoryStateOnSurface tsos;
 
     if( track.get()->eta() < -corner ) {
       tsos = myAP.propagate( fts, *lendcap);
     }
     else if( fabs(track.get()->eta()) < corner ) {
       tsos = myAP.propagate( fts, *barrel);
     }
     else if( track.get()->eta() > corner ) {
       tsos = myAP.propagate( fts, *rendcap);
     }
 
     if( tsos.isValid() )
       outerTrkPosition.SetXYZ( tsos.globalPosition().x(), tsos.globalPosition().y(), tsos.globalPosition().z() );
 
     else 
       outerTrkPosition.SetPtEtaPhi( 999., -10., 2 * TMath::Pi() );
 
     return outerTrkPosition;
   }

Member Data Documentation

edm::Handle<reco::BeamSpot> cms::MuonTCMETValueMapProducer::beamSpot_h

private

Definition at line 53 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), isGoodTrack(), and produce().

edm::InputTag cms::MuonTCMETValueMapProducer::beamSpotInputTag_

private

Definition at line 57 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

class MagneticField* cms::MuonTCMETValueMapProducer::bField

private

Definition at line 60 of file MuonTCMETValueMapProducer.h.

Referenced by produce(), and propagateTrack().

double cms::MuonTCMETValueMapProducer::d0cuta_

private

Definition at line 91 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::d0cutb_

private

Definition at line 92 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

bool cms::MuonTCMETValueMapProducer::hasValidVertex

private

Definition at line 69 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and produce().

double cms::MuonTCMETValueMapProducer::maxchi2_

private

Definition at line 84 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::maxchi2_tight_

private

Definition at line 88 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::maxd0cut_

private

Definition at line 87 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::maxeta_

private

Definition at line 83 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::maxpt_

private

Definition at line 82 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::maxpt_eta20_

private

Definition at line 79 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::maxpt_eta25_

private

Definition at line 78 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::maxPtErr_

private

Definition at line 86 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::maxPtErr_tight_

private

Definition at line 90 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

int cms::MuonTCMETValueMapProducer::maxTrackAlgo_

private

Definition at line 80 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::minhits_

private

Definition at line 85 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::minhits_tight_

private

Definition at line 89 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::minpt_

private

Definition at line 81 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

double cms::MuonTCMETValueMapProducer::muon_dptrel_

private

Definition at line 104 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

edm::Handle<reco::MuonCollection> cms::MuonTCMETValueMapProducer::muon_h

private

Definition at line 52 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodCaloMuon(), and produce().

double cms::MuonTCMETValueMapProducer::muonchi2_

private

Definition at line 100 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::muond0_

private

Definition at line 102 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::muonDeltaR_

private

Definition at line 103 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodCaloMuon(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::muoneta_

private

Definition at line 99 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

bool cms::MuonTCMETValueMapProducer::muonGlobal_

private

Definition at line 66 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::muonhits_

private

Definition at line 101 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

edm::InputTag cms::MuonTCMETValueMapProducer::muonInputTag_

private

Definition at line 56 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

int cms::MuonTCMETValueMapProducer::muonMinValidStaHits_

private

Definition at line 97 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::muonpt_

private

Definition at line 98 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

bool cms::MuonTCMETValueMapProducer::muonTracker_

private

Definition at line 67 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodMuon(), and MuonTCMETValueMapProducer().

int cms::MuonTCMETValueMapProducer::nLayers_

private

Definition at line 72 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

int cms::MuonTCMETValueMapProducer::nLayersTight_

private

Definition at line 73 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

class TH2D* cms::MuonTCMETValueMapProducer::response_function

private

Definition at line 64 of file MuonTCMETValueMapProducer.h.

Referenced by beginJob(), MuonTCMETValueMapProducer(), and produce().

int cms::MuonTCMETValueMapProducer::rfType_

private

Definition at line 71 of file MuonTCMETValueMapProducer.h.

Referenced by beginJob(), and MuonTCMETValueMapProducer().

TCMETAlgo* cms::MuonTCMETValueMapProducer::tcmetAlgo_

private

Definition at line 105 of file MuonTCMETValueMapProducer.h.

Referenced by beginJob(), MuonTCMETValueMapProducer(), and ~MuonTCMETValueMapProducer().

std::vector<int> cms::MuonTCMETValueMapProducer::trkAlgos_

private

Definition at line 95 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

std::vector<int> cms::MuonTCMETValueMapProducer::trkQuality_

private

Definition at line 94 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

bool cms::MuonTCMETValueMapProducer::useCaloMuons_

private

Definition at line 68 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

bool cms::MuonTCMETValueMapProducer::usePvtxd0_

private

Definition at line 93 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

const reco::VertexCollection* cms::MuonTCMETValueMapProducer::vertexColl

private

Definition at line 62 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), isValidVertex(), and produce().

edm::Handle<reco::VertexCollection> cms::MuonTCMETValueMapProducer::VertexHandle

private

Definition at line 54 of file MuonTCMETValueMapProducer.h.

Referenced by produce().

edm::InputTag cms::MuonTCMETValueMapProducer::vertexInputTag_

private

Definition at line 58 of file MuonTCMETValueMapProducer.h.

Referenced by MuonTCMETValueMapProducer(), and produce().

double cms::MuonTCMETValueMapProducer::vertexMaxDZ_

private

Definition at line 77 of file MuonTCMETValueMapProducer.h.

Referenced by isGoodTrack(), and MuonTCMETValueMapProducer().

int cms::MuonTCMETValueMapProducer::vertexNdof_

private

Definition at line 74 of file MuonTCMETValueMapProducer.h.

Referenced by isValidVertex(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::vertexRho_

private

Definition at line 76 of file MuonTCMETValueMapProducer.h.

Referenced by isValidVertex(), and MuonTCMETValueMapProducer().

double cms::MuonTCMETValueMapProducer::vertexZ_

private

Definition at line 75 of file MuonTCMETValueMapProducer.h.

Referenced by isValidVertex(), and MuonTCMETValueMapProducer().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation