#include <SignCaloSpecificAlgo.h>

Public Types
typedef math::XYZTLorentzVector	LorentzVector
typedef math::XYZPoint	Point
typedef std::vector< const reco::Candidate * >	TowerCollection
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_

Detailed Description

Definition at line 35 of file SignCaloSpecificAlgo.h.

Member Typedef Documentation

typedef math::XYZTLorentzVector SignCaloSpecificAlgo::LorentzVector

Definition at line 43 of file SignCaloSpecificAlgo.h.

typedef math::XYZPoint SignCaloSpecificAlgo::Point

Definition at line 44 of file SignCaloSpecificAlgo.h.

typedef std::vector<const reco::Candidate*> SignCaloSpecificAlgo::TowerCollection

Definition at line 45 of file SignCaloSpecificAlgo.h.

Constructor & Destructor Documentation

SignCaloSpecificAlgo::SignCaloSpecificAlgo ( )

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

{
}

Member Function Documentation

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

Definition at line 168 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().

{ 

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

References significance_.

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

{return significance_;}

TMatrixD SignCaloSpecificAlgo::getSignificanceMatrix ( void ) const [inline]

Definition at line 49 of file SignCaloSpecificAlgo.h.

References matrix_.

Referenced by cms::METProducer::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 67 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(), 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 47 of file SignCaloSpecificAlgo.h.

References makeMuonMisalignmentScenario::matrix, and matrix_.

{matrix_=matrix;}

void SignCaloSpecificAlgo::usePreviousSignif ( const std::vector< double > & values )

Definition at line 52 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;
}

Member Data Documentation

TMatrixD SignCaloSpecificAlgo::matrix_ [private]

Definition at line 58 of file SignCaloSpecificAlgo.h.

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

double SignCaloSpecificAlgo::significance_ [private]

Definition at line 57 of file SignCaloSpecificAlgo.h.

Referenced by calculateBaseCaloMET(), and getSignificance().

SignCaloSpecificAlgo Class Reference

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation