#include <JetPlusTrackProducer.cc>

Inheritance diagram for JetPlusTrackProducer:

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

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

	~JetPlusTrackProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
std::string	alias

edm::EDGetTokenT< edm::View< reco::CaloJet > >	input_jets_token_

edm::EDGetTokenT< reco::VertexCollection >	input_vertex_token_

JetPlusTrackCorrector *	mJPTalgo

ZSPJPTJetCorrector *	mZSPalgo

double	ptCUT

edm::InputTag	src

edm::InputTag	srcPVs_

bool	useZSP

bool	vectorial_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 41 of file JetPlusTrackProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 57 of file JetPlusTrackProducer.cc.

References SiStripOfflineCRack_cfg::alias, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), RecoJPTJets_cff::ptCUT, and TrackRefitter_38T_cff::src.

                                                                          {
   //register your products
   src = iConfig.getParameter<edm::InputTag>("src");
   alias = iConfig.getUntrackedParameter<string>("alias");
   srcPVs_ = iConfig.getParameter<edm::InputTag>("srcPVs");
   vectorial_ = iConfig.getParameter<bool>("VectorialCorrection");
   useZSP = iConfig.getParameter<bool>("UseZSP");
   ptCUT = iConfig.getParameter<double>("ptCUT");
   mJPTalgo = new JetPlusTrackCorrector(iConfig, consumesCollector());
   if (useZSP)
     mZSPalgo = new ZSPJPTJetCorrector(iConfig);
 
   produces<reco::JPTJetCollection>().setBranchAlias(alias);
 
   input_jets_token_ = consumes<edm::View<reco::CaloJet> >(src);
   input_vertex_token_ = consumes<reco::VertexCollection>(srcPVs_);
 }

JetPlusTrackProducer::~JetPlusTrackProducer ( )

override

Definition at line 75 of file JetPlusTrackProducer.cc.

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

Member Function Documentation

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

override

Definition at line 85 of file JetPlusTrackProducer.cc.

                                                                                 {
   using namespace edm;
 
   //  std::cout<<" RecoJets::JetPlusTrackProducer::produce "<<std::endl;
 
   // get stuff from Event
   edm::Handle<edm::View<reco::CaloJet> > jets_h;
   iEvent.getByToken(input_jets_token_, jets_h);
 
   //  auto pOut = std::make_unique<reco::CaloJetCollection>();
   auto pOut = std::make_unique<reco::JPTJetCollection>();
 
   for (unsigned i = 0; i < jets_h->size(); ++i) {
     const reco::CaloJet* oldjet = &(*(jets_h->refAt(i)));
 
     reco::CaloJet corrected = *oldjet;
 
     // ZSP corrections
 
     double factorZSP = 1.;
     if (useZSP)
       factorZSP = mZSPalgo->correction(corrected, iEvent, iSetup);
 
     corrected.scaleEnergy(factorZSP);
 
     // JPT corrections
 
     double scaleJPT = 1.;
 
     math::XYZTLorentzVector p4;
 
     jpt::MatchedTracks pions;
     jpt::MatchedTracks muons;
     jpt::MatchedTracks elecs;
     bool ok = false;
 
     if (!vectorial_) {
       scaleJPT = mJPTalgo->correction(corrected, *oldjet, iEvent, iSetup, pions, muons, elecs, ok);
       p4 = math::XYZTLorentzVector(corrected.px() * scaleJPT,
                                    corrected.py() * scaleJPT,
                                    corrected.pz() * scaleJPT,
                                    corrected.energy() * scaleJPT);
     } else {
       scaleJPT = mJPTalgo->correction(corrected, *oldjet, iEvent, iSetup, p4, pions, muons, elecs, ok);
     }
 
     reco::JPTJet::Specific specific;
 
     if (ok) {
       specific.pionsInVertexInCalo = pions.inVertexInCalo_;
       specific.pionsInVertexOutCalo = pions.inVertexOutOfCalo_;
       specific.pionsOutVertexInCalo = pions.outOfVertexInCalo_;
       specific.muonsInVertexInCalo = muons.inVertexInCalo_;
       specific.muonsInVertexOutCalo = muons.inVertexOutOfCalo_;
       specific.muonsOutVertexInCalo = muons.outOfVertexInCalo_;
       specific.elecsInVertexInCalo = elecs.inVertexInCalo_;
       specific.elecsInVertexOutCalo = elecs.inVertexOutOfCalo_;
       specific.elecsOutVertexInCalo = elecs.outOfVertexInCalo_;
     }
 
     // Fill JPT Specific
     edm::RefToBase<reco::Jet> myjet = (edm::RefToBase<reco::Jet>)jets_h->refAt(i);
     specific.theCaloJetRef = myjet;
     specific.mZSPCor = factorZSP;
     specific.mResponseOfChargedWithEff = (float)mJPTalgo->getResponseOfChargedWithEff();
     specific.mResponseOfChargedWithoutEff = (float)mJPTalgo->getResponseOfChargedWithoutEff();
     specific.mSumPtOfChargedWithEff = (float)mJPTalgo->getSumPtWithEff();
     specific.mSumPtOfChargedWithoutEff = (float)mJPTalgo->getSumPtWithoutEff();
     specific.mSumEnergyOfChargedWithEff = (float)mJPTalgo->getSumEnergyWithEff();
     specific.mSumEnergyOfChargedWithoutEff = (float)mJPTalgo->getSumEnergyWithoutEff();
     specific.mChargedHadronEnergy = (float)mJPTalgo->getSumEnergyWithoutEff();
 
     // Fill Charged Jet shape parameters
     double deR2Tr = 0.;
     double deEta2Tr = 0.;
     double dePhi2Tr = 0.;
     double Zch = 0.;
     double Pout2 = 0.;
     double Pout = 0.;
     double denominator_tracks = 0.;
     int ntracks = 0;
 
     for (reco::TrackRefVector::const_iterator it = pions.inVertexInCalo_.begin(); it != pions.inVertexInCalo_.end();
          it++) {
       double deR = deltaR((*it)->eta(), (*it)->phi(), p4.eta(), p4.phi());
       double deEta = (*it)->eta() - p4.eta();
       double dePhi = deltaPhi((*it)->phi(), p4.phi());
       if ((**it).ptError() / (**it).pt() < 0.1) {
         deR2Tr = deR2Tr + deR * deR * (*it)->pt();
         deEta2Tr = deEta2Tr + deEta * deEta * (*it)->pt();
         dePhi2Tr = dePhi2Tr + dePhi * dePhi * (*it)->pt();
         denominator_tracks = denominator_tracks + (*it)->pt();
         Zch = Zch + (*it)->pt();
 
         Pout2 = Pout2 + (**it).p() * (**it).p() - (Zch * p4.P()) * (Zch * p4.P());
         ntracks++;
       }
     }
 
     for (reco::TrackRefVector::const_iterator it = muons.inVertexInCalo_.begin(); it != muons.inVertexInCalo_.end();
          it++) {
       double deR = deltaR((*it)->eta(), (*it)->phi(), p4.eta(), p4.phi());
       double deEta = (*it)->eta() - p4.eta();
       double dePhi = deltaPhi((*it)->phi(), p4.phi());
       if ((**it).ptError() / (**it).pt() < 0.1) {
         deR2Tr = deR2Tr + deR * deR * (*it)->pt();
         deEta2Tr = deEta2Tr + deEta * deEta * (*it)->pt();
         dePhi2Tr = dePhi2Tr + dePhi * dePhi * (*it)->pt();
         denominator_tracks = denominator_tracks + (*it)->pt();
         Zch = Zch + (*it)->pt();
 
         Pout2 = Pout2 + (**it).p() * (**it).p() - (Zch * p4.P()) * (Zch * p4.P());
         ntracks++;
       }
     }
     for (reco::TrackRefVector::const_iterator it = elecs.inVertexInCalo_.begin(); it != elecs.inVertexInCalo_.end();
          it++) {
       double deR = deltaR((*it)->eta(), (*it)->phi(), p4.eta(), p4.phi());
       double deEta = (*it)->eta() - p4.eta();
       double dePhi = deltaPhi((*it)->phi(), p4.phi());
       if ((**it).ptError() / (**it).pt() < 0.1) {
         deR2Tr = deR2Tr + deR * deR * (*it)->pt();
         deEta2Tr = deEta2Tr + deEta * deEta * (*it)->pt();
         dePhi2Tr = dePhi2Tr + dePhi * dePhi * (*it)->pt();
         denominator_tracks = denominator_tracks + (*it)->pt();
         Zch = Zch + (*it)->pt();
 
         Pout2 = Pout2 + (**it).p() * (**it).p() - (Zch * p4.P()) * (Zch * p4.P());
         ntracks++;
       }
     }
     for (reco::TrackRefVector::const_iterator it = pions.inVertexOutOfCalo_.begin();
          it != pions.inVertexOutOfCalo_.end();
          it++) {
       Zch = Zch + (*it)->pt();
     }
     for (reco::TrackRefVector::const_iterator it = muons.inVertexOutOfCalo_.begin();
          it != muons.inVertexOutOfCalo_.end();
          it++) {
       Zch = Zch + (*it)->pt();
     }
     for (reco::TrackRefVector::const_iterator it = elecs.inVertexOutOfCalo_.begin();
          it != elecs.inVertexOutOfCalo_.end();
          it++) {
       Zch = Zch + (*it)->pt();
     }
 
     if (mJPTalgo->getSumPtForBeta() > 0.)
       Zch = Zch / mJPTalgo->getSumPtForBeta();
 
     //     std::cout<<" Zch "<< Zch<<" "<<mJPTalgo->getSumPtForBeta()<<std::endl;
 
     if (ntracks > 0) {
       Pout = sqrt(fabs(Pout2)) / ntracks;
     }
     if (denominator_tracks != 0) {
       deR2Tr = deR2Tr / denominator_tracks;
       deEta2Tr = deEta2Tr / denominator_tracks;
       dePhi2Tr = dePhi2Tr / denominator_tracks;
     }
 
     specific.R2momtr = deR2Tr;
     specific.Eta2momtr = deEta2Tr;
     specific.Phi2momtr = dePhi2Tr;
     specific.Pout = Pout;
     specific.Zch = Zch;
 
     //       std::cout<<" Moments for charged component "<<deR2_Tr<<" "<<deEta2_Tr<<" "<<dePhi2_Tr<<std::endl;
 
     // Create JPT jet
 
     reco::Particle::Point vertex_ = reco::Jet::Point(0, 0, 0);
 
     // If we add primary vertex
     edm::Handle<reco::VertexCollection> pvCollection;
     iEvent.getByToken(input_vertex_token_, pvCollection);
     if (pvCollection.isValid() && !pvCollection->empty())
       vertex_ = pvCollection->begin()->position();
 
     reco::JPTJet fJet(p4, vertex_, specific, corrected.getJetConstituents());
 
     //   fJet.printJet();
 
     // Output module
     if (fJet.pt() > ptCUT)
       pOut->push_back(fJet);
   }
 
   iEvent.put(std::move(pOut));
 }

Member Data Documentation

std::string JetPlusTrackProducer::alias

private

Definition at line 53 of file JetPlusTrackProducer.h.

edm::EDGetTokenT<edm::View<reco::CaloJet> > JetPlusTrackProducer::input_jets_token_

private

Definition at line 58 of file JetPlusTrackProducer.h.

edm::EDGetTokenT<reco::VertexCollection> JetPlusTrackProducer::input_vertex_token_

private

Definition at line 59 of file JetPlusTrackProducer.h.

JetPlusTrackCorrector* JetPlusTrackProducer::mJPTalgo

private

Definition at line 49 of file JetPlusTrackProducer.h.

ZSPJPTJetCorrector* JetPlusTrackProducer::mZSPalgo

private

Definition at line 50 of file JetPlusTrackProducer.h.

double JetPlusTrackProducer::ptCUT

private

Definition at line 56 of file JetPlusTrackProducer.h.

edm::InputTag JetPlusTrackProducer::src

private

Definition at line 51 of file JetPlusTrackProducer.h.

edm::InputTag JetPlusTrackProducer::srcPVs_

private

Definition at line 52 of file JetPlusTrackProducer.h.

bool JetPlusTrackProducer::useZSP

private

Definition at line 55 of file JetPlusTrackProducer.h.

bool JetPlusTrackProducer::vectorial_

private

Definition at line 54 of file JetPlusTrackProducer.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation