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 ()
TMatrixD	getSignificanceMatrix () const
	SignCaloSpecificAlgo ()
void	usePreviousSignif (const std::vector< double > &values)
void	usePreviousSignif (const TMatrixD &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
TMatrixD	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_(2,2)
{
  matrix_(0,0)=matrix_(1,0)=matrix_(0,1)=matrix_(1,1)=0.;
}

SignCaloSpecificAlgo::~SignCaloSpecificAlgo

(

)

Definition at line 31 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 151 of file SignCaloSpecificAlgo.cc.

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

Referenced by cms::CaloMETProducer::produce().

{ 

  //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 43 of file SignCaloSpecificAlgo.h.

References significance_.

Referenced by cms::CaloMETProducer::produce().

{return significance_;}

TMatrixD SignCaloSpecificAlgo::getSignificanceMatrix ( void  ) const

inline

Definition at line 44 of file SignCaloSpecificAlgo.h.

References matrix_.

Referenced by cms::CaloMETProducer::produce().

{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 50 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(), edm::View< T >::end(), 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(), and groupFilesInBlocks::temp.

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=1;
        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=1;
        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 35 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 TMatrixD &  matrix )

inline

Definition at line 42 of file SignCaloSpecificAlgo.h.

References makeMuonMisalignmentScenario::matrix, and matrix_.

{ matrix_ = matrix; }

Member Data Documentation

TMatrixD SignCaloSpecificAlgo::matrix_

private

Definition at line 53 of file SignCaloSpecificAlgo.h.

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

double SignCaloSpecificAlgo::significance_

private

Definition at line 52 of file SignCaloSpecificAlgo.h.

Referenced by calculateBaseCaloMET(), and getSignificance().