Concrete class for ROOT based IO of Correlations. More...

#include <AlignmentCorrelationsIORoot.h>

Inheritance diagram for AlignmentCorrelationsIORoot:

Private Member Functions
	AlignmentCorrelationsIORoot ()
	constructor More...

int	close (void) override
	close IO More...

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

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

align::Correlations	read (const align::Alignables &alivec, int &ierr) override
	read correlations More...

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

int	write (const align::Correlations &cor, bool validCheck) override
	write correlations More...

Private Attributes
align::ID	Ali1Id
	correlation tree More...

align::StructureType	Ali1ObjId

align::ID	Ali2Id

align::StructureType	Ali2ObjId

double	CorMatrix [nParMax *nParMax]

int	corSize

Friends
class	AlignmentIORoot

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 AlignmentCorrelationsIO
virtual	~AlignmentCorrelationsIO ()
	destructor More...

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 Correlations.

Definition at line 10 of file AlignmentCorrelationsIORoot.h.

Constructor & Destructor Documentation

◆ AlignmentCorrelationsIORoot()

AlignmentCorrelationsIORoot::AlignmentCorrelationsIORoot ( )

private

constructor

Definition at line 13 of file AlignmentCorrelationsIORoot.cc.

References AlignmentIORootBase::treename, and AlignmentIORootBase::treetxt.

                                                          {
   treename = "AlignmentCorrelations";
   treetxt = "Correlations";
 }

Member Function Documentation

◆ close()

int AlignmentCorrelationsIORoot::close ( void )

inlineoverrideprivatevirtual

close IO

Implements AlignmentCorrelationsIO.

Definition at line 23 of file AlignmentCorrelationsIORoot.h.

References AlignmentIORootBase::closeRoot().

Referenced by esMonitoring.AsyncLineReaderMixin::handle_close(), esMonitoring.FDJsonServer::handle_close(), AlignmentIORoot::readCorrelations(), and AlignmentIORoot::writeCorrelations().

23 { return closeRoot(); };

AlignmentIORootBase::closeRoot

int closeRoot(void)

close IO

Definition: AlignmentIORootBase.cc:92

◆ createBranches()

void AlignmentCorrelationsIORoot::createBranches ( void )

overrideprivatevirtual

create root branches

Implements AlignmentIORootBase.

Definition at line 20 of file AlignmentCorrelationsIORoot.cc.

References Ali1Id, Ali1ObjId, Ali2Id, Ali2ObjId, CorMatrix, and corSize.

                                                      {
   tree->Branch("Ali1Id", &Ali1Id, "Ali1Id/i");
   tree->Branch("Ali2Id", &Ali2Id, "Ali2Id/i");
   tree->Branch("Ali1ObjId", &Ali1ObjId, "Ali1ObjId/I");
   tree->Branch("Ali2ObjId", &Ali2ObjId, "Ali2ObjId/I");
   tree->Branch("corSize", &corSize, "corSize/I");
   tree->Branch("CorMatrix", &CorMatrix, "CorMatrix[corSize]/D");
 }

◆ open()

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

inlineoverrideprivatevirtual

open IO

Implements AlignmentCorrelationsIO.

Definition at line 18 of file AlignmentCorrelationsIORoot.h.

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

Referenced by AlignmentIORoot::readCorrelations(), and AlignmentIORoot::writeCorrelations().

                                                                          {
     return openRoot(filename, iteration, writemode);
   };

◆ read()

align::Correlations AlignmentCorrelationsIORoot::read	(	const align::Alignables &	alivec,
		int &	ierr
	)

overrideprivatevirtual

read correlations

Implements AlignmentCorrelationsIO.

Definition at line 71 of file AlignmentCorrelationsIORoot.cc.

References Ali1Id, Ali1ObjId, Ali2Id, Ali2ObjId, cuy::col, CorMatrix, mps_splice::entry, contentValuesCheck::nfound, and AlignmentIORootBase::nParMax.

Referenced by edmIntegrityCheck.PublishToFileSystem::get(), and AlignmentIORoot::readCorrelations().

                                                                                             {
   align::Correlations theMap;
 
   // create ID map for all Alignables in alivec
   align::Alignables::const_iterator it1;
   std::map<std::pair<unsigned int, int>, Alignable*> idAlis;
   for (it1 = alivec.begin(); it1 != alivec.end(); ++it1)
     idAlis[std::make_pair((*it1)->id(), (*it1)->alignableObjectId())] = (*it1);
 
   std::map<std::pair<unsigned int, int>, Alignable*>::const_iterator aliSearch1;
   std::map<std::pair<unsigned int, int>, Alignable*>::const_iterator aliSearch2;
   int nfound = 0;
   double maxEntry = tree->GetEntries();
   for (int entry = 0; entry < maxEntry; entry++) {
     tree->GetEntry(entry);
     aliSearch1 = idAlis.find(std::make_pair(Ali1Id, Ali1ObjId));
     aliSearch2 = idAlis.find(std::make_pair(Ali2Id, Ali2ObjId));
     if (aliSearch1 != idAlis.end() && aliSearch2 != idAlis.end()) {
       // Alignables for this pair found
       nfound++;
       Alignable* myAli1 = (*aliSearch1).second;
       Alignable* myAli2 = (*aliSearch2).second;
       // FIXME: instead of nParMax in the next few lines one should probably
       //        use something like sqrt(corSize) - but take care of rounding!
       //        I have no time to test... :-( GF
       AlgebraicMatrix mat(nParMax, nParMax);
       for (int row = 0; row < nParMax; row++)
         for (int col = 0; col < nParMax; col++)
           mat[row][col] = CorMatrix[row + col * nParMax];
       theMap[std::make_pair(myAli1, myAli2)] = mat;
     }
   }
 
   edm::LogInfo("AlignmentCorrelationsIORoot") << "Read correlations: all,read: " << alivec.size() << "," << nfound;
 
   ierr = 0;
   return theMap;
 }

◆ setBranchAddresses()

void AlignmentCorrelationsIORoot::setBranchAddresses ( void )

overrideprivatevirtual

set root branches

Implements AlignmentIORootBase.

Definition at line 31 of file AlignmentCorrelationsIORoot.cc.

References Ali1Id, Ali1ObjId, Ali2Id, Ali2ObjId, CorMatrix, and corSize.

                                                          {
   tree->SetBranchAddress("corSize", &corSize);
   tree->SetBranchAddress("Ali1Id", &Ali1Id);
   tree->SetBranchAddress("Ali2Id", &Ali2Id);
   tree->SetBranchAddress("Ali1ObjId", &Ali1ObjId);
   tree->SetBranchAddress("Ali2ObjId", &Ali2ObjId);
   tree->SetBranchAddress("CorMatrix", &CorMatrix);
 }

◆ write()

int AlignmentCorrelationsIORoot::write	(	const align::Correlations &	cor,
		bool	validCheck
	)

overrideprivatevirtual

write correlations

Implements AlignmentCorrelationsIO.

Definition at line 42 of file AlignmentCorrelationsIORoot.cc.

References Ali1Id, Ali1ObjId, Ali2Id, Ali2ObjId, Alignable::alignableObjectId(), Alignable::alignmentParameters(), cuy::col, CorMatrix, corSize, Alignable::id(), AlignmentParameters::isValid(), and ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

Referenced by pkg.AbstractPkg::generate(), querying.connection::write_and_commit(), and AlignmentIORoot::writeCorrelations().

                                                                                     {
   int icount = 0;
 
   for (align::Correlations::const_iterator it = cor.begin(); it != cor.end(); ++it) {
     AlgebraicMatrix mat = (*it).second;
     std::pair<Alignable*, Alignable*> Pair = (*it).first;
     Alignable* ali1 = Pair.first;
     Alignable* ali2 = Pair.second;
     if ((ali1->alignmentParameters()->isValid() && ali2->alignmentParameters()->isValid()) || !(validCheck)) {
       Ali1ObjId = ali1->alignableObjectId();
       Ali2ObjId = ali2->alignableObjectId();
       Ali1Id = ali1->id();
       Ali2Id = ali2->id();
       int maxColumn = mat.num_row();
       corSize = maxColumn * maxColumn;
       for (int row = 0; row < maxColumn; row++)
         for (int col = 0; col < maxColumn; col++)
           CorMatrix[row + col * maxColumn] = mat[row][col];
       tree->Fill();
       icount++;
     }
   }
   edm::LogInfo("AlignmentCorrelationsIORoot") << "Writing correlations: all,written: " << cor.size() << "," << icount;
   return 0;
 }