StripCompactDigiSimLinksProducer Class Reference

Inheritance diagram for StripCompactDigiSimLinksProducer:

Public Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &) override
	StripCompactDigiSimLinksProducer (const edm::ParameterSet &iConfig)
	~StripCompactDigiSimLinksProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	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
	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
uint32_t	maxHoleSize_
edm::InputTag	src_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
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)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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 26 of file StripCompactDigiSimLinksProducer.cc.

Constructor & Destructor Documentation

StripCompactDigiSimLinksProducer::StripCompactDigiSimLinksProducer

(

const edm::ParameterSet &

iConfig

)

Definition at line 38 of file StripCompactDigiSimLinksProducer.cc.

                                                                                                 :
    src_(iConfig.getParameter<edm::InputTag>("src")),
    maxHoleSize_(iConfig.getParameter<uint32_t>("maxHoleSize"))
{
    produces<StripCompactDigiSimLinks>();
}

StripCompactDigiSimLinksProducer::~StripCompactDigiSimLinksProducer

(

)

Definition at line 45 of file StripCompactDigiSimLinksProducer.cc.

{
}

Member Function Documentation

void StripCompactDigiSimLinksProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &    )

overridevirtual

Implements edm::EDProducer.

Definition at line 50 of file StripCompactDigiSimLinksProducer.cc.

References edm::DetSet< T >::begin(), dtNoiseDBValidation_cfg::cerr, DEBUG, DEFINE_FWK_MODULE, edm::DetSet< T >::detId(), edm::DetSet< T >::end(), spr::find(), edm::Event::getByLabel(), combine::key, maxHoleSize_, convertSQLitetoXML_cfg::output, edm::Event::put(), edm::DetSet< T >::size(), findQualityFiles::size, alcazmumu_cfi::src, src_, and swap().

{
    using namespace edm;
    Handle<DetSetVector<StripDigiSimLink> > src;
    iEvent.getByLabel(src_, src);

    StripCompactDigiSimLinks::Filler output;

    int previousStrip;     // previous strip with at least one link (might not be the strip of the previous link there are overlapping clusters)
    int previousLinkStrip; // strip of the previous link (can be the same as the one of this link if there are overlapping clusters)
    std::vector<StripCompactDigiSimLinks::key_type> thisStripSignals;      // particles on this strip
    std::vector<StripCompactDigiSimLinks::key_type> previousStripSignals;  // particles on the previous strip

    foreach(const DetSet<StripDigiSimLink> &det, *src) {
        DEBUG(std::cerr << "\n\nProcessing detset " << det.detId() << ", size = " << det.size() << std::endl;)
        previousStrip     = -2; // previous strip with at least one link (might not be the strip of the previous link there are overlapping clusters)
        previousLinkStrip = -2; // strip of the previous link (can be the same as the one of this link if there are overlapping clusters)
        thisStripSignals.clear();
        previousStripSignals.clear();
        for (DetSet<StripDigiSimLink>::const_iterator it = det.begin(), ed = det.end(); it != ed; ++it) {
            DEBUG(std::cerr << "  processing digiSimLink on strip " << it->channel() << " left by particle " << it->SimTrackId() << ", event " << it->eventId().rawId() << std::endl;)
            if (int(it->channel()) != previousLinkStrip) { 
                previousStrip = previousLinkStrip;
                DEBUG(std::cerr << "   strip changed!" << std::endl;)
                swap(thisStripSignals, previousStripSignals); 
                thisStripSignals.clear();
            }
            DEBUG(std::cerr << "    previous strip " << previousStrip << ", previous link strip " << previousLinkStrip << std::endl;)
            //DEBUG(std::cerr << "    on previous strip: "; foreach(StripCompactDigiSimLinks::key_type k, previousStripSignals) { std::cerr << "(ev " << k.first.rawId() << ", id " << k.second << ") "; } std::cerr << std::endl;)
            //DEBUG(std::cerr << "    on this strip: ";     foreach(StripCompactDigiSimLinks::key_type k,     thisStripSignals) { std::cerr << "(ev " << k.first.rawId() << ", id " << k.second << ") "; } std::cerr << std::endl;)
            StripCompactDigiSimLinks::key_type key(it->eventId(), it->SimTrackId());
            bool alreadyClusterized = false;
            if (int(it->channel()) == previousStrip+1) {
                DEBUG(std::cerr << "  on next strip" << std::endl;)
                if (std::find(previousStripSignals.begin(), previousStripSignals.end(), key) != previousStripSignals.end()) {
                    alreadyClusterized = true;
                    DEBUG(std::cerr << "    and part of previous cluster" << std::endl;)
                }
            }
            if (!alreadyClusterized) {
                DEBUG(std::cerr << "   clusterize!" << std::endl;)
                unsigned int size = 1;
                int myLastStrip = it->channel(); // last known strip with signals from this particle
                for (DetSet<StripDigiSimLink>::const_iterator it2 = it+1; it2 < ed; ++it2) {
                    DEBUG(std::cerr << "     digiSimLink on strip " << it2->channel() << " left by particle " << it2->SimTrackId() << ", event " << it2->eventId().rawId() << std::endl;)
                    if ((it2->channel() - myLastStrip) > maxHoleSize_+1) {
                        DEBUG(std::cerr << "        found hole of size " << (it2->channel() - myLastStrip) << ", stopping." << std::endl;)
                        break;
                    }
                    if ((it2->eventId() == key.first) && (it2->SimTrackId() == key.second)) {
                        size++;
                        DEBUG(std::cerr << "        extending cluster, now size = " << size << std::endl;)
                        myLastStrip = it2->channel();
                    }
                }
                output.insert(key, StripCompactDigiSimLinks::HitRecord(det.detId(), it->channel(), size));
            }
            if (int(it->channel()) != previousLinkStrip) {
                previousLinkStrip = it->channel();
            }
            thisStripSignals.push_back(key);
            DEBUG(std::cerr << "    ending state " << previousStrip << ", previous link strip " << previousLinkStrip << std::endl;)
            //DEBUG(std::cerr << "    on previous strip: "; foreach(StripCompactDigiSimLinks::key_type k, previousStripSignals) { std::cerr << "(ev " << k.first.rawId() << ", id " << k.second << ") "; } std::cerr << std::endl;)
            //DEBUG(std::cerr << "    on this strip: ";     foreach(StripCompactDigiSimLinks::key_type k,     thisStripSignals) { std::cerr << "(ev " << k.first.rawId() << ", id " << k.second << ") "; } std::cerr << std::endl;)
            DEBUG(std::cerr << std::endl;)
        }
    }
   
    std::auto_ptr< StripCompactDigiSimLinks > ptr(new StripCompactDigiSimLinks(output));
    iEvent.put(ptr);
}

Member Data Documentation

uint32_t StripCompactDigiSimLinksProducer::maxHoleSize_

private

Definition at line 35 of file StripCompactDigiSimLinksProducer.cc.

Referenced by produce().

edm::InputTag StripCompactDigiSimLinksProducer::src_

private

Definition at line 34 of file StripCompactDigiSimLinksProducer.cc.

Referenced by produce().