Public Member Functions | List of all members
protoana::ProtoDUNEShowerUtils Class Reference

#include <ProtoDUNEShowerUtils.h>

Public Member Functions

 ProtoDUNEShowerUtils ()
 
 ~ProtoDUNEShowerUtils ()
 
const std::vector< const recob::Hit * > GetRecoShowerHits (const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
 Get the hits from a given reco shower. More...
 
const std::vector< art::Ptr< recob::Hit > > GetRecoShowerArtHits (const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
 
unsigned int GetNumberRecoShowerHits (const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
 Get the number of hits from a given reco shower. More...
 
std::vector< const recob::PCAxis * > GetRecoShowerPCAxis (const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
 Get the associated PCAxis object (from a principal component analysis) More...
 
std::vector< double > EstimateEnergyFromHitCharge (detinfo::DetectorClocksData const &clockData, detinfo::DetectorPropertiesData const &detProp, const std::vector< const recob::Hit * > &hits, calo::CalorimetryAlg caloAlg)
 
int GetShowerIndex (const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
 If the shower.ID() isn't filled we must find the actual shower index ourselves. More...
 
std::vector< anab::Calorimetry > GetRecoShowerCalorimetry (const recob::Shower &shower, art::Event const &evt, const std::string showerModule, const std::string caloModule) const
 Get shower calo info. More...
 

Detailed Description

Definition at line 29 of file ProtoDUNEShowerUtils.h.

Constructor & Destructor Documentation

protoana::ProtoDUNEShowerUtils::ProtoDUNEShowerUtils ( )

Definition at line 8 of file ProtoDUNEShowerUtils.cxx.

8  {
9 
10 }
protoana::ProtoDUNEShowerUtils::~ProtoDUNEShowerUtils ( )

Definition at line 12 of file ProtoDUNEShowerUtils.cxx.

12  {
13 
14 }

Member Function Documentation

std::vector< double > protoana::ProtoDUNEShowerUtils::EstimateEnergyFromHitCharge ( detinfo::DetectorClocksData const &  clockData,
detinfo::DetectorPropertiesData const &  detProp,
const std::vector< const recob::Hit * > &  hits,
calo::CalorimetryAlg  caloAlg 
)

Definition at line 78 of file ProtoDUNEShowerUtils.cxx.

81 {
82  double kGeVtoElectrons { 4.237e7 }; // obtained from utils class.. Copied for now, should use class (although this is a physical constant, so hopefully doesn't change).
83  double recombination { 1/0.63 };
84 
85  std::vector<double> showerEnergy = {0,0,0};
86 
87  // Find the total charge on each plane
88  for ( size_t h{0} ; h < hits.size() ; h++ ) {
89  const recob::Hit* hit = hits[h];
90  const int plane = hit->WireID().Plane;
91  showerEnergy[ plane ] += ( caloAlg.ElectronsFromADCArea( hit->Integral(), plane) * caloAlg.LifetimeCorrection(clockData, detProp, hit->PeakTime()) ) / kGeVtoElectrons;
92  }
93 
94  showerEnergy[0] *= recombination;
95  showerEnergy[1] *= recombination;
96  showerEnergy[2] *= recombination;
97 
98  // caloAlg.ElectronsFromADCArea( hit->Integral(), plane) -> Does hit->Integral()/AreaConstants(plane)
99  // AreaConstants(plane) is defined in calorimetry_pdune.fcl. Although that fcl file has a typo.
100  // These probably need tuning for protodune data.
101 
102  return showerEnergy;
103 }
recob::Hit::WireID
geo::WireID WireID() const
Definition: Hit.h:233
recob::Hit::Integral
float Integral() const
Integral under the calibrated signal waveform of the hit, in tick x ADC units.
Definition: Hit.h:224
calo::CalorimetryAlg::ElectronsFromADCArea
double ElectronsFromADCArea(double area, unsigned short plane) const
Definition: CalorimetryAlg.h:153
reco_momentum_tuples.h
h
Definition: reco_momentum_tuples.py:66
recombination
double recombination
Definition: doAna.cpp:19
geo::PlaneID::Plane
PlaneID_t Plane
Index of the plane within its TPC.
Definition: geo_types.h:493
hit
Detector simulation of raw signals on wires.
Definition: NumberOfHitsFilter_module.cc:25
recob::Hit::PeakTime
float PeakTime() const
Time of the signal peak, in tick units.
Definition: Hit.h:218
recob::Hit
2D representation of charge deposited in the TDC/wire plane
Definition: Hit.h:48
calo::CalorimetryAlg::LifetimeCorrection
double LifetimeCorrection(detinfo::DetectorClocksData const &clock_data, detinfo::DetectorPropertiesData const &det_prop, double time, double T0=0) const
Definition: CalorimetryAlg.cxx:246
unsigned int protoana::ProtoDUNEShowerUtils::GetNumberRecoShowerHits ( const recob::Shower &  shower,
art::Event const &  evt,
const std::string  showerModule 
) const

Get the number of hits from a given reco shower.

Definition at line 55 of file ProtoDUNEShowerUtils.cxx.

55  {
56 
57  return GetRecoShowerHits(shower,evt,showerModule).size();
58 
59 }
protoana::ProtoDUNEShowerUtils::GetRecoShowerHits
const std::vector< const recob::Hit * > GetRecoShowerHits(const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
Get the hits from a given reco shower.
Definition: ProtoDUNEShowerUtils.cxx:17
tca::evt
TCEvent evt
Definition: DataStructs.cxx:7
const std::vector< art::Ptr< recob::Hit > > protoana::ProtoDUNEShowerUtils::GetRecoShowerArtHits ( const recob::Shower &  shower,
art::Event const &  evt,
const std::string  showerModule 
) const

Definition at line 40 of file ProtoDUNEShowerUtils.cxx.

41  {
42 
43  auto recoShowers = evt.getValidHandle<std::vector<recob::Shower> >(showerModule);
44  art::FindManyP<recob::Hit> findHits(recoShowers,evt,showerModule);
45 
46  // Shower.ID is sometimes at a default value - make sure we get the correct one
47  int actualIndex = GetShowerIndex(shower,evt,showerModule);
48 
49  std::vector<art::Ptr<recob::Hit>> showerHits = findHits.at(actualIndex);
50  return showerHits;
51 
52 }
protoana::ProtoDUNEShowerUtils::GetShowerIndex
int GetShowerIndex(const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
If the shower.ID() isn&#39;t filled we must find the actual shower index ourselves.
Definition: ProtoDUNEShowerUtils.cxx:106
tca::evt
TCEvent evt
Definition: DataStructs.cxx:7
std::vector< anab::Calorimetry > protoana::ProtoDUNEShowerUtils::GetRecoShowerCalorimetry ( const recob::Shower &  shower,
art::Event const &  evt,
const std::string  showerModule,
const std::string  caloModule 
) const

Get shower calo info.

Definition at line 130 of file ProtoDUNEShowerUtils.cxx.

130  {
131 
132  auto recoShowers = evt.getValidHandle<std::vector<recob::Shower> >(showerModule);
133  std::vector<anab::Calorimetry> caloInfo;
134 
135  try{
136  const art::FindManyP<anab::Calorimetry> findCalorimetry(recoShowers,evt,caloModule);
137  int actualIndex = GetShowerIndex(shower,evt,showerModule);
138  std::vector<art::Ptr<anab::Calorimetry>> theseCalos = findCalorimetry.at(actualIndex);
139 
140  for( auto calo : theseCalos){
141  caloInfo.push_back(*calo);
142  }
143  }
144  catch(...){
145  std::cerr << "No shower calorimetry object found... returning empty vector" << std::endl;
146  }
147 
148  return caloInfo;
149 
150 
151 }
protoana::ProtoDUNEShowerUtils::GetShowerIndex
int GetShowerIndex(const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
If the shower.ID() isn&#39;t filled we must find the actual shower index ourselves.
Definition: ProtoDUNEShowerUtils.cxx:106
tca::evt
TCEvent evt
Definition: DataStructs.cxx:7
calo
calorimetry
Definition: CaloChecker_module.cc:43
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030
const std::vector< const recob::Hit * > protoana::ProtoDUNEShowerUtils::GetRecoShowerHits ( const recob::Shower &  shower,
art::Event const &  evt,
const std::string  showerModule 
) const

Get the hits from a given reco shower.

Definition at line 17 of file ProtoDUNEShowerUtils.cxx.

17  {
18 
19  auto recoShowers = evt.getValidHandle<std::vector<recob::Shower> >(showerModule);
20  art::FindManyP<recob::Hit> findHits(recoShowers,evt,showerModule);
21 
22  // Shower.ID is sometimes at a default value - make sure we get the correct one
23  int actualIndex = GetShowerIndex(shower,evt,showerModule);
24 
25  std::vector<art::Ptr<recob::Hit>> inputHits = findHits.at(actualIndex);
26  std::vector<const recob::Hit*> showerHits;
27 
28  for(const art::Ptr<recob::Hit> hit : inputHits){
29 
30  showerHits.push_back(hit.get());
31 
32  }
33 
34  return showerHits;
35 
36 }
protoana::ProtoDUNEShowerUtils::GetShowerIndex
int GetShowerIndex(const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
If the shower.ID() isn&#39;t filled we must find the actual shower index ourselves.
Definition: ProtoDUNEShowerUtils.cxx:106
hit
Detector simulation of raw signals on wires.
Definition: NumberOfHitsFilter_module.cc:25
tca::evt
TCEvent evt
Definition: DataStructs.cxx:7
std::vector< const recob::PCAxis * > protoana::ProtoDUNEShowerUtils::GetRecoShowerPCAxis ( const recob::Shower &  shower,
art::Event const &  evt,
const std::string  showerModule 
) const

Get the associated PCAxis object (from a principal component analysis)

Definition at line 62 of file ProtoDUNEShowerUtils.cxx.

62  {
63 
64  auto recoShowers = evt.getValidHandle<std::vector<recob::Shower> >(showerModule);
65  art::FindManyP<recob::PCAxis> findPCA(recoShowers,evt,showerModule);
66 
67  // Shower.ID is sometimes at a default value - make sure we get the correct one
68  int actualIndex = GetShowerIndex(shower,evt,showerModule);
69 
70  std::vector<const recob::PCAxis*> pcaVec;
71  for(auto const pca : findPCA.at(actualIndex)){
72  pcaVec.push_back(pca.get());
73  }
74 
75  return pcaVec;
76 }
protoana::ProtoDUNEShowerUtils::GetShowerIndex
int GetShowerIndex(const recob::Shower &shower, art::Event const &evt, const std::string showerModule) const
If the shower.ID() isn&#39;t filled we must find the actual shower index ourselves.
Definition: ProtoDUNEShowerUtils.cxx:106
tca::evt
TCEvent evt
Definition: DataStructs.cxx:7
int protoana::ProtoDUNEShowerUtils::GetShowerIndex ( const recob::Shower &  shower,
art::Event const &  evt,
const std::string  showerModule 
) const

If the shower.ID() isn't filled we must find the actual shower index ourselves.

Definition at line 106 of file ProtoDUNEShowerUtils.cxx.

106  {
107 
108  if(shower.ID() != -999) return shower.ID();
109 
110  auto recoShowers = evt.getValidHandle<std::vector<recob::Shower> >(showerModule);
111 
112  // Iterate through all showers to find the matching one to our shower
113  int actualIndex = shower.ID();
114  if(shower.ID() < 0){
115  for(unsigned int s = 0; s < recoShowers->size(); ++s){
116  const recob::Shower thisShower = (*recoShowers)[s];
117  // Can't compare actual objects so look at a property
118  if(fabs(thisShower.Length() - shower.Length()) < 1e-5){
119  actualIndex = s;
120  continue;
121  }
122  }
123  }
124 
125  return actualIndex;
126 
127 }
recob::Shower::Length
double Length() const
Definition: Shower.h:201
QCString::size
uint size() const
Definition: qcstring.h:201
e
const double e
Definition: gUpMuFluxGen.cxx:165
tca::evt
TCEvent evt
Definition: DataStructs.cxx:7
recob::Shower::ID
int ID() const
Definition: Shower.h:187
s
static QCString * s
Definition: config.cpp:1042
The documentation for this class was generated from the following files:
Generated by   doxygen 1.8.11