FakeTBEventHeaderProducer.cc

Go to the documentation of this file.

/*
 * \file FakeTBEventHeaderProducer.cc
 *
 *
 */

#include "SimG4CMS/EcalTestBeam/interface/FakeTBEventHeaderProducer.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"

using namespace cms;
using namespace std;


FakeTBEventHeaderProducer::FakeTBEventHeaderProducer(const edm::ParameterSet& ps) {
  produces<EcalTBEventHeader>();
  ecalTBInfoLabel_ = ps.getUntrackedParameter<string>("EcalTBInfoLabel","SimEcalTBG4Object");

}

 
FakeTBEventHeaderProducer::~FakeTBEventHeaderProducer() 
{
}

 void FakeTBEventHeaderProducer::produce(edm::Event & event, const edm::EventSetup& eventSetup)
{
  auto_ptr<EcalTBEventHeader> product(new EcalTBEventHeader());

  // get the vertex information from the event

  const PEcalTBInfo* theEcalTBInfo=0;
  edm::Handle<PEcalTBInfo> EcalTBInfo;
  event.getByLabel(ecalTBInfoLabel_,EcalTBInfo);
  if (EcalTBInfo.isValid()){
    theEcalTBInfo = EcalTBInfo.product(); 
  } else {
    edm::LogError("FakeTBEventHeaderProducer") << "Error! can't get the product " << ecalTBInfoLabel_.c_str() ;
  }
  
  if (!theEcalTBInfo)
    return;
  
  // 64 bits event ID in CMSSW converted to EcalTBEventHeader ID
  int evtid = (int)event.id().event();
  product->setEventNumber(evtid);
  product->setRunNumber(event.id().run());
  product->setBurstNumber(1);
  product->setTriggerMask(0x1);
  product->setCrystalInBeam(EBDetId(1,theEcalTBInfo->nCrystal(),EBDetId::SMCRYSTALMODE));
  
//   LogDebug("FakeTBHeader") << (*product);
//   LogDebug("FakeTBHeader") << (*product).eventType();
//   LogDebug("FakeTBHeader") << (*product).crystalInBeam();
  event.put(product);
  
}