#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

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

52 {}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 49 of file HcalGainsCheck.cc.

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

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

 {
   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 f, and rootfile.

 {
   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(), and endJob().

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