L1Analysis::L1AnalysisRecoMuon2 Class Reference

#include <L1AnalysisRecoMuon2.h>

Public Member Functions
L1AnalysisRecoMuon2DataFormat *	getData ()
void	init (const edm::EventSetup &eventSetup)
	L1AnalysisRecoMuon2 (const edm::ParameterSet &pset, edm::ConsumesCollector)
void	Reset ()
void	SetMuon (const edm::Event &event, const edm::EventSetup &setup, const edm::Handle< reco::MuonCollection > muons, const edm::Handle< reco::VertexCollection > vertices, double METx, double METy, unsigned maxMuon)
	~L1AnalysisRecoMuon2 ()

Private Attributes
PropagateToMuon	muPropagator1st_
PropagateToMuon	muPropagator2nd_
L1AnalysisRecoMuon2DataFormat	recoMuon_

Detailed Description

Definition at line 41 of file L1AnalysisRecoMuon2.h.

Constructor & Destructor Documentation

L1Analysis::L1AnalysisRecoMuon2::L1AnalysisRecoMuon2	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector	iC
	)

Definition at line 10 of file L1AnalysisRecoMuon2.cc.

    : muPropagator1st_(pset.getParameter<edm::ParameterSet>("muProp1st"), iC),
      muPropagator2nd_(pset.getParameter<edm::ParameterSet>("muProp2nd"), iC) {}

L1Analysis::L1AnalysisRecoMuon2::~L1AnalysisRecoMuon2

(

)

Definition at line 14 of file L1AnalysisRecoMuon2.cc.

{}

Member Function Documentation

L1AnalysisRecoMuon2DataFormat* L1Analysis::L1AnalysisRecoMuon2::getData ( void  )

inline

Definition at line 60 of file L1AnalysisRecoMuon2.h.

References recoMuon_.

{ return &recoMuon_; }

void L1Analysis::L1AnalysisRecoMuon2::init

(

const edm::EventSetup &

eventSetup

)

Definition at line 94 of file L1AnalysisRecoMuon2.cc.

                                                                        {
  muPropagator1st_.init(eventSetup);
  muPropagator2nd_.init(eventSetup);
}

void L1Analysis::L1AnalysisRecoMuon2::Reset ( )

inline

Definition at line 61 of file L1AnalysisRecoMuon2.h.

References recoMuon_, and L1Analysis::L1AnalysisRecoMuon2DataFormat::Reset().

{ recoMuon_.Reset(); }

void L1Analysis::L1AnalysisRecoMuon2::SetMuon	(	const edm::Event &	event,
		const edm::EventSetup &	setup,
		const edm::Handle< reco::MuonCollection >	muons,
		const edm::Handle< reco::VertexCollection >	vertices,
		double	METx,
		double	METy,
		unsigned	maxMuon
	)

Definition at line 16 of file L1AnalysisRecoMuon2.cc.

References PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalPosition(), isLooseMuonCustom(), isMediumMuonCustom(), isTightMuonCustom(), TrajectoryStateOnSurface::isValid(), edm::HandleBase::isValid(), SiStripPI::max, dumpRecoGeometry_cfg::Muon, and PV3DBase< T, PVType, FrameType >::phi().

                                                                {
  recoMuon_.nMuons = 0;

  for (reco::MuonCollection::const_iterator it = muons->begin(); it != muons->end() && recoMuon_.nMuons < maxMuon;
       ++it) {
    recoMuon_.e.push_back(it->energy());
    recoMuon_.pt.push_back(it->pt());
    recoMuon_.et.push_back(it->et());
    recoMuon_.eta.push_back(it->eta());
    recoMuon_.phi.push_back(it->phi());
    recoMuon_.charge.push_back(it->charge());

    //check isLooseMuon
    bool flagLoose = isLooseMuonCustom(*it);
    recoMuon_.isLooseMuon.push_back(flagLoose);

    //check isMediumMuon
    bool flagMedium = isMediumMuonCustom(*it);
    recoMuon_.isMediumMuon.push_back(flagMedium);

    //check isTightMuon
    bool flagTight = false;
    if (vertices.isValid())
      flagTight = isTightMuonCustom(*it, (*vertices)[0]);
    recoMuon_.isTightMuon.push_back(flagTight);

    double iso = (it->pfIsolationR04().sumChargedHadronPt +
                  max(0.,
                      it->pfIsolationR04().sumNeutralHadronEt + it->pfIsolationR04().sumPhotonEt -
                          0.5 * it->pfIsolationR04().sumPUPt)) /
                 it->pt();
    recoMuon_.iso.push_back(iso);

    double MET_local = TMath::Sqrt(METx * METx + METy * METy);
    recoMuon_.met.push_back(MET_local);

    TLorentzVector METP4;
    METP4.SetPxPyPzE(METx, METy, 0, MET_local);

    TLorentzVector Muon;
    Muon.SetPtEtaPhiE(it->pt(), it->eta(), it->phi(), it->energy());

    double scalSum = MET_local + Muon.Pt();
    TLorentzVector vecSum(Muon);
    vecSum += METP4;
    double vecSumPt = vecSum.Pt();

    recoMuon_.mt.push_back(TMath::Sqrt(scalSum * scalSum - vecSumPt * vecSumPt));

    recoMuon_.nMuons++;

    // extrapolation of track coordinates
    TrajectoryStateOnSurface stateAtMuSt1 = muPropagator1st_.extrapolate(*it);
    if (stateAtMuSt1.isValid()) {
      recoMuon_.etaSt1.push_back(stateAtMuSt1.globalPosition().eta());
      recoMuon_.phiSt1.push_back(stateAtMuSt1.globalPosition().phi());
    } else {
      recoMuon_.etaSt1.push_back(-9999);
      recoMuon_.phiSt1.push_back(-9999);
    }

    TrajectoryStateOnSurface stateAtMuSt2 = muPropagator2nd_.extrapolate(*it);
    if (stateAtMuSt2.isValid()) {
      recoMuon_.etaSt2.push_back(stateAtMuSt2.globalPosition().eta());
      recoMuon_.phiSt2.push_back(stateAtMuSt2.globalPosition().phi());
    } else {
      recoMuon_.etaSt2.push_back(-9999);
      recoMuon_.phiSt2.push_back(-9999);
    }
  }
}

Member Data Documentation

PropagateToMuon L1Analysis::L1AnalysisRecoMuon2::muPropagator1st_

private

Definition at line 66 of file L1AnalysisRecoMuon2.h.

PropagateToMuon L1Analysis::L1AnalysisRecoMuon2::muPropagator2nd_

private

Definition at line 67 of file L1AnalysisRecoMuon2.h.

L1AnalysisRecoMuon2DataFormat L1Analysis::L1AnalysisRecoMuon2::recoMuon_

private

Definition at line 64 of file L1AnalysisRecoMuon2.h.

Referenced by getData(), and Reset().