HcalCableMapper Class Reference

Inheritance diagram for HcalCableMapper:

Public Member Functions
virtual void	analyze (edm::Event const &e, edm::EventSetup const &c)
virtual void	endJob ()
	HcalCableMapper (edm::ParameterSet const &conf)
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Types
typedef std::map< HcalDetId, HcalElectronicsId >	IdMap
typedef std::map< HcalDetId, SampleSet >	PathSet
typedef std::vector < HcalQIESample >	SampleSet

Private Member Functions
void	process (const PathSet &ps, const IdMap &im)
template<class DigiCollection >
void	record (const DigiCollection &digis)

Private Attributes
std::map< HcalDetId, std::vector< SampleSet > >	fullHistory_
edm::InputTag	hbheLabel_
edm::InputTag	hfLabel_
edm::InputTag	hoLabel_
IdMap	IdSet

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 16 of file HcalCableMapper.cc.

Member Typedef Documentation

typedef std::map<HcalDetId,HcalElectronicsId> HcalCableMapper::IdMap

private

Definition at line 26 of file HcalCableMapper.cc.

typedef std::map<HcalDetId,SampleSet> HcalCableMapper::PathSet

private

Definition at line 25 of file HcalCableMapper.cc.

typedef std::vector<HcalQIESample> HcalCableMapper::SampleSet

private

Definition at line 23 of file HcalCableMapper.cc.

Constructor & Destructor Documentation

HcalCableMapper::HcalCableMapper ( edm::ParameterSet const &  conf )

explicit

Definition at line 57 of file HcalCableMapper.cc.

                                                            :
  hbheLabel_(conf.getParameter<edm::InputTag>("hbheLabel")),
  hoLabel_(conf.getParameter<edm::InputTag>("hoLabel")),
  hfLabel_(conf.getParameter<edm::InputTag>("hfLabel")){
  
}

Member Function Documentation

void HcalCableMapper::analyze ( edm::Event const &  e, edm::EventSetup const &  c  )

virtual

Implements edm::EDAnalyzer.

Definition at line 160 of file HcalCableMapper.cc.

References edm::Event::getByLabel(), hbheLabel_, hfLabel_, hoLabel_, and record().

                                                                       {
   
  edm::Handle<HBHEDigiCollection> hbhe;
  e.getByLabel(hbheLabel_,hbhe);

  edm::Handle<HFDigiCollection> hf;
  e.getByLabel(hfLabel_,hf);
  edm::Handle<HODigiCollection> ho;
  e.getByLabel(hoLabel_,ho);
  
   
  record(*hbhe);
  record(*hf);  
  record(*ho);
}

void HcalCableMapper::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 179 of file HcalCableMapper.cc.

References trackerHits::c, fullHistory_, i, IdSet, j, gen::k, GetRecoTauVFromDQM_MC_cff::kk, process(), and alignCSCRings::s.

                            {
  
  
  std::vector<SampleSet>::iterator j; 
  int c [128];
  int k,ii,kk;
  int c_max=0;
  
  std::map<HcalDetId,std::vector<SampleSet> >::iterator i;
  
  PathSet consensus;
 
  for (i=fullHistory_.begin(); i!=fullHistory_.end(); i++) {
    //i.first --> id
    //i.second --> vector<SampleSet>
    SampleSet s;
    for (k=0; k<10; k++) {
      for (ii=0; ii<128; ii++) c[ii]=0;

      for (j=i->second.begin();j!=i->second.end();j++){//word number
    if (int(j->size())>k) 
      c[(*j)[k].adc()]++;
    
    
      }//j loop
      //sort c-array
      for (kk=0;kk<128;kk++){  
    if (c[kk]>c[c_max]){
      c_max = kk;
    }
      }    
      
      s.push_back(((c_max&0x7F)));

      c_max=0;
    }//k-loop    
  consensus[i->first]=s;
  
  }//i loop
 
  process(consensus,IdSet);
    
   



}//end of endjob 

void HcalCableMapper::process ( const PathSet &  ps, const IdMap &  im  )

private

Definition at line 68 of file HcalCableMapper.cc.

References trackerHits::c, gather_cfg::cout, HcalElectronicsId::dccid(), createTree::dd, HcalDetId::depth(), det_names, HcalElectronicsId::fiberChanId(), HcalElectronicsId::fiberIndex(), errorMatrix2Lands::header, HcalElectronicsId::htrSlot(), HcalElectronicsId::htrTopBottom(), HcalDetId::ietaAbs(), HcalDetId::iphi(), HcalElectronicsId::readoutVMECrateId(), HcalElectronicsId::spigot(), and HcalDetId::subdet().

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::dumpPython(), endJob(), ConfigBuilder.ConfigBuilder.PrintAllModules::leave(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::open(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::setProcess(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::setProperty().

                                                               {
  
  PathSet::const_iterator iii;
  IdMap::const_iterator ij;
  
  for (iii=ps.begin();iii!=ps.end();iii++){
    
    
    SampleSet ss = iii->second;
    const HcalDetId dd = iii->first;
    
    ij=im.find(dd);
    HcalElectronicsId eid=ij->second;

    int header = ((ss[0].adc())&0x7F);
    int ieta  = ((ss[1].adc())&0x3F);
    int z_ieta = (((ss[1].adc())>>6)&0x1);
    int iphi = ((ss[2].adc())&0x7F);
    int depth = ((ss[3].adc())&0x7);
    int det = (((ss[3].adc())>>3)&0xF);
    int spigot = ((ss[4].adc())&0xF);
    int fiber = (((ss[4].adc())>>4)&0x7);
    int crate  = ((ss[5].adc())&0x1F);
    int fiber_chan= (((ss[5].adc())>>5)&0x3);
    int G_Dcc  = ((ss[6].adc())&0x3F);
    int H_slot  = ((ss[7].adc())&0x1F);
    int TB = (((ss[7].adc())>>5)&0x1);
    int RBX_7 = (((ss[7].adc())>>6)&0x1);
    int RBX = ((ss[8].adc())&0x7F);
    int RM = ((ss[9].adc())&0x3);
    int RM_card= (((ss[9].adc())>>2)&0x3);
    int RM_chan= (((ss[9].adc())>>4)&0x7);
    string eta_sign;
    std::string det_name;
    if (det>4 || det<0) {
      char c[20];
      snprintf(c,20,"Det=%d",det);
      det_name=c;
    } else det_name=det_names[det];
    
    
    if (z_ieta==1){
      eta_sign = "+";
    }else{eta_sign = "-";}
    string is_header;
    if (header == 0x75){
     
      //NO SHIFT
      if((spigot==eid.spigot())&&(fiber+1==eid.fiberIndex())&&(fiber_chan==eid.fiberChanId())&&(H_slot==eid.htrSlot())&&(G_Dcc==eid.dccid())&&(crate==eid.readoutVMECrateId())&&(iphi==dd.iphi())&&(depth==dd.depth())&&(ieta==dd.ietaAbs())&&(TB==eid.htrTopBottom())&&(det==dd.subdet())){//&&(z_ieta==dd.zside())
    std::cout <<"Pathway match"<<std::endl;
      }else{
    
    is_header=" Header found";
    
    std::cout <<" Digi ID: " << dd << is_header<< " ieta: "<< eta_sign << ieta << " iphi: "<< iphi << " Depth: " << depth << " Detector: " << det_name << " Spigot: "<< spigot<<"/"<<eid.spigot() << " Fiber: " << fiber+1<<"/"<<eid.fiberIndex() << " Fiber Channel: "<< fiber_chan <<"/"<<eid.fiberChanId()<< " Crate: " << crate<<"/"<<eid.readoutVMECrateId() << " Global Dcc: " << G_Dcc <<"/"<<eid.dccid() << " HTR Slot: " << H_slot <<"/ " <<eid.htrSlot()<< " Top/Bottom: " << TB<<"/"<< eid.htrTopBottom() << " RBX: " << (RBX_7*128+RBX) << " RM: " << RM+1 << " RM Card: " << RM_card+1 << " RM Channel: " << RM_chan << std::endl;
      }    
    }else if (ieta+64==0x75){
      
      ieta  = ((ss[2].adc())&0x3F);
      z_ieta = (((ss[2].adc())>>6)&0x1);
      iphi = ((ss[3].adc())&0x7F);
      depth = ((ss[4].adc())&0x7);
      det = (((ss[4].adc())>>3)&0xF);
      spigot = ((ss[5].adc())&0xF);
      fiber = (((ss[5].adc())>>4)&0x7);
      crate  = ((ss[6].adc())&0x1F);
      fiber_chan= (((ss[6].adc())>>5)&0x3);
      G_Dcc  = ((ss[7].adc())&0x3F);
      H_slot  = ((ss[8].adc())&0x1F);
      TB = (((ss[8].adc())>>5)&0x1);
      RBX_7 = (((ss[8].adc())>>6)&0x1);
      RBX = ((ss[9].adc())&0x7F);
      

      //SHIFT
        if((spigot==eid.spigot())&&(fiber+1==eid.fiberIndex())&&(fiber_chan==eid.fiberChanId())&&(H_slot==eid.htrSlot())&&(G_Dcc==eid.dccid())&&(TB==eid.htrTopBottom())&&(crate==eid.readoutVMECrateId())&&(iphi==dd.iphi())&&(depth==dd.depth())&&(det==dd.subdet())&&(ieta==dd.ietaAbs())){//&&(z_ieta==dd.zside())

    std::cout <<"Pathway match (SHIFT)"<<std::endl;
      }else{


      is_header=" DATA SHIFT";
      
     std::cout <<" Digi ID: " << dd << is_header<< " ieta: "<< eta_sign << ieta << " iphi: "<< iphi << " Depth: " << depth << " Detector: " << det_name << " Spigot: "<< spigot<<"/"<<eid.spigot() << " Fiber: " << fiber+1<<"/"<<eid.fiberIndex() << " Fiber Channel: "<< fiber_chan <<"/"<<eid.fiberChanId()<< " Crate: " << crate<<"/"<<eid.readoutVMECrateId() << " Global Dcc: " << G_Dcc <<"/"<<eid.dccid() << " HTR Slot: " << H_slot <<"/ " <<eid.htrSlot()<< " Top/Bottom: " << TB<<"/"<< eid.htrTopBottom() << " RBX: " << (RBX_7*128+RBX) << std::endl;
      
    }
    }else { std::cout<<" Digi ID: " <<dd << " +NO HEADER+  " << " RBX: " << (RBX_7*128+RBX) << std::endl;
    }
  }
}

template<class DigiCollection >

void HcalCableMapper::record ( const DigiCollection &  digis )

inlineprivate

Definition at line 38 of file HcalCableMapper.cc.

References i, and lumiQueryAPI::q.

Referenced by analyze().

                                           {

    for (typename DigiCollection::const_iterator digi=digis.begin(); digi!=digis.end(); digi++) {

      SampleSet q;
      for (int i=0; i<digi->size(); i++)
    q.push_back(digi->sample(i));
      
      if (fullHistory_.find(digi->id())==fullHistory_.end()) fullHistory_.insert(std::pair<HcalDetId,std::vector<SampleSet> >(digi->id(),std::vector<SampleSet>()));
      if (IdSet.find(digi->id())==IdSet.end()) IdSet.insert(std::pair<HcalDetId,HcalElectronicsId>(digi->id(),HcalElectronicsId())); 
      fullHistory_[digi->id()].push_back(q);
      IdSet[digi->id()]=digi->elecId();
    }
  }

Member Data Documentation

std::map<HcalDetId,std::vector<SampleSet> > HcalCableMapper::fullHistory_

private

Definition at line 33 of file HcalCableMapper.cc.

Referenced by endJob().

edm::InputTag HcalCableMapper::hbheLabel_

private

Definition at line 35 of file HcalCableMapper.cc.

Referenced by analyze().

edm::InputTag HcalCableMapper::hfLabel_

private

Definition at line 35 of file HcalCableMapper.cc.

Referenced by analyze().

edm::InputTag HcalCableMapper::hoLabel_

private

Definition at line 35 of file HcalCableMapper.cc.

Referenced by analyze().

IdMap HcalCableMapper::IdSet

private

Definition at line 34 of file HcalCableMapper.cc.

Referenced by endJob().