HcalGainsCheck Class Reference

#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
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

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
typedef WorkerT< EDAnalyzer >	WorkerType
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::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

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

{}

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();
    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().