#include <PiZeroProcess.h>

Public Member Functions
	PiZeroProcess (const simb::MCParticle &mcp, const art::Event &evt, const detinfo::DetectorClocksData &clockData, const detinfo::DetectorPropertiesData &detProp, std::string showerLabel)

	PiZeroProcess (const recob::Shower &shower, const art::Event &evt, const detinfo::DetectorClocksData &clockData, const detinfo::DetectorPropertiesData &detProp, std::string showerLabel)

	~PiZeroProcess ()

void	reset ()

const simb::MCParticle *	pi0 () const

const ShowerProcess *	showerProcess1 () const

const ShowerProcess *	showerProcess2 () const

const simb::MCParticle *	photon1 () const

const simb::MCParticle *	photon2 () const

const recob::Shower *	shower1 () const

const recob::Shower *	shower2 () const

const recob::Track *	track1 () const

const recob::Track *	track2 () const

const Cone *	cone1 () const

const Cone *	cone2 () const

const art::Event *	evt () const

std::string	showerLabel () const

bool	allMCSet () const

bool	allRecoSet () const

void	setPi0 (const simb::MCParticle &newPi0)

void	createBranches (TTree *tree)

void	fillMCInfo ()

void	fillRecoInfo (calo::CalorimetryAlg caloAlg)

Public Attributes
bool	_haveMCInfo

bool	_haveRecoInfo

Private Attributes
std::unique_ptr< ShowerProcess >	_shProcess1

std::unique_ptr< ShowerProcess >	_shProcess2

const simb::MCParticle *	_pi0 = 0x0

const art::Event &	_evt

const detinfo::DetectorClocksData &	_clockData

const detinfo::DetectorPropertiesData &	_detProp

const std::string	_showerLabel

art::ServiceHandle< cheat::ParticleInventoryService >	pi_serv

protoana::ProtoDUNETruthUtils	truthUtils

Detailed Description

Definition at line 95 of file PiZeroProcess.h.

Constructor & Destructor Documentation

pizero::PiZeroProcess::PiZeroProcess	(	const simb::MCParticle &	mcp,
		const art::Event &	evt,
		const detinfo::DetectorClocksData &	clockData,
		const detinfo::DetectorPropertiesData &	detProp,
		std::string	showerLabel
	)

Definition at line 191 of file PiZeroProcess.h.

                                                    :
   _evt(evt), _clockData(clockData), _detProp{detProp}, _showerLabel(showerLabel)
 {
   if(mcp.PdgCode() == 111) {
     // Got a pi0.
     setPi0(mcp);
   } else if(mcp.PdgCode() == 22) {
     // Got a photon. Find the mother and run through the above process again.
     if(mcp.Mother() == 0) return;
     const simb::MCParticle* mcpmom = pi_serv->TrackIdToParticle_P(mcp.Mother());
     if(mcpmom->PdgCode() == 111) setPi0(*mcpmom);
   }
 
   // Quit if the system did not manage to find appropriate MCParticles.
   if(!allMCSet()) return;
 }

pizero::PiZeroProcess::PiZeroProcess	(	const recob::Shower &	shower,
		const art::Event &	evt,
		const detinfo::DetectorClocksData &	clockData,
		const detinfo::DetectorPropertiesData &	detProp,
		std::string	showerLabel
	)

Definition at line 212 of file PiZeroProcess.h.

                                                    :
   _evt(evt), _clockData(clockData), _detProp{detProp}, _showerLabel(showerLabel)
 {
   // Get the event's showers.
   auto showerHandle = _evt.getHandle<std::vector<recob::Shower> >(_showerLabel);
   if (!showerHandle) return;
   if(showerHandle->size() < 2) return;
 
   // For each shower, find the one which is closest to intersecting.
   double min_dist = 100000;
   const recob::Shower* shMatch = 0x0;
   for(const recob::Shower& potMatch : *showerHandle) {
     if(&potMatch == &shower) continue;
 
     const double this_dist = pizero::ClosestDistance(&shower, &potMatch);
     if(this_dist < min_dist) {
       min_dist = this_dist;
       shMatch = &potMatch;
     }
   }
   TVector3 reco_pi0_location = pizero::ClosestPoint(&shower, shMatch);
 
   // Find closest MC pi0 match to the reconstructed pi0 vertex.
   double reco_mismatch = 100000;
   const simb::MCParticle* min_pi0 = 0;
   for(auto const& p : pi_serv->ParticleList()) {
     const simb::MCParticle* true_pi0 = p.second;
     // std::cout << "Particle PDG: " << true_pi0->PdgCode() << '\n';
     if(true_pi0->PdgCode() != 111) continue;
 
     const double this_dist = (true_pi0->Position().Vect() - reco_pi0_location).Mag();
     if(this_dist < reco_mismatch) {
       reco_mismatch = this_dist;
       min_pi0 = true_pi0;
     }
   }
   if(min_pi0 == 0) return;
 
   // Fill internal data fields with the found pointers.
   setPi0(*min_pi0);
   // Check whether we have all necessary pointers.
   if(!allMCSet()) {
     _pi0 = 0x0;
     _shProcess1 = 0x0;
     _shProcess2 = 0x0;
     // _photon1 = 0x0;
     // _photon2 = 0x0;
     // _shower1 = 0x0;
     // _shower2 = 0x0;
     return;
   }
   // // Use truth utilities to see which shower matches the first photon best.
   // for(auto const& p : truthUtils.GetRecoShowerListFromMCParticle(*_photon1, _evt, _showerLabel)) {
   //   if(p.first->Length() - shower.Length() < 1e-5) {
   //     _shower1 = &shower;
   //     _shower2 = shMatch;
   //     break;
   //   } else if(p.first->Length() - shMatch->Length() < 1e-5) {
   //     _shower1 = shMatch;
   //     _shower2 = &shower;
   //     break;
   //   }
   // }
 }

pizero::PiZeroProcess::~PiZeroProcess ( )

inline

Definition at line 128 of file PiZeroProcess.h.

128 {};

Member Function Documentation

bool pizero::PiZeroProcess::allMCSet ( ) const

Definition at line 281 of file PiZeroProcess.h.

                                    {
   return (pi0() != 0x0) && (photon1() != 0x0) && (photon2() != 0x0);
 }

bool pizero::PiZeroProcess::allRecoSet ( ) const

Definition at line 284 of file PiZeroProcess.h.

                                      {
   return (shower1() != 0x0) && (shower2() != 0x0);
 }

const Cone* pizero::PiZeroProcess::cone1 ( ) const

inline

Definition at line 162 of file PiZeroProcess.h.

                             {
     return _shProcess1? _shProcess1->cone(): 0x0;
   }

const Cone* pizero::PiZeroProcess::cone2 ( ) const

inline

Definition at line 165 of file PiZeroProcess.h.

                             {
     return _shProcess2? _shProcess2->cone(): 0x0;
   }

void pizero::PiZeroProcess::createBranches ( TTree * tree )

const art::Event* pizero::PiZeroProcess::evt ( ) const

inline

Definition at line 168 of file PiZeroProcess.h.

168 { return &_evt; }

pizero::PiZeroProcess::_evt

const art::Event & _evt

Definition: PiZeroProcess.h:107

void pizero::PiZeroProcess::fillMCInfo ( )

void pizero::PiZeroProcess::fillRecoInfo ( calo::CalorimetryAlg caloAlg )

const simb::MCParticle* pizero::PiZeroProcess::photon1 ( ) const

inline

Definition at line 144 of file PiZeroProcess.h.

                                         {
     return _shProcess1? _shProcess1->mcparticle(): 0x0;
   }

const simb::MCParticle* pizero::PiZeroProcess::photon2 ( ) const

inline

Definition at line 147 of file PiZeroProcess.h.

                                         {
     return _shProcess2? _shProcess2->mcparticle(): 0x0;
   }

const simb::MCParticle* pizero::PiZeroProcess::pi0 ( ) const

inline

Definition at line 137 of file PiZeroProcess.h.

137 { return _pi0; }

pizero::PiZeroProcess::_pi0

const simb::MCParticle * _pi0

Definition: PiZeroProcess.h:103

void pizero::PiZeroProcess::reset ( )

void pizero::PiZeroProcess::setPi0 ( const simb::MCParticle & newPi0 )

Definition at line 289 of file PiZeroProcess.h.

                                                        {
   _pi0 = &newPi0;
   if (_pi0->NumberDaughters() == 2) { // Skip rare decays for now.
     // Find the associated photons through particle inventory.
     const simb::MCParticle* photon1 = pi_serv->TrackIdToParticle_P(_pi0->Daughter(0));
     const simb::MCParticle* photon2 = pi_serv->TrackIdToParticle_P(_pi0->Daughter(1));
     // Set photon1 to be the more energetic one.
     if(photon1->E() < photon2->E()) std::swap(photon1, photon2);
     // Fill the showerProcesses.
     _shProcess1 = std::make_unique<ShowerProcess>(*photon1, _evt, _clockData, _detProp);
     _shProcess2 = std::make_unique<ShowerProcess>(*photon2, _evt, _clockData, _detProp);
   }
 }

const recob::Shower* pizero::PiZeroProcess::shower1 ( ) const

inline

Definition at line 150 of file PiZeroProcess.h.

                                      {
     return _shProcess1? _shProcess1->shower(): 0x0;
   }

const recob::Shower* pizero::PiZeroProcess::shower2 ( ) const

inline

Definition at line 153 of file PiZeroProcess.h.

                                      {
     return _shProcess2? _shProcess2->shower(): 0x0;
   }

std::string pizero::PiZeroProcess::showerLabel ( ) const

inline

Definition at line 169 of file PiZeroProcess.h.

169 { return _showerLabel; }

pizero::PiZeroProcess::_showerLabel

const std::string _showerLabel

Definition: PiZeroProcess.h:110

const ShowerProcess* pizero::PiZeroProcess::showerProcess1 ( ) const

inline

Definition at line 138 of file PiZeroProcess.h.

                                               {
     return &*_shProcess1;
   }

const ShowerProcess* pizero::PiZeroProcess::showerProcess2 ( ) const

inline

Definition at line 141 of file PiZeroProcess.h.

                                               {
     return &*_shProcess2;
   }

const recob::Track* pizero::PiZeroProcess::track1 ( ) const

inline

Definition at line 156 of file PiZeroProcess.h.

                                    {
     return _shProcess1? _shProcess1->track(): 0x0;
   }

const recob::Track* pizero::PiZeroProcess::track2 ( ) const

inline

Definition at line 159 of file PiZeroProcess.h.

                                    {
     return _shProcess2? _shProcess2->track(): 0x0;
   }

Member Data Documentation

const detinfo::DetectorClocksData& pizero::PiZeroProcess::_clockData

private

Definition at line 108 of file PiZeroProcess.h.

const detinfo::DetectorPropertiesData& pizero::PiZeroProcess::_detProp

private

Definition at line 109 of file PiZeroProcess.h.

const art::Event& pizero::PiZeroProcess::_evt

private

Definition at line 107 of file PiZeroProcess.h.

bool pizero::PiZeroProcess::_haveMCInfo

Definition at line 128 of file PiZeroProcess.h.

bool pizero::PiZeroProcess::_haveRecoInfo

Definition at line 131 of file PiZeroProcess.h.

const simb::MCParticle* pizero::PiZeroProcess::_pi0 = 0x0

private

Definition at line 103 of file PiZeroProcess.h.

const std::string pizero::PiZeroProcess::_showerLabel

private

Definition at line 110 of file PiZeroProcess.h.

std::unique_ptr<ShowerProcess> pizero::PiZeroProcess::_shProcess1

private

Definition at line 98 of file PiZeroProcess.h.

std::unique_ptr<ShowerProcess> pizero::PiZeroProcess::_shProcess2

private

Definition at line 99 of file PiZeroProcess.h.

art::ServiceHandle<cheat::ParticleInventoryService> pizero::PiZeroProcess::pi_serv

private

Definition at line 113 of file PiZeroProcess.h.

protoana::ProtoDUNETruthUtils pizero::PiZeroProcess::truthUtils

private

Definition at line 114 of file PiZeroProcess.h.

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

protoduneana/protoduneana/singlephase/EMTaskForce/PiZeroProcess.h

Public Member Functions

Public Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation