FittedEntriesManager Class Reference

#include <FittedEntriesManager.h>

Public Member Functions
void	AddFittedEntriesSet (FittedEntriesSet *fents)
	FittedEntriesManager ()
std::vector< FittedEntriesSet * >	getFittedEntriesSets () const
void	MakeHistos ()
	~FittedEntriesManager ()

Static Public Member Functions
static FittedEntriesManager *	getInstance ()

Private Member Functions
ALIstring	createFileName (const ALIstring &optoName, const ALIstring &entryName)
void	dumpEntriesSubstraction (std::ofstream &fout, FittedEntriesSet &fes, ALIint order1, ALIint order2)
void	GetDifferentBetweenLasers ()

Private Attributes
std::vector< FittedEntriesSet * >	theFittedEntriesSets

Static Private Attributes
static FittedEntriesManager *	instance = nullptr

Detailed Description

Definition at line 16 of file FittedEntriesManager.h.

Constructor & Destructor Documentation

FittedEntriesManager::FittedEntriesManager ( )

inline

Definition at line 19 of file FittedEntriesManager.h.

Referenced by getInstance().

{};

FittedEntriesManager::~FittedEntriesManager ( )

inline

Definition at line 20 of file FittedEntriesManager.h.

{};

Member Function Documentation

void FittedEntriesManager::AddFittedEntriesSet

(

FittedEntriesSet *

fents

)

Definition at line 32 of file FittedEntriesManager.cc.

References theFittedEntriesSets.

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

{ theFittedEntriesSets.push_back(fents); }

ALIstring FittedEntriesManager::createFileName ( const ALIstring &  optoName, const ALIstring &  entryName  )

private

Definition at line 164 of file FittedEntriesManager.cc.

References gather_cfg::cout, ALIUtils::debug, lut2db_cfg::filename, and pickleFileParser::slash.

Referenced by MakeHistos().

                                                                                                    {
  //  std::cout << "in createFileName " << optoName << " " << entryName << std::endl;
  ALIstring filename;
  //-  std::cout << "o" << optoName << " e " << entryName << std::endl;
  /*  ALIint last_slash =  optoName.rfind('/');
  ALIint size = optoName.length();
  //-   std::cout << last_slash << " " << size << "optoname " << optoName << std::endl;

  filename = optoName.substr( last_slash+1, size );
  */
  //----- substitute '/' by '.' in opto name
  filename = optoName.substr(2, optoName.size());  // skip the first 's/'
  ALIint slash = -1;
  for (;;) {
    slash = filename.find('/', slash + 1);
    if (slash == -1)
      break;
    filename[slash] = '.';
  }

  //----- Check if there is a ' ' in entry (should not happen now)
  ALIint space = entryName.rfind(' ');
  filename.append(".");
  ALIstring en = entryName;
  if (space != -1)
    en[space] = '_';

  //----- Merge opto and entry names
  // now it is not used as filename   filename.append( en + ".out");
  filename.append(en);
  if (ALIUtils::debug >= 3)
    std::cout << "filename " << filename << std::endl;

  return filename;
}

void FittedEntriesManager::dumpEntriesSubstraction ( std::ofstream &  fout, FittedEntriesSet &  fes, ALIint  order1, ALIint  order2  )

private

Definition at line 204 of file FittedEntriesManager.cc.

References SplitLinear::begin, gather_cfg::cout, ALIUtils::debug, FittedEntriesSet::FittedEntries(), FittedEntry::getSigma(), FittedEntry::getValue(), mathSSE::sqrt(), and hgcalPerformanceValidation::val.

                                                                  {
  //---------- Found order of entry1 and entry2 in FittedEntriesSet fes
  // (the order will be the same for every FittedEntriesSet
  //std::vector< FittedEntriesSet* >::const_iterator vfescite = theFittedEntriesSets.begin();
  /* std::vector< FittedEntry* >::const_iterator vfecite;
  ALIint order1, order2;
  ALIint jj=0;
  for( vfecite = (fes.FittedEntries()).begin(); vfecite != (fes.FittedEntries()).end(); vfecite++) {
    ALIstring entryTemp = (*vfecite)->OptOName() + "/" + (*vfecite)->EntryName();
    if(entryTemp == entry1) order1 = jj;
    if(entryTemp == entry2) order2 = jj;
    jj++;
  }
  */

  FittedEntry* fe1 = *((fes.FittedEntries()).begin() + order1);
  FittedEntry* fe2 = *((fes.FittedEntries()).begin() + order2);
  //-------- Substract values of entry1 and entry2 (with errors)
  ALIdouble val1 = fe1->getValue();
  ALIdouble val2 = fe2->getValue();
  ALIdouble sig1 = fe1->getSigma();
  ALIdouble sig2 = fe2->getSigma();
  ALIdouble val = val1 - val2;
  ALIdouble sig = sqrt(sig1 * sig1 + sig2 * sig2);
  //-  std::cout << "CHECK " << val1 << " "<< val2 << " "<< sig1 << " "<< sig2 << std::endl;
  fout << std::setprecision(6) << std::setw(8) << val << " +- " << sig << "  ";
  if (ALIUtils::debug >= 3)
    std::cout << " FITTEDENTRY:" << std::setprecision(5) << std::setw(8) << val << " +- " << sig << std::endl;
}

void FittedEntriesManager::GetDifferentBetweenLasers ( )

private

Definition at line 129 of file FittedEntriesManager.cc.

References gather_cfg::cout, ALIUtils::debug, M_PI, Model::OptOList(), and LightRay::startLightRay().

Referenced by MakeHistos().

                                                     {
  std::vector<OpticalObject*> optoList = Model::OptOList();
  std::vector<OpticalObject*>::const_iterator ite;
  std::map<ALIstring, LightRay*> lrays;

  for (ite = optoList.begin(); ite != optoList.end(); ++ite) {
    if ((*ite)->type() == "laser") {
      LightRay* lightray = new LightRay;
      lightray->startLightRay(*ite);
      lrays[(*ite)->parent()->name()] = lightray;
    }
  }

  std::map<ALIstring, LightRay*>::const_iterator lite1, lite2;
  for (lite1 = lrays.begin(); lite1 != lrays.end(); ++lite1) {
    lite2 = lite1;
    ++lite2;
    for (; lite2 != lrays.end(); ++lite2) {
      if (lite1 == lite2)
        continue;
      CLHEP::Hep3Vector dirdiff = ((*lite1).second->direction() - (*lite2).second->direction());
      if (ALIUtils::debug >= 0) {
        std::cout << "LASER DIFF " << (*lite1).first << " & " << (*lite2).first << " " << dirdiff.mag() * 180. / M_PI
                  << "o " << dirdiff.mag() << " rad " << dirdiff << std::endl;

        (*lite1).second->dumpData(ALIstring(" laser ") + (*lite1).first);
        (*lite2).second->dumpData(ALIstring(" laser ") + (*lite2).first);
      }
    }
  }
}

std::vector<FittedEntriesSet*> FittedEntriesManager::getFittedEntriesSets ( ) const

inline

Definition at line 26 of file FittedEntriesManager.h.

References theFittedEntriesSets.

{ return theFittedEntriesSets; }

FittedEntriesManager * FittedEntriesManager::getInstance ( )

static

Definition at line 22 of file FittedEntriesManager.cc.

References FittedEntriesManager(), and instance.

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

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

void FittedEntriesManager::MakeHistos

(

)

Definition at line 38 of file FittedEntriesManager.cc.

References AddFittedEntriesSet(), SplitLinear::begin, gather_cfg::cout, createFileName(), ALIUtils::debug, ALIUtils::dumpDimensions(), dataset::end, lut2db_cfg::filename, groupFilesInBlocks::fout, GetDifferentBetweenLasers(), GlobalOptionMgr::getInstance(), GlobalOptionMgr::GlobalOptions(), cuy::ii, findQualityFiles::jj, findQualityFiles::size, and theFittedEntriesSets.

Referenced by Fit::startFit().

                                      {
  //----------- Get
  //----------- Loop entries
  //-  vfescite = theFittedEntriesSets.begin();
  //  std::vector< FittedEntry* > vfe = theFittedEntriesSets.begin()->FittedEntries();
  //  std::vector< FittedEntry* >::const_iterator vfecite2;
  //-- Number of fitted entries (equal for all Fitted Entries Sets )
  if (ALIUtils::debug >= 5)
    std::cout << "No sets2 " << theFittedEntriesSets.size() << " No ent "
              << ((*(theFittedEntriesSets.begin()))->FittedEntries()).size() << std::endl;
  std::ofstream fout;
  std::ofstream fout2;
  fout.open("fittedEntries.out");
  fout2.open("longFittedEntries.out");
  //---------- Dump dimensions
  ALIUtils::dumpDimensions(fout);
  ALIUtils::dumpDimensions(fout2);

  AddFittedEntriesSet(new FittedEntriesSet(theFittedEntriesSets));  //add a new set that averages all the others

  //---------- Loop sets of entries
  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
      fout << "  ";
      ALIint ii = 0;
      for (vfecite = ((*vfescite)->FittedEntries()).begin(); vfecite != ((*vfescite)->FittedEntries()).end();
           ++vfecite) {
        ALIstring filename = createFileName((*vfecite)->getOptOName(), (*vfecite)->getEntryName());
        fout << ii << ": " << std::setw(13) << filename << " ";
        ;
        if (ALIUtils::debug >= 3)
          std::cout << ii << ": " << std::setw(13) << filename << " = " << (*vfecite)->getName() << std::endl;
        if (ALIUtils::debug >= 3)
          std::cout << filename << " ";
        if (ALIUtils::debug >= 3)
          std::cout << "OPENING FITTED ENTRIES file " << filename << std::endl;
        ii++;
      }
      //      fout << std::setw(17) << "2:-4:";
      fout << std::endl;
      if (ALIUtils::debug >= 3)
        std::cout << std::endl;
    }
    //----- Dump entry set number
    fout << jj << " ";
    fout2 << jj << " ";
    //----- Dump measurements date
    GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
    if (gomgr->GlobalOptions()["DumpDateInFittedEntries"] >= 1) {
      fout << (*vfescite)->getDate() << " " << (*vfescite)->getTime() << " ";
      fout2 << (*vfescite)->getDate() << " " << (*vfescite)->getTime() << " ";
    }

    ALIint ii = 0;
    for (vfecite = ((*vfescite)->FittedEntries()).begin(); vfecite != ((*vfescite)->FittedEntries()).end(); ++vfecite) {
      //      std::cout << ii << *vfescite << " FITTEDENTRY: "   << vfecite << " " <<*vfecite << " " << (*vfecite)->Value() << std::endl;
      //      if( ii == 2 || ii == 4 ) {
      fout << std::setprecision(8) << std::setw(10) << (*vfecite)->getValue() << " " << (*vfecite)->getSigma() << "  ";
      //- fout << std::setw(9) << std::setprecision(6) << (*vfecite)->getValue()  << " +- " << (*vfecite)->getSigma() << "  ";
      //      }
      if (ALIUtils::debug >= 3)
        std::cout << " FITTEDENTRY:" << std::setprecision(5) << std::setw(8) << (*vfecite)->getValue() << " +- "
                  << (*vfecite)->getSigma() << std::endl;

      ALIstring filename = createFileName((*vfecite)->getOptOName(), (*vfecite)->getEntryName());
      fout2 << std::setprecision(8) << std::setw(10) << filename << " " << (*vfecite)->getValue() << " "
            << (*vfecite)->getSigma() << "  ";
      ii++;
    }
    //    dumpEntriesSubstraction( fout, *(*vfescite), 2, 4);
    fout << std::endl;
    fout2 << std::endl;
    if (ALIUtils::debug >= 3)
      std::cout << std::endl;
    jj++;
  }
  fout.close();
  fout2.close();

  GetDifferentBetweenLasers();
}

Member Data Documentation

FittedEntriesManager * FittedEntriesManager::instance = nullptr

staticprivate

Definition at line 33 of file FittedEntriesManager.h.

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

std::vector<FittedEntriesSet*> FittedEntriesManager::theFittedEntriesSets

private

Definition at line 34 of file FittedEntriesManager.h.

Referenced by AddFittedEntriesSet(), getFittedEntriesSets(), and MakeHistos().