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>("lut_dump_flag");
  my_lut_btic = pset.getUntrackedParameter<int>("lut_btic");
}

DTTrigProd::~DTTrigProd(){

  if (my_trig) delete my_trig;

}

// void DTTrigProd::beginJob(const EventSetup & iEventSetup){


//   // get DTConfigManager
//   // ESHandle< DTConfigManager > confManager ;
//   // iEventSetup.get< DTConfigManagerRcd >().get( confManager ) ;
//   // my_BXoffset = confManager->getBXOffset();

//   if (my_debug)
//     cout << "[DTTrigProd] DTTrig istance Created" << endl;

// }

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

   if(my_debug)
    cout << "DTTrigProd::beginRun  " << iRun.id().run() << endl;

   using namespace edm;

   ESHandle< DTConfigManager > dtConfig ;
   iEventSetup.get< DTConfigManagerRcd >().get( dtConfig ) ;

   if(my_debug)
   {
    cout << "DTConfigManagerRcd : Print some Config stuff" << endl;
    DTBtiId btiid(1,1,1,1,1);
    DTTracoId tracoid(1,1,1,1);
    DTChamberId chid(1,1,1);
    DTSectCollId scid(1,1);
    cout    << "BtiMap & TracoMap Size for chamber (1,1,1):" << dtConfig->getDTConfigBtiMap(chid).size() 
        << " " << dtConfig->getDTConfigTracoMap(chid).size() << endl;

    dtConfig->getDTConfigBti(btiid)->print();
    dtConfig->getDTConfigTraco(tracoid)->print();
    dtConfig->getDTConfigTSTheta(chid)->print();
    dtConfig->getDTConfigTSPhi(chid)->print();
        // 100511 SV LUTs will NOT be configurated from cfg if not found in DB
        //dtConfig->getDTConfigLUTs(chid)->print(); 
  }


  if (!my_trig) {
    my_trig = new DTTrig(my_params);
    my_trig->createTUs(iEventSetup);
    if (my_debug)
      cout << "[DTTrigProd] TU's Created" << endl;

    //SV 090903 check luts
    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){

  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);

}