JetCaloCellQualityUtils.h

Go to the documentation of this file.

// This file is -*- c++ -*-
 
/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef JETUTILS_JETCALOCELLQUALHELPER_H
#define JETUTILS_JETCALOCELLQUALHELPER_H
 
#include "JetUtils/JetCaloQualityUtils.h"
 
 
 
class CaloCell;
 
namespace jet {

  class JetCaloCellCalculator : virtual public JetCaloCalculator {
  public:
    virtual bool processCell(const CaloCell* /*cell*/, weight_t /*w*/)  =0;
    virtual double operator()(const xAOD::Jet* jet, xAOD::JetConstitScale s=xAOD::UncalibratedJetConstituent);
 
    virtual JetCaloCellCalculator* clone()const = 0;
 
  };
 

#define DECLARECELLCALC( baseclass, derivedclass, inits ) \
  class derivedclass : public baseclass ,  public JetCaloCellCalculator  {                            \
  public:                                                               \
    derivedclass() :  baseclass() {inits;}              \
    virtual bool processCell(const CaloCell* , weight_t )  ;             \
    JetCaloCellCalculator* clone()const {return new derivedclass(*this);}    \
  }
    
 
  // below we declare cell based calo calculators inheriting the cluster based calculators.
  DECLARECELLCALC(JetCalcnLeadingCells, JetCalcnLeadingCells_fromCells, m_id = JetAttribute::N90Cells; );
  DECLARECELLCALC(JetCalcOutOfTimeEnergyFraction, JetCalcOutOfTimeEnergyFraction_fromCells, );
  DECLARECELLCALC(JetCalcTimeCells, JetCalcTimeCells_fromCells, );
  DECLARECELLCALC(JetCalcAverageLArQualityF, JetCalcAverageLArQualityF_fromCells, );
  DECLARECELLCALC(JetCalcQualityHEC, JetCalcQualityHEC_fromCells, );
  DECLARECELLCALC(JetCalcQuality, JetCalcQuality_fromCells, );
  DECLARECELLCALC(JetCalcNegativeEnergy, JetCalcNegativeEnergy_fromCells, );
  DECLARECELLCALC(JetCalcCentroid, JetCalcCentroid_fromCells, );
  //  DECLARECELLCALC(JetCalcBadCellsFrac, JetCalcBadCellsFrac_fromCells, );
 

  class JetCaloCellCalculations : public JetCaloCalculations  {
  public:
    
    JetCaloCellCalculations(bool ownCalc=true) : JetCaloCalculations(ownCalc) {}

    virtual std::vector<double> process(const xAOD::Jet* jet) const ;    
 
    void addCellCalculator(JetCaloCellCalculator* c);
 
  protected:
    std::vector<JetCaloCellCalculator*> m_cellcalculators;
 
  };
 
 
 
}
 
#endif