|
|
|
#include <CocoaDBMgr.h>
Public Member Functions | |
| bool | DumpCocoaResults () |
| ~CocoaDBMgr () | |
Static Public Member Functions | |
| static CocoaDBMgr * | getInstance () |
Private Member Functions | |
| std::pair< Alignments *, AlignmentErrors * > | BuildAlignments (bool bDT) |
| OpticalAlignments * | BuildOpticalAlignments () |
| CocoaDBMgr () | |
| AlignTransformError * | GetAlignInfoErrorFromOptO (OpticalObject *opto) |
| AlignTransform * | GetAlignInfoFromOptO (OpticalObject *opto) |
| double | GetEntryError (const Entry *entry1, const Entry *entry2) |
| double | GetEntryError (const Entry *entry) |
| OpticalAlignInfo | GetOptAlignInfoFromOptO (OpticalObject *opto) |
Static Private Attributes | |
| static CocoaDBMgr * | instance = 0 |
Definition at line 37 of file CocoaDBMgr.h.
| CocoaDBMgr::CocoaDBMgr | ( | ) | [private] |
| CocoaDBMgr::~CocoaDBMgr | ( | ) | [inline] |
Definition at line 45 of file CocoaDBMgr.h.
{ }
| std::pair< Alignments *, AlignmentErrors * > CocoaDBMgr::BuildAlignments | ( | bool | bDT | ) | [private] |
Definition at line 286 of file CocoaDBMgr.cc.
References gather_cfg::cout, ALIUtils::debug, GetAlignInfoErrorFromOptO(), GetAlignInfoFromOptO(), Alignments::m_align, AlignmentErrors::m_alignError, and Model::OptOList().
Referenced by DumpCocoaResults().
{
Alignments* alignments = new Alignments;
AlignmentErrors* alignmentErrors = new AlignmentErrors;
//read
static std::vector< OpticalObject* > optolist = Model::OptOList();
static std::vector< OpticalObject* >::const_iterator ite;
for(ite = optolist.begin(); ite != optolist.end(); ite++ ){
if( (*ite)->type() == "system" ) continue;
std::cout << "CocoaDBMgr::BuildAlignments getCmsswID " << (*ite) << std::endl;
std::cout << "CocoaDBMgr::BuildAlignments getCmsswID " << (*ite)->getCmsswID() << std::endl;
//check CMSSW ID
if( (*ite)->getCmsswID() > 0 ) { //put the numbers of DT or CSC objects
std::cout << " cal fill alignments " << std::endl;
alignments->m_align.push_back( *(GetAlignInfoFromOptO( *ite )));
std::cout << " fill alignments " << std::endl;
// AlignTransformError* err =
//GetAlignInfoErrorFromOptO( *ite );
alignmentErrors->m_alignError.push_back(*(GetAlignInfoErrorFromOptO( *ite )));
std::cout << "CocoaDBMgr::BuildAlignments add alignmentError " << alignmentErrors->m_alignError.size() << std::endl;
}
}
if(ALIUtils::debug >= 4) std::cout << "CocoaDBMgr::BuildAlignments end with n alignment " << alignments->m_align.size() << " n alignmentError " << alignmentErrors->m_alignError.size() << std::endl;
return std::pair< Alignments*,AlignmentErrors*>(alignments,alignmentErrors);
}
| OpticalAlignments * CocoaDBMgr::BuildOpticalAlignments | ( | ) | [private] |
Definition at line 265 of file CocoaDBMgr.cc.
References gather_cfg::cout, data, ALIUtils::debug, GetOptAlignInfoFromOptO(), OpticalAlignments::opticalAlignments_, and Model::OptOList().
Referenced by DumpCocoaResults().
{
OpticalAlignments* optalign= new OpticalAlignments;
static std::vector< OpticalObject* > optolist = Model::OptOList();
static std::vector< OpticalObject* >::const_iterator ite;
for(ite = optolist.begin(); ite != optolist.end(); ite++ ){
if( (*ite)->type() == "system" ) continue;
OpticalAlignInfo data = GetOptAlignInfoFromOptO( *ite );
optalign->opticalAlignments_.push_back(data);
if(ALIUtils::debug >= 5) {
std::cout << "@@@@ OPTALIGNINFO TO BE WRITTEN TO DB "
<< data
<< std::endl;
}
}
return optalign;
}
| bool CocoaDBMgr::DumpCocoaResults | ( | ) |
Definition at line 52 of file CocoaDBMgr.cc.
References BuildAlignments(), BuildOpticalAlignments(), gather_cfg::cout, ALIUtils::debug, Exception, getInstance(), GlobalOptionMgr::GlobalOptions(), edm::Service< T >::isAvailable(), Alignments::m_align, AlignmentErrors::m_alignError, Fit::nEvent, OpticalAlignments::opticalAlignments_, and TRUE.
{
edm::Service<cond::service::PoolDBOutputService> myDbService;
GlobalOptionMgr* gomgr = GlobalOptionMgr::getInstance();
int nrcd;
cond::Time_t appendTime = Fit::nEvent+1;
if(gomgr->GlobalOptions()["writeDBOptAlign"] > 0 ) {
//----- Build OpticalAlignments
OpticalAlignments* optalign = BuildOpticalAlignments();
//--- Dump OpticalAlignments
nrcd = optalign->opticalAlignments_.size();
if( !myDbService.isAvailable() ){
throw cms::Exception("CocoaDBMgr::DumpCocoaResults DB not available");
}
// try{
if ( myDbService->isNewTagRequest( "OpticalAlignmentsRcd" ) ) {
std::cout << " new OA to DB " << "begin " << myDbService->beginOfTime() << " current " << myDbService->currentTime() << " end " << myDbService->endOfTime() << std::endl;
myDbService->createNewIOV<OpticalAlignments>(optalign,
myDbService->beginOfTime(),
myDbService->endOfTime(),
// myDbService->endOfTime(),
"OpticalAlignmentsRcd");
} else {
std::cout << " old OA to DB " << " current " << myDbService->currentTime() << " end " << myDbService->endOfTime() << std::endl;
myDbService->appendSinceTime<OpticalAlignments>( optalign,
// myDbService->endOfTime(),
appendTime,
// myDbService->currentTime(),
"OpticalAlignmentsRcd");
}
/* }catch(const cond::Exception& er) {
std::cout<<er.what()<<std::endl;
}catch(const std::exception& er){
std::cout<<"caught std::exception "<<er.what()<<std::endl;
}catch(...){
std::cout<<"Funny error"<<std::endl;
} */
if(ALIUtils::debug >= 2) std::cout << "OpticalAlignmentsRcd WRITTEN TO DB : "<< nrcd << std::endl;
}
if( gomgr->GlobalOptions()["writeDBAlign"] > 0) {
// Build DT alignments and errors
std::pair< Alignments*,AlignmentErrors*> dtali = BuildAlignments(1);
Alignments* dt_Alignments = dtali.first;
AlignmentErrors* dt_AlignmentErrors = dtali.second;
// Dump DT alignments and errors
nrcd = dt_Alignments->m_align.size();
if ( myDbService->isNewTagRequest( "DTAlignmentRcd" ) ) {
myDbService->createNewIOV<Alignments>(&(*dt_Alignments),
myDbService->beginOfTime(),
myDbService->endOfTime(),
"DTAlignmentRcd");
} else {
myDbService->appendSinceTime<Alignments>( &(*dt_Alignments),
appendTime,
// myDbService->currentTime(),
"DTAlignmentRcd");
}
if(ALIUtils::debug >= 2) std::cout << "DTAlignmentRcd WRITTEN TO DB : "<< nrcd << std::endl;
nrcd = dt_AlignmentErrors->m_alignError.size();
if ( myDbService->isNewTagRequest( "DTAlignmentErrorRcd" ) ) {
myDbService->createNewIOV<AlignmentErrors>(&(*dt_AlignmentErrors),
myDbService->beginOfTime(),
myDbService->endOfTime(),
"DTAlignmentErrorRcd");
} else {
myDbService->appendSinceTime<AlignmentErrors>( &(*dt_AlignmentErrors),
appendTime,
"DTAlignmentErrorRcd");
}
if(ALIUtils::debug >= 2) std::cout << "DTAlignmentErrorRcd WRITTEN TO DB : "<< nrcd << std::endl;
// Build CSC alignments and errors
std::pair< Alignments*,AlignmentErrors*> cscali = BuildAlignments(0);
Alignments* csc_Alignments = cscali.first;
AlignmentErrors* csc_AlignmentErrors = cscali.second;
// Dump CSC alignments and errors
nrcd = csc_Alignments->m_align.size();
if ( myDbService->isNewTagRequest( "CSCAlignmentRcd" ) ) {
myDbService->createNewIOV<Alignments>(&(*csc_Alignments),
myDbService->beginOfTime(),
myDbService->endOfTime(),
"CSCAlignmentRcd");
} else {
myDbService->appendSinceTime<Alignments>( &(*csc_Alignments),
appendTime,
"CSCAlignmentRcd");
}
if(ALIUtils::debug >= 2) std::cout << "CSCAlignmentRcd WRITTEN TO DB : "<< nrcd << std::endl;
nrcd = csc_AlignmentErrors->m_alignError.size();
if ( myDbService->isNewTagRequest( "CSCAlignmentErrorRcd" ) ) {
myDbService->createNewIOV<AlignmentErrors>(&(*csc_AlignmentErrors),
myDbService->beginOfTime(),
myDbService->endOfTime(),
"CSCAlignmentErrorRcd");
} else {
myDbService->appendSinceTime<AlignmentErrors>( &(*csc_AlignmentErrors),
appendTime,
"CSCAlignmentErrorRcd");
}
if(ALIUtils::debug >= 2) std::cout << "CSCAlignmentErrorRcd WRITTEN TO DB : "<< nrcd << std::endl;
//? gives unreadable error??? std::cout << "@@@@ OPTICALALIGNMENTS WRITTEN TO DB " << *optalign << std::endl;
return TRUE;
}
return TRUE;
}
| AlignTransformError * CocoaDBMgr::GetAlignInfoErrorFromOptO | ( | OpticalObject * | opto | ) | [private] |
Definition at line 334 of file CocoaDBMgr.cc.
References asHepMatrix(), OpticalObject::CoordinateEntryList(), gather_cfg::cout, ALIUtils::debug, OpticalObject::getCmsswID(), GetEntryError(), GlobalErrorBase< T, ErrorWeightType >::matrix(), and OpticalObject::name().
Referenced by BuildAlignments().
{
if(ALIUtils::debug >= 3) std::cout << "@@@ CocoaDBMgr::GetAlignInfoErrorFromOptO " << opto->name() << std::endl;
align::ID cmsswID = opto->getCmsswID();
GlobalError gerr(1.,
0.,
1.,
0.,
0.,
1.);
//double(dx*dx), 0., double(dy*dy), 0., 0., double(dz*dz) ) ;
CLHEP::HepSymMatrix errms = asHepMatrix(gerr.matrix());
AlignTransformError* alignError = new AlignTransformError( errms, cmsswID );
return alignError;
CLHEP::HepMatrix errm(3,3);
const std::vector< Entry* > theCoordinateEntryVector = opto->CoordinateEntryList();
std::cout << "@@@ CocoaDBMgr::GetAlignInfoFromOptOfill errm " << opto->name() << std::endl;
errm(0,0) = GetEntryError( theCoordinateEntryVector[0] ) / 100.; // in cm
errm(1,1) = GetEntryError( theCoordinateEntryVector[1] ) / 100.; // in cm
errm(2,2) = GetEntryError( theCoordinateEntryVector[2] ) / 100.; // in cm
errm(0,1) = GetEntryError( theCoordinateEntryVector[0], theCoordinateEntryVector[1] ) / 100.; // in cm
errm(0,2) = GetEntryError( theCoordinateEntryVector[0], theCoordinateEntryVector[2] ) / 100.; // in cm
errm(1,2) = GetEntryError( theCoordinateEntryVector[1], theCoordinateEntryVector[2] ) / 100.; // in cm
// errm(1,0) = errm(0,1);
// errm(2,0) = errm(0,2);
// errm(2,1) = errm(1,2);
std::cout << "@@@ CocoaDBMgr::GetAlignInfoFromOptO errm filled" << opto->name() << std::endl;
// CLHEP::HepSymMatrix errms(3);
// errms.assign(errm);
std::cout << "@@@ CocoaDBMgr::GetAlignInfoFromOptO errms filled " << opto->name() << std::endl;
// AlignTransformError* alignError = new AlignTransformError( errms, cmsswID );
// AlignTransformError* alignError = 0;
std::cout << alignError << "@@@ CocoaDBMgr::GetAlignInfoFromOptO error built " << opto->name() << std::endl;
//t return alignError;
return (AlignTransformError*)(0);
}
| AlignTransform * CocoaDBMgr::GetAlignInfoFromOptO | ( | OpticalObject * | opto | ) | [private] |
Definition at line 316 of file CocoaDBMgr.cc.
References OpticalObject::centreGlob(), gather_cfg::cout, ALIUtils::debug, OpticalObject::getCmsswID(), OpticalObject::name(), OpticalObject::rmGlob(), and makeMuonMisalignmentScenario::rot.
Referenced by BuildAlignments().
{
if(ALIUtils::debug >= 3) std::cout << "@@@ CocoaDBMgr::GetAlignInfoFromOptO " << opto->name() << std::endl;
AlignTransform::Translation trans = opto->centreGlob();
AlignTransform::Rotation rot = opto->rmGlob();
align::ID cmsswID = opto->getCmsswID();
std::cout << "@@@ CocoaDBMgr::GetAlignInfoFromOptO buildalign" << opto->name() << std::endl;
AlignTransform* align = new AlignTransform( trans, rot, cmsswID );
std::cout << "@@@ CocoaDBMgr::GetAlignInfoFromOptO alig built " << opto->name() << std::endl;
return align;
// return dd;
}
Definition at line 252 of file CocoaDBMgr.cc.
References Entry::fitPos(), Fit::GetAtWAMatrix(), GetEntryError(), Entry::quality(), and mathSSE::sqrt().
{
if( entry1 == entry2 ) return GetEntryError( entry1 );
if( entry1->quality() > 0 && entry2->quality() > 0 ) {
return sqrt(Fit::GetAtWAMatrix()->Mat()->me[entry1->fitPos()][entry2->fitPos()]);
} else { //entries not fitted, correlation is 0
return 0.;
}
}
| double CocoaDBMgr::GetEntryError | ( | const Entry * | entry | ) | [private] |
Definition at line 241 of file CocoaDBMgr.cc.
References Entry::fitPos(), Fit::GetAtWAMatrix(), Entry::quality(), Entry::sigma(), and mathSSE::sqrt().
Referenced by GetAlignInfoErrorFromOptO(), GetEntryError(), and GetOptAlignInfoFromOptO().
| CocoaDBMgr * CocoaDBMgr::getInstance | ( | ) | [static] |
Definition at line 38 of file CocoaDBMgr.cc.
References CocoaDBMgr(), and instance.
Referenced by DumpCocoaResults().
{
if(!instance) {
instance = new CocoaDBMgr;
}
return instance;
}
| OpticalAlignInfo CocoaDBMgr::GetOptAlignInfoFromOptO | ( | OpticalObject * | opto | ) | [private] |
Definition at line 176 of file CocoaDBMgr.cc.
References OpticalAlignInfo::angx_, OpticalAlignInfo::angy_, OpticalAlignInfo::angz_, OpticalObject::centreGlob(), OpticalObject::CoordinateEntryList(), gather_cfg::cout, data, OpticalAlignParam::dim_type_, OpticalAlignParam::error_, OpticalAlignInfo::extraEntries_, OpticalObject::ExtraEntryList(), Fit::GetAtWAMatrix(), OpticalObject::getCmsswID(), GetEntryError(), OpticalObject::getLocalRotationAngles(), OpticalAlignInfo::ID_, M_PI, MatrixMeschach::Mat(), OpticalObject::name(), OpticalAlignInfo::name_, OpticalAlignParam::name_, OpticalObject::parent(), OpticalAlignParam::quality_, OpticalObject::rmGlob(), OpticalObject::type(), OpticalAlignInfo::type_, OpticalAlignParam::value_, OpticalAlignInfo::x_, OpticalAlignInfo::y_, and OpticalAlignInfo::z_.
Referenced by BuildOpticalAlignments().
{
std::cout << " CocoaDBMgr::GetOptAlignInfoFromOptO " << opto->name() << std::endl;
OpticalAlignInfo data;
data.ID_=opto->getCmsswID();
data.type_=opto->type();
data.name_=opto->name();
//----- Centre in local coordinates
CLHEP::Hep3Vector centreLocal = opto->centreGlob() - opto->parent()->centreGlob();
CLHEP::HepRotation parentRmGlobInv = inverseOf( opto->parent()->rmGlob() );
centreLocal = parentRmGlobInv * centreLocal;
const std::vector< Entry* > theCoordinateEntryVector = opto->CoordinateEntryList();
std::cout << " CocoaDBMgr::GetOptAlignInfoFromOptO starting coord " <<std::endl;
data.x_.value_= centreLocal.x() / 100.; // in cm
std::cout << " matrix " << Fit::GetAtWAMatrix() << std::endl;
std::cout << " matrix " << Fit::GetAtWAMatrix()->Mat() << " " << theCoordinateEntryVector[0]->fitPos() << std::endl;
data.x_.error_= GetEntryError( theCoordinateEntryVector[0] ) / 100.; // in cm
data.y_.value_= centreLocal.y() / 100.; // in cm
std::cout << " matrix " << Fit::GetAtWAMatrix()->Mat() << " " << theCoordinateEntryVector[1]->fitPos() << std::endl;
data.y_.error_= GetEntryError( theCoordinateEntryVector[1] ) / 100.; // in cm
data.z_.value_= centreLocal.z() / 100.; // in cm
std::cout << " matrix " << Fit::GetAtWAMatrix()->Mat() << " " << theCoordinateEntryVector[2]->fitPos() << std::endl;
data.z_.error_= GetEntryError( theCoordinateEntryVector[2] ) / 100.; // in cm
//----- angles in local coordinates
std::vector<double> anglocal = opto->getLocalRotationAngles( theCoordinateEntryVector );
data.angx_.value_= anglocal[0] *180./M_PI; // in deg
std::cout << " matrix " << Fit::GetAtWAMatrix()->Mat() << theCoordinateEntryVector[3]->fitPos() << std::endl;
data.angx_.error_= GetEntryError( theCoordinateEntryVector[3] ) * 180./M_PI; // in deg;
data.angy_.value_= anglocal[1] * 180./M_PI; // in deg
std::cout << " matrix " << Fit::GetAtWAMatrix()->Mat() << theCoordinateEntryVector[4]->fitPos() << std::endl;
data.angy_.error_= GetEntryError( theCoordinateEntryVector[4] ) * 180./M_PI; // in deg;;
data.angz_.value_= anglocal[2] *180./M_PI; // in deg
std::cout << " matrix " << Fit::GetAtWAMatrix()->Mat() << theCoordinateEntryVector[5]->fitPos() << std::endl;
data.angz_.error_= GetEntryError( theCoordinateEntryVector[5] ) * 180./M_PI; // in deg;
const std::vector< Entry* > theExtraEntryVector = opto->ExtraEntryList(); std::cout << " CocoaDBMgr::GetOptAlignInfoFromOptO starting entry " << std::endl;
std::vector< Entry* >::const_iterator ite;
for( ite = theExtraEntryVector.begin(); ite != theExtraEntryVector.end(); ite++ ) {
OpticalAlignParam extraEntry;
extraEntry.name_ = (*ite)->name();
extraEntry.dim_type_ = (*ite)->type();
extraEntry.value_ = (*ite)->value();
extraEntry.error_ = (*ite)->sigma();
extraEntry.quality_ = (*ite)->quality();
data.extraEntries_.push_back( extraEntry );
std::cout << " CocoaDBMgr::GetOptAlignInfoFromOptO done extra entry " << extraEntry.name_ << std::endl;
}
return data;
}
CocoaDBMgr * CocoaDBMgr::instance = 0 [static, private] |
Definition at line 62 of file CocoaDBMgr.h.
Referenced by getInstance().
1.7.3