/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/JetUtils/JetUtils/JetCaloCellQualityUtils.h

// This file is -*- c++ -*-
#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

---