/data/refman/pasoursint/CMSSW_4_1_8_patch9/src/RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h

Go to the documentation of this file.

#ifndef RecoLocalCalo_EcalRecAlgos_EcalSeverityLevelAlgo_hh
#define RecoLocalCalo_EcalRecAlgos_EcalSeverityLevelAlgo_hh

#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "CondFormats/EcalObjects/interface/EcalChannelStatus.h"

class EcalSeverityLevelAlgo {
        public:
                // give the severity level from the EcalRecHit flags + the DB information stored in EcalChannelStatus
                // Levels of severity:
                // - kGood        --> good channel
                // - kProblematic --> problematic (e.g. noisy)
                // - kRecovered   --> recovered (e.g. an originally dead or saturated)
                // - kTime        --> the channel is out of time (e.g. spike)
                // - kWeird       --> weird (e.g. spike)
                // - kBad         --> bad, not suitable to be used in the reconstruction

                enum EcalSeverityLevel { kGood=0, kProblematic, kRecovered, kTime, kWeird, kBad };

                enum SpikeId { kE1OverE9=0, kSwissCross, kSwissCrossBordersIncluded };

                static int severityLevel( const DetId,
                                          const EcalRecHitCollection &,
                                          const EcalChannelStatus &,
                                          float recHitEtThreshold = 5.,
                                          SpikeId spId = kSwissCross,
                                          float spIdThreshold = 0.95,
                                          float recHitEnergyThresholdForTiming = 2.,
                                          float recHitEnergyThresholdForEE = 15,
                                          float spIdThresholdIEta85 = 0.999
                                          );

                static float spikeFromNeighbours( const DetId id,
                                                  const EcalRecHitCollection &,
                                                  float recHitEtThreshold,
                                                  SpikeId spId
                                                  );

                static float E1OverE9( const DetId id, const EcalRecHitCollection &, float recHitEtThreshold = 0. );

                static float E2overE9( const DetId id, const EcalRecHitCollection &, float recHitEtThreshold = 10.0 , 
                                       float recHitEtThreshold2 = 1.0 , bool avoidIeta85=false, bool KillSecondHit=true);


                static float swissCross( const DetId id, const EcalRecHitCollection &, float recHitEtThreshold = 0. , bool avoidIeta85=true);

                static bool spikeFromTiming( const EcalRecHit &, float recHitEnergyThreshold );

        private:

                static int severityLevel( uint32_t rhFlag, uint16_t dbStatus );
                static int severityLevel( const EcalRecHit &, const EcalChannelStatus & );

                static uint16_t retrieveDBStatus( const DetId , const EcalChannelStatus &chStatus );

                static float recHitE( const DetId id, const EcalRecHitCollection &recHits );
                static float recHitApproxEt( const DetId id, const EcalRecHitCollection &recHits );
                static float recHitE( const DetId id, const EcalRecHitCollection & recHits, int dEta, int dPhi );
};

#endif