CSCDigiValidator Class Reference

Inheritance diagram for CSCDigiValidator:

Public Member Functions
	CSCDigiValidator (const edm::ParameterSet &)
	~CSCDigiValidator () override
Public Member Functions inherited from edm::one::EDFilter<>
	EDFilter ()=default
	EDFilter (const EDFilter &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDFilter &	operator= (const EDFilter &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept final
Public Member Functions inherited from edm::one::EDFilterBase
	EDFilterBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const noexcept
bool	wantsStreamRuns () const noexcept
	~EDFilterBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
std::vector< bool > const &	recordProvenanceList () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob () override
void	endJob () override
bool	filter (edm::Event &, const edm::EventSetup &) override
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::ESGetToken< CSCChamberMap, CSCChamberMapRcd >	cham_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
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::one::EDFilterBase
typedef EDFilterBase	ModuleType
Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDFilterBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
template<Transition Tr = Transition::Event>
auto	produces () noexcept
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 48 of file CSCDigiValidator.cc.

Constructor & Destructor Documentation

CSCDigiValidator()

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

explicit

Definition at line 99 of file CSCDigiValidator.cc.

References al1_token, al2_token, cd1_token, cd2_token, cham_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"));
  cham_token = esConsumes<CSCChamberMap, CSCChamberMapRcd>();
}

~CSCDigiValidator()

CSCDigiValidator::~CSCDigiValidator ( )

override

Definition at line 119 of file CSCDigiValidator.cc.

{}

Member Function Documentation

beginJob()

void CSCDigiValidator::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDFilterBase.

Definition at line 869 of file CSCDigiValidator.cc.

{}

endJob()

void CSCDigiValidator::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDFilterBase.

Definition at line 871 of file CSCDigiValidator.cc.

{}

filter()

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

overrideprivatevirtual

Implements edm::one::EDFilterBase.

Definition at line 121 of file CSCDigiValidator.cc.

References a, al1_token, al2_token, CustomPhysics_cfi::amplitude, b, DummyCfis::c, cham_token, haddnano::cl, cl1_token, cl2_token, fileCollector::cmp, co1_token, co2_token, gather_cfg::cout, cuy::cv, TauDecayModes::dec, MillePedeFileConverter_cfg::e, CSCTriggerContainer< T >::get(), edm::EventSetup::getHandle(), mps_fire::i, iEvent, sistrip::SpyUtilities::isValid(), isotrackApplyRegressor::k, WZElectronSkims53X_cff::max, SiStripPI::min, Hcal_Conditions_forGlobalTag_cff::pedestal, edm::ESHandle< T >::product(), CSCTriggerContainer< T >::push_back(), relabelCompDigis(), relabelStripDigis(), alignCSCRings::s, sanitizeWireDigis(), sd1_token, sd2_token, edm::second(), findQualityFiles::size, DigiDM_cff::strips, pfDeepBoostedJetPreprocessParams_cfi::sv, tr1_token, tr2_token, DiMuonV_cfg::tracks, ts1_token, ts2_token, w(), wd1_token, wd2_token, DigiDM_cff::wires, zeroSupCompDigis(), zeroSupStripDigis(), and SurfaceOrientation::zplus.

                                                                           {
  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.getHandle(cham_token);
  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();

  //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;
}

relabelCompDigis()

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 802 of file CSCDigiValidator.cc.

References b, MillePedeFileConverter_cfg::e, mps_fire::i, CSCDetId::ring(), CSCDetId::station(), and nano_mu_digi_cff::strip.

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;
}

relabelStripDigis()

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 779 of file CSCDigiValidator.cc.

References b, MillePedeFileConverter_cfg::e, mps_fire::i, and nano_mu_digi_cff::strip.

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;
}

sanitizeStripDigis()

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

private

Definition at line 825 of file CSCDigiValidator.cc.

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

  return _r;
}

sanitizeWireDigis()

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

private

Definition at line 753 of file CSCDigiValidator.cc.

References b, ztail::d, MillePedeFileConverter_cfg::e, mps_fire::i, submitPVValidationJobs::t, and nano_mu_digi_cff::wire.

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;
}

zeroSupCompDigis()

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

private

Definition at line 852 of file CSCDigiValidator.cc.

References b, MillePedeFileConverter_cfg::e, and mps_fire::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;
}

zeroSupStripDigis()

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

private

Definition at line 832 of file CSCDigiValidator.cc.

References a, b, dqmiodumpmetadata::counts, MillePedeFileConverter_cfg::e, mps_fire::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;
}

Member Data Documentation

al1_token

edm::EDGetTokenT<CSCALCTDigiCollection> CSCDigiValidator::al1_token

private

Definition at line 82 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

al2_token

edm::EDGetTokenT<CSCALCTDigiCollection> CSCDigiValidator::al2_token

private

Definition at line 91 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

alct1

edm::InputTag CSCDigiValidator::alct1

private

Definition at line 76 of file CSCDigiValidator.cc.

alct2

edm::InputTag CSCDigiValidator::alct2

private

Definition at line 77 of file CSCDigiValidator.cc.

cd1_token

edm::EDGetTokenT<CSCComparatorDigiCollection> CSCDigiValidator::cd1_token

private

Definition at line 81 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator().

cd2_token

edm::EDGetTokenT<CSCComparatorDigiCollection> CSCDigiValidator::cd2_token

private

Definition at line 90 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator().

cham_token

edm::ESGetToken<CSCChamberMap, CSCChamberMapRcd> CSCDigiValidator::cham_token

private

Definition at line 96 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

cl1_token

edm::EDGetTokenT<CSCCLCTDigiCollection> CSCDigiValidator::cl1_token

private

Definition at line 83 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

cl2_token

edm::EDGetTokenT<CSCCLCTDigiCollection> CSCDigiValidator::cl2_token

private

Definition at line 92 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

clct1

edm::InputTag CSCDigiValidator::clct1

private

Definition at line 76 of file CSCDigiValidator.cc.

clct2

edm::InputTag CSCDigiValidator::clct2

private

Definition at line 77 of file CSCDigiValidator.cc.

co1_token

edm::EDGetTokenT<CSCCorrelatedLCTDigiCollection> CSCDigiValidator::co1_token

private

Definition at line 84 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

co2_token

edm::EDGetTokenT<CSCCorrelatedLCTDigiCollection> CSCDigiValidator::co2_token

private

Definition at line 93 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

comp1

edm::InputTag CSCDigiValidator::comp1

private

Definition at line 76 of file CSCDigiValidator.cc.

comp2

edm::InputTag CSCDigiValidator::comp2

private

Definition at line 77 of file CSCDigiValidator.cc.

csctf1

edm::InputTag CSCDigiValidator::csctf1

private

Definition at line 76 of file CSCDigiValidator.cc.

csctf2

edm::InputTag CSCDigiValidator::csctf2

private

Definition at line 77 of file CSCDigiValidator.cc.

csctfstubs1

edm::InputTag CSCDigiValidator::csctfstubs1

private

Definition at line 76 of file CSCDigiValidator.cc.

csctfstubs2

edm::InputTag CSCDigiValidator::csctfstubs2

private

Definition at line 77 of file CSCDigiValidator.cc.

lct1

edm::InputTag CSCDigiValidator::lct1

private

Definition at line 76 of file CSCDigiValidator.cc.

lct2

edm::InputTag CSCDigiValidator::lct2

private

Definition at line 77 of file CSCDigiValidator.cc.

sd1_token

edm::EDGetTokenT<CSCStripDigiCollection> CSCDigiValidator::sd1_token

private

Definition at line 80 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

sd2_token

edm::EDGetTokenT<CSCStripDigiCollection> CSCDigiValidator::sd2_token

private

Definition at line 89 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

strip1

edm::InputTag CSCDigiValidator::strip1

private

Definition at line 76 of file CSCDigiValidator.cc.

strip2

edm::InputTag CSCDigiValidator::strip2

private

Definition at line 77 of file CSCDigiValidator.cc.

tr1_token

edm::EDGetTokenT<L1CSCTrackCollection> CSCDigiValidator::tr1_token

private

Definition at line 85 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

tr2_token

edm::EDGetTokenT<L1CSCTrackCollection> CSCDigiValidator::tr2_token

private

Definition at line 94 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

ts1_token

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

private

Definition at line 86 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

ts2_token

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

private

Definition at line 95 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

wd1_token

edm::EDGetTokenT<CSCWireDigiCollection> CSCDigiValidator::wd1_token

private

Definition at line 79 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

wd2_token

edm::EDGetTokenT<CSCWireDigiCollection> CSCDigiValidator::wd2_token

private

Definition at line 88 of file CSCDigiValidator.cc.

Referenced by CSCDigiValidator(), and filter().

wire1

edm::InputTag CSCDigiValidator::wire1

private

Definition at line 76 of file CSCDigiValidator.cc.

wire2

edm::InputTag CSCDigiValidator::wire2

private

Definition at line 77 of file CSCDigiValidator.cc.