/data/refman/pasoursint/CMSSW_5_3_10/src/CalibTracker/SiPixelLorentzAngle/src/SiPixelLorentzAngleDB.cc

Go to the documentation of this file.

#include <memory>
#include <string>
#include <iostream>
#include <fstream>
#include "CondCore/DBOutputService/interface/PoolDBOutputService.h"
#include "CalibTracker/SiPixelLorentzAngle/interface/SiPixelLorentzAngleDB.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelLorentzAngle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "FWCore/ServiceRegistry/interface/Service.h"


#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementPoint.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementError.h"
#include "DataFormats/GeometrySurface/interface/GloballyPositioned.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"

using namespace std;
using namespace edm;

  //Constructor

SiPixelLorentzAngleDB::SiPixelLorentzAngleDB(edm::ParameterSet const& conf) : 
  conf_(conf){
        magneticField_ = conf_.getParameter<double>("magneticField");
        recordName_ = conf_.getUntrackedParameter<std::string>("record","SiPixelLorentzAngleRcd");
        bPixLorentzAnglePerTesla_ = (float)conf_.getParameter<double>("bPixLorentzAnglePerTesla");
        fPixLorentzAnglePerTesla_ = (float)conf_.getParameter<double>("fPixLorentzAnglePerTesla");
        useFile_ = conf_.getParameter<bool>("useFile");         
        fileName_ = conf_.getParameter<string>("fileName");

}

  //BeginJob

void SiPixelLorentzAngleDB::beginJob(){
  
}
// Virtual destructor needed.

SiPixelLorentzAngleDB::~SiPixelLorentzAngleDB() {  

}  

// Analyzer: Functions that gets called by framework every event

void SiPixelLorentzAngleDB::analyze(const edm::Event& e, const edm::EventSetup& es)
{

        SiPixelLorentzAngle* LorentzAngle = new SiPixelLorentzAngle();
           
        
        edm::ESHandle<TrackerGeometry> pDD;
        es.get<TrackerDigiGeometryRecord>().get( pDD );
        edm::LogInfo("SiPixelLorentzAngle") <<" There are "<<pDD->detUnits().size() <<" detectors"<<std::endl;
        
        for(TrackerGeometry::DetUnitContainer::const_iterator it = pDD->detUnits().begin(); it != pDD->detUnits().end(); it++){
    
                if( dynamic_cast<PixelGeomDetUnit*>((*it))!=0){
                        DetId detid=(*it)->geographicalId();
                        
                        // fill bpix values for LA 
                        if(detid.subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel)) {
                                
                                if(!useFile_){
                                        if ( ! LorentzAngle->putLorentzAngle(detid.rawId(),bPixLorentzAnglePerTesla_) )
                                        edm::LogError("SiPixelLorentzAngleDB")<<"[SiPixelLorentzAngleDB::analyze] detid already exists"<<std::endl;
                                } else {
                                        cout << "method for reading file not implemented yet" << endl;
                                }
                                
                        
                        // fill bpix values for LA 
                        } else if(detid.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap)) {
                                
                                if ( ! LorentzAngle->putLorentzAngle(detid.rawId(),fPixLorentzAnglePerTesla_) )
                                        edm::LogError("SiPixelLorentzAngleDB")<<"[SiPixelLorentzAngleDB::analyze] detid already exists"<<std::endl;
                        } else {
                                edm::LogError("SiPixelLorentzAngleDB")<<"[SiPixelLorentzAngleDB::analyze] detid is Pixel but neither bpix nor fpix"<<std::endl;
                        }
                
                }
                        
        }      
        

        edm::Service<cond::service::PoolDBOutputService> mydbservice;
        if( mydbservice.isAvailable() ){
                try{
                        if( mydbservice->isNewTagRequest(recordName_) ){
                                mydbservice->createNewIOV<SiPixelLorentzAngle>(LorentzAngle,
                                                                               mydbservice->beginOfTime(),
                                                                               mydbservice->endOfTime(),
                                                                               recordName_);
                        } else {
                                mydbservice->appendSinceTime<SiPixelLorentzAngle>(LorentzAngle,
                                                                                  mydbservice->currentTime(),
                                                                                  recordName_);
                        }
                }catch(const cond::Exception& er){
                        edm::LogError("SiPixelLorentzAngleDB")<<er.what()<<std::endl;
                }catch(const std::exception& er){
                        edm::LogError("SiPixelLorentzAngleDB")<<"caught std::exception "<<er.what()<<std::endl;
                }catch(...){
                        edm::LogError("SiPixelLorentzAngleDB")<<"Funny error"<<std::endl;
                }
        }else{
                edm::LogError("SiPixelLorentzAngleDB")<<"Service is unavailable"<<std::endl;
        }
   

}

void SiPixelLorentzAngleDB::endJob(){


}