CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
JetIDProducer.cc
Go to the documentation of this file.
3 
4 #include <vector>
5 
6 //
7 // constants, enums and typedefs
8 //
9 
10 
11 //
12 // static data member definitions
13 //
14 
15 //
16 // constructors and destructor
17 //
19  src_ ( iConfig.getParameter<edm::InputTag>("src") ),
20  helper_ ( iConfig ),
21  muHelper_ ( iConfig )
22 {
23  produces< reco::JetIDValueMap >();
24 }
25 
26 
28 {
29 }
30 
31 
32 //
33 // member functions
34 //
35 
36 // ------------ method called to produce the data ------------
37 void
39 {
40 
41  // get the input jets
43  iEvent.getByLabel( src_, h_jets );
44 
45  // allocate the jet--->jetid value map
46  std::auto_ptr<reco::JetIDValueMap> jetIdValueMap( new reco::JetIDValueMap );
47  // instantiate the filler with the map
48  reco::JetIDValueMap::Filler filler(*jetIdValueMap);
49 
50  // allocate the vector of ids
51  size_t njets = h_jets->size();
52  std::vector<reco::JetID> ids (njets);
53 
54  // loop over the jets
55  for ( edm::View<reco::CaloJet>::const_iterator jetsBegin = h_jets->begin(),
56  jetsEnd = h_jets->end(),
57  ijet = jetsBegin;
58  ijet != jetsEnd; ++ijet ) {
59 
60  // get the id from each jet
61  helper_.calculate( iEvent, *ijet );
62 
63  muHelper_.calculate( iEvent, iSetup, *ijet );
64 
65  ids[ijet-jetsBegin].fHPD = helper_.fHPD();
66  ids[ijet-jetsBegin].fRBX = helper_.fRBX();
67  ids[ijet-jetsBegin].n90Hits = helper_.n90Hits();
68  ids[ijet-jetsBegin].fSubDetector1 = helper_.fSubDetector1();
69  ids[ijet-jetsBegin].fSubDetector2 = helper_.fSubDetector2();
70  ids[ijet-jetsBegin].fSubDetector3 = helper_.fSubDetector3();
71  ids[ijet-jetsBegin].fSubDetector4 = helper_.fSubDetector4();
72  ids[ijet-jetsBegin].restrictedEMF = helper_.restrictedEMF();
73  ids[ijet-jetsBegin].nHCALTowers = helper_.nHCALTowers();
74  ids[ijet-jetsBegin].nECALTowers = helper_.nECALTowers();
75  ids[ijet-jetsBegin].approximatefHPD = helper_.approximatefHPD();
76  ids[ijet-jetsBegin].approximatefRBX = helper_.approximatefRBX();
77  ids[ijet-jetsBegin].hitsInN90 = helper_.hitsInN90();
78 
79  ids[ijet-jetsBegin].numberOfHits2RPC = muHelper_.numberOfHits2RPC();
80  ids[ijet-jetsBegin].numberOfHits3RPC = muHelper_.numberOfHits3RPC();
81  ids[ijet-jetsBegin].numberOfHitsRPC = muHelper_.numberOfHitsRPC();
82 
83  ids[ijet-jetsBegin].fEB = helper_.fEB ();
84  ids[ijet-jetsBegin].fEE = helper_.fEE ();
85  ids[ijet-jetsBegin].fHB = helper_.fHB ();
86  ids[ijet-jetsBegin].fHE = helper_.fHE ();
87  ids[ijet-jetsBegin].fHO = helper_.fHO ();
88  ids[ijet-jetsBegin].fLong = helper_.fLong ();
89  ids[ijet-jetsBegin].fShort = helper_.fShort();
90  ids[ijet-jetsBegin].fLS = helper_.fLSbad ();
91  ids[ijet-jetsBegin].fHFOOT = helper_.fHFOOT();
92  }
93 
94  // set up the map
95  filler.insert( h_jets, ids.begin(), ids.end() );
96 
97  // fill the vals
98  filler.fill();
99 
100  // write map to the event
101  iEvent.put( jetIdValueMap );
102 }
103 
104 // ------------ method called once each job just before starting event loop ------------
105 void
107 {
108 }
109 
110 // ------------ method called once each job just after ending the event loop ------------
111 void
113 }
114 
115 //define this as a plug-in
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:81
double fHE() const
Definition: JetIDHelper.h:43
double fHB() const
Definition: JetIDHelper.h:42
double fLSbad() const
Definition: JetIDHelper.h:47
int nHCALTowers() const
Definition: JetIDHelper.h:51
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
double approximatefRBX() const
Definition: JetIDHelper.h:55
double fLong() const
Definition: JetIDHelper.h:45
void insert(const H &h, I begin, I end)
Definition: ValueMap.h:53
int nECALTowers() const
Definition: JetIDHelper.h:52
double fHPD() const
Definition: JetIDHelper.h:33
double fSubDetector4() const
Definition: JetIDHelper.h:39
reco::helper::JetIDHelper helper_
Definition: JetIDProducer.h:59
double restrictedEMF() const
Definition: JetIDHelper.h:50
double fSubDetector1() const
Definition: JetIDHelper.h:36
void calculate(const edm::Event &event, const edm::EventSetup &isetup, const reco::Jet &jet, const int iDbg=0)
edm::InputTag src_
Definition: JetIDProducer.h:58
int iEvent
Definition: GenABIO.cc:243
double approximatefHPD() const
Definition: JetIDHelper.h:54
double fEB() const
Definition: JetIDHelper.h:40
JetIDProducer(const edm::ParameterSet &)
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:94
double fEE() const
Definition: JetIDHelper.h:41
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
virtual void endJob()
double fSubDetector3() const
Definition: JetIDHelper.h:38
double fShort() const
Definition: JetIDHelper.h:46
reco::helper::JetMuonHitsIDHelper muHelper_
Definition: JetIDProducer.h:60
double fSubDetector2() const
Definition: JetIDHelper.h:37
double fHO() const
Definition: JetIDHelper.h:44
double fHFOOT() const
Definition: JetIDHelper.h:48
virtual void produce(edm::Event &, const edm::EventSetup &)
virtual void beginJob()
void calculate(const edm::Event &event, const reco::CaloJet &jet, const int iDbg=0)
Definition: JetIDHelper.cc:86
double fRBX() const
Definition: JetIDHelper.h:34