#include <SiPixelLorentzAngleDB.h>

Inheritance diagram for SiPixelLorentzAngleDB:

Public Member Functions
virtual void	analyze (const edm::Event &e, const edm::EventSetup &c)

virtual void	beginJob ()

virtual void	endJob ()

	SiPixelLorentzAngleDB (const edm::ParameterSet &conf)

virtual	~SiPixelLorentzAngleDB ()

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Types
typedef std::vector < edm::ParameterSet >	Parameters

Private Member Functions
unsigned int	HVgroup (unsigned int panel, unsigned int module)

Private Attributes
Parameters	BPixParameters_

edm::ParameterSet	conf_

std::vector< std::pair < uint32_t, float > >	detid_la

std::string	fileName_

Parameters	FPixParameters_

double	magneticField_

std::string	recordName_

bool	useFile_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 23 of file SiPixelLorentzAngleDB.h.

Member Typedef Documentation

typedef std::vector< edm::ParameterSet > SiPixelLorentzAngleDB::Parameters

private

Definition at line 51 of file SiPixelLorentzAngleDB.h.

Constructor & Destructor Documentation

SiPixelLorentzAngleDB::SiPixelLorentzAngleDB ( const edm::ParameterSet & conf )

explicit

Definition at line 26 of file SiPixelLorentzAngleDB.cc.

References BPixParameters_, conf_, fileName_, FPixParameters_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), magneticField_, recordName_, and useFile_.

                                                                         : 
   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");
 
         BPixParameters_ = conf_.getUntrackedParameter<Parameters>("BPixParameters");
         FPixParameters_ = conf_.getUntrackedParameter<Parameters>("FPixParameters");
 }

SiPixelLorentzAngleDB::~SiPixelLorentzAngleDB ( )

virtual

Definition at line 46 of file SiPixelLorentzAngleDB.cc.

                                               {  
 
 }  

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 52 of file SiPixelLorentzAngleDB.cc.

References cond::service::PoolDBOutputService::appendSinceTime(), cond::service::PoolDBOutputService::beginOfTime(), PXFDetId::blade(), BPixParameters_, gather_cfg::cout, cond::service::PoolDBOutputService::createNewIOV(), cond::service::PoolDBOutputService::currentTime(), cond::rpcobgas::detid, PXFDetId::disk(), cond::service::PoolDBOutputService::endOfTime(), cppFunctionSkipper::exception, FPixParameters_, edm::EventSetup::get(), HVgroup(), edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), PXBDetId::ladder(), PXBDetId::layer(), PXBDetId::module(), PXFDetId::module(), PXFDetId::panel(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, SiPixelLorentzAngle::putLorentzAngle(), DetId::rawId(), recordName_, PXFDetId::side(), DetId::subdetId(), useFile_, and cms::Exception::what().

 {
 
     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)) {
                 
                 PXBDetId pxdetid = PXBDetId(detid);
                 cout << " hp:barrel:" << "  layer=" << pxdetid.layer() << "  ladder=" << pxdetid.ladder() << "  module=" << pxdetid.module() << endl;
 
            if(!useFile_){
 /*hp
             if ( ! LorentzAngle->putLorentzAngle(detid.rawId(),bPixLorentzAnglePerTesla_) )
             edm::LogError("SiPixelLorentzAngleDB")<<"[SiPixelLorentzAngleDB::analyze] detid already exists"<<std::endl;
 */
 //hp
                 for(Parameters::iterator it = BPixParameters_.begin(); it != BPixParameters_.end(); ++it) {
 //                           cout << " PSet: " << *it << ", module = " << it->getParameter<unsigned int>("module") << endl;
                            if( it->getParameter<unsigned int>("module") == pxdetid.module() && it->getParameter<unsigned int>("layer") == pxdetid.layer() )
                            {
                               float lorentzangle = (float)it->getParameter<double>("angle");
                               LorentzAngle->putLorentzAngle(detid.rawId(),lorentzangle);
                            }
                         }
 
            } else {
 //          cout << "method for reading file not implemented yet" << endl;
             edm::LogError("SiPixelLorentzAngleDB")<<"[SiPixelLorentzAngleDB::analyze] method for reading file not implemented yet" << std::endl;
            }
             
            // fill fpix values for LA 
         } else if(detid.subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap)) {
                 
 //hp
                       PXFDetId pxdetid = PXFDetId(detid);
                       cout << " hp:endcap:" << "  side=" << pxdetid.side() << "  disk=" << pxdetid.disk() << "  blade=" << pxdetid.blade() << "  panel=" << pxdetid.panel() << "  module=" << pxdetid.module() << endl;
 /*hp
               if ( ! LorentzAngle->putLorentzAngle(detid.rawId(),fPixLorentzAnglePerTesla_) )  edm::LogError("SiPixelLorentzAngleDB")<<"[SiPixelLorentzAngleDB::analyze] detid already exists"<<std::endl;
 */
 
 //hp
                 for(Parameters::iterator it = FPixParameters_.begin(); it != FPixParameters_.end(); ++it) {
 //                           cout << " PSet: " << *it << ", module = " << it->getParameter<unsigned int>("module") << endl;
                            if( it->getParameter<unsigned int>("side") == pxdetid.side() && it->getParameter<unsigned int>("disk") == pxdetid.disk() && it->getParameter<unsigned int>("HVgroup") == HVgroup(pxdetid.panel(),pxdetid.module()) )
                            {
                               float lorentzangle = (float)it->getParameter<double>("angle");
                               LorentzAngle->putLorentzAngle(detid.rawId(),lorentzangle);
                            }
                         }
 
         } 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::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 41 of file SiPixelLorentzAngleDB.cc.

                                     {
   
 }

void SiPixelLorentzAngleDB::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 172 of file SiPixelLorentzAngleDB.cc.

                                   {
 
 }

unsigned int SiPixelLorentzAngleDB::HVgroup	(	unsigned int	panel,
		unsigned int	module
	)

private

Definition at line 150 of file SiPixelLorentzAngleDB.cc.

References gather_cfg::cout.

Referenced by analyze().

                                                                                   {
 
    if( 1 == panel && ( 1 == module || 2 == module ))  {
       return 1;
    }
    else if( 1 == panel && ( 3 == module || 4 == module ))  {
       return 2;
    }
    else if( 2 == panel && 1 == module )  {
       return 1;
    }
    else if( 2 == panel && ( 2 == module || 3 == module ))  {
       return 2;
    }
    else {
       cout << " *** error *** in SiPixelLorentzAngleDB::HVgroup(...), panel = " << panel << ", module = " << module << endl;
       return 0;
    }
    
 }