#include <Digi2Raw2Digi.h>

Inheritance diagram for Digi2Raw2Digi:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob ()

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 &)

virtual void	endJob ()

	~Digi2Raw2Digi ()

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
	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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
DQMStore *	dbe_

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

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 16 of file Digi2Raw2Digi.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 219 of file Digi2Raw2Digi.cc.

References DQMStore::book1D(), dbe_, timingPdfMaker::histo, inputTag1_, inputTag2_, meStatus, cppFunctionSkipper::operator, outputFile_, DQMStore::setCurrentFolder(), 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")),
     dbe_(0)
 {
 
   // 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";
   }
 
   // DQM service initialization
   dbe_ = edm::Service<DQMStore>().operator->();
    
   if ( dbe_ ) {
     dbe_->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    = dbe_->book1D(histo, histo, 5, 0., 5.);
 
 }

Digi2Raw2Digi::~Digi2Raw2Digi ( )

Definition at line 268 of file Digi2Raw2Digi.cc.

268 {}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 271 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::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 263 of file Digi2Raw2Digi.cc.

263 {}

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

void Digi2Raw2Digi::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 264 of file Digi2Raw2Digi.cc.

References dbe_, outputFile_, and DQMStore::save().

                            {
   if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
 }

Member Data Documentation

DQMStore* Digi2Raw2Digi::dbe_

private

Definition at line 39 of file Digi2Raw2Digi.h.

Referenced by Digi2Raw2Digi(), and endJob().

edm::InputTag Digi2Raw2Digi::inputTag1_

private

Definition at line 26 of file Digi2Raw2Digi.h.

Referenced by Digi2Raw2Digi().

edm::InputTag Digi2Raw2Digi::inputTag2_

private

Definition at line 27 of file Digi2Raw2Digi.h.

Referenced by Digi2Raw2Digi().

MonitorElement* Digi2Raw2Digi::meStatus

private

Definition at line 41 of file Digi2Raw2Digi.h.

Referenced by compare(), and Digi2Raw2Digi().

std::string Digi2Raw2Digi::outputFile_

private

Definition at line 38 of file Digi2Raw2Digi.h.

Referenced by Digi2Raw2Digi(), and endJob().

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

private

Definition at line 29 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

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

private

Definition at line 30 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

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

private

Definition at line 33 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

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

private

Definition at line 34 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

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

private

Definition at line 31 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

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

private

Definition at line 32 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

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

private

Definition at line 35 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

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

private

Definition at line 36 of file Digi2Raw2Digi.h.

Referenced by analyze(), and Digi2Raw2Digi().

int Digi2Raw2Digi::unsuppressed

private

Definition at line 43 of file Digi2Raw2Digi.h.

Referenced by analyze(), and compare().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation