LArSoft  v10_04_05
Liquid Argon Software toolkit - https://larsoft.org/
ShowerRecoTools::ShowerStartPositionCheater Class Reference
Inheritance diagram for ShowerRecoTools::ShowerStartPositionCheater:
ShowerRecoTools::IShowerTool

Public Member Functions

 ShowerStartPositionCheater (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::LArPandoraShowerAlgGetLArPandoraShowerAlg () 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
 
const art::InputTag fPFParticleLabel
 
const art::InputTag fHitModuleLabel
 
const std::string fShowerStartPositionOutputLabel
 
const std::string fTrueParticleOutputLabel
 

Detailed Description

Definition at line 21 of file ShowerStartPositionCheater_tool.cc.

Constructor & Destructor Documentation

ShowerRecoTools::ShowerStartPositionCheater::ShowerStartPositionCheater ( const fhicl::ParameterSet pset)

Definition at line 43 of file ShowerStartPositionCheater_tool.cc.

44  : IShowerTool(pset.get<fhicl::ParameterSet>("BaseTools"))
45  , fLArPandoraShowerCheatingAlg(pset.get<fhicl::ParameterSet>("LArPandoraShowerCheatingAlg"))
46  , fPFParticleLabel(pset.get<art::InputTag>("PFParticleLabel"))
47  , fHitModuleLabel(pset.get<art::InputTag>("HitModuleLabel"))
48  , fShowerStartPositionOutputLabel(pset.get<std::string>("ShowerStartPositionOutputLabel"))
49  , fTrueParticleOutputLabel(pset.get<std::string>("TrueParticleOutputLabel"))
50  {}
T get(std::string const &key) const
Definition: ParameterSet.h:314
shower::LArPandoraShowerCheatingAlg fLArPandoraShowerCheatingAlg
IShowerTool(const fhicl::ParameterSet &pset)
Definition: IShowerTool.h:31

Member Function Documentation

virtual int ShowerRecoTools::IShowerTool::AddAssociations ( const art::Ptr< recob::PFParticle > &  pfpPtr,
art::Event Event,
reco::shower::ShowerElementHolder ShowerEleHolder 
)
inlinevirtualinherited
template<class T , class A , class B >
void ShowerRecoTools::IShowerTool::AddSingle ( A &  a,
B &  b,
std::string  Name 
)
inlineprotectedinherited

Definition at line 150 of file IShowerTool.h.

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

151  {
152  UniquePtrs->AddSingle<T>(a, b, Name);
153  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:83
void AddSingle(A &a, B &b, const std::string &Name)
int ShowerRecoTools::ShowerStartPositionCheater::CalculateElement ( const art::Ptr< recob::PFParticle > &  pfparticle,
art::Event Event,
reco::shower::ShowerElementHolder ShowerEleHolder 
)
overridevirtual

Implements ShowerRecoTools::IShowerTool.

Definition at line 52 of file ShowerStartPositionCheater_tool.cc.

References util::abs(), DEFINE_ART_CLASS_TOOL, simb::MCParticle::E(), fLArPandoraShowerCheatingAlg, fPFParticleLabel, fShowerStartPositionOutputLabel, fTrueParticleOutputLabel, shower::LArPandoraShowerCheatingAlg::GetTrueChain(), shower::LArPandoraShowerCheatingAlg::GetTrueParticleMap(), art::ProductRetriever::getValidHandle(), hits(), art::Ptr< T >::key(), simb::MCParticle::NumberTrajectoryPoints(), simb::MCParticle::PdgCode(), simb::MCParticle::Position(), reco::shower::ShowerElementHolder::SetElement(), geo::vect::toPoint(), and shower::LArPandoraShowerCheatingAlg::TrueParticleIDFromTrueChain().

56  {
57 
58  //Could store these in the shower element holder and just calculate once?
59  std::map<int, const simb::MCParticle*> trueParticles =
61  std::map<int, std::vector<int>> showersMothers =
63 
64  //Get the hits from the shower:
65  auto const pfpHandle = Event.getValidHandle<std::vector<recob::PFParticle>>(fPFParticleLabel);
66 
67  //Get the clusters
68  auto const clusHandle = Event.getValidHandle<std::vector<recob::Cluster>>(fPFParticleLabel);
69 
70  art::FindManyP<recob::Cluster> fmc(pfpHandle, Event, fPFParticleLabel);
71  std::vector<art::Ptr<recob::Cluster>> clusters = fmc.at(pfparticle.key());
72 
73  //Get the hit association
74  art::FindManyP<recob::Hit> fmhc(clusHandle, Event, fPFParticleLabel);
75 
76  std::vector<art::Ptr<recob::Hit>> showerHits;
77  for (auto const& cluster : clusters) {
78 
79  //Get the hits
80  std::vector<art::Ptr<recob::Hit>> hits = fmhc.at(cluster.key());
81  showerHits.insert(showerHits.end(), hits.begin(), hits.end());
82  }
83 
84  //Get the true particle from the shower
85  auto const clockData =
87 
88  std::pair<int, double> ShowerTrackInfo =
90  clockData, showersMothers, showerHits, 2);
91 
92  if (ShowerTrackInfo.first == -99999) {
93  mf::LogError("ShowerStartPositionCheater") << "True Shower Not Found";
94  return 1;
95  }
96 
97  const simb::MCParticle* trueParticle = trueParticles[ShowerTrackInfo.first];
98  if (!trueParticle) {
99  mf::LogError("ShowerDirectionCheater") << "True shower not found, returning";
100  return 1;
101  }
102 
103  ShowerEleHolder.SetElement(trueParticle, fTrueParticleOutputLabel);
104 
105  geo::Point_t trueStartPos{-999, -999, -999};
106  // If the true particle is a photon, we need to be smarter.
107  // Select the first traj point where tne photon loses energy
108  if (abs(trueParticle->PdgCode()) == 22) {
109  double initialEnergy = trueParticle->E();
110  unsigned int TrajPoints = trueParticle->NumberTrajectoryPoints();
111  for (unsigned int trajPoint = 0; trajPoint < TrajPoints; trajPoint++) {
112  if (trueParticle->E(trajPoint) < initialEnergy) {
113  trueStartPos = geo::vect::toPoint(trueParticle->Position(trajPoint).Vect());
114  break;
115  }
116  }
117  }
118  else {
119  trueStartPos = geo::vect::toPoint(trueParticle->Position().Vect());
120  }
121 
122  geo::Point_t trueStartPosErr = {-999, -999, -999};
123  ShowerEleHolder.SetElement(trueStartPos, trueStartPosErr, fShowerStartPositionOutputLabel);
124 
125  return 0;
126  }
double E(const int i=0) const
Definition: MCParticle.h:234
unsigned int NumberTrajectoryPoints() const
Definition: MCParticle.h:219
const TLorentzVector & Position(const int i=0) const
Definition: MCParticle.h:220
int PdgCode() const
Definition: MCParticle.h:213
std::map< int, const simb::MCParticle * > GetTrueParticleMap() const
void SetElement(T &dataproduct, const std::string &Name, bool checktag=false)
constexpr auto abs(T v)
Returns the absolute value of the argument.
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Cluster finding and building.
void hits()
Definition: readHits.C:15
key_type key() const noexcept
Definition: Ptr.h:166
std::map< int, std::vector< int > > GetTrueChain(std::map< int, const simb::MCParticle * > &trueParticles) const
shower::LArPandoraShowerCheatingAlg fLArPandoraShowerCheatingAlg
Point_t toPoint(Point const &p)
Convert the specified point into a geo::Point_t.
std::pair< int, double > TrueParticleIDFromTrueChain(detinfo::DetectorClocksData const &clockData, std::map< int, std::vector< int >> const &ShowersMothers, std::vector< art::Ptr< recob::Hit >> const &hits, int planeid) const
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Point_t
Type for representation of position in physical 3D space.
Definition: geo_vectors.h:180
ValidHandle< PROD > getValidHandle(InputTag const &tag) const
const shower::LArPandoraShowerAlg& ShowerRecoTools::IShowerTool::GetLArPandoraShowerAlg ( ) const
inlineprotectedinherited

Definition at line 80 of file IShowerTool.h.

Referenced by ShowerRecoTools::ShowerLengthPercentile::CalculateElement(), ShowerRecoTools::Shower3DCylinderTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPCAPropergationStartPosition::CalculateElement(), ShowerRecoTools::Shower2DLinearRegressionTrackHitFinder::CalculateElement(), ShowerRecoTools::ShowerPFPVertexStartPosition::CalculateElement(), ShowerRecoTools::ShowerUnidirectiondEdx::CalculateElement(), ShowerRecoTools::ShowerDirectionCheater::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().

81  {
82  return fLArPandoraShowerAlg;
83  };
shower::LArPandoraShowerAlg fLArPandoraShowerAlg
Definition: IShowerTool.h:90
template<class T >
art::Ptr<T> ShowerRecoTools::IShowerTool::GetProducedElementPtr ( std::string  Name,
reco::shower::ShowerElementHolder ShowerEleHolder,
int  iter = -1 
)
inlineprotectedinherited

Definition at line 100 of file IShowerTool.h.

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

103  {
104 
105  //Check the element has been set
106  bool check_element = ShowerEleHolder.CheckElement(Name);
107  if (!check_element) {
108  throw cet::exception("IShowerTool") << "tried to get a element that does not exist. Failed "
109  "at making the art ptr for Element: "
110  << Name << std::endl;
111  }
112 
113  //Check the unique ptr has been set.
114  bool check_ptr = UniquePtrs->CheckUniqueProduerPtr(Name);
115  if (!check_ptr) {
116  throw cet::exception("IShowerTool")
117  << "tried to get a ptr that does not exist. Failed at making the art ptr for Element"
118  << Name;
119  }
120 
121  //Check if the user has defined an index if not just use the current shower index/
122  int index;
123  if (iter != -1) { index = iter; }
124  else {
125  index = ShowerEleHolder.GetShowerNumber();
126  }
127 
128  //Make the ptr
129  return UniquePtrs->GetArtPtr<T>(Name, index);
130  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:83
bool CheckUniqueProduerPtr(const std::string &Name) const
bool CheckElement(const std::string &Name) const
art::Ptr< T > GetArtPtr(const std::string &Name, const int &iter) const
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
int ShowerRecoTools::IShowerTool::GetVectorPtrSize ( std::string  Name)
inlineprotectedinherited
void ShowerRecoTools::IShowerTool::InitaliseProducerPtr ( reco::shower::ShowerProducedPtrsHolder uniqueproducerPtrs)
inlineinherited

Definition at line 66 of file IShowerTool.h.

67  {
68  UniquePtrs = &uniqueproducerPtrs;
69  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:83
virtual void ShowerRecoTools::IShowerTool::InitialiseProducers ( )
inlinevirtualinherited
template<class T >
void ShowerRecoTools::IShowerTool::InitialiseProduct ( std::string  Name,
std::string  InstanceName = "" 
)
inlineprotectedinherited

Definition at line 135 of file IShowerTool.h.

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

136  {
137 
138  if (collectorPtr == nullptr) {
139  mf::LogWarning("IShowerTool") << "The art::ProducesCollector ptr has not been set";
140  return;
141  }
142 
143  collectorPtr->produces<T>(InstanceName);
144  UniquePtrs->SetShowerUniqueProduerPtr(type<T>(), Name, InstanceName);
145  }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:83
art::ProducesCollector * collectorPtr
Definition: IShowerTool.h:95
void produces(std::string const &instanceName={}, Persistable const persistable=Persistable::Yes)
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
int SetShowerUniqueProduerPtr(type< T >, const std::string &Name, const std::string &Instance="")
void ShowerRecoTools::IShowerTool::PrintPtr ( std::string  Name)
inlineprotectedinherited

Definition at line 160 of file IShowerTool.h.

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

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

160 { UniquePtrs->PrintPtr(Name); }
reco::shower::ShowerProducedPtrsHolder * UniquePtrs
Definition: IShowerTool.h:83
void PrintPtr(const std::string &Name) const
void ShowerRecoTools::IShowerTool::PrintPtrs ( )
inlineprotectedinherited
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 44 of file IShowerTool.h.

48  {
49 
50  int calculation_status = CalculateElement(pfparticle, Event, ShowerEleHolder);
51  if (calculation_status != 0) return calculation_status;
52  if (fRunEventDisplay) {
54  pfparticle, Event, ShowerEleHolder, evd_display_name_append);
55  }
56  return calculation_status;
57  }
void DebugEVD(art::Ptr< recob::PFParticle > const &pfparticle, art::Event const &Event, const reco::shower::ShowerElementHolder &ShowerEleHolder, std::string const &evd_disp_name_append="") const
virtual int CalculateElement(const art::Ptr< recob::PFParticle > &pfparticle, art::Event &Event, reco::shower::ShowerElementHolder &ShowerEleHolder)=0
const shower::LArPandoraShowerAlg & GetLArPandoraShowerAlg() const
Definition: IShowerTool.h:80
void ShowerRecoTools::IShowerTool::SetPtr ( art::ProducesCollector collector)
inlineinherited

Definition at line 63 of file IShowerTool.h.

63 { collectorPtr = collector; }
art::ProducesCollector * collectorPtr
Definition: IShowerTool.h:95

Member Data Documentation

const art::InputTag ShowerRecoTools::ShowerStartPositionCheater::fHitModuleLabel
private

Definition at line 37 of file ShowerStartPositionCheater_tool.cc.

shower::LArPandoraShowerCheatingAlg ShowerRecoTools::ShowerStartPositionCheater::fLArPandoraShowerCheatingAlg
private

Definition at line 33 of file ShowerStartPositionCheater_tool.cc.

Referenced by CalculateElement().

const art::InputTag ShowerRecoTools::ShowerStartPositionCheater::fPFParticleLabel
private

Definition at line 36 of file ShowerStartPositionCheater_tool.cc.

Referenced by CalculateElement().

const std::string ShowerRecoTools::ShowerStartPositionCheater::fShowerStartPositionOutputLabel
private

Definition at line 39 of file ShowerStartPositionCheater_tool.cc.

Referenced by CalculateElement().

const std::string ShowerRecoTools::ShowerStartPositionCheater::fTrueParticleOutputLabel
private

Definition at line 40 of file ShowerStartPositionCheater_tool.cc.

Referenced by CalculateElement().


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