AnomalousECALVariables Class Reference

#include <AnomalousECALVariables.h>

Public Member Functions
	AnomalousECALVariables ()
	AnomalousECALVariables (const std::vector< BoundaryInformation > &p_enNeighboursGap_EB, const std::vector< BoundaryInformation > &p_enNeighboursGap_EE, const std::vector< BoundaryInformation > &p_boundaryInfoDeadCells_EB, const std::vector< BoundaryInformation > &p_boundaryInfoDeadCells_EE)
bool	isDeadEcalCluster (double maxBoundaryEnergy=10, const std::vector< int > &limitDeadCellToChannelStatusEB=std::vector< int >(), const std::vector< int > &limitDeadCellToChannelStatusEE=std::vector< int >()) const
bool	isGapEcalCluster (double maxGapEnergyEB=10, double maxGapEnergyEE=10) const
	~AnomalousECALVariables ()

Public Attributes
std::vector< BoundaryInformation >	v_boundaryInfoDeadCells_EB
std::vector< BoundaryInformation >	v_boundaryInfoDeadCells_EE
std::vector< BoundaryInformation >	v_enNeighboursGap_EB
std::vector< BoundaryInformation >	v_enNeighboursGap_EE

Detailed Description

Definition at line 20 of file AnomalousECALVariables.h.

Constructor & Destructor Documentation

AnomalousECALVariables() [1/2]

AnomalousECALVariables::AnomalousECALVariables ( )

inline

Definition at line 22 of file AnomalousECALVariables.h.

References v_boundaryInfoDeadCells_EB, v_boundaryInfoDeadCells_EE, v_enNeighboursGap_EB, and v_enNeighboursGap_EE.

                           {
    //energy next to ECAL Gap
    v_enNeighboursGap_EB.reserve(50);
    v_enNeighboursGap_EE.reserve(50);
    v_enNeighboursGap_EB.clear();
    v_enNeighboursGap_EE.clear();

    //energy around dead cells
    v_boundaryInfoDeadCells_EB = std::vector<BoundaryInformation>();
    v_boundaryInfoDeadCells_EE = std::vector<BoundaryInformation>();
    v_boundaryInfoDeadCells_EB.reserve(50);
    v_boundaryInfoDeadCells_EE.reserve(50);
    v_boundaryInfoDeadCells_EB.clear();
    v_boundaryInfoDeadCells_EE.clear();
  };

AnomalousECALVariables() [2/2]

AnomalousECALVariables::AnomalousECALVariables ( const std::vector< BoundaryInformation > &  p_enNeighboursGap_EB, const std::vector< BoundaryInformation > &  p_enNeighboursGap_EE, const std::vector< BoundaryInformation > &  p_boundaryInfoDeadCells_EB, const std::vector< BoundaryInformation > &  p_boundaryInfoDeadCells_EE  )

inline

Definition at line 38 of file AnomalousECALVariables.h.

References v_boundaryInfoDeadCells_EB, v_boundaryInfoDeadCells_EE, v_enNeighboursGap_EB, and v_enNeighboursGap_EE.

                                                                                           {
    v_boundaryInfoDeadCells_EB = std::vector<BoundaryInformation>();
    v_boundaryInfoDeadCells_EE = std::vector<BoundaryInformation>();
    v_boundaryInfoDeadCells_EB.reserve(50);
    v_boundaryInfoDeadCells_EE.reserve(50);
    v_boundaryInfoDeadCells_EB.clear();
    v_boundaryInfoDeadCells_EE.clear();
    v_boundaryInfoDeadCells_EB = p_boundaryInfoDeadCells_EB;
    v_boundaryInfoDeadCells_EE = p_boundaryInfoDeadCells_EE;

    v_enNeighboursGap_EB = p_enNeighboursGap_EB;
    v_enNeighboursGap_EE = p_enNeighboursGap_EE;
  };

~AnomalousECALVariables()

AnomalousECALVariables::~AnomalousECALVariables ( )

inline

Definition at line 55 of file AnomalousECALVariables.h.

References v_boundaryInfoDeadCells_EB, v_boundaryInfoDeadCells_EE, v_enNeighboursGap_EB, and v_enNeighboursGap_EE.

                            {
    //cout<<"destructor AnomalousECAL"<<endl;
    v_enNeighboursGap_EB.clear();
    v_enNeighboursGap_EE.clear();
    v_boundaryInfoDeadCells_EB.clear();
    v_boundaryInfoDeadCells_EE.clear();
  };

Member Function Documentation

isDeadEcalCluster()

bool AnomalousECALVariables::isDeadEcalCluster ( double  maxBoundaryEnergy = 10, const std::vector< int > &  limitDeadCellToChannelStatusEB = std::vector<int>(), const std::vector< int > &  limitDeadCellToChannelStatusEE = std::vector<int>()  ) const

inline

Definition at line 67 of file AnomalousECALVariables.h.

References funct::abs(), BoundaryInformation::boundaryEnergy, BoundaryInformation::boundaryET, BoundaryInformation::channelStatus, callgraph::cs, mps_fire::i, createfilelist::int, EcalDeadCellBoundaryEnergyFilter_cfi::limitDeadCellToChannelStatusEB, EcalDeadCellBoundaryEnergyFilter_cfi::limitDeadCellToChannelStatusEE, mps_update::status, v_boundaryInfoDeadCells_EB, and v_boundaryInfoDeadCells_EE.

                                                                                                        {
    float highestEnergyDepositAroundDeadCell = 0;

    for (int i = 0; i < (int)v_boundaryInfoDeadCells_EB.size(); ++i) {
      BoundaryInformation bInfo = v_boundaryInfoDeadCells_EB[i];

      //check if channel limitation rejectsbInfo
      bool passChannelLimitation = false;
      std::vector<int> status = bInfo.channelStatus;

      for (int cs = 0; cs < (int)limitDeadCellToChannelStatusEB.size(); ++cs) {
        int channelAllowed = limitDeadCellToChannelStatusEB[cs];

        for (std::vector<int>::iterator st_it = status.begin(); st_it != status.end(); ++st_it) {
          if (channelAllowed == *st_it || (channelAllowed < 0 && abs(channelAllowed) <= *st_it)) {
            passChannelLimitation = true;
            break;
          }
        }
      }

      if (passChannelLimitation || limitDeadCellToChannelStatusEB.empty()) {
        if (bInfo.boundaryEnergy > highestEnergyDepositAroundDeadCell) {
          highestEnergyDepositAroundDeadCell = bInfo.boundaryET;
          //highestEnergyDepositAroundDeadCell = bInfo.boundaryEnergy;
        }
      }
    }

    for (int i = 0; i < (int)v_boundaryInfoDeadCells_EE.size(); ++i) {
      BoundaryInformation bInfo = v_boundaryInfoDeadCells_EE[i];

      //check if channel limitation rejectsbInfo
      bool passChannelLimitation = false;
      std::vector<int> status = bInfo.channelStatus;

      for (int cs = 0; cs < (int)limitDeadCellToChannelStatusEE.size(); ++cs) {
        int channelAllowed = limitDeadCellToChannelStatusEE[cs];

        for (std::vector<int>::iterator st_it = status.begin(); st_it != status.end(); ++st_it) {
          if (channelAllowed == *st_it || (channelAllowed < 0 && abs(channelAllowed) <= *st_it)) {
            passChannelLimitation = true;
            break;
          }
        }
      }

      if (passChannelLimitation || limitDeadCellToChannelStatusEE.empty()) {
        if (bInfo.boundaryEnergy > highestEnergyDepositAroundDeadCell) {
          highestEnergyDepositAroundDeadCell = bInfo.boundaryET;
          //highestEnergyDepositAroundDeadCell = bInfo.boundaryEnergy;
        }
      }
    }

    if (highestEnergyDepositAroundDeadCell > maxBoundaryEnergy) {
      //            cout << "<<<<<<<<<< List of EB  Boundary objects <<<<<<<<<<" << endl;
      //            for (int i = 0; i < (int) v_boundaryInfoDeadCells_EB.size(); ++i) {
      //               BoundaryInformation bInfo = v_boundaryInfoDeadCells_EB[i];
      //               cout << "no of neighbouring RecHits:" << bInfo.recHits.size() << endl;
      //               cout << "no of neighbouring DetIds:" << bInfo.detIds.size() << endl;
      //               cout << "boundary energy:" << bInfo.boundaryEnergy << endl;
      //               cout << "Channel stati: ";
      //               for (std::vector<int>::iterator it = bInfo.channelStatus.begin(); it != bInfo.channelStatus.end(); ++it) {
      //                  cout << *it << " ";
      //               }
      //               cout << endl;
      //            }
      //            cout << "<<<<<<<<<< List of EE  Boundary objects <<<<<<<<<<" << endl;
      //            for (int i = 0; i < (int) v_boundaryInfoDeadCells_EE.size(); ++i) {
      //               BoundaryInformation bInfo = v_boundaryInfoDeadCells_EE[i];
      //               cout << "no of neighbouring RecHits:" << bInfo.recHits.size() << endl;
      //               cout << "no of neighbouring DetIds:" << bInfo.detIds.size() << endl;
      //               cout << "boundary energy:" << bInfo.boundaryEnergy << endl;
      //               cout << "Channel stati: ";
      //               for (std::vector<int>::iterator it = bInfo.channelStatus.begin(); it != bInfo.channelStatus.end(); ++it) {
      //                  cout << *it << " ";
      //               }
      //               cout << endl;
      //            }
      return true;
    } else
      return false;
  }

isGapEcalCluster()

bool AnomalousECALVariables::isGapEcalCluster ( double  maxGapEnergyEB = 10, double  maxGapEnergyEE = 10  ) const

inline

Definition at line 154 of file AnomalousECALVariables.h.

References BoundaryInformation::boundaryEnergy, BoundaryInformation::boundaryET, mps_fire::i, createfilelist::int, v_enNeighboursGap_EB, and v_enNeighboursGap_EE.

                                                                                      {
    float highestEnergyDepositAlongGapEB = 0;

    for (int i = 0; i < (int)v_enNeighboursGap_EB.size(); ++i) {
      BoundaryInformation gapInfo = v_enNeighboursGap_EB[i];

      if (gapInfo.boundaryEnergy > highestEnergyDepositAlongGapEB) {
        highestEnergyDepositAlongGapEB = gapInfo.boundaryET;
        //highestEnergyDepositAlongGapEB = gapInfo.boundaryEnergy;
      }
    }

    float highestEnergyDepositAlongGapEE = 0;

    for (int i = 0; i < (int)v_enNeighboursGap_EE.size(); ++i) {
      BoundaryInformation gapInfo = v_enNeighboursGap_EE[i];

      if (gapInfo.boundaryEnergy > highestEnergyDepositAlongGapEE) {
        highestEnergyDepositAlongGapEE = gapInfo.boundaryET;
        //highestEnergyDepositAlongGapEE = gapInfo.boundaryEnergy;
      }
    }

    if (highestEnergyDepositAlongGapEB > maxGapEnergyEB || highestEnergyDepositAlongGapEE > maxGapEnergyEE) {
      //            cout << "<<<<<<<<<< List of EB Gap objects <<<<<<<<<<" << endl;
      //            for (int i = 0; i < (int) v_enNeighboursGap_EB.size(); ++i) {
      //               BoundaryInformation gapInfo = v_enNeighboursGap_EB[i];
      //               cout << "no of neighbouring RecHits:" << gapInfo.recHits.size() << endl;
      //               cout << "no of neighbouring DetIds:" << gapInfo.detIds.size() << endl;
      //               cout << "gap energy:" << gapInfo.boundaryEnergy << endl;
      //            }
      //            cout << "<<<<<<<<<< List of EE Gap objects <<<<<<<<<<" << endl;
      //            for (int i = 0; i < (int) v_enNeighboursGap_EE.size(); ++i) {
      //               BoundaryInformation gapInfo = v_enNeighboursGap_EE[i];
      //               cout << "no of neighbouring RecHits:" << gapInfo.recHits.size() << endl;
      //               cout << "no of neighbouring DetIds:" << gapInfo.detIds.size() << endl;
      //               cout << "gap energy:" << gapInfo.boundaryEnergy << endl;
      //            }
      return true;
    } else
      return false;
  }

Member Data Documentation

v_boundaryInfoDeadCells_EB

std::vector<BoundaryInformation> AnomalousECALVariables::v_boundaryInfoDeadCells_EB

Definition at line 200 of file AnomalousECALVariables.h.

Referenced by AnomalousECALVariables(), isDeadEcalCluster(), and ~AnomalousECALVariables().

v_boundaryInfoDeadCells_EE

std::vector<BoundaryInformation> AnomalousECALVariables::v_boundaryInfoDeadCells_EE

Definition at line 201 of file AnomalousECALVariables.h.

Referenced by AnomalousECALVariables(), isDeadEcalCluster(), and ~AnomalousECALVariables().

v_enNeighboursGap_EB

std::vector<BoundaryInformation> AnomalousECALVariables::v_enNeighboursGap_EB

Definition at line 197 of file AnomalousECALVariables.h.

Referenced by AnomalousECALVariables(), isGapEcalCluster(), and ~AnomalousECALVariables().

v_enNeighboursGap_EE

std::vector<BoundaryInformation> AnomalousECALVariables::v_enNeighboursGap_EE

Definition at line 198 of file AnomalousECALVariables.h.

Referenced by AnomalousECALVariables(), isGapEcalCluster(), and ~AnomalousECALVariables().