---

SurveyInputTrackerFromDB.cc

Go to the documentation of this file.

#include "Alignment/CommonAlignment/interface/SurveyDet.h"
#include "Alignment/TrackerAlignment/interface/AlignableTracker.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeomBuilderFromGeometricDet.h"

#include "Alignment/SurveyAnalysis/plugins/SurveyInputTrackerFromDB.h"

SurveyInputTrackerFromDB::SurveyInputTrackerFromDB(const edm::ParameterSet& cfg):
  textFileName( cfg.getParameter<std::string>("textFileName") )
{
}

void SurveyInputTrackerFromDB::beginJob(const edm::EventSetup& setup )
{

//  std::cout << "***************ENTERING BEGIN JOB******************" << std::endl;
        
  //Get map from textreader
  SurveyInputTextReader dataReader;
  dataReader.readFile( textFileName );
  uIdMap = dataReader.UniqueIdMap();

  edm::ESHandle<GeometricDet>  geom;
  setup.get<IdealGeometryRecord>().get(geom);    
  TrackerGeometry* tracker = TrackerGeomBuilderFromGeometricDet().build(&*geom);
  
  addComponent( new AlignableTracker( tracker ) );
  addSurveyInfo( detector() );
//  std::cout << "*************END BEGIN JOB***************" << std::endl;
}

void SurveyInputTrackerFromDB::addSurveyInfo(Alignable* ali)
{
  const align::Alignables& comp = ali->components();
        //std::cout << "in addSurveyInfo" << std::endl;
  unsigned int nComp = comp.size();
        //std::cout << "alOjbId: " << ali->alignableObjectId();
        //std::cout << " , detId: " << ali->geomDetId().rawId() << std::endl;
        
  for (unsigned int i = 0; i < nComp; ++i) addSurveyInfo(comp[i]);
        
  align::ErrorMatrix error;
        
  //if (ali->alignableObjectId() != AlignableObjectId::AlignableDetUnit){
                
  //    std::cout << "UniqueId: " << id.first << ", " << id.second << std::endl;
  SurveyInputTextReader::MapType::const_iterator it
    = uIdMap.find(std::make_pair(ali->id(), ali->alignableObjectId()));
  //    std::cout << "itID: " << it->first.first << ", " << it->first.second << std::endl;
  if (it != uIdMap.end()){
                
    const align::Scalars& parameters = (it)->second;
                
    //move the surface
    //displacement
    align::LocalVector lvector (parameters[0], parameters[1], parameters[2]);
    align::GlobalVector gvector = ali->surface().toGlobal(lvector);
    ali->move(gvector);
    //rotation
    Basic3DVector<align::Scalar> rot_aa(parameters[3], parameters[4], parameters[5]);
    align::RotationType rotation(rot_aa, rot_aa.mag());
    ali->rotateInLocalFrame(rotation);
                
    //sets the errors for the hierarchy level
    double* errorData = error.Array();
    for (unsigned int i = 0; i < 21; ++i){errorData[i] = parameters[i+6];}

    ali->setSurvey( new SurveyDet(ali->surface(), error*(1e-6)) );
  }
  else {
    error = ROOT::Math::SMatrixIdentity();
    ali->setSurvey( new SurveyDet(ali->surface(), error * 1e-6) );
  }
                
}
// Plug in to framework

#include "FWCore/Framework/interface/MakerMacros.h"


DEFINE_FWK_MODULE(SurveyInputTrackerFromDB);

---