HIPUserVariablesIORoot Class Reference

#include <HIPUserVariablesIORoot.h>

Inheritance diagram for HIPUserVariablesIORoot:

Public Types
typedef std::vector< Alignable * >	Alignables

Public Member Functions
	HIPUserVariablesIORoot ()
std::vector < AlignmentUserVariables * >	readHIPUserVariables (const Alignables &alivec, const char *filename, int iter, int &ierr)
void	writeHIPUserVariables (const Alignables &alivec, const char *filename, int iter, bool validCheck, int &ierr)

Private Types
typedef std::map< std::pair < int, int >, int >	treemaptype

Private Member Functions
int	close (void)
void	createBranches (void)
	create root branches More...
int	findEntry (unsigned int detId, int comp)
int	open (const char *filename, int iteration, bool writemode)
AlignmentUserVariables *	readOne (Alignable *ali, int &ierr)
void	setBranchAddresses (void)
	set root branches More...
int	writeOne (Alignable *ali)

Private Attributes
double	AlignableChi2
unsigned int	AlignableNdof
unsigned int	Id
double	Jtve [nparmax]
double	Jtvj [nparmax *(nparmax+1)/2]
bool	newopen
int	Nhit
int	Npare
int	Nparj
int	ObjId
double	Par [nparmax]
double	ParError [nparmax]
treemaptype	treemap

Static Private Attributes
static const int	nparmax =6

Additional Inherited Members
Protected Member Functions inherited from AlignmentIORootBase
	AlignmentIORootBase ()
	constructor More...
int	closeRoot (void)
	close IO More...
int	openRoot (const char *filename, int iteration, bool writemode)
	open IO More...
int	testFile (const char *filename, const TString &tname)
	test if file is existing and if so, what the highest iteration is More...
TString	treeName (int iter, const TString &tname)
	compose tree name More...
virtual	~AlignmentIORootBase ()
	destructor More...
Protected Member Functions inherited from AlignmentUserVariablesIO
std::vector < AlignmentUserVariables * >	read (const align::Alignables &alivec, int &ierr)
int	write (const align::Alignables &alivec, bool validCheck)
virtual	~AlignmentUserVariablesIO ()
Protected Attributes inherited from AlignmentIORootBase
bool	bWrite
TTree *	tree
TString	treename
TString	treetxt
Static Protected Attributes inherited from AlignmentIORootBase
static const int	itermax = 1000
static const int	nParMax = 20

Detailed Description

concrete class for ROOT based IO of AlignmentUserVariables

Definition at line 9 of file HIPUserVariablesIORoot.h.

Member Typedef Documentation

typedef std::vector<Alignable*> HIPUserVariablesIORoot::Alignables

Definition at line 15 of file HIPUserVariablesIORoot.h.

typedef std::map< std::pair<int,int> , int > HIPUserVariablesIORoot::treemaptype

private

Definition at line 67 of file HIPUserVariablesIORoot.h.

Constructor & Destructor Documentation

HIPUserVariablesIORoot::HIPUserVariablesIORoot

(

)

constructor

Definition at line 16 of file HIPUserVariablesIORoot.cc.

References i, Jtve, Jtvj, nparmax, Par, ParError, AlignmentIORootBase::treename, and AlignmentIORootBase::treetxt.

                                               :
  ObjId(0), Id(0),Nhit(0), Nparj(0), Npare(0),
  AlignableChi2(0.), AlignableNdof(0)
{
  treename = "T9";
  treetxt = "HIP User Variables";
  
  for (int i=0;i<nparmax*(nparmax+1)/2;++i) 
    Jtvj[i] = 0.;
  for (int i=0;i<nparmax;++i){ 
    Jtve[i] = 0.;
    Par[i] = 0.; 
    ParError[i] = 0.;}
}

Member Function Documentation

int HIPUserVariablesIORoot::close ( void  )

inlineprivatevirtual

close IO

Implements AlignmentUserVariablesIO.

Definition at line 43 of file HIPUserVariablesIORoot.h.

References AlignmentIORootBase::closeRoot().

Referenced by lumiQTWidget.ApplicationWindow::fileQuit(), esMonitoring.AsyncLineReaderMixin::handle_close(), esMonitoring.FDJsonServer::handle_close(), Vispa.Gui.BoxContentDialog.BoxContentDialog::keyPressEvent(), Vispa.Gui.FindDialog.FindDialog::keyPressEvent(), readHIPUserVariables(), and writeHIPUserVariables().

{return closeRoot();};

void HIPUserVariablesIORoot::createBranches ( void  )

privatevirtual

create root branches

Implements AlignmentIORootBase.

Definition at line 33 of file HIPUserVariablesIORoot.cc.

References AlignableChi2, AlignableNdof, Id, Jtve, Jtvj, Nhit, Npare, Nparj, ObjId, Par, ParError, and AlignmentIORootBase::tree.

{
  tree->Branch("Id",        &Id,        "Id/i");
  tree->Branch("ObjId",     &ObjId,     "ObjId/I");

  tree->Branch("Nhit",      &Nhit,      "Nhit/I");
  tree->Branch("Nparj",     &Nparj,     "Nparj/I");
  tree->Branch("Jtvj",      &Jtvj,      "Jtvj[Nparj]/D");
  tree->Branch("Npare",     &Npare,     "Npare/I");
  tree->Branch("Jtve",      &Jtve,      "Jtve[Npare]/D");
  tree->Branch("AlignableChi2",          &AlignableChi2, "AlignableChi2/D");
  tree->Branch("AlignableNdof",          &AlignableNdof, "AlignableNdof/i");
  tree->Branch("Par",         &Par, "Par[Npare]/D");
  tree->Branch("ParError",         &ParError, "ParError[Npare]/D");
}

int HIPUserVariablesIORoot::findEntry ( unsigned int  detId, int  comp  )

private

Definition at line 70 of file HIPUserVariablesIORoot.cc.

References ev, Id, newopen, ObjId, mps_fire::result, AlignmentIORootBase::tree, and treemap.

Referenced by QGLikelihoodCalculator.QGLikelihoodCalculator::computeQGLikelihood(), and readOne().

{
  if (newopen) { // we're here for the first time
    edm::LogInfo("Alignment") <<"[HIPUserVariablesIORoot::findEntry] fill map ...";
    treemap.erase(treemap.begin(),treemap.end());
    for (int ev = 0;ev<tree->GetEntries();ev++) {
      tree->GetEntry(ev); 
      treemap[std::make_pair(Id,ObjId)]=ev;
    }
    newopen=false;
  }
  
  // now we have filled the map
  treemaptype::iterator imap = treemap.find(std::make_pair(detId,comp));
  int result=-1;
  if (imap != treemap.end()) result=(*imap).second;
  return result;


  //double noAliPar = tree->GetEntries();
  //for (int ev = 0;ev<noAliPar;ev++) {
  //  tree->GetEntry(ev); 
  //  if(Id==detId&&comp==ObjId) return (ev);
  //}
  //return(-1);
}

int HIPUserVariablesIORoot::open ( const char *  filename, int  iteration, bool  writemode  )

inlineprivatevirtual

open IO

Implements AlignmentUserVariablesIO.

Definition at line 39 of file HIPUserVariablesIORoot.h.

References newopen, and AlignmentIORootBase::openRoot().

Referenced by Vispa.Plugins.ConfigEditor.ConfigEditorTabController.ConfigEditorTabController::importConfig(), readHIPUserVariables(), and writeHIPUserVariables().

    {newopen=true; return openRoot(filename,iteration,writemode);};

std::vector< AlignmentUserVariables * > HIPUserVariablesIORoot::readHIPUserVariables	(	const Alignables &	alivec,
		const char *	filename,
		int	iter,
		int &	ierr
	)

read user variables

Definition at line 203 of file HIPUserVariablesIORoot.cc.

References close(), open(), AlignmentUserVariablesIO::read(), and mps_fire::result.

Referenced by HIPAlignmentAlgorithm::collector().

{
  std::vector<AlignmentUserVariables*> result;
  ierr=0;
  int iret;
  iret = open(filename,iter,false);
  if (iret!=0) { ierr=-1; return result;}
  result = read(alivec,iret);
  if (iret!=0) { ierr=-2; return result;}
  iret = close();
  if (iret!=0) { ierr=-3; return result;}

  return result;
}

AlignmentUserVariables * HIPUserVariablesIORoot::readOne ( Alignable *  ali, int &  ierr  )

privatevirtual

read AlignmentParameters of one Alignable

Implements AlignmentUserVariablesIO.

Definition at line 142 of file HIPUserVariablesIORoot.cc.

References HIPUserVariables::alichi2, HIPUserVariables::alierr, AlignableChi2, AlignableNdof, Alignable::alignableObjectId(), HIPUserVariables::alindof, HIPUserVariables::alipar, cuy::col, KineDebug3::count(), mps_splice::entry, findEntry(), Alignable::id(), HIPUserVariables::jtve, Jtve, HIPUserVariables::jtvj, Jtvj, HIPUserVariables::nhit, Nhit, np, Npare, Par, ParError, and AlignmentIORootBase::tree.

{
  ierr=0;
  HIPUserVariables* uvar;

  int entry = findEntry(ali->id(), ali->alignableObjectId());
  if(entry!=-1) {
    tree->GetEntry(entry);

    int np=Npare;
    AlgebraicVector jtve(np,0);
    AlgebraicSymMatrix jtvj(np,0);
    AlgebraicVector alipar(np,0);
    AlgebraicVector alierr(np,0);
    int count=0;
    for(int row=0;row<np;row++) {
      jtve[row]=Jtve[row];
      alipar[row]=Par[row];
      alierr[row]=ParError[row];
      for(int col=0; col < np;col++) {
    if(row-1<col) {jtvj[row][col]=Jtvj[count];count++;}
      }
    } 

    uvar = new HIPUserVariables(np);
    uvar->jtvj=jtvj;
    uvar->jtve=jtve;
    uvar->nhit=Nhit;
    uvar->alipar=alipar;
    uvar->alierr=alierr;
    //Chi2n
    uvar->alichi2=AlignableChi2;
    uvar->alindof=AlignableNdof;

    return uvar;
  }

  //  ierr=-1;
  return 0 ;
}

void HIPUserVariablesIORoot::setBranchAddresses ( void  )

privatevirtual

set root branches

Implements AlignmentIORootBase.

Definition at line 51 of file HIPUserVariablesIORoot.cc.

References AlignableChi2, AlignableNdof, Id, Jtve, Jtvj, Nhit, Npare, Nparj, ObjId, Par, ParError, and AlignmentIORootBase::tree.

{
  tree->SetBranchAddress("Id",        &Id);
  tree->SetBranchAddress("ObjId",     &ObjId);

  tree->SetBranchAddress("Nhit",      &Nhit);
  tree->SetBranchAddress("Nparj",     &Nparj);
  tree->SetBranchAddress("Jtvj",      &Jtvj);
  tree->SetBranchAddress("Npare",     &Npare);
  tree->SetBranchAddress("Jtve",      &Jtve);
  tree->SetBranchAddress("AlignableChi2",     &AlignableChi2);
  tree->SetBranchAddress("AlignableNdof",      &AlignableNdof);
  tree->SetBranchAddress("Par",     &Par);
  tree->SetBranchAddress("ParError",     &ParError);
}

void HIPUserVariablesIORoot::writeHIPUserVariables	(	const Alignables &	alivec,
		const char *	filename,
		int	iter,
		bool	validCheck,
		int &	ierr
	)

write user variables

Definition at line 187 of file HIPUserVariablesIORoot.cc.

References close(), open(), and AlignmentUserVariablesIO::write().

Referenced by HIPAlignmentAlgorithm::terminate().

{
  ierr=0;
  int iret;
  iret = open(filename,iter,true);
  if (iret!=0) { ierr=-1; return;}
  iret = write(alivec,validCheck);
  if (iret!=0) { ierr=-2; return;}
  iret = close();
  if (iret!=0) { ierr=-3; return;}
}

int HIPUserVariablesIORoot::writeOne ( Alignable *  ali )

privatevirtual

write AlignmentParameters of one Alignable

Implements AlignmentUserVariablesIO.

Definition at line 99 of file HIPUserVariablesIORoot.cc.

References HIPUserVariables::alichi2, HIPUserVariables::alierr, AlignableChi2, AlignableNdof, Alignable::alignableObjectId(), Alignable::alignmentParameters(), HIPUserVariables::alindof, HIPUserVariables::alipar, cuy::col, KineDebug3::count(), Id, Alignable::id(), HIPUserVariables::jtve, Jtve, HIPUserVariables::jtvj, Jtvj, HIPUserVariables::nhit, Nhit, np, Npare, Nparj, ObjId, Par, ParError, AlignmentIORootBase::tree, and AlignmentParameters::userVariables().

{
  AlignmentParameters* ap=ali->alignmentParameters();

  if ((ap->userVariables())==0) { 
    edm::LogError("Alignment") <<"UserVariables not found!"; 
    return -1; 
  }

  HIPUserVariables* uvar = 
    dynamic_cast<HIPUserVariables*>(ap->userVariables());

  AlgebraicSymMatrix jtvj = uvar->jtvj;
  AlgebraicVector jtve = uvar->jtve;
  AlgebraicVector alipar = uvar->alipar;
  AlgebraicVector alierr = uvar->alierr;
  int nhit=uvar->nhit;
  int np=jtve.num_row();

  Nhit=nhit;
  Npare=np;
  Nparj=np*(np+1)/2;
  int count=0;
  for(int row=0;row<np;row++){
    Jtve[row]=jtve[row];
    Par[row]=alipar[row];
    ParError[row]=alierr[row];
    for(int col=0;col<np;col++){
      if(row-1<col){Jtvj[count]=jtvj[row][col];count++;}
    }
  }
  Id = ali->id();
  ObjId = ali->alignableObjectId();

  //Chi^2 of alignable
  AlignableChi2= uvar->alichi2 ;
  AlignableNdof= uvar->alindof ;
  tree->Fill();
  return 0;
}

Member Data Documentation

double HIPUserVariablesIORoot::AlignableChi2

private

Definition at line 61 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), readOne(), setBranchAddresses(), and writeOne().

unsigned int HIPUserVariablesIORoot::AlignableNdof

private

Definition at line 62 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), readOne(), setBranchAddresses(), and writeOne().

unsigned int HIPUserVariablesIORoot::Id

private

Definition at line 57 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), findEntry(), setBranchAddresses(), and writeOne().

double HIPUserVariablesIORoot::Jtve[nparmax]

private

Definition at line 60 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), HIPUserVariablesIORoot(), readOne(), setBranchAddresses(), and writeOne().

double HIPUserVariablesIORoot::Jtvj[nparmax *(nparmax+1)/2]

private

Definition at line 59 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), HIPUserVariablesIORoot(), readOne(), setBranchAddresses(), and writeOne().

bool HIPUserVariablesIORoot::newopen

private

Definition at line 66 of file HIPUserVariablesIORoot.h.

Referenced by findEntry(), and open().

int HIPUserVariablesIORoot::Nhit

private

Definition at line 58 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), readOne(), setBranchAddresses(), and writeOne().

int HIPUserVariablesIORoot::Npare

private

Definition at line 58 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), readOne(), setBranchAddresses(), and writeOne().

int HIPUserVariablesIORoot::Nparj

private

Definition at line 58 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), setBranchAddresses(), and writeOne().

const int HIPUserVariablesIORoot::nparmax =6

staticprivate

Definition at line 53 of file HIPUserVariablesIORoot.h.

Referenced by HIPUserVariablesIORoot().

int HIPUserVariablesIORoot::ObjId

private

alignment parameter tree

Definition at line 56 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), findEntry(), setBranchAddresses(), and writeOne().

double HIPUserVariablesIORoot::Par[nparmax]

private

Definition at line 63 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), HIPUserVariablesIORoot(), readOne(), setBranchAddresses(), and writeOne().

double HIPUserVariablesIORoot::ParError[nparmax]

private

Definition at line 64 of file HIPUserVariablesIORoot.h.

Referenced by createBranches(), HIPUserVariablesIORoot(), readOne(), setBranchAddresses(), and writeOne().

treemaptype HIPUserVariablesIORoot::treemap

private

Definition at line 68 of file HIPUserVariablesIORoot.h.

Referenced by findEntry().