#include <DataMixingSiStripRawWorker.h>

Public Member Functions
void	addSiStripPileups (const int bcr, edm::EventPrincipal *, unsigned int EventId)

void	addSiStripSignals (const edm::Event &e)

	DataMixingSiStripRawWorker ()

	DataMixingSiStripRawWorker (const edm::ParameterSet &ps)

void	putSiStrip (edm::Event &e)

virtual	~DataMixingSiStripRawWorker ()

Private Types
typedef std::vector< SiStripDigi >	OneDetectorMap

typedef std::map< uint32_t, OneDetectorMap >	SiGlobalIndex

Private Attributes
const edm::DetSetVector < SiStripDigi > *	digicollection_

std::string	label_

const edm::DetSetVector < SiStripRawDigi > *	rawdigicollection_

Selector *	sel_

SiGlobalIndex	SiHitStorage_

edm::InputTag	Sistripdigi_collectionSig_

std::string	SiStripDigiCollectionDM_

edm::InputTag	SistripLabelSig_

edm::InputTag	SiStripPileInputTag_

std::string	SiStripRawDigiSource_

edm::InputTag	SiStripRawInputTag_

Detailed Description

Definition at line 37 of file DataMixingSiStripRawWorker.h.

Member Typedef Documentation

typedef std::vector<SiStripDigi> edm::DataMixingSiStripRawWorker::OneDetectorMap

private

Definition at line 66 of file DataMixingSiStripRawWorker.h.

typedef std::map<uint32_t, OneDetectorMap> edm::DataMixingSiStripRawWorker::SiGlobalIndex

private

Definition at line 67 of file DataMixingSiStripRawWorker.h.

Constructor & Destructor Documentation

DataMixingSiStripRawWorker::DataMixingSiStripRawWorker ( )

Definition at line 26 of file DataMixingSiStripRawWorker.cc.

26 { sel_=0;}

edm::DataMixingSiStripRawWorker::sel_

Selector * sel_

Definition: DataMixingSiStripRawWorker.h:75

DataMixingSiStripRawWorker::DataMixingSiStripRawWorker ( const edm::ParameterSet & ps )

explicit

standard constructor

Definition at line 29 of file DataMixingSiStripRawWorker.cc.

References edm::ParameterSet::getParameter(), label_, sel_, SiHitStorage_, Sistripdigi_collectionSig_, SiStripDigiCollectionDM_, SistripLabelSig_, SiStripPileInputTag_, SiStripRawDigiSource_, and SiStripRawInputTag_.

                                                                                   : 
                                 label_(ps.getParameter<std::string>("Label"))
 
   {                                                         
 
     // get the subdetector names
     //    this->getSubdetectorNames();  //something like this may be useful to check what we are supposed to do...
 
     // create input selector
     if (label_.size()>0){
       sel_=new Selector( ModuleLabelSelector(label_));
     }
     else {
       sel_=new Selector( MatchAllSelector());
     }
 
     // declare the products to produce
 
     Sistripdigi_collectionSig_   = ps.getParameter<edm::InputTag>("SistripdigiCollectionSig");
     SistripLabelSig_   = ps.getParameter<edm::InputTag>("SistripLabelSig");
 
     SiStripPileInputTag_ = ps.getParameter<edm::InputTag>("SiStripPileInputTag");
     SiStripRawInputTag_ = ps.getParameter<edm::InputTag>("SiStripRawInputTag");
 
     SiStripDigiCollectionDM_  = ps.getParameter<std::string>("SiStripDigiCollectionDM");
 
     SiStripRawDigiSource_ = ps.getParameter<std::string>("SiStripRawDigiSource");
 
     // clear local storage for this event                                                                     
     SiHitStorage_.clear();
 
   }

DataMixingSiStripRawWorker::~DataMixingSiStripRawWorker ( )

virtual

Default destructor

Definition at line 64 of file DataMixingSiStripRawWorker.cc.

References sel_.

                                                           { 
     delete sel_;
     sel_=0;
   }  

Member Function Documentation

void DataMixingSiStripRawWorker::addSiStripPileups	(	const int	bcr,
		edm::EventPrincipal *	ep,
		unsigned int	EventId
	)

Definition at line 92 of file DataMixingSiStripRawWorker.cc.

References gather_cfg::cout, digicollection_, edm::EventPrincipal::id(), LogDebug, rawdigicollection_, SiStripPileInputTag_, SiStripRawDigiSource_, and SiStripRawInputTag_.

Referenced by edm::DataMixingModule::addPileups().

                                                                                                             {
     
     LogDebug("DataMixingSiStripRawWorker") << "\n===============> adding pileups from event  "
                        << ep->id() << " for bunchcrossing " << bcr;
 
     boost::shared_ptr<Wrapper<edm::DetSetVector<SiStripDigi> > const>    pSSD;
     boost::shared_ptr<Wrapper<edm::DetSetVector<SiStripRawDigi> > const> pSSRD;
     
     if (SiStripRawDigiSource_=="SIGNAL") {
       pSSD = getProductByTag<edm::DetSetVector<SiStripDigi> >(*ep, SiStripPileInputTag_);
       digicollection_ = const_cast< edm::DetSetVector<SiStripDigi> * >(pSSD->product());
     } else if (SiStripRawDigiSource_=="PILEUP") {
       pSSRD = getProductByTag<edm::DetSetVector<SiStripRawDigi> >(*ep, SiStripRawInputTag_ );
       rawdigicollection_ = const_cast< edm::DetSetVector<SiStripRawDigi> * >(pSSRD->product());
     } else {
       std::cout << "you shouldn't be here" << std::endl;
     }
 
   } // end of addSiStripPileups

void DataMixingSiStripRawWorker::addSiStripSignals ( const edm::Event & e )

Definition at line 71 of file DataMixingSiStripRawWorker.cc.

References gather_cfg::cout, digicollection_, edm::Event::getByLabel(), edm::InputTag::label(), edm::Handle< T >::product(), rawdigicollection_, Sistripdigi_collectionSig_, SistripLabelSig_, SiStripRawDigiSource_, and SiStripRawInputTag_.

Referenced by edm::DataMixingModule::addSignals().

                                                                       { 
 
     
     edm::Handle< edm::DetSetVector<SiStripDigi> >      hSSD;
     edm::Handle< edm::DetSetVector<SiStripRawDigi> >   hSSRD;
     
     if (SiStripRawDigiSource_=="SIGNAL") {
       e.getByLabel(SiStripRawInputTag_,hSSRD);
       rawdigicollection_ = hSSRD.product();
     } else if (SiStripRawDigiSource_=="PILEUP") {
       e.getByLabel(Sistripdigi_collectionSig_.label(),SistripLabelSig_.label(),hSSD);
       digicollection_ =  hSSD.product();
     } else {
       std::cout << "you shouldn't be here" << std::endl;
     }
     
 
   } // end of addSiStripSignals

void DataMixingSiStripRawWorker::putSiStrip ( edm::Event & e )

Definition at line 113 of file DataMixingSiStripRawWorker.cc.

References begin, edm::DetSetVector< T >::begin(), digicollection_, end, edm::DetSetVector< T >::end(), edm::DetSetVector< T >::insert(), LogDebug, edm::DetSet< T >::push_back(), edm::Event::put(), rawdigicollection_, SiHitStorage_, and SiStripDigiCollectionDM_.

Referenced by edm::DataMixingModule::put().

                                                          {
 
 
     //------------------
     //  (1) Fill a map from the Digi collection
     //
 
     // fill in maps of SiStripDigis
     OneDetectorMap LocalMap;
     
     //loop on all detsets (detectorIDs) inside the input collection
     edm::DetSetVector<SiStripDigi>::const_iterator DSViter=digicollection_->begin();
     for (; DSViter!=digicollection_->end();DSViter++){
       
 #ifdef DEBUG
       LogDebug("DataMixingSiStripRawWorker")  << "Processing DetID " << DSViter->id;
 #endif
       
       LocalMap.clear();
       LocalMap.reserve((DSViter->data).size());
       LocalMap.insert(LocalMap.end(),(DSViter->data).begin(),(DSViter->data).end());    
       
       SiHitStorage_.insert( SiGlobalIndex::value_type( DSViter->id, LocalMap ) );
     }
 
 
     //------------------
     //  (2) Loop over the input RawDigi collection and add the Digis from the map
     //
 
     // collection of RawDigis to put back in the event
     std::vector< edm::DetSet<SiStripRawDigi> > vSiStripRawDigi;
 
     //loop on all detsets (detectorIDs) inside the SiStripRawDigis collection
     edm::DetSetVector<SiStripRawDigi>::const_iterator rawDSViter=rawdigicollection_->begin();
     for (; rawDSViter!=rawdigicollection_->end();rawDSViter++){
 
       // Make empty collection with this detID
       edm::DetSet<SiStripRawDigi> SSRD(rawDSViter->id); 
 
       // find local map (if it exists) for this detector ID
       SiGlobalIndex::const_iterator itest;
       itest = SiHitStorage_.find(rawDSViter->id);
 
       // if detID already has digis in existing map, add them to rawdigis
       if(itest!=SiHitStorage_.end()) {  
 
 #ifdef DEBUG
     LogDebug("DataMixingSiStripRawWorker")  << "Pileups: Processing DetID " << rawDSViter->id;
 #endif
 
     // get the map from storage
     LocalMap = itest->second;
     OneDetectorMap::const_iterator iLocal  = LocalMap.begin();
 
     // loop on all strips in rawdigi detset
     int currentstrip=0;
     edm::DetSet<SiStripRawDigi>::const_iterator iRawDigi = rawDSViter->begin();
     while( iRawDigi != rawDSViter->end() ) {
 
       int ADCSum = iRawDigi->adc();
 
       // if current strip exists in map, add ADC values
       if(iLocal->strip() == currentstrip) {
         ADCSum += iLocal->adc();
         iLocal++;
       }
 
       // put ADC sum in DetSet and go to next strip
       SSRD.push_back( SiStripRawDigi(ADCSum) );
       iRawDigi++;
       currentstrip++;
 
     }
 
     // copy combined digi+rawdigi into rawdigi DetSetVector
     vSiStripRawDigi.push_back(SSRD);
 
       // otherwise, just copy the rawdigis from the background event to the output
       } else {
     vSiStripRawDigi.push_back(*rawDSViter);
       }
 
     }
 
 
     //------------------
     //  (3) Put the new RawDigi collection back into the event
     //
 
     // make new raw digi collection
     std::auto_ptr< edm::DetSetVector<SiStripRawDigi> > MySiStripRawDigis(new edm::DetSetVector<SiStripRawDigi>(vSiStripRawDigi) );
 
     // put collection
     e.put( MySiStripRawDigis, SiStripDigiCollectionDM_ );
 
     // clear local storage for this event
     SiHitStorage_.clear();
   }