#include <UserCode/CSCDigiValidator/src/CSCDigiValidator.cc>

Inheritance diagram for CSCDigiValidator:

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

	~CSCDigiValidator ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

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

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

virtual bool	filter (edm::Event &, const edm::EventSetup &)

std::vector< CSCComparatorDigi >	relabelCompDigis (const CSCChamberMap *m, CSCDetId _id, std::vector< CSCComparatorDigi >::const_iterator b, std::vector< CSCComparatorDigi >::const_iterator e)

std::vector< CSCStripDigi >	relabelStripDigis (const CSCChamberMap *, CSCDetId, std::vector< CSCStripDigi >::const_iterator, std::vector< CSCStripDigi >::const_iterator)

std::vector< CSCStripDigi >	sanitizeStripDigis (std::vector< CSCStripDigi >::const_iterator, std::vector< CSCStripDigi >::const_iterator)

std::vector< CSCWireDigi >	sanitizeWireDigis (std::vector< CSCWireDigi >::const_iterator, std::vector< CSCWireDigi >::const_iterator)

std::vector< CSCComparatorDigi >	zeroSupCompDigis (std::vector< CSCComparatorDigi >::const_iterator, std::vector< CSCComparatorDigi >::const_iterator)

std::vector< CSCStripDigi >	zeroSupStripDigis (std::vector< CSCStripDigi >::const_iterator, std::vector< CSCStripDigi >::const_iterator)

Private Attributes
edm::EDGetTokenT < CSCALCTDigiCollection >	al1_token

edm::EDGetTokenT < CSCALCTDigiCollection >	al2_token

edm::InputTag	alct1

edm::InputTag	alct2

edm::EDGetTokenT < CSCComparatorDigiCollection >	cd1_token

edm::EDGetTokenT < CSCComparatorDigiCollection >	cd2_token

edm::EDGetTokenT < CSCCLCTDigiCollection >	cl1_token

edm::EDGetTokenT < CSCCLCTDigiCollection >	cl2_token

edm::InputTag	clct1

edm::InputTag	clct2

edm::EDGetTokenT < CSCCorrelatedLCTDigiCollection >	co1_token

edm::EDGetTokenT < CSCCorrelatedLCTDigiCollection >	co2_token

edm::InputTag	comp1

edm::InputTag	comp2

edm::InputTag	csctf1

edm::InputTag	csctf2

edm::InputTag	csctfstubs1

edm::InputTag	csctfstubs2

edm::InputTag	lct1

edm::InputTag	lct2

bool	reorderStrips

edm::EDGetTokenT < CSCStripDigiCollection >	sd1_token

edm::EDGetTokenT < CSCStripDigiCollection >	sd2_token

edm::InputTag	strip1

edm::InputTag	strip2

edm::EDGetTokenT < L1CSCTrackCollection >	tr1_token

edm::EDGetTokenT < L1CSCTrackCollection >	tr2_token

edm::EDGetTokenT < CSCTriggerContainer < csctf::TrackStub > >	ts1_token

edm::EDGetTokenT < CSCTriggerContainer < csctf::TrackStub > >	ts2_token

edm::EDGetTokenT < CSCWireDigiCollection >	wd1_token

edm::EDGetTokenT < CSCWireDigiCollection >	wd2_token

edm::InputTag	wire1

edm::InputTag	wire2

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

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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 35 of file CSCDigiValidator.h.

Constructor & Destructor Documentation

CSCDigiValidator::CSCDigiValidator ( const edm::ParameterSet & iConfig )

explicit

Definition at line 26 of file CSCDigiValidator.cc.

References al1_token, al2_token, cd1_token, cd2_token, cl1_token, cl2_token, co1_token, co2_token, edm::ParameterSet::getParameter(), sd1_token, sd2_token, tr1_token, tr2_token, ts1_token, ts2_token, wd1_token, and wd2_token.

 {
   wd1_token = consumes<CSCWireDigiCollection>( iConfig.getParameter<edm::InputTag>("inputWire") );
   wd2_token = consumes<CSCWireDigiCollection>( iConfig.getParameter<edm::InputTag>("repackWire") );
   sd1_token = consumes<CSCStripDigiCollection>( iConfig.getParameter<edm::InputTag>("inputStrip") );
   sd2_token = consumes<CSCStripDigiCollection>( iConfig.getParameter<edm::InputTag>("inputStrip") );
   cd1_token = consumes<CSCComparatorDigiCollection>( iConfig.getParameter<edm::InputTag>("inputComp") );
   cd2_token = consumes<CSCComparatorDigiCollection>( iConfig.getParameter<edm::InputTag>("RepackComp") );
   al1_token = consumes<CSCALCTDigiCollection>( iConfig.getParameter<edm::InputTag>("inputALCT") );
   al2_token = consumes<CSCALCTDigiCollection>( iConfig.getParameter<edm::InputTag>("repackALCT") );
   cl1_token = consumes<CSCCLCTDigiCollection>( iConfig.getParameter<edm::InputTag>("inputCLCT") );
   cl2_token = consumes<CSCCLCTDigiCollection>( iConfig.getParameter<edm::InputTag>("repackCLCT") );
   co1_token = consumes<CSCCorrelatedLCTDigiCollection>( iConfig.getParameter<edm::InputTag>("inputCorrLCT") );
   co2_token = consumes<CSCCorrelatedLCTDigiCollection>( iConfig.getParameter<edm::InputTag>("repackCorrLCT") );
   tr1_token = consumes<L1CSCTrackCollection>( iConfig.getParameter<edm::InputTag>("inputCSCTF") );
   tr2_token = consumes<L1CSCTrackCollection>( iConfig.getParameter<edm::InputTag>("repackCSCTF") );
   ts1_token = consumes<CSCTriggerContainer<csctf::TrackStub> >( iConfig.getParameter<edm::InputTag>("inputCSCTFStubs") );
   ts2_token = consumes<CSCTriggerContainer<csctf::TrackStub> >( iConfig.getParameter<edm::InputTag>("repackCSCTFStubs") );
 
   //  reorderStrips(iConfig.getUntrackedParameter<bool>("applyStripReordering",true))
 
 }

CSCDigiValidator::~CSCDigiValidator ( )

Definition at line 49 of file CSCDigiValidator.cc.

50 {

51 }

Member Function Documentation

void CSCDigiValidator::beginJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 988 of file CSCDigiValidator.cc.

989 {

990 }

void CSCDigiValidator::endJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 993 of file CSCDigiValidator.cc.

994 {

995 }

bool CSCDigiValidator::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDFilter.

Definition at line 54 of file CSCDigiValidator.cc.

References a, al1_token, al2_token, b, trackerHits::c, GetRecoTauVFromDQM_MC_cff::cl, cl1_token, cl2_token, co1_token, co2_token, gather_cfg::cout, CSCDetId, cuy::cv, dt, alignCSCRings::e, edm::EventSetup::get(), edm::Event::getByToken(), i, gen::k, python.multivaluedict::map(), max(), bookConverter::min, EcalCondDBWriter_cfi::pedestal, relabelCompDigis(), relabelStripDigis(), alignCSCRings::s, sanitizeWireDigis(), sd1_token, sd2_token, edm::second(), findQualityFiles::size, RecoTauPiZeroBuilderPlugins_cfi::strips, tr1_token, tr2_token, testEve_cfg::tracks, ts1_token, ts2_token, w(), wd1_token, wd2_token, zeroSupCompDigis(), zeroSupStripDigis(), and SurfaceOrientation::zplus.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

 {
   bool _err = false;
   using namespace edm;
 
   typedef std::map<CSCDetId,
     std::pair<std::vector<CSCWireDigi>,std::vector<CSCWireDigi> > > 
     matchingDetWireCollection;
   typedef std::map<CSCDetId,
     std::pair<std::vector<CSCStripDigi>,std::vector<CSCStripDigi> > > 
     matchingDetStripCollection;
   typedef std::map<CSCDetId,
     std::pair<std::vector<CSCComparatorDigi>,std::vector<CSCComparatorDigi> > > 
     matchingDetComparatorCollection;
   typedef std::map<CSCDetId,
     std::pair<std::vector<CSCCLCTDigi>,std::vector<CSCCLCTDigi> > > 
     matchingDetCLCTCollection;
   typedef std::map<CSCDetId,
     std::pair<std::vector<CSCALCTDigi>,std::vector<CSCALCTDigi> > > 
     matchingDetALCTCollection;
   typedef std::map<CSCDetId,
     std::pair<std::vector<CSCCorrelatedLCTDigi>,std::vector<CSCCorrelatedLCTDigi> > > 
     matchingDetLCTCollection;
 
   edm::ESHandle<CSCChamberMap> hcham;
   iSetup.get<CSCChamberMapRcd>().get(hcham); 
   const CSCChamberMap* theMapping = hcham.product();
     
   Handle<CSCWireDigiCollection> _wi,_swi;
   Handle<CSCStripDigiCollection> _st,_sst;
   Handle<CSCComparatorDigiCollection> _cmp,_scmp;
   Handle<CSCCLCTDigiCollection> _clct, _sclct;
   Handle<CSCALCTDigiCollection> _alct, _salct;
   Handle<CSCCorrelatedLCTDigiCollection> _lct, _slct;
   Handle<L1CSCTrackCollection> _trk, _strk;
   Handle<CSCTriggerContainer<csctf::TrackStub> > _dt, _sdt;
 
    // get wire digis before and after unpacking
    iEvent.getByToken( wd1_token, _swi );
    iEvent.getByToken( wd2_token, _wi );
 
    //get strip digis before and after unpacking
    iEvent.getByToken( sd1_token, _sst );
    iEvent.getByToken( sd2_token, _st );
 
    //get comparator digis before and after unpacking
    iEvent.getByToken( co1_token, _scmp );
    iEvent.getByToken( co2_token, _cmp );
 
    //get clcts
    iEvent.getByToken( cl1_token, _sclct );
    iEvent.getByToken( cl2_token, _clct );
 
    //get alcts
    iEvent.getByToken( al1_token, _salct );
    iEvent.getByToken( al2_token, _alct );
 
    //get corr lcts
    iEvent.getByToken( co1_token, _slct );
    iEvent.getByToken( co2_token, _lct );
 
    //get l1 tracks
    iEvent.getByToken( tr1_token, _strk );
    iEvent.getByToken( tr2_token, _trk );
 
    //get DT stubs for L1 Tracks
    iEvent.getByToken( ts1_token, _sdt );
    iEvent.getByToken( ts2_token, _dt );
 
 
    CSCWireDigiCollection::DigiRangeIterator 
      wi = _wi->begin(), swi= _swi->begin();
    CSCStripDigiCollection::DigiRangeIterator 
      st = _st->begin(), sst = _sst->begin();
    CSCComparatorDigiCollection::DigiRangeIterator 
      cmp = _cmp->begin(), scmp = _scmp->begin();
    CSCCLCTDigiCollection::DigiRangeIterator
      clct = _clct->begin(), sclct = _sclct->begin();
    CSCALCTDigiCollection::DigiRangeIterator
      alct = _alct->begin(), salct = _salct->begin();
    CSCCorrelatedLCTDigiCollection::DigiRangeIterator
      lct = _lct->begin(), slct = _slct->begin();
    L1CSCTrackCollection::const_iterator 
      trk = _trk->begin(), strk = _strk->begin();
    std::vector<csctf::TrackStub>::const_iterator
      dt = _dt->get().begin(), sdt = _sdt->get().begin();
      // WARNING 5_0_X
      dt++; dt--; sdt++; sdt--;
 
    //per detID, create lists of various digi types
    matchingDetWireCollection wires;
    matchingDetStripCollection strips;
    matchingDetComparatorCollection comps;
    matchingDetCLCTCollection clcts;
    matchingDetALCTCollection alcts;
    matchingDetLCTCollection lcts,trackstubs;
    CSCTriggerContainer<csc::L1Track> tracks, simtracks;
    
    //wires
    for(;wi != _wi->end();++wi)
      {
        CSCWireDigiCollection::const_iterator 
      b=(*wi).second.first,e=(*wi).second.second;
        std::vector<CSCWireDigi>::iterator 
      beg=wires[(*wi).first].first.end();
        wires[(*wi).first].first.insert(beg,b,e);
      }
    for(;swi != _swi->end();++swi)
      {
        CSCWireDigiCollection::const_iterator 
      b=(*swi).second.first,e=(*swi).second.second;
        //convert sim ring 4(ME1/a) to ring 1
        CSCDetId _id = (*swi).first;
        if((*swi).first.ring() == 4)
      _id = CSCDetId((*swi).first.endcap(),(*swi).first.station(),
             1, (*swi).first.chamber(),(*swi).first.layer());
 
        std::vector<CSCWireDigi>::iterator 
      beg=wires[_id].second.end();
        
        wires[_id].second.insert(beg,b,e);
        // automatically combine wire digis after each insertion
        wires[_id].second = sanitizeWireDigis(wires[_id].second.begin(),
                          wires[_id].second.end());
      }
    
    //strips
    for(;st != _st->end();++st)
      {
        CSCStripDigiCollection::const_iterator 
      b=(*st).second.first,e=(*st).second.second;
        std::vector<CSCStripDigi>::iterator 
      beg=strips[(*st).first].first.end();
               
        //need to remove strips with no active ADCs
        std::vector<CSCStripDigi> zs = zeroSupStripDigis(b,e);
 
        strips[(*st).first].first.insert(beg,zs.begin(),zs.end());
      }
    for(;sst != _sst->end();++sst)
      {
        CSCStripDigiCollection::const_iterator 
      b=(*sst).second.first,e=(*sst).second.second;
        // conversion of ring 4->1 not necessary here
        CSCDetId _id = (*sst).first;
        //if((*sst).first.ring() == 4)    
        //    _id = CSCDetId((*sst).first.endcap(),(*sst).first.station(),
        //       1, (*sst).first.chamber(),(*sst).first.layer());
 
        std::vector<CSCStripDigi>::iterator 
      beg=strips[_id].second.end();       
        
        std::vector<CSCStripDigi> relab = relabelStripDigis(theMapping,(*sst).first,b,e);
         
        strips[_id].second.insert(beg,relab.begin(),relab.end());              
        //strips[_id].second.insert(beg,b,e);  
      }
    
    //comparators
    for(;cmp != _cmp->end();++cmp)
      {
        CSCComparatorDigiCollection::const_iterator 
      b=(*cmp).second.first,e=(*cmp).second.second;
        std::vector<CSCComparatorDigi>::iterator 
      beg=comps[(*cmp).first].first.end();
 
        comps[(*cmp).first].first.insert(beg,b,e);
      }
    for(;scmp != _scmp->end();++scmp)
      {
        CSCComparatorDigiCollection::const_iterator 
      b=(*scmp).second.first,e=(*scmp).second.second;
        // convert sim ring 4 (ME1/a) to ring 1
        CSCDetId _id = (*scmp).first;
        if((*scmp).first.ring() == 4)
          _id = CSCDetId((*scmp).first.endcap(),(*scmp).first.station(),
             1, (*scmp).first.chamber(),(*scmp).first.layer());
        
        std::vector<CSCComparatorDigi>::iterator 
      beg=comps[_id].second.begin();
 
        if((*scmp).first.ring()==4)
      beg=comps[_id].second.end();
 
        std::vector<CSCComparatorDigi> zs = 
      zeroSupCompDigis(b,e);           
 
        std::vector<CSCComparatorDigi> relab = 
      relabelCompDigis(theMapping,(*scmp).first,
               zs.begin(),
               zs.end());
 
        comps[_id].second.insert(beg,relab.begin(),relab.end());
      }
 
    //CLCTs
    for(;clct != _clct->end();++clct)
      {
        CSCCLCTDigiCollection::const_iterator 
      b=(*clct).second.first,e=(*clct).second.second;
        std::vector<CSCCLCTDigi>::iterator 
      beg=clcts[(*clct).first].first.end();
 
        clcts[(*clct).first].first.insert(beg,b,e);
      }
    for(;sclct != _sclct->end();++sclct)
      {
        CSCCLCTDigiCollection::const_iterator 
      b=(*sclct).second.first,e=(*sclct).second.second;
        // convert sim ring 4 (ME1/a) to ring 1
        CSCDetId _id = (*sclct).first;
        if((*sclct).first.ring() == 4)
          _id = CSCDetId((*sclct).first.endcap(),(*sclct).first.station(),
             1, (*sclct).first.chamber(),(*sclct).first.layer());
        
        std::vector<CSCCLCTDigi>::iterator 
      beg=clcts[_id].second.begin();
 
        if((*sclct).first.ring()==4)
      beg=clcts[_id].second.end();
        
        clcts[_id].second.insert(beg,b,e);
      }
 
    //ALCTs
    for(;alct != _alct->end();++alct)
      {
        CSCALCTDigiCollection::const_iterator 
      b=(*alct).second.first,e=(*alct).second.second;
        std::vector<CSCALCTDigi>::iterator 
      beg=alcts[(*alct).first].first.end();
 
        alcts[(*alct).first].first.insert(beg,b,e);
      }
    for(;salct != _salct->end();++salct)
      {
        CSCALCTDigiCollection::const_iterator 
      b=(*salct).second.first,e=(*salct).second.second;
        // convert sim ring 4 (ME1/a) to ring 1
        CSCDetId _id = (*salct).first;
        if((*salct).first.ring() == 4)
          _id = CSCDetId((*salct).first.endcap(),(*salct).first.station(),
             1, (*salct).first.chamber(),(*salct).first.layer());
        
        std::vector<CSCALCTDigi>::iterator 
      beg=alcts[_id].second.begin();
 
        if((*salct).first.ring()==4)
      beg=alcts[_id].second.end();
        
        alcts[_id].second.insert(beg,b,e);
      }
 
    // Correlated LCTs
    for(;lct != _lct->end();++lct)
      {
        CSCCorrelatedLCTDigiCollection::const_iterator 
      b=(*lct).second.first,e=(*lct).second.second;
        std::vector<CSCCorrelatedLCTDigi>::iterator 
      beg=lcts[(*lct).first].first.end();
 
        lcts[(*lct).first].first.insert(beg,b,e);
      }
    for(;slct != _slct->end();++slct)
      {
        CSCCorrelatedLCTDigiCollection::const_iterator 
      b=(*slct).second.first,e=(*slct).second.second;
        // convert sim ring 4 (ME1/a) to ring 1
        CSCDetId _id = (*slct).first;
        if((*slct).first.ring() == 4)
          _id = CSCDetId((*slct).first.endcap(),(*slct).first.station(),
             1, (*slct).first.chamber(),(*slct).first.layer());
        
        std::vector<CSCCorrelatedLCTDigi>::iterator 
      beg=lcts[_id].second.begin();
 
        if((*slct).first.ring()==4)
      beg=lcts[_id].second.end();
        
        lcts[_id].second.insert(beg,b,e);
      }
    // remove attached LCT digis from tracks, should be put into their own collection and checked separately
    for(; trk != _trk->end(); ++trk)
      {
        tracks.push_back(trk->first);
        
      }
    for(;strk != _strk->end(); ++strk)
      {
        simtracks.push_back(strk->first);
      }
 
    //now loop through each set and process if there are differences!
    matchingDetWireCollection::const_iterator w;
    matchingDetStripCollection::const_iterator s;
    matchingDetComparatorCollection::const_iterator c;
    matchingDetCLCTCollection::const_iterator cl;
    matchingDetALCTCollection::const_iterator al;
    matchingDetLCTCollection::const_iterator lc;
    
    for(w = wires.begin(); w != wires.end(); ++w)
      {
        if(w->second.first.size() != w->second.second.size())
      {
        std::cout << "Major error! # of wire digis in detID: " << w->first 
              << " is not equal between sim and unpacked!" << std::endl;
        //eventually do more in this case!
 
        std::vector<CSCWireDigi> a = w->second.second;
        std::vector<CSCWireDigi> b = w->second.first;
        std::cout << "SIM OUTPUT:" << std::endl;
        for(std::vector<CSCWireDigi>::const_iterator i = a.begin(); i != a.end(); ++i)
          i->print();
        std::cout << "UNPACKER OUTPUT:" << std::endl;
        for(std::vector<CSCWireDigi>::const_iterator i = b.begin(); i != b.end(); ++i)
          i->print();
          
      }
        int max = std::min(w->second.first.size(),w->second.second.size());
        std::vector<CSCWireDigi> cv = w->second.first;
        std::vector<CSCWireDigi> sv = w->second.second;
        for(int i = 0; i < max; ++i)
      {
        if(sv[i].getWireGroup() != cv[i].getWireGroup())
          {
            std::cout << "In detId: " << w->first << std::endl;
            std::cout << "Wire Groups do not match: " << sv[i].getWireGroup() 
              << " != " << cv[i].getWireGroup() << std::endl;           
          }
        if(sv[i].getTimeBin() != cv[i].getTimeBin())
          {
            std::cout << "In detId: " << w->first << std::endl;
            std::cout << "First Time Bins do not match: " << sv[i].getTimeBin()
              << " != " << cv[i].getTimeBin() << std::endl;         
          }     
        if(sv[i].getTimeBinWord() != cv[i].getTimeBinWord())
          {
            std::cout << "In detId: " << w->first << std::endl;
            std::cout << "Time Bin Words do not match: " << sv[i].getTimeBinWord()
              << " != " << cv[i].getTimeBinWord() << std::endl;
          }     
      }
      }
    for(s = strips.begin(); s != strips.end(); ++s)
      {
        if(s->second.first.size() != s->second.second.size())
      {
        std::cout << "Major error! # of strip digis in detID: " << s->first 
              << " is not equal between sim and unpacked!" << std::endl;
        //eventually do more in this case!
 
        std::vector<CSCStripDigi> a = s->second.second;
        std::vector<CSCStripDigi> b = s->second.first;
        std::cout << "SIM OUTPUT:" << std::endl;
        for(std::vector<CSCStripDigi>::const_iterator i = a.begin(); i != a.end(); ++i)
          i->print();
        std::cout << "UNPACKER OUTPUT:" << std::endl;
        for(std::vector<CSCStripDigi>::const_iterator i = b.begin(); i != b.end(); ++i)
          i->print();         
      }
        int max = std::min(s->second.first.size(),s->second.second.size());
        std::vector<CSCStripDigi> cv = s->second.first;
        std::vector<CSCStripDigi> sv = s->second.second;
        for(int i = 0; i < max; ++i)
      {
        bool me1a = s->first.station()==1 && s->first.ring()==4;
        bool me1b = s->first.station()==1 && s->first.ring()==1;
        bool zplus = s->first.endcap()==1;
        int k=i;
        
        if(me1a && zplus) k=max-i-1;
        if(me1b && !zplus) k=max-i-1;
 
        if(sv[k].getStrip() != cv[i].getStrip())
          {
            std::cout << "In detId: " << s->first << std::endl;
            std::cout << "Strips do not match: " << sv[k].getStrip() 
              << " != " << cv[i].getStrip() << std::endl;           
          }
        if(sv[k].getADCCounts().size() != cv[i].getADCCounts().size())
          {
            std::cout << "In detId: " << s->first << std::endl;
            std::cout << "ADC Readouts not of equal size!" << std::endl;
            std::cout << sv[k].getADCCounts().size() << ' ' 
              << cv[i].getADCCounts().size() << std::endl;
          }
        else
          {
            std::vector<int> sADC = sv[k].getADCCounts();
            std::vector<int> uADC = cv[i].getADCCounts();
 
            for(unsigned iadc = 0; iadc < sADC.size(); ++iadc)
          if(sADC[iadc] != uADC[iadc])
            {
              std::cout << "In detId: " << s->first << std::endl;
              std::cout << "ADC counts not equal at index: " << iadc << std::endl
                    << std::hex <<sADC[iadc] << " != " << uADC[iadc] << std::dec
                    << std::endl;    
            }     
          }
        if(sv[k].getADCOverflow().size() != cv[i].getADCOverflow().size())
          {
            std::cout << "In detId: " << s->first << std::endl;
            std::cout << "ADC Overflows not of equal size!" << std::endl;
            std::cout << sv[k].getADCOverflow().size() << ' ' 
              << cv[i].getADCOverflow().size() << std::endl;
          }
        else
          {
            std::vector<uint16_t> sADC = sv[k].getADCOverflow();
            std::vector<uint16_t> uADC = cv[i].getADCOverflow();
 
            for(unsigned iadc = 0; iadc < sADC.size(); ++iadc)
          if(sADC[iadc] != uADC[iadc])
            {
              std::cout << "In detId: " << s->first << std::endl;
              std::cout << "ADC overflows not equal at index: " << iadc << std::endl
                    << std::hex <<sADC[iadc] << " != " << uADC[iadc] << std::dec
                    << std::endl;         
            }
          }
        if(sv[k].getOverlappedSample().size() != cv[i].getOverlappedSample().size())
          {
            std::cout << "In detId: " << s->first << std::endl;
            std::cout << "Overlapped Samples not of equal size!" << std::endl;
            std::cout << sv[k].getOverlappedSample().size() << ' ' 
              << cv[i].getOverlappedSample().size() << std::endl;
          }
        else
          {
            std::vector<uint16_t> sADC = sv[k].getOverlappedSample();
            std::vector<uint16_t> uADC = cv[i].getOverlappedSample();
 
            for(unsigned iadc = 0; iadc < sADC.size(); ++iadc)
          if(sADC[iadc] != uADC[iadc])
            {
              std::cout << "In detId: " << s->first << std::endl;
              std::cout << "Overlapped Samples not equal at index: " << iadc << std::endl
                    << std::hex <<sADC[iadc] << " != " << uADC[iadc] << std::dec
                    << std::endl;         
            }
          }
        if(sv[k].getErrorstat().size() != cv[i].getErrorstat().size())
          {
            std::cout << "In detId: " << s->first << std::endl;
            std::cout << "Errorstat not of equal size!" << std::endl;
            std::cout << sv[k].getErrorstat().size() << ' ' 
              << cv[i].getErrorstat().size() << std::endl;
          }
        else
          {
            std::vector<uint16_t> sADC = sv[k].getErrorstat();
            std::vector<uint16_t> uADC = cv[i].getErrorstat();
 
            for(unsigned iadc = 0; iadc < sADC.size(); ++iadc)
          if(sADC[iadc] != uADC[iadc])
            {
              std::cout << "In detId: " << s->first << std::endl;
              std::cout << "Errorstat not equal at index: " << iadc << std::endl
                    << std::hex <<sADC[iadc] << " != " << uADC[iadc] << std::dec
                    << std::endl;         
            }
          }
        if(sv[k].pedestal() != cv[i].pedestal())
          {
            std::cout << "In detId: " << s->first << std::endl;
            std::cout << "Pedestals not equal: " << sv[k].pedestal() << " != " 
              << cv[i].pedestal() << std::endl;
          }
        if(sv[k].amplitude() != cv[i].amplitude())
          {
            std::cout << "In detId: " << s->first << std::endl;
            std::cout << "Amplitudes not equal: " << sv[k].amplitude() << " != " 
              << cv[i].amplitude() << std::endl;
          }
      }
      }
    for(c = comps.begin(); c != comps.end(); ++c)
      {
        if(c->second.first.size() != c->second.second.size())
      {
        std::cout << "Major error! # of comparator digis in detID: " << c->first
              << " is not equal between sim and unpacked!" << std::endl;
        //eventually do more in this case!
 
        std::vector<CSCComparatorDigi> a = c->second.second;
        std::vector<CSCComparatorDigi> b = c->second.first;
        std::cout << "SIM OUTPUT:" << std::endl;
        for(std::vector<CSCComparatorDigi>::const_iterator i = a.begin(); i != a.end(); ++i)
          i->print();
        std::cout << "UNPACKER OUTPUT:" << std::endl;
        for(std::vector<CSCComparatorDigi>::const_iterator i = b.begin(); i != b.end(); ++i)
          i->print();         
      }
        int max = std::min(c->second.first.size(),c->second.second.size());
        std::vector<CSCComparatorDigi> cv = c->second.first;
        std::vector<CSCComparatorDigi> sv = c->second.second;
        for(int i = 0; i < max; ++i)
      {     
        if(sv[i].getStrip() != cv[i].getStrip())
          {
            std::cout << "In detId: " << s->first << std::endl;
            std::cout << "Comparator strips do not match: " << sv[i].getStrip() 
              << " != " << cv[i].getStrip() << std::endl;
          }
        if(sv[i].getComparator() != cv[i].getComparator())
          {         
            std::cout << "In detId: " << c->first << std::endl;
            std::cout << "Comparators do not match: " << sv[i].getComparator()
              << " != " << cv[i].getComparator() << std::endl;         
          }
        if(sv[i].getTimeBinWord() != cv[i].getTimeBinWord())
          {
            std::cout << "In detId: " << c->first << std::endl;
            std::cout << "Comparator time bins words do not match: " << sv[i].getTimeBinWord()
              << " != " << cv[i].getTimeBinWord() << std::endl;
          }
      }
      }
    for(cl = clcts.begin(); cl != clcts.end(); ++cl)
      {
        if(cl->second.first.size() != cl->second.second.size())
      {
        std::cout << "Major error! # of CLCT digis in detID: " << cl->first
              << " is not equal between sim and unpacked!" << std::endl;
        //eventually do more in this case!
 
        std::vector<CSCCLCTDigi> a = cl->second.second;
        std::vector<CSCCLCTDigi> b = cl->second.first;
        std::cout << "SIM OUTPUT:" << std::endl;
        for(std::vector<CSCCLCTDigi>::const_iterator i = a.begin(); i != a.end(); ++i)
          i->print();
        std::cout << "UNPACKER OUTPUT:" << std::endl;
        for(std::vector<CSCCLCTDigi>::const_iterator i = b.begin(); i != b.end(); ++i)
          i->print();         
      }
        int max = std::min(cl->second.first.size(),cl->second.second.size());
        std::vector<CSCCLCTDigi> cv = cl->second.first;
        std::vector<CSCCLCTDigi> sv = cl->second.second;
        for(int i = 0; i < max; ++i)
      {
        if(cv[i].getKeyStrip() != sv[i].getKeyStrip())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT key strips do not match: " << sv[i].getKeyStrip()
              << " != " << cv[i].getKeyStrip() << std::endl;
          }
        if(cv[i].getStrip() != sv[i].getStrip())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT strips do not match: " << sv[i].getStrip()
              << " != " << cv[i].getStrip() << std::endl;
          }     
        if(cv[i].isValid() != sv[i].isValid())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT Valid bits do not match: " << sv[i].isValid() 
              << " != " << cv[i].isValid() << std::endl;
          }
        if(cv[i].getQuality() != sv[i].getQuality())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT qualities do not match: " << sv[i].getQuality()
              << " != " << cv[i].getQuality() << std::endl;
          }
        if(cv[i].getPattern() != sv[i].getPattern())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT patterns do not match: " << sv[i].getPattern()
              << " != " << cv[i].getPattern() << std::endl;
          }
        if(cv[i].getStripType() != sv[i].getStripType())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT strip types do not match: " << sv[i].getStripType()
              << " != " << cv[i].getStripType() << std::endl;
          }
        if(cv[i].getBend() != sv[i].getBend())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT bends do not match: " << sv[i].getBend()
              << " != " << cv[i].getBend() << std::endl;
          }
        if(cv[i].getCFEB() != sv[i].getCFEB())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT CFEBs do not match: " << sv[i].getCFEB()
              << " != " << cv[i].getCFEB() << std::endl;
          }
        if(((short)cv[i].getBX()) != ((short)sv[i].getBX()) - 4)
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT BX do not match: " << sv[i].getBX() - 4
              << " != " << cv[i].getBX() << std::endl;
          }
        if(cv[i].getFullBX() != sv[i].getFullBX())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT Full BX do not match: " << sv[i].getFullBX()
              << " != " << cv[i].getFullBX() << std::endl;
          }
        if(cv[i].getTrknmb() != sv[i].getTrknmb())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "CLCT Track numbers do not match: " << sv[i].getTrknmb()
              << " != " << cv[i].getTrknmb() << std::endl;
          }     
      }
      }
    for(al = alcts.begin(); al != alcts.end(); ++al)
      {
        if(al->second.first.size() != al->second.second.size())
      {
        std::cout << "Major error! # of ALCT digis in detID: " << al->first
              << " is not equal between sim and unpacked!" << std::endl;
        //eventually do more in this case!
 
        std::vector<CSCALCTDigi> a = al->second.second;
        std::vector<CSCALCTDigi> b = al->second.first;
        std::cout << "SIM OUTPUT:" << std::endl;
        for(std::vector<CSCALCTDigi>::const_iterator i = a.begin(); i != a.end(); ++i)
          i->print();
        std::cout << "UNPACKER OUTPUT:" << std::endl;
        for(std::vector<CSCALCTDigi>::const_iterator i = b.begin(); i != b.end(); ++i)
          i->print();         
      }
        int max = std::min(al->second.first.size(),al->second.second.size());
        std::vector<CSCALCTDigi> cv = al->second.first;
        std::vector<CSCALCTDigi> sv = al->second.second;
        for(int i = 0; i < max; ++i)
      {
        if(cv[i].getKeyWG() != sv[i].getKeyWG())
          {
            std::cout << "In detId: " << al->first << std::endl;
            std::cout << "ALCT key wire groups do not match: " << sv[i].getKeyWG()
              << " != " << cv[i].getKeyWG() << std::endl;
          }
        if(cv[i].isValid() != sv[i].isValid())
          {
            std::cout << "In detId: " << al->first << std::endl;
            std::cout << "ALCT Valid bits do not match: " << sv[i].isValid() 
              << " != " << cv[i].isValid() << std::endl;
          }
        if(cv[i].getQuality() != sv[i].getQuality())
          {
            std::cout << "In detId: " << al->first << std::endl;
            std::cout << "ALCT qualities do not match: " << sv[i].getQuality()
              << " != " << cv[i].getQuality() << std::endl;
          }
        if(cv[i].getAccelerator() != sv[i].getAccelerator())
          {
            std::cout << "In detId: " << al->first << std::endl;
            std::cout << "ALCT accelerator bits do not match: " << sv[i].getAccelerator()
              << " != " << cv[i].getAccelerator() << std::endl;
          }
        if(cv[i].getCollisionB() != sv[i].getCollisionB())
          {
            std::cout << "In detId: " << al->first << std::endl;
            std::cout << "ALCT CollisionB flags do not match: " << sv[i].getCollisionB()
              << " != " << cv[i].getCollisionB() << std::endl;
          }
        if((cv[i].getBX()) != (sv[i].getBX()))
          {
            std::cout << "In detId: " << al->first << std::endl;
            std::cout << "ALCT BX do not match: " << sv[i].getBX()
              << " != " << cv[i].getBX() << std::endl;
          }
        if(cv[i].getFullBX() != sv[i].getFullBX())
          {
            std::cout << "In detId: " << cl->first << std::endl;
            std::cout << "ALCT Full BX do not match: " << sv[i].getFullBX()
              << " != " << cv[i].getFullBX() << std::endl;
          }
      }
      }
    for(lc = lcts.begin(); lc != lcts.end(); ++lc)
      {
        if(lc->second.first.size() != lc->second.second.size())
      {
        std::cout << "Major error! # of Correlated LCT digis in detID: " << lc->first
              << " is not equal between sim and unpacked!" << std::endl;
        //eventually do more in this case!
 
        std::vector<CSCCorrelatedLCTDigi> a = lc->second.second;
        std::vector<CSCCorrelatedLCTDigi> b = lc->second.first;
        std::cout << "SIM OUTPUT:" << std::endl;
        for(std::vector<CSCCorrelatedLCTDigi>::const_iterator i = a.begin(); i != a.end(); ++i)
          i->print();
        std::cout << "UNPACKER OUTPUT:" << std::endl;
        for(std::vector<CSCCorrelatedLCTDigi>::const_iterator i = b.begin(); i != b.end(); ++i)
          i->print();         
      }
        int max = std::min(lc->second.first.size(),lc->second.second.size());
        std::vector<CSCCorrelatedLCTDigi> cv = lc->second.first;
        std::vector<CSCCorrelatedLCTDigi> sv = lc->second.second;
        for(int i = 0; i < max; ++i)
      {
        if(cv[i].getStrip() != sv[i].getStrip())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT strips do not match: " << sv[i].getStrip()
              << " != " << cv[i].getStrip() << std::endl;
          }   
        if(cv[i].getKeyWG() != sv[i].getKeyWG())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT key wire groups do not match: " << sv[i].getKeyWG()
              << " != " << cv[i].getKeyWG() << std::endl;
          }
        if(cv[i].isValid() != sv[i].isValid())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT Valid bits do not match: " << sv[i].isValid() 
              << " != " << cv[i].isValid() << std::endl;
          }
        if(cv[i].getQuality() != sv[i].getQuality())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT qualities do not match: " << sv[i].getQuality()
              << " != " << cv[i].getQuality() << std::endl;
          }
        if(cv[i].getPattern() != sv[i].getPattern())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT ALCT patterns do not match: " << sv[i].getPattern()
              << " != " << cv[i].getPattern() << std::endl;
          }
        if(cv[i].getCLCTPattern() != sv[i].getCLCTPattern())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT CLCT patterns do not match: " << sv[i].getCLCTPattern()
              << " != " << cv[i].getCLCTPattern() << std::endl;
          }
        if(cv[i].getStripType() != sv[i].getStripType())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT strip types do not match: " << sv[i].getStripType()
              << " != " << cv[i].getStripType() << std::endl;
          }
        if(cv[i].getBend() != sv[i].getBend())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT bends do not match: " << sv[i].getBend()
              << " != " << cv[i].getBend() << std::endl;
          }
        if(cv[i].getMPCLink() != sv[i].getMPCLink())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT MPC Links do not match: " << sv[i].getMPCLink()
              << " != " << cv[i].getMPCLink() << std::endl;
          }
        if((cv[i].getBX()) != (sv[i].getBX()-6))
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT BX do not match: " << sv[i].getBX()-6
              << " != " << cv[i].getBX() << std::endl;
          }
        if(cv[i].getCSCID() != sv[i].getCSCID())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT CSCIDs do not match: " << sv[i].getCSCID()
              << " != " << cv[i].getCSCID() << std::endl;
          }
        if(cv[i].getBX0() != sv[i].getBX0())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT BX0s do not match: " << sv[i].getBX0()
              << " != " << cv[i].getBX0() << std::endl;
          }
        if(cv[i].getSyncErr() != sv[i].getSyncErr())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT SyncErrs do not match: " << sv[i].getSyncErr()
              << " != " << cv[i].getSyncErr() << std::endl;
          }
        if(cv[i].getTrknmb() != sv[i].getTrknmb())
          {
            std::cout << "In detId: " << lc->first << std::endl;
            std::cout << "Correlated LCT Track numbers do not match: " << sv[i].getTrknmb()
              << " != " << cv[i].getTrknmb() << std::endl;
          }     
      }
      }
    if(tracks.get().size() != simtracks.get().size())
      {
        std::cout << "Major error! # of L1 Tracks is not equal between sim and unpacked!" << std::endl;
        std::vector<csc::L1Track> a = simtracks.get();
        std::vector<csc::L1Track> b = tracks.get();
        std::cout << "SIM OUTPUT:" << std::endl;
        for(std::vector<csc::L1Track>::const_iterator i = a.begin(); i != a.end(); ++i)
      i->print();
        std::cout << "UNPACKER OUTPUT:" << std::endl;
        for(std::vector<csc::L1Track>::const_iterator i = b.begin(); i != b.end(); ++i)
      i->print();
      }
    
 
    return _err;
 }

std::vector< CSCComparatorDigi > CSCDigiValidator::relabelCompDigis	(	const CSCChamberMap *	m,
		CSCDetId	_id,
		std::vector< CSCComparatorDigi >::const_iterator	b,
		std::vector< CSCComparatorDigi >::const_iterator	e
	)

private

Definition at line 915 of file CSCDigiValidator.cc.

References alignCSCRings::e, i, CSCDetId::ring(), and CSCDetId::station().

Referenced by filter().

 {
   std::vector<CSCComparatorDigi> _r; // the vector of comp digis with appropriate strip #'s
 
   bool me1a = _id.station()==1 && _id.ring()==4;
   //bool zplus = _id.endcap()==1;
   //bool me1b = _id.station()==1 && _id.ring()==1;
 
   for(std::vector<CSCComparatorDigi>::const_iterator i = b; i != e; ++i)
   {
     int strip=i->getStrip();
     
     //    if(me1a&&zplus) strip=17-strip;
     //    if(me1b&&!zplus) strip=65-strip;
     if(me1a) strip+=64;
     
     _r.push_back(CSCComparatorDigi(strip,i->getComparator(),
                    i->getTimeBinWord()));
   }
   return _r;
 }

std::vector< CSCStripDigi > CSCDigiValidator::relabelStripDigis	(	const CSCChamberMap *	m,
		CSCDetId	_id,
		std::vector< CSCStripDigi >::const_iterator	b,
		std::vector< CSCStripDigi >::const_iterator	e
	)

private

Definition at line 890 of file CSCDigiValidator.cc.

References alignCSCRings::e, and i.

Referenced by filter().

 {
   std::vector<CSCStripDigi> _r; // the vector of strip digis with appropriate strip #'s
 
   //bool me1a = _id.station()==1 && _id.ring()==4;
   //bool zplus = _id.endcap()==1;
   //bool me1b = _id.station()==1 && _id.ring()==1;
 
   for(std::vector<CSCStripDigi>::const_iterator i = b; i != e; ++i)
   {
     int strip=i->getStrip();
     
     //if(me1a&&zplus) strip=17-strip;
     //if(me1b&&!zplus) strip=(65-strip-1)%(m->dmb(_id)*16) + 1;
     //if(me1a) strip+=64;
     
     _r.push_back(CSCStripDigi(strip,i->getADCCounts(),i->getADCOverflow(),
                  i->getOverlappedSample(),i->getErrorstat()));
   }
   return _r;
 }

std::vector< CSCStripDigi > CSCDigiValidator::sanitizeStripDigis	(	std::vector< CSCStripDigi >::const_iterator	b,
		std::vector< CSCStripDigi >::const_iterator	e
	)

private

Definition at line 940 of file CSCDigiValidator.cc.

 {
   std::vector<CSCStripDigi> _r; // vector of digis in proper order
 
   return _r;
 }

std::vector< CSCWireDigi > CSCDigiValidator::sanitizeWireDigis	(	std::vector< CSCWireDigi >::const_iterator	b,
		std::vector< CSCWireDigi >::const_iterator	e
	)

private

Definition at line 858 of file CSCDigiValidator.cc.

References alignCSCRings::e, i, and edmStreamStallGrapher::t.

Referenced by filter().

 {
   typedef std::map<int,std::vector<CSCWireDigi> > wire2digi;
 
   std::vector<CSCWireDigi> _r; // the resulting vector of wire digis
   wire2digi _wr2digis; // map of wires to a set of digis
 
   for(std::vector<CSCWireDigi>::const_iterator i = b; i != e; ++i)
     _wr2digis[i->getWireGroup()].push_back(*i);
 
   for(wire2digi::const_iterator i = _wr2digis.begin(); i != _wr2digis.end(); ++i)
     {
       int wire = i->first;
       unsigned tbin = 0x0;
 
       for(std::vector<CSCWireDigi>::const_iterator d = i->second.begin();
       d != i->second.end(); ++d)
     {
       std::vector<int> binson = d->getTimeBinsOn();
       for(std::vector<int>::const_iterator t = binson.begin(); 
           t != binson.end(); ++t)
         tbin |= 1<<(*t);
     }
       
       _r.push_back(CSCWireDigi(wire,tbin));
     }
 
   return _r;
 }

std::vector< CSCComparatorDigi > CSCDigiValidator::zeroSupCompDigis	(	std::vector< CSCComparatorDigi >::const_iterator	b,
		std::vector< CSCComparatorDigi >::const_iterator	e
	)

private

Definition at line 970 of file CSCDigiValidator.cc.

References alignCSCRings::e, and i.

Referenced by filter().

 {
   std::vector<CSCComparatorDigi> _r;
 
   for(std::vector<CSCComparatorDigi>::const_iterator i = b; i != e; ++i)
     {
       bool present = false;
 
       if(i->getTimeBin() < 10) present=true;
 
       if(present) _r.push_back(*i);
     }
 
   return _r;
 }

std::vector< CSCStripDigi > CSCDigiValidator::zeroSupStripDigis	(	std::vector< CSCStripDigi >::const_iterator	b,
		std::vector< CSCStripDigi >::const_iterator	e
	)

private

Definition at line 949 of file CSCDigiValidator.cc.

References a, alignCSCRings::e, i, and beamvalidation::nonzero().

Referenced by filter().

 {
   std::vector<CSCStripDigi> _r; // zero-suppressed strip digis
   std::vector<int> counts;
 
   for(std::vector<CSCStripDigi>::const_iterator i = b; i != e; ++i)
     {
       bool nonzero=false;
       counts = i->getADCCounts();
       for(std::vector<int>::const_iterator a = counts.begin(); a != counts.end(); ++a)
     if((*a) != 0) nonzero = true;
 
       if(nonzero) _r.push_back(*i);
     }
 
   return _r;
 }