NtupleManager Class Reference

#include <NtupleManager.h>

Public Member Functions
void	BookNtuple ()
void	FillChi2 ()
void	FillFitParameters (MatrixMeschach *AtWAMatrix)
void	FillMeasurements ()
void	FillNtupleTree ()
void	FillOptObjects (MatrixMeschach *AtWAMatrix)
void	InitNtuple ()
	NtupleManager ()
void	WriteNtuple ()
	~NtupleManager ()

Static Public Member Functions
static NtupleManager *	getInstance ()

Private Member Functions
void	GetGlobalAngles (const CLHEP::HepRotation &rmGlob, double *theta)

Private Attributes
double	Chi2CalibratedParameters
double	Chi2Measurements
TClonesArray *	CloneCopsMeas
TClonesArray *	CloneDistancemeter1DimMeas
TClonesArray *	CloneDistancemeterMeas
TClonesArray *	CloneFitParam
TClonesArray *	CloneOptObject
TClonesArray *	CloneSensor2DMeas
TClonesArray *	CloneTiltmeterMeas
TTree *	CocoaTree
CopsMeas *	CopsMeasA
Distancemeter1DimMeas *	Distancemeter1DimMeasA
DistancemeterMeas *	DistancemeterMeasA
FitParam *	FitParamA
int	NCops
int	NDegreesOfFreedom
int	NDistancemeter
int	NDistancemeter1Dim
int	NFitParameters
int	NOptObjects
int	NSensor2D
int	NTiltmeter
OptObject *	OptObjectA
Sensor2DMeas *	Sensor2DMeasA
TFile *	theRootFile
TiltmeterMeas *	TiltmeterMeasA

Static Private Attributes
static NtupleManager *	instance = nullptr

Detailed Description

Definition at line 20 of file NtupleManager.h.

Constructor & Destructor Documentation

NtupleManager()

NtupleManager::NtupleManager ( )

inline

Definition at line 23 of file NtupleManager.h.

Referenced by getInstance().

{};

~NtupleManager()

NtupleManager::~NtupleManager ( )

inline

Definition at line 24 of file NtupleManager.h.

{};

Member Function Documentation

BookNtuple()

void NtupleManager::BookNtuple

(

)

Definition at line 39 of file NtupleManager.cc.

References Chi2CalibratedParameters, Chi2Measurements, CloneCopsMeas, CloneDistancemeter1DimMeas, CloneDistancemeterMeas, CloneFitParam, CloneOptObject, CloneSensor2DMeas, CloneTiltmeterMeas, CocoaTree, NCops, NDegreesOfFreedom, NDistancemeter, NDistancemeter1Dim, NFitParameters, NOptObjects, NSensor2D, NTiltmeter, and theRootFile.

Referenced by Fit::startFit().

                               {
  theRootFile = new TFile("report.root", "RECREATE", "Simple ROOT Ntuple");

  CocoaTree = new TTree("CocoaTree", "CocoaTree");

  CocoaTree->Branch("Chi2Measurements", &Chi2Measurements, "Chi2Measurements/D");
  CocoaTree->Branch("Chi2CalibratedParameters", &Chi2CalibratedParameters, "Chi2CalibratedParameters/D");
  CocoaTree->Branch("NDegreesOfFreedom", &NDegreesOfFreedom, "NDegreesOfFreedom/I");

  CloneFitParam = new TClonesArray("FitParam");
  CocoaTree->Branch("FitParameters", &CloneFitParam, 32000, 2);
  CocoaTree->Branch("NFitParameters", &NFitParameters, "NFitParameters/I");

  CloneOptObject = new TClonesArray("OptObject");
  CocoaTree->Branch("OptObjects", &CloneOptObject, 32000, 2);
  CocoaTree->Branch("NOptObjects", &NOptObjects, "NOptObjects/I");

  CloneSensor2DMeas = new TClonesArray("Sensor2DMeas");
  CocoaTree->Branch("Sensor2DMeasurements", &CloneSensor2DMeas, 32000, 2);
  CocoaTree->Branch("NSensor2D", &NSensor2D, "NSensor2D/I");

  CloneDistancemeterMeas = new TClonesArray("DistancemeterMeas");
  CocoaTree->Branch("DistancemeterMeasurements", &CloneDistancemeterMeas, 32000, 2);
  CocoaTree->Branch("NDistancemeter", &NDistancemeter, "NDistancemeter/I");

  CloneDistancemeter1DimMeas = new TClonesArray("Distancemeter1DimMeas");
  CocoaTree->Branch("Distancemeter1DimMeasurements", &CloneDistancemeter1DimMeas, 32000, 2);
  CocoaTree->Branch("NDistancemeter1Dim", &NDistancemeter1Dim, "NDistancemeter1Dim/I");

  CloneTiltmeterMeas = new TClonesArray("TiltmeterMeas");
  CocoaTree->Branch("TiltmeterMeasurements", &CloneTiltmeterMeas, 32000, 2);
  CocoaTree->Branch("NTiltmeter", &NTiltmeter, "NTiltmeter/I");

  CloneCopsMeas = new TClonesArray("CopsMeas");
  CocoaTree->Branch("CopsMeasurements", &CloneCopsMeas, 32000, 2);
  CocoaTree->Branch("NCops", &NCops, "NCops/I");

  theRootFile->Add(CocoaTree);

  //   FitParametersTree = new TTree("FitParametersTree","FitParametersTree");
  //   FitParametersTree->Branch("NFitParameters",&NFitParameters,"NFitParameters/I");
  //   BookFitParameters = false;
  //   theRootFile->Add(FitParametersTree);

  //   MeasurementsTree = new TTree("MeasurementsTree","MeasurementsTree");
  //   MeasurementsTree->Branch("NMeasurements",&NMeasurements,"NMeasurements/I");
  //   BookMeasurements = false;
  //   theRootFile->Add(MeasurementsTree);
}

FillChi2()

void NtupleManager::FillChi2

(

)

Definition at line 119 of file NtupleManager.cc.

References Chi2CalibratedParameters, Chi2Measurements, Model::EntryList(), cuy::ii, Model::MeasurementList(), and NDegreesOfFreedom.

Referenced by Fit::fitNextEvent().

                             {
  double chi2meas = 0;
  double chi2cal = 0;
  ALIint nMeas = 0, nUnk = 0;

  //----- Calculate the chi2 of measurements
  std::vector<Measurement*>::const_iterator vmcite;
  for (vmcite = Model::MeasurementList().begin(); vmcite != Model::MeasurementList().end(); ++vmcite) {
    for (ALIuint ii = 0; ii < ALIuint((*vmcite)->dim()); ii++) {
      nMeas++;
      double c2 = ((*vmcite)->value(ii) - (*vmcite)->valueSimulated(ii)) / (*vmcite)->sigma(ii);
      chi2meas += c2 * c2;
    }
  }

  //----- Calculate the chi2 of calibrated parameters
  std::vector<Entry*>::iterator veite;
  for (veite = Model::EntryList().begin(); veite != Model::EntryList().end(); ++veite) {
    if ((*veite)->quality() == 2)
      nUnk++;
    if ((*veite)->quality() == 1) {
      //       std::cout << " " << (*veite)->valueDisplacementByFitting() << " "
      //        << (*veite)->value << " " << (*veite)->sigma() << std::endl;
      double c2 = (*veite)->valueDisplacementByFitting() / (*veite)->sigma();
      chi2cal += c2 * c2;
    }
  }

  Chi2Measurements = chi2meas;
  Chi2CalibratedParameters = chi2cal;
  NDegreesOfFreedom = nMeas - nUnk;
}

FillFitParameters()

void NtupleManager::FillFitParameters

(

MatrixMeschach *

AtWAMatrix

)

Definition at line 154 of file NtupleManager.cc.

References CloneOptObject, gather_cfg::cout, Model::EntryList(), FitParamA, FitParam::FittedSigma, FitParam::FittedValue, cms::cuda::for(), FittedEntry::getEntryName(), FittedEntry::getName(), FittedEntry::getOptOName(), FittedEntry::getOrder(), FittedEntry::getOrigSigma(), FittedEntry::getOrigValue(), FittedEntry::getQuality(), FittedEntry::getSigma(), FittedEntry::getValue(), if(), cuy::ii, FitParam::InitialSigma, FitParam::InitialValue, MatrixMeschach::Mat(), FitParam::Name, OptObject::Name, NFitParameters, NOptObjects, FitParam::OptObjectIndex, FitParam::Quality, and mathSSE::sqrt().

Referenced by Fit::fitNextEvent().

                                                                {
  //   double ParValue[1000], ParError[1000];
  int theMinEntryQuality = 1;
  int ii = 0;
  std::vector<Entry*>::const_iterator vecite;
  for (vecite = Model::EntryList().begin(); vecite != Model::EntryList().end(); ++vecite) {
    //--- Only for good quality parameters (='unk')
    if ((*vecite)->quality() >= theMinEntryQuality) {
      ALIint ipos = (*vecite)->fitPos();
      FittedEntry* fe = new FittedEntry((*vecite), ipos, sqrt(AtWAMatrix->Mat()->me[ipos][ipos]));
      //       if (!BookFitParameters) {
      //    CocoaTree->Branch("NFitParameters",&NFitParameters,"NFitParameters/I:");
      //    ALIstring partype = fe->getName() + "/D";
      //    FitParametersTree->Branch(fe->getName().c_str(), &ParValue[ii], partype.c_str());
      //    ALIstring parerrname = fe->getName() + "_err";
      //    ALIstring parerrtype = parerrname + "/D";
      //    FitParametersTree->Branch(parerrname.c_str(), &ParError[ii], parerrtype.c_str());
      //       }
      //       ParValue[ii] = fe->getValue();
      //       ParError[ii] = fe->getSigma();
      std::cout << "EEE " << (*vecite)->ValueDimensionFactor() << " " << (*vecite)->SigmaDimensionFactor() << " "
                << fe->getOptOName() << " " << fe->getEntryName() << " " << fe->getName() << " " << fe->getOrder()
                << " " << fe->getQuality() << " " << (*vecite)->type() << " " << std::endl;
      FitParamA = new ((*CloneFitParam)[ii]) FitParam();
      FitParamA->Name = fe->getName();
      if (fe->getQuality() == 1)
        FitParamA->Quality = "Calibrated";
      else if (fe->getQuality() == 2)
        FitParamA->Quality = "Unknown";
      for (int no = 0; no < NOptObjects; no++) {
        OptObject* optobj = (OptObject*)CloneOptObject->At(no);
        if (optobj->Name == fe->getOptOName())
          FitParamA->OptObjectIndex = no;
      }
      float DF = 1.;
      if ((*vecite)->type() == "centre" || (*vecite)->type() == "length")
        DF = 1000.;
      FitParamA->InitialValue = DF * fe->getOrigValue() * (*vecite)->ValueDimensionFactor();
      FitParamA->InitialSigma = DF * fe->getOrigSigma() * (*vecite)->SigmaDimensionFactor();
      FitParamA->FittedValue = DF * fe->getValue() * (*vecite)->ValueDimensionFactor();
      FitParamA->FittedSigma = DF * fe->getSigma() * (*vecite)->SigmaDimensionFactor();
      ii++;
    }
  }
  //   BookFitParameters = true;
  NFitParameters = ii;
  //   FitParametersTree->Fill();

  /*  
  //---------- Loop sets of entries
  std::vector< FittedEntriesSet* > theFittedEntriesSets;
  std::vector< FittedEntriesSet* >::const_iterator vfescite;
  std::vector< FittedEntry* >::const_iterator vfecite;
  ALIint jj = 1;
  for( vfescite = theFittedEntriesSets.begin(); vfescite != theFittedEntriesSets.end(); vfescite++) {
    //---------- Loop entries
    if( vfescite == theFittedEntriesSets.begin() ) {
    //----- dump entries names if first set 
      ALIint ii = 0;
      for( vfecite = ((*vfescite)->FittedEntries()).begin(); vfecite != ((*vfescite)->FittedEntries()).end(); vfecite++) {
    ALIstring partype = (*vfecite)->getName() + "/D";
    FitParametersTree->Branch((*vfecite)->getName().c_str(), &ParValue[ii], partype.c_str());
    ALIstring parerrname = (*vfecite)->getName() + "_err";
    ALIstring parerrtype = parerrname + "/D";
    FitParametersTree->Branch(parerrname.c_str(), &ParError[ii], parerrtype.c_str());
        ii++;
      }
    }
    ALIint ii = 0;
    for( vfecite = ((*vfescite)->FittedEntries()).begin(); vfecite != ((*vfescite)->FittedEntries()).end(); vfecite++) {
      ParValue[ii] = (*vfecite)->getValue();
      ParError[ii] = (*vfecite)->getSigma();
      ii++;
    }
    NFitParameters = ii;
    FitParametersTree->Fill();
    jj++;
  }
  */
}

FillMeasurements()

void NtupleManager::FillMeasurements

(

)

Definition at line 287 of file NtupleManager.cc.

References TiltmeterMeas::AngError, TiltmeterMeas::Angle, CloneOptObject, CopsMeasA, d1, createTree::dd, DistancemeterMeas::DisError, Distancemeter1DimMeas::DisError, DistancemeterMeas::Distance, Distancemeter1DimMeas::Distance, Distancemeter1DimMeasA, DistancemeterMeasA, mps_fire::i, dqmdumpme::last, Model::MeasurementList(), OptObject::Name, Sensor2DMeas::Name, DistancemeterMeas::Name, Distancemeter1DimMeas::Name, TiltmeterMeas::Name, CopsMeas::Name, NCops, NDistancemeter, NDistancemeter1Dim, NOptObjects, NSensor2D, NTiltmeter, Sensor2DMeas::OptObjectIndex, DistancemeterMeas::OptObjectIndex, Distancemeter1DimMeas::OptObjectIndex, TiltmeterMeas::OptObjectIndex, CopsMeas::OptObjectIndex, Sensor2DMeas::PosError, CopsMeas::PosError, Sensor2DMeas::Position, CopsMeas::Position, Sensor2DMeasA, TiltmeterMeas::SimulatedAngle, DistancemeterMeas::SimulatedDistance, Distancemeter1DimMeas::SimulatedDistance, Sensor2DMeas::SimulatedPosition, CopsMeas::SimulatedPosition, contentValuesCheck::ss, TiltmeterMeasA, and groupFilesInBlocks::tt.

Referenced by Fit::fitNextEvent().

                                     {
  //---------- Loop Measurements
  int ss = 0, dd = 0, d1 = 0, tt = 0, cc = 0;
  std::vector<Measurement*>::const_iterator vmcite;
  for (vmcite = Model::MeasurementList().begin(); vmcite != Model::MeasurementList().end(); ++vmcite) {
    std::vector<ALIstring> optonamelist = (*vmcite)->OptONameList();
    int last = optonamelist.size() - 1;
    ALIstring LastOptOName = optonamelist[last];
    int optoind = -999;
    for (int no = 0; no < NOptObjects; no++) {
      OptObject* optobj = (OptObject*)CloneOptObject->At(no);
      if (optobj->Name == LastOptOName)
        optoind = no;
    }
    //std::cout << "DimSens " << (*vmcite)->type() << " " << (*vmcite)->sigma(0) << " " << LastOptOName << " " << optoind << std::endl;
    if ((*vmcite)->type() == "SENSOR2D") {
      Sensor2DMeasA = new ((*CloneSensor2DMeas)[ss]) Sensor2DMeas();
      Sensor2DMeasA->Name = (*vmcite)->name();
      Sensor2DMeasA->OptObjectIndex = optoind;
      for (ALIuint i = 0; i < (*vmcite)->dim(); i++) {
        Sensor2DMeasA->Position[i] = 1000. * (*vmcite)->value()[i];
        Sensor2DMeasA->PosError[i] = 1000. * (*vmcite)->sigma()[i];
        Sensor2DMeasA->SimulatedPosition[i] = 1000. * (*vmcite)->valueSimulated(i);
      }
      ss++;
    }
    if ((*vmcite)->type() == "DISTANCEMETER") {
      DistancemeterMeasA = new ((*CloneDistancemeterMeas)[dd]) DistancemeterMeas();
      DistancemeterMeasA->Name = (*vmcite)->name();
      DistancemeterMeasA->OptObjectIndex = optoind;
      DistancemeterMeasA->Distance = 1000. * (*vmcite)->value()[0];
      DistancemeterMeasA->DisError = 1000. * (*vmcite)->sigma()[0];
      DistancemeterMeasA->SimulatedDistance = 1000. * (*vmcite)->valueSimulated(0);
      dd++;
    }
    if ((*vmcite)->type() == "DISTANCEMETER1DIM") {
      Distancemeter1DimMeasA = new ((*CloneDistancemeter1DimMeas)[d1]) Distancemeter1DimMeas();
      Distancemeter1DimMeasA->Name = (*vmcite)->name();
      Distancemeter1DimMeasA->OptObjectIndex = optoind;
      Distancemeter1DimMeasA->Distance = 1000. * (*vmcite)->value()[0];
      Distancemeter1DimMeasA->DisError = 1000. * (*vmcite)->sigma()[0];
      Distancemeter1DimMeasA->SimulatedDistance = 1000. * (*vmcite)->valueSimulated(0);
      d1++;
    }
    if ((*vmcite)->type() == "TILTMETER") {
      TiltmeterMeasA = new ((*CloneTiltmeterMeas)[tt]) TiltmeterMeas();
      TiltmeterMeasA->Name = (*vmcite)->name();
      TiltmeterMeasA->OptObjectIndex = optoind;
      TiltmeterMeasA->Angle = (*vmcite)->value()[0];
      TiltmeterMeasA->AngError = (*vmcite)->sigma()[0];
      TiltmeterMeasA->SimulatedAngle = (*vmcite)->valueSimulated(0);
      tt++;
    }
    if ((*vmcite)->type() == "COPS") {
      CopsMeasA = new ((*CloneCopsMeas)[cc]) CopsMeas();
      CopsMeasA->Name = (*vmcite)->name();
      CopsMeasA->OptObjectIndex = optoind;
      for (ALIuint i = 0; i < (*vmcite)->dim(); i++) {
        CopsMeasA->Position[i] = 1000. * (*vmcite)->value()[i];
        CopsMeasA->PosError[i] = 1000. * (*vmcite)->sigma()[i];
        CopsMeasA->SimulatedPosition[i] = 1000. * (*vmcite)->valueSimulated(i);
      }
      cc++;
    }
  }
  NSensor2D = ss;
  NDistancemeter = dd;
  NDistancemeter1Dim = d1;
  NTiltmeter = tt;
  NCops = cc;
  //   MeasurementsTree->Fill();
}

FillNtupleTree()

void NtupleManager::FillNtupleTree

(

)

Definition at line 108 of file NtupleManager.cc.

References CocoaTree.

Referenced by Fit::fitNextEvent().

{ CocoaTree->Fill(); }

FillOptObjects()

void NtupleManager::FillOptObjects

(

MatrixMeschach *

AtWAMatrix

)

Definition at line 237 of file NtupleManager.cc.

References OptObject::AnglesGlobal, OptObject::AnglesLocal, OptObject::CentreGlobal, OptObject::CentreLocal, GetGlobalAngles(), mps_fire::i, cuy::ii, OptObject::Name, NOptObjects, OptObjectA, Model::OptOList(), OptObject::Parent, createTree::pp, theta(), OptObject::Type, XCoor, YCoor, and ZCoor.

Referenced by Fit::fitNextEvent().

                                                             {
  int ii = 0;
  std::vector<OpticalObject*>::const_iterator vecobj;
  for (vecobj = Model::OptOList().begin(); vecobj != Model::OptOList().end(); ++vecobj) {
    OptObjectA = new ((*CloneOptObject)[ii]) OptObject();

    OptObjectA->Name = (*vecobj)->name();
    OptObjectA->Type = (*vecobj)->type();

    if (!(*vecobj)->parent()) {
      OptObjectA->Parent = ii;
      ii++;
      continue;
    }

    int pp = 0;
    std::vector<OpticalObject*>::const_iterator vecobj2;
    for (vecobj2 = Model::OptOList().begin(); vecobj2 != Model::OptOList().end(); ++vecobj2) {
      if ((*vecobj2)->name() == (*vecobj)->parent()->name()) {
        OptObjectA->Parent = pp;
        continue;
      }
      pp++;
    }

    OptObjectA->CentreGlobal[0] = 1000. * (*vecobj)->centreGlobal().x();
    OptObjectA->CentreGlobal[1] = 1000. * (*vecobj)->centreGlobal().y();
    OptObjectA->CentreGlobal[2] = 1000. * (*vecobj)->centreGlobal().z();

    OptObjectA->CentreLocal[0] = 1000. * (*vecobj)->centreLocal().x();
    OptObjectA->CentreLocal[1] = 1000. * (*vecobj)->centreLocal().y();
    OptObjectA->CentreLocal[2] = 1000. * (*vecobj)->centreLocal().z();

    OptObjectA->AnglesLocal[0] = (*vecobj)->getEntryRMangle(XCoor);
    OptObjectA->AnglesLocal[1] = (*vecobj)->getEntryRMangle(YCoor);
    OptObjectA->AnglesLocal[2] = (*vecobj)->getEntryRMangle(ZCoor);

    double theta[3];
    GetGlobalAngles((*vecobj)->rmGlob(), theta);
    for (int i = 0; i < 3; i++)
      OptObjectA->AnglesGlobal[i] = theta[i];

    ii++;
  }

  NOptObjects = ii;
}

GetGlobalAngles()

void NtupleManager::GetGlobalAngles ( const CLHEP::HepRotation &  rmGlob, double *  theta  )

private

Definition at line 362 of file NtupleManager.cc.

References alpha, HLT_2022v12_cff::beta, funct::cos(), gather_cfg::cout, MillePedeFileConverter_cfg::e, CustomPhysics_cfi::gamma, M_PI, funct::sin(), theta(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, geometryCSVtoXML::xz, geometryCSVtoXML::yy, geometryCSVtoXML::yz, and geometryCSVtoXML::zz.

Referenced by FillOptObjects().

                                                                                 {
  double xx = rmGlob.xx();
  if (fabs(xx) < 1.e-08)
    xx = 0.;
  double xy = rmGlob.xy();
  if (fabs(xy) < 1.e-08)
    xy = 0.;
  double xz = rmGlob.xz();
  if (fabs(xz) < 1.e-08)
    xz = 0.;
  double yx = rmGlob.yx();
  if (fabs(yx) < 1.e-08)
    yx = 0.;
  double yy = rmGlob.yy();
  if (fabs(yy) < 1.e-08)
    yy = 0.;
  double yz = rmGlob.yz();
  if (fabs(yz) < 1.e-08)
    yz = 0.;
  double zx = rmGlob.zx();
  if (fabs(zx) < 1.e-08)
    zx = 0.;
  double zy = rmGlob.zy();
  if (fabs(zy) < 1.e-08)
    zy = 0.;
  double zz = rmGlob.zz();
  if (fabs(zz) < 1.e-08)
    zz = 0.;

  double beta = asin(-zx);

  double alpha, gamma;
  if (fabs(zx) != 1.) {
    double sinalpha = zy / cos(beta);
    double cosalpha = zz / cos(beta);
    if (cosalpha >= 0)
      alpha = asin(sinalpha);
    else
      alpha = M_PI - asin(sinalpha);
    if (alpha > M_PI)
      alpha -= 2 * M_PI;

    double singamma = yx / cos(beta);
    double cosgamma = xx / cos(beta);
    if (cosgamma >= 0)
      gamma = asin(singamma);
    else
      gamma = M_PI - asin(singamma);
    if (gamma > M_PI)
      gamma -= 2 * M_PI;

  } else {
    alpha = 0.;

    double singamma = yz / sin(beta);
    double cosgamma = yy;
    if (cosgamma >= 0)
      gamma = asin(singamma);
    else
      gamma = M_PI - asin(singamma);
    if (gamma > M_PI)
      gamma -= 2 * M_PI;
  }

  int GotGlobalAngles = 0;
  if (fabs(xy - (sin(alpha) * sin(beta) * cos(gamma) - sin(gamma) * cos(alpha))) > 1.e-08)
    GotGlobalAngles += 1;
  if (fabs(xz - (cos(alpha) * sin(beta) * cos(gamma) + sin(gamma) * sin(alpha))) > 1.e-08)
    GotGlobalAngles += 10;
  if (fabs(yy - (sin(alpha) * sin(beta) * sin(gamma) + cos(gamma) * cos(alpha))) > 1.e-08)
    GotGlobalAngles += 100;
  if (fabs(yz - (cos(alpha) * sin(beta) * sin(gamma) - cos(gamma) * sin(alpha))) > 1.e-08)
    GotGlobalAngles += 1000;
  if (GotGlobalAngles > 0)
    std::cout << "NtupleManager Warning: cannot get global rotation: " << GotGlobalAngles << std::endl;

  theta[0] = alpha;
  theta[1] = beta;
  theta[2] = gamma;
}

getInstance()

NtupleManager * NtupleManager::getInstance ( )

static

Definition at line 29 of file NtupleManager.cc.

References instance, and NtupleManager().

Referenced by Fit::fitNextEvent(), and Fit::startFit().

                                          {
  if (!instance) {
    instance = new NtupleManager;
  }
  return instance;
}

InitNtuple()

void NtupleManager::InitNtuple

(

)

Definition at line 91 of file NtupleManager.cc.

References Chi2CalibratedParameters, Chi2Measurements, CloneFitParam, NCops, NDegreesOfFreedom, NDistancemeter, NDistancemeter1Dim, NFitParameters, NOptObjects, NSensor2D, and NTiltmeter.

Referenced by Fit::fitNextEvent().

                               {
  CloneFitParam->Clear();

  Chi2Measurements = 0.;
  Chi2CalibratedParameters = 0.;
  NDegreesOfFreedom = 0;
  NFitParameters = 0;
  NOptObjects = 0;
  NSensor2D = 0;
  NDistancemeter = 0;
  NDistancemeter1Dim = 0;
  NTiltmeter = 0;
  NCops = 0;
}

WriteNtuple()

void NtupleManager::WriteNtuple

(

)

Definition at line 112 of file NtupleManager.cc.

References theRootFile.

Referenced by Fit::startFit().

                                {
  theRootFile->Write();
  theRootFile->Close();
}

Member Data Documentation

Chi2CalibratedParameters

double NtupleManager::Chi2CalibratedParameters

private

Definition at line 64 of file NtupleManager.h.

Referenced by BookNtuple(), FillChi2(), and InitNtuple().

Chi2Measurements

double NtupleManager::Chi2Measurements

private

Definition at line 64 of file NtupleManager.h.

Referenced by BookNtuple(), FillChi2(), and InitNtuple().

CloneCopsMeas

TClonesArray* NtupleManager::CloneCopsMeas

private

Definition at line 58 of file NtupleManager.h.

Referenced by BookNtuple().

CloneDistancemeter1DimMeas

TClonesArray* NtupleManager::CloneDistancemeter1DimMeas

private

Definition at line 54 of file NtupleManager.h.

Referenced by BookNtuple().

CloneDistancemeterMeas

TClonesArray* NtupleManager::CloneDistancemeterMeas

private

Definition at line 52 of file NtupleManager.h.

Referenced by BookNtuple().

CloneFitParam

TClonesArray* NtupleManager::CloneFitParam

private

Definition at line 46 of file NtupleManager.h.

Referenced by BookNtuple(), and InitNtuple().

CloneOptObject

TClonesArray* NtupleManager::CloneOptObject

private

Definition at line 48 of file NtupleManager.h.

Referenced by BookNtuple(), FillFitParameters(), and FillMeasurements().

CloneSensor2DMeas

TClonesArray* NtupleManager::CloneSensor2DMeas

private

Definition at line 50 of file NtupleManager.h.

Referenced by BookNtuple().

CloneTiltmeterMeas

TClonesArray* NtupleManager::CloneTiltmeterMeas

private

Definition at line 56 of file NtupleManager.h.

Referenced by BookNtuple().

CocoaTree

TTree* NtupleManager::CocoaTree

private

Definition at line 42 of file NtupleManager.h.

Referenced by BookNtuple(), and FillNtupleTree().

CopsMeasA

CopsMeas* NtupleManager::CopsMeasA

private

Definition at line 59 of file NtupleManager.h.

Referenced by FillMeasurements().

Distancemeter1DimMeasA

Distancemeter1DimMeas* NtupleManager::Distancemeter1DimMeasA

private

Definition at line 55 of file NtupleManager.h.

Referenced by FillMeasurements().

DistancemeterMeasA

DistancemeterMeas* NtupleManager::DistancemeterMeasA

private

Definition at line 53 of file NtupleManager.h.

Referenced by FillMeasurements().

FitParamA

FitParam* NtupleManager::FitParamA

private

Definition at line 47 of file NtupleManager.h.

Referenced by FillFitParameters().

instance

NtupleManager * NtupleManager::instance = nullptr

staticprivate

Definition at line 36 of file NtupleManager.h.

Referenced by getInstance(), and production_tasks.Task::getname().

NCops

int NtupleManager::NCops

private

Definition at line 72 of file NtupleManager.h.

Referenced by BookNtuple(), FillMeasurements(), and InitNtuple().

NDegreesOfFreedom

int NtupleManager::NDegreesOfFreedom

private

Definition at line 65 of file NtupleManager.h.

Referenced by BookNtuple(), FillChi2(), and InitNtuple().

NDistancemeter

int NtupleManager::NDistancemeter

private

Definition at line 69 of file NtupleManager.h.

Referenced by BookNtuple(), FillMeasurements(), and InitNtuple().

NDistancemeter1Dim

int NtupleManager::NDistancemeter1Dim

private

Definition at line 70 of file NtupleManager.h.

Referenced by BookNtuple(), FillMeasurements(), and InitNtuple().

NFitParameters

int NtupleManager::NFitParameters

private

Definition at line 66 of file NtupleManager.h.

Referenced by BookNtuple(), FillFitParameters(), and InitNtuple().

NOptObjects

int NtupleManager::NOptObjects

private

Definition at line 67 of file NtupleManager.h.

Referenced by BookNtuple(), FillFitParameters(), FillMeasurements(), FillOptObjects(), and InitNtuple().

NSensor2D

int NtupleManager::NSensor2D

private

Definition at line 68 of file NtupleManager.h.

Referenced by BookNtuple(), FillMeasurements(), and InitNtuple().

NTiltmeter

int NtupleManager::NTiltmeter

private

Definition at line 71 of file NtupleManager.h.

Referenced by BookNtuple(), FillMeasurements(), and InitNtuple().

OptObjectA

OptObject* NtupleManager::OptObjectA

private

Definition at line 49 of file NtupleManager.h.

Referenced by FillOptObjects().

Sensor2DMeasA

Sensor2DMeas* NtupleManager::Sensor2DMeasA

private

Definition at line 51 of file NtupleManager.h.

Referenced by FillMeasurements().

theRootFile

TFile* NtupleManager::theRootFile

private

Definition at line 40 of file NtupleManager.h.

Referenced by BookNtuple(), and WriteNtuple().

TiltmeterMeasA

TiltmeterMeas* NtupleManager::TiltmeterMeasA

private

Definition at line 57 of file NtupleManager.h.

Referenced by FillMeasurements().