DTTrigProd.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//
//--------------------------------------------------

// This class's header
#include "L1Trigger/DTTrigger/interface/DTTrigProd.h"

// Framework related classes
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Run.h"

#include "L1TriggerConfig/DTTPGConfig/interface/DTConfigManager.h"
#include "L1TriggerConfig/DTTPGConfig/interface/DTConfigManagerRcd.h"

// Data Formats classes
#include "L1Trigger/DTSectorCollector/interface/DTSectCollPhSegm.h"
#include "L1Trigger/DTSectorCollector/interface/DTSectCollThSegm.h"
#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambPhContainer.h"
#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambPhDigi.h"
#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambThContainer.h"
#include "DataFormats/L1DTTrackFinder/interface/L1MuDTChambThDigi.h"

// Collaborating classes
#include <iostream>

using namespace edm;
using namespace std;

// DataFormats interface
typedef vector<DTSectCollPhSegm>  SectCollPhiColl;
typedef SectCollPhiColl::const_iterator SectCollPhiColl_iterator;
typedef vector<DTSectCollThSegm>  SectCollThetaColl;
typedef SectCollThetaColl::const_iterator SectCollThetaColl_iterator;

DTTrigProd::DTTrigProd(const ParameterSet& pset) : my_trig(0) {
  
  produces<L1MuDTChambPhContainer>();
  produces<L1MuDTChambThContainer>();

  my_debug = pset.getUntrackedParameter<bool>("debug");
  my_DTTFnum = pset.getParameter<bool>("DTTFSectorNumbering");
  my_params = pset;

  my_lut_dump_flag = pset.getUntrackedParameter<bool>("lutDumpFlag");
  my_lut_btic = pset.getUntrackedParameter<int>("lutBtic");

}

DTTrigProd::~DTTrigProd(){

  if (my_trig) delete my_trig;

}

void DTTrigProd::beginRun(edm::Run& iRun, const edm::EventSetup& iEventSetup) {

  if(my_debug)
    cout << "DTTrigProd::beginRun  " << iRun.id().run() << endl;
  
  ESHandle< DTConfigManager > dtConfig ;
  iEventSetup.get< DTConfigManagerRcd >().get( dtConfig ) ;

  my_CCBValid = dtConfig->CCBConfigValidity();

  if (!my_trig) {
    my_trig = new DTTrig(my_params);
    my_trig->createTUs(iEventSetup);
    if (my_debug)
      cout << "[DTTrigProd] TU's Created" << endl;
    
    if(my_lut_dump_flag) {
      cout << "Dumping luts...." << endl;
      my_trig->dumpLuts(my_lut_btic, dtConfig.product());
    }   
  }


  
}


void DTTrigProd::produce(Event & iEvent, const EventSetup& iEventSetup){

  // SV check if CCB configuration is valid, otherwise just produce empty collections
  if(!my_CCBValid)  {
    if (my_debug)
      cout << "[DTTrigProd] CCB configuration is not valid for this run, emulator will not run ! " << endl;
  }  else  {
    my_trig->triggerReco(iEvent,iEventSetup);
    my_BXoffset = my_trig->getBXOffset();
  
    if (my_debug)
      cout << "[DTTrigProd] Trigger algorithm run for " <<iEvent.id() << endl;
  
    // Convert Phi Segments
    SectCollPhiColl myPhiSegments;
    myPhiSegments = my_trig->SCPhTrigs();
    vector<L1MuDTChambPhDigi> outPhi;

    SectCollPhiColl_iterator SCPCend = myPhiSegments.end();
    for (SectCollPhiColl_iterator it=myPhiSegments.begin();it!=SCPCend;++it){
      int step = (*it).step() - my_BXoffset; // Shift correct BX to 0 (needed for DTTF data processing)
      int sc_sector = (*it).SCId().sector();
      if (my_DTTFnum == true) sc_sector--; // Modified for DTTF numbering [0-11]
        outPhi.push_back(L1MuDTChambPhDigi(step,
                       (*it).ChamberId().wheel(),
                       sc_sector,
                       (*it).ChamberId().station(),
                       (*it).phi(),
                       (*it).phiB(),
                       (*it).code(),
                       !(*it).isFirst(),
                       0
                       ));
    }

    // Convert Theta Segments
    SectCollThetaColl myThetaSegments;
    myThetaSegments = my_trig->SCThTrigs();
    vector<L1MuDTChambThDigi> outTheta;
  
    SectCollThetaColl_iterator SCTCend = myThetaSegments.end();
    for (SectCollThetaColl_iterator it=myThetaSegments.begin();it!=SCTCend;++it){
      int pos[7], qual[7];
      for (int i=0; i<7; i++){
        pos[i] =(*it).position(i);
        qual[i]=(*it).quality(i);
      }
      int step =(*it).step() - my_BXoffset; // Shift correct BX to 0 (needed for DTTF data processing)
      int sc_sector =  (*it).SCId().sector();
      if (my_DTTFnum == true) sc_sector--; // Modified for DTTF numbering [0-11]
      outTheta.push_back(L1MuDTChambThDigi( step,
                     (*it).ChamberId().wheel(),
                     sc_sector,
                     (*it).ChamberId().station(),
                     pos,
                     qual
                     ));
    }

   // Write everything into the event
   std::auto_ptr<L1MuDTChambPhContainer> resultPhi (new L1MuDTChambPhContainer);
   resultPhi->setContainer(outPhi);
   iEvent.put(resultPhi);
   std::auto_ptr<L1MuDTChambThContainer> resultTheta (new L1MuDTChambThContainer);
   resultTheta->setContainer(outTheta);
   iEvent.put(resultTheta);
  }
}