showerAna::ShowerAnalysis Class Reference

Inheritance diagram for showerAna::ShowerAnalysis:

Public Member Functions
	ShowerAnalysis (const fhicl::ParameterSet &pset)
	~ShowerAnalysis ()
void	analyze (const art::Event &evt)
void	MakeDataProducts ()
void	FillData (const std::map< int, std::shared_ptr< ShowerParticle > > &particles)
void	FillPi0Data (const std::map< int, std::shared_ptr< ShowerParticle > > &particles, const std::vector< int > &pi0Decays)
int	FindTrueParticle (detinfo::DetectorClocksData const &clockData, const std::vector< art::Ptr< recob::Hit > > &hits)
int	FindParticleID (detinfo::DetectorClocksData const &clockData, const art::Ptr< recob::Hit > &hit)
std::vector< art::Ptr< recob::Hit > >	FindTrueHits (detinfo::DetectorClocksData const &clockData, const std::vector< art::Ptr< recob::Hit > > &hits, int trueParticle)
Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDAnalyzer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept
	Analyzer (fhicl::ParameterSet const &pset)
template<typename Config >
	Analyzer (Table< Config > const &config)
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Observer
	~Observer () noexcept
	Observer (Observer const &)=delete
	Observer (Observer &&)=delete
Observer &	operator= (Observer const &)=delete
Observer &	operator= (Observer &&)=delete
void	registerProducts (ProductDescriptions &, ModuleDescription const &)
void	fillDescriptions (ModuleDescription const &)
fhicl::ParameterSetID	selectorConfig () const
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
std::string	fShowerModuleLabel
std::string	fClusterModuleLabel
std::string	fHitsModuleLabel
double	fElectrondEdx
double	fPhotondEdx
double	fElectrondEdxWidth
double	fPhotondEdxWidth
art::ServiceHandle< cheat::BackTrackerService >	bt_serv
art::ServiceHandle< cheat::ParticleInventoryService >	pi_serv
art::ServiceHandle< art::TFileService >	tfs
art::ServiceHandle< geo::Geometry >	geom
TTree *	fTree
TH1D *	hClusterCompleteness
TH1D *	hLargestClusterCompleteness
TH1D *	hClusterPurity
TH1D *	hLargestClusterPurity
TH2D *	hClusterCompletenessEnergy
TH2D *	hLargestClusterCompletenessEnergy
TH2D *	hClusterCompletenessDirection
TH2D *	hLargestClusterCompletenessDirection
TH1D *	hShowerCompleteness
TH1D *	hLargestShowerCompleteness
TH1D *	hShowerPurity
TH1D *	hLargestShowerPurity
TH2D *	hShowerCompletenessEnergy
TH2D *	hLargestShowerCompletenessEnergy
TH2D *	hShowerCompletenessDirection
TH2D *	hLargestShowerCompletenessDirection
TH1D *	hShowerEnergy
TH1D *	hShowerDepositedEnergy
TH1D *	hShowerDirection
TH1D *	hShowerdEdx
TH1D *	hShowerStart
TH1D *	hShowerReconstructed
TH1D *	hNumShowersReconstructed
TH1D *	hNumShowersReconstructedNonZero
TH2D *	hShowerdEdxEnergy
TH2D *	hNumShowersReconstructedEnergy
TProfile *	hShowerReconstructedEnergy
TProfile *	hShowerdEdxEnergyProfile
TProfile *	hNumShowersReconstructedEnergyProfile
TProfile *	hShowerEnergyCompleteness
TProfile *	hShowerStartEnergy
TH1D *	hElectronPull
TH1D *	hPhotonPull
TH1D *	hPi0MassPeakReconEnergyReconAngle
TH1D *	hPi0MassPeakReconEnergyTrueAngle
TH1D *	hPi0MassPeakTrueEnergyReconAngle
TH1D *	hPi0MassPeakTrueEnergyTrueAngle

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >
using	ModuleType = EDAnalyzer
Protected Member Functions inherited from art::Observer
std::string const &	processName () const
bool	wantAllEvents () const noexcept
bool	wantEvent (ScheduleID id, Event const &e) const
Handle< TriggerResults >	getTriggerResults (Event const &e) const
	Observer (fhicl::ParameterSet const &config)
	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 194 of file ShowerAnalysis_module.cc.

Constructor & Destructor Documentation

showerAna::ShowerAnalysis::ShowerAnalysis

(

const fhicl::ParameterSet &

pset

)

Definition at line 243 of file ShowerAnalysis_module.cc.

                                                                     : EDAnalyzer(pset) {
  fShowerModuleLabel  = pset.get<std::string>("ShowerModuleLabel");
  fClusterModuleLabel = pset.get<std::string>("ClusterModuleLabel");
  fHitsModuleLabel    = pset.get<std::string>("HitsModuleLabel");
  fElectrondEdx       = pset.get<double>("ElectrondEdx",2.1);
  fPhotondEdx         = pset.get<double>("PhotondEdx",4.2);
  fElectrondEdxWidth  = pset.get<double>("ElectrondEdxWidth",0.410);
  fPhotondEdxWidth    = pset.get<double>("PhotondEdxWidth",1.217);
  this->MakeDataProducts();
}

showerAna::ShowerAnalysis::~ShowerAnalysis

(

)

Definition at line 254 of file ShowerAnalysis_module.cc.

                                         {
}

Member Function Documentation

void showerAna::ShowerAnalysis::analyze

(

const art::Event &

evt

)

Definition at line 257 of file ShowerAnalysis_module.cc.

                                                         {

  // Get showers out of event
  std::vector<art::Ptr<recob::Shower> > showers;
  auto showerHandle = evt.getHandle<std::vector<recob::Shower> >(fShowerModuleLabel);
  if (showerHandle)
    art::fill_ptr_vector(showers, showerHandle);

  // Get clusters out of event
  std::vector<art::Ptr<recob::Cluster> > clusters;
  auto clusterHandle = evt.getHandle<std::vector<recob::Cluster> >(fClusterModuleLabel);
  if (clusterHandle)
    art::fill_ptr_vector(clusters, clusterHandle);

  // Get hits out of event
  std::vector<art::Ptr<recob::Hit> > hits;
  auto hitHandle = evt.getHandle<std::vector<recob::Hit> >(fHitsModuleLabel);
  if (hitHandle)
    art::fill_ptr_vector(hits, hitHandle);

  // Get space points out of event
  std::vector<art::Ptr<recob::SpacePoint> > spacePoints;
  auto spacePointHandle = evt.getHandle<std::vector<recob::SpacePoint> >(fShowerModuleLabel);
  if (spacePointHandle)
    art::fill_ptr_vector(spacePoints, spacePointHandle);

  // Get associations out of event
  art::FindManyP<recob::Hit> fmhc(clusterHandle, evt, fClusterModuleLabel);
  art::FindManyP<recob::Hit> fmhs(showerHandle, evt, fShowerModuleLabel);
  art::FindManyP<recob::Hit> fmhsp(spacePointHandle, evt, fShowerModuleLabel);

  // Map all the true and reconstructed information for each particle
  std::map<int,std::shared_ptr<ShowerParticle> > particles;

  // Keep an eye out for pi0s!
  bool isPi0 = false;
  std::vector<int> pi0Decays;

  // Fill true properties
  const sim::ParticleList& trueParticles = pi_serv->ParticleList();
  for (sim::ParticleList::const_iterator particleIt = trueParticles.begin(); particleIt != trueParticles.end(); ++particleIt) {

    const simb::MCParticle* trueParticle = particleIt->second;
    int mother = trueParticle->Mother();
    if (mother != 0 and trueParticles.at(mother)->PdgCode() == 111) {
      isPi0 = true;
      pi0Decays.push_back(particleIt->first);
    }

    std::map<int,double> depositedEnergy;
    const std::vector<art::Ptr< sim::SimChannel >>& simChannels = bt_serv->SimChannels();
    for (std::vector<art::Ptr< sim::SimChannel >>::const_iterator channelIt = simChannels.begin(); channelIt != simChannels.end(); ++channelIt) {
      int plane = geom->View((*channelIt)->Channel());
      auto const & tdcidemap = (*channelIt)->TDCIDEMap();
      for (auto const& tdcIt : tdcidemap) {
        auto const& idevec = tdcIt.second;
        for (auto const& ideIt : idevec) {
          if (TMath::Abs(ideIt.trackID) != trueParticle->TrackId())
            continue;
          depositedEnergy[plane] += ideIt.energy / 1000;
        }
      }
    }

    std::shared_ptr<ShowerParticle> particle = std::make_shared<ShowerParticle>(trueParticle->TrackId());
    particle->SetEnergy(trueParticle->E());
    particle->SetDepositedEnergy(depositedEnergy);
    particle->SetDirection(trueParticle->Momentum().Vect().Unit());
    particle->SetStart(trueParticle->Position().Vect());
    particle->SetEnd(trueParticle->EndPosition().Vect());
    particle->SetPDG(trueParticle->PdgCode());

    particles[particleIt->first] = std::move(particle);
  }

  // Fill recon properties
  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
  for (std::vector<art::Ptr<recob::Hit> >::iterator hitIt = hits.begin(); hitIt != hits.end(); ++hitIt) {
    int trueParticle = FindParticleID(clockData, *hitIt);
    if (particles.count(trueParticle))
      particles[trueParticle]->AddAssociatedHit(*hitIt);
  }

  for (std::vector<art::Ptr<recob::Cluster> >::iterator clusterIt = clusters.begin(); clusterIt != clusters.end(); ++clusterIt) {
    std::vector<art::Ptr<recob::Hit> > hits = fmhc.at(clusterIt->key());
    int trueParticle = FindTrueParticle(clockData, hits);
    std::vector<art::Ptr<recob::Hit> > trueHits = FindTrueHits(clockData, hits, trueParticle);
    if (particles.count(trueParticle))
      particles[trueParticle]->AddAssociatedCluster(*clusterIt, hits, trueHits);
  }

  for (std::vector<art::Ptr<recob::Shower> >::iterator showerIt = showers.begin(); showerIt != showers.end(); ++showerIt) {
    std::vector<art::Ptr<recob::Hit> > hits = fmhs.at(showerIt->key());
    int trueParticle = FindTrueParticle(clockData, hits);
    std::vector<art::Ptr<recob::Hit> > trueHits = FindTrueHits(clockData, hits, trueParticle);
    if (particles.count(trueParticle))
      particles[trueParticle]->AddAssociatedShower(*showerIt, hits, trueHits);
  }

  // Fill output data products
  this->FillData(particles);
  if (isPi0 and pi0Decays.size() == 2)
    this->FillPi0Data(particles, pi0Decays);

  return;

}

void showerAna::ShowerAnalysis::FillData

(

const std::map< int, std::shared_ptr< ShowerParticle > > &

particles

)

Definition at line 365 of file ShowerAnalysis_module.cc.

                                                                                                  {

  // Look at each particle
  for (std::map<int,std::shared_ptr<ShowerParticle> >::const_iterator particleIt = particles.begin(); particleIt != particles.end(); ++particleIt) {

    std::shared_ptr<ShowerParticle> particle = particleIt->second;

    // No point making these plots for particles which don't shower!
    if (TMath::Abs(particle->PDG()) != 11 and TMath::Abs(particle->PDG()) != 22)
      continue;

    // Cluster plots
    for (int cluster = 0; cluster < particle->NumClusters(); ++cluster) {
      hClusterCompleteness              ->Fill(particle->ClusterCompleteness(cluster));
      hClusterPurity                    ->Fill(particle->ClusterPurity(cluster));
      hClusterCompletenessEnergy        ->Fill(particle->Energy(), particle->ClusterCompleteness(cluster));
      hClusterCompletenessDirection     ->Fill(particle->Direction().Angle(TVector3(0,1,0)), particle->ClusterCompleteness(cluster));
      if (particle->LargestCluster(cluster)) {
        hLargestClusterCompleteness             ->Fill(particle->ClusterCompleteness(cluster));
        hLargestClusterPurity                   ->Fill(particle->ClusterPurity(cluster));
        hLargestClusterCompletenessEnergy       ->Fill(particle->Energy(), particle->ClusterCompleteness(cluster));
        hLargestClusterCompletenessDirection    ->Fill(particle->Direction().Angle(TVector3(0,1,0)), particle->ClusterCompleteness(cluster));
      }
    }

    // Shower plots
    for (int shower = 0; shower < particle->NumShowers(); ++shower) {
      hShowerCompleteness               ->Fill(particle->ShowerCompleteness(shower));
      hShowerPurity                     ->Fill(particle->ShowerPurity(shower));
      hShowerCompletenessEnergy         ->Fill(particle->Energy(), particle->ShowerCompleteness(shower));
      hShowerCompletenessDirection      ->Fill(particle->Direction().Angle(TVector3(0,1,0)), particle->ShowerCompleteness(shower));
      if (particle->LargestShower(shower)) {
        hLargestShowerCompleteness              ->Fill(particle->ShowerCompleteness(shower));
        hLargestShowerPurity                    ->Fill(particle->ShowerPurity(shower));
        hLargestShowerCompletenessEnergy        ->Fill(particle->Energy(), particle->ShowerCompleteness(shower));
        hLargestShowerCompletenessDirection     ->Fill(particle->Direction().Angle(TVector3(0,1,0)), particle->ShowerCompleteness(shower));
        hShowerEnergy                           ->Fill(particle->ShowerEnergy() / particle->Energy());
        hShowerDepositedEnergy                  ->Fill(particle->ShowerEnergy() / particle->DepositedEnergy());
        hShowerEnergyCompleteness               ->Fill(particle->Energy(), particle->ShowerEnergy()/particle->DepositedEnergy());
        hShowerDirection                        ->Fill(particle->Direction().Dot(particle->ShowerDirection()));
        hShowerdEdx                             ->Fill(particle->ShowerdEdx());
        if (particle->PDG() == 22) {
          hShowerStart                          ->Fill((particle->End() - particle->ShowerStart()).Mag());
          hShowerStartEnergy                    ->Fill(particle->Energy(), (particle->End() - particle->ShowerStart()).Mag());
        }
        else {
          hShowerStart                          ->Fill((particle->Start() - particle->ShowerStart()).Mag());
          hShowerStartEnergy                    ->Fill(particle->Energy(), (particle->Start() - particle->ShowerStart()).Mag());
        }
      }
    }

    hShowerReconstructed->Fill(particle->NumShowers() != 0);
    hNumShowersReconstructed->Fill(particle->NumShowers());
    hShowerReconstructedEnergy->Fill(particle->Energy(), particle->NumShowers() != 0);
    hNumShowersReconstructedEnergy->Fill(particle->Energy(), particle->NumShowers());
    hNumShowersReconstructedEnergyProfile->Fill(particle->Energy(), particle->NumShowers());
    if (particle->NumShowers()) {
      hNumShowersReconstructedNonZero->Fill(particle->NumShowers());
      hShowerdEdxEnergy->Fill(particle->Energy(), particle->ShowerdEdx());
      hShowerdEdxEnergyProfile->Fill(particle->Energy(), particle->ShowerdEdx());
      hElectronPull->Fill((particle->ShowerdEdx() - fElectrondEdx)/fElectrondEdxWidth);
      hPhotonPull->Fill((particle->ShowerdEdx() - fPhotondEdx)/fPhotondEdxWidth);
    }

  }

  return;

}

void showerAna::ShowerAnalysis::FillPi0Data	(	const std::map< int, std::shared_ptr< ShowerParticle > > &	particles,
		const std::vector< int > &	pi0Decays
	)

Fill data products related to pi0s We only call this function if there are exactly two decay products, so we can assume the decay was pi0 -> gamma gamma

Definition at line 436 of file ShowerAnalysis_module.cc.

                                                                                                                                      {

  /// Fill data products related to pi0s
  /// We only call this function if there are exactly two decay products, so we can assume the decay was pi0 -> gamma gamma

  std::shared_ptr<ShowerParticle> photon1 = particles.at(pi0Decays.at(0));
  std::shared_ptr<ShowerParticle> photon2 = particles.at(pi0Decays.at(1));

  // Make sure there is at least one shower and it was fully reconstructed
  if (photon1->NumShowers() == 0 or photon2->NumShowers() == 0 or
      photon1->ShowerStart() == TVector3(0,0,0) or photon2->ShowerStart() == TVector3(0,0,0)) {
    std::cout << "Event doesn't have enough info for pi0 mass determination" << std::endl;
    return;
  }

  double reconOpeningAngle = TMath::ASin(TMath::Sin(photon1->ShowerDirection().Angle(photon2->ShowerDirection())));
  double trueOpeningAngle = photon1->Direction().Angle(photon2->Direction());

  // Reconstruct the pi0 mass
  double massReconEnergyReconAngle = TMath::Sqrt(4 * photon1->ShowerEnergy() * photon2->ShowerEnergy() * TMath::Power(TMath::Sin(reconOpeningAngle/2),2));
  double massReconEnergyTrueAngle = TMath::Sqrt(4 * photon1->ShowerEnergy() * photon2->ShowerEnergy() * TMath::Power(TMath::Sin(trueOpeningAngle/2),2));
  double massTrueEnergyReconAngle = TMath::Sqrt(4 * photon1->Energy() * photon2->Energy() * TMath::Power(TMath::Sin(reconOpeningAngle/2),2));
  double massTrueEnergyTrueAngle = TMath::Sqrt(4 * photon1->Energy() * photon2->Energy() * TMath::Power(TMath::Sin(trueOpeningAngle/2),2));

  hPi0MassPeakReconEnergyReconAngle->Fill(massReconEnergyReconAngle);
  hPi0MassPeakReconEnergyTrueAngle->Fill(massReconEnergyTrueAngle);
  hPi0MassPeakTrueEnergyReconAngle->Fill(massTrueEnergyReconAngle);
  hPi0MassPeakTrueEnergyTrueAngle->Fill(massTrueEnergyTrueAngle);

  return;
  
}

int showerAna::ShowerAnalysis::FindParticleID	(	detinfo::DetectorClocksData const &	clockData,
		const art::Ptr< recob::Hit > &	hit
	)

Returns the true track ID associated with this hit (if more than one, returns the one with highest energy)

Definition at line 496 of file ShowerAnalysis_module.cc.

                                                                           {

  /// Returns the true track ID associated with this hit (if more than one, returns the one with highest energy)

  double particleEnergy = 0;
  int likelyTrackID = 0;
  std::vector<sim::TrackIDE> trackIDs = bt_serv->HitToTrackIDEs(clockData, hit);
  for (unsigned int idIt = 0; idIt < trackIDs.size(); ++idIt) {
    if (trackIDs.at(idIt).energy > particleEnergy) {
      particleEnergy = trackIDs.at(idIt).energy;
      likelyTrackID = TMath::Abs(trackIDs.at(idIt).trackID);
    }
  }

  return likelyTrackID;

}

std::vector< art::Ptr< recob::Hit > > showerAna::ShowerAnalysis::FindTrueHits	(	detinfo::DetectorClocksData const &	clockData,
		const std::vector< art::Ptr< recob::Hit > > &	hits,
		int	trueParticle
	)

Definition at line 515 of file ShowerAnalysis_module.cc.

                                                                                                                                       {

  std::vector<art::Ptr<recob::Hit> > trueHits;
  for (std::vector<art::Ptr<recob::Hit> >::const_iterator hitIt = hits.begin(); hitIt != hits.end(); ++hitIt)
    if (FindParticleID(clockData, *hitIt) == trueParticle)
      trueHits.push_back(*hitIt);

  return trueHits;

}

int showerAna::ShowerAnalysis::FindTrueParticle	(	detinfo::DetectorClocksData const &	clockData,
		const std::vector< art::Ptr< recob::Hit > > &	hits
	)

Returns the true particle most likely associated with this shower

Definition at line 469 of file ShowerAnalysis_module.cc.

                                                                                                {

  /// Returns the true particle most likely associated with this shower

  // Make a map of the tracks which are associated with this shower and the charge each contributes
  std::map<int,double> trackMap;
  for (std::vector<art::Ptr<recob::Hit> >::const_iterator showerHitIt = showerHits.begin(); showerHitIt != showerHits.end(); ++showerHitIt) {
    art::Ptr<recob::Hit> hit = *showerHitIt;
    int trackID = FindParticleID(clockData, hit);
    trackMap[trackID] += hit->Integral();
  }

  // Pick the track with the highest charge as the 'true track'
  double highestCharge = 0;
  int showerTrack = 0;
  for (std::map<int,double>::iterator trackIt = trackMap.begin(); trackIt != trackMap.end(); ++trackIt) {
    if (trackIt->second > highestCharge) {
      highestCharge = trackIt->second;
      showerTrack  = trackIt->first;
    }
  }

  return showerTrack;

}

void showerAna::ShowerAnalysis::MakeDataProducts

(

)

Definition at line 527 of file ShowerAnalysis_module.cc.

                                               {

  fTree = tfs->make<TTree>("ShowerAnalysis","ShowerAnalysis");

  // Cluster
  hClusterCompleteness = tfs->make<TH1D>("ClusterCompleteness","Completeness of all clusters;Completeness;",101,0,1.01);
  hLargestClusterCompleteness = tfs->make<TH1D>("LargestClusterCompleteness","Completeness of largest cluster;Completeness;",101,0,1.01);
  hClusterPurity = tfs->make<TH1D>("ClusterPurity","Purity of all clusters;Purity;",101,0,1.01);
  hLargestClusterPurity = tfs->make<TH1D>("LargestClusterPurity","Purity of largest cluster;Purity;",101,0,1.01);
  hClusterCompletenessEnergy = tfs->make<TH2D>("ClusterCompletenessEnergy","Completeness of all clusters vs Energy;Energy (GeV);Completeness;",100,0,10,101,0,1.01);
  hLargestClusterCompletenessEnergy = tfs->make<TH2D>("LargestClusterCompletenessEnergy","Completeness of largest cluster vs Energy;Energy (GeV);Completeness;",100,0,10,101,0,1.01);
  hClusterCompletenessDirection = tfs->make<TH2D>("ClusterCompletenessDirection","Completeness of all clusters vs Direction;Direction;Completeness;",100,0,5,101,0,1.01);
  hLargestClusterCompletenessDirection = tfs->make<TH2D>("LargestClusterCompletenessDirection","Completeness of largest cluster vs Direction;Direction;Completeness;",100,0,5,101,0,1.01);

  // Shower
  hShowerCompleteness = tfs->make<TH1D>("ShowerCompleteness","Completeness of all showers;Completeness;",101,0,1.01);
  hLargestShowerCompleteness = tfs->make<TH1D>("LargestShowerCompleteness","Completeness of largest shower;Completeness;",101,0,1.01);
  hShowerPurity = tfs->make<TH1D>("ShowerPurity","Purity of all showers;Purity;",101,0,1.01);
  hLargestShowerPurity = tfs->make<TH1D>("LargestShowerPurity","Purity of largest shower;Purity;",101,0,1.01);
  hShowerCompletenessEnergy = tfs->make<TH2D>("ShowerCompletenessEnergy","Completeness of all showers vs Energy;Energy (GeV);Completeness;",100,0,10,101,0,1.01);
  hLargestShowerCompletenessEnergy = tfs->make<TH2D>("LargestShowerCompletenessEnergy","Completeness of largest shower vs Energy;Energy (GeV);Completeness;",100,0,10,101,0,1.01);
  hShowerCompletenessDirection = tfs->make<TH2D>("ShowerCompletenessDirection","Completeness of all showers vs Direction;Direction;Completeness;",100,0,5,101,0,1.01);
  hLargestShowerCompletenessDirection = tfs->make<TH2D>("LargestShowerCompletenessDirection","Completeness of largest shower vs Direction;Direction;Completeness;",100,0,5,101,0,1.01);
  hShowerEnergy = tfs->make<TH1D>("ShowerEnergy","Shower energy;Recon Energy/True Energy;",120,0,1.2);
  hShowerDepositedEnergy = tfs->make<TH1D>("ShowerDepositedEnergy","Shower deposited energy;Recon Energy/True (Deposited) Energy;",120,0,1.2);
  hShowerDirection = tfs->make<TH1D>("ShowerDirection","Shower direction;True Direction.(Recon Direction);",202,-1.01,1.01);
  hShowerdEdx = tfs->make<TH1D>("ShowerdEdx","dEdx of Shower;dE/dx (MeV/cm)",50,0,10);
  hShowerStart = tfs->make<TH1D>("ShowerStart","Distance from reconstructed shower start to true shower start;True Start - Recon Start (cm);",100,0,20);
  hShowerReconstructed = tfs->make<TH1D>("ShowerReconstructed","% of showering particles with reconstructed shower",101,0,1.01);
  hShowerEnergyCompleteness = tfs->make<TProfile>("ShowerEnergyCompleteness","Shower energy completeness vs true deposited energy;True (Deposited) Energy (GeV);Energy Completeness;",100,0,5);
  hShowerStartEnergy = tfs->make<TProfile>("ShowerStartEnergy","Shower start distance vs true deposited energy;True (Deposited) Energy (GeV);True Start - Recon Start (cm);",100,0,5);
  hNumShowersReconstructed = tfs->make<TH1D>("NumShowersReconstructed","Number of showers reconstructed for each showering particle;Number of Showers",10,0,10);
  hNumShowersReconstructedNonZero = tfs->make<TH1D>("NumShowersReconstructedNonZero","Number of showers reconstructed for each showering particle;Number of Showers",9,1,10);
  for (int n_showers = 0; n_showers < 9; ++n_showers) {
    hNumShowersReconstructed->GetXaxis()->SetBinLabel(n_showers+1,std::to_string(n_showers).c_str());
    hNumShowersReconstructedNonZero->GetXaxis()->SetBinLabel(n_showers+1,std::to_string(n_showers+1).c_str());
  }
  hNumShowersReconstructed->GetXaxis()->SetBinLabel(10,"9 or more");
  hNumShowersReconstructedNonZero->GetXaxis()->SetBinLabel(9,"9 or more");
  hShowerReconstructedEnergy = tfs->make<TProfile>("ShowerReconstructedEnergyProfile","% of showering particles with reconstructed shower vs true energy;True Energy (GeV);Fraction of particles with reconstructed shower;",100,0,5);
  hNumShowersReconstructedEnergy = tfs->make<TH2D>("NumShowersReconstructedEnergy","Number of showers reconstructed per showering particle vs true energy;True Energy (GeV);Average Number of Showers;",100,0,5,10,0,10);
  hNumShowersReconstructedEnergyProfile = tfs->make<TProfile>("NumShowersReconstructedEnergyProfile",";True Energy (GeV);Number of Showers;",100,0,5);
  hShowerdEdxEnergy = tfs->make<TH2D>("ShowerdEdxEnergy","Shower dE/dx vs true energy;True Energy (GeV);dE/dx (MeV/cm);",100,0,5,50,0,10);
  hShowerdEdxEnergyProfile = tfs->make<TProfile>("ShowerdEdxEnergyProfile",";True Energy (GeV);dE/dx (MeV/cm);",100,0,5);
  hElectronPull = tfs->make<TH1D>("ElectronPull","Electron pull;Electron pull;",100,-10,10);
  hPhotonPull = tfs->make<TH1D>("PhotonPull","Photon pull;Photon pull;",100,-10,10);

  // pi0
  hPi0MassPeakReconEnergyReconAngle = tfs->make<TH1D>("Pi0MassPeakReconEnergyReconAngle",";Invariant Mass (GeV);",40,0,0.5);
  hPi0MassPeakTrueEnergyReconAngle = tfs->make<TH1D>("Pi0MassPeakTrueEnergyReconAngle",";Invariant Mass (GeV);",40,0,0.5);
  hPi0MassPeakReconEnergyTrueAngle = tfs->make<TH1D>("Pi0MassPeakReconEnergyTrueAngle",";Invariant Mass (GeV);",40,0,0.5);
  hPi0MassPeakTrueEnergyTrueAngle = tfs->make<TH1D>("Pi0MassPeakTrueEnergyTrueAngle",";Invariant Mass (GeV);",40,0,0.5);

}

Member Data Documentation

art::ServiceHandle<cheat::BackTrackerService> showerAna::ShowerAnalysis::bt_serv

private

Definition at line 218 of file ShowerAnalysis_module.cc.

std::string showerAna::ShowerAnalysis::fClusterModuleLabel

private

Definition at line 214 of file ShowerAnalysis_module.cc.

double showerAna::ShowerAnalysis::fElectrondEdx

private

Definition at line 216 of file ShowerAnalysis_module.cc.

double showerAna::ShowerAnalysis::fElectrondEdxWidth

private

Definition at line 216 of file ShowerAnalysis_module.cc.

std::string showerAna::ShowerAnalysis::fHitsModuleLabel

private

Definition at line 214 of file ShowerAnalysis_module.cc.

double showerAna::ShowerAnalysis::fPhotondEdx

private

Definition at line 216 of file ShowerAnalysis_module.cc.

double showerAna::ShowerAnalysis::fPhotondEdxWidth

private

Definition at line 216 of file ShowerAnalysis_module.cc.

std::string showerAna::ShowerAnalysis::fShowerModuleLabel

private

Definition at line 214 of file ShowerAnalysis_module.cc.

TTree* showerAna::ShowerAnalysis::fTree

private

Definition at line 223 of file ShowerAnalysis_module.cc.

art::ServiceHandle<geo::Geometry> showerAna::ShowerAnalysis::geom

private

Definition at line 221 of file ShowerAnalysis_module.cc.

TH1D* showerAna::ShowerAnalysis::hClusterCompleteness

private

Definition at line 226 of file ShowerAnalysis_module.cc.

TH2D * showerAna::ShowerAnalysis::hClusterCompletenessDirection

private

Definition at line 227 of file ShowerAnalysis_module.cc.

TH2D* showerAna::ShowerAnalysis::hClusterCompletenessEnergy

private

Definition at line 227 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hClusterPurity

private

Definition at line 226 of file ShowerAnalysis_module.cc.

TH1D* showerAna::ShowerAnalysis::hElectronPull

private

Definition at line 233 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hLargestClusterCompleteness

private

Definition at line 226 of file ShowerAnalysis_module.cc.

TH2D * showerAna::ShowerAnalysis::hLargestClusterCompletenessDirection

private

Definition at line 227 of file ShowerAnalysis_module.cc.

TH2D * showerAna::ShowerAnalysis::hLargestClusterCompletenessEnergy

private

Definition at line 227 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hLargestClusterPurity

private

Definition at line 226 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hLargestShowerCompleteness

private

Definition at line 228 of file ShowerAnalysis_module.cc.

TH2D * showerAna::ShowerAnalysis::hLargestShowerCompletenessDirection

private

Definition at line 229 of file ShowerAnalysis_module.cc.

TH2D * showerAna::ShowerAnalysis::hLargestShowerCompletenessEnergy

private

Definition at line 229 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hLargestShowerPurity

private

Definition at line 228 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hNumShowersReconstructed

private

Definition at line 230 of file ShowerAnalysis_module.cc.

TH2D * showerAna::ShowerAnalysis::hNumShowersReconstructedEnergy

private

Definition at line 231 of file ShowerAnalysis_module.cc.

TProfile * showerAna::ShowerAnalysis::hNumShowersReconstructedEnergyProfile

private

Definition at line 232 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hNumShowersReconstructedNonZero

private

Definition at line 230 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hPhotonPull

private

Definition at line 233 of file ShowerAnalysis_module.cc.

TH1D* showerAna::ShowerAnalysis::hPi0MassPeakReconEnergyReconAngle

private

Definition at line 236 of file ShowerAnalysis_module.cc.

TH1D* showerAna::ShowerAnalysis::hPi0MassPeakReconEnergyTrueAngle

private

Definition at line 237 of file ShowerAnalysis_module.cc.

TH1D* showerAna::ShowerAnalysis::hPi0MassPeakTrueEnergyReconAngle

private

Definition at line 238 of file ShowerAnalysis_module.cc.

TH1D* showerAna::ShowerAnalysis::hPi0MassPeakTrueEnergyTrueAngle

private

Definition at line 239 of file ShowerAnalysis_module.cc.

TH1D* showerAna::ShowerAnalysis::hShowerCompleteness

private

Definition at line 228 of file ShowerAnalysis_module.cc.

TH2D * showerAna::ShowerAnalysis::hShowerCompletenessDirection

private

Definition at line 229 of file ShowerAnalysis_module.cc.

TH2D* showerAna::ShowerAnalysis::hShowerCompletenessEnergy

private

Definition at line 229 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hShowerdEdx

private

Definition at line 230 of file ShowerAnalysis_module.cc.

TH2D* showerAna::ShowerAnalysis::hShowerdEdxEnergy

private

Definition at line 231 of file ShowerAnalysis_module.cc.

TProfile * showerAna::ShowerAnalysis::hShowerdEdxEnergyProfile

private

Definition at line 232 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hShowerDepositedEnergy

private

Definition at line 230 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hShowerDirection

private

Definition at line 230 of file ShowerAnalysis_module.cc.

TH1D* showerAna::ShowerAnalysis::hShowerEnergy

private

Definition at line 230 of file ShowerAnalysis_module.cc.

TProfile * showerAna::ShowerAnalysis::hShowerEnergyCompleteness

private

Definition at line 232 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hShowerPurity

private

Definition at line 228 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hShowerReconstructed

private

Definition at line 230 of file ShowerAnalysis_module.cc.

TProfile* showerAna::ShowerAnalysis::hShowerReconstructedEnergy

private

Definition at line 232 of file ShowerAnalysis_module.cc.

TH1D * showerAna::ShowerAnalysis::hShowerStart

private

Definition at line 230 of file ShowerAnalysis_module.cc.

TProfile * showerAna::ShowerAnalysis::hShowerStartEnergy

private

Definition at line 232 of file ShowerAnalysis_module.cc.

art::ServiceHandle<cheat::ParticleInventoryService> showerAna::ShowerAnalysis::pi_serv

private

Definition at line 219 of file ShowerAnalysis_module.cc.

art::ServiceHandle<art::TFileService> showerAna::ShowerAnalysis::tfs

private

Definition at line 220 of file ShowerAnalysis_module.cc.

---

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

• duneana/duneana/ShowerAna/ShowerAnalysis_module.cc