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

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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 HLT_25ns14e33_v1_cff::InputTag, 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.

265 {}

Member Function Documentation

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

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 268 of file Digi2Raw2Digi.cc.

References 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(), timingPdfMaker::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, HLT_25ns14e33_v1_cff::depth, MonitorElement::Fill(), edm::Event::getByToken(), i, fjr2json::lumi, match(), meStatus, unsuppressed, and ecaldqm::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;
   }
 }