#include <significanceAlgo.h>

Public Member Functions
const void	addObjects (const std::vector< metsig::SigInputObj > &EventVec)

const void	addSignifMatrix (const reco::METCovMatrix &input)

reco::METCovMatrix	getSignifMatrix () const

const void	setSignifMatrix (const reco::METCovMatrix &input, const double &met_r, const double &met_phi, const double &met_set)
	reset the signficance matrix (this is the most likely case), so that the vector sum can be continued More...

const double	significance (double &met_r, double &met_phi, double &met_set)

	significanceAlgo ()

const void	subtractObjects (const std::vector< metsig::SigInputObj > &EventVec)

	~significanceAlgo ()

Private Attributes
double	set_worker_

reco::METCovMatrix	signifmatrix_

double	xmet_

double	ymet_

Detailed Description

Definition at line 86 of file significanceAlgo.h.

Constructor & Destructor Documentation

metsig::significanceAlgo::significanceAlgo ( )

Definition at line 27 of file significanceAlgo.cc.

References signifmatrix_.

                                         :
   //  eventVec_(0),
   set_worker_(0),
   xmet_(0),
   ymet_(0)
 
 {
   //  std::cout << "in constructor ! " << std::endl;
   signifmatrix_(0,0)=signifmatrix_(1,0)=signifmatrix_(0,1)=signifmatrix_(1,1)=0;
 }

metsig::significanceAlgo::~significanceAlgo ( )

Definition at line 60 of file significanceAlgo.cc.

60 {

61 }

Member Function Documentation

const void metsig::significanceAlgo::addObjects ( const std::vector< metsig::SigInputObj > & EventVec )

Definition at line 116 of file significanceAlgo.cc.

References funct::cos(), hgcalPlots::obj, set_worker_, signifmatrix_, funct::sin(), xmet_, and ymet_.

Referenced by metsig::SignPFSpecificAlgo::addPFCandidate(), metsig::SignPFSpecificAlgo::addPFJets(), SignCaloSpecificAlgo::calculateBaseCaloMET(), computePFMEtSignificance(), PFMEtSignInterfaceBase::operator()(), and pat::PATMHTProducer::produce().

 { 
   reco::METCovMatrix v_tot = signifmatrix_;
   //--- Loop over physics objects in the event ---//
   //  for(unsigned int objnum=1; objnum < EventVec.size(); objnum++ ) {
   for(std::vector<SigInputObj>::const_iterator obj = eventVec.begin(); obj!= eventVec.end(); ++obj){
     double et_tmp     = obj->get_energy();
     double phi_tmp    = obj->get_phi();
     double sigma_et   = obj->get_sigma_e();
     double sigma_tan  = obj->get_sigma_tan();
 
     double cosphi=cos(phi_tmp);
     double sinphi=sin(phi_tmp);
     double xval = et_tmp * cosphi;
     double yval = et_tmp * sinphi;
     xmet_ -= xval;
     ymet_ -= yval;
     set_worker_ += et_tmp;
 
     double sigma0_2=sigma_et*sigma_et;
     double sigma1_2=sigma_tan*sigma_tan;
 
     v_tot(0,0)+= sigma0_2*cosphi*cosphi + sigma1_2*sinphi*sinphi;
     v_tot(0,1)+= cosphi*sinphi*(sigma0_2 - sigma1_2);
     v_tot(1,0)+= cosphi*sinphi*(sigma0_2 - sigma1_2);
     v_tot(1,1)+= sigma1_2*cosphi*cosphi + sigma0_2*sinphi*sinphi;
     
   }
   signifmatrix_=v_tot;
 }

const void metsig::significanceAlgo::addSignifMatrix ( const reco::METCovMatrix & input )

Definition at line 41 of file significanceAlgo.cc.

References input, and signifmatrix_.

Referenced by SignCaloSpecificAlgo::calculateBaseCaloMET().

                                                                     {
   signifmatrix_+=input;
   return;
 }

reco::METCovMatrix metsig::significanceAlgo::getSignifMatrix ( ) const

inline

Definition at line 96 of file significanceAlgo.h.

References signifmatrix_.

Referenced by SignCaloSpecificAlgo::calculateBaseCaloMET(), computePFMEtSignificance(), metsig::SignPFSpecificAlgo::getSignifMatrix(), and PFMEtSignInterfaceBase::operator()().

96 {return signifmatrix_;}

metsig::significanceAlgo::signifmatrix_

reco::METCovMatrix signifmatrix_

Definition: significanceAlgo.h:102

const void metsig::significanceAlgo::setSignifMatrix	(	const reco::METCovMatrix &	input,
		const double &	met_r,
		const double &	met_phi,
		const double &	met_set
	)

reset the signficance matrix (this is the most likely case), so that the vector sum can be continued

Definition at line 49 of file significanceAlgo.cc.

References funct::cos(), input, set_worker_, signifmatrix_, funct::sin(), xmet_, and ymet_.

                                                                                                                                       {
   signifmatrix_=input;
   set_worker_=met_set;
   xmet_=met_r*cos(met_phi);
   ymet_=met_r*sin(met_phi);
   
   return;
 }

const double metsig::significanceAlgo::significance	(	double &	met_r,
		double &	met_phi,
		double &	met_set
	)

Definition at line 149 of file significanceAlgo.cc.

References set_worker_, signifmatrix_, mathSSE::sqrt(), xmet_, and ymet_.

Referenced by SignCaloSpecificAlgo::calculateBaseCaloMET(), and pat::PATMHTProducer::produce().

 {
   
   if(signifmatrix_(0,0)==0 && signifmatrix_(1,1)==0 && signifmatrix_(1,0)==0 && signifmatrix_(0,1)==0){
     //edm::LogWarning("SignCaloSpecificAlgo") << "Event Vector is empty!  Return significance -1";
     return(-1);
   } 
 
   //--- Temporary variables ---//
  
   reco::METCovMatrix v_tot;
   ROOT::Math::SVector<double, 2> metvec;
 
  //--- Initialize sum of rotated covariance matrices ---//  
   v_tot=signifmatrix_;
   //  std::cout << "INPUT:\n"<< v_tot(0,0) << "," << v_tot(0,1) << "\n" << v_tot(1,0) << "," << v_tot(1,1) << std::endl;
 
 
 
   //--- Calculate magnitude and angle of MET, store in returned variables ---//
   met_r = sqrt(xmet_*xmet_ + ymet_*ymet_);
   met_set = set_worker_;
   met_phi= TMath::ATan2(ymet_, xmet_);
 
   // one other option: if particles cancel there could be small numbers.
   // this check fixes this, added by F.Blekman
   double det=0;
   v_tot.Det(det);
   if(fabs(det)<0.000001)
     return -1;
 
 
   // save matrix into object:
   //  std::cout << "SAVED:\n"<< v_tot(0,0) << "," << v_tot(0,1) << "\n" << v_tot(1,0) << "," << v_tot(1,1) << std::endl;
   v_tot.Invert();
   //  std::cout << "INVERTED:\n"<< v_tot(0,0) << "," << v_tot(0,1) << "\n" << v_tot(1,0) << "," << v_tot(1,1) << std::endl;
   
 
 
   metvec(0) = xmet_; metvec(1) = ymet_;
   double lnSignificance = ROOT::Math::Dot(metvec, (v_tot * metvec) );
 
   //  v_tot.Invert();
   //  std::cout << "INVERTED AGAIN:\n"<< v_tot(0,0) << "," << v_tot(0,1) << "\n" << v_tot(1,0) << "," << v_tot(1,1) << std::endl;
   return lnSignificance;
 }

const void metsig::significanceAlgo::subtractObjects ( const std::vector< metsig::SigInputObj > & EventVec )

Definition at line 82 of file significanceAlgo.cc.

References funct::cos(), hgcalPlots::obj, set_worker_, signifmatrix_, funct::sin(), xmet_, and ymet_.

 { 
   reco::METCovMatrix v_tot = signifmatrix_;
   //--- Loop over physics objects in the event ---//
   //  for(unsigned int objnum=1; objnum < EventVec.size(); objnum++ ) {
   for(std::vector<SigInputObj>::const_iterator obj = eventVec.begin(); obj!= eventVec.end(); ++obj){
     double et_tmp     = obj->get_energy();
     double phi_tmp    = obj->get_phi();
     double sigma_et   = obj->get_sigma_e();
     double sigma_tan  = obj->get_sigma_tan();
 
     double cosphi=cos(phi_tmp);
     double sinphi=sin(phi_tmp);
     double xval = et_tmp * cosphi;
     double yval = et_tmp * sinphi;
     xmet_ += xval;
     ymet_ += yval;
     set_worker_ -= et_tmp;
 
     double sigma0_2=sigma_et*sigma_et;
     double sigma1_2=sigma_tan*sigma_tan;
 
     v_tot(0,0)-= sigma0_2*cosphi*cosphi + sigma1_2*sinphi*sinphi;
     v_tot(0,1)-= cosphi*sinphi*(sigma0_2 - sigma1_2);
     v_tot(1,0)-= cosphi*sinphi*(sigma0_2 - sigma1_2);
     v_tot(1,1)-= sigma1_2*cosphi*cosphi + sigma0_2*sinphi*sinphi;
     
   }
   signifmatrix_=v_tot;
 }