---

CSCTFTrackBuilder.cc

Go to the documentation of this file.

#include <L1Trigger/CSCTrackFinder/src/CSCTFTrackBuilder.h>

#include <FWCore/MessageLogger/interface/MessageLogger.h>

#include <DataFormats/MuonDetId/interface/CSCDetId.h>
#include <DataFormats/MuonDetId/interface/CSCTriggerNumbering.h>
#include <DataFormats/L1DTTrackFinder/interface/L1MuDTChambPhContainer.h>

#include <DataFormats/L1CSCTrackFinder/interface/CSCTriggerContainer.h>
#include <L1Trigger/CSCTrackFinder/interface/CSCTFSectorProcessor.h>

#include "CondFormats/L1TObjects/interface/L1MuCSCTFConfiguration.h"
#include "CondFormats/DataRecord/interface/L1MuCSCTFConfigurationRcd.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include <sstream>
#include <stdlib.h>

CSCTFTrackBuilder::CSCTFTrackBuilder(const edm::ParameterSet& pset, bool TMB07,
                                     const L1MuTriggerScales* scales,
                                     const L1MuTriggerPtScale* ptScale ){
  my_dtrc = new CSCTFDTReceiver();

  for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
    {
      for(int s = CSCTriggerNumbering::minTriggerSectorId();
          s <= CSCTriggerNumbering::maxTriggerSectorId(); ++s)
        {
          // All SPs work with the same configuration (impossible to make it more exclusive in this framework)
          my_SPs[e-1][s-1] = new CSCTFSectorProcessor(e, s, pset, TMB07,
                                                      scales, ptScale);
        }
    }
}

void CSCTFTrackBuilder::initialize(const edm::EventSetup& c){
 //my_dtrc->initialize(c);
        for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
        {
                for(int s = CSCTriggerNumbering::minTriggerSectorId(); s <= CSCTriggerNumbering::maxTriggerSectorId(); ++s)
                {
                        my_SPs[e-1][s-1]->initialize(c);
                }
        }
}

CSCTFTrackBuilder::~CSCTFTrackBuilder()
{
  delete my_dtrc;
  my_dtrc = NULL;

  for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
    {
      for(int s = CSCTriggerNumbering::minTriggerSectorId();
          s <= CSCTriggerNumbering::maxTriggerSectorId(); ++s)
        {
          delete my_SPs[e-1][s-1];
          my_SPs[e-1][s-1] = NULL;
        }
    }
}

void CSCTFTrackBuilder::buildTracks(const CSCCorrelatedLCTDigiCollection* lcts, const L1MuDTChambPhContainer* dttrig,
                                    L1CSCTrackCollection* trkcoll, CSCTriggerContainer<csctf::TrackStub>* stubs_to_dt)
{
  std::vector<csc::L1Track> trks;
  CSCTriggerContainer<csctf::TrackStub> stub_list;

  CSCCorrelatedLCTDigiCollection::DigiRangeIterator Citer;

  for(Citer = lcts->begin(); Citer != lcts->end(); Citer++)
    {
      CSCCorrelatedLCTDigiCollection::const_iterator Diter = (*Citer).second.first;
      CSCCorrelatedLCTDigiCollection::const_iterator Dend = (*Citer).second.second;

          for(; Diter != Dend; Diter++)
        {
          csctf::TrackStub theStub((*Diter),(*Citer).first);
          stub_list.push_back(theStub);
        }
    }

  // Now we append the track stubs the the DT Sector Collector
  // after processing from the DT Receiver.

  CSCTriggerContainer<csctf::TrackStub> dtstubs = my_dtrc->process(dttrig);
  stub_list.push_many(dtstubs);

  // run each sector processor in the TF
  for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
  {
      for(int s = CSCTriggerNumbering::minTriggerSectorId();
        s <= CSCTriggerNumbering::maxTriggerSectorId(); ++s)
      {
         CSCTriggerContainer<csctf::TrackStub> current_e_s = stub_list.get(e, s);
         if(my_SPs[e-1][s-1]->run(current_e_s))
         {
           std::vector<csc::L1Track> theTracks = my_SPs[e-1][s-1]->tracks().get();
           trks.insert(trks.end(), theTracks.begin(), theTracks.end());
         }
         stubs_to_dt->push_many(my_SPs[e-1][s-1]->dtStubs()); // send stubs whether or not we find a track!!!
      }
  }

  // Now to combine tracks with their track stubs and send them off.
  trkcoll->resize(trks.size());
  std::vector<csc::L1Track>::const_iterator titr = trks.begin();
  L1CSCTrackCollection::iterator tcitr = trkcoll->begin();

  for(; titr != trks.end(); titr++)
    {
      tcitr->first = (*titr);
      std::vector<csctf::TrackStub> possible_stubs = my_SPs[titr->endcap()-1][titr->sector()-1]->filteredStubs();
      std::vector<csctf::TrackStub>::const_iterator tkstbs = possible_stubs.begin();

      int me1ID = titr->me1ID();
      int me2ID = titr->me2ID();
      int me3ID = titr->me3ID();
      int me4ID = titr->me4ID();
      int mb1ID = titr->mb1ID();

      for(; tkstbs != possible_stubs.end(); tkstbs++)
        {
          switch(tkstbs->station())
            {
            case 1:
              if((tkstbs->getMPCLink()
                  +(3*(CSCTriggerNumbering::triggerSubSectorFromLabels(CSCDetId(tkstbs->getDetId().rawId())) - 1))) == me1ID && me1ID != 0)
                {
                  tcitr->second.insertDigi(CSCDetId(tkstbs->getDetId().rawId()), *(tkstbs->getDigi()));
                }
              break;
            case 2:
              if(tkstbs->getMPCLink() == me2ID && me2ID != 0)
                {
                  tcitr->second.insertDigi(CSCDetId(tkstbs->getDetId().rawId()), *(tkstbs->getDigi()));
                }
              break;
            case 3:
              if(tkstbs->getMPCLink() == me3ID && me3ID != 0)
                {
                  tcitr->second.insertDigi(CSCDetId(tkstbs->getDetId().rawId()), *(tkstbs->getDigi()));
                }
              break;
            case 4:
              if(tkstbs->getMPCLink() == me4ID && me4ID != 0)
                {
                  tcitr->second.insertDigi(CSCDetId(tkstbs->getDetId().rawId()), *(tkstbs->getDigi()));
                }
              break;
            case 5:
              if(tkstbs->getMPCLink() == mb1ID && mb1ID != 0)
              {
              }
              break;
            default:
              edm::LogWarning("CSCTFTrackBuilder::buildTracks()") << "SERIOUS ERROR: STATION " << tkstbs->station() << " NOT IN RANGE [1,5]\n";
            };
        }
      tcitr++; // increment to next track in the collection
    }

}

---