#include <HcalGainsCheck.h>

Inheritance diagram for HcalGainsCheck:

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

virtual void	beginJob ()

virtual void	endJob ()

	HcalGainsCheck (edm::ParameterSet const &ps)

	~HcalGainsCheck ()

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 Attributes
TH1F *	diffUpRefCap0

TH1F *	diffUpRefCap1

TH1F *	diffUpRefCap2

TH1F *	diffUpRefCap3

std::string	dumprefs

std::string	dumpupdate

bool	emapflag

double	epsilon

TFile *	f

TH1F *	gainsRefCap0

TH1F *	gainsRefCap0vsEta

TH1F *	gainsRefCap1

TH1F *	gainsRefCap1vsEta

TH1F *	gainsRefCap2

TH1F *	gainsRefCap2vsEta

TH1F *	gainsRefCap3

TH1F *	gainsRefCap3vsEta

TH1F *	gainsUpCap0

TH1F *	gainsUpCap0vsEta

TH1F *	gainsUpCap1

TH1F *	gainsUpCap1vsEta

TH1F *	gainsUpCap2

TH1F *	gainsUpCap2vsEta

TH1F *	gainsUpCap3

TH1F *	gainsUpCap3vsEta

TH2F *	ocMapRef

TH2F *	ocMapUp

std::string	outfile

TH1F *	ratioUpRefCap0

TH1F *	ratioUpRefCap1

TH1F *	ratioUpRefCap2

TH1F *	ratioUpRefCap3

std::string	rootfile

bool	validategainsflag

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 47 of file HcalGainsCheck.h.

Constructor & Destructor Documentation

HcalGainsCheck::HcalGainsCheck ( edm::ParameterSet const & ps )

Definition at line 3 of file HcalGainsCheck.cc.

References dumprefs, dumpupdate, emapflag, epsilon, edm::ParameterSet::getUntrackedParameter(), outfile, rootfile, AlCaHLTBitMon_QueryRunRegistry::string, and validategainsflag.

 {
   rootfile = ps.getUntrackedParameter<std::string>("rootfile","null");
   outfile = ps.getUntrackedParameter<std::string>("outFile","null");
   dumpupdate = ps.getUntrackedParameter<std::string>("dumpUpdateGainsTo","null");
   dumprefs = ps.getUntrackedParameter<std::string>("dumpRefGainsTo","null");
   emapflag = ps.getUntrackedParameter<bool>("checkEmap",false);
   validategainsflag = ps.getUntrackedParameter<bool>("validateGains",false);
   epsilon = ps.getUntrackedParameter<double>("deltaG",1000000);
 }

HcalGainsCheck::~HcalGainsCheck ( )

inline

Definition at line 52 of file HcalGainsCheck.h.

References analyze(), beginJob(), endJob(), and ev.

52 {}

Member Function Documentation

void HcalGainsCheck::analyze	(	const edm::Event &	ev,
		const edm::EventSetup &	es
	)

Definition at line 49 of file HcalGainsCheck.cc.

References HcalCondObjectContainer< Item >::addValues(), HcalElectronicsMap::allPrecisionId(), gather_cfg::cout, diffUpRefCap0, diffUpRefCap1, diffUpRefCap2, diffUpRefCap3, HcalDbASCIIIO::dumpObject(), dumprefs, dumpupdate, emapflag, epsilon, Exception, spr::find(), gainsRefCap0, gainsRefCap1, gainsRefCap2, gainsRefCap3, gainsUpCap0, gainsUpCap1, gainsUpCap2, gainsUpCap3, edm::EventSetup::get(), HcalCondObjectContainer< Item >::getAllChannels(), HcalGain::getValue(), HcalCondObjectContainer< Item >::getValues(), outfile, edm::ESHandle< T >::product(), ratioUpRefCap0, ratioUpRefCap1, ratioUpRefCap2, ratioUpRefCap3, HcalCondObjectContainerBase::topo(), and validategainsflag.

Referenced by ~HcalGainsCheck().

 {
   using namespace edm::eventsetup;
   bool epsilonflag = false;
   bool notequalsflag = false;
   // get new gains
   edm::ESHandle<HcalGains> newGains;
   es.get<HcalGainsRcd>().get("update",newGains);
   const HcalGains* myNewGains = newGains.product();
 
   // get reference gains
   edm::ESHandle<HcalGains> refGains;
   es.get<HcalGainsRcd>().get("reference",refGains);
   const HcalGains* myRefGains = refGains.product();
 
   // get e-map from reference
   edm::ESHandle<HcalElectronicsMap> refEMap;
   es.get<HcalElectronicsMapRcd>().get("reference",refEMap);
   const HcalElectronicsMap* myRefEMap = refEMap.product();
 
 
   // dump gains:
    if(dumpupdate.compare("null")!=0){
    std::ofstream outStream(dumpupdate.c_str());
    std::cout << "--- Dumping Gains - update ---" << std::endl;
    HcalDbASCIIIO::dumpObject (outStream, (*myNewGains) );
    }
    if(dumprefs.compare("null")!=0){
     std::ofstream outStream2(dumprefs.c_str());
     std::cout << "--- Dumping Gains - reference ---" << std::endl;
     HcalDbASCIIIO::dumpObject (outStream2, (*myRefGains) );
    }
     // get the list of all channels
     std::vector<DetId> listNewChan = myNewGains->getAllChannels();
     std::vector<DetId> listRefChan = myRefGains->getAllChannels();
     
     std::vector<DetId>::const_iterator cell;
 
     //plots: occupancy map, value map, difference, ratio, gains:
     for (std::vector<DetId>::const_iterator it = listRefChan.begin(); it!=listRefChan.end(); it++)
       {
     HcalGenericDetId myId(*it);
     //  ocMapRef->Fill(myId->);
 
     float valCap0 = myRefGains->getValues(*it)->getValue(0);
     float valCap1 = myRefGains->getValues(*it)->getValue(1);
     float valCap2 = myRefGains->getValues(*it)->getValue(2);
     float valCap3 = myRefGains->getValues(*it)->getValue(3);
 
     gainsRefCap0->Fill(valCap0);
     gainsRefCap1->Fill(valCap1);
     gainsRefCap2->Fill(valCap2);
     gainsRefCap3->Fill(valCap3);
 
     cell = std::find(listNewChan.begin(), listNewChan.end(), (*it));
     if (cell != listNewChan.end() ) //found
       {
         float valCap0up = myNewGains->getValues(*it)->getValue(0);
         float valCap1up = myNewGains->getValues(*it)->getValue(1);
         float valCap2up = myNewGains->getValues(*it)->getValue(2);
         float valCap3up = myNewGains->getValues(*it)->getValue(3);
         
         diffUpRefCap0->Fill(valCap0up - valCap0);
         diffUpRefCap1->Fill(valCap1up - valCap1);
         diffUpRefCap2->Fill(valCap2up - valCap2);
         diffUpRefCap3->Fill(valCap3up - valCap3);
 
         if(fabs(valCap0up - valCap0) > epsilon) epsilonflag = true;
             if(fabs(valCap1up - valCap1) > epsilon) epsilonflag = true;
             if(fabs(valCap2up - valCap2) > epsilon) epsilonflag = true;
             if(fabs(valCap3up - valCap3) > epsilon) epsilonflag = true;
 
             if(valCap0up != valCap0) notequalsflag = true;
             if(valCap1up != valCap1) notequalsflag = true;
             if(valCap2up != valCap2) notequalsflag = true;
             if(valCap3up != valCap3) notequalsflag = true;
 
         ratioUpRefCap0->Fill(valCap0up / valCap0);
         ratioUpRefCap1->Fill(valCap1up / valCap1);
         ratioUpRefCap2->Fill(valCap2up / valCap2);
         ratioUpRefCap3->Fill(valCap3up / valCap3);
       }
       }
     for (std::vector<DetId>::const_iterator it = listNewChan.begin(); it!=listNewChan.end(); it++)
       {
     float valCap0 = myNewGains->getValues(*it)->getValue(0);
     float valCap1 = myNewGains->getValues(*it)->getValue(1);
     float valCap2 = myNewGains->getValues(*it)->getValue(2);
     float valCap3 = myNewGains->getValues(*it)->getValue(3);
 
     gainsUpCap0->Fill(valCap0);
     gainsUpCap1->Fill(valCap1);
     gainsUpCap2->Fill(valCap2);
     gainsUpCap3->Fill(valCap3);
       }
     
     if(epsilon != 1000000){
        if(epsilonflag) throw cms::Exception("DataDoesNotMatch") << "Values differ by more than deltaG" << std::endl;
     }else{
        std::cout << "These gains do not differ by more than deltaG" << std::endl;
     }
 
     if(validategainsflag){
        if(notequalsflag) throw cms::Exception("DataDoesNotMatch") << "Values do not match" << std::endl;
     }else{
        std::cout << "These gains are identical" << std::endl;
     }
 
     // go through list of valid channels from reference, look up if conditions exist for update
     // push back into new vector the corresponding updated conditions,
     // or if it doesn't exist, the reference
 
     if(outfile.compare("null")!=0){
       HcalGains *resultGains = new HcalGains(refGains->topo());
     for (std::vector<DetId>::const_iterator it = listRefChan.begin(); it != listRefChan.end(); it++)
       {
     DetId mydetid = *it;
     HcalGenericDetId myId(*it);
     cell = std::find(listNewChan.begin(), listNewChan.end(), mydetid);
     if (cell == listNewChan.end()) // not present in new list, take old conditions
       {
         const HcalGain* item = myRefGains->getValues(mydetid);
         std::cout << "o";
         resultGains->addValues(*item);
       }
     else // present in new list, take new conditions
       {
         const HcalGain* item = myNewGains->getValues(mydetid);
         std::cout << "n";
         resultGains->addValues(*item);
       }
       }
     std::cout << std::endl;
     
     std::vector<DetId> listResult = resultGains->getAllChannels();
     // get the e-map list of channels
     if(emapflag){
     std::vector<HcalGenericDetId> listEMap = myRefEMap->allPrecisionId();
     // look up if emap channels are all present in pedestals, if not then cerr
     for (std::vector<HcalGenericDetId>::const_iterator it = listEMap.begin(); it != listEMap.end(); it++)
       {
       DetId mydetid = DetId(it->rawId());
     if (std::find(listResult.begin(), listResult.end(), mydetid ) == listResult.end()  )
       {
         std::cout << "Conditions not found for DetId = " << HcalGenericDetId(it->rawId()) << std::endl;
       }
       }
     }
    
     // dump the resulting list of pedestals into a file
     //    std::ostringstream filename3;
     //    filename3 << "test_combined.txt";
     std::ofstream outStream3(outfile.c_str());
     std::cout << "--- Dumping Gains - the combined ones ---" << std::endl;
     HcalDbASCIIIO::dumpObject (outStream3, (*resultGains) );
 
     }
     // const float* values = myped->getValues (channelID);
     //    if (values) std::cout << "pedestals for channel " << channelID << ": "
     //            << values [0] << '/' << values [1] << '/' << values [2] << '/' << values [3] << std::endl; 
 
 }

void HcalGainsCheck::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 14 of file HcalGainsCheck.cc.

References diffUpRefCap0, diffUpRefCap1, diffUpRefCap2, diffUpRefCap3, f, gainsRefCap0, gainsRefCap1, gainsRefCap2, gainsRefCap3, gainsUpCap0, gainsUpCap1, gainsUpCap2, gainsUpCap3, ocMapRef, ocMapUp, ratioUpRefCap0, ratioUpRefCap1, ratioUpRefCap2, ratioUpRefCap3, and rootfile.

Referenced by ~HcalGainsCheck().

 {
   f = new TFile(rootfile.c_str(),"RECREATE");
 
   //book histos:
   ocMapUp = new TH2F("ocMapUp","occupancy_map_updated_gains",83,-41.5,41.5,72,0.5,72.5);
   ocMapRef = new TH2F("ocMapRef","occupancy_map_updated_gains",83,-41.5,41.5,72,0.5,72.5);
 //  valMapUp;
 //  valMapRef;
 
   diffUpRefCap0 = new TH1F("diffUpRefCap0","difference_update_reference_Cap0",100,-0.5,0.5);
   ratioUpRefCap0 = new TH1F("ratioUpRefCap0", "ration_update_reference_Cap0",100,0.5,1.5);
   gainsUpCap0 = new TH1F("gainsUpCap0","gains_update_Cap0",100,0.0,0.6);
   gainsRefCap0 = new TH1F("gainsRefCap0","gains_reference_Cap0",100,0.0,0.6);
 
   diffUpRefCap1 = new TH1F("diffUpRefCap1","difference_update_reference_Cap1",100,-0.5,0.5);
   ratioUpRefCap1 = new TH1F("ratioUpRefCap1", "ration_update_reference_Cap1",100,0.5,1.5);
   gainsUpCap1 = new TH1F("gainsUpCap1","gains_update_Cap1",100,0.0,0.6);
   gainsRefCap1 = new TH1F("gainsRefCap1","gains_reference_Cap1",100,0.0,0.6);
 
   diffUpRefCap2 = new TH1F("diffUpRefCap2","difference_update_reference_Cap2",100,-0.5,0.5);
   ratioUpRefCap2 = new TH1F("ratioUpRefCap2", "ration_update_reference_Cap2",100,0.5,1.5);
   gainsUpCap2 = new TH1F("gainsUpCap2","gains_update_Cap2",100,0.0,0.6);
   gainsRefCap2 = new TH1F("gainsRefCap2","gains_reference_Cap2",100,0.0,0.6);
 
   diffUpRefCap3 = new TH1F("diffUpRefCap3","difference_update_reference_Cap3",100,-0.5,0.5);
   ratioUpRefCap3 = new TH1F("ratioUpRefCap3", "ration_update_reference_Cap3",100,0.5,1.5);
   gainsUpCap3 = new TH1F("gainsUpCap3","gains_update_Cap3",100,0.0,0.6);
   gainsRefCap3 = new TH1F("gainsRefCap3","gains_reference_Cap3",100,0.0,0.6);
 
   //  gainsUpCap0vsEta = new TGraph("gainsUpCap0vsEta","gains_update_Cap0_vsEta",100,-41,0.6);
   //  gainsRefCap0vsEta = new TGraph("gainsRefCap0vsEta","gains_reference_Cap0_vsEta",100,0.0,0.6);
 }

void HcalGainsCheck::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 215 of file HcalGainsCheck.cc.

References DEFINE_FWK_MODULE, f, and rootfile.

Referenced by ~HcalGainsCheck().

 {
   if(rootfile.compare("null")!=0){
   f->Write();
 }
   f->Close();
 
 }

Member Data Documentation

TH1F* HcalGainsCheck::diffUpRefCap0

private

Definition at line 77 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::diffUpRefCap1

private

Definition at line 78 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::diffUpRefCap2

private

Definition at line 79 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::diffUpRefCap3

private

Definition at line 80 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

std::string HcalGainsCheck::dumprefs

private

Definition at line 63 of file HcalGainsCheck.h.

Referenced by analyze(), and HcalGainsCheck().

std::string HcalGainsCheck::dumpupdate

private

Definition at line 62 of file HcalGainsCheck.h.

Referenced by analyze(), and HcalGainsCheck().

bool HcalGainsCheck::emapflag

private

Definition at line 66 of file HcalGainsCheck.h.

Referenced by analyze(), and HcalGainsCheck().

double HcalGainsCheck::epsilon

private

Definition at line 68 of file HcalGainsCheck.h.

Referenced by analyze(), utils.BinToBin::checkBinningMatches(), utils.BinToBin1percent::checkBinningMatches(), utils.BinToBin::do_test(), and HcalGainsCheck().

TFile* HcalGainsCheck::f

private

Definition at line 70 of file HcalGainsCheck.h.

Referenced by beginJob(), Vispa.Views.RootCanvasView.RootCanvasView::createGraph(), Vispa.Views.RootCanvasView.RootCanvasView::createLegoPlot(), endJob(), and ztail.Decoder::initial_synchronize().

TH1F* HcalGainsCheck::gainsRefCap0

private

Definition at line 89 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::gainsRefCap0vsEta

private

Definition at line 97 of file HcalGainsCheck.h.

TH1F* HcalGainsCheck::gainsRefCap1

private

Definition at line 90 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::gainsRefCap1vsEta

private

Definition at line 98 of file HcalGainsCheck.h.

TH1F* HcalGainsCheck::gainsRefCap2

private

Definition at line 91 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::gainsRefCap2vsEta

private

Definition at line 99 of file HcalGainsCheck.h.

TH1F* HcalGainsCheck::gainsRefCap3

private

Definition at line 92 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::gainsRefCap3vsEta

private

Definition at line 100 of file HcalGainsCheck.h.

TH1F* HcalGainsCheck::gainsUpCap0

private

Definition at line 85 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::gainsUpCap0vsEta

private

Definition at line 93 of file HcalGainsCheck.h.

TH1F* HcalGainsCheck::gainsUpCap1

private

Definition at line 86 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::gainsUpCap1vsEta

private

Definition at line 94 of file HcalGainsCheck.h.

TH1F* HcalGainsCheck::gainsUpCap2

private

Definition at line 87 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::gainsUpCap2vsEta

private

Definition at line 95 of file HcalGainsCheck.h.

TH1F* HcalGainsCheck::gainsUpCap3

private

Definition at line 88 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::gainsUpCap3vsEta

private

Definition at line 96 of file HcalGainsCheck.h.

TH2F* HcalGainsCheck::ocMapRef

private

Definition at line 73 of file HcalGainsCheck.h.

Referenced by beginJob().

TH2F* HcalGainsCheck::ocMapUp

private

Definition at line 72 of file HcalGainsCheck.h.

Referenced by beginJob().

std::string HcalGainsCheck::outfile

private

Definition at line 65 of file HcalGainsCheck.h.

Referenced by analyze(), and HcalGainsCheck().

TH1F* HcalGainsCheck::ratioUpRefCap0

private

Definition at line 81 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::ratioUpRefCap1

private

Definition at line 82 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::ratioUpRefCap2

private

Definition at line 83 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

TH1F* HcalGainsCheck::ratioUpRefCap3

private

Definition at line 84 of file HcalGainsCheck.h.

Referenced by analyze(), and beginJob().

std::string HcalGainsCheck::rootfile

private

Definition at line 64 of file HcalGainsCheck.h.

Referenced by beginJob(), endJob(), and HcalGainsCheck().

bool HcalGainsCheck::validategainsflag

private

Definition at line 67 of file HcalGainsCheck.h.

Referenced by analyze(), and HcalGainsCheck().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation