covmatrix.h

Go to the documentation of this file.

/*
 * To change this license header, choose License Headers in Project Properties.
 * To change this template file, choose Tools | Templates
 * and open the template in the editor.
 */
 
/*
 * File:   covmatrix.h
 * Author: dglazier
 *
 * Created on April 27, 2017, 10:01 AM
 */
 
#ifndef CLAS12_COVMATRIX_H
#define CLAS12_COVMATRIX_H
 
#include <hipo4/bank.h>
#include <vector>
#include <map>
#include <memory>
 
namespace clas12 {
  
  using CovMatrix =std::vector<std::vector<float > >; 
  using IndMatrix =std::vector<std::vector<int > >; 
 
  class covmatrix : public hipo::bank {
 
 
  public:
 
    //covmatrix() { }
    //vmatrix(const char *bankName, hipo::reader &r) : hipo::bank(bankName,r){
    //
    covmatrix(hipo::schema __schema);
 
  
    //get the pindex = corresponding particle index
    int getPindex(int index) { return getInt(_pindex_order,index);}
    int getPindex() { return getInt(_pindex_order,_index);}
    int getIndex() const {return _index;}
    int getIndex(int pindex);
    void setIndex(int ind) const{_index=ind;}
    void setEntry(int ind) const{_index=ind;}
    void setBankEntry(int ind){_index=ind;}
 
    const  CovMatrix* matrix();
    float getC11()  const noexcept{ return getFloat(_morder[0][0],_index);}
    float getC12()  const noexcept{ return getFloat(_morder[0][1],_index);}
    float getC13()  const noexcept{ return getFloat(_morder[0][2],_index);}
    float getC14()  const noexcept{ return getFloat(_morder[0][3],_index);}
    float getC15()  const noexcept{ return getFloat(_morder[0][4],_index);}
    float getC22()  const noexcept{ return getFloat(_morder[1][1],_index);}
    float getC23()  const noexcept{ return getFloat(_morder[1][2],_index);}
    float getC24()  const noexcept{ return getFloat(_morder[1][3],_index);}
    float getC25()  const noexcept{ return getFloat(_morder[1][4],_index);}
    float getC33()  const noexcept{ return getFloat(_morder[2][2],_index);}
    float getC34()  const noexcept{ return getFloat(_morder[2][3],_index);}
    float getC35()  const noexcept{ return getFloat(_morder[2][4],_index);}
    float getC44()  const noexcept{ return getFloat(_morder[3][3],_index);}
    float getC45()  const noexcept{ return getFloat(_morder[3][4],_index);}
    float getC55()  const noexcept{ return getFloat(_morder[4][4],_index);}

    
    void notify() final{
      bank::notify();
      scanIndex();
    }
    void scanIndex();
    void print();
    
    private :
    
    CovMatrix _matrix={5,std::vector<float>(5,0)};
    IndMatrix _morder={5,std::vector<int>(5,-1)};
    int  _pindex_order=-1;
    mutable int _index=-1;
    //std::map<int,int> _rmap;
    std::vector<int> _rvec;
  };
 
  using covmat_ptr=clas12::covmatrix*;
  using covmat_uptr=std::unique_ptr<clas12::covmatrix>;
}
 
#endif /* UTILS_H */