---

HPDNoiseLibraryReader Class Reference

#include <SimCalorimetry/HcalSimAlgos/interface/HPDNoiseLibraryReader.h>

List of all members.

Public Member Functions
std::vector< std::pair < HcalDetId, const float * > >	getBiasedNoisyHcalDetIds ()
std::vector< std::pair < HcalDetId, const float * > >	getNoisyHcalDetIds ()
	HPDNoiseLibraryReader (const edm::ParameterSet &)
	~HPDNoiseLibraryReader ()
Static Public Member Functions
static void	initializeServices ()
Public Attributes
std::string	theHPDName
std::vector< std::string >	theNames
std::vector< float >	theNoiseRate
std::vector< int >	theNoisyPhi
CLHEP::RandFlat *	theRandFlat
CLHEP::RandGaussQ *	theRandGaussQ
HPDNoiseReader *	theReader
HcalTopology	theTopology
Protected Member Functions
void	setRandomEngine (CLHEP::HepRandomEngine &engine)
void	setRandomEngine ()
Private Member Functions
bool	applyNoise (int iphi)
void	clearPhi ()
void	fillRate ()
void	getBiasedNoisyPhis ()
HPDNoiseData *	getNoiseData (int iphi)
void	getNoisyPhis ()
std::string	itos (int i)

---

Detailed Description

Definition at line 42 of file HPDNoiseLibraryReader.h.

---

Constructor & Destructor Documentation

HPDNoiseLibraryReader::HPDNoiseLibraryReader

(

const edm::ParameterSet &

iConfig

)

Definition at line 22 of file HPDNoiseLibraryReader.cc.

References HPDNoiseReader::allNames(), fillRate(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), setRandomEngine(), theHPDName, theNames, and theReader.

: theNoiseRate(0), theNoisyPhi(0), theRandFlat(0), theRandGaussQ(0)
{
    setRandomEngine();
    
    ParameterSet pSet     = iConfig.getParameter<edm::ParameterSet>("HPDNoiseLibrary");
    FileInPath   filepath = pSet.getParameter<edm::FileInPath>("FileName");
    theHPDName            = pSet.getUntrackedParameter<string> ("HPDName", "HPD");
    string pName          = filepath.fullPath();
    if (pName.find(".") == 0) pName.erase(0,2);
    theReader = new HPDNoiseReader(pName);
    theNames = theReader->allNames();   // all 72x2 HPDs

    fillRate();
}

HPDNoiseLibraryReader::~HPDNoiseLibraryReader

(

)

Definition at line 38 of file HPDNoiseLibraryReader.cc.

References theRandFlat, and theRandGaussQ.

{
    if (theRandFlat)delete theRandFlat;
    if (theRandGaussQ)delete theRandGaussQ;
}

---

Member Function Documentation

bool HPDNoiseLibraryReader::applyNoise

(

int

iphi

)

[private]

Definition at line 205 of file HPDNoiseLibraryReader.cc.

References find(), and theNoisyPhi.

Referenced by getNoiseData().

{

    bool isAccepted = false;
    vector < int >::iterator phi_iter;
    phi_iter = find(theNoisyPhi.begin(), theNoisyPhi.end(), iphi);
    if (phi_iter != theNoisyPhi.end()){
        isAccepted = true;
    }
    return isAccepted;
}

void HPDNoiseLibraryReader::clearPhi

(

)

[private]

Definition at line 223 of file HPDNoiseLibraryReader.cc.

References theNoisyPhi.

Referenced by getBiasedNoisyPhis(), and getNoisyPhis().

                                     {
    theNoisyPhi.clear();
}

void HPDNoiseLibraryReader::fillRate

(

)

[private]

Definition at line 82 of file HPDNoiseLibraryReader.cc.

References TestMuL1L2Filter_cff::cerr, lat::endl(), HPDNoiseReader::getHandle(), i, HPDNoiseReader::rate(), theNames, theNoiseRate, theReader, and HPDNoiseReader::valid().

Referenced by HPDNoiseLibraryReader().

{
    for (size_t i = 0; i < theNames.size(); ++i) {
        HPDNoiseReader::Handle hpdObj = theReader->getHandle(theNames[i]);
        if (theReader->valid(hpdObj)) {
            theNoiseRate.push_back(theReader->rate(hpdObj));
        } else {
            std::cerr << "HPD Handle Object is not valid!" << endl;
        }
    }
}

vector< pair< HcalDetId, const float * > > HPDNoiseLibraryReader::getBiasedNoisyHcalDetIds

(

)

Definition at line 184 of file HPDNoiseLibraryReader.cc.

References data, getBiasedNoisyPhis(), HPDNoiseData::getDataFrame(), HPDNoiseDataFrame::getFrame(), getNoiseData(), i, HPDNoiseDataFrame::id(), int, HLT_VtxMuL3::result, HPDNoiseData::size(), and theNoisyPhi.

{
    
    vector< pair<HcalDetId, const float* > > result;
    // decide which phi are noisy by using noise rates 
    // and random numbers (to be called for each event)
    // at least one Phi is always noisy.
    getBiasedNoisyPhis();
    for (int i = 0; i < int(theNoisyPhi.size()); ++i) {
        int iphi = theNoisyPhi[i];
        HPDNoiseData* data;
        data = getNoiseData(iphi);
        for(unsigned int i=0; i<data->size();++i)
        {
            pair < HcalDetId, const float* >tmp_pair( data->getDataFrame(i).id(), data->getDataFrame(i).getFrame());
            result.push_back(tmp_pair);
        }
    }
    return result;
}

void HPDNoiseLibraryReader::getBiasedNoisyPhis

(

)

[private]

Definition at line 144 of file HPDNoiseLibraryReader.cc.

References clearPhi(), i, theNoiseRate, theNoisyPhi, and theRandFlat.

Referenced by getBiasedNoisyHcalDetIds().

{

    clearPhi();
    double rndm[144];

    theRandFlat->shootArray(144, rndm);
    for (int i = 0; i < 144; ++i) {
        if (rndm[i] < theNoiseRate[i]) {
            theNoisyPhi.push_back(i + 1);
        }
    }
    // make sure one HPD is always noisy
    if (theNoisyPhi.size() == 0) {
        int iPhi = (theRandFlat->fireInt(144)) + 1; // integer from interval [0-144[ + 1

        theNoisyPhi.push_back(iPhi);
    }
}

HPDNoiseData * HPDNoiseLibraryReader::getNoiseData

(

int

iphi

)

[private]

Definition at line 94 of file HPDNoiseLibraryReader.cc.

References applyNoise(), TestMuL1L2Filter_cff::cerr, data, lat::endl(), HPDNoiseReader::getEntry(), HPDNoiseReader::getHandle(), itos(), name, theHPDName, theRandFlat, theReader, HPDNoiseReader::totalEntries(), and HPDNoiseReader::valid().

Referenced by getBiasedNoisyHcalDetIds(), and getNoisyHcalDetIds().

{
    
    
    HPDNoiseData* data;
    // make sure that iphi from HcalDetId is found noisy at first.
    // In other words, be sure that iphi is in the collection of
    // noisy Phis
    if (!(applyNoise(iphi))) return data;
    
    int zside = 1;
    if(iphi>72){
       iphi = iphi-72;
       zside = -1;
    }
    std::string name;
    if (zside == 1) {
        name = "ZPlus" + theHPDName + itos(iphi);
    } else if (zside == -1){
        name = "ZMinus" + theHPDName + itos(iphi);
    }else {
        cerr << " ZSide Calculation Error." << endl;
    }
    HPDNoiseReader::Handle hpdObj = theReader->getHandle(name);
    if (theReader->valid(hpdObj)) {
        // randomly select one entry from library for this HPD
        unsigned long entry = theRandFlat->fireInt( theReader->totalEntries(hpdObj));
        theReader->getEntry(hpdObj, entry, &data);
    }else{
        std::cerr << " HPD Name in the library is not valid." << std::endl;
    }
    return data;
}

vector< pair< HcalDetId, const float * > > HPDNoiseLibraryReader::getNoisyHcalDetIds

(

)

Definition at line 164 of file HPDNoiseLibraryReader.cc.

References data, HPDNoiseData::getDataFrame(), HPDNoiseDataFrame::getFrame(), getNoiseData(), getNoisyPhis(), i, HPDNoiseDataFrame::id(), int, HLT_VtxMuL3::result, HPDNoiseData::size(), and theNoisyPhi.

{
    
    vector< pair<HcalDetId, const float* > > result;
    // decide which phi are noisy by using noise rates 
    // and random numbers (to be called for each event)
    getNoisyPhis();
    for (int i = 0; i < int(theNoisyPhi.size()); ++i) {
        int iphi = theNoisyPhi[i];
        HPDNoiseData* data;
        data = getNoiseData(iphi);
        for(unsigned int i=0; i<data->size();++i)
        {
            pair < HcalDetId, const float* >tmp_pair( data->getDataFrame(i).id(), data->getDataFrame(i).getFrame());
            result.push_back(tmp_pair);
        }
    }
    return result;
}

void HPDNoiseLibraryReader::getNoisyPhis

(

)

[private]

Definition at line 129 of file HPDNoiseLibraryReader.cc.

References clearPhi(), i, theNoiseRate, theNoisyPhi, and theRandFlat.

Referenced by getNoisyHcalDetIds().

{

    clearPhi();
    double rndm[144];

    theRandFlat->shootArray(144, rndm);

    for (int i = 0; i < 144; ++i) {
        if (rndm[i] < theNoiseRate[i]) {
            theNoisyPhi.push_back(i + 1);
        }
    }
}

void HPDNoiseLibraryReader::initializeServices

(

)

[static]

Definition at line 44 of file HPDNoiseLibraryReader.cc.

References edmplugin::standard::config(), edmplugin::PluginManager::configure(), edm::ServiceRegistry::createServicesFromConfig(), and edmplugin::PluginManager::isAvailable().

{
    if (not edmplugin::PluginManager::isAvailable()) {
        edmplugin::PluginManager::configure(edmplugin::standard::config());
    }

    std::string config =
      "process CorrNoise = {"
          "service = RandomNumberGeneratorService" 
             "{" 
                  "untracked uint32 sourceSeed = 123456789" 
              "}" 
      "}";

    // create the services
    edm::ServiceToken tempToken = edm::ServiceRegistry::createServicesFromConfig(config);

    // make the services available
    edm::ServiceRegistry::Operate operate(tempToken);
}

string HPDNoiseLibraryReader::itos

(

int

i

)

[private]

Definition at line 216 of file HPDNoiseLibraryReader.cc.

References s.

Referenced by getNoiseData().

{
    stringstream s;
    s << i;
    return s.str();
}

void HPDNoiseLibraryReader::setRandomEngine

(

CLHEP::HepRandomEngine &

engine

)

[protected]

Definition at line 74 of file HPDNoiseLibraryReader.cc.

References theRandFlat, and theRandGaussQ.

{
    if(theRandGaussQ) delete theRandGaussQ;
    if(theRandFlat) delete theRandFlat;
    theRandGaussQ = new CLHEP::RandGaussQ(engine);
    theRandFlat = new CLHEP::RandFlat(engine);
}

void HPDNoiseLibraryReader::setRandomEngine

(

)

[protected]

Definition at line 64 of file HPDNoiseLibraryReader.cc.

References Exception, and edm::Service< T >::isAvailable().

Referenced by HPDNoiseLibraryReader().

{
    edm::Service < edm::RandomNumberGenerator > rng;
    if (!rng.isAvailable()) {
        throw cms::Exception("Configuration") << "HcalHPDNoiseLibrary requires the RandomNumberGeneratorService\n"
          "which is not present in the configuration file.  You must add the service\n"
          "in the configuration file or remove the modules that require it.";
    }
    setRandomEngine(rng->getEngine());
}

---

Member Data Documentation

std::string HPDNoiseLibraryReader::theHPDName

Definition at line 85 of file HPDNoiseLibraryReader.h.

Referenced by getNoiseData(), and HPDNoiseLibraryReader().

std::vector<std::string> HPDNoiseLibraryReader::theNames

Definition at line 84 of file HPDNoiseLibraryReader.h.

Referenced by fillRate(), and HPDNoiseLibraryReader().

std::vector<float> HPDNoiseLibraryReader::theNoiseRate

Definition at line 79 of file HPDNoiseLibraryReader.h.

Referenced by fillRate(), getBiasedNoisyPhis(), and getNoisyPhis().

std::vector<int> HPDNoiseLibraryReader::theNoisyPhi

Definition at line 80 of file HPDNoiseLibraryReader.h.

Referenced by applyNoise(), clearPhi(), getBiasedNoisyHcalDetIds(), getBiasedNoisyPhis(), getNoisyHcalDetIds(), and getNoisyPhis().

CLHEP::RandFlat* HPDNoiseLibraryReader::theRandFlat

Definition at line 81 of file HPDNoiseLibraryReader.h.

Referenced by getBiasedNoisyPhis(), getNoiseData(), getNoisyPhis(), setRandomEngine(), and ~HPDNoiseLibraryReader().

CLHEP::RandGaussQ* HPDNoiseLibraryReader::theRandGaussQ

Definition at line 82 of file HPDNoiseLibraryReader.h.

Referenced by setRandomEngine(), and ~HPDNoiseLibraryReader().

HPDNoiseReader* HPDNoiseLibraryReader::theReader

Definition at line 83 of file HPDNoiseLibraryReader.h.

Referenced by fillRate(), getNoiseData(), and HPDNoiseLibraryReader().

HcalTopology HPDNoiseLibraryReader::theTopology

Definition at line 76 of file HPDNoiseLibraryReader.h.

---

The documentation for this class was generated from the following files:

• SimCalorimetry/HcalSimAlgos/interface/HPDNoiseLibraryReader.h
• SimCalorimetry/HcalSimAlgos/src/HPDNoiseLibraryReader.cc

---