Helper for MultiVar2D. More...

Public Member Functions
	MultiVar2DFunc (const _MultiVar< T > &a, const Binning binsa, const _MultiVar< T > &b, const Binning binsb)

std::vector< double >	operator() (const T *sr) const

Protected Attributes
const _MultiVar< T >	fA

const Binning	fBinsA

const _MultiVar< T >	fB

const Binning	fBinsB

Detailed Description

template<class T>
class ana::MultiVar2DFunc< T >

Helper for MultiVar2D.

Definition at line 19 of file MultiVar.cxx.

Constructor & Destructor Documentation

template<class T >

ana::MultiVar2DFunc< T >::MultiVar2DFunc	(	const _MultiVar< T > &	a,
		const Binning	binsa,
		const _MultiVar< T > &	b,
		const Binning	binsb
	)

inline

Definition at line 22 of file MultiVar.cxx.

       : fA(a), fBinsA(binsa),
         fB(b), fBinsB(binsb)
     {
     }

Member Function Documentation

template<class T >

std::vector<double> ana::MultiVar2DFunc< T >::operator() ( const T * sr ) const

inline

Definition at line 29 of file MultiVar.cxx.

     {
       std::vector<double> returnVec;
 
       const std::vector<double> vaVec = fA(sr);
       const std::vector<double> vbVec = fB(sr);
 
       if(vaVec.size() != vbVec.size())
         {
           std::cout << "MultiVars need to be same size, these two are size " 
                     << vaVec.size() << " " << vbVec.size() << "respectively." << std::endl;
           std::abort();
         }
       
       for(unsigned n = 0; n < vaVec.size(); ++n){
         const double va = vaVec.at(n);
         const double vb = vbVec.at(n);
         // Since there are no overflow/underflow bins, check the range
         if(va < fBinsA.Min() || vb < fBinsB.Min()){ returnVec.push_back(-1); continue;}
         if(va > fBinsA.Max() || vb > fBinsB.Max()){ returnVec.push_back(fBinsA.NBins() * fBinsB.NBins()); continue;}
 
         // FindBin uses root convention, first bin is bin 1, bin 0 is underflow
         const int ia = fBinsA.FindBin(va) - 1;
         const int ib = fBinsB.FindBin(vb) - 1;
 
         const int i = ia*fBinsB.NBins()+ib;
 
         returnVec.push_back(i+.5);
       }
       return returnVec;
     }