#include <SiStripZeroSuppression.h>

Inheritance diagram for SiStripZeroSuppression:

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

	SiStripZeroSuppression (const edm::ParameterSet &)

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Types
typedef std::vector < edm::InputTag > ::const_iterator	tag_iterator_t

typedef token_v::const_iterator	token_iterator_t

typedef edm::EDGetTokenT < edm::DetSetVector < SiStripRawDigi > >	token_t

typedef std::vector< token_t >	token_v

Private Member Functions
void	CollectionMergedZeroSuppression (edm::Event &)

void	formatRawDigis (edm::DetSetVector< SiStripRawDigi >::const_iterator, edm::DetSet< SiStripRawDigi > &)

void	MergeCollectionsZeroSuppression (edm::Event &)

void	processRaw (const edm::InputTag &, const edm::DetSetVector< SiStripRawDigi > &)

void	StandardZeroSuppression (edm::Event &)

void	storeBaseline (uint32_t, const std::vector< std::pair< short, float > > &)

void	storeBaselinePoints (uint32_t)

void	storeCMN (uint32_t, const std::vector< std::pair< short, float > > &)

void	storeExtraOutput (uint32_t, int16_t)

Private Attributes
std::auto_ptr < SiStripRawProcessingAlgorithms >	algorithms

edm::EDGetTokenT < edm::DetSetVector < SiStripRawDigi > >	DigisToMergeVR

edm::EDGetTokenT < edm::DetSetVector < SiStripDigi > >	DigisToMergeZS

bool	fixCM

std::vector< edm::InputTag >	inputTags

token_v	inputTokens

bool	mergeCollections

std::vector< edm::DetSet < SiStripProcessedRawDigi > >	output_apvcm

std::vector< edm::DetSet < SiStripDigi > >	output_base

std::vector< edm::DetSet < SiStripRawDigi > >	output_base_raw

std::vector< edm::DetSet < SiStripProcessedRawDigi > >	output_baseline

std::vector< edm::DetSet < SiStripDigi > >	output_baseline_points

bool	produceBaselinePoints

bool	produceCalculatedBaseline

bool	produceRawDigis

bool	storeCM

bool	storeInZScollBadAPV

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 16 of file SiStripZeroSuppression.h.

Member Typedef Documentation

typedef std::vector<edm::InputTag>::const_iterator SiStripZeroSuppression::tag_iterator_t

private

Definition at line 40 of file SiStripZeroSuppression.h.

typedef token_v::const_iterator SiStripZeroSuppression::token_iterator_t

private

Definition at line 49 of file SiStripZeroSuppression.h.

typedef edm::EDGetTokenT< edm::DetSetVector<SiStripRawDigi> > SiStripZeroSuppression::token_t

private

Definition at line 47 of file SiStripZeroSuppression.h.

typedef std::vector<token_t> SiStripZeroSuppression::token_v

private

Definition at line 48 of file SiStripZeroSuppression.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 15 of file SiStripZeroSuppression.cc.

References DigisToMergeVR, DigisToMergeZS, fixCM, edm::ParameterSet::getParameter(), beam_dqm_sourceclient-live_cfg::inputTag, inputTags, inputTokens, mergeCollections, produceBaselinePoints, produceCalculatedBaseline, produceRawDigis, storeCM, storeInZScollBadAPV, GlobalPosition_Frontier_DevDB_cff::tag, and edm::vector_transform().

   : inputTags(conf.getParameter<std::vector<edm::InputTag> >("RawDigiProducersList")),
     algorithms(SiStripRawProcessingFactory::create(conf.getParameter<edm::ParameterSet>("Algorithms"))),
     storeCM(conf.getParameter<bool>("storeCM")),
     mergeCollections(conf.getParameter<bool>("mergeCollections"))
     
 {
    
     produceRawDigis = conf.getParameter<bool>("produceRawDigis");
     produceCalculatedBaseline = conf.getParameter<bool>("produceCalculatedBaseline");
     produceBaselinePoints = conf.getParameter<bool>("produceBaselinePoints");
     storeInZScollBadAPV = conf.getParameter<bool>("storeInZScollBadAPV");
     fixCM = conf.getParameter<bool>("fixCM");  
   
   if(mergeCollections){
     storeCM = false;
     produceRawDigis = false;
     DigisToMergeZS = consumes< edm::DetSetVector<SiStripDigi> >(conf.getParameter<edm::InputTag>("DigisToMergeZS"));
     DigisToMergeVR = consumes< edm::DetSetVector<SiStripRawDigi> >(conf.getParameter<edm::InputTag>("DigisToMergeVR"));
     produces< edm::DetSetVector<SiStripDigi> > ("ZeroSuppressed");
   }
   
   for(tag_iterator_t inputTag = inputTags.begin(); inputTag != inputTags.end(); ++inputTag ){
     produces< edm::DetSetVector<SiStripDigi> > (inputTag->instance());
     if(produceRawDigis) produces< edm::DetSetVector<SiStripRawDigi> > (inputTag->instance());
     if(produceCalculatedBaseline) produces< edm::DetSetVector<SiStripProcessedRawDigi> > ("BADAPVBASELINE"+inputTag->instance());
     if(produceBaselinePoints) produces< edm::DetSetVector<SiStripDigi> > ("BADAPVBASELINEPOINTS"+inputTag->instance());
     if(storeCM) produces< edm::DetSetVector<SiStripProcessedRawDigi> > ("APVCM"+inputTag->instance());
 //     tokens consumes<reco::BeamSpot>(
   } 
   inputTokens = edm::vector_transform( inputTags, [this](edm::InputTag const & tag) { return consumes< edm::DetSetVector<SiStripRawDigi> >(tag);} );
   
   
   
 }

Member Function Documentation

void SiStripZeroSuppression::CollectionMergedZeroSuppression ( edm::Event & e )

inlineprivate

Definition at line 415 of file SiStripZeroSuppression.cc.

References edm::Event::getByLabel(), beam_dqm_sourceclient-live_cfg::inputTag, inputTags, convertSQLitetoXML_cfg::output, output_base, processRaw(), and edm::Event::put().

                                                                               {
 
   for(tag_iterator_t inputTag = inputTags.begin(); inputTag != inputTags.end(); ++inputTag ) {
 
     edm::Handle< edm::DetSetVector<SiStripDigi> > inputdigi;
     edm::Handle< edm::DetSetVector<SiStripRawDigi> > inputraw;
     e.getByLabel(*inputTag,inputdigi);
     e.getByLabel(*inputTag,inputraw);
     
     std::vector<edm::DetSet<SiStripDigi> > outputdigi; 
     std::vector<edm::DetSet<SiStripRawDigi> > outputraw;  
        
     
     if (inputraw->size())   
       processRaw(*inputTag, *inputraw);
     
     
     for ( std::vector<edm::DetSet<SiStripDigi> >::const_iterator itinputdigi = inputdigi->begin(); itinputdigi !=inputdigi->end(); ++itinputdigi) {
       output_base.push_back(*itinputdigi);  
     }
     
     std::auto_ptr< edm::DetSetVector<SiStripDigi> > output(new edm::DetSetVector<SiStripDigi>(output_base) );
     e.put( output, inputTag->instance() );
     
   }
 
 }

void SiStripZeroSuppression::formatRawDigis	(	edm::DetSetVector< SiStripRawDigi >::const_iterator	rawDigis,
		edm::DetSet< SiStripRawDigi > &	outRawDigis
	)

inlineprivate

Definition at line 148 of file SiStripZeroSuppression.cc.

References algorithms, edm::DetSetVector< T >::begin(), edm::DetSetVector< T >::end(), and edm::DetSet< T >::push_back().

Referenced by processRaw().

                                                                          {
      
       const std::vector<bool>& apvf = algorithms->GetAPVFlags();
       edm::DetSet<SiStripRawDigi>::const_iterator itRawDigis = rawDigis->begin(); 
          
       uint32_t strip=0;
       for (; itRawDigis != rawDigis->end(); ++itRawDigis){
     int16_t APVn = strip/128;
         if(apvf[APVn]) outRawDigis.push_back(*itRawDigis); 
         else outRawDigis.push_back(SiStripRawDigi(0));
         ++strip;
        }
           
 }

void SiStripZeroSuppression::MergeCollectionsZeroSuppression ( edm::Event & e )

inlineprivate

Definition at line 265 of file SiStripZeroSuppression.cc.

References ecalMGPA::adc(), algorithms, edm::DetSet< T >::begin(), edm::HandleBase::clear(), KineDebug3::count(), gather_cfg::cout, DigisToMergeVR, DigisToMergeZS, edm::DetSet< T >::end(), edm::Event::getByToken(), convertSQLitetoXML_cfg::output, edm::Event::put(), and edm::DetSet< T >::size().

Referenced by produce().

                                                                               {
     
     std::cout<< "starting Merging" << std::endl;
     edm::Handle< edm::DetSetVector<SiStripDigi> > inputdigi;
     edm::Handle< edm::DetSetVector<SiStripRawDigi> > inputraw;
     e.getByToken(DigisToMergeZS,inputdigi);
     e.getByToken(DigisToMergeVR,inputraw);
     
     std::cout << inputdigi->size() << " " << inputraw->size() << std::endl;
     
     if (inputraw->size()){
         
         std::vector<edm::DetSet<SiStripDigi> > outputdigi; 
         outputdigi.clear();
         
         //std::cout << "copying the input ZS to the output ZS collection" << std::endl;
         for ( edm::DetSetVector<SiStripDigi>::const_iterator Digis = inputdigi->begin(); Digis != inputdigi->end(); ++Digis)  outputdigi.push_back(*Digis);
        
             
         std::cout << "looping over the raw data collection" << std::endl;
         for ( edm::DetSetVector<SiStripRawDigi>::const_iterator rawDigis = inputraw->begin(); rawDigis != inputraw->end(); ++rawDigis) {
            
             edm::DetSet<SiStripRawDigi>::const_iterator itRawDigis = rawDigis->begin();
             uint16_t nAPV = rawDigis->size()/128;
             uint32_t rawDetId = rawDigis->id;
           
             std::vector<bool> restoredAPV;
             restoredAPV.clear();
             restoredAPV.insert(restoredAPV.begin(), nAPV, false); 
           
           
             bool isModuleRestored = false;
             for( uint16_t strip =0; strip < rawDigis->size();++strip){
                 if(itRawDigis[strip].adc()!=0){
                   restoredAPV[strip/128] = true;       
                   isModuleRestored = true;
                 }
             }
  
             
             if(isModuleRestored){
                 std::cout << "apply the ZS to the raw data collection" << std::endl;
                 edm::DetSet<SiStripDigi> suppressedDigis(rawDetId);
                 std::vector<int16_t> processedRawDigis(rawDigis->size());
                                 algorithms->SuppressVirginRawData(*rawDigis, suppressedDigis);
                
                 if(suppressedDigis.size()){   
                     std::cout << "looking for the detId with the new ZS in the collection of the zero suppressed data" << std::endl; 
                     std::vector<edm::DetSet<SiStripDigi> >::iterator zsModule = outputdigi.begin();
                     //std::vector<edm::DetSet<SiStripDigi> >::iterator LastLowerIdDigis = zsModule;
                     
                     uint32_t zsDetId = zsModule->id;
                     bool isModuleInZscollection = false;
                     while((zsDetId <= rawDetId)&&(zsModule != outputdigi.end())&&(!isModuleInZscollection)){
                         //std::cout << rawDetId << " ==== " <<  zsDetId << std::endl;
                         if( zsDetId == rawDetId){
                             isModuleInZscollection = true;
                         }else{
                             ++zsModule;
                             zsDetId = zsModule->id;
                         }
                     }         
                     std::cout << "after the look " << rawDetId << " ==== " <<  zsDetId << std::endl;
                     std::cout << "exiting looking for the detId with the new ZS in the collection of the zero suppressed data" << std::endl; 
         
                     //creating the map containing the digis (in rawdigi format) merged
                     std::vector<uint16_t> MergedRawDigis;
                     MergedRawDigis.clear();
                     MergedRawDigis.insert(MergedRawDigis.begin(), nAPV*128, 0);
                     
                     uint32_t count=0; // to be removed...
                     
                     edm::DetSet<SiStripDigi> newDigiToIndert(rawDetId);
                     if(!isModuleInZscollection){
                       std::cout << "WE HAVE A PROBLEM, THE MODULE IS NTOT FOUND" << std::endl;
                       outputdigi.insert(zsModule, newDigiToIndert);
                         --zsModule; 
                       std::cout << "new module id -1 " << zsModule->id << std::endl;
                         ++zsModule;
                       std::cout << "new module id " << zsModule->id << std::endl;
                         ++zsModule;
                         std::cout << "new module id +1 " << zsModule->id << std::endl;
                         --zsModule;
                         
                     } else {
                         std::cout << "inserting only the digis for not restored APVs" << std::endl;
                         std::cout << "size : " << zsModule->size() << std::endl;
                         edm::DetSet<SiStripDigi>::iterator itZsModule = zsModule->begin(); 
                         for(; itZsModule != zsModule->end(); ++itZsModule){
                             uint16_t adc = itZsModule->adc();
                             uint16_t strip = itZsModule->strip();
                             if(!restoredAPV[strip/128]){
                                 MergedRawDigis[strip] = adc;
                                 ++count;
                                 std::cout << "original count: "<< count << " strip: " << strip << " adc: " << adc << std::endl;  
                             }
                         } 
                         
                     }
                                          
                     std::cout << "size of digis to keep: " << count << std::endl;               
                     std::cout << "inserting only the digis for the restored APVs" << std::endl;
                     std::cout << "size : " << suppressedDigis.size() << std::endl;
                     edm::DetSet<SiStripDigi>::iterator itSuppDigi = suppressedDigis.begin();
                     for(; itSuppDigi != suppressedDigis.end(); ++itSuppDigi){
                       uint16_t adc = itSuppDigi->adc();
                       uint16_t strip = itSuppDigi->strip();
                       if(restoredAPV[strip/128]){
                         MergedRawDigis[strip] = adc;
                         std::cout << "new suppressed strip: " << strip << " adc: " << adc << std::endl;
                       }
                     } 
                   
                   
                   
                     std::cout << "suppressing the raw digis" << std::endl;
                     zsModule->clear();
                     for(uint16_t strip=0; strip < MergedRawDigis.size(); ++strip){
                       uint16_t adc = MergedRawDigis[strip];
                       if(adc) zsModule->push_back(SiStripDigi(strip, adc));
                     }
                     std::cout << "size zsModule after the merging: " << zsModule->size() << std::endl;
                     if((count + suppressedDigis.size()) != zsModule->size()) std::cout << "WE HAVE A PROBLEM!!!! THE NUMBER OF DIGIS IS NOT RIGHT==============" << std::endl;
                     std::cout << "exiting suppressing the raw digis" << std::endl;
                 }//if new ZS digis size
             } //if module restored
         }//loop over raw data collection
         
         uint32_t oldid =0;
         for(edm::DetSetVector<SiStripDigi>::const_iterator dg = outputdigi.begin(); dg != outputdigi.end(); ++dg){
             uint32_t iddg = dg->id;
             if(iddg < oldid){
                 std::cout<< "NOT IN THE RIGHT ORGER" << std:: endl;
                 std::cout<< "======================="<< std:: endl;
             }
             oldid = iddg; 
         }
         
         
         std::cout << "write the output vector" << std::endl;
         std::auto_ptr< edm::DetSetVector<SiStripDigi> > output(new edm::DetSetVector<SiStripDigi>(outputdigi) );
         e.put( output, "ZeroSuppressed" );  
         
         
     }//if inputraw.size
 
 
    
 }

void SiStripZeroSuppression::processRaw	(	const edm::InputTag &	inputTag,
		const edm::DetSetVector< SiStripRawDigi > &	input
	)

inlineprivate

Definition at line 105 of file SiStripZeroSuppression.cc.

References algorithms, edm::DetSetVector< T >::begin(), edm::DetSetVector< T >::end(), Exception, formatRawDigis(), edm::InputTag::instance(), output_apvcm, output_base, output_base_raw, output_baseline, output_baseline_points, produceBaselinePoints, produceCalculatedBaseline, produceRawDigis, edm::DetSet< T >::size(), storeCM, storeExtraOutput(), and storeInZScollBadAPV.

Referenced by CollectionMergedZeroSuppression(), and StandardZeroSuppression().

                                                                                        {
 
   output_apvcm.clear();
   output_baseline.clear();
   output_baseline_points.clear();
   output_base.clear(); 
   output_base_raw.clear();
 
   if(storeCM) output_apvcm.reserve(16000);
   if(produceCalculatedBaseline) output_baseline.reserve(16000);
   if(produceBaselinePoints) output_baseline_points.reserve(16000);
   if(produceRawDigis) output_base_raw.reserve(16000);
   output_base.reserve(16000);    
  
   
   for ( edm::DetSetVector<SiStripRawDigi>::const_iterator 
       rawDigis = input.begin(); rawDigis != input.end(); ++rawDigis) {
         
       edm::DetSet<SiStripDigi> suppressedDigis(rawDigis->id);
       int16_t nAPVflagged = 0;
         
       if ( "ProcessedRaw" == inputTag.instance()) nAPVflagged = algorithms->SuppressProcessedRawData(*rawDigis, suppressedDigis);
       else if ( "VirginRaw" == inputTag.instance()) nAPVflagged = algorithms->SuppressVirginRawData(*rawDigis, suppressedDigis); 
       else     
       throw cms::Exception("Unknown input type") 
     << inputTag.instance() << " unknown.  SiStripZeroZuppression can only process types \"VirginRaw\" and \"ProcessedRaw\" ";
 
       //here storing the output
       this->storeExtraOutput(rawDigis->id, nAPVflagged);
       if (suppressedDigis.size() && (storeInZScollBadAPV || nAPVflagged ==0)) 
     output_base.push_back(suppressedDigis); 
          
       if (produceRawDigis && nAPVflagged > 0){  
     edm::DetSet<SiStripRawDigi> outRawDigis(rawDigis->id);
     this->formatRawDigis(rawDigis, outRawDigis);
     output_base_raw.push_back(outRawDigis);
       }
          
   }
   
 }

void SiStripZeroSuppression::produce	(	edm::Event &	e,
		const edm::EventSetup &	es
	)

virtual

Implements edm::stream::EDProducerBase.

Definition at line 52 of file SiStripZeroSuppression.cc.

References algorithms, mergeCollections, MergeCollectionsZeroSuppression(), and StandardZeroSuppression().

                                               {
     
   algorithms->initialize(es, e);
      
   if(mergeCollections){
     this->MergeCollectionsZeroSuppression(e);
   }else{ 
     this->StandardZeroSuppression(e);
   }
 }

void SiStripZeroSuppression::StandardZeroSuppression ( edm::Event & e )

inlineprivate

Definition at line 63 of file SiStripZeroSuppression.cc.

References edm::Event::getByToken(), input, beam_dqm_sourceclient-live_cfg::inputTag, inputTags, inputTokens, convertSQLitetoXML_cfg::output, output_apvcm, output_base, output_base_raw, output_baseline, output_baseline_points, processRaw(), produceBaselinePoints, produceCalculatedBaseline, produceRawDigis, edm::Event::put(), and storeCM.

Referenced by produce().

                                                                       {
 
   token_iterator_t inputToken = inputTokens.begin();
   for(tag_iterator_t inputTag = inputTags.begin(); inputTag != inputTags.end(); ++inputTag,++inputToken ) {
 
     edm::Handle< edm::DetSetVector<SiStripRawDigi> > input;
     e.getByToken(*inputToken,input);
 
     if (input->size())
       processRaw(*inputTag, *input);
     
       std::auto_ptr< edm::DetSetVector<SiStripDigi> > output(new edm::DetSetVector<SiStripDigi>(output_base) );
       e.put( output, inputTag->instance() );
         
       if(produceRawDigis){
     std::auto_ptr< edm::DetSetVector<SiStripRawDigi> > outputraw(new edm::DetSetVector<SiStripRawDigi>(output_base_raw) );
     e.put(outputraw, inputTag->instance() );
       }
     
       if(produceCalculatedBaseline){
     std::auto_ptr< edm::DetSetVector<SiStripProcessedRawDigi> > outputbaseline(new edm::DetSetVector<SiStripProcessedRawDigi>(output_baseline) );
     e.put(outputbaseline, "BADAPVBASELINE"+inputTag->instance() );
       }
   
       if(produceBaselinePoints){
     std::auto_ptr< edm::DetSetVector<SiStripDigi> > outputbaselinepoints(new edm::DetSetVector<SiStripDigi>(output_baseline_points) );
     e.put(outputbaselinepoints, "BADAPVBASELINEPOINTS"+inputTag->instance() );
       }
   
       if(storeCM){
     std::auto_ptr< edm::DetSetVector<SiStripProcessedRawDigi> > outputAPVCM(new edm::DetSetVector<SiStripProcessedRawDigi>(output_apvcm) );
     e.put( outputAPVCM,"APVCM"+inputTag->instance());
       }
     
   }
 }

void SiStripZeroSuppression::storeBaseline	(	uint32_t	id,
		const std::vector< std::pair< short, float > > &	vmedians
	)

inlineprivate

Definition at line 178 of file SiStripZeroSuppression.cc.

References algorithms, i, output_baseline, edm::DetSet< T >::push_back(), and edm::DetSet< T >::size().

Referenced by storeExtraOutput().

                                                                                                         {
   
   std::map< uint16_t, std::vector < int16_t> >& baselinemap = algorithms->GetBaselineMap(); 
 
   edm::DetSet<SiStripProcessedRawDigi> baselineDetSet(id);
   std::map< uint16_t, std::vector < int16_t> >::iterator itBaselineMap;
   
   for(size_t i=0;i<vmedians.size();++i){
     uint16_t APVn = vmedians[i].first;
     float median = vmedians[i].second;
     itBaselineMap = baselinemap.find(APVn);
     if(itBaselineMap==baselinemap.end()){
       for(size_t strip=0; strip < 128; ++strip)  baselineDetSet.push_back(SiStripProcessedRawDigi(median));
     } else {
        for(size_t strip=0; strip < 128; ++strip) baselineDetSet.push_back(SiStripProcessedRawDigi((itBaselineMap->second)[strip]));
     }
     
   }
   
   if(baselineDetSet.size())
     output_baseline.push_back(baselineDetSet);
   
 }

void SiStripZeroSuppression::storeBaselinePoints ( uint32_t id )

inlineprivate

Definition at line 203 of file SiStripZeroSuppression.cc.

References algorithms, output_baseline_points, edm::DetSet< T >::push_back(), and edm::DetSet< T >::size().

Referenced by storeExtraOutput().

                                                            {
 
     std::map< uint16_t, std::map< uint16_t, int16_t> >&  BasPointVec = algorithms->GetSmoothedPoints();
     std::map< uint16_t, std::map< uint16_t, int16_t> >::iterator itBasPointVect = BasPointVec.begin() ;
     std::map< uint16_t, int16_t>::iterator itBaselinePointMap;
     
     edm::DetSet<SiStripDigi> baspointDetSet(id);
           
     for(; itBasPointVect != BasPointVec.end(); ++itBasPointVect){
         uint16_t APVn= itBasPointVect->first;
         itBaselinePointMap =itBasPointVect->second.begin();
         for(;itBaselinePointMap != itBasPointVect->second.end(); ++itBaselinePointMap){
                   uint16_t bpstrip = (itBaselinePointMap->first) + APVn*128;
                   int16_t  bp = itBaselinePointMap->second;
                   baspointDetSet.push_back(SiStripDigi(bpstrip,bp+128));
                
           } 
       }    
 
     
     if(baspointDetSet.size())
     output_baseline_points.push_back(baspointDetSet);
   
 }

void SiStripZeroSuppression::storeCMN	(	uint32_t	id,
		const std::vector< std::pair< short, float > > &	vmedians
	)

inlineprivate

Definition at line 229 of file SiStripZeroSuppression.cc.

References algorithms, plotBeamSpotDB::first, fixCM, i, output_apvcm, edm::DetSet< T >::push_back(), edm::second(), and edm::DetSet< T >::size().

Referenced by storeExtraOutput().

                                                                                                    {
     
   edm::DetSet<SiStripProcessedRawDigi> apvDetSet(id);
   short apvNb=0;
   
   std::vector<bool> apvf;
   apvf.clear();
   apvf.insert(apvf.begin(), 6, false);
 
   if(fixCM){
     std::vector<bool>& apvFlagged = algorithms->GetAPVFlags();
     for(uint16_t it=0; it< apvFlagged.size(); ++it) apvf[it] = apvFlagged[it];
   }
 
   for(size_t i=0;i<vmedians.size();++i){
     if(vmedians[i].first>apvNb){
       for(int i=0;i<vmedians[i].first-apvNb;++i){
     apvDetSet.push_back(SiStripProcessedRawDigi(-999.));
     apvNb++;
       }
     }
 
     if(fixCM&&apvf[vmedians[i].first]){
       apvDetSet.push_back(SiStripProcessedRawDigi(-999.));
     }else{
       apvDetSet.push_back(SiStripProcessedRawDigi(vmedians[i].second));
     }
     apvNb++;
   }
    
   if(apvDetSet.size())
     output_apvcm.push_back(apvDetSet);
   
 }

void SiStripZeroSuppression::storeExtraOutput	(	uint32_t	id,
		int16_t	nAPVflagged
	)

inlineprivate

Definition at line 166 of file SiStripZeroSuppression.cc.

References algorithms, produceBaselinePoints, produceCalculatedBaseline, storeBaseline(), storeBaselinePoints(), storeCM, and storeCMN().

Referenced by processRaw().

                                                                              {
 
       const std::vector< std::pair<short,float> >& vmedians = algorithms->getAPVsCM();
       if(storeCM) this->storeCMN(id, vmedians);
       if(nAPVflagged > 0){
     if(produceCalculatedBaseline) this->storeBaseline(id, vmedians);
     if(produceBaselinePoints) this->storeBaselinePoints(id);
       }
 }