CaloTrkProcessing Class Reference

#include <CaloTrkProcessing.h>

Inheritance diagram for CaloTrkProcessing:

Classes
struct	Detector

Public Member Functions
	CaloTrkProcessing (const std::string &aSDname, const edm::EventSetup &es, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *)
void	clearHits () override
void	EndOfEvent (G4HCofThisEvent *) override
void	fillHits (edm::PCaloHitContainer &, const std::string &) override
void	Initialize (G4HCofThisEvent *) override
bool	ProcessHits (G4Step *, G4TouchableHistory *) override
uint32_t	setDetUnitId (const G4Step *step) override
	~CaloTrkProcessing () override
Public Member Functions inherited from SensitiveCaloDetector
virtual void	reset ()
	SensitiveCaloDetector (const std::string &iname, const edm::EventSetup &es, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)
Public Member Functions inherited from SensitiveDetector
void	EndOfEvent (G4HCofThisEvent *eventHC) override
const std::vector< std::string > &	getNames () const
void	Initialize (G4HCofThisEvent *eventHC) override
bool	isCaloSD () const
	SensitiveDetector (const std::string &iname, const edm::EventSetup &es, const SensitiveDetectorCatalog &, edm::ParameterSet const &p, bool calo)
	~SensitiveDetector () override
Public Member Functions inherited from Observer< const BeginOfEvent * >
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const G4Step * >
	Observer ()
void	slotForUpdate (const G4Step * iT)
virtual	~Observer ()

Private Member Functions
void	detectorLevel (const G4VTouchable *, int &, int *, G4String *) const
int	detLevels (const G4VTouchable *) const
G4LogicalVolume *	detLV (const G4VTouchable *, int) const
int	isItCalo (const G4VTouchable *, const std::vector< Detector > &)
int	isItInside (const G4VTouchable *, int, int)
void	update (const BeginOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const G4Step *) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
std::vector< Detector >	detectors_
bool	doFineCalo_
double	eMin_
double	eMinFine_
double	eMinFinePhoton_
std::vector< Detector >	fineDetectors_
int	lastTrackID_
bool	putHistory_
bool	testBeam_

Additional Inherited Members
Protected Types inherited from SensitiveDetector
enum	coordinates { WorldCoordinates, LocalCoordinates }
Protected Member Functions inherited from SensitiveDetector
TrackInformation *	cmsTrackInformation (const G4Track *aTrack)
Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point) const
Local3DPoint	FinalStepPosition (const G4Step *step, coordinates) const
Local3DPoint	InitialStepPosition (const G4Step *step, coordinates) const
Local3DPoint	LocalPostStepPosition (const G4Step *step) const
Local3DPoint	LocalPreStepPosition (const G4Step *step) const
void	NaNTrap (const G4Step *step) const
void	setNames (const std::vector< std::string > &)

Detailed Description

Definition at line 22 of file CaloTrkProcessing.h.

Constructor & Destructor Documentation

CaloTrkProcessing::CaloTrkProcessing	(	const std::string &	aSDname,
		const edm::EventSetup &	es,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *
	)

Definition at line 25 of file CaloTrkProcessing.cc.

References CaloSimulationParameters::caloNames_, hgcalTestNeighbor_cfi::detector, detectors_, doFineCalo_, eMin_, eMinFine_, eMinFinePhoton_, Exception, CaloSimulationParameters::fCaloNames_, fineDetectors_, CaloSimulationParameters::fLevels_, CaloTrkProcessing::Detector::fromDetL, CaloTrkProcessing::Detector::fromDets, CaloTrkProcessing::Detector::fromLevels, edm::EventSetup::get(), edm::ParameterSet::getParameter(), mps_fire::i, CaloSimulationParameters::insideLevel_, CaloSimulationParameters::insideNames_, edm::ESHandleBase::isValid(), dqmdumpme::k, CaloTrkProcessing::Detector::level, CaloSimulationParameters::levels_, CaloTrkProcessing::Detector::lv, MeV, Skims_PA_cff::name, CaloTrkProcessing::Detector::name, CaloSimulationParameters::neighbours_, edm::ESHandle< T >::product(), putHistory_, and testBeam_.

    : SensitiveCaloDetector(name, es, clg, p), lastTrackID_(-1) {
  //Initialise the parameter set
  edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("CaloTrkProcessing");
  testBeam_ = m_p.getParameter<bool>("TestBeam");
  eMin_ = m_p.getParameter<double>("EminTrack") * MeV;
  putHistory_ = m_p.getParameter<bool>("PutHistory");
  doFineCalo_ = m_p.getParameter<bool>("DoFineCalo");
  eMinFine_ = m_p.getParameter<double>("EminFineTrack") * MeV;
  eMinFinePhoton_ = m_p.getParameter<double>("EminFinePhoton") * MeV;

  edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: Initailised with TestBeam = " << testBeam_ << " Emin = " << eMin_
                              << ":" << eMinFine_ << ":" << eMinFinePhoton_ << " MeV and Flags " << putHistory_
                              << " (History), " << doFineCalo_ << " (Special Calorimeter)";

  // Get pointers to HcalDDDConstant and HcalSimulationParameters
  edm::ESHandle<CaloSimulationParameters> csps;
  es.get<HcalParametersRcd>().get(csps);
  if (csps.isValid()) {
    const CaloSimulationParameters* csp = csps.product();
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: " << csp->caloNames_.size() << " entries for caloNames:";
    for (unsigned int i = 0; i < csp->caloNames_.size(); i++)
      edm::LogVerbatim("CaloSim") << " (" << i << ") " << csp->caloNames_[i];
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: " << csp->levels_.size() << " entries for levels:";
    for (unsigned int i = 0; i < csp->levels_.size(); i++)
      edm::LogVerbatim("CaloSim") << " (" << i << ") " << csp->levels_[i];
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: " << csp->neighbours_.size() << " entries for neighbours:";
    for (unsigned int i = 0; i < csp->neighbours_.size(); i++)
      edm::LogVerbatim("CaloSim") << " (" << i << ") " << csp->neighbours_[i];
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: " << csp->insideNames_.size() << " entries for insideNames:";
    for (unsigned int i = 0; i < csp->insideNames_.size(); i++)
      edm::LogVerbatim("CaloSim") << " (" << i << ") " << csp->insideNames_[i];
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: " << csp->insideLevel_.size() << " entries for insideLevel:";
    for (unsigned int i = 0; i < csp->insideLevel_.size(); i++)
      edm::LogVerbatim("CaloSim") << " (" << i << ") " << csp->insideLevel_[i];
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: " << csp->fCaloNames_.size()
                                << " entries for fineCalorimeterNames:";
    for (unsigned int i = 0; i < csp->fCaloNames_.size(); i++)
      edm::LogVerbatim("CaloSim") << " (" << i << ") " << csp->fCaloNames_[i];
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: " << csp->fLevels_.size()
                                << " entries for fineCalorimeterLevels:";
    for (unsigned int i = 0; i < csp->fLevels_.size(); i++)
      edm::LogVerbatim("CaloSim") << " (" << i << ") " << csp->fLevels_[i];
#endif

    if (csp->caloNames_.size() < csp->neighbours_.size()) {
      edm::LogError("CaloSim") << "CaloTrkProcessing: # of Calorimeter bins " << csp->caloNames_.size()
                               << " does not match with " << csp->neighbours_.size() << " ==> illegal ";
      throw cms::Exception("Unknown", "CaloTrkProcessing")
          << "Calorimeter array size does not match with size of neighbours\n";
    }

    const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
    std::vector<G4LogicalVolume*>::const_iterator lvcite;
    int istart = 0;
    for (unsigned int i = 0; i < csp->caloNames_.size(); i++) {
      G4LogicalVolume* lv = nullptr;
      G4String name = static_cast<G4String>(csp->caloNames_[i]);
      for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) {
        if ((*lvcite)->GetName() == name) {
          lv = (*lvcite);
          break;
        }
      }
      if (lv != nullptr) {
        CaloTrkProcessing::Detector detector;
        detector.name = name;
        detector.lv = lv;
        detector.level = csp->levels_[i];
        if (istart + csp->neighbours_[i] > static_cast<int>(csp->insideNames_.size())) {
          edm::LogError("CaloSim") << "CaloTrkProcessing: # of InsideNames bins " << csp->insideNames_.size()
                                   << " too few compaerd to " << istart + csp->neighbours_[i]
                                   << " requested ==> illegal ";
          throw cms::Exception("Unknown", "CaloTrkProcessing")
              << "InsideNames array size does not match with list of neighbours\n";
        }
        std::vector<std::string> inside;
        std::vector<G4LogicalVolume*> insideLV;
        std::vector<int> insideLevels;
        for (int k = 0; k < csp->neighbours_[i]; k++) {
          lv = nullptr;
          name = static_cast<G4String>(csp->insideNames_[istart + k]);
          for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) {
            if ((*lvcite)->GetName() == name) {
              lv = (*lvcite);
              break;
            }
          }
          inside.push_back(name);
          insideLV.push_back(lv);
          insideLevels.push_back(csp->insideLevel_[istart + k]);
        }
        detector.fromDets = inside;
        detector.fromDetL = insideLV;
        detector.fromLevels = insideLevels;
        detectors_.emplace_back(detector);
      }
      istart += csp->neighbours_[i];
    }

    for (unsigned int i = 0; i < csp->fCaloNames_.size(); i++) {
      G4LogicalVolume* lv = nullptr;
      G4String name = static_cast<G4String>(csp->fCaloNames_[i]);
      for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) {
        if ((*lvcite)->GetName() == name) {
          lv = (*lvcite);
          break;
        }
      }
      if (lv != nullptr) {
        CaloTrkProcessing::Detector detector;
        detector.name = name;
        detector.lv = lv;
        detector.level = csp->fLevels_[i];
        detector.fromDets.clear();
        detector.fromDetL.clear();
        detector.fromLevels.clear();
        fineDetectors_.emplace_back(detector);
      }
    }
  } else {
    edm::LogError("HcalSim") << "CaloTrkProcessing: Cannot find CaloSimulationParameters";
    throw cms::Exception("Unknown", "CaloSD") << "Cannot find CaloSimulationParameters\n";
  }

  edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: with " << detectors_.size() << " calorimetric volumes";
  for (unsigned int i = 0; i < detectors_.size(); i++) {
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: Calorimeter volume " << i << " " << detectors_[i].name << " LV "
                                << detectors_[i].lv << " at level " << detectors_[i].level << " with "
                                << detectors_[i].fromDets.size() << " neighbours";
    for (unsigned int k = 0; k < detectors_[i].fromDets.size(); k++)
      edm::LogVerbatim("CaloSim") << "                   Element " << k << " " << detectors_[i].fromDets[k] << " LV "
                                  << detectors_[i].fromDetL[k] << " at level " << detectors_[i].fromLevels[k];
  }

  doFineCalo_ = !(fineDetectors_.empty());
  edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: with " << fineDetectors_.size() << " special calorimetric volumes";
  for (unsigned int i = 0; i < detectors_.size(); i++)
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: Calorimeter volume " << i << " " << detectors_[i].name << " LV "
                                << detectors_[i].lv << " at level " << detectors_[i].level;
}

CaloTrkProcessing::~CaloTrkProcessing ( )

override

Definition at line 172 of file CaloTrkProcessing.cc.

{}

Member Function Documentation

void CaloTrkProcessing::clearHits ( )

inlineoverridevirtual

Implements SensitiveDetector.

Definition at line 33 of file CaloTrkProcessing.h.

{}

void CaloTrkProcessing::detectorLevel ( const G4VTouchable *  touch, int &  level, int *  copyno, G4String *  name  ) const

private

Definition at line 372 of file CaloTrkProcessing.cc.

References mps_fire::i, cuy::ii, personalPlayback::level, MetAnalyzer::pv(), AlCaHLTBitMon_QueryRunRegistry::string, and susybsm::HSCParticleType::unknown.

Referenced by isItCalo(), and isItInside().

                                                                                                              {
  static const std::string unknown("Unknown");
  //Get name and copy numbers
  if (level > 0) {
    for (int ii = 0; ii < level; ii++) {
      int i = level - ii - 1;
      G4VPhysicalVolume* pv = touch->GetVolume(i);
      if (pv != nullptr)
        name[ii] = pv->GetName();
      else
        name[ii] = unknown;
      copyno[ii] = touch->GetReplicaNumber(i);
    }
  }
}

int CaloTrkProcessing::detLevels ( const G4VTouchable *  touch ) const

private

Definition at line 351 of file CaloTrkProcessing.cc.

Referenced by isItCalo(), and isItInside().

                                                                {
  //Return number of levels
  if (touch)
    return ((touch->GetHistoryDepth()) + 1);
  else
    return 0;
}

G4LogicalVolume * CaloTrkProcessing::detLV ( const G4VTouchable *  touch, int  currentlevel  ) const

private

Definition at line 359 of file CaloTrkProcessing.cc.

References cuy::ii, and personalPlayback::level.

Referenced by isItCalo(), and isItInside().

                                                                                           {
  G4LogicalVolume* lv = nullptr;
  if (touch) {
    int level = ((touch->GetHistoryDepth()) + 1);
    if (level > 0 && level >= currentlevel) {
      int ii = level - currentlevel;
      lv = touch->GetVolume(ii)->GetLogicalVolume();
      return lv;
    }
  }
  return lv;
}

void CaloTrkProcessing::EndOfEvent ( G4HCofThisEvent *  )

inlineoverride

Definition at line 36 of file CaloTrkProcessing.h.

{}

void CaloTrkProcessing::fillHits ( edm::PCaloHitContainer &  , const std::string &    )

inlineoverridevirtual

Implements SensitiveCaloDetector.

Definition at line 37 of file CaloTrkProcessing.h.

{}

void CaloTrkProcessing::Initialize ( G4HCofThisEvent *  )

inlineoverride

Definition at line 32 of file CaloTrkProcessing.h.

{}

int CaloTrkProcessing::isItCalo ( const G4VTouchable *  touch, const std::vector< Detector > &  detectors  )

private

Definition at line 261 of file CaloTrkProcessing.cc.

References detectorLevel(), detLevels(), detLV(), testProducerWithPsetDescEmpty_cfi::i2, personalPlayback::level, jets_cff::levels, hfnoseParametersInitialization_cfi::name2, convertSQLiteXML::ok, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by update().

                                                                                               {
  int lastLevel = -1;
  G4LogicalVolume* lv = nullptr;
  for (unsigned int it = 0; it < detectors.size(); it++) {
    if (lastLevel != detectors[it].level) {
      lastLevel = detectors[it].level;
      lv = detLV(touch, lastLevel);
#ifdef EDM_ML_DEBUG
      std::string name1 = "Unknown";
      if (lv != 0)
        name1 = lv->GetName();
      edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: volume " << name1 << " at Level " << lastLevel;
      int levels = detLevels(touch);
      if (levels > 0) {
        G4String name2[20];
        int copyno2[20];
        detectorLevel(touch, levels, copyno2, name2);
        for (int i2 = 0; i2 < levels; i2++)
          edm::LogVerbatim("CaloSim") << " " << i2 << " " << name2[i2] << " " << copyno2[i2];
      }
#endif
    }
    bool ok = (lv == detectors[it].lv);
    if (ok)
      return it;
  }
  return -1;
}

int CaloTrkProcessing::isItInside ( const G4VTouchable *  touch, int  idcal, int  idin  )

private

Definition at line 290 of file CaloTrkProcessing.cc.

References detectorLevel(), detectors_, detLevels(), detLV(), testProducerWithPsetDescEmpty_cfi::i2, globals_cff::id1, globals_cff::id2, jets_cff::levels, hfnoseParametersInitialization_cfi::name2, convertSQLiteXML::ok, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by update().

                                                                                {
  int lastLevel = -1;
  G4LogicalVolume* lv = nullptr;
  int id1, id2;
  if (idcal < 0) {
    id1 = 0;
    id2 = static_cast<int>(detectors_.size());
  } else {
    id1 = idcal;
    id2 = id1 + 1;
  }
  for (int it1 = id1; it1 < id2; it1++) {
    if (idin < 0) {
      for (unsigned int it2 = 0; it2 < detectors_[it1].fromDets.size(); it2++) {
        if (lastLevel != detectors_[it1].fromLevels[it2]) {
          lastLevel = detectors_[it1].fromLevels[it2];
          lv = detLV(touch, lastLevel);
#ifdef EDM_ML_DEBUG
          std::string name1 = "Unknown";
          if (lv != 0)
            name1 = lv->GetName();
          edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: volume " << name1 << " at Level " << lastLevel;
          int levels = detLevels(touch);
          if (levels > 0) {
            G4String name2[20];
            int copyno2[20];
            detectorLevel(touch, levels, copyno2, name2);
            for (int i2 = 0; i2 < levels; i2++)
              edm::LogVerbatim("CaloSim") << " " << i2 << " " << name2[i2] << " " << copyno2[i2];
          }
#endif
        }
        bool ok = (lv == detectors_[it1].fromDetL[it2]);
        if (ok)
          return it2;
      }
    } else {
      lastLevel = detectors_[it1].fromLevels[idin];
      lv = detLV(touch, lastLevel);
#ifdef EDM_ML_DEBUG
      std::string name1 = "Unknown";
      if (lv != 0)
        name1 = lv->GetName();
      edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: volume " << name1 << " at Level " << lastLevel;
      int levels = detLevels(touch);
      if (levels > 0) {
        G4String name2[20];
        int copyno2[20];
        detectorLevel(touch, levels, copyno2, name2);
        for (int i2 = 0; i2 < levels; i2++)
          edm::LogVerbatim("CaloSim") << " " << i2 << " " << name2[i2] << " " << copyno2[i2];
      }
#endif
      bool ok = (lv == detectors_[it1].fromDetL[idin]);
      if (ok)
        return idin;
    }
  }
  return -1;
}

bool CaloTrkProcessing::ProcessHits ( G4Step *  , G4TouchableHistory *    )

inlineoverridevirtual

Implements SensitiveDetector.

Definition at line 34 of file CaloTrkProcessing.h.

{ return true; }

uint32_t CaloTrkProcessing::setDetUnitId ( const G4Step *  step )

inlineoverridevirtual

Implements SensitiveDetector.

Definition at line 35 of file CaloTrkProcessing.h.

{ return 0; }

void CaloTrkProcessing::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 174 of file CaloTrkProcessing.cc.

References lastTrackID_.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{ lastTrackID_ = -1; }

void CaloTrkProcessing::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 176 of file CaloTrkProcessing.cc.

References TrackInformation::caloIDChecked(), TkAlMuonSelectors_cfi::cut, detectors_, doFineCalo_, eMin_, eMinFine_, eMinFinePhoton_, Exception, fineDetectors_, TrackInformation::getIDCaloVolume(), TrackInformation::getIDLastVolume(), TrackInformation::getIDonCaloSurface(), triggerObjects_cff::id, TrackInformation::isInHistory(), isItCalo(), isItInside(), lastTrackID_, MeV, putHistory_, TrackInformation::putInHistory(), TrackInformation::setCaloIDChecked(), TrackInformation::setIDonCaloSurface(), and testBeam_.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                  {
  // define if you are at the surface of CALO

  G4Track* theTrack = aStep->GetTrack();
  int id = theTrack->GetTrackID();

  TrackInformation* trkInfo = dynamic_cast<TrackInformation*>(theTrack->GetUserInformation());

  if (trkInfo == nullptr) {
    edm::LogError("CaloSim") << "CaloTrkProcessing: No trk info !!!! abort ";
    throw cms::Exception("Unknown", "CaloTrkProcessing") << "cannot get trkInfo for Track " << id << "\n";
  }

  if (testBeam_) {
    if (trkInfo->getIDonCaloSurface() == 0) {
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("CaloSim") << "CaloTrkProcessing set IDonCaloSurface to " << id << " at step Number "
                                  << theTrack->GetCurrentStepNumber();
#endif
      trkInfo->setIDonCaloSurface(id, 0, 0, theTrack->GetDefinition()->GetPDGEncoding(), theTrack->GetMomentum().mag());
      lastTrackID_ = id;
      if (theTrack->GetKineticEnergy() / MeV > eMin_)
        trkInfo->putInHistory();
    }
  } else {
    if (putHistory_) {
      trkInfo->putInHistory();
      //      trkInfo->setAncestor();
    }
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("CaloSim") << "CaloTrkProcessing Entered for " << id << " at stepNumber "
                                << theTrack->GetCurrentStepNumber() << " IDonCaloSur.. "
                                << trkInfo->getIDonCaloSurface() << " CaloCheck " << trkInfo->caloIDChecked();
#endif
    if (trkInfo->getIDonCaloSurface() != 0) {
      if (trkInfo->caloIDChecked() == false) {
        G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
        const G4VTouchable* post_touch = postStepPoint->GetTouchable();

        if (isItInside(post_touch, trkInfo->getIDCaloVolume(), trkInfo->getIDLastVolume()) > 0) {
          trkInfo->setIDonCaloSurface(0, -1, -1, 0, 0);
        } else {
          trkInfo->setCaloIDChecked(true);
        }
      }
    } else {
      G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
      const G4VTouchable* post_touch = postStepPoint->GetTouchable();
      int ical = isItCalo(post_touch, detectors_);
      if (ical >= 0) {
        G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
        const G4VTouchable* pre_touch = preStepPoint->GetTouchable();
        int inside = isItInside(pre_touch, ical, -1);
        if (inside >= 0 || (theTrack->GetCurrentStepNumber() == 1)) {
          trkInfo->setIDonCaloSurface(
              id, ical, inside, theTrack->GetDefinition()->GetPDGEncoding(), theTrack->GetMomentum().mag());
          trkInfo->setCaloIDChecked(true);
          lastTrackID_ = id;
          if (theTrack->GetKineticEnergy() / MeV > eMin_)
            trkInfo->putInHistory();
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: set ID on Calo " << ical << " surface (Inside " << inside
                                      << ") to " << id << " of a Track with Kinetic Energy "
                                      << theTrack->GetKineticEnergy() / MeV << " MeV";
#endif
        }
      }
    }
  }
  if (doFineCalo_ && (!trkInfo->isInHistory())) {
    const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
    if (isItCalo(touch, fineDetectors_) >= 0) {
      int pdg = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
      double cut = (pdg == 22) ? eMinFinePhoton_ : eMinFine_;
      if (aStep->GetTrack()->GetKineticEnergy() / MeV > cut)
        trkInfo->putInHistory();
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("CaloSim") << "CaloTrkProcessing: the track with PDGID " << pdg << " and kinetic energy "
                                  << aStep->GetTrack()->GetKineticEnergy() / MeV << " is tested against " << cut
                                  << " to be put in history";
#endif
    }
  }
}

Member Data Documentation

std::vector<Detector> CaloTrkProcessing::detectors_

private

Definition at line 63 of file CaloTrkProcessing.h.

Referenced by CaloTrkProcessing(), isItInside(), and update().

bool CaloTrkProcessing::doFineCalo_

private

Definition at line 60 of file CaloTrkProcessing.h.

Referenced by CaloTrkProcessing(), and update().

double CaloTrkProcessing::eMin_

private

Definition at line 61 of file CaloTrkProcessing.h.

Referenced by CaloTrkProcessing(), and update().

double CaloTrkProcessing::eMinFine_

private

Definition at line 61 of file CaloTrkProcessing.h.

Referenced by CaloTrkProcessing(), and update().

double CaloTrkProcessing::eMinFinePhoton_

private

Definition at line 61 of file CaloTrkProcessing.h.

Referenced by CaloTrkProcessing(), and update().

std::vector<Detector> CaloTrkProcessing::fineDetectors_

private

Definition at line 63 of file CaloTrkProcessing.h.

Referenced by CaloTrkProcessing(), and update().

int CaloTrkProcessing::lastTrackID_

private

Definition at line 62 of file CaloTrkProcessing.h.

Referenced by update().

bool CaloTrkProcessing::putHistory_

private

Definition at line 60 of file CaloTrkProcessing.h.

Referenced by CaloTrkProcessing(), and update().

bool CaloTrkProcessing::testBeam_

private

Definition at line 60 of file CaloTrkProcessing.h.

Referenced by CaloTrkProcessing(), and update().