|
|
|
#include <HcalGainsCheck.h>
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 () | |
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 |
Definition at line 47 of file HcalGainsCheck.h.
| 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.
{}
| 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, CastorDbASCIIIO::dumpObject(), dumprefs, dumpupdate, emapflag, epsilon, 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.
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(), and HcalGainsCheck().
TFile* HcalGainsCheck::f [private] |
Definition at line 70 of file HcalGainsCheck.h.
Referenced by beginJob(), 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().
1.7.3