HIPUserVariablesIORoot.cc

Go to the documentation of this file.

#include "TTree.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "Alignment/CommonAlignment/interface/Alignable.h"
#include "Alignment/CommonAlignment/interface/AlignmentParameters.h"
#include "Alignment/HIPAlignmentAlgorithm/interface/HIPUserVariables.h"

// this class's header
#include "Alignment/HIPAlignmentAlgorithm/interface/HIPUserVariablesIORoot.h"

// ----------------------------------------------------------------------------
// constructor

HIPUserVariablesIORoot::HIPUserVariablesIORoot()
    : 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.;
  }
}

// ----------------------------------------------------------------------------

void HIPUserVariablesIORoot::createBranches(void) {
  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");
}

// ----------------------------------------------------------------------------

void HIPUserVariablesIORoot::setBranchAddresses(void) {
  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);
}

// ----------------------------------------------------------------------------
// find tree entry based on detID and typeID

int HIPUserVariablesIORoot::findEntry(unsigned int detId, int comp) {
  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::writeOne(Alignable* ali) {
  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;
}

// ----------------------------------------------------------------------------

AlignmentUserVariables* HIPUserVariablesIORoot::readOne(Alignable* ali, int& ierr) {
  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;
}

//-----------------------------------------------------------------------------

void HIPUserVariablesIORoot::writeHIPUserVariables(
    const Alignables& alivec, const char* filename, int iter, bool validCheck, int& ierr) {
  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;
  }
}

//-----------------------------------------------------------------------------

std::vector<AlignmentUserVariables*> HIPUserVariablesIORoot::readHIPUserVariables(const Alignables& alivec,
                                                                                  const char* filename,
                                                                                  int iter,
                                                                                  int& ierr) {
  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;
}