SignCaloSpecificAlgo Class Reference

#include <RecoMET/METAlgorithms/interface/SignCaloSpecificAlgo.h>

Public Member Functions
void	calculateBaseCaloMET (edm::Handle< edm::View< reco::Candidate > > towers, const CommonMETData &met, const metsig::SignAlgoResolutions &resolutions, bool noHF, double globalthreshold)
double	getSignificance ()
reco::METCovMatrix	getSignificanceMatrix () const
	SignCaloSpecificAlgo ()
void	usePreviousSignif (const std::vector< double > &values)
void	usePreviousSignif (const reco::METCovMatrix &matrix)
	~SignCaloSpecificAlgo ()

Private Member Functions
std::vector< metsig::SigInputObj >	makeVectorOutOfCaloTowers (edm::Handle< edm::View< reco::Candidate > > towers, const metsig::SignAlgoResolutions &resolutions, bool noHF, double globalthreshold)

Private Attributes
reco::METCovMatrix	matrix_
double	significance_

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 34 of file SignCaloSpecificAlgo.h.

Constructor & Destructor Documentation

SignCaloSpecificAlgo::SignCaloSpecificAlgo

(

)

Definition at line 25 of file SignCaloSpecificAlgo.cc.

References matrix_.

                                           : significance_(0.) {
  matrix_(0, 0) = matrix_(1, 0) = matrix_(0, 1) = matrix_(1, 1) = 0.;
}

SignCaloSpecificAlgo::~SignCaloSpecificAlgo

(

)

Definition at line 28 of file SignCaloSpecificAlgo.cc.

{}

Member Function Documentation

void SignCaloSpecificAlgo::calculateBaseCaloMET	(	edm::Handle< edm::View< reco::Candidate > >	towers,
		const CommonMETData &	met,
		const metsig::SignAlgoResolutions &	resolutions,
		bool	noHF,
		double	globalthreshold
	)

Definition at line 152 of file SignCaloSpecificAlgo.cc.

References metsig::significanceAlgo::addObjects(), metsig::significanceAlgo::addSignifMatrix(), metsig::significanceAlgo::getSignifMatrix(), makeVectorOutOfCaloTowers(), matrix_, metsig::significanceAlgo::significance(), significance_, and HLT_FULL_cff::towers.

                                                                        {
  //retreive calo tower information from candidates
  //start with the first element of the candidate list

  // use this container to calculate the significance. SigInputObj are objects that contain both directional and uncertainty information and are used as input to the significance calculation

  std::vector<metsig::SigInputObj> signInputVec = makeVectorOutOfCaloTowers(towers, resolutions, noHF, globalThreshold);

  // now run the significance algorithm.

  double sign_calo_met_total = 0;
  double sign_calo_met_phi = 0;
  double sign_calo_met_set = 0;
  metsig::significanceAlgo signifalgo;
  // check the caloMET, if significance was already run continue with the matrix that is stored..
  signifalgo.addSignifMatrix(matrix_);
  signifalgo.addObjects(signInputVec);
  matrix_ = signifalgo.getSignifMatrix();
  significance_ = signifalgo.significance(sign_calo_met_total, sign_calo_met_phi, sign_calo_met_set);
  // cleanup everything:
  signInputVec.clear();
  // and return
}

double SignCaloSpecificAlgo::getSignificance ( )

inline

Definition at line 41 of file SignCaloSpecificAlgo.h.

References significance_.

{ return significance_; }

reco::METCovMatrix SignCaloSpecificAlgo::getSignificanceMatrix ( void  ) const

inline

Definition at line 42 of file SignCaloSpecificAlgo.h.

References matrix_.

{ return matrix_; }

std::vector< metsig::SigInputObj > SignCaloSpecificAlgo::makeVectorOutOfCaloTowers ( edm::Handle< edm::View< reco::Candidate > >  towers, const metsig::SignAlgoResolutions &  resolutions, bool  noHF, double  globalthreshold  )

private

Definition at line 43 of file SignCaloSpecificAlgo.cc.

References metsig::caloEB, metsig::caloEE, metsig::caloHB, metsig::caloHE, metsig::caloHF, metsig::caloHO, CaloTower::constituent(), CaloTower::constituentsSize(), DetId::Ecal, EcalBarrel, EcalEndcap, CaloTower::emEt(), metsig::ET, CaloTower::et(), reco::LeafCandidate::eta(), CaloTower::hadEt(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, CaloTower::outerEt(), metsig::PHI, reco::LeafCandidate::phi(), AlCaHLTBitMon_QueryRunRegistry::string, HcalDetId::subdet(), groupFilesInBlocks::temp, and HLT_FULL_cff::towers.

Referenced by calculateBaseCaloMET().

                            {
  edm::View<Candidate>::const_iterator towerCand = towers->begin();
  std::vector<metsig::SigInputObj> signInputVec;
  //iterate over all CaloTowers and record information
  for (; towerCand != towers->end(); towerCand++) {
    const Candidate* candidate = &(*towerCand);
    if (candidate) {
      const CaloTower* calotower = dynamic_cast<const CaloTower*>(candidate);
      if (calotower) {
        double sign_tower_et = calotower->et();
        if (sign_tower_et < globalThreshold)
          continue;
        bool wasused = false;
        double sign_tower_phi = calotower->phi();
        double sign_tower_sigma_et = 0;
        double sign_tower_sigma_phi = 0;
        std::string sign_tower_type = "";

        bool hadIsDone = false;
        bool emIsDone = false;
        int cell = calotower->constituentsSize();

        while (--cell >= 0 && (!hadIsDone || !emIsDone)) {
          DetId id = calotower->constituent(cell);
          if (!hadIsDone && id.det() == DetId::Hcal) {
            HcalSubdetector subdet = HcalDetId(id).subdet();
            if (subdet == HcalBarrel) {
              sign_tower_type = "hadcalotower";
              sign_tower_et = calotower->hadEt();
              sign_tower_sigma_et =
                  resolutions.eval(metsig::caloHB, metsig::ET, sign_tower_et, calotower->phi(), calotower->eta());
              sign_tower_sigma_phi =
                  resolutions.eval(metsig::caloHB, metsig::PHI, sign_tower_et, calotower->phi(), calotower->eta());
            } else if (subdet == HcalOuter) {
              sign_tower_type = "hadcalotower";
              sign_tower_et = calotower->outerEt();
              sign_tower_sigma_et =
                  resolutions.eval(metsig::caloHO, metsig::ET, sign_tower_et, calotower->phi(), calotower->eta());
              sign_tower_sigma_phi =
                  resolutions.eval(metsig::caloHO, metsig::PHI, sign_tower_et, calotower->phi(), calotower->eta());
            } else if (subdet == HcalEndcap) {
              sign_tower_type = "hadcalotower";
              sign_tower_et = calotower->hadEt();
              sign_tower_sigma_et =
                  resolutions.eval(metsig::caloHE, metsig::ET, sign_tower_et, calotower->phi(), calotower->eta());
              sign_tower_sigma_phi =
                  resolutions.eval(metsig::caloHE, metsig::PHI, sign_tower_et, calotower->phi(), calotower->eta());
            } else if (subdet == HcalForward) {
              sign_tower_type = "hadcalotower";
              sign_tower_et = calotower->et();
              sign_tower_sigma_et =
                  resolutions.eval(metsig::caloHF, metsig::ET, sign_tower_et, calotower->phi(), calotower->eta());
              sign_tower_sigma_phi =
                  resolutions.eval(metsig::caloHF, metsig::PHI, sign_tower_et, calotower->phi(), calotower->eta());
            } else {
              edm::LogWarning("SignCaloSpecificAlgo")
                  << " HCAL tower cell not assigned to an HCAL subdetector!!!" << std::endl;
            }
            // and book!
            metsig::SigInputObj temp(
                sign_tower_type, sign_tower_et, sign_tower_phi, sign_tower_sigma_et, sign_tower_sigma_phi);
            if (!noHF || subdet != HcalForward)
              signInputVec.push_back(temp);

            wasused = true;
            hadIsDone = true;
          } else if (!emIsDone && id.det() == DetId::Ecal) {
            EcalSubdetector subdet = EcalSubdetector(id.subdetId());

            if (subdet == EcalBarrel) {
              sign_tower_type = "emcalotower";
              sign_tower_et = calotower->emEt();
              sign_tower_sigma_et =
                  resolutions.eval(metsig::caloEB, metsig::ET, sign_tower_et, calotower->phi(), calotower->eta());
              sign_tower_sigma_phi =
                  resolutions.eval(metsig::caloEB, metsig::PHI, sign_tower_et, calotower->phi(), calotower->eta());
            } else if (subdet == EcalEndcap) {
              sign_tower_type = "emcalotower";
              sign_tower_et = calotower->emEt();
              sign_tower_sigma_et =
                  resolutions.eval(metsig::caloEE, metsig::ET, sign_tower_et, calotower->phi(), calotower->eta());
              sign_tower_sigma_phi =
                  resolutions.eval(metsig::caloEE, metsig::PHI, sign_tower_et, calotower->phi(), calotower->eta());

            } else {
              edm::LogWarning("SignCaloSpecificAlgo")
                  << " ECAL tower cell not assigned to an ECAL subdetector!!!" << std::endl;
            }
            metsig::SigInputObj temp(
                sign_tower_type, sign_tower_et, sign_tower_phi, sign_tower_sigma_et, sign_tower_sigma_phi);
            signInputVec.push_back(temp);
            wasused = true;
            emIsDone = true;
          }
        }
        if (wasused == 0)
          edm::LogWarning("SignCaloSpecificAlgo") << "found non-assigned cell, " << std::endl;
      }
    }
  }
  return signInputVec;
}

void SignCaloSpecificAlgo::usePreviousSignif

(

const std::vector< double > &

values

)

Definition at line 30 of file SignCaloSpecificAlgo.cc.

References matrix_.

                                                                            {
  if (values.size() != 4)
    return;
  matrix_(0, 0) = values[0];
  matrix_(0, 1) = values[1];
  matrix_(1, 0) = values[2];
  matrix_(1, 1) = values[3];
  return;
}

void SignCaloSpecificAlgo::usePreviousSignif ( const reco::METCovMatrix &  matrix )

inline

Definition at line 40 of file SignCaloSpecificAlgo.h.

References makeMuonMisalignmentScenario::matrix, and matrix_.

{ matrix_ = matrix; }

Member Data Documentation

reco::METCovMatrix SignCaloSpecificAlgo::matrix_

private

Definition at line 57 of file SignCaloSpecificAlgo.h.

Referenced by calculateBaseCaloMET(), getSignificanceMatrix(), SignCaloSpecificAlgo(), and usePreviousSignif().

double SignCaloSpecificAlgo::significance_

private

Definition at line 56 of file SignCaloSpecificAlgo.h.

Referenced by calculateBaseCaloMET(), and getSignificance().