L1TMuonCaloSumProducer.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    L1TMuonCaloSumProducer
// Class:      L1TMuonCaloSumProducer
//
//
// Original Author:  Joschka Philip Lingemann,40 3-B01,+41227671598,
//         Created:  Thu Oct  3 10:12:30 CEST 2013
// $Id$
//
//

// system include files
#include <memory>
#include <fstream>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/L1TMuon/interface/MuonCaloSumFwd.h"
#include "DataFormats/L1TMuon/interface/MuonCaloSum.h"
#include "DataFormats/L1TCalorimeter/interface/CaloTower.h"

#include "TMath.h"
#include "TRandom3.h"

//
// class declaration
//
using namespace l1t;

class L1TMuonCaloSumProducer : public edm::stream::EDProducer<> {
public:
  explicit L1TMuonCaloSumProducer(const edm::ParameterSet&);

  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

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

  edm::EDGetTokenT<CaloTowerBxCollection> m_caloTowerToken;
  edm::InputTag m_caloLabel;
};

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
L1TMuonCaloSumProducer::L1TMuonCaloSumProducer(const edm::ParameterSet& iConfig) {
  //register your inputs:
  m_caloLabel = iConfig.getParameter<edm::InputTag>("caloStage2Layer2Label");
  m_caloTowerToken = consumes<CaloTowerBxCollection>(m_caloLabel);
  //register your products
  produces<MuonCaloSumBxCollection>("TriggerTowerSums");
  produces<MuonCaloSumBxCollection>("TriggerTower2x2s");
}

//
// member functions
//

// ------------ method called to produce the data  ------------
void L1TMuonCaloSumProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  using namespace edm;
  std::unique_ptr<MuonCaloSumBxCollection> towerSums(std::make_unique<MuonCaloSumBxCollection>());
  std::unique_ptr<MuonCaloSumBxCollection> tower2x2s(std::make_unique<MuonCaloSumBxCollection>());

  edm::Handle<CaloTowerBxCollection> caloTowers;

  if (iEvent.getByToken(m_caloTowerToken, caloTowers)) {
    int detamax = 4;
    int dphimax = 4;

    const int iFirstBx = caloTowers->getFirstBX();
    const int iLastBx = caloTowers->getLastBX();

    // set BX range for sums
    towerSums->setBXRange(iFirstBx, iLastBx);
    tower2x2s->setBXRange(iFirstBx, iLastBx);

    for (int bx = iFirstBx; bx <= iLastBx; ++bx) {
      std::map<int, MuonCaloSum> sums;
      std::map<int, MuonCaloSum> regs;

      for (auto it = caloTowers->begin(bx); it != caloTowers->end(bx); ++it) {
        const CaloTower& twr = *it;
        int hwEta = twr.hwEta();
        if (std::abs(hwEta) > 27) {
          continue;
        }
        int hwPt = twr.hwPt();
        if (hwPt < 1) {
          continue;
        }
        int hwPhi = twr.hwPhi();

        // calculating tower2x2s
        int ieta2x2 = (hwEta + 27) / 2;
        int iphi2x2 = hwPhi / 2;
        int muon_idx = iphi2x2 * 28 + ieta2x2;
        if (regs.count(muon_idx) == 0) {
          regs[muon_idx] = MuonCaloSum(hwPt, iphi2x2, ieta2x2, muon_idx);
        } else {
          regs.at(muon_idx).setEtBits(regs.at(muon_idx).etBits() + hwPt);
        }

        // std::cout << "iphi; phi " << hwPhi << "; " << phi << " .. ieta; eta" << hwEta << "; " << twr.eta() << std::endl;

        // calculating towerSums
        int ietamax = hwEta + detamax + 1;
        for (int ieta = hwEta - detamax; ieta < ietamax; ++ieta) {
          if (std::abs(ieta) > 27) {
            continue;
          }
          int ietamu = (ieta + 27) / 2;
          int iphimax = hwPhi + dphimax + 1;
          for (int iphi = hwPhi - dphimax; iphi < iphimax; ++iphi) {
            int iphiwrapped = iphi;
            if (iphiwrapped < 0) {
              iphiwrapped += 72;
            } else if (iphiwrapped > 71) {
              iphiwrapped -= 72;
            }
            int iphimu = iphiwrapped / 2;
            int idxmu = iphimu * 28 + ietamu;
            if (sums.count(idxmu) == 0) {
              sums[idxmu] = MuonCaloSum(hwPt, iphimu, ietamu, idxmu);
            } else {
              sums.at(idxmu).setEtBits(sums.at(idxmu).etBits() + hwPt);
            }
          }
        }
      }

      // fill towerSums output collection for this BX
      for (auto it = sums.begin(); it != sums.end(); ++it) {
        if (it->second.etBits() > 0) {
          MuonCaloSum sum = MuonCaloSum(it->second);
          // convert Et to correct scale:
          if (sum.etBits() > 31) {
            sum.setEtBits(31);
          }
          towerSums->push_back(bx, sum);
        }
      }
      // fill tower2x2s output collection for this BX
      for (auto it = regs.begin(); it != regs.end(); ++it) {
        if (it->second.etBits() > 0) {
          tower2x2s->push_back(bx, it->second);
        }
      }
    }
  } else {
    LogWarning("GlobalMuon") << "CaloTowers not found. Producing empty collections." << std::endl;
  }

  iEvent.put(std::move(towerSums), "TriggerTowerSums");
  iEvent.put(std::move(tower2x2s), "TriggerTower2x2s");
}

// ------------ method fills 'descriptions' with the allowed parameters for the module  ------------
void L1TMuonCaloSumProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

//define this as a plug-in
DEFINE_FWK_MODULE(L1TMuonCaloSumProducer);