CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
JetUnpacker.cc
Go to the documentation of this file.
2 
4 
5 #include "L1TObjectCollections.h"
6 
7 namespace l1t {
8  namespace stage2 {
9  class JetUnpacker : public Unpacker {
10  public:
11  virtual bool unpack(const Block& block, UnpackerCollections *coll) override;
12  };
13  }
14 }
15 
16 // Implementation
17 
18 namespace l1t {
19 namespace stage2 {
20  bool
22  {
23 
24  LogDebug("L1T") << "Block ID = " << block.header().getID() << " size = " << block.header().getSize();
25 
26  int nBX = int(ceil(block.header().getSize() / 6.)); // Since there are 12 jets reported per event (see CMS IN-2013/005)
27 
28  // Find the first and last BXs
29  int firstBX = -(ceil((double)nBX/2.)-1);
30  int lastBX;
31  if (nBX % 2 == 0) {
32  lastBX = ceil((double)nBX/2.);
33  } else {
34  lastBX = ceil((double)nBX/2.)-1;
35  }
36 
37  auto res_ = static_cast<L1TObjectCollections*>(coll)->getJets();
38  res_->setBXRange(firstBX, lastBX);
39 
40  LogDebug("L1T") << "nBX = " << nBX << " first BX = " << firstBX << " lastBX = " << lastBX;
41 
42  // Initialise index
43  int unsigned i = 0;
44 
45  // Loop over multiple BX and then number of jets filling jet collection
46  for (int bx=firstBX; bx<=lastBX; bx++){
47  for (unsigned nJet=0; nJet < 6 && nJet < block.header().getSize(); nJet++){
48  uint32_t raw_data = block.payload()[i++];
49 
50  if (raw_data == 0)
51  continue;
52 
53  l1t::Jet jet = l1t::Jet();
54 
55  jet.setHwPt(raw_data & 0x7FF);
56 
57  int abs_eta = (raw_data >> 11) & 0x7F;
58  if ((raw_data >> 18) & 0x1) {
59  jet.setHwEta(-1 * (128-abs_eta));
60  } else {
61  jet.setHwEta(abs_eta);
62  }
63 
64  jet.setHwPhi((raw_data >> 19) & 0xFF);
65  jet.setHwQual((raw_data >> 27) & 0x7); // Assume 3 bits for now? Leaves 2 bits spare
66 
67  LogDebug("L1T") << "Jet: eta " << jet.hwEta() << " phi " << jet.hwPhi() << " pT " << jet.hwPt() << " qual " << jet.hwQual() << " bx " << bx;
68 
69  res_->push_back(bx,jet);
70  }
71  }
72 
73  return true;
74  }
75 }
76 }
77 
#define LogDebug(id)
int i
Definition: DBlmapReader.cc:9
unsigned int getID() const
Definition: Block.h:22
BlockHeader header() const
Definition: Block.h:56
void setHwQual(int qual)
Definition: L1Candidate.cc:64
int hwPhi() const
Definition: L1Candidate.cc:79
std::vector< uint32_t > payload() const
Definition: Block.h:57
#define DEFINE_L1T_UNPACKER(type)
Definition: Unpacker.h:31
Definition: Jet.h:13
int hwEta() const
Definition: L1Candidate.cc:74
int hwQual() const
Definition: L1Candidate.cc:89
virtual bool unpack(const Block &block, UnpackerCollections *coll) override
Definition: JetUnpacker.cc:21
void setHwPhi(int phi)
Definition: L1Candidate.cc:54
JetCorrectorParametersCollection coll
Definition: classes.h:10
int hwPt() const
Definition: L1Candidate.cc:69
void setHwPt(int pt)
Definition: L1Candidate.cc:44
unsigned int getSize() const
Definition: Block.h:23
void setHwEta(int eta)
Definition: L1Candidate.cc:49