SimG4HcalValidation Class Reference

#include <SimG4HcalValidation.h>

Inheritance diagram for SimG4HcalValidation:

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override
	SimG4HcalValidation (const edm::ParameterSet &p)
	~SimG4HcalValidation () override
Public Member Functions inherited from SimProducer
void	registerProducts (edm::ProducesCollector producesCollector)
	SimProducer ()
Public Member Functions inherited from SimWatcher
	SimWatcher ()
virtual	~SimWatcher ()
Public Member Functions inherited from Observer< const BeginOfJob * >
	Observer ()
void	slotForUpdate (const BeginOfJob * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfRun * >
	Observer ()
void	slotForUpdate (const BeginOfRun * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfEvent * >
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const EndOfEvent * >
	Observer ()
void	slotForUpdate (const EndOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const G4Step * >
	Observer ()
void	slotForUpdate (const G4Step * iT)
virtual	~Observer ()

Private Member Functions
void	clear ()
void	collectEnergyRdir (const double, const double)
void	fetchHits (PHcalValidInfoLayer &)
void	fill (const EndOfEvent *ev)
double	getHcalScale (std::string, int) const
void	init ()
void	jetAnalysis (PHcalValidInfoJets &)
void	layerAnalysis (PHcalValidInfoLayer &)
void	nxNAnalysis (PHcalValidInfoNxN &)
const SimG4HcalValidation &	operator= (const SimG4HcalValidation &)=delete
	SimG4HcalValidation (const SimG4HcalValidation &)=delete
void	update (const BeginOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfJob *job) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const G4Step *step) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
bool	applySampling
double	coneSize
unsigned int	count
std::vector< double >	dEta
std::vector< double >	dPhi
double	edepd [5]
double	edepEB
double	edepEE
double	edepHB
double	edepHE
double	edepHO
double	edepl [20]
double	een
double	ehitThreshold
double	enEcal
double	enHcal
float	eta0
bool	hcalOnly
double	hen
double	hhitThreshold
std::vector< CaloHit >	hitcache
double	hoen
int	infolevel
SimG4HcalHitJetFinder *	jetf
float	jetThreshold
std::string	labelJets
std::string	labelLayer
std::string	labelNxN
std::vector< std::string >	names
HcalNumberingFromDDD *	numberingFromDDD
HcalTestNumberingScheme *	org
float	phi0
std::vector< float >	scaleHB
std::vector< float >	scaleHE
std::vector< float >	scaleHF
float	timeLowlim
float	timeUplim
double	vhitec
double	vhithc

Additional Inherited Members
Protected Member Functions inherited from SimProducer
template<class T >
void	produces ()
template<class T >
void	produces (const std::string &instanceName)

Detailed Description

Definition at line 34 of file SimG4HcalValidation.h.

Constructor & Destructor Documentation

SimG4HcalValidation() [1/2]

SimG4HcalValidation::SimG4HcalValidation

(

const edm::ParameterSet &

p

)

Definition at line 36 of file SimG4HcalValidation.cc.

    : jetf(nullptr), numberingFromDDD(nullptr), org(nullptr) {
  edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("SimG4HcalValidation");
  infolevel = m_Anal.getParameter<int>("InfoLevel");
  hcalOnly = m_Anal.getParameter<bool>("HcalClusterOnly");
  applySampling = m_Anal.getParameter<bool>("HcalSampling");
  coneSize = m_Anal.getParameter<double>("ConeSize");
  ehitThreshold = m_Anal.getParameter<double>("EcalHitThreshold");
  hhitThreshold = m_Anal.getParameter<double>("HcalHitThreshold");
  timeLowlim = m_Anal.getParameter<double>("TimeLowLimit");
  timeUplim = m_Anal.getParameter<double>("TimeUpLimit");
  jetThreshold = m_Anal.getParameter<double>("JetThreshold");
  eta0 = m_Anal.getParameter<double>("Eta0");
  phi0 = m_Anal.getParameter<double>("Phi0");
  names = m_Anal.getParameter<std::vector<std::string>>("Names");
  labelLayer = m_Anal.getUntrackedParameter<std::string>("LabelLayerInfo", "HcalInfoLayer");
  labelNxN = m_Anal.getUntrackedParameter<std::string>("LabelNxNInfo", "HcalInfoNxN");
  labelJets = m_Anal.getUntrackedParameter<std::string>("LabelJetsInfo", "HcalInfoJets");
 
  produces<PHcalValidInfoLayer>(labelLayer);
  if (infolevel > 0)
    produces<PHcalValidInfoNxN>(labelNxN);
  if (infolevel > 1)
    produces<PHcalValidInfoJets>(labelJets);
 
  edm::LogVerbatim("ValidHcal") << "HcalTestAnalysis:: Initialised as observer of begin/end events and "
                                << "of G4step with Parameter values: \n\tInfoLevel     = " << infolevel
                                << "\n\thcalOnly      = " << hcalOnly << "\n\tapplySampling = " << applySampling
                                << "\n\tconeSize      = " << coneSize << "\n\tehitThreshold = " << ehitThreshold
                                << "\n\thhitThreshold = " << hhitThreshold << "\n\tttimeLowlim   = " << timeLowlim
                                << "\n\tttimeUplim    = " << timeUplim << "\n\teta0          = " << eta0
                                << "\n\tphi0          = " << phi0 << "\nLabels (Layer): " << labelLayer
                                << " (NxN): " << labelNxN << " (Jets): " << labelJets;
 
  init();
}

References applySampling, coneSize, ehitThreshold, eta0, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalOnly, hhitThreshold, infolevel, init(), jetThreshold, labelJets, labelLayer, labelNxN, names, AlCaHLTBitMon_ParallelJobs::p, phi0, AlCaHLTBitMon_QueryRunRegistry::string, timeLowlim, and timeUplim.

~SimG4HcalValidation()

SimG4HcalValidation::~SimG4HcalValidation ( )

override

Definition at line 73 of file SimG4HcalValidation.cc.

                                          {
  edm::LogVerbatim("ValidHcal") << "\n -------->  Total number of selected entries"
                                << " : " << count << "\nPointers:: JettFinder " << jetf << ", Numbering Scheme " << org
                                << " and FromDDD " << numberingFromDDD;
  if (jetf) {
    edm::LogVerbatim("ValidHcal") << "Delete Jetfinder";
    delete jetf;
    jetf = nullptr;
  }
  if (numberingFromDDD) {
    edm::LogVerbatim("ValidHcal") << "Delete HcalNumberingFromDDD";
    delete numberingFromDDD;
    numberingFromDDD = nullptr;
  }
}

References count, jetf, numberingFromDDD, and org.

SimG4HcalValidation() [2/2]

SimG4HcalValidation::SimG4HcalValidation ( const SimG4HcalValidation &  )

privatedelete

Member Function Documentation

clear()

void SimG4HcalValidation::clear ( void  )

private

Definition at line 677 of file SimG4HcalValidation.cc.

{ hitcache.erase(hitcache.begin(), hitcache.end()); }

References hitcache.

Referenced by update().

collectEnergyRdir()

void SimG4HcalValidation::collectEnergyRdir ( const double  eta0, const double  phi0  )

private

Definition at line 680 of file SimG4HcalValidation.cc.

                                                                                {
  std::vector<CaloHit> *hits = &hitcache;
  std::vector<CaloHit>::iterator hit_itr;
 
  double sume = 0., sumh = 0., sumho = 0.;
 
  for (hit_itr = hits->begin(); hit_itr < hits->end(); hit_itr++) {
    double e = hit_itr->e();
    double eta = hit_itr->eta();
    double phi = hit_itr->phi();
 
    int type = hit_itr->det();
 
    double r = jetf->rDist(eta0, phi0, eta, phi);
    if (r < coneSize) {
      if (type == 10 || type == 11 || type == 12) {
        sume += e;
      } else {
        sumh += e;
        if (type == static_cast<int>(HcalBarrel) && hit_itr->layer() > 17)
          sumho += e;
      }
    }
  }
 
  een = sume;
  hen = sumh;
  hoen = sumho;
}

References coneSize, MillePedeFileConverter_cfg::e, een, PVValHelper::eta, eta0, HcalBarrel, hen, hitcache, hfClusterShapes_cfi::hits, hoen, jetf, phi, phi0, alignCSCRings::r, and SimG4HcalHitJetFinder::rDist().

Referenced by jetAnalysis(), and nxNAnalysis().

fetchHits()

void SimG4HcalValidation::fetchHits ( PHcalValidInfoLayer &  product )

private

Definition at line 610 of file SimG4HcalValidation.cc.

                                                                {
  LogDebug("ValidHcal") << "Enter SimG4HcalValidation::fetchHits with " << hitcache.size() << " hits";
  int nHit = hitcache.size();
  int hit = 0;
  int i;
  std::vector<CaloHit>::iterator itr;
  std::vector<CaloHit *> lhits(nHit);
  for (i = 0, itr = hitcache.begin(); itr != hitcache.end(); i++, itr++) {
    uint32_t unitID = itr->id();
    int subdet, zside, group, ieta, iphi, lay;
    HcalTestNumbering::unpackHcalIndex(unitID, subdet, zside, group, ieta, iphi, lay);
    subdet = itr->det();
    lay = itr->layer();
    group = (subdet & 15) << 20;
    group += ((lay - 1) & 31) << 15;
    group += (zside & 1) << 14;
    group += (ieta & 127) << 7;
    group += (iphi & 127);
    itr->setId(group);
    lhits[i] = &hitcache[i];
    LogDebug("ValidHcal") << "SimG4HcalValidation::fetchHits:Original " << i << " " << hitcache[i] << "\n"
                          << "SimG4HcalValidation::fetchHits:Copied   " << i << " " << *lhits[i];
  }
  sort(lhits.begin(), lhits.end(), CaloHitIdMore());
  std::vector<CaloHit *>::iterator k1, k2;
  for (i = 0, k1 = lhits.begin(); k1 != lhits.end(); i++, k1++)
    LogDebug("ValidHcal") << "SimG4HcalValidation::fetchHits:Sorted " << i << " " << **k1;
  int nHits = 0;
  for (i = 0, k1 = lhits.begin(); k1 != lhits.end(); i++, k1++) {
    double ehit = (**k1).e();
    double t = (**k1).t();
    uint32_t unitID = (**k1).id();
    int jump = 0;
    LogDebug("ValidHcal") << "SimG4HcalValidation::fetchHits:Start " << i << " U/T/E"
                          << " 0x" << std::hex << unitID << std::dec << " " << t << " " << ehit;
    for (k2 = k1 + 1; k2 != lhits.end() && (t - (**k2).t()) < 1 && (t - (**k2).t()) > -1 && unitID == (**k2).id();
         k2++) {
      ehit += (**k2).e();
      LogDebug("ValidHcal") << "\t + " << (**k2).e();
      jump++;
    }
    LogDebug("ValidHcal") << "\t = " << ehit << " in " << jump;
 
    double eta = (*k1)->eta();
    double phi = (*k1)->phi();
    int lay = ((unitID >> 15) & 31) + 1;
    int subdet = (unitID >> 20) & 15;
    int zside = (unitID >> 14) & 1;
    int ieta = (unitID >> 7) & 127;
    int iphi = (unitID)&127;
 
    // All hits in cache
    product.fillHits(nHits, lay, subdet, eta, phi, ehit, t);
    nHits++;
 
    LogDebug("ValidHcal") << "SimG4HcalValidation::fetchHits:Hit " << nHits << " " << i << " ID 0x" << std::hex
                          << unitID << "   det " << std::dec << subdet << " " << lay << " " << zside << " " << ieta
                          << " " << iphi << " Time " << t << " E " << ehit;
 
    i += jump;
    k1 += jump;
  }
 
  LogDebug("ValidHcal") << "SimG4HcalValidation::fetchHits called with " << nHit << " hits"
                        << " and writes out " << nHits << '(' << hit << ") hits";
}

References TauDecayModes::dec, PVValHelper::eta, PHcalValidInfoLayer::fillHits(), watchdog::group, hitcache, mps_fire::i, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, LogDebug, phi, OrderedSet::t, HcalTestNumbering::unpackHcalIndex(), and ecaldqm::zside().

Referenced by layerAnalysis().

fill()

void SimG4HcalValidation::fill ( const EndOfEvent *  ev )

private

Definition at line 265 of file SimG4HcalValidation.cc.

                                                    {
  LogDebug("ValidHcal") << "SimG4HcalValidation:Fill event " << (*evt)()->GetEventID();
 
  // access to the G4 hit collections
  G4HCofThisEvent *allHC = (*evt)()->GetHCofThisEvent();
 
  int nhc = 0, j = 0;
 
  // Hcal
  int HCHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names[0]);
  CaloG4HitCollection *theHCHC = (CaloG4HitCollection *)allHC->GetHC(HCHCid);
  LogDebug("ValidHcal") << "SimG4HcalValidation :: Hit Collection for " << names[0] << " of ID " << HCHCid
                        << " is obtained at " << theHCHC;
  if (HCHCid >= 0 && theHCHC != nullptr) {
    int nhits = theHCHC->entries();
    for (j = 0; j < nhits; j++) {
      CaloG4Hit *aHit = (*theHCHC)[j];
 
      double e = aHit->getEnergyDeposit() / GeV;
      double time = aHit->getTimeSlice();
 
      math::XYZPoint pos = aHit->getPosition();
      double theta = pos.theta();
      double eta = -log(tan(theta * 0.5));
      double phi = pos.phi();
 
      uint32_t unitID = aHit->getUnitID();
      int subdet, zside, layer, etaIndex, phiIndex, lay;
      org->unpackHcalIndex(unitID, subdet, zside, layer, etaIndex, phiIndex, lay);
 
      // some logic to separate HO ...
      layer--;
      std::string det = "HB";
      if (subdet == static_cast<int>(HcalForward)) {
        det = "HF";
        uint16_t depth = aHit->getDepth();
        if (depth != 0) {
          layer += 2;
          lay += 2;
        }
        if (layer == 1)
          vhithc += e;
        else if (layer == 0)
          vhitec += e;
        else
          edm::LogVerbatim("ValidHcal") << "SimG4HcalValidation::HitPMT " << subdet << " " << (2 * zside - 1) * etaIndex
                                        << " " << phiIndex << " " << layer + 1 << " R " << pos.rho() << " Phi "
                                        << convertRadToDeg(phi) << " Edep " << e << " Time " << time;
      } else if (subdet == static_cast<int>(HcalEndcap)) {
        if (etaIndex <= 20) {
          det = "HES";
        } else {
          det = "HED";
        }
      }
      LogDebug("ValidHcal") << "SimG4HcalValidation::debugFill Hcal " << det << " layer " << std::setw(2) << layer
                            << " lay " << std::setw(2) << lay << " time " << std::setw(6) << time << " theta "
                            << std::setw(8) << theta << " eta " << std::setw(8) << eta << " phi " << std::setw(8) << phi
                            << " e " << std::setw(8) << e << " ID 0x" << std::hex << unitID << " ID dec " << std::dec
                            << (int)unitID;
 
      // if desired, apply sampling factors in HCAL !!!
      if (applySampling)
        e *= getHcalScale(det, layer);
 
      //    filter on time & energy
      if (time >= timeLowlim && time <= timeUplim && e > hhitThreshold) {
        enHcal += e;
        CaloHit ahit(subdet, lay, e, eta, phi, time, unitID);
        hitcache.push_back(ahit);  // fill cache
        ++nhc;
      }
    }
  }
  LogDebug("ValidHcal") << "SimG4HcalValidation:: HCAL hits : " << nhc;
 
  if (!hcalOnly) {  //--------------------------  ECAL hits --------------------
    int ndets = names.size();
    if (ndets > 3)
      ndets = 3;
    for (int idty = 1; idty < ndets; idty++) {
      std::string det = "EB";
      if (idty == 2)
        det = "EE";
      else if (idty > 2)
        det = "ES";
 
      int nec = 0;
      int ECHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names[idty]);
      CaloG4HitCollection *theECHC = (CaloG4HitCollection *)allHC->GetHC(ECHCid);
      LogDebug("ValidHcal") << "SimG4HcalValidation:: Hit Collection for " << names[idty] << " of ID " << ECHCid
                            << " is obtained at " << theECHC;
      int theechc_entries = theECHC->entries();
      if (ECHCid >= 0 && theECHC != nullptr) {
        for (j = 0; j < theechc_entries; j++) {
          CaloG4Hit *aHit = (*theECHC)[j];
 
          double e = aHit->getEnergyDeposit() / GeV;
          double time = aHit->getTimeSlice();
 
          math::XYZPoint pos = aHit->getPosition();
          double theta = pos.theta();
          double eta = -log(tan(theta / 2.));
          double phi = pos.phi();
          HcalNumberingFromDDD::HcalID id = numberingFromDDD->unitID(eta, phi, 1, 1);
          uint32_t unitID = org->getUnitID(id);
          int subdet, zside, layer, ieta, iphi, lay;
          org->unpackHcalIndex(unitID, subdet, zside, layer, ieta, iphi, lay);
          subdet = idty + 9;
          layer = 0;
          unitID = org->packHcalIndex(subdet, zside, layer, ieta, iphi, lay);
 
          //    filter on time & energy
          if (time >= timeLowlim && time <= timeUplim && e > ehitThreshold) {
            enEcal += e;
            CaloHit ahit(subdet, lay, e, eta, phi, time, unitID);
            hitcache.push_back(ahit);  // fill cache
            ++nec;
          }
 
          LogDebug("ValidHcal") << "SimG4HcalValidation::debugFill Ecal " << det << " layer " << std::setw(2) << layer
                                << " lay " << std::setw(2) << lay << " time " << std::setw(6) << time << " theta "
                                << std::setw(8) << theta << " eta " << std::setw(8) << eta << " phi " << std::setw(8)
                                << phi << " e " << std::setw(8) << e << " ID 0x" << std::hex << unitID << " ID dec "
                                << std::dec << (int)unitID;
        }
      }
 
      LogDebug("ValidHcal") << "SimG4HcalValidation:: " << det << " hits : " << nec;
    }
  }  // end of if(!hcalOnly)
}

References applySampling, angle_units::operators::convertRadToDeg(), TauDecayModes::dec, LEDCalibrationChannels::depth, MillePedeFileConverter_cfg::e, ehitThreshold, enEcal, enHcal, PVValHelper::eta, CaloG4Hit::getDepth(), CaloG4Hit::getEnergyDeposit(), getHcalScale(), CaloG4Hit::getPosition(), CaloG4Hit::getTimeSlice(), HcalTestNumberingScheme::getUnitID(), CaloG4Hit::getUnitID(), GeV, HcalEndcap, HcalForward, hcalOnly, hhitThreshold, hitcache, LEDCalibrationChannels::ieta, createfilelist::int, LEDCalibrationChannels::iphi, dqmiolumiharvest::j, dqm-mbProfile::log, LogDebug, names, nhits, numberingFromDDD, org, HcalTestNumberingScheme::packHcalIndex(), phi, AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), ntuplemaker::time, timeLowlim, HcalNumberingFromDDD::unitID(), HcalTestNumberingScheme::unpackHcalIndex(), vhitec, vhithc, and ecaldqm::zside().

Referenced by update().

getHcalScale()

double SimG4HcalValidation::getHcalScale ( std::string  det, int  layer  ) const

private

Definition at line 711 of file SimG4HcalValidation.cc.

                                                                       {
  double tmp = 0.;
 
  if (det == "HB") {
    tmp = scaleHB[layer];
  } else if (det == "HES" || det == "HED") {
    tmp = scaleHE[layer];
  } else if (det == "HF") {
    tmp = scaleHF[layer];
  }
 
  return tmp;
}

References scaleHB, scaleHE, scaleHF, and createJobs::tmp.

Referenced by fill().

init()

void SimG4HcalValidation::init ( void  )

private

Definition at line 108 of file SimG4HcalValidation.cc.

                               {
  float sHB[4] = {117., 117., 178., 217.};
  float sHE[3] = {178., 178., 178.};
  float sHF[3] = {2.84, 2.09, 0.};  //
 
  float deta[4] = {0.0435, 0.1305, 0.2175, 0.3045};
  float dphi[4] = {0.0436, 0.1309, 0.2182, 0.3054};
 
  int i = 0;
  for (i = 0; i < 4; i++) {
    scaleHB.push_back(sHB[i]);
  }
  for (i = 0; i < 3; i++) {
    scaleHE.push_back(sHE[i]);
  }
  for (i = 0; i < 3; i++) {
    scaleHF.push_back(sHF[i]);
  }
 
  // window steps;
  for (i = 0; i < 4; i++) {
    dEta.push_back(deta[i]);
    dPhi.push_back(dphi[i]);
  }
 
  // jetfinder conse size setting
  jetf = new SimG4HcalHitJetFinder(coneSize);
 
  // counter
  count = 0;
}

References coneSize, count, dEta, dPhi, mps_fire::i, jetf, scaleHB, scaleHE, and scaleHF.

Referenced by SimG4HcalValidation().

jetAnalysis()

void SimG4HcalValidation::jetAnalysis ( PHcalValidInfoJets &  product )

private

Definition at line 488 of file SimG4HcalValidation.cc.

                                                                 {
  std::vector<CaloHit> *hhit = &hitcache;
 
  jetf->setInput(hhit);
 
  // zeroing cluster list, perfom clustering, fill cluster list & return pntr
  std::vector<SimG4HcalHitCluster> *result = jetf->getClusters(hcalOnly);
 
  std::vector<SimG4HcalHitCluster>::iterator clus_itr;
 
  LogDebug("ValidHcal") << "\n ---------- Final list of " << (*result).size() << " clusters ---------------";
  for (clus_itr = result->begin(); clus_itr < result->end(); clus_itr++)
    LogDebug("ValidHcal") << (*clus_itr);
 
  std::vector<double> enevec, etavec, phivec;
 
  if (!(*result).empty()) {
    sort((*result).begin(), (*result).end());
 
    clus_itr = result->begin();  // first cluster only
    double etac = clus_itr->eta();
    double phic = clus_itr->phi();
 
    double ecal_collect = 0.;  // collect Ecal energy in the cone
    if (!hcalOnly) {
      ecal_collect = clus_itr->collectEcalEnergyR();
    } else {
      collectEnergyRdir(etac, phic);
      ecal_collect = een;
    }
    LogDebug("ValidHcal") << " JetAnalysis ===> ecal_collect  = " << ecal_collect;
 
    // eta-phi deviation of the cluster from nominal (eta0, phi0) values
    double dist = jetf->rDist(eta0, phi0, etac, phic);
    LogDebug("ValidHcal") << " JetAnalysis ===> eta phi deviation  = " << dist;
    product.fillEtaPhiProfileJet(eta0, phi0, etac, phic, dist);
 
    LogDebug("ValidHcal") << " JetAnalysis ===> after fillEtaPhiProfileJet";
 
    std::vector<CaloHit> *hits = clus_itr->getHits();
    std::vector<CaloHit>::iterator hit_itr;
 
    double ee = 0., he = 0., hoe = 0., etot = 0.;
 
    // cycle over all hits in the FIRST cluster
    for (hit_itr = hits->begin(); hit_itr < hits->end(); hit_itr++) {
      double e = hit_itr->e();
      double t = hit_itr->t();
      double r = jetf->rDist(&(*clus_itr), &(*hit_itr));
 
      // energy collection
      etot += e;
      if (hit_itr->det() == 10 || hit_itr->det() == 11 || hit_itr->det() == 12)
        ee += e;
      if (hit_itr->det() == static_cast<int>(HcalBarrel) || hit_itr->det() == static_cast<int>(HcalEndcap) ||
          hit_itr->det() == static_cast<int>(HcalForward)) {
        he += e;
        if (hit_itr->det() == static_cast<int>(HcalBarrel) && hit_itr->layer() > 17)
          hoe += e;
      }
 
      if (hit_itr->det() == static_cast<int>(HcalBarrel) || hit_itr->det() == static_cast<int>(HcalEndcap) ||
          hit_itr->det() == static_cast<int>(HcalForward)) {
        product.fillTProfileJet(he, r, t);
      }
    }
 
    product.fillEcollectJet(ee, he, hoe, etot);
 
    LogDebug("ValidHcal") << " JetAnalysis ===> after fillEcollectJet: "
                          << "ee/he/hoe/etot " << ee << "/" << he << "/" << hoe << "/" << etot;
 
    // Loop over clusters
    for (clus_itr = result->begin(); clus_itr < result->end(); clus_itr++) {
      if (clus_itr->e() > jetThreshold) {
        enevec.push_back(clus_itr->e());
        etavec.push_back(clus_itr->eta());
        phivec.push_back(clus_itr->phi());
      }
    }
    product.fillJets(enevec, etavec, phivec);
 
    LogDebug("ValidHcal") << " JetAnalysis ===> after fillJets\n"
                          << " JetAnalysis ===> (*result).size() " << (*result).size();
 
    // Di-jet mass
    if (etavec.size() > 1) {
      if (etavec[0] > -2.5 && etavec[0] < 2.5 && etavec[1] > -2.5 && etavec[1] < 2.5) {
        LogDebug("ValidHcal") << " JetAnalysis ===> Di-jet mass enter\n"
                              << " JetAnalysis ===> Di-jet vectors ";
        for (unsigned int i = 0; i < enevec.size(); i++)
          LogDebug("ValidHcal") << "   e, eta, phi = " << enevec[i] << " " << etavec[i] << " " << phivec[i];
 
        double et0 = enevec[0] / cosh(etavec[0]);
        double px0 = et0 * cos(phivec[0]);
        double py0 = et0 * sin(phivec[0]);
        double pz0 = et0 * sinh(etavec[0]);
        double et1 = enevec[1] / cosh(etavec[1]);
        double px1 = et1 * cos(phivec[1]);
        double py1 = et1 * sin(phivec[1]);
        double pz1 = et1 * sinh(etavec[1]);
 
        double dijetmass2 = (enevec[0] + enevec[1]) * (enevec[0] + enevec[1]) - (px1 + px0) * (px1 + px0) -
                            (py1 + py0) * (py1 + py0) - (pz1 + pz0) * (pz1 + pz0);
 
        LogDebug("ValidHcal") << " JetAnalysis ===> Di-jet massSQ " << dijetmass2;
 
        double dijetmass;
        if (dijetmass2 >= 0.)
          dijetmass = sqrt(dijetmass2);
        else
          dijetmass = -sqrt(-1. * dijetmass2);
 
        product.fillDiJets(dijetmass);
 
        LogDebug("ValidHcal") << " JetAnalysis ===> after fillDiJets";
      }
    }
  }
}

References collectEnergyRdir(), funct::cos(), MillePedeFileConverter_cfg::e, een, eta0, PHcalValidInfoJets::fillDiJets(), PHcalValidInfoJets::fillEcollectJet(), PHcalValidInfoJets::fillEtaPhiProfileJet(), PHcalValidInfoJets::fillJets(), PHcalValidInfoJets::fillTProfileJet(), SimG4HcalHitJetFinder::getClusters(), HcalBarrel, HcalEndcap, HcalForward, hcalOnly, es_hardcode_cfi::he, hitcache, hfClusterShapes_cfi::hits, electrons_cff::hoe, mps_fire::i, jetf, jetThreshold, LogDebug, phi0, alignCSCRings::r, SimG4HcalHitJetFinder::rDist(), mps_fire::result, SimG4HcalHitJetFinder::setInput(), funct::sin(), mathSSE::sqrt(), and OrderedSet::t.

Referenced by produce().

layerAnalysis()

void SimG4HcalValidation::layerAnalysis ( PHcalValidInfoLayer &  product )

private

Definition at line 398 of file SimG4HcalValidation.cc.

                                                                    {
  int i = 0;
  LogDebug("ValidHcal") << "\n ===>>> SimG4HcalValidation: Energy deposit "
                        << "in MeV\n at EB : " << std::setw(6) << edepEB / MeV << "\n at EE : " << std::setw(6)
                        << edepEE / MeV << "\n at HB : " << std::setw(6) << edepHB / MeV
                        << "\n at HE : " << std::setw(6) << edepHE / MeV << "\n at HO : " << std::setw(6)
                        << edepHO / MeV << "\n ---- SimG4HcalValidation: Energy deposit in";
  for (i = 0; i < 5; i++)
    LogDebug("ValidHcal") << " Depth " << std::setw(2) << i << " E " << std::setw(8) << edepd[i] / MeV << " MeV";
  LogDebug("ValidHcal") << " ---- SimG4HcalValidation: Energy deposit in"
                        << "layers";
  for (i = 0; i < 20; i++)
    LogDebug("ValidHcal") << " Layer " << std::setw(2) << i << " E " << std::setw(8) << edepl[i] / MeV << " MeV";
 
  product.fillLayers(edepl, edepd, edepHO, edepHB + edepHE, edepEB + edepEE);
 
  // Hits in HF
  product.fillHF(vhitec, vhithc, enEcal, enHcal);
  LogDebug("ValidHcal") << "SimG4HcalValidation::HF hits " << vhitec << " in EC and " << vhithc << " in HC\n"
                        << "                     HB/HE   " << enEcal << " in EC and " << enHcal << " in HC";
 
  // Another HCAL hist to porcess and store separately (a bit more complicated)
  fetchHits(product);
 
  LogDebug("ValidHcal") << " layerAnalysis ===> after fetchHits";
}

References edepd, edepEB, edepEE, edepHB, edepHE, edepHO, edepl, enEcal, enHcal, fetchHits(), PHcalValidInfoLayer::fillHF(), PHcalValidInfoLayer::fillLayers(), mps_fire::i, LogDebug, MeV, vhitec, and vhithc.

Referenced by produce().

nxNAnalysis()

void SimG4HcalValidation::nxNAnalysis ( PHcalValidInfoNxN &  product )

private

Definition at line 426 of file SimG4HcalValidation.cc.

                                                                {
  std::vector<CaloHit> *hits = &hitcache;
  std::vector<CaloHit>::iterator hit_itr;
 
  LogDebug("ValidHcal") << "SimG4HcalValidation::NxNAnalysis : entrance ";
 
  collectEnergyRdir(eta0, phi0);  // HCAL and ECAL energy in  SimHitCache
                                  // around (eta0,phi0)
 
  LogDebug("ValidHcal") << " NxNAnalysis : coolectEnergyRdir - Ecal " << een << "   Hcal " << hen;
 
  double etot = 0.;  // total e deposited     in "cluster"
  double ee = 0.;    // ECAL  e deposited     in "cluster"
  double he = 0.;    // HCAL  e deposited     in "cluster"
  double hoe = 0.;   // HO    e deposited     in "cluster"
 
  int max = dEta.size();  // 4
 
  for (hit_itr = hits->begin(); hit_itr < hits->end(); hit_itr++) {
    double e = hit_itr->e();
    double t = hit_itr->t();
    double eta = hit_itr->eta();
    double phi = hit_itr->phi();
    int type = hit_itr->det();
    int layer = hit_itr->layer();
 
    // NxN calulation
 
    if (fabs(eta0 - eta) <= dEta[max - 1] && fabs(phi0 - phi) <= dPhi[max - 1]) {
      etot += e;
      if (type == 10 || type == 11 || type == 12)
        ee += e;
      if (type == static_cast<int>(HcalBarrel) || type == static_cast<int>(HcalEndcap) ||
          type == static_cast<int>(HcalForward)) {
        he += e;
        if (type == static_cast<int>(HcalBarrel) && layer > 17)
          hoe += e;
 
        // which concrete i-th square ?
        for (int i = 0; i < max; i++) {
          if ((eta0 - eta) <= dEta[i] && (eta0 - eta) >= -dEta[i] && (phi0 - phi) <= dPhi[i] &&
              (phi0 - phi) >= -dPhi[i]) {
            LogDebug("ValidHcal") << "SimG4HcalValidation:: nxNAnalysis eta0,"
                                  << " phi0 = " << eta0 << " " << phi0 << "   type, layer = " << type << "," << layer
                                  << "  eta, phi = " << eta << " " << phi;
 
            product.fillTProfileNxN(e, i, t);
            break;
          }
        }
        // here comes break ...
      }
    }
  }
 
  product.fillEcollectNxN(ee, he, hoe, etot);
  product.fillHvsE(een, hen, hoen, een + hen);
 
  LogDebug("ValidHcal") << " nxNAnalysis ===> after fillHvsE";
}

References collectEnergyRdir(), dEta, dPhi, MillePedeFileConverter_cfg::e, een, PVValHelper::eta, eta0, PHcalValidInfoNxN::fillEcollectNxN(), PHcalValidInfoNxN::fillHvsE(), PHcalValidInfoNxN::fillTProfileNxN(), HcalBarrel, HcalEndcap, HcalForward, es_hardcode_cfi::he, hen, hitcache, hfClusterShapes_cfi::hits, electrons_cff::hoe, hoen, mps_fire::i, LogDebug, SiStripPI::max, phi, phi0, and OrderedSet::t.

Referenced by produce().

operator=()

const SimG4HcalValidation& SimG4HcalValidation::operator= ( const SimG4HcalValidation &  )

privatedelete

produce()

void SimG4HcalValidation::produce ( edm::Event &  e, const edm::EventSetup &    )

overridevirtual

Implements SimProducer.

Definition at line 89 of file SimG4HcalValidation.cc.

                                                                    {
  std::unique_ptr<PHcalValidInfoLayer> productLayer(new PHcalValidInfoLayer);
  layerAnalysis(*productLayer);
  e.put(std::move(productLayer), labelLayer);
 
  if (infolevel > 0) {
    std::unique_ptr<PHcalValidInfoNxN> productNxN(new PHcalValidInfoNxN);
    nxNAnalysis(*productNxN);
    e.put(std::move(productNxN), labelNxN);
  }
 
  if (infolevel > 1) {
    std::unique_ptr<PHcalValidInfoJets> productJets(new PHcalValidInfoJets);
    jetAnalysis(*productJets);
    e.put(std::move(productJets), labelJets);
  }
}

References MillePedeFileConverter_cfg::e, infolevel, jetAnalysis(), labelJets, labelLayer, labelNxN, layerAnalysis(), eostools::move(), and nxNAnalysis().

update() [1/5]

void SimG4HcalValidation::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 182 of file SimG4HcalValidation.cc.

                                                        {
  int i = 0;
  edepEB = edepEE = edepHB = edepHE = edepHO = 0.;
  for (i = 0; i < 5; i++)
    edepd[i] = 0.;
  for (i = 0; i < 20; i++)
    edepl[i] = 0.;
  vhitec = vhithc = enEcal = enHcal = 0;
  // Cache reset
  clear();
 
  int iev = (*evt)()->GetEventID();
  LogDebug("ValidHcal") << "SimG4HcalValidation: =====> Begin of event = " << iev;
}

References clear(), edepd, edepEB, edepEE, edepHB, edepHE, edepHO, edepl, enEcal, enHcal, mps_fire::i, LogDebug, vhitec, and vhithc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

update() [2/5]

void SimG4HcalValidation::update ( const BeginOfJob *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 140 of file SimG4HcalValidation.cc.

                                                      {
  // Numbering From DDD
  edm::ESHandle<HcalDDDSimConstants> hdc;
  (*job)()->get<HcalSimNumberingRecord>().get(hdc);
  const HcalDDDSimConstants *hcons = hdc.product();
  edm::LogVerbatim("ValidHcal") << "HcalTestAnalysis:: Initialise "
                                << "HcalNumberingFromDDD";
  numberingFromDDD = new HcalNumberingFromDDD(hcons);
 
  // Numbering scheme
  org = new HcalTestNumberingScheme(false);
}

References numberingFromDDD, org, and edm::ESHandle< T >::product().

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

update() [3/5]

void SimG4HcalValidation::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 153 of file SimG4HcalValidation.cc.

                                                      {
  int irun = (*run)()->GetRunID();
 
  edm::LogVerbatim("ValidHcal") << " =====> Begin of Run = " << irun;
 
  std::string sdname = names[0];
  G4SDManager *sd = G4SDManager::GetSDMpointerIfExist();
  if (sd != nullptr) {
    G4VSensitiveDetector *aSD = sd->FindSensitiveDetector(sdname);
    if (aSD == nullptr) {
      edm::LogWarning("ValidHcal") << "SimG4HcalValidation::beginOfRun: No SD"
                                   << " with name " << sdname << " in this "
                                   << "Setup";
    } else {
      HCalSD *theCaloSD = dynamic_cast<HCalSD *>(aSD);
      edm::LogVerbatim("ValidHcal") << "SimG4HcalValidation::beginOfRun: Finds SD with name " << theCaloSD->GetName()
                                    << " in this Setup";
      if (org) {
        theCaloSD->setNumberingScheme(org);
        edm::LogVerbatim("ValidHcal") << "SimG4HcalValidation::beginOfRun: set a new numbering scheme";
      }
    }
  } else {
    edm::LogWarning("ValidHcal") << "SimG4HcalValidation::beginOfRun: Could "
                                 << "not get SD Manager!";
  }
}

References names, org, sd, HCalSD::setNumberingScheme(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

update() [4/5]

void SimG4HcalValidation::update ( const EndOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 256 of file SimG4HcalValidation.cc.

                                                      {
  count++;
 
  // Fill hits cache for jetfinding etc.
  fill(evt);
  LogDebug("ValidHcal") << "SimG4HcalValidation:: ---  after Fill ";
}

References count, fill(), and LogDebug.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

update() [5/5]

void SimG4HcalValidation::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 198 of file SimG4HcalValidation.cc.

                                                    {
  if (aStep != nullptr) {
    G4VPhysicalVolume *curPV = aStep->GetPreStepPoint()->GetPhysicalVolume();
    G4String name = curPV->GetName();
    name.assign(name, 0, 3);
    double edeposit = aStep->GetTotalEnergyDeposit();
    int layer = -1, depth = -1;
    if (name == "EBR") {
      depth = 0;
      edepEB += edeposit;
    } else if (name == "EFR") {
      depth = 0;
      edepEE += edeposit;
    } else if (name == "HBS") {
      layer = (curPV->GetCopyNo() / 10) % 100;
      depth = (curPV->GetCopyNo()) % 10 + 1;
      if (depth > 0 && depth < 4 && layer >= 0 && layer < 17) {
        edepHB += edeposit;
      } else {
        edm::LogWarning("ValidHcal") << "SimG4HcalValidation:Error " << curPV->GetName() << curPV->GetCopyNo();
        depth = -1;
        layer = -1;
      }
    } else if (name == "HES") {
      layer = (curPV->GetCopyNo() / 10) % 100;
      depth = (curPV->GetCopyNo()) % 10 + 1;
      if (depth > 0 && depth < 3 && layer >= 0 && layer < 19) {
        edepHE += edeposit;
      } else {
        edm::LogWarning("ValidHcal") << "SimG4HcalValidation:Error " << curPV->GetName() << curPV->GetCopyNo();
        depth = -1;
        layer = -1;
      }
    } else if (name == "HTS") {
      layer = (curPV->GetCopyNo() / 10) % 100;
      depth = (curPV->GetCopyNo()) % 10 + 1;
      if (depth > 3 && depth < 5 && layer >= 17 && layer < 20) {
        edepHO += edeposit;
      } else {
        edm::LogWarning("ValidHcal") << "SimG4HcalValidation:Error " << curPV->GetName() << curPV->GetCopyNo();
        depth = -1;
        layer = -1;
      }
    }
    if (depth >= 0 && depth < 5)
      edepd[depth] += edeposit;
    if (layer >= 0 && layer < 20)
      edepl[layer] += edeposit;
 
    if (layer >= 0 && layer < 20) {
      LogDebug("ValidHcal") << "SimG4HcalValidation:: G4Step: " << name << " Layer " << std::setw(3) << layer
                            << " Depth " << std::setw(2) << depth << " Edep " << std::setw(6) << edeposit / MeV
                            << " MeV";
    }
  }
}

References LEDCalibrationChannels::depth, edepd, edepEB, edepEE, edepHB, edepHE, edepHO, edepl, LogDebug, MeV, and Skims_PA_cff::name.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

Member Data Documentation

applySampling

bool SimG4HcalValidation::applySampling

private

Definition at line 91 of file SimG4HcalValidation.h.

Referenced by fill(), and SimG4HcalValidation().

coneSize

double SimG4HcalValidation::coneSize

private

Definition at line 89 of file SimG4HcalValidation.h.

Referenced by collectEnergyRdir(), init(), and SimG4HcalValidation().

count

unsigned int SimG4HcalValidation::count

private

Definition at line 100 of file SimG4HcalValidation.h.

Referenced by init(), update(), and ~SimG4HcalValidation().

dEta

std::vector<double> SimG4HcalValidation::dEta

private

Definition at line 96 of file SimG4HcalValidation.h.

Referenced by init(), and nxNAnalysis().

dPhi

std::vector<double> SimG4HcalValidation::dPhi

private

Definition at line 97 of file SimG4HcalValidation.h.

Referenced by init(), and nxNAnalysis().

edepd

double SimG4HcalValidation::edepd[5]

private

Definition at line 102 of file SimG4HcalValidation.h.

Referenced by layerAnalysis(), and update().

edepEB

double SimG4HcalValidation::edepEB

private

Definition at line 101 of file SimG4HcalValidation.h.

Referenced by layerAnalysis(), and update().

edepEE

double SimG4HcalValidation::edepEE

private

Definition at line 101 of file SimG4HcalValidation.h.

Referenced by layerAnalysis(), and update().

edepHB

double SimG4HcalValidation::edepHB

private

Definition at line 101 of file SimG4HcalValidation.h.

Referenced by layerAnalysis(), and update().

edepHE

double SimG4HcalValidation::edepHE

private

Definition at line 101 of file SimG4HcalValidation.h.

Referenced by layerAnalysis(), and update().

edepHO

double SimG4HcalValidation::edepHO

private

Definition at line 101 of file SimG4HcalValidation.h.

Referenced by layerAnalysis(), and update().

edepl

double SimG4HcalValidation::edepl[20]

private

Definition at line 102 of file SimG4HcalValidation.h.

Referenced by layerAnalysis(), and update().

een

double SimG4HcalValidation::een

private

Definition at line 103 of file SimG4HcalValidation.h.

Referenced by collectEnergyRdir(), jetAnalysis(), and nxNAnalysis().

ehitThreshold

double SimG4HcalValidation::ehitThreshold

private

Definition at line 89 of file SimG4HcalValidation.h.

Referenced by fill(), and SimG4HcalValidation().

enEcal

double SimG4HcalValidation::enEcal

private

Definition at line 104 of file SimG4HcalValidation.h.

Referenced by fill(), layerAnalysis(), and update().

enHcal

double SimG4HcalValidation::enHcal

private

Definition at line 104 of file SimG4HcalValidation.h.

Referenced by fill(), layerAnalysis(), and update().

eta0

float SimG4HcalValidation::eta0

private

Definition at line 90 of file SimG4HcalValidation.h.

Referenced by collectEnergyRdir(), jetAnalysis(), nxNAnalysis(), and SimG4HcalValidation().

hcalOnly

bool SimG4HcalValidation::hcalOnly

private

Definition at line 91 of file SimG4HcalValidation.h.

Referenced by fill(), jetAnalysis(), and SimG4HcalValidation().

hen

double SimG4HcalValidation::hen

private

Definition at line 103 of file SimG4HcalValidation.h.

Referenced by collectEnergyRdir(), and nxNAnalysis().

hhitThreshold

double SimG4HcalValidation::hhitThreshold

private

Definition at line 89 of file SimG4HcalValidation.h.

Referenced by fill(), and SimG4HcalValidation().

hitcache

std::vector<CaloHit> SimG4HcalValidation::hitcache

private

Definition at line 80 of file SimG4HcalValidation.h.

Referenced by clear(), collectEnergyRdir(), fetchHits(), fill(), jetAnalysis(), and nxNAnalysis().

hoen

double SimG4HcalValidation::hoen

private

Definition at line 103 of file SimG4HcalValidation.h.

Referenced by collectEnergyRdir(), and nxNAnalysis().

infolevel

int SimG4HcalValidation::infolevel

private

Definition at line 92 of file SimG4HcalValidation.h.

Referenced by produce(), and SimG4HcalValidation().

jetf

SimG4HcalHitJetFinder* SimG4HcalValidation::jetf

private

Definition at line 71 of file SimG4HcalValidation.h.

Referenced by collectEnergyRdir(), init(), jetAnalysis(), and ~SimG4HcalValidation().

jetThreshold

float SimG4HcalValidation::jetThreshold

private

Definition at line 90 of file SimG4HcalValidation.h.

Referenced by jetAnalysis(), and SimG4HcalValidation().

labelJets

std::string SimG4HcalValidation::labelJets

private

Definition at line 93 of file SimG4HcalValidation.h.

Referenced by produce(), and SimG4HcalValidation().

labelLayer

std::string SimG4HcalValidation::labelLayer

private

Definition at line 93 of file SimG4HcalValidation.h.

Referenced by produce(), and SimG4HcalValidation().

labelNxN

std::string SimG4HcalValidation::labelNxN

private

Definition at line 93 of file SimG4HcalValidation.h.

Referenced by produce(), and SimG4HcalValidation().

names

std::vector<std::string> SimG4HcalValidation::names

private

Definition at line 88 of file SimG4HcalValidation.h.

Referenced by fill(), SimG4HcalValidation(), and update().

numberingFromDDD

HcalNumberingFromDDD* SimG4HcalValidation::numberingFromDDD

private

Definition at line 74 of file SimG4HcalValidation.h.

Referenced by fill(), update(), and ~SimG4HcalValidation().

org

HcalTestNumberingScheme* SimG4HcalValidation::org

private

Definition at line 77 of file SimG4HcalValidation.h.

Referenced by fill(), update(), and ~SimG4HcalValidation().

phi0

float SimG4HcalValidation::phi0

private

Definition at line 90 of file SimG4HcalValidation.h.

Referenced by collectEnergyRdir(), jetAnalysis(), nxNAnalysis(), and SimG4HcalValidation().

scaleHB

std::vector<float> SimG4HcalValidation::scaleHB

private

Definition at line 83 of file SimG4HcalValidation.h.

Referenced by getHcalScale(), and init().

scaleHE

std::vector<float> SimG4HcalValidation::scaleHE

private

Definition at line 84 of file SimG4HcalValidation.h.

Referenced by getHcalScale(), and init().

scaleHF

std::vector<float> SimG4HcalValidation::scaleHF

private

Definition at line 85 of file SimG4HcalValidation.h.

Referenced by getHcalScale(), and init().

timeLowlim

float SimG4HcalValidation::timeLowlim

private

Definition at line 90 of file SimG4HcalValidation.h.

Referenced by fill(), and SimG4HcalValidation().

timeUplim

float SimG4HcalValidation::timeUplim

private

Definition at line 90 of file SimG4HcalValidation.h.

Referenced by SimG4HcalValidation().

vhitec

double SimG4HcalValidation::vhitec

private

Definition at line 104 of file SimG4HcalValidation.h.

Referenced by fill(), layerAnalysis(), and update().

vhithc

double SimG4HcalValidation::vhithc

private

Definition at line 104 of file SimG4HcalValidation.h.

Referenced by fill(), layerAnalysis(), and update().