|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
|
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/RecoCandidate/interface/CaloRecHitCandidate.h"
namespace reco {
namespace modules {
template <typename HitCollection>
class CaloRecHitCandidateProducer : public edm::global::EDProducer<> {
public:
/// constructor
CaloRecHitCandidateProducer(const edm::ParameterSet &cfg)
: srcToken_(consumes<HitCollection>(cfg.template getParameter<edm::InputTag>("src"))) {}
/// destructor
~CaloRecHitCandidateProducer() override {}
private:
/// process one event
void produce(edm::StreamID, edm::Event &, const edm::EventSetup &) const override;
/// source collection tag
const edm::EDGetTokenT<HitCollection> srcToken_;
};
} // namespace modules
} // namespace reco
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Math/interface/Vector3D.h"
namespace reco {
namespace modules {
template <typename HitCollection>
void CaloRecHitCandidateProducer<HitCollection>::produce(edm::StreamID,
edm::Event &evt,
const edm::EventSetup &) const {
using namespace edm;
using namespace reco;
using namespace std;
Handle<HitCollection> hits;
evt.getByToken(srcToken_, hits);
unique_ptr<CandidateCollection> cands(new CandidateCollection);
size_t size = hits->size();
cands->reserve(size);
for (size_t idx = 0; idx != size; ++idx) {
const CaloRecHit &hit = (*hits)[idx];
/// don't know how to set eta and phi
double eta = 0, phi = 0, energy = hit.energy();
math::RhoEtaPhiVector p(1, eta, phi);
p *= (energy / p.r());
CaloRecHitCandidate *c = new CaloRecHitCandidate(Candidate::LorentzVector(p.x(), p.y(), p.z(), energy));
c->setCaloRecHit(Ptr<CaloRecHit>(hits, idx));
cands->push_back(c);
}
evt.put(std::move(cands));
}
} // namespace modules
} // namespace reco
#include "FWCore/Framework/interface/MakerMacros.h"
#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
typedef reco::modules::CaloRecHitCandidateProducer<HFRecHitCollection> HFRecHitCandidateProducer;
typedef reco::modules::CaloRecHitCandidateProducer<ZDCRecHitCollection> ZDCRecHitCandidateProducer;
typedef reco::modules::CaloRecHitCandidateProducer<HBHERecHitCollection> HBHERecHitCandidateProducer;
typedef reco::modules::CaloRecHitCandidateProducer<HORecHitCollection> HORecHitCandidateProducer;
DEFINE_FWK_MODULE(HFRecHitCandidateProducer);
DEFINE_FWK_MODULE(ZDCRecHitCandidateProducer);
DEFINE_FWK_MODULE(HBHERecHitCandidateProducer);
DEFINE_FWK_MODULE(HORecHitCandidateProducer);
|