Digi2Raw2Digi Class Reference

#include <Digi2Raw2Digi.h>

Inheritance diagram for Digi2Raw2Digi:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
template<class Digi >
void	compare (const edm::Event &, const edm::EventSetup &, const edm::EDGetTokenT< edm::SortedCollection< Digi > > &tok1, const edm::EDGetTokenT< edm::SortedCollection< Digi > > &tok2)
	Digi2Raw2Digi (const edm::ParameterSet &)
	~Digi2Raw2Digi ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	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 &	itemsToGetFromEvent () 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 Attributes
edm::InputTag	inputTag1_
edm::InputTag	inputTag2_
MonitorElement *	meStatus
std::string	outputFile_
edm::EDGetTokenT< edm::SortedCollection< HBHEDataFrame > >	tok_hbhe1_
edm::EDGetTokenT< edm::SortedCollection< HBHEDataFrame > >	tok_hbhe2_
edm::EDGetTokenT< edm::SortedCollection< HFDataFrame > >	tok_hf1_
edm::EDGetTokenT< edm::SortedCollection< HFDataFrame > >	tok_hf2_
edm::EDGetTokenT< edm::SortedCollection< HODataFrame > >	tok_ho1_
edm::EDGetTokenT< edm::SortedCollection< HODataFrame > >	tok_ho2_
edm::EDGetTokenT< edm::SortedCollection< ZDCDataFrame > >	tok_zdc1_
edm::EDGetTokenT< edm::SortedCollection< ZDCDataFrame > >	tok_zdc2_
int	unsuppressed

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 Digi2Raw2Digi.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 219 of file Digi2Raw2Digi.cc.

References inputTag1_, inputTag2_, outputFile_, tok_hbhe1_, tok_hbhe2_, tok_hf1_, tok_hf2_, tok_ho1_, tok_ho2_, tok_zdc1_, and tok_zdc2_.

  : inputTag1_(iConfig.getParameter<edm::InputTag>("digiLabel1")),
    inputTag2_(iConfig.getParameter<edm::InputTag>("digiLabel2")),
    outputFile_(iConfig.getUntrackedParameter<std::string>("outputFile"))
{

  // register for data access
  tok_hbhe1_ = consumes<edm::SortedCollection<HBHEDataFrame> >(inputTag1_);
  tok_hbhe2_ = consumes<edm::SortedCollection<HBHEDataFrame> >(inputTag2_);
  tok_ho1_ = consumes<edm::SortedCollection<HODataFrame> >(inputTag1_);
  tok_ho2_ = consumes<edm::SortedCollection<HODataFrame> >(inputTag2_);
  tok_hf1_ = consumes<edm::SortedCollection<HFDataFrame> >(inputTag1_);
  tok_hf2_ = consumes<edm::SortedCollection<HFDataFrame> >(inputTag2_);
  tok_zdc1_ = consumes<edm::SortedCollection<ZDCDataFrame> >(edm::InputTag("simHcalUnsuppressedDigis"));
  tok_zdc2_ = consumes<edm::SortedCollection<ZDCDataFrame> >(inputTag2_);
  

  // DQM ROOT output
  if ( outputFile_.size() != 0 ) {
    edm::LogInfo("OutputInfo")
      << " Hcal RecHit Task histograms will be saved to '" 
      << outputFile_.c_str() << "'";
  } else {
    edm::LogInfo("OutputInfo") 
      << " Hcal RecHit Task histograms will NOT be saved";
  }

   
}

Digi2Raw2Digi::~Digi2Raw2Digi

(

)

Definition at line 265 of file Digi2Raw2Digi.cc.

{}

Member Function Documentation

void Digi2Raw2Digi::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Definition at line 268 of file Digi2Raw2Digi.cc.

References DEFINE_FWK_MODULE, iEvent, tok_hbhe1_, tok_hbhe2_, tok_hf1_, tok_hf2_, tok_ho1_, tok_ho2_, tok_zdc1_, tok_zdc2_, and unsuppressed.

{
  unsuppressed = 0;
  
  //  std::cout << "=== HBHE ==================" << std::endl; 
  compare<HBHEDataFrame>(iEvent,iSetup,tok_hbhe1_,tok_hbhe2_);

  //  std::cout << "=== HO ====================" << std::endl; 
  compare<HODataFrame>(iEvent,iSetup,tok_ho1_,tok_ho2_);

  //  std::cout << "=== HF ====================" << std::endl; 
  compare<HFDataFrame>(iEvent,iSetup,tok_hf1_,tok_hf2_);  


  //  std::cout << "=== ZDC ===================" << std::endl; 
  unsuppressed = 1;
  compare<ZDCDataFrame>(iEvent,iSetup,tok_zdc1_,tok_zdc2_);
  
  
  //  std::cout << "=== CASTOR ================" << std::endl; 
  //  compare<CastorDataFrame>(iEvent,iSetup); 
  
}

void Digi2Raw2Digi::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  run, edm::EventSetup const &  es  )

virtual

Implements DQMEDAnalyzer.

Definition at line 249 of file Digi2Raw2Digi.cc.

References DQMStore::IBooker::book1D(), trackerHits::histo, meStatus, and DQMStore::IBooker::setCurrentFolder().

{

  ibooker.setCurrentFolder("Digi2Raw2DigiV/Digi2Raw2DigiTask");
 

  // const char * sub = hcalselector_.c_str();
  char histo[100];

  sprintf (histo, "Digi2Raw2Digi_status") ;
  // bins: 1)full match 2)ID match, not content 3) no match 
  // 4) number of events with diff number of Digis
  meStatus    = ibooker.book1D(histo, histo, 5, 0., 5.);

}

template<class Digi >

void Digi2Raw2Digi::compare	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const edm::EDGetTokenT< edm::SortedCollection< Digi > > &	tok1,
		const edm::EDGetTokenT< edm::SortedCollection< Digi > > &	tok2
	)

Definition at line 39 of file Digi2Raw2Digi.cc.

References ecalMGPA::adc(), gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, MonitorElement::Fill(), edm::Event::getByToken(), mps_fire::i, csvLumiCalc::lumi, match(), meStatus, vertexPlots::section, unsuppressed, and HcalDigiParam_cfi::zside.

                                                                                                                                                                                               {

  typename edm::Handle<edm::SortedCollection<Digi> > digiCollection1;
  typename edm::SortedCollection<Digi>::const_iterator digiItr1;
  typename edm::Handle<edm::SortedCollection<Digi> > digiCollection2;
  typename edm::SortedCollection<Digi>::const_iterator digiItr2;
  
  if(unsuppressed) {  // ZDC
     iEvent.getByToken (tok1, digiCollection1); 
  }
  else iEvent.getByToken (tok1, digiCollection1);
  
  iEvent.getByToken (tok2, digiCollection2);
  
  int size1 = 0;
  int size2 = 0;
  
  for (digiItr1=digiCollection1->begin();digiItr1!=digiCollection1->end();digiItr1++) {
    size1++;
  }

  for (digiItr2=digiCollection2->begin();digiItr2!=digiCollection2->end();digiItr2++) {
    size2++;
  }

  //std::cout << "Digi collections   size1 = "<< size1 
  //    << "   size2 = " << size2 << std::endl;
  

  // CYCLE over first DIGI collection ======================================
  
  for (digiItr1=digiCollection1->begin();digiItr1!=digiCollection1->end();digiItr1++) {
    HcalGenericDetId HcalGenDetId(digiItr1->id());
    int tsize =  (*digiItr1).size();
    int match = 0;

    if(HcalGenDetId.isHcalZDCDetId()){
      //for zdc
      HcalZDCDetId element(digiItr1->id());
      int zside =  element.zside();
      int section = element.section();
      int channel = element.channel();
      int gsub = HcalGenDetId.genericSubdet();
 
    if(section==3){// lumi section not reconstructed
      size2++;
      match = 1;
      goto lumi;
    }
      

    //std::cout<< " Zdc genSubdet="<< gsub << " zside=" <<zside
    //       << " section= "<< section << " channel " <<channel
    //      <<std::endl; 

      for (digiItr2=digiCollection2->begin();digiItr2!=digiCollection2->end();digiItr2++) {
    HcalZDCDetId element2(digiItr2->id());
    
    //int zside2 =  element2.zside();
    //int section2 = element2.section();
    //int channel2 = element2.channel();
    //int gsub2 = HcalGenDetId.genericSubdet();

    //std::cout<< " Zdc genSubdet="<<gsub2 
    //   << " zside=" <<zside2
    //   << " section= "<<section2
    //   <<" channel "<<channel2 
    //   <<std::endl; 

    if(element == element2) {
      match = 1;
      int identical = 1; 
      for (int i=0; i<tsize; i++) {
        double adc  =  (*digiItr1)[i].adc();     
        int capid   =  (*digiItr1)[i].capid();
        //    std::cout << std::endl << "     capid1=" << capid 
        //                << " adc1=" << adc 
        //          << std::endl;
        double adc2 =  (*digiItr2)[i].adc();     
        int capid2  =  (*digiItr2)[i].capid();
        //    std::cout << "     capid2=" << capid2 << " adc2=" << adc2 
        //          << std::endl;
        if( capid != capid2 || adc !=  adc2) {
          std::cout << "===> PROBLEM !!!  gebsubdet=" << gsub 
            << " zside=" <<zside
            << " section= "<< section << " channel " <<channel
            << std::endl;
          std::cout << "     capid1["<< i << "]=" << capid 
            << " adc1["<< i << "]=" << adc 
            << "     capid2["<< i << "]=" << capid2 
            << " adc2["<< i << "]=" << adc2
            << std::endl;   
          identical = 0;
          meStatus->Fill(1.); 
          break;
        }
      } // end of DataFrames array
      if(identical) meStatus->Fill(0.); 
      break; // matched HcalZDCID  is processed,
      //  go to next (primary collection) cell  
    } 
      } // end of cycle over 2d DIGI collection 
    lumi:
      if (!match) {
    meStatus->Fill(2.); 
    std::cout << "===> PROBLEM !!!  gsubdet=" << gsub
          << " zside=" <<zside
          << " section= "<< section << " channel " <<channel
          << " HcalZDCId match is not found !!!" 
          << std::endl;
      }
  
    }
    else{
      //for Hcal subdetectors
      HcalDetId cell(digiItr1->id()); 
      int depth = cell.depth();
      int iphi  = cell.iphi()-1;
      int ieta  = cell.ieta();
      int sub   = cell.subdet();
      //    if(ieta > 0) ieta--;
      //    std::cout << " Cell subdet=" << sub << "  ieta=" << ieta 
      //          << "  inphi=" << iphi << "  depth=" << depth << std::endl;
      
      // CYCLE over second DIGI collection ======================================
      for (digiItr2=digiCollection2->begin();digiItr2!=digiCollection2->end();digiItr2++) {
  
    HcalDetId cell2(digiItr2->id());
    
    if( cell == cell2) {
      match = 1;
      int identical = 1; 
      for (int i=0; i<tsize; i++) {
        double adc  =  (*digiItr1)[i].adc();     
        int capid   =  (*digiItr1)[i].capid();
        //    std::cout << std::endl << "     capid1=" << capid 
        //                << " adc1=" << adc 
        //          << std::endl;
        double adc2 =  (*digiItr2)[i].adc();     
        int capid2  =  (*digiItr2)[i].capid();
        //    std::cout << "     capid2=" << capid2 << " adc2=" << adc2 
        //          << std::endl;
        if( capid != capid2 || adc !=  adc2) {
          std::cout << "===> PROBLEM !!!  subdet=" << sub << "  ieta="
            << ieta  << "  inphi=" << iphi << "  depth=" << depth
            << std::endl;
          std::cout << "     capid1["<< i << "]=" << capid 
            << " adc1["<< i << "]=" << adc 
            << "     capid2["<< i << "]=" << capid2 
            << " adc2["<< i << "]=" << adc2
            << std::endl;   
            identical = 0;
            meStatus->Fill(1.); 
            break;
        }
      } // end of DataFrames array
      if(identical) meStatus->Fill(0.); 
      break; // matched HcalID  is processed,
      //  go to next (primary collection) cell  
    } 
      } // end of cycle over 2d DIGI collection 
      if (!match) {
    meStatus->Fill(2.); 
    std::cout << "===> PROBLEM !!!  subdet=" << sub << "  ieta="
          << ieta  << "  inphi=" << iphi << "  depth=" << depth
          << " HcalID match is not found !!!" 
          << std::endl;
      }
    }
  } // end of cycle over 1st DIGI collection    
  
  if (size1 != size2) {
    meStatus->Fill(3.); 
    std::cout << "===> PROBLEM !!!  Different size of Digi collections : "
          << size1  << " and " << size2 
          << std::endl;
  }
}

Member Data Documentation

edm::InputTag Digi2Raw2Digi::inputTag1_

private

Definition at line 29 of file Digi2Raw2Digi.h.

Referenced by Digi2Raw2Digi().

edm::InputTag Digi2Raw2Digi::inputTag2_

private

Definition at line 30 of file Digi2Raw2Digi.h.

Referenced by Digi2Raw2Digi().

MonitorElement* Digi2Raw2Digi::meStatus

private

Definition at line 43 of file Digi2Raw2Digi.h.

Referenced by bookHistograms(), and compare().

std::string Digi2Raw2Digi::outputFile_

private

Definition at line 41 of file Digi2Raw2Digi.h.

Referenced by Digi2Raw2Digi().

edm::EDGetTokenT<edm::SortedCollection<HBHEDataFrame> > Digi2Raw2Digi::tok_hbhe1_

private

Definition at line 32 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

edm::EDGetTokenT<edm::SortedCollection<HBHEDataFrame> > Digi2Raw2Digi::tok_hbhe2_

private

Definition at line 33 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

edm::EDGetTokenT<edm::SortedCollection<HFDataFrame> > Digi2Raw2Digi::tok_hf1_

private

Definition at line 36 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

edm::EDGetTokenT<edm::SortedCollection<HFDataFrame> > Digi2Raw2Digi::tok_hf2_

private

Definition at line 37 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

edm::EDGetTokenT<edm::SortedCollection<HODataFrame> > Digi2Raw2Digi::tok_ho1_

private

Definition at line 34 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

edm::EDGetTokenT<edm::SortedCollection<HODataFrame> > Digi2Raw2Digi::tok_ho2_

private

Definition at line 35 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

edm::EDGetTokenT<edm::SortedCollection<ZDCDataFrame> > Digi2Raw2Digi::tok_zdc1_

private

Definition at line 38 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

edm::EDGetTokenT<edm::SortedCollection<ZDCDataFrame> > Digi2Raw2Digi::tok_zdc2_

private

Definition at line 39 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

int Digi2Raw2Digi::unsuppressed

private

Definition at line 45 of file Digi2Raw2Digi.h.

Referenced by analyze(), and compare().