|
|
|
#include <FittedEntriesManager.h>
Definition at line 17 of file FittedEntriesManager.h.
| FittedEntriesManager::FittedEntriesManager | ( | ) | [inline] |
| FittedEntriesManager::~FittedEntriesManager | ( | ) | [inline] |
Definition at line 23 of file FittedEntriesManager.h.
{ };
| FittedEntriesManager::FittedEntriesManager | ( | ) | [inline] |
Definition at line 22 of file FittedEntriesManager.h.
{ };
| FittedEntriesManager::~FittedEntriesManager | ( | ) | [inline] |
Definition at line 23 of file FittedEntriesManager.h.
{ };
| void FittedEntriesManager::AddFittedEntriesSet | ( | FittedEntriesSet * | fents | ) |
Definition at line 33 of file FittedEntriesManager.cc.
References theFittedEntriesSets.
Referenced by MakeHistos().
{
theFittedEntriesSets.push_back( fents );
}
| void FittedEntriesManager::AddFittedEntriesSet | ( | FittedEntriesSet * | 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;
}
| ALIstring FittedEntriesManager::createFileName | ( | const ALIstring & | optoName, |
| const ALIstring & | entryName | ||
| ) | [private] |
| void FittedEntriesManager::dumpEntriesSubstraction | ( | std::ofstream & | fout, |
| FittedEntriesSet & | fes, | ||
| ALIint | order1, | ||
| ALIint | order2 | ||
| ) | [private] |
| 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] |
| 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; }
| 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 MakeHistos().
{
if(!instance) {
instance = new FittedEntriesManager;
}
return instance;
}
| static FittedEntriesManager* FittedEntriesManager::getInstance | ( | ) | [static] |
| 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(), 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();
}
| void FittedEntriesManager::MakeHistos | ( | ) |
FittedEntriesManager * FittedEntriesManager::instance = 0 [static, private] |
Definition at line 37 of file FittedEntriesManager.h.
Referenced by getInstance().
std::vector< FittedEntriesSet * > FittedEntriesManager::theFittedEntriesSets [private] |
Definition at line 38 of file FittedEntriesManager.h.
Referenced by AddFittedEntriesSet(), getFittedEntriesSets(), and MakeHistos().
1.7.3