Inheritance diagram for sistrip::FEDEmulatorModule:

Public Member Functions
	FEDEmulatorModule (const edm::ParameterSet &)

	~FEDEmulatorModule ()

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Member Functions
virtual void	beginRun (edm::Run &, const edm::EventSetup &)

virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
std::auto_ptr < SiStripRawProcessingAlgorithms >	algorithms_
	object for zero-suppression More...

bool	byModule_

sistrip::FEDEmulator	fedEmulator_

edm::InputTag	spyReorderedDigisTag_

edm::InputTag	spyVirginRawDigisTag_

SpyUtilities	utility_

Static Private Attributes
static const char *	messageLabel_ = "SiStripFEDEmulatorModule"

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

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

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

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 58 of file SiStripFEDEmulatorModule.cc.

Constructor & Destructor Documentation

sistrip::FEDEmulatorModule::FEDEmulatorModule ( const edm::ParameterSet & iConfig )

explicit

Definition at line 99 of file SiStripFEDEmulatorModule.cc.

References byModule_, fedEmulator_, and sistrip::FEDEmulator::initialise().

                                                                     : 
     spyReorderedDigisTag_(iConfig.getParameter<edm::InputTag>("SpyReorderedDigisTag")),
     spyVirginRawDigisTag_(iConfig.getParameter<edm::InputTag>("SpyVirginRawDigisTag")),
     byModule_(iConfig.getParameter<bool>("ByModule")),
     algorithms_(SiStripRawProcessingFactory::create(iConfig.getParameter<edm::ParameterSet>("Algorithms"))) 
   {
 
     fedEmulator_.initialise(byModule_);
 
     if (!byModule_) { //if not by module
       //the medians will be produced by fed id/channel
       produces<std::map<uint32_t,std::vector<uint32_t> > >("Medians");
       produces<edm::DetSetVector<SiStripRawDigi> >("PedestalsOrdered");
       produces<edm::DetSetVector<SiStripProcessedRawDigi> >("NoisesOrdered");
       produces<edm::DetSetVector<SiStripRawDigi> >("PedSubtrDigisOrdered");
       produces<edm::DetSetVector<SiStripRawDigi> >("CMSubtrDigisOrdered");
     }
     else { //by module
       produces<edm::DetSetVector<SiStripRawDigi> >("ModulePedestals");
       produces<edm::DetSetVector<SiStripProcessedRawDigi> >("ModuleNoises");
       produces<edm::DetSetVector<SiStripRawDigi> >("PedSubtrModuleDigis");
       produces<std::map<uint32_t,std::vector<uint32_t> > >("ModuleMedians");
       produces<edm::DetSetVector<SiStripRawDigi> >("CMSubtrModuleDigis");
       produces<edm::DetSetVector<SiStripDigi> >("ZSModuleDigis");
     }//end of by module check
 
   }//end of FEDEmulatorModule constructor

sistrip::FEDEmulatorModule::~FEDEmulatorModule ( )

Definition at line 127 of file SiStripFEDEmulatorModule.cc.

128 {

129 }

Member Function Documentation

void sistrip::FEDEmulatorModule::beginRun	(	edm::Run &	,
		const edm::EventSetup &
	)

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 281 of file SiStripFEDEmulatorModule.cc.

   {
 
   }

void sistrip::FEDEmulatorModule::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 133 of file SiStripFEDEmulatorModule.cc.

References algorithms_, FedChannelConnection::apvPairNumber(), edm::DetSetVector< T >::begin(), byModule_, SiStripFedCabling::connection(), gather_cfg::cout, FedChannelConnection::detId(), edm::DetSetVector< T >::end(), fedEmulator_, sistrip::SpyUtilities::fedIndex(), edm::Event::getByLabel(), sistrip::SpyUtilities::getCabling(), sistrip::SpyUtilities::getNoiseHandle(), sistrip::SpyUtilities::getPedestalHandle(), sistrip::FEDEmulator::initialiseModule(), python.multivaluedict::map(), FedChannelConnection::nApvPairs(), edm::Handle< T >::product(), edm::Event::put(), sistrip::FEDEmulator::retrieveNoises(), sistrip::FEDEmulator::retrievePedestals(), edm::DetSet< T >::size(), edm::DetSetVector< T >::size(), spyReorderedDigisTag_, spyVirginRawDigisTag_, sistrip::STRIPS_PER_FEDCH, sistrip::FEDEmulator::subtractCM(), sistrip::FEDEmulator::subtractPedestals(), utility_, cms::Exception::what(), and sistrip::FEDEmulator::zeroSuppress().

Referenced by python.JSONExport.JsonExport::export(), and python.HTMLExport.HTMLExport::export().

   {
     //update cabling and pedestals
     const SiStripFedCabling* lCabling = utility_.getCabling( iSetup );
     edm::ESHandle<SiStripPedestals> lPedsHandle = utility_.getPedestalHandle(iSetup);
     edm::ESHandle<SiStripNoises> lNoiseHandle = utility_.getNoiseHandle(iSetup);
   
     //initialise the algorithms object for the zero suppression
     algorithms_->initialize(iSetup);
 
     //retrieve the digis
     edm::Handle<edm::DetSetVector<SiStripRawDigi> > lDigisHandle;
     try { //to get the digis from the event
       if (!byModule_) {
     iEvent.getByLabel(spyReorderedDigisTag_, lDigisHandle);
       }
       else { //digis supplied by module
     iEvent.getByLabel(spyVirginRawDigisTag_, lDigisHandle);
       }//end of by module check
     } catch (const cms::Exception& e) {
       std::cout << e.what() ;
       return;
     } //end of get digis try
 
     const edm::DetSetVector<SiStripRawDigi> * lInputDigis = lDigisHandle.product();
 
     unsigned int lNDigis = lInputDigis->size();
 
     //define output containers
     //reserve space, will push_back elements
     std::vector<edm::DetSetVector<SiStripRawDigi>::detset > pedsData;
     pedsData.reserve(lNDigis);
     std::vector<edm::DetSetVector<SiStripProcessedRawDigi>::detset > noiseData;
     noiseData.reserve(lNDigis);
     std::vector<edm::DetSetVector<SiStripRawDigi>::detset > pedSubtrData;
     pedSubtrData.reserve(lNDigis);
     std::vector<edm::DetSetVector<SiStripRawDigi>::detset > cmSubtrData;
     cmSubtrData.reserve(lNDigis);
     //zero suppressed contained - no fixed size (could be empty) so no need to reserve elements.
     std::vector< edm::DetSet<SiStripDigi> > zsData;
     
     //this is a map: no reserve/resize
     std::map<uint32_t,std::vector<uint32_t> > medsData;
 
     edm::DetSetVector<SiStripRawDigi>::const_iterator inputChannel = lInputDigis->begin();
 
     for ( ; inputChannel!=lInputDigis->end(); ++inputChannel){//loop on input channels
       uint32_t lDetId = inputChannel->detId(); //either fedIndex or detId
 
       pedsData.push_back(    edm::DetSetVector<SiStripRawDigi>::detset(lDetId) );
       noiseData.push_back(   edm::DetSetVector<SiStripProcessedRawDigi>::detset(lDetId) );
       pedSubtrData.push_back(edm::DetSetVector<SiStripRawDigi>::detset(lDetId) );
       cmSubtrData.push_back( edm::DetSetVector<SiStripRawDigi>::detset(lDetId) );
 
       unsigned int lNStrips = inputChannel->size();
 
       //define output digi containers
       std::vector<SiStripRawDigi>& pedsDetSetData = pedsData.back().data;
       pedsDetSetData.reserve(lNStrips);
       std::vector<SiStripProcessedRawDigi>& noiseDetSetData = noiseData.back().data;
       noiseDetSetData.reserve(lNStrips);
       std::vector<SiStripRawDigi>& pedSubtrDetSetData = pedSubtrData.back().data;
       pedSubtrDetSetData.reserve(lNStrips);
       std::vector<SiStripRawDigi>& cmSubtrDetSetData = cmSubtrData.back().data;
       cmSubtrDetSetData.reserve(lNStrips);
       //zero suppressed - slightly different procedure as not fixed size
       edm::DetSet<SiStripDigi> zsDetSetData(lDetId);
       
       //determine the number of APV pairs in the channel
       uint32_t lNPairs = static_cast<uint32_t>(lNStrips*1./sistrip::STRIPS_PER_FEDCH);
       uint32_t lPair = 0;
       
       std::vector<uint32_t> medsDetSetData;
       medsDetSetData.reserve(lNPairs*2); //2*number of pairs per module. If not by module, lNPairs = 1...
       
       if (!byModule_) { //the input is not stored by module
     //need to retrieve the proper detId from cabling
     uint16_t lFedId = 0;
     uint16_t lFedChannel = 0;
     sistrip::SpyUtilities::fedIndex(lDetId, lFedId, lFedChannel);
             
     const FedChannelConnection & lConnection = lCabling->connection(lFedId,lFedChannel);
     lDetId = lConnection.detId();
     lNPairs = lConnection.nApvPairs();
     lPair = lConnection.apvPairNumber();
       }//end of by module check
 
       fedEmulator_.initialiseModule(lDetId,lNPairs,lPair);
 
       //get the pedestal values
       //stored by module in the database
       fedEmulator_.retrievePedestals(lPedsHandle);
       fedEmulator_.retrieveNoises(lNoiseHandle);
       
       //last option: fill medians from these ped subtr data
       //if want something else, need to call a method to fill
       //the data member medians_ of the class fedEmulator.
       fedEmulator_.subtractPedestals(inputChannel,
                      pedsDetSetData,
                      noiseDetSetData,
                      pedSubtrDetSetData,
                      medsDetSetData,
                      true);
       
       fedEmulator_.subtractCM(pedSubtrDetSetData,cmSubtrDetSetData);
 
       //fill the median map
       medsData[inputChannel->detId()] = medsDetSetData;
 
       //zero suppress the digis
       fedEmulator_.zeroSuppress(cmSubtrDetSetData, zsDetSetData, algorithms_);
       if (zsDetSetData.size()) zsData.push_back( zsDetSetData );
       
     }//loop on input channels
 
 
     std::auto_ptr<edm::DetSetVector<SiStripRawDigi> > lPeds(new edm::DetSetVector<SiStripRawDigi>(pedsData,true));
     std::auto_ptr<edm::DetSetVector<SiStripProcessedRawDigi> > lNoises(new edm::DetSetVector<SiStripProcessedRawDigi>(noiseData,true));
 
     std::auto_ptr<edm::DetSetVector<SiStripRawDigi> > lOutputPedSubtr(new edm::DetSetVector<SiStripRawDigi>(pedSubtrData,true));
 
     std::auto_ptr<edm::DetSetVector<SiStripRawDigi> > lOutputCMSubtr(new edm::DetSetVector<SiStripRawDigi>(cmSubtrData,true));
 
     std::auto_ptr<std::map<uint32_t,std::vector<uint32_t> > > lMedians(new std::map<uint32_t,std::vector<uint32_t> >(medsData));
   
     //zero suppressed digis
     std::auto_ptr< edm::DetSetVector<SiStripDigi> > lOutputZS(new edm::DetSetVector<SiStripDigi>(zsData));
     
     if (!byModule_) {
       iEvent.put(lMedians,"Medians");
       iEvent.put(lPeds,"PedestalsOrdered");
       iEvent.put(lNoises,"NoisesOrdered");
       iEvent.put(lOutputPedSubtr,"PedSubtrDigisOrdered");
       iEvent.put(lOutputCMSubtr,"CMSubtrDigisOrdered");
     }
     else {
       iEvent.put(lPeds,"ModulePedestals");
       iEvent.put(lNoises,"ModuleNoises");
       iEvent.put(lOutputPedSubtr,"PedSubtrModuleDigis");
       iEvent.put(lMedians,"ModuleMedians");
       iEvent.put(lOutputCMSubtr,"CMSubtrModuleDigis");
       iEvent.put(lOutputZS,"ZSModuleDigis");
     }
 
   }//produce method