FastHFShowerLibrary Class Reference

#include <FastHFShowerLibrary.h>

Public Member Functions
	FastHFShowerLibrary (edm::ParameterSet const &, edm::ConsumesCollector &&)
const std::map< CaloHitID, float > &	getHitsMap ()
void const	initHFShowerLibrary (const edm::EventSetup &)
void	modifyDepth (HcalNumberingFromDDD::HcalID &id)
void	recoHFShowerLibrary (const FSimTrack &myTrack)
void	SetRandom (const RandomEngineAndDistribution *)
	~FastHFShowerLibrary ()

Private Attributes
bool	applyFidCut
const edm::ParameterSet	fast
const HcalDDDSimConstants *	hcalConstants
const edm::ESGetToken< HcalDDDSimConstants, HcalSimNumberingRecord >	hcalDDDSimConstantsESToken_
const edm::ESGetToken< HcalSimulationConstants, HcalSimNumberingRecord >	hcalSimulationConstantsESToken_
std::unique_ptr< HFShowerLibrary >	hfshower
std::map< CaloHitID, float >	hitMap
std::string	name
std::unique_ptr< HcalNumberingFromDDD >	numberingFromDDD
HcalNumberingScheme	numberingScheme

Detailed Description

Definition at line 39 of file FastHFShowerLibrary.h.

Constructor & Destructor Documentation

FastHFShowerLibrary()

FastHFShowerLibrary::FastHFShowerLibrary	(	edm::ParameterSet const &	p,
		edm::ConsumesCollector &&	iC
	)

Definition at line 41 of file FastHFShowerLibrary.cc.

References applyFidCut, edm::ParameterSet::getParameter(), and AlCaHLTBitMon_ParallelJobs::p.

    : fast(p), hcalDDDSimConstantsESToken_(iC.esConsumes()), hcalSimulationConstantsESToken_(iC.esConsumes()) {
  edm::ParameterSet m_HS = p.getParameter<edm::ParameterSet>("HFShowerLibrary");
  applyFidCut = m_HS.getParameter<bool>("ApplyFiducialCut");
}

~FastHFShowerLibrary()

FastHFShowerLibrary::~FastHFShowerLibrary ( )

inline

Definition at line 43 of file FastHFShowerLibrary.h.

{}

Member Function Documentation

getHitsMap()

const std::map<CaloHitID, float>& FastHFShowerLibrary::getHitsMap ( )

inline

Definition at line 49 of file FastHFShowerLibrary.h.

References hitMap.

Referenced by CalorimetryManager::HDShowerSimulation(), and CalorimetryManager::reconstructTrack().

{ return hitMap; };

initHFShowerLibrary()

void const FastHFShowerLibrary::initHFShowerLibrary

(

const edm::EventSetup &

iSetup

)

Definition at line 47 of file FastHFShowerLibrary.cc.

References bJpsiMuMuTrigSettings_cff::decays, fast, edm::EventSetup::getData(), hcalConstants, hcalDDDSimConstantsESToken_, HcalSimulationConstants::hcalsimpar(), hcalSimulationConstantsESToken_, hfshower, initializeOnce, name, numberingFromDDD, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                               {
  edm::LogInfo("FastCalorimetry") << "initHFShowerLibrary::initialization";

  hcalConstants = &iSetup.getData(hcalDDDSimConstantsESToken_);
  const HcalSimulationConstants* hsps = &iSetup.getData(hcalSimulationConstantsESToken_);

  std::string name = "HcalHits";
  numberingFromDDD = std::make_unique<HcalNumberingFromDDD>(hcalConstants);
  hfshower = std::make_unique<HFShowerLibrary>(name, hcalConstants, hsps->hcalsimpar(), fast);

  //only one thread can be allowed to setup the G4 physics table.
  std::call_once(initializeOnce, []() {
    // Geant4 particles
    G4DecayPhysics decays;
    decays.ConstructParticle();
    G4ParticleTable* partTable = G4ParticleTable::GetParticleTable();
    partTable->SetReadiness();
  });
  //G4ParticleTable* partTable = G4ParticleTable::GetParticleTable();
  //hfshower->initRun(partTable, hcalConstants);  // init particle code
}

modifyDepth()

void FastHFShowerLibrary::modifyDepth

(

HcalNumberingFromDDD::HcalID &

id

)

Definition at line 131 of file FastHFShowerLibrary.cc.

References hcalRecHitTable_cff::depth, hcalConstants, HcalForward, hcalRecHitTable_cff::ieta, and HcalDDDSimConstants::maxHFDepth().

Referenced by recoHFShowerLibrary().

                                                                    {
  if (id.subdet == HcalForward) {
    int ieta = (id.zside == 0) ? -id.etaR : id.etaR;
    if (hcalConstants->maxHFDepth(ieta, id.phis) > 2) {
      if (id.depth <= 2) {
        if (G4UniformRand() > 0.5)
          id.depth += 2;
      }
    }
  }
}

recoHFShowerLibrary()

void FastHFShowerLibrary::recoHFShowerLibrary

(

const FSimTrack &

myTrack

)

Definition at line 76 of file FastHFShowerLibrary.cc.

References applyFidCut, hcalRecHitTable_cff::depth, RawParticle::e(), HcalNumberingScheme::getUnitID(), hfshower, hitMap, hfClusterShapes_cfi::hits, mps_fire::i, FSimTrack::id(), modifyDepth(), numberingFromDDD, numberingScheme, convertSQLiteXML::ok, FSimTrack::onVFcal(), HFFibreFiducial::PMTNumber(), hcalRecHitTable_cff::time, createJobs::tmp, CoreSimTrack::type(), RawParticle::Vect(), bphysicsOniaDQM_cfi::vertex, RawParticle::vertex(), and FSimTrack::vfcalEntrance().

Referenced by CalorimetryManager::HDShowerSimulation(), and CalorimetryManager::reconstructTrack().

                                                                      {
#ifdef DebugLog
  edm::LogInfo("FastCalorimetry") << "FastHFShowerLibrary: recoHFShowerLibrary ";
#endif

  if (!myTrack.onVFcal()) {
#ifdef DebugLog
    edm::LogInfo("FastCalorimetry") << "FastHFShowerLibrary: we should not be here ";
#endif
  }

  hitMap.clear();
  double eGen = 1000. * myTrack.vfcalEntrance().e();  // energy in [MeV]
  double delZv = (myTrack.vfcalEntrance().vertex().Z() > 0.0) ? 50.0 : -50.0;
  G4ThreeVector vertex(10. * myTrack.vfcalEntrance().vertex().X(),
                       10. * myTrack.vfcalEntrance().vertex().Y(),
                       10. * myTrack.vfcalEntrance().vertex().Z() + delZv);  // in [mm]

  G4ThreeVector direction(
      myTrack.vfcalEntrance().Vect().X(), myTrack.vfcalEntrance().Vect().Y(), myTrack.vfcalEntrance().Vect().Z());

  bool ok;
  double weight = 1.0;  // rad. damage
  int parCode = myTrack.type();
  double tSlice = 0.1 * vertex.mag() / 29.98;

  std::vector<HFShowerLibrary::Hit> hits =
      hfshower->fillHits(vertex, direction, parCode, eGen, ok, weight, tSlice, false);

  for (unsigned int i = 0; i < hits.size(); ++i) {
    G4ThreeVector pos = hits[i].position;
    int depth = hits[i].depth;
    double time = hits[i].time;
    if (!applyFidCut || (HFFibreFiducial::PMTNumber(pos) > 0)) {
      //    if (!applyFidCut || (applyFidCut && HFFibreFiducial::PMTNumber(pos)>0)) {
      int det = 5;
      int lay = 1;
      uint32_t id = 0;
      HcalNumberingFromDDD::HcalID tmp =
          numberingFromDDD->unitID(det, math::XYZVectorD(pos.x(), pos.y(), pos.z()), depth, lay);
      modifyDepth(tmp);
      id = numberingScheme.getUnitID(tmp);

      CaloHitID current_id(id, time, myTrack.id());
      std::map<CaloHitID, float>::iterator cellitr;
      cellitr = hitMap.find(current_id);
      if (cellitr == hitMap.end()) {
        hitMap.insert(std::pair<CaloHitID, float>(current_id, 1.0));
      } else {
        cellitr->second += 1.0;
      }
    }  // end of isItinFidVolume check
  }    // end loop over hits
}

SetRandom()

void FastHFShowerLibrary::SetRandom

(

const RandomEngineAndDistribution *

rnd

)

Definition at line 69 of file FastHFShowerLibrary.cc.

References LogDebug, and RandomEngineAndDistribution::theEngine().

Referenced by CalorimetryManager::HDShowerSimulation(), and CalorimetryManager::initialize().

                                                                          {
  // define Geant4 engine per thread
  G4Random::setTheEngine(&(rnd->theEngine()));
  LogDebug("FastHFShowerLibrary::recoHFShowerLibrary")
      << "Begin of event " << G4UniformRand() << "  " << rnd->theEngine().name() << "  " << rnd->theEngine();
}

Member Data Documentation

applyFidCut

bool FastHFShowerLibrary::applyFidCut

private

Definition at line 62 of file FastHFShowerLibrary.h.

Referenced by FastHFShowerLibrary(), and recoHFShowerLibrary().

fast

const edm::ParameterSet FastHFShowerLibrary::fast

private

Definition at line 54 of file FastHFShowerLibrary.h.

Referenced by initHFShowerLibrary().

hcalConstants

const HcalDDDSimConstants* FastHFShowerLibrary::hcalConstants

private

Definition at line 57 of file FastHFShowerLibrary.h.

Referenced by initHFShowerLibrary(), and modifyDepth().

hcalDDDSimConstantsESToken_

const edm::ESGetToken<HcalDDDSimConstants, HcalSimNumberingRecord> FastHFShowerLibrary::hcalDDDSimConstantsESToken_

private

Definition at line 65 of file FastHFShowerLibrary.h.

Referenced by initHFShowerLibrary().

hcalSimulationConstantsESToken_

const edm::ESGetToken<HcalSimulationConstants, HcalSimNumberingRecord> FastHFShowerLibrary::hcalSimulationConstantsESToken_

private

Definition at line 66 of file FastHFShowerLibrary.h.

Referenced by initHFShowerLibrary().

hfshower

std::unique_ptr<HFShowerLibrary> FastHFShowerLibrary::hfshower

private

Definition at line 55 of file FastHFShowerLibrary.h.

Referenced by initHFShowerLibrary(), and recoHFShowerLibrary().

hitMap

std::map<CaloHitID, float> FastHFShowerLibrary::hitMap

private

Definition at line 60 of file FastHFShowerLibrary.h.

Referenced by getHitsMap(), and recoHFShowerLibrary().

name

std::string FastHFShowerLibrary::name

private

Definition at line 63 of file FastHFShowerLibrary.h.

Referenced by ElectronMVAID.ElectronMVAID::__call__(), FWLite.ElectronMVAID::__call__(), dirstructure.Directory::__create_pie_image(), DisplayManager.DisplayManager::__del__(), dqm_interfaces.DirID::__eq__(), dirstructure.Directory::__get_full_path(), dirstructure.Comparison::__get_img_name(), dataset.Dataset::__getDataType(), dataset.Dataset::__getFileInfoList(), dirstructure.Comparison::__make_image(), core.autovars.NTupleVariable::__repr__(), core.autovars.NTupleObjectType::__repr__(), core.autovars.NTupleObject::__repr__(), core.autovars.NTupleCollection::__repr__(), dirstructure.Directory::__repr__(), dqm_interfaces.DirID::__repr__(), dirstructure.Comparison::__repr__(), config.Service::__setattr__(), config.CFG::__str__(), counter.Counter::__str__(), average.Average::__str__(), FWLite.WorkingPoints::_reformat_cut_definitions(), core.autovars.NTupleObjectType::addSubObjects(), core.autovars.NTupleObjectType::addVariables(), core.autovars.NTupleObjectType::allVars(), dataset.CMSDataset::buildListOfFiles(), dataset.LocalDataset::buildListOfFiles(), dataset.CMSDataset::buildListOfFilesDBS(), dirstructure.Directory::calcStats(), validation.Sample::digest(), python.rootplot.utilities.Hist::divide(), python.rootplot.utilities.Hist::divide_wilson(), DisplayManager.DisplayManager::Draw(), TreeCrawler.Package::dump(), core.autovars.NTupleVariable::fillBranch(), core.autovars.NTupleObject::fillBranches(), core.autovars.NTupleCollection::fillBranchesScalar(), core.autovars.NTupleCollection::fillBranchesVector(), core.autovars.NTupleCollection::get_cpp_declaration(), core.autovars.NTupleCollection::get_cpp_wrapper_class(), core.autovars.NTupleCollection::get_py_wrapper_class(), utils.StatisticalTest::get_status(), production_tasks.Task::getname(), dataset.CMSDataset::getPrimaryDatasetEntries(), dataset.PrivateDataset::getPrimaryDatasetEntries(), initHFShowerLibrary(), VIDSelectorBase.VIDSelectorBase::initialize(), core.autovars.NTupleVariable::makeBranch(), core.autovars.NTupleObject::makeBranches(), core.autovars.NTupleCollection::makeBranchesScalar(), core.autovars.NTupleCollection::makeBranchesVector(), dirstructure.Directory::print_report(), dataset.BaseDataset::printInfo(), dataset.Dataset::printInfo(), production_tasks.MonitorJobs::run(), python.rootplot.utilities.Hist::TGraph(), python.rootplot.utilities.Hist::TH1F(), counter.Counter::write(), and average.Average::write().

numberingFromDDD

std::unique_ptr<HcalNumberingFromDDD> FastHFShowerLibrary::numberingFromDDD

private

Definition at line 56 of file FastHFShowerLibrary.h.

Referenced by initHFShowerLibrary(), and recoHFShowerLibrary().

numberingScheme

HcalNumberingScheme FastHFShowerLibrary::numberingScheme

private

Definition at line 58 of file FastHFShowerLibrary.h.

Referenced by recoHFShowerLibrary().