|
|
|
#include <RecoMET/METAlgorithms/interface/SignCaloSpecificAlgo.h>
Public Member Functions | |
| void | calculateBaseCaloMET (edm::Handle< edm::View< reco::Candidate > > towers, CommonMETData met, const metsig::SignAlgoResolutions &resolutions, bool noHF, double globalthreshold) |
| double | getSignificance () |
| TMatrixD | getSignificanceMatrix () const |
| SignCaloSpecificAlgo () | |
| void | usePreviousSignif (const TMatrixD &matrix) |
| void | usePreviousSignif (const std::vector< double > &values) |
| ~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_ |
Description: <one line="" class="" summary>="">
Implementation: <Notes on="" implementation>="">
Definition at line 35 of file SignCaloSpecificAlgo.h.
| SignCaloSpecificAlgo::SignCaloSpecificAlgo | ( | ) |
Definition at line 26 of file SignCaloSpecificAlgo.cc.
References matrix_.
| SignCaloSpecificAlgo::~SignCaloSpecificAlgo | ( | ) |
Definition at line 32 of file SignCaloSpecificAlgo.cc.
{
}
| void SignCaloSpecificAlgo::calculateBaseCaloMET | ( | edm::Handle< edm::View< reco::Candidate > > | towers, |
| 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(), and significance_.
Referenced by cms::METProducer::produce_CaloMET().
{
//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 44 of file SignCaloSpecificAlgo.h.
References significance_.
Referenced by cms::METProducer::produce_CaloMET().
{return significance_;}
| TMatrixD SignCaloSpecificAlgo::getSignificanceMatrix | ( | void | ) | const [inline] |
Definition at line 45 of file SignCaloSpecificAlgo.h.
References matrix_.
Referenced by cms::METProducer::produce_CaloMET().
{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 51 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 TMatrixD & | matrix | ) | [inline] |
Definition at line 43 of file SignCaloSpecificAlgo.h.
References makeMuonMisalignmentScenario::matrix, and matrix_.
| void SignCaloSpecificAlgo::usePreviousSignif | ( | const std::vector< double > & | values | ) |
TMatrixD SignCaloSpecificAlgo::matrix_ [private] |
Definition at line 54 of file SignCaloSpecificAlgo.h.
Referenced by calculateBaseCaloMET(), getSignificanceMatrix(), SignCaloSpecificAlgo(), and usePreviousSignif().
double SignCaloSpecificAlgo::significance_ [private] |
Definition at line 53 of file SignCaloSpecificAlgo.h.
Referenced by calculateBaseCaloMET(), and getSignificance().
1.7.3