TRTAlign::HistoChisqProjection Class Reference

#include <TRTAlignHistoChisqProjection.h>

Collaboration diagram for TRTAlign::HistoChisqProjection:

Public Member Functions
	HistoChisqProjection (const char *name, int dim, int nbins, float xmin, float xmax)
	~HistoChisqProjection ()
	HistoChisqProjection (const HistoChisqProjection &)=delete
HistoChisqProjection &	operator= (const HistoChisqProjection &)=delete
void	add (const CLHEP::HepMatrix &deriv, const CLHEP::HepSymMatrix &weight, const CLHEP::HepSymMatrix &variance, const CLHEP::HepVector &residual)
void	add (const CLHEP::HepMatrix &deriv, const CLHEP::HepSymMatrix &weight, const CLHEP::HepSymMatrix &variance, const CLHEP::HepVector &residual, const CLHEP::HepSymMatrix &secondweight)
CLHEP::HepVector	getFittedHalfDChisqDX (std::ostream &os=std::cout) const
CLHEP::HepVector	getNormalHalfDChisqDX () const
void	setDirectory (TDirectory *dir)

Private Member Functions
void	addfirst (const CLHEP::HepMatrix &derivative, const CLHEP::HepSymMatrix &weight, const CLHEP::HepSymMatrix &variance, const CLHEP::HepVector &residual)
void	addsecond (const CLHEP::HepMatrix &derivative, const CLHEP::HepSymMatrix &weight)

Private Attributes
int	m_dim
TH2 *	m_h2
TH2 *	m_unweightedh2
std::vector< double >	m_integrals
CLHEP::HepVector	m_halfdChisqdX
CLHEP::HepSymMatrix	m_halfd2ChisqdX2

Detailed Description

Definition at line 23 of file TRTAlignHistoChisqProjection.h.

Constructor & Destructor Documentation

HistoChisqProjection() [1/2]

TRTAlign::HistoChisqProjection::HistoChisqProjection	(	const char *	name,
		int	dim,
		int	nbins,
		float	xmin,
		float	xmax
	)

Definition at line 14 of file TRTAlignHistoChisqProjection.cxx.

    : m_dim(dim), m_integrals(dim,0),m_halfdChisqdX(dim,0)
  {
    m_h2 = new TH2F("h2",name,dim,-0.5,dim-0.5,nbins,xmin,xmax) ;
    m_unweightedh2 = new TH2F("unweightedh2",name,dim,-0.5,dim-0.5,nbins,xmin,xmax) ;
    m_h2->Sumw2() ;
  } ;

~HistoChisqProjection()

TRTAlign::HistoChisqProjection::~HistoChisqProjection

(

)

HistoChisqProjection() [2/2]

TRTAlign::HistoChisqProjection::HistoChisqProjection ( const HistoChisqProjection &  )

delete

Member Function Documentation

add() [1/2]

void TRTAlign::HistoChisqProjection::add ( const CLHEP::HepMatrix &  deriv, const CLHEP::HepSymMatrix &  weight, const CLHEP::HepSymMatrix &  variance, const CLHEP::HepVector &  residual  )

inline

Definition at line 33 of file TRTAlignHistoChisqProjection.h.

                                                                                  {
      addfirst( deriv, weight, variance, residual ) ;
      addsecond( deriv, weight ) ;
    }

add() [2/2]

void TRTAlign::HistoChisqProjection::add ( const CLHEP::HepMatrix &  deriv, const CLHEP::HepSymMatrix &  weight, const CLHEP::HepSymMatrix &  variance, const CLHEP::HepVector &  residual, const CLHEP::HepSymMatrix &  secondweight  )

inline

Definition at line 40 of file TRTAlignHistoChisqProjection.h.

                                                 {
      addfirst( deriv, weight, variance, residual ) ;
      addsecond( deriv, secondweight ) ;
    }

addfirst()

void TRTAlign::HistoChisqProjection::addfirst ( const CLHEP::HepMatrix &  derivative, const CLHEP::HepSymMatrix &  weight, const CLHEP::HepSymMatrix &  variance, const CLHEP::HepVector &  residual  )

private

Definition at line 22 of file TRTAlignHistoChisqProjection.cxx.

  {
    int dimr = weight.num_row() ;
    CLHEP::HepMatrix derivtimesweight(deriv*weight) ;
    
    for(int irow=1; irow<=dimr; ++irow) {
      double thissigma = sqrt( variance.fast(irow,irow) ) ;
      
      for(int ipar=0; ipar<m_dim; ++ipar) {
    double thisweight = derivtimesweight(ipar+1, irow ) * thissigma ;
    double thispull   = residual(irow) / thissigma ;
    if( thisweight< 0 ) {
      thisweight*=-1 ;
      thispull*=-1 ;
    }
    m_h2->Fill( double(ipar), thispull, thisweight ) ;
    m_unweightedh2->Fill( double(ipar), thispull ) ;
    m_integrals[ipar] += thisweight ;
    m_halfdChisqdX(ipar+1) += thisweight*thispull ;
      }
    }
  }

addsecond()

void TRTAlign::HistoChisqProjection::addsecond ( const CLHEP::HepMatrix &  derivative, const CLHEP::HepSymMatrix &  weight  )

private

Definition at line 46 of file TRTAlignHistoChisqProjection.cxx.

  {
    m_halfd2ChisqdX2 += weight.similarity(deriv) ;
  }

getFittedHalfDChisqDX()

CLHEP::HepVector TRTAlign::HistoChisqProjection::getFittedHalfDChisqDX

(

std::ostream &

os = std::cout

)

const

Definition at line 67 of file TRTAlignHistoChisqProjection.cxx.

  {
    CLHEP::HepVector halfDChisqDX(m_dim,0) ;
    // perform a fit to each hist
    TF1 f1("mygaus",gausfunc,m_h2->GetYaxis()->GetXmin(),m_h2->GetYaxis()->GetXmax(),3) ;
    for(int ipar=1; ipar<=m_dim; ++ipar) {
      char hisname[256] ;
      sprintf(hisname,"tmph1_%d",ipar) ;
      TH1* tmph1 = m_h2->ProjectionY(hisname,ipar,ipar,"e") ; // the 'e' triggers the error calculation (quite important!)
      binwidth = tmph1->GetXaxis()->GetBinWidth(1) ;
      tmph1->SetDirectory(m_h2->GetDirectory()) ;
      f1.SetParameter(0,tmph1->Integral());
      f1.SetParLimits(0,0.5*tmph1->Integral(),2*tmph1->Integral()) ;
      f1.SetParameter(1,0) ;
      f1.SetParameter(2,tmph1->GetRMS()) ;
      tmph1->Fit(&f1,"Q0I") ;
      tmph1->GetListOfFunctions()->Add(f1.Clone()) ;
      halfDChisqDX(ipar) = f1.GetParameter(1) * m_integrals[ipar-1] ;
      os << "fitresult: " << f1.GetParameter(1) << " weight: " << m_integrals[ipar-1] << " integrated mean: "
     << m_halfdChisqdX(ipar) / m_integrals[ipar-1] << std::endl ;
    }
    return halfDChisqDX ;
  }

getNormalHalfDChisqDX()

CLHEP::HepVector TRTAlign::HistoChisqProjection::getNormalHalfDChisqDX ( ) const

inline

Definition at line 48 of file TRTAlignHistoChisqProjection.h.

{ return m_halfdChisqdX ; } 

operator=()

HistoChisqProjection& TRTAlign::HistoChisqProjection::operator= ( const HistoChisqProjection &  )

delete

setDirectory()

void TRTAlign::HistoChisqProjection::setDirectory

(

TDirectory *

dir

)

Definition at line 91 of file TRTAlignHistoChisqProjection.cxx.

  { 
    m_h2->SetDirectory(dir) ; 
    m_unweightedh2->SetDirectory(dir) ; 
  }

Member Data Documentation

m_dim

int TRTAlign::HistoChisqProjection::m_dim

private

Definition at line 59 of file TRTAlignHistoChisqProjection.h.

m_h2

TH2* TRTAlign::HistoChisqProjection::m_h2

private

Definition at line 60 of file TRTAlignHistoChisqProjection.h.

m_halfd2ChisqdX2

CLHEP::HepSymMatrix TRTAlign::HistoChisqProjection::m_halfd2ChisqdX2

private

Definition at line 64 of file TRTAlignHistoChisqProjection.h.

m_halfdChisqdX

CLHEP::HepVector TRTAlign::HistoChisqProjection::m_halfdChisqdX

private

Definition at line 63 of file TRTAlignHistoChisqProjection.h.

m_integrals

std::vector<double> TRTAlign::HistoChisqProjection::m_integrals

private

Definition at line 62 of file TRTAlignHistoChisqProjection.h.

m_unweightedh2

TH2* TRTAlign::HistoChisqProjection::m_unweightedh2

private

Definition at line 61 of file TRTAlignHistoChisqProjection.h.

---

The documentation for this class was generated from the following files:

• TRTAlignHistoChisqProjection.h
• TRTAlignHistoChisqProjection.cxx