#include <DTVDriftWriter.h>

Inheritance diagram for DTVDriftWriter:

Public Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup)

virtual void	beginRun (const edm::Run &run, const edm::EventSetup &setup)

	DTVDriftWriter (const edm::ParameterSet &pset)

virtual void	endJob ()

virtual	~DTVDriftWriter ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
edm::ESHandle< DTGeometry >	dtGeom_

std::string	granularity_

const DTMtime *	mTimeMap_

dtCalibration::DTVDriftBaseAlgo *	vDriftAlgo_

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

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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 24 of file DTVDriftWriter.h.

Constructor & Destructor Documentation

DTVDriftWriter::DTVDriftWriter ( const edm::ParameterSet & pset )

Definition at line 33 of file DTVDriftWriter.cc.

References edm::hlt::Exception, reco::get(), edm::ParameterSet::getParameter(), granularity_, and vDriftAlgo_.

                                                       :
   granularity_( pset.getUntrackedParameter<string>("calibGranularity","bySL") ) {
 
   LogVerbatim("Calibration") << "[DTVDriftWriter]Constructor called!";
 
   if(granularity_ != "bySL")
      throw cms::Exception("Configuration") << "[DTVDriftWriter] Check parameter calibGranularity: " << granularity_ << " option not available.";
 
   // Get the concrete algo from the factory
   string algoName = pset.getParameter<string>("vDriftAlgo");
   ParameterSet algoPSet = pset.getParameter<ParameterSet>("vDriftAlgoConfig");
   vDriftAlgo_ = DTVDriftPluginFactory::get()->create(algoName,algoPSet);
 }

DTVDriftWriter::~DTVDriftWriter ( )

virtual

Definition at line 47 of file DTVDriftWriter.cc.

References vDriftAlgo_.

                                {
   LogVerbatim("Calibration") << "[DTVDriftWriter]Destructor called!";
   delete vDriftAlgo_;
 }

Member Function Documentation

virtual void DTVDriftWriter::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

inlinevirtual

Implements edm::EDAnalyzer.

Definition at line 31 of file DTVDriftWriter.h.

31 {}

void DTVDriftWriter::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 52 of file DTVDriftWriter.cc.

References dtGeom_, edm::EventSetup::get(), mTimeMap_, dtCalibration::DTVDriftBaseAlgo::setES(), and vDriftAlgo_.

                                                                            {
   // Get the map of ttrig from the Setup
   ESHandle<DTMtime> mTime;
   setup.get<DTMtimeRcd>().get(mTime);
   mTimeMap_ = &*mTime;
 
   // Get geometry from Event Setup
   setup.get<MuonGeometryRecord>().get(dtGeom_);
   // Pass EventSetup to concrete implementation
   vDriftAlgo_->setES(setup); 
 }

void DTVDriftWriter::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 64 of file DTVDriftWriter.cc.

References DTVelocityUnits::cm_per_ns, dtCalibration::DTVDriftBaseAlgo::compute(), dtGeom_, alignCSCRings::e, cms::Exception::explainSelf(), DTMtime::get(), granularity_, mTimeMap_, GlobalPosition_Frontier_DevDB_cff::record, dtDQMClient_cfg::resolution, dtCalibration::DTVDriftData::resolution, DTMtime::set(), ntuplemaker::status, dtCalibration::DTVDriftData::vdrift, and vDriftAlgo_.

                             {
   // Create the object to be written to DB
   DTMtime* mTimeNewMap = new DTMtime();
 
   if(granularity_ == "bySL") {    
      // Get all the sls from the geometry
      const vector<DTSuperLayer*>& superLayers = dtGeom_->superLayers(); 
      vector<DTSuperLayer*>::const_iterator sl = superLayers.begin();
      vector<DTSuperLayer*>::const_iterator sl_end = superLayers.end();
      for(; sl != sl_end; ++sl){
         DTSuperLayerId slId = (*sl)->id();
         // Get original value from DB
         float vDrift = 0., resolution = 0.;
         // vdrift is cm/ns , resolution is cm
         int status = mTimeMap_->get(slId,vDrift,resolution,DTVelocityUnits::cm_per_ns);
 
         // Compute vDrift
         try{
            dtCalibration::DTVDriftData vDriftData = vDriftAlgo_->compute(slId);
            float vDriftNew = vDriftData.vdrift;
            float resolutionNew = vDriftData.resolution; 
            // vdrift is cm/ns , resolution is cm
            mTimeNewMap->set(slId,
                             vDriftNew,
                     resolutionNew,
                     DTVelocityUnits::cm_per_ns);
            LogVerbatim("Calibration") << "vDrift for: " << slId
                                       << " Mean " << vDriftNew
                                       << " Resolution " << resolutionNew;
         } catch(cms::Exception& e){
            LogError("Calibration") << e.explainSelf();
            // Go back to original value in case of error
            if(!status){  
               mTimeNewMap->set(slId,
                                vDrift,
                                resolution,
                                DTVelocityUnits::cm_per_ns);
               LogVerbatim("Calibration") << "Keep original vDrift for: " << slId
                                          << " Mean " << vDrift
                                          << " Resolution " << resolution;
            }
         }
      } // End of loop on superlayers
   }
 
   // Write the vDrift object to DB
   LogVerbatim("Calibration") << "[DTVDriftWriter]Writing vdrift object to DB!";
   string record = "DTMtimeRcd";
   DTCalibDBUtils::writeToDB<DTMtime>(record, mTimeNewMap);
 }