HcalCableMapper Class Reference

Inheritance diagram for HcalCableMapper:

Public Member Functions
virtual void	analyze (edm::Event const &e, edm::EventSetup const &c) override
virtual void	endJob () override
	HcalCableMapper (edm::ParameterSet const &conf)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
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
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::vector< ModuleDescription const * > &modules, 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
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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_
IdMap	IdSet
edm::EDGetTokenT< HBHEDigiCollection >	tok_hbhe_
edm::EDGetTokenT< HFDigiCollection >	tok_hf_
edm::EDGetTokenT< HODigiCollection >	tok_ho_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 15 of file HcalCableMapper.cc.

Member Typedef Documentation

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

private

Definition at line 25 of file HcalCableMapper.cc.

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

private

Definition at line 24 of file HcalCableMapper.cc.

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

private

Definition at line 22 of file HcalCableMapper.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 58 of file HcalCableMapper.cc.

References edm::ParameterSet::getParameter().

                                                            {
  tok_hbhe_ = consumes<HBHEDigiCollection>(conf.getParameter<edm::InputTag>("hbheLabel"));
  tok_ho_ = consumes<HODigiCollection>(conf.getParameter<edm::InputTag>("hoLabel"));
  tok_hf_ = consumes<HFDigiCollection>(conf.getParameter<edm::InputTag>("hfLabel"));
  
}

Member Function Documentation

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

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 161 of file HcalCableMapper.cc.

References edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, photonIsolationHIProducer_cfi::hf, photonIsolationHIProducer_cfi::ho, and record.

                                                                       {
   
  edm::Handle<HBHEDigiCollection> hbhe;
  e.getByToken(tok_hbhe_,hbhe);

  edm::Handle<HFDigiCollection> hf;
  e.getByToken(tok_hf_,hf);
  edm::Handle<HODigiCollection> ho;
  e.getByToken(tok_ho_,ho);
  
   
  record(*hbhe);
  record(*hf);  
  record(*ho);
}

void HcalCableMapper::endJob ( void  )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 180 of file HcalCableMapper.cc.

References EnergyCorrector::c, DEFINE_FWK_MODULE, mps_fire::i, cuy::ii, gen::k, GetRecoTauVFromDQM_MC_cff::kk, LaserDQM_cfg::process, and alignCSCRings::s.

Referenced by o2olib.O2ORunMgr::executeJob().

                            {
  
  
  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 69 of file HcalCableMapper.cc.

References EnergyCorrector::c, gather_cfg::cout, HcalElectronicsId::dccid(), createTree::dd, particleFlowClusterECALTimeSelected_cfi::depth, HcalDetId::depth(), det_names, runTauDisplay::eid, HcalElectronicsId::fiberChanId(), HcalElectronicsId::fiberIndex(), RecoTauValidation_cfi::header, HcalElectronicsId::htrSlot(), HcalElectronicsId::htrTopBottom(), HcalDetId::ietaAbs(), HcalDetId::iphi(), HcalElectronicsId::readoutVMECrateId(), HcalElectronicsId::spigot(), AlCaHLTBitMon_QueryRunRegistry::string, and HcalDetId::subdet().

Referenced by ConfigBuilder.ConfigBuilder::addExtraStream(), ConfigBuilder.ConfigBuilder::completeInputCommand(), ConfigBuilder.ConfigBuilder::doNotInlineEventContent(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::dumpPython(), ConfigBuilder.ConfigBuilder.PrintAllModules::leave(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::open(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), ConfigBuilder.ConfigBuilder::prepare(), ConfigBuilder.ConfigBuilder::prepare_ALCA(), ConfigBuilder.ConfigBuilder::prepare_DQM(), ConfigBuilder.ConfigBuilder::prepare_HLT(), ConfigBuilder.ConfigBuilder::prepare_LHE(), ConfigBuilder.ConfigBuilder::prepare_PATFILTER(), ConfigBuilder.ConfigBuilder::prepare_VALIDATION(), ConfigBuilder.ConfigBuilder::renameHLTprocessInSequence(), ConfigBuilder.ConfigBuilder::renameInputTagsInSequence(), ConfigBuilder.ConfigBuilder::scheduleSequence(), 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 39 of file HcalCableMapper.cc.

References mps_fire::i, and lumiQueryAPI::q.

                                           {

    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 32 of file HcalCableMapper.cc.

IdMap HcalCableMapper::IdSet

private

Definition at line 33 of file HcalCableMapper.cc.

edm::EDGetTokenT<HBHEDigiCollection> HcalCableMapper::tok_hbhe_

private

Definition at line 34 of file HcalCableMapper.cc.

edm::EDGetTokenT<HFDigiCollection> HcalCableMapper::tok_hf_

private

Definition at line 36 of file HcalCableMapper.cc.

edm::EDGetTokenT<HODigiCollection> HcalCableMapper::tok_ho_

private

Definition at line 35 of file HcalCableMapper.cc.