L1AnalysisCaloTP.cc

Go to the documentation of this file.

#include "L1Trigger/L1TNtuples/interface/L1AnalysisCaloTP.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "CondFormats/L1TObjects/interface/L1CaloEcalScale.h"
#include "CondFormats/L1TObjects/interface/L1CaloHcalScale.h"


L1Analysis::L1AnalysisCaloTP::L1AnalysisCaloTP():verbose_(false)
{
}

L1Analysis::L1AnalysisCaloTP::L1AnalysisCaloTP(bool verbose)
{
  verbose_ = verbose;
}

L1Analysis::L1AnalysisCaloTP::~L1AnalysisCaloTP()
{

}

void L1Analysis::L1AnalysisCaloTP::SetHCAL( const HcalTrigPrimDigiCollection& hcalTPs ) {

  if (verbose_) edm::LogInfo("L1Ntuple") << "HCAL TPs : " << hcalTPs.size() << std::endl;

  for (unsigned i=0; i<hcalTPs.size(); ++i) {

    short ieta = (short) hcalTPs[i].id().ieta(); 
    unsigned short absIeta = (unsigned short) abs(ieta);
    short sign = ieta/absIeta;

    unsigned short cal_iphi = (unsigned short) hcalTPs[i].id().iphi();
    unsigned short iphi = (72 + 18 - cal_iphi) % 72;
    if (absIeta >= 29) {  // special treatment for HF
      iphi = iphi/4;
    }

    unsigned short compEt = hcalTPs[i].SOI_compressedEt();
    double et = 0.;
    if (hcalScale_!=nullptr) et = hcalScale_->et( compEt, absIeta, sign );

    unsigned short fineGrain = (unsigned short) hcalTPs[i].SOI_fineGrain();

    tp_.hcalTPieta.push_back( ieta );
    tp_.hcalTPCaliphi.push_back( cal_iphi );
    tp_.hcalTPiphi.push_back( iphi );
    tp_.hcalTPet.push_back( et );
    tp_.hcalTPcompEt.push_back( compEt );
    tp_.hcalTPfineGrain.push_back( fineGrain );
    tp_.nHCALTP++;

  }

}

void L1Analysis::L1AnalysisCaloTP::SetECAL( const EcalTrigPrimDigiCollection& ecalTPs ) {
  
  if (verbose_) edm::LogInfo("L1Ntuple") << "ECAL TPs : " << ecalTPs.size() << std::endl;

  for (unsigned i=0; i<ecalTPs.size(); ++i) {

    short ieta = (short) ecalTPs[i].id().ieta(); 
    unsigned short absIeta = (unsigned short) abs(ieta);
    short sign = ieta/absIeta;
    
    unsigned short cal_iphi = (unsigned short) ecalTPs[i].id().iphi(); 
    unsigned short iphi = (72 + 18 - cal_iphi) % 72; // transform TOWERS (not regions) into local rct (intuitive) phi bins
    
    unsigned short compEt = ecalTPs[i].compressedEt();
    double et = 0.;
    if (ecalScale_!=nullptr) et = ecalScale_->et( compEt, absIeta, sign );
    
    unsigned short fineGrain = (unsigned short) ecalTPs[i].fineGrain();  // 0 or 1
    
    tp_.ecalTPieta.push_back( ieta );
    tp_.ecalTPCaliphi.push_back( cal_iphi );
    tp_.ecalTPiphi.push_back( iphi );
    tp_.ecalTPet.push_back( et );
    tp_.ecalTPcompEt.push_back( compEt );
    tp_.ecalTPfineGrain.push_back( fineGrain );
    tp_.nECALTP++;
    
  }  

  
}