CMS 3D CMS Logo

JetTracksAssociatorAtCaloFace.cc
Go to the documentation of this file.
1 // \class JetTracksAssociatorAtCaloFace JetTracksAssociatorAtCaloFace.cc
2 //
3 // Original Author: Andrea Rizzi
4 // Created: Wed Apr 12 11:12:49 CEST 2006
5 // Accommodated for Jet Package by: Fedor Ratnikov Jul. 30, 2007
6 //
8 
9 
16 
17 
18 
20 
22  : firstRun(true),
23  dR_(fConfig.getParameter<double>("coneSize"))
24 {
25  mJets = consumes<edm::View <reco::Jet> >(fConfig.getParameter<edm::InputTag> ("jets"));
26  mExtrapolations = consumes<std::vector<reco::TrackExtrapolation> >(fConfig.getParameter<edm::InputTag> ("extrapolations")),
27 
28  produces<reco::JetTracksAssociation::Container> ();
29 }
30 
32 
33 
34  // Get geometry
35  if ( firstRun ) {
36  fSetup.get<CaloGeometryRecord>().get(pGeo);
37  firstRun = false;
38  }
39 
40  if ( !pGeo.isValid() ) {
41  throw cms::Exception("InvalidInput") << "Did not get geometry" << std::endl;
42  }
43 
44  // get stuff from Event
46  fEvent.getByToken (mJets, jets_h);
48  fEvent.getByToken (mExtrapolations, extrapolations_h);
49 
50  auto jetTracks = std::make_unique<reco::JetTracksAssociation::Container>(reco::JetRefBaseProd(jets_h));
51 
52  // Check to make sure we have inputs
53  if ( jets_h->size() == 0 ){
54  // store output regardless the size of the inputs
55  fEvent.put(std::move(jetTracks));
56  return;
57  }
58  // Check to make sure the inputs are calo jets
59  reco::CaloJet const * caloJet0 = dynamic_cast<reco::CaloJet const *>( & (jets_h->at(0)) );
60  // Disallowed non-CaloJet inputs
61  if ( caloJet0 == 0 ) {
62  throw cms::Exception("InvalidInput") << " Jet-track association is only defined for CaloJets.";
63  }
64 
65 
66  // format inputs
67  std::vector <edm::RefToBase<reco::Jet> > allJets;
68  allJets.reserve (jets_h->size());
69  for (unsigned i = 0; i < jets_h->size(); ++i) allJets.push_back (jets_h->refAt(i));
70  mAssociator.produce (&*jetTracks, allJets, *extrapolations_h, *pGeo, dR_ );
71 
72 
73  // store output
74  fEvent.put(std::move(jetTracks));
75 }
76 
77 
std::vector< reco::JetBaseRef > allJets(const Container &)
fill list of all jets associated with values. Return # of jets in the list
T getParameter(std::string const &) const
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
Jets made from CaloTowers.
Definition: CaloJet.h:29
virtual void produce(edm::Event &, const edm::EventSetup &)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
JetTracksAssociationXtrpCalo mAssociator
edm::ESHandle< CaloGeometry > pGeo
edm::RefToBaseProd< reco::Jet > JetRefBaseProd
Definition: JetCollection.h:14
void produce(Association *, JetRefs const &, std::vector< reco::TrackExtrapolation > const &, CaloGeometry const &, double dR)
Associates tracks to jets.
edm::EDGetTokenT< std::vector< reco::TrackExtrapolation > > mExtrapolations
const T & get() const
Definition: EventSetup.h:55
firstRun
Definition: dataset.py:933
bool isValid() const
Definition: ESHandle.h:47
def move(src, dest)
Definition: eostools.py:510
edm::EDGetTokenT< edm::View< reco::Jet > > mJets