FastL1CaloSim.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    L1CaloTriggerProducer
// Class:      FastL1CaloSim
// 
//
// Original Author:  Chi Nhan Nguyen
//         Created:  Mon Feb 19 13:25:24 CST 2007
// $Id: FastL1CaloSim.cc,v 1.13 2009/03/28 14:40:52 chinhan Exp $
//
//

#include "FastSimulation/L1CaloTriggerProducer/interface/FastL1CaloSim.h"


//
// constructors and destructor
//
FastL1CaloSim::FastL1CaloSim(const edm::ParameterSet& iConfig)
{
  //register your products
  /* old labels
  produces<l1extra::L1EtMissParticle>("MET");
  produces<l1extra::L1JetParticleCollection>("TauJets");
  produces<l1extra::L1JetParticleCollection>("CenJets");
  produces<l1extra::L1JetParticleCollection>("ForJets");
  produces<l1extra::L1EmParticleCollection>("Egammas");
  produces<l1extra::L1EmParticleCollection>("isoEgammas");
  */

  //produces<l1extra::L1EtMissParticle>();
  produces<l1extra::L1EtMissParticleCollection>("MET");
  produces<l1extra::L1JetParticleCollection>("Tau");
  produces<l1extra::L1JetParticleCollection>("Central");
  produces<l1extra::L1JetParticleCollection>("Forward");
  produces<l1extra::L1EmParticleCollection>("NonIsolated");
  produces<l1extra::L1EmParticleCollection>("Isolated");
  produces<l1extra::L1MuonParticleCollection>(); // muon is dummy for L1extraParticleMap!

  m_AlgorithmSource = iConfig.getParameter<std::string>("AlgorithmSource");

  // No BitInfos for release versions
  m_DoBitInfo = iConfig.getParameter<bool>("DoBitInfo");
  if (m_DoBitInfo)
    produces<FastL1BitInfoCollection>("L1BitInfos");

  m_FastL1GlobalAlgo = new FastL1GlobalAlgo(iConfig);
}

FastL1CaloSim::~FastL1CaloSim()
{
  delete m_FastL1GlobalAlgo;
}


//
// member functions
//

// ------------ method called to produce the data  ------------
void
FastL1CaloSim::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  //edm::LogInfo("FastL1CaloSim::produce()");

  if (m_AlgorithmSource == "RecHits") {
    m_FastL1GlobalAlgo->FillL1Regions(iEvent, iSetup);
    m_FastL1GlobalAlgo->FillEgammas(iEvent);
  } else if (m_AlgorithmSource == "TrigPrims") {
    m_FastL1GlobalAlgo->FillL1RegionsTP(iEvent,iSetup);
    m_FastL1GlobalAlgo->FillEgammasTP(iEvent);
  } else {
    std::cerr<<"AlgorithmSource not valid: "<<m_AlgorithmSource<<std::endl;
    return;
  }
  m_FastL1GlobalAlgo->FillMET(iEvent); // using CaloTowers
  //m_FastL1GlobalAlgo->FillMET();     // using Regions
  m_FastL1GlobalAlgo->FillJets(iSetup);
  
  if (m_DoBitInfo)
    m_FastL1GlobalAlgo->FillBitInfos();

  //std::auto_ptr<l1extra::L1EtMissParticle> METResult(new l1extra::L1EtMissParticle);
  std::auto_ptr<l1extra::L1EtMissParticleCollection> METResult(new l1extra::L1EtMissParticleCollection);
  std::auto_ptr<l1extra::L1JetParticleCollection> TauJetResult(new l1extra::L1JetParticleCollection);
  std::auto_ptr<l1extra::L1JetParticleCollection> CenJetResult(new l1extra::L1JetParticleCollection);
  std::auto_ptr<l1extra::L1JetParticleCollection> ForJetResult(new l1extra::L1JetParticleCollection);
  std::auto_ptr<l1extra::L1EmParticleCollection> EgammaResult(new l1extra::L1EmParticleCollection);
  std::auto_ptr<l1extra::L1EmParticleCollection> isoEgammaResult(new l1extra::L1EmParticleCollection);
  // muon is dummy!
  std::auto_ptr<l1extra::L1MuonParticleCollection> muonDummy(new l1extra::L1MuonParticleCollection);
  //
  //*METResult = m_FastL1GlobalAlgo->getMET();
  for (int i=0; i<(int)m_FastL1GlobalAlgo->getMET().size(); i++) {
    METResult->push_back(m_FastL1GlobalAlgo->getMET().at(i));
  }
  for (int i=0; i<std::min(4,(int)m_FastL1GlobalAlgo->getTauJets().size()); i++) {
    TauJetResult->push_back(m_FastL1GlobalAlgo->getTauJets().at(i));
  }
  for (int i=0; i<std::min(4,(int)m_FastL1GlobalAlgo->getCenJets().size()); i++) {
    CenJetResult->push_back(m_FastL1GlobalAlgo->getCenJets().at(i));
  }
  for (int i=0; i<std::min(4,(int)m_FastL1GlobalAlgo->getForJets().size()); i++) {
    ForJetResult->push_back(m_FastL1GlobalAlgo->getForJets().at(i));
  }
  for (int i=0; i<std::min(4,(int)m_FastL1GlobalAlgo->getEgammas().size()); i++) {
    EgammaResult->push_back(m_FastL1GlobalAlgo->getEgammas().at(i));
  }
  for (int i=0; i<std::min(4,(int)m_FastL1GlobalAlgo->getisoEgammas().size()); i++) {
    isoEgammaResult->push_back(m_FastL1GlobalAlgo->getisoEgammas().at(i));
  }

  if (m_DoBitInfo) {
  std::auto_ptr<FastL1BitInfoCollection> L1BitInfoResult(new FastL1BitInfoCollection);
    for (int i=0; i<(int)m_FastL1GlobalAlgo->getBitInfos().size(); i++) {
      L1BitInfoResult->push_back(m_FastL1GlobalAlgo->getBitInfos().at(i));
    }
    iEvent.put(L1BitInfoResult,"L1BitInfos");
  }

  // put the collections into the event
  /* old labels
  iEvent.put(METResult,"MET");
  iEvent.put(TauJetResult,"TauJets");
  iEvent.put(CenJetResult,"CenJets");
  iEvent.put(ForJetResult,"ForJets");
  iEvent.put(EgammaResult,"Egammas");
  iEvent.put(isoEgammaResult,"isoEgammas");
  */
  iEvent.put(METResult);
  iEvent.put(TauJetResult,"Tau");
  iEvent.put(CenJetResult,"Central");
  iEvent.put(ForJetResult,"Forward");
  iEvent.put(EgammaResult,"NonIsolated");
  iEvent.put(isoEgammaResult,"Isolated");
  iEvent.put(muonDummy);

}

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