EcalExclusiveTrigFilter.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    EcalExclusiveTrigFilter
// Class:      EcalExclusiveTrigFilter
// 
//
// Original Author:  Seth COOPER
//         Created:  Thu May 22 11:40:12 CEST 2008
//
//

#include "CaloOnlineTools/EcalTools/plugins/EcalExclusiveTrigFilter.h"

using namespace cms;
using namespace edm;
using namespace std;

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
EcalExclusiveTrigFilter::EcalExclusiveTrigFilter(const edm::ParameterSet& iConfig) :
  l1GTReadoutRecTag_ (iConfig.getUntrackedParameter<std::string>("L1GlobalReadoutRecord","gtDigis"))
{
   //now do what ever initialization is needed

}


EcalExclusiveTrigFilter::~EcalExclusiveTrigFilter()
{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}


//
// member functions
//

// ------------ method called on each new Event  ------------
bool
EcalExclusiveTrigFilter::filter(edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  // get the GMTReadoutCollection
  Handle<L1MuGMTReadoutCollection> gmtrc_handle; 
  iEvent.getByLabel(l1GTReadoutRecTag_,gmtrc_handle);
  L1MuGMTReadoutCollection const* gmtrc = gmtrc_handle.product();
  if (!(gmtrc_handle.isValid())) 
  {
    LogWarning("EcalExclusiveTrigFilter") << "l1MuGMTReadoutCollection" << " not available";
    //return;
  }

  // get hold of L1GlobalReadoutRecord
  Handle<L1GlobalTriggerReadoutRecord> L1GTRR;
  iEvent.getByLabel(l1GTReadoutRecTag_,L1GTRR);
  bool isEcalL1 = false;
  const unsigned int sizeOfDecisionWord(L1GTRR->decisionWord().size());
  if (!(L1GTRR.isValid()))
  {
    LogWarning("EcalExclusiveTrigFilter") << l1GTReadoutRecTag_ << " not available";
    //return;
  }
  else if(sizeOfDecisionWord<128)
  {
    LogWarning("EcalExclusiveTrigFilter") << "size of L1 decisionword is " << sizeOfDecisionWord
      << "; L1 Ecal triggering bits not available";
  }
  else
  {
    l1Names_.resize(sizeOfDecisionWord);
    l1Accepts_.resize(sizeOfDecisionWord);
    for (unsigned int i=0; i!=sizeOfDecisionWord; ++i) {
      l1Accepts_[i]=0;
      l1Names_[i]="NameNotAvailable";
    }
    for (unsigned int i=0; i!=sizeOfDecisionWord; ++i) {
      if (L1GTRR->decisionWord()[i])
      {
        l1Accepts_[i]++;
        //cout << "L1A bit: " << i << endl;
      }
    }

    if(l1Accepts_[14] || l1Accepts_[15] || l1Accepts_[16] || l1Accepts_[17]
        || l1Accepts_[18] || l1Accepts_[19] || l1Accepts_[20])
      isEcalL1 = true;
    if(l1Accepts_[73] || l1Accepts_[74] || l1Accepts_[75] || l1Accepts_[76]
        || l1Accepts_[77] || l1Accepts_[78])
      isEcalL1 = true;
  }

  bool isRPCL1 = false;
  bool isDTL1 = false;
  bool isCSCL1 = false;
  bool isHCALL1 = false;

  std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
  std::vector<L1MuGMTReadoutRecord>::const_iterator igmtrr;

    for(igmtrr=gmt_records.begin(); igmtrr!=gmt_records.end(); igmtrr++) {

      std::vector<L1MuRegionalCand>::const_iterator iter1;
      std::vector<L1MuRegionalCand> rmc;

      //DT triggers
      int idt = 0;
      rmc = igmtrr->getDTBXCands();
      for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
        if ( !(*iter1).empty() ) {
          idt++;
        }
      }

      //if(idt>0) std::cout << "Found " << idt << " valid DT candidates in bx wrt. L1A = " 
      //  << igmtrr->getBxInEvent() << std::endl;
      if(igmtrr->getBxInEvent()==0 && idt>0) isDTL1 = true;

      //RPC triggers
      int irpcb = 0;
      rmc = igmtrr->getBrlRPCCands();
      for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
        if ( !(*iter1).empty() ) {
          irpcb++;
        }
      }

      //if(irpcb>0) std::cout << "Found " << irpcb << " valid RPC candidates in bx wrt. L1A = " 
      //  << igmtrr->getBxInEvent() << std::endl;
      if(igmtrr->getBxInEvent()==0 && irpcb>0) isRPCL1 = true;

      //CSC Triggers
      int icsc = 0;
      rmc = igmtrr->getCSCCands();
      for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
        if ( !(*iter1).empty() ) {
          icsc++;
        }
      }
      //if(icsc>0) std::cout << "Found " << icsc << " valid CSC candidates in bx wrt. L1A = " 
      //    //  << igmtrr->getBxInEvent() << std::endl;
      if(igmtrr->getBxInEvent()==0 && icsc>0) isCSCL1 = true;
    }

    L1GlobalTriggerReadoutRecord const* gtrr = L1GTRR.product();

    for(int ibx=-1; ibx<=1; ibx++) {
      bool hcal_top = false;
      bool hcal_bot = false;
      const L1GtPsbWord psb = gtrr->gtPsbWord(0xbb0d,ibx);
      std::vector<int> valid_phi;
      if((psb.aData(4)&0x3f) >= 1) {valid_phi.push_back( (psb.aData(4)>>10)&0x1f ); }
      if((psb.bData(4)&0x3f) >= 1) {valid_phi.push_back( (psb.bData(4)>>10)&0x1f ); }
      if((psb.aData(5)&0x3f) >= 1) {valid_phi.push_back( (psb.aData(5)>>10)&0x1f ); }
      if((psb.bData(5)&0x3f) >= 1) {valid_phi.push_back( (psb.bData(5)>>10)&0x1f ); }
      std::vector<int>::const_iterator iphi;
      for(iphi=valid_phi.begin(); iphi!=valid_phi.end(); iphi++) {
        //std::cout << "Found HCAL mip with phi=" << *iphi << " in bx wrt. L1A = " << ibx << std::endl;
        if(*iphi<9) hcal_top=true;
        if(*iphi>8) hcal_bot=true;
      }
      if(ibx==0 && hcal_top && hcal_bot) isHCALL1=true;
    }

    std::cout << "**** Trigger Source ****" << std::endl;
    if(isDTL1) std::cout << "DT" << std::endl;
    if(isRPCL1) std::cout << "RPC" << std::endl;
    if(isCSCL1) cout << "CSC" << endl;
    if(isHCALL1) std::cout << "HCAL" << std::endl;
    if(isEcalL1) std::cout << "ECAL" << std::endl;
    std::cout << "************************" << std::endl;

    return(isEcalL1 && !isDTL1 && !isRPCL1 && !isCSCL1 && !isHCALL1);

}