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

Detailed Description

Definition at line 17 of file FittedEntriesManager.h.

Constructor & Destructor Documentation

FittedEntriesManager::FittedEntriesManager ( )

inline

Definition at line 22 of file FittedEntriesManager.h.

Referenced by getInstance().

{ };

FittedEntriesManager::~FittedEntriesManager ( )

inline

Definition at line 23 of file FittedEntriesManager.h.

{ };

Member Function Documentation

void FittedEntriesManager::AddFittedEntriesSet

(

FittedEntriesSet *

fents

)

Definition at line 33 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 163 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 begin, gather_cfg::cout, ALIUtils::debug, FittedEntriesSet::FittedEntries(), FittedEntry::getSigma(), FittedEntry::getValue(), and mathSSE::sqrt().

{
  
  //---------- 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 126 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 29 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 44 of file FittedEntriesManager.cc.

References AddFittedEntriesSet(), begin, gather_cfg::cout, createFileName(), ALIUtils::debug, ALIUtils::dumpDimensions(), 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 = 0

staticprivate

Definition at line 37 of file FittedEntriesManager.h.

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

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

private

Definition at line 38 of file FittedEntriesManager.h.

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