metsig::significanceAlgo Class Reference

#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 Member Functions
void	rotateMatrix (Double_t theta, reco::METCovMatrix &v)

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.

                                         {
}

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(), getGTfromDQMFile::obj, and funct::sin().

Referenced by 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.

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()().

{return signifmatrix_;}

void metsig::significanceAlgo::rotateMatrix ( Double_t  theta, reco::METCovMatrix &  v  )

private

Definition at line 66 of file significanceAlgo.cc.

References funct::cos(), alignCSCRings::r, and funct::sin().

{
  // I suggest not using this to rotate trivial matrices.
  reco::METCovMatrix r;
  reco::METCovMatrix rInv;

  r(0,0) = cos(theta); r(0,1) = sin(theta); r(1,0) = -sin(theta); r(1,1) = cos(theta);
  rInv = r;
  rInv.Invert();
  //-- Rotate v --//
  v = rInv * v * r;
}

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, and funct::sin().

                                                                                                                                      {
  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 mathSSE::sqrt().

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(), getGTfromDQMFile::obj, and funct::sin().

{ 
  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;
}

Member Data Documentation

double metsig::significanceAlgo::set_worker_

private

Definition at line 104 of file significanceAlgo.h.

reco::METCovMatrix metsig::significanceAlgo::signifmatrix_

private

Definition at line 102 of file significanceAlgo.h.

Referenced by getSignifMatrix(), and significanceAlgo().

double metsig::significanceAlgo::xmet_

private

Definition at line 105 of file significanceAlgo.h.

double metsig::significanceAlgo::ymet_

private

Definition at line 106 of file significanceAlgo.h.