#include <HIPUserVariablesIORoot.h>

Inheritance diagram for HIPUserVariablesIORoot:

Public Types
using	Alignables = align::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) override

void	createBranches (void) override
	create root branches More...

int	findEntry (unsigned int detId, int comp)

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

AlignmentUserVariables *	readOne (Alignable *ali, int &ierr) override

void	setBranchAddresses (void) override
	set root branches More...

int	writeOne (Alignable *ali) override

Private Attributes
double	AlignableChi2

unsigned int	AlignableNdof

int	DataType

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 = 19

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
const static int	itermax = 1000

const static int	nParMax = 20

Detailed Description

concrete class for ROOT based IO of AlignmentUserVariables

Definition at line 10 of file HIPUserVariablesIORoot.h.

Member Typedef Documentation

◆ Alignables

using HIPUserVariablesIORoot::Alignables = align::Alignables

Definition at line 12 of file HIPUserVariablesIORoot.h.

◆ treemaptype

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

private

Definition at line 65 of file HIPUserVariablesIORoot.h.

Constructor & Destructor Documentation

◆ HIPUserVariablesIORoot()

HIPUserVariablesIORoot::HIPUserVariablesIORoot ( )

constructor

Definition at line 15 of file HIPUserVariablesIORoot.cc.

     : ObjId(0), Id(0), Nhit(0), Nparj(0), Npare(0), DataType(-1), 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.;
   }
 }

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

Member Function Documentation

◆ close()

int HIPUserVariablesIORoot::close ( void )

inlineoverrideprivatevirtual

close IO

Implements AlignmentUserVariablesIO.

Definition at line 40 of file HIPUserVariablesIORoot.h.

40 { return closeRoot(); }

References AlignmentIORootBase::closeRoot().

Referenced by esMonitoring.AsyncLineReaderMixin::handle_close(), esMonitoring.FDJsonServer::handle_close(), readHIPUserVariables(), and writeHIPUserVariables().

◆ createBranches()

void HIPUserVariablesIORoot::createBranches ( void )

overrideprivatevirtual

create root branches

Implements AlignmentIORootBase.

Definition at line 31 of file HIPUserVariablesIORoot.cc.

                                                 {
   tree->Branch("Id", &Id, "Id/i");
   tree->Branch("ObjId", &ObjId, "ObjId/I");
  
   tree->Branch("Nhit", &Nhit, "Nhit/I");
   tree->Branch("DataType", &DataType, "DataType/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");
 }

References AlignableChi2, AlignableNdof, DataType, Id, Jtve, Jtvj, Nhit, Npare, Nparj, ObjId, Par, and ParError.

◆ findEntry()

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

private

Definition at line 68 of file HIPUserVariablesIORoot.cc.

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

References AlCaHLTBitMon_QueryRunRegistry::comp, ev, Id, newopen, ObjId, mps_fire::result, and treemap.

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

◆ open()

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

inlineoverrideprivatevirtual

open IO

Implements AlignmentUserVariablesIO.

Definition at line 34 of file HIPUserVariablesIORoot.h.

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

References corrVsCorr::filename, align_cfg::iteration, newopen, and AlignmentIORootBase::openRoot().

Referenced by readHIPUserVariables(), and writeHIPUserVariables().

◆ readHIPUserVariables()

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

read user variables

Definition at line 207 of file HIPUserVariablesIORoot.cc.

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

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

Referenced by HIPAlignmentAlgorithm::collector().

◆ readOne()

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

overrideprivatevirtual

read AlignmentParameters of one Alignable

Implements AlignmentUserVariablesIO.

Definition at line 140 of file HIPUserVariablesIORoot.cc.

                                                                                  {
   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 = row; col < np; col++) {
         jtvj[row][col] = Jtvj[count];
         count++;
       }
     }
  
     uvar = new HIPUserVariables(np);
     uvar->jtvj = jtvj;
     uvar->jtve = jtve;
     uvar->nhit = Nhit;
     uvar->datatype = DataType;
     uvar->alipar = alipar;
     uvar->alierr = alierr;
     //Chi2n
     uvar->alichi2 = AlignableChi2;
     uvar->alindof = AlignableNdof;
  
     return uvar;
   }
  
   //  ierr=-1;
   return nullptr;
 }

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

◆ setBranchAddresses()

void HIPUserVariablesIORoot::setBranchAddresses ( void )

overrideprivatevirtual

set root branches

Implements AlignmentIORootBase.

Definition at line 49 of file HIPUserVariablesIORoot.cc.

                                                     {
   tree->SetBranchAddress("Id", &Id);
   tree->SetBranchAddress("ObjId", &ObjId);
  
   tree->SetBranchAddress("Nhit", &Nhit);
   tree->SetBranchAddress("DataType", &DataType);
   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);
 }

References AlignableChi2, AlignableNdof, DataType, Id, Jtve, Jtvj, Nhit, Npare, Nparj, ObjId, Par, and ParError.

◆ writeHIPUserVariables()

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

write user variables

Definition at line 184 of file HIPUserVariablesIORoot.cc.

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

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

Referenced by HIPAlignmentAlgorithm::terminate().

◆ writeOne()

int HIPUserVariablesIORoot::writeOne ( Alignable * ali )

overrideprivatevirtual

write AlignmentParameters of one Alignable

Implements AlignmentUserVariablesIO.

Definition at line 96 of file HIPUserVariablesIORoot.cc.

                                                    {
   AlignmentParameters* ap = ali->alignmentParameters();
  
   if ((ap->userVariables()) == nullptr) {
     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 datatype = uvar->datatype;
   int np = jtve.num_row();
  
   Nhit = nhit;
   DataType = datatype;
   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 = row; col < np; 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;
 }

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

Member Data Documentation

◆ AlignableChi2

double HIPUserVariablesIORoot::AlignableChi2

private

Definition at line 59 of file HIPUserVariablesIORoot.h.

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

◆ AlignableNdof

unsigned int HIPUserVariablesIORoot::AlignableNdof

private

Definition at line 60 of file HIPUserVariablesIORoot.h.

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

◆ DataType

int HIPUserVariablesIORoot::DataType

private

Definition at line 56 of file HIPUserVariablesIORoot.h.

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

◆ Id

unsigned int HIPUserVariablesIORoot::Id

private

Definition at line 54 of file HIPUserVariablesIORoot.h.

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

◆ Jtve

double HIPUserVariablesIORoot::Jtve[nparmax]

private

Definition at line 58 of file HIPUserVariablesIORoot.h.

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

◆ Jtvj

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

private

Definition at line 57 of file HIPUserVariablesIORoot.h.

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

◆ newopen

bool HIPUserVariablesIORoot::newopen

private

Definition at line 64 of file HIPUserVariablesIORoot.h.

Referenced by findEntry(), and open().

◆ Nhit

int HIPUserVariablesIORoot::Nhit

private

Definition at line 55 of file HIPUserVariablesIORoot.h.

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

◆ Npare

int HIPUserVariablesIORoot::Npare

private

Definition at line 55 of file HIPUserVariablesIORoot.h.

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

◆ Nparj

int HIPUserVariablesIORoot::Nparj

private

Definition at line 55 of file HIPUserVariablesIORoot.h.

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

◆ nparmax

const int HIPUserVariablesIORoot::nparmax = 19

staticprivate

Definition at line 50 of file HIPUserVariablesIORoot.h.

Referenced by HIPUserVariablesIORoot().

◆ ObjId

int HIPUserVariablesIORoot::ObjId

private

alignment parameter tree

Definition at line 53 of file HIPUserVariablesIORoot.h.

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

◆ Par

double HIPUserVariablesIORoot::Par[nparmax]

private

Definition at line 61 of file HIPUserVariablesIORoot.h.

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

◆ ParError

double HIPUserVariablesIORoot::ParError[nparmax]

private

Definition at line 62 of file HIPUserVariablesIORoot.h.

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

◆ treemap

treemaptype HIPUserVariablesIORoot::treemap

private

Definition at line 66 of file HIPUserVariablesIORoot.h.

Referenced by findEntry().

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ Alignables

◆ treemaptype

Constructor & Destructor Documentation

◆ HIPUserVariablesIORoot()

Member Function Documentation

◆ close()

◆ createBranches()

◆ findEntry()

◆ open()

◆ readHIPUserVariables()

◆ readOne()

◆ setBranchAddresses()

◆ writeHIPUserVariables()

◆ writeOne()

Member Data Documentation

◆ AlignableChi2

◆ AlignableNdof

◆ DataType

◆ Id

◆ Jtve

◆ Jtvj

◆ newopen

◆ Nhit

◆ Npare

◆ Nparj

◆ nparmax

◆ ObjId

◆ Par

◆ ParError

◆ treemap