ShowerRecoTools::ShowerDirectionCheater Class Reference

Inheritance diagram for ShowerRecoTools::ShowerDirectionCheater:

Public Member Functions
	ShowerDirectionCheater (const fhicl::ParameterSet &pset)
int	CalculateElement (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder) override

Private Member Functions
int	RunShowerTool (const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder, std::string evd_display_name_append="")
virtual void	InitialiseProducers ()
void	SetPtr (art::ProducesCollector *collector)
void	InitaliseProducerPtr (reco::shower::ShowerProducedPtrsHolder &uniqueproducerPtrs)
virtual int	AddAssociations (const art::Ptr< recob::PFParticle > &pfpPtr, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder)
const shower::LArPandoraShowerAlg &	GetLArPandoraShowerAlg () const
template<class T >
art::Ptr< T >	GetProducedElementPtr (std::string Name, reco::shower::ShowerElementHolder &ShowerEleHolder, int iter=-1)
template<class T >
void	InitialiseProduct (std::string Name, std::string InstanceName="")
template<class T , class A , class B >
void	AddSingle (A &a, B &b, std::string Name)
int	GetVectorPtrSize (std::string Name)
void	PrintPtrs ()
void	PrintPtr (std::string Name)

Private Attributes
shower::LArPandoraShowerCheatingAlg	fLArPandoraShowerCheatingAlg
art::ServiceHandle< art::TFileService >	tfs
const art::InputTag	fPFParticleLabel
const unsigned int	fNSegments
const bool	fRMSFlip
const bool	fVertexFlip
TTree *	Tree
float	vertexDotProduct
float	rmsGradient
const std::string	fShowerStartPositionInputLabel
const std::string	fTrueParticleInputLabel
const std::string	fShowerDirectionOutputLabel

Detailed Description

Definition at line 24 of file ShowerDirectionCheater_tool.cc.

Constructor & Destructor Documentation

ShowerRecoTools::ShowerDirectionCheater::ShowerDirectionCheater

(

const fhicl::ParameterSet &

pset

)

Definition at line 59 of file ShowerDirectionCheater_tool.cc.

References fRMSFlip, fVertexFlip, rmsGradient, tfs, Tree, and vertexDotProduct.

    : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
    , fLArPandoraShowerCheatingAlg(pset.get<fhicl::ParameterSet>("LArPandoraShowerCheatingAlg"))
    , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
    , fNSegments(pset.get<unsigned int>("NSegments"))
    , fRMSFlip(pset.get<bool>("RMSFlip"))
    , fVertexFlip(pset.get<bool>("VertexFlip"))
    , fShowerStartPositionInputLabel(pset.get<std::string>("ShowerStartPositionInputLabel"))
    , fTrueParticleInputLabel(pset.get<std::string>("TrueParticleInputLabel"))
    , fShowerDirectionOutputLabel(pset.get<std::string>("ShowerDirectionOutputLabel"))
  {
    if (fVertexFlip || fRMSFlip) {
      Tree = tfs->make<TTree>("DebugTreeDirCheater", "DebugTree from shower direction cheater");
      if (fVertexFlip) Tree->Branch("vertexDotProduct", &vertexDotProduct);
      if (fRMSFlip) Tree->Branch("rmsGradient", &rmsGradient);
    }
  }

Member Function Documentation

virtual int ShowerRecoTools::IShowerTool::AddAssociations ( const art::Ptr< recob::PFParticle > &  pfpPtr, art::Event &  Event, reco::shower::ShowerElementHolder &  ShowerEleHolder  )

inlinevirtualinherited

Reimplemented in ShowerRecoTools::ShowerPCADirection, ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder, ShowerRecoTools::ShowerExampleTool, and ShowerRecoTools::ShowerSkeletonTool.

Definition at line 74 of file IShowerTool.h.

    {
      return 0;
    }

template<class T , class A , class B >

void ShowerRecoTools::IShowerTool::AddSingle ( A &  a, B &  b, std::string  Name  )

inlineprotectedinherited

Definition at line 152 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::AddSingle().

    {
      UniquePtrs->AddSingle<T>(a, b, Name);
    }

int ShowerRecoTools::ShowerDirectionCheater::CalculateElement ( const art::Ptr< recob::PFParticle > &  pfparticle, art::Event &  Event, reco::shower::ShowerElementHolder &  ShowerEleHolder  )

overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 77 of file ShowerDirectionCheater_tool.cc.

References reco::shower::ShowerElementHolder::CheckElement(), DEFINE_ART_CLASS_TOOL, fLArPandoraShowerCheatingAlg, fNSegments, fPFParticleLabel, fRMSFlip, fShowerDirectionOutputLabel, fShowerStartPositionInputLabel, fTrueParticleInputLabel, fVertexFlip, reco::shower::ShowerElementHolder::GetElement(), ShowerRecoTools::IShowerTool::GetLArPandoraShowerAlg(), shower::LArPandoraShowerCheatingAlg::GetTrueChain(), shower::LArPandoraShowerCheatingAlg::GetTrueParticleMap(), art::ProductRetriever::getValidHandle(), hits(), art::Ptr< T >::key(), simb::MCParticle::Px(), simb::MCParticle::Py(), simb::MCParticle::Pz(), rmsGradient, shower::LArPandoraShowerAlg::RMSShowerGradient(), reco::shower::ShowerElementHolder::SetElement(), shower::LArPandoraShowerAlg::ShowerCentre(), Tree, shower::LArPandoraShowerCheatingAlg::TrueParticleIDFromTrueChain(), and vertexDotProduct.

  {

    const simb::MCParticle* trueParticle;

    //Get the hits from the shower:
    auto const pfpHandle = Event.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleLabel);

    auto const clockData =
      art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(Event);
    auto const detProp =
      art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(Event, clockData);

    if (ShowerEleHolder.CheckElement(fTrueParticleInputLabel)) {
      ShowerEleHolder.GetElement(fTrueParticleInputLabel, trueParticle);
    }
    else {

      //Could store these in the shower element holder and just calculate once?
      std::map<int, const simb::MCParticle*> trueParticles =
        fLArPandoraShowerCheatingAlg.GetTrueParticleMap();
      std::map<int, std::vector<int>> showersMothers =
        fLArPandoraShowerCheatingAlg.GetTrueChain(trueParticles);

      //Get the clusters
      auto const clusHandle = Event.getValidHandle<std::vector<recob::Cluster>>(fPFParticleLabel);
      art::FindManyP<recob::Cluster> fmc(pfpHandle, Event, fPFParticleLabel);
      std::vector<art::Ptr<recob::Cluster>> clusters = fmc.at(pfparticle.key());

      //Get the hit association
      art::FindManyP<recob::Hit> fmhc(clusHandle, Event, fPFParticleLabel);

      std::vector<art::Ptr<recob::Hit>> showerHits;
      for (auto const& cluster : clusters) {
        //Get the hits
        std::vector<art::Ptr<recob::Hit>> hits = fmhc.at(cluster.key());
        showerHits.insert(showerHits.end(), hits.begin(), hits.end());
      }

      //Get the true particle from the shower
      std::pair<int, double> ShowerTrackInfo =
        fLArPandoraShowerCheatingAlg.TrueParticleIDFromTrueChain(
          clockData, showersMothers, showerHits, 2);

      if (ShowerTrackInfo.first == -99999) {
        mf::LogError("ShowerDirectionCheater") << "True shower not found, returning";
        return 1;
      }
      trueParticle = trueParticles[ShowerTrackInfo.first];
      ShowerEleHolder.SetElement(trueParticle, fTrueParticleInputLabel);
    }

    if (!trueParticle) {
      mf::LogError("ShowerDirectionCheater") << "True shower not found, returning";
      return 1;
    }

    auto trueDir = geo::Vector_t{trueParticle->Px(), trueParticle->Py(), trueParticle->Pz()}.Unit();

    geo::Vector_t trueDirErr = {-999, -999, -999};
    ShowerEleHolder.SetElement(trueDir, trueDirErr, fShowerDirectionOutputLabel);

    if (fRMSFlip || fVertexFlip) {

      // Reset the tree values to defaults
      rmsGradient = std::numeric_limits<float>::lowest();
      vertexDotProduct = std::numeric_limits<float>::lowest();

      //Get the SpacePoints and hits
      art::FindManyP<recob::SpacePoint> fmspp(pfpHandle, Event, fPFParticleLabel);

      if (!fmspp.isValid()) {
        throw cet::exception("ShowerDirectionCheater")
          << "Trying to get the spacepoint and failed. Something is not configured correctly. "
             "Stopping ";
      }

      auto const spHandle = Event.getValidHandle<std::vector<recob::SpacePoint>>(fPFParticleLabel);
      art::FindManyP<recob::Hit> fmh(spHandle, Event, fPFParticleLabel);
      if (!fmh.isValid()) {
        throw cet::exception("ShowerDirectionCheater")
          << "Spacepoint and hit association not valid. Stopping.";
      }
      std::vector<art::Ptr<recob::SpacePoint>> spacePoints = fmspp.at(pfparticle.key());

      if (spacePoints.size() < 3) {
        mf::LogWarning("ShowerDirectionCheater")
          << spacePoints.size() << " spacepoints in shower, not calculating direction" << std::endl;
        return 1;
      }

      //Get Shower Centre
      float TotalCharge;

      auto const ShowerCentre = IShowerTool::GetLArPandoraShowerAlg().ShowerCentre(
        clockData, detProp, spacePoints, fmh, TotalCharge);

      //Check if we are pointing the correct direction or not, First try the start position
      if (ShowerEleHolder.CheckElement(fShowerStartPositionInputLabel) && fVertexFlip) {

        //Get the General direction as the vector between the start position and the centre
        geo::Point_t StartPositionVec = {-999, -999, -999};
        ShowerEleHolder.GetElement(fShowerStartPositionInputLabel, StartPositionVec);

        auto const GeneralDir = (ShowerCentre - StartPositionVec).Unit();

        //Dot product
        vertexDotProduct = trueDir.Dot(GeneralDir);
      }

      if (fRMSFlip) {
        //Otherwise Check against the RMS of the shower. Method adapated from EMShower Thanks Mike.
        rmsGradient = IShowerTool::GetLArPandoraShowerAlg().RMSShowerGradient(
          spacePoints, ShowerCentre, trueDir, fNSegments);
      }
      Tree->Fill();
    }

    return 0;
  }

const shower::LArPandoraShowerAlg& ShowerRecoTools::IShowerTool::GetLArPandoraShowerAlg ( ) const

inlineprotectedinherited

Definition at line 82 of file IShowerTool.h.

Referenced by ShowerRecoTools::ShowerLengthPercentile::CalculateElement(), ShowerRecoTools::ShowerUnidirectiondEdx::CalculateElement(), ShowerRecoTools::Shower3DCylinderTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPCAPropergationStartPosition::CalculateElement(), ShowerRecoTools::Shower2DLinearRegressionTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPFPVertexStartPosition::CalculateElement(), CalculateElement(), ShowerRecoTools::ShowerPCADirection::CalculateElement(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerTrajPointdEdx::CalculateElement(), ShowerRecoTools::ShowerPCADirection::CalculateShowerPCA(), ShowerRecoTools::Shower2DLinearRegressionTrackHitFinder::FindInitialTrackHits(), ShowerRecoTools::Shower3DCylinderTrackHitFinder::FindTrackSpacePoints(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::FitSegmentAndCalculateResidual(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::PruneFrontOfSPSPool(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::PruneTrack(), ShowerRecoTools::ShowerIncrementalTrackHitFinder::RunTestOfIncrementalSpacePointFinder(), ShowerRecoTools::ShowerTrackPCADirection::ShowerPCAVector(), and ShowerRecoTools::ShowerIncrementalTrackHitFinder::ShowerPCAVector().

    {
      return fLArPandoraShowerAlg;
    };

template<class T >

art::Ptr<T> ShowerRecoTools::IShowerTool::GetProducedElementPtr ( std::string  Name, reco::shower::ShowerElementHolder &  ShowerEleHolder, int  iter = -1  )

inlineprotectedinherited

Definition at line 102 of file IShowerTool.h.

References reco::shower::ShowerElementHolder::CheckElement(), reco::shower::ShowerProducedPtrsHolder::CheckUniqueProduerPtr(), reco::shower::ShowerProducedPtrsHolder::GetArtPtr(), and reco::shower::ShowerElementHolder::GetShowerNumber().

    {

      //Check the element has been set
      bool check_element = ShowerEleHolder.CheckElement(Name);
      if (!check_element) {
        throw cet::exception("IShowerTool") << "tried to get a element that does not exist. Failed "
                                               "at making the art ptr for Element: "
                                            << Name << std::endl;
      }

      //Check the unique ptr has been set.
      bool check_ptr = UniquePtrs->CheckUniqueProduerPtr(Name);
      if (!check_ptr) {
        throw cet::exception("IShowerTool")
          << "tried to get a ptr that does not exist. Failed at making the art ptr for Element"
          << Name;
      }

      //Check if the user has defined an index if not just use the current shower index/
      int index;
      if (iter != -1) { index = iter; }
      else {
        index = ShowerEleHolder.GetShowerNumber();
      }

      //Make the ptr
      return UniquePtrs->GetArtPtr<T>(Name, index);
    }

int ShowerRecoTools::IShowerTool::GetVectorPtrSize ( std::string  Name )

inlineprotectedinherited

Definition at line 158 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::GetVectorPtrSize().

Referenced by ShowerRecoTools::ShowerExampleTool::AddAssociations(), ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder::AddAssociations(), and ShowerRecoTools::ShowerPCADirection::AddAssociations().

{ return UniquePtrs->GetVectorPtrSize(Name); }

void ShowerRecoTools::IShowerTool::InitaliseProducerPtr ( reco::shower::ShowerProducedPtrsHolder &  uniqueproducerPtrs )

inlineinherited

Definition at line 68 of file IShowerTool.h.

    {
      UniquePtrs = &uniqueproducerPtrs;
    }

virtual void ShowerRecoTools::IShowerTool::InitialiseProducers ( )

inlinevirtualinherited

Reimplemented in ShowerRecoTools::ShowerPCADirection, ShowerRecoTools::ShowerPandoraSlidingFitTrackFinder, and ShowerRecoTools::ShowerExampleTool.

Definition at line 62 of file IShowerTool.h.

{}

template<class T >

void ShowerRecoTools::IShowerTool::InitialiseProduct ( std::string  Name, std::string  InstanceName = ""  )

inlineprotectedinherited

Definition at line 137 of file IShowerTool.h.

References art::ProducesCollector::produces(), and reco::shower::ShowerProducedPtrsHolder::SetShowerUniqueProduerPtr().

    {

      if (collectorPtr == nullptr) {
        mf::LogWarning("IShowerTool") << "The art::ProducesCollector ptr has not been set";
        return;
      }

      collectorPtr->produces<T>(InstanceName);
      UniquePtrs->SetShowerUniqueProduerPtr(type<T>(), Name, InstanceName);
    }

void ShowerRecoTools::IShowerTool::PrintPtr ( std::string  Name )

inlineprotectedinherited

Definition at line 162 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::PrintPtr().

Referenced by ShowerRecoTools::ShowerExampleTool::CalculateElement().

{ UniquePtrs->PrintPtr(Name); }

void ShowerRecoTools::IShowerTool::PrintPtrs ( )

inlineprotectedinherited

Definition at line 160 of file IShowerTool.h.

References reco::shower::ShowerProducedPtrsHolder::PrintPtrs().

Referenced by ShowerRecoTools::ShowerExampleTool::CalculateElement().

{ UniquePtrs->PrintPtrs(); }

int ShowerRecoTools::IShowerTool::RunShowerTool ( const art::Ptr< recob::PFParticle > &  pfparticle, art::Event &  Event, reco::shower::ShowerElementHolder &  ShowerEleHolder, std::string  evd_display_name_append = ""  )

inlineinherited

Definition at line 46 of file IShowerTool.h.

    {

      int calculation_status = CalculateElement(pfparticle, Event, ShowerEleHolder);
      if (calculation_status != 0) return calculation_status;
      if (fRunEventDisplay) {
        IShowerTool::GetLArPandoraShowerAlg().DebugEVD(
          pfparticle, Event, ShowerEleHolder, evd_display_name_append);
      }
      return calculation_status;
    }

void ShowerRecoTools::IShowerTool::SetPtr ( art::ProducesCollector *  collector )

inlineinherited

Definition at line 65 of file IShowerTool.h.

{ collectorPtr = collector; }

Member Data Documentation

shower::LArPandoraShowerCheatingAlg ShowerRecoTools::ShowerDirectionCheater::fLArPandoraShowerCheatingAlg

private

Definition at line 36 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement().

const unsigned int ShowerRecoTools::ShowerDirectionCheater::fNSegments

private

Definition at line 44 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement().

const art::InputTag ShowerRecoTools::ShowerDirectionCheater::fPFParticleLabel

private

Definition at line 42 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement().

const bool ShowerRecoTools::ShowerDirectionCheater::fRMSFlip

private

Definition at line 45 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement(), and ShowerDirectionCheater().

const std::string ShowerRecoTools::ShowerDirectionCheater::fShowerDirectionOutputLabel

private

Definition at line 56 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement().

const std::string ShowerRecoTools::ShowerDirectionCheater::fShowerStartPositionInputLabel

private

Definition at line 54 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement().

const std::string ShowerRecoTools::ShowerDirectionCheater::fTrueParticleInputLabel

private

Definition at line 55 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement().

const bool ShowerRecoTools::ShowerDirectionCheater::fVertexFlip

private

Definition at line 47 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement(), and ShowerDirectionCheater().

float ShowerRecoTools::ShowerDirectionCheater::rmsGradient

private

Definition at line 52 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement(), and ShowerDirectionCheater().

art::ServiceHandle<art::TFileService> ShowerRecoTools::ShowerDirectionCheater::tfs

private

Definition at line 39 of file ShowerDirectionCheater_tool.cc.

Referenced by ShowerDirectionCheater().

TTree* ShowerRecoTools::ShowerDirectionCheater::Tree

private

Definition at line 50 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement(), and ShowerDirectionCheater().

float ShowerRecoTools::ShowerDirectionCheater::vertexDotProduct

private

Definition at line 51 of file ShowerDirectionCheater_tool.cc.

Referenced by CalculateElement(), and ShowerDirectionCheater().

---

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

• larpandora/v09_22_05/source/larpandora/LArPandoraEventBuilding/LArPandoraShower/Tools/Cheating/ShowerDirectionCheater_tool.cc