CMS 3D CMS Logo

Public Member Functions | Private Member Functions | Private Attributes

CSCDigiValidator Class Reference

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

Inheritance diagram for CSCDigiValidator:
edm::EDFilter edm::ProducerBase edm::ProductRegistryHelper

List of all members.

Public Member Functions

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

Private Member Functions

virtual void beginJob ()
virtual void endJob ()
virtual bool filter (edm::Event &, const edm::EventSetup &)
std::vector< CSCComparatorDigirelabelCompDigis (const CSCChamberMap *m, CSCDetId _id, std::vector< CSCComparatorDigi >::const_iterator b, std::vector< CSCComparatorDigi >::const_iterator e)
std::vector< CSCStripDigirelabelStripDigis (const CSCChamberMap *, CSCDetId, std::vector< CSCStripDigi >::const_iterator, std::vector< CSCStripDigi >::const_iterator)
std::vector< CSCStripDigisanitizeStripDigis (std::vector< CSCStripDigi >::const_iterator, std::vector< CSCStripDigi >::const_iterator)
std::vector< CSCWireDigisanitizeWireDigis (std::vector< CSCWireDigi >::const_iterator, std::vector< CSCWireDigi >::const_iterator)
std::vector< CSCComparatorDigizeroSupCompDigis (std::vector< CSCComparatorDigi >::const_iterator, std::vector< CSCComparatorDigi >::const_iterator)
std::vector< CSCStripDigizeroSupStripDigis (std::vector< CSCStripDigi >::const_iterator, std::vector< CSCStripDigi >::const_iterator)

Private Attributes

edm::InputTag alct1
edm::InputTag alct2
edm::InputTag clct1
edm::InputTag clct2
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::InputTag strip1
edm::InputTag strip2
edm::InputTag wire1
edm::InputTag wire2

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 39 of file CSCDigiValidator.h.


Constructor & Destructor Documentation

CSCDigiValidator::CSCDigiValidator ( const edm::ParameterSet iConfig) [explicit]

Definition at line 55 of file CSCDigiValidator.cc.

                                                                 :
  wire1(iConfig.getUntrackedParameter<edm::InputTag>("inputWire")),
  strip1(iConfig.getUntrackedParameter<edm::InputTag>("inputStrip")),
  comp1(iConfig.getUntrackedParameter<edm::InputTag>("inputComp")),
  clct1(iConfig.getUntrackedParameter<edm::InputTag>("inputCLCT")),
  alct1(iConfig.getUntrackedParameter<edm::InputTag>("inputALCT")),
  lct1(iConfig.getUntrackedParameter<edm::InputTag>("inputCorrLCT")),
  csctf1(iConfig.getUntrackedParameter<edm::InputTag>("inputCSCTF")),
  csctfstubs1(iConfig.getUntrackedParameter<edm::InputTag>("inputCSCTFStubs")),
  wire2(iConfig.getUntrackedParameter<edm::InputTag>("repackWire")),
  strip2(iConfig.getUntrackedParameter<edm::InputTag>("repackStrip")),
  comp2(iConfig.getUntrackedParameter<edm::InputTag>("repackComp")),
  clct2(iConfig.getUntrackedParameter<edm::InputTag>("repackCLCT")),
  alct2(iConfig.getUntrackedParameter<edm::InputTag>("repackALCT")),
  lct2(iConfig.getUntrackedParameter<edm::InputTag>("repackCorrLCT")),
  csctf2(iConfig.getUntrackedParameter<edm::InputTag>("repackCSCTF")),
  csctfstubs2(iConfig.getUntrackedParameter<edm::InputTag>("repackCSCTFStubs"))
  //  reorderStrips(iConfig.getUntrackedParameter<bool>("applyStripReordering",true))
{
   //now do what ever initialization is needed
  //  produces<std::map<std::string,unsigned> >(); // # of errors enumerated by error type
}
CSCDigiValidator::~CSCDigiValidator ( )

Definition at line 79 of file CSCDigiValidator.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

void CSCDigiValidator::beginJob ( void  ) [private, virtual]

Reimplemented from edm::EDFilter.

Definition at line 1030 of file CSCDigiValidator.cc.

{
}
void CSCDigiValidator::endJob ( void  ) [private, virtual]

Reimplemented from edm::EDFilter.

Definition at line 1036 of file CSCDigiValidator.cc.

{
}
bool CSCDigiValidator::filter ( edm::Event iEvent,
const edm::EventSetup iSetup 
) [private, virtual]

Implements edm::EDFilter.

Definition at line 94 of file CSCDigiValidator.cc.

References a, alct1, alct2, b, trackerHits::c, clct1, clct2, comp1, comp2, gather_cfg::cout, CSCDetId, csctf1, csctf2, csctfstubs1, csctfstubs2, dt, edm::EventSetup::get(), edm::Event::getByLabel(), i, gen::k, lct1, lct2, python::multivaluedict::map(), max(), min, relabelCompDigis(), relabelStripDigis(), asciidump::s, sanitizeWireDigis(), edm::second(), findQualityFiles::size, strip1, strip2, RecoTauPiZeroBuilderPlugins_cfi::strips, testEve_cfg::tracks, wire1, wire2, 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;

  // std::auto_ptr<std::map<std::string,unsigned> >
  //    errors(new std::map<std::string,unsigned>);

  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.getByLabel(wire1,_swi);
   iEvent.getByLabel(wire2,_wi);
   
   //get strip digis before and after unpacking
   iEvent.getByLabel(strip1,_sst);
   iEvent.getByLabel(strip2,_st);

   //get comparator digis before and after unpacking
   iEvent.getByLabel(comp1,_scmp);
   iEvent.getByLabel(comp2,_cmp);

   //get clcts
   iEvent.getByLabel(clct1,_sclct);
   iEvent.getByLabel(clct2,_clct);

   //get alcts
   iEvent.getByLabel(alct1,_salct);
   iEvent.getByLabel(alct2,_alct);

   //get lcts
   iEvent.getByLabel(lct1,_slct);
   iEvent.getByLabel(lct2,_lct);

   //get l1 tracks
   iEvent.getByLabel(csctf1,_strk);
   iEvent.getByLabel(csctfstubs1,_sdt);
   iEvent.getByLabel(csctf2,_trk);
   iEvent.getByLabel(csctfstubs2,_dt);
   
   //get DT stubs for L1 Tracks

   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();

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

   //   iEvent.put(errors);
   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 956 of file CSCDigiValidator.cc.

References i, CSCDetId::ring(), CSCDetId::station(), and 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;
}
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 931 of file CSCDigiValidator.cc.

References i, and 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;
}
std::vector< CSCStripDigi > CSCDigiValidator::sanitizeStripDigis ( std::vector< CSCStripDigi >::const_iterator  b,
std::vector< CSCStripDigi >::const_iterator  e 
) [private]

Definition at line 981 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 899 of file CSCDigiValidator.cc.

References i, and matplotRender::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 1011 of file CSCDigiValidator.cc.

References i, and edm::productstatus::present().

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

References a, 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

Definition at line 71 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 72 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 71 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 72 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 71 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 72 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 71 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 72 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 71 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 72 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 71 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 72 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 74 of file CSCDigiValidator.h.

Definition at line 71 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 72 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 71 of file CSCDigiValidator.h.

Referenced by filter().

Definition at line 72 of file CSCDigiValidator.h.

Referenced by filter().