LArSoft  v06_85_00
Liquid Argon Software toolkit - http://larsoft.org/
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | List of all members
calo::Calorimetry Class Reference
Inheritance diagram for calo::Calorimetry:
art::EDProducer art::ProducerBase art::Consumer art::EngineCreator art::ProductRegistryHelper

Public Types

using ModuleType = EDProducer
 
using WorkerType = WorkerT< EDProducer >
 
template<typename UserConfig , typename KeysToIgnore = void>
using Table = ProducerBase::Table< UserConfig, KeysToIgnore >
 

Public Member Functions

 Calorimetry (fhicl::ParameterSet const &pset)
 
virtual ~Calorimetry ()
 
void beginJob ()
 
void produce (art::Event &evt)
 
template<typename PROD , BranchType B = InEvent>
ProductID getProductID (std::string const &instanceName={}) const
 
template<typename PROD , BranchType B>
ProductID getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const
 
bool modifiesEvent () const
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< T > consumes (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< T > consumesView (InputTag const &it)
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< T > mayConsume (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< T > mayConsumeView (InputTag const &it)
 
base_engine_t & createEngine (seed_t seed)
 
base_engine_t & createEngine (seed_t seed, std::string const &kind_of_engine_to_make)
 
base_engine_t & createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
 
seed_t get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)
 

Static Public Member Functions

static cet::exempt_ptr< Consumer > non_module_context ()
 

Protected Member Functions

CurrentProcessingContext const * currentContext () const
 
void validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
 
void prepareForJob (fhicl::ParameterSet const &pset)
 
void showMissingConsumes () const
 

Private Member Functions

void ReadCaloTree ()
 
bool BeginsOnBoundary (art::Ptr< recob::Track > lar_track)
 
bool EndsOnBoundary (art::Ptr< recob::Track > lar_track)
 
void GetPitch (art::Ptr< recob::Hit > hit, std::vector< double > trkx, std::vector< double > trky, std::vector< double > trkz, std::vector< double > trkw, std::vector< double > trkx0, double *xyz3d, double &pitch, double TickT0)
 

Private Attributes

std::string fTrackModuleLabel
 
std::string fSpacePointModuleLabel
 
std::string fT0ModuleLabel
 
bool fUseArea
 
bool fFlipTrack_dQdx
 
CalorimetryAlg caloAlg
 
int fnsps
 
std::vector< int > fwire
 
std::vector< double > ftime
 
std::vector< double > fstime
 
std::vector< double > fetime
 
std::vector< double > fMIPs
 
std::vector< double > fdQdx
 
std::vector< double > fdEdx
 
std::vector< double > fResRng
 
std::vector< double > fpitch
 
std::vector< TVector3 > fXYZ
 

Detailed Description

Definition at line 75 of file Calorimetry_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer
inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>
using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>
inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>
inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

calo::Calorimetry::Calorimetry ( fhicl::ParameterSet const &  pset)
explicit

Definition at line 125 of file Calorimetry_module.cc.

126  : fTrackModuleLabel(pset.get< std::string >("TrackModuleLabel") ),
127  fSpacePointModuleLabel (pset.get< std::string >("SpacePointModuleLabel") ),
128  fT0ModuleLabel (pset.get< std::string >("T0ModuleLabel") ),
129  fUseArea(pset.get< bool >("UseArea") ),
130  fFlipTrack_dQdx(pset.get< bool >("FlipTrack_dQdx",true)),
131  caloAlg(pset.get< fhicl::ParameterSet >("CaloAlg"))
132 {
133  produces< std::vector<anab::Calorimetry> >();
134  produces< art::Assns<recob::Track, anab::Calorimetry> >();
135 }
calo::Calorimetry::fT0ModuleLabel
std::string fT0ModuleLabel
Definition: Calorimetry_module.cc:98
calo::Calorimetry::fSpacePointModuleLabel
std::string fSpacePointModuleLabel
Definition: Calorimetry_module.cc:97
calo::Calorimetry::caloAlg
CalorimetryAlg caloAlg
Definition: Calorimetry_module.cc:101
calo::Calorimetry::fTrackModuleLabel
std::string fTrackModuleLabel
Definition: Calorimetry_module.cc:96
calo::Calorimetry::~Calorimetry ( )
virtual

Definition at line 138 of file Calorimetry_module.cc.

139 {
140 
141 }

Member Function Documentation

void calo::Calorimetry::beginJob ( )
virtual

Reimplemented from art::EDProducer.

Definition at line 144 of file Calorimetry_module.cc.

145 {
146  return;
147 }
bool calo::Calorimetry::BeginsOnBoundary ( art::Ptr< recob::Track >  lar_track)
private
template<typename T , BranchType = InEvent>
ProductToken<T> art::Consumer::consumes ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it)
inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

148 {
149  if (!moduleContext_)
150  return ProductToken<T>::invalid();
151 
152  consumables_[BT].emplace_back(ConsumableType::Product,
153  TypeID{typeid(T)},
154  it.label(),
155  it.instance(),
156  it.process());
157  return ProductToken<T>{it};
158 }
art::ProductToken::invalid
static ProductToken< T > invalid()
Definition: ProductToken.h:47
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
template<typename T , art::BranchType BT>
void art::Consumer::consumesMany ( )
inherited

Definition at line 162 of file Consumer.h.

163 {
164  if (!moduleContext_)
165  return;
166 
167  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
168 }
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it)
inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

173 {
174  if (!moduleContext_)
175  return ViewToken<T>::invalid();
176 
177  consumables_[BT].emplace_back(ConsumableType::ViewElement,
178  TypeID{typeid(T)},
179  it.label(),
180  it.instance(),
181  it.process());
182  return ViewToken<T>{it};
183 }
art::ViewToken::invalid
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed)
inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

27 {
28  return rng()->createEngine(placeholder_schedule_id(), seed);
29 }
seed
long seed
Definition: chem4.cc:68
art::EngineCreator::rng
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
Definition: EngineCreator.cc:58
EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make 
)
inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

34 {
35  return rng()->createEngine(
36  placeholder_schedule_id(), seed, kind_of_engine_to_make);
37 }
seed
long seed
Definition: chem4.cc:68
art::EngineCreator::rng
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
Definition: EngineCreator.cc:58
EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make,
label_t const &  engine_label 
)
inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

43 {
44  return rng()->createEngine(
45  placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
46 }
seed
long seed
Definition: chem4.cc:68
art::EngineCreator::rng
static art::ServiceHandle< art::RandomNumberGenerator > & rng()
Definition: EngineCreator.cc:58
CurrentProcessingContext const * art::EDProducer::currentContext ( ) const
protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

121  {
122  return current_context_.get();
123  }
art::EDProducer::current_context_
CPC_exempt_ptr current_context_
Definition: EDProducer.h:116
bool calo::Calorimetry::EndsOnBoundary ( art::Ptr< recob::Track >  lar_track)
private
EngineCreator::seed_t EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset,
char const  key[] = "seed",
seed_t const  implicit_seed = -1 
)
inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

52 {
53  auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
54  return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
55 }
void calo::Calorimetry::GetPitch ( art::Ptr< recob::Hit >  hit,
std::vector< double >  trkx,
std::vector< double >  trky,
std::vector< double >  trkz,
std::vector< double >  trkw,
std::vector< double >  trkx0,
double *  xyz3d,
double &  pitch,
double  TickT0 
)
private

Definition at line 573 of file Calorimetry_module.cc.

References geo::CryostatID::Cryostat, DEFINE_ART_MODULE, recob::Hit::PeakTime(), geo::PlaneID::Plane, geo::GeometryCore::Plane(), geo::PlaneGeo::ThetaZ(), geo::TPCID::TPC, geo::WireID::Wire, recob::Hit::WireID(), and geo::GeometryCore::WirePitch().

Referenced by produce().

573  {
574  //Get 3d coordinates and track pitch for each hit
575  //Find 5 nearest space points and determine xyz and curvature->track pitch
576 
577  //std::cout << "Start of get pitch" << std::endl;
578 
579  // Get services
580  art::ServiceHandle<geo::Geometry> geom;
581  auto const* dp = lar::providerFrom<detinfo::DetectorPropertiesService>();
582 
583  //save distance to each spacepoint sorted by distance
584  std::map<double,size_t> sptmap;
585  //save the sign of distance
586  std::map<size_t, int> sptsignmap;
587 
588  double wire_pitch = geom->WirePitch(0);
589 
590  double t0 = hit->PeakTime() - TickT0;
591  double x0 = dp->ConvertTicksToX(t0, hit->WireID().Plane, hit->WireID().TPC, hit->WireID().Cryostat);
592  double w0 = hit->WireID().Wire;
593 
594  for (size_t i = 0; i<trkx.size(); ++i){
595  double distance = pow((trkw[i]-w0)*wire_pitch,2)+pow(trkx0[i]-x0,2);
596  if (distance>0) distance = sqrt(distance);
597  //std::cout << "Dis " << distance << ", sqaured " << distance*distance << " = " << wire_pitch*wire_pitch <<"("<<trkw[i]<<"-"<<w0<<")^2 + ("<<trkx0[i]<<"-"<<x0<<")^2"<<std::endl;
598  sptmap.insert(std::pair<double,size_t>(distance,i));
599  if (w0-trkw[i]>0) sptsignmap.insert(std::pair<size_t,int>(i,1));
600  else sptsignmap.insert(std::pair<size_t,int>(i,-1));
601  }
602 
603  //x,y,z vs distance
604  std::vector<double> vx;
605  std::vector<double> vy;
606  std::vector<double> vz;
607  std::vector<double> vs;
608 
609  double kx = 0, ky = 0, kz = 0;
610 
611  int np = 0;
612  for (auto isp = sptmap.begin(); isp!=sptmap.end(); isp++){
613 // const double *xyz = new double[3];
614 // xyz = isp->second->XYZ();
615  double xyz[3];
616  xyz[0] = trkx[isp->second];
617  xyz[1] = trky[isp->second];
618  xyz[2] = trkz[isp->second];
619 
620  double distancesign = sptsignmap[isp->second];
621  //std::cout<<np<<" "<<xyz[0]<<" "<<xyz[1]<<" "<<xyz[2]<<" "<<(*isp).first<<std::endl;
622  if (np==0&&isp->first>30){//hit not on track
623  xyz3d[0] = -1000;
624  xyz3d[1] = -1000;
625  xyz3d[2] = -1000;
626  pitch = -1;
627  return;
628  }
629  //std::cout<<np<<" "<<xyz[0]<<" "<<xyz[1]<<" "<<xyz[2]<<" "<<(*isp).first<<" Plane " << hit->WireID().Plane << " TPC " << hit->WireID().TPC << std::endl;
630  if (np<5) {
631  vx.push_back(xyz[0]);
632  vy.push_back(xyz[1]);
633  vz.push_back(xyz[2]);
634  vs.push_back(isp->first*distancesign);
635  }
636  else {
637  break;
638  }
639  np++;
640  //delete [] xyz;
641  }
642  //std::cout<<"np="<<np<<std::endl;
643  if (np>=2){//at least two points
644  //std::cout << "At least 2 points.."<<std::endl;
645  TGraph *xs = new TGraph(np,&vs[0],&vx[0]);
646  //for (int i = 0; i<np; i++) std::cout<<i<<" "<<vs[i]<<" "<<vx[i]<<" "<<vy[i]<<" "<<vz[i]<<std::endl;
647  try{
648  if (np>2){
649  xs->Fit("pol2","Q");
650  }
651  else{
652  xs->Fit("pol1","Q");
653  }
654  TF1 *pol = 0;
655  if (np>2) pol = (TF1*) xs->GetFunction("pol2");
656  else pol = (TF1*) xs->GetFunction("pol1");
657  xyz3d[0] = pol->Eval(0);
658  kx = pol->GetParameter(1);
659  //std::cout<<"X fit "<<xyz3d[0]<<" "<<kx<<std::endl;
660  }
661  catch(...){
662  mf::LogWarning("Calorimetry::GetPitch") <<"Fitter failed";
663  xyz3d[0] = vx[0];
664  }
665  delete xs;
666  TGraph *ys = new TGraph(np,&vs[0],&vy[0]);
667  try{
668  if (np>2){
669  ys->Fit("pol2","Q");
670  }
671  else{
672  ys->Fit("pol1","Q");
673  }
674  TF1 *pol = 0;
675  if (np>2) pol = (TF1*) ys->GetFunction("pol2");
676  else pol = (TF1*) ys->GetFunction("pol1");
677  xyz3d[1] = pol->Eval(0);
678  ky = pol->GetParameter(1);
679  //std::cout<<"Y fit "<<xyz3d[1]<<" "<<ky<<std::endl;
680  }
681  catch(...){
682  mf::LogWarning("Calorimetry::GetPitch") <<"Fitter failed";
683  xyz3d[1] = vy[0];
684  }
685  delete ys;
686  TGraph *zs = new TGraph(np,&vs[0],&vz[0]);
687  try{
688  if (np>2){
689  zs->Fit("pol2","Q");
690  }
691  else{
692  zs->Fit("pol1","Q");
693  }
694  TF1 *pol = 0;
695  if (np>2) pol = (TF1*) zs->GetFunction("pol2");
696  else pol = (TF1*) zs->GetFunction("pol1");
697  xyz3d[2] = pol->Eval(0);
698  kz = pol->GetParameter(1);
699  //std::cout<<"Z fit "<<xyz3d[2]<<" "<<kz<<std::endl;
700  }
701  catch(...){
702  mf::LogWarning("Calorimetry::GetPitch") <<"Fitter failed";
703  xyz3d[2] = vz[0];
704  }
705  delete zs;
706  }
707  else if (np){
708  xyz3d[0] = vx[0];
709  xyz3d[1] = vy[0];
710  xyz3d[2] = vz[0];
711  }
712  else{
713  xyz3d[0] = -1000;
714  xyz3d[1] = -1000;
715  xyz3d[2] = -1000;
716  pitch = -1;
717  return;
718  }
719  pitch = -1;
720  if (kx*kx+ky*ky+kz*kz){
721  double tot = sqrt(kx*kx+ky*ky+kz*kz);
722  kx /= tot;
723  ky /= tot;
724  kz /= tot;
725  //get pitch
726  double wirePitch = geom->WirePitch(hit->WireID().Plane,hit->WireID().TPC,hit->WireID().Cryostat);
727  double angleToVert = geom->Plane(hit->WireID().Plane,hit->WireID().TPC,hit->WireID().Cryostat).Wire(0).ThetaZ(false) - 0.5*TMath::Pi();
728  double cosgamma = TMath::Abs(TMath::Sin(angleToVert)*ky+TMath::Cos(angleToVert)*kz);
729  if (cosgamma>0) pitch = wirePitch/cosgamma;
730 
731  }
732  //std::cout << "At end of get pitch " << xyz3d[0] << " " << xyz3d[1] << " " << xyz3d[2] << " " << x0 << " " << std::endl;
733 }
t0
code to link reconstructed objects back to the MC truth information
Definition: DirectHitParticleAssns_tool.cc:18
geo::GeometryCore::Plane
PlaneGeo const & Plane(unsigned int const p, unsigned int const tpc=0, unsigned int const cstat=0) const
Returns the specified wire.
Definition: GeometryCore.h:2938
recob::Hit::WireID
geo::WireID WireID() const
Initial tdc tick for hit.
Definition: Hit.h:234
geo::CryostatID::Cryostat
CryostatID_t Cryostat
Index of cryostat.
Definition: geo_types.h:130
geo::WireID::Wire
WireID_t Wire
Index of the wire within its plane.
Definition: geo_types.h:313
geo::GeometryCore::WirePitch
geo::Length_t WirePitch(geo::PlaneID const &planeid) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:751
geo::PlaneGeo::ThetaZ
double ThetaZ() const
Angle of the wires from positive z axis; .
Definition: PlaneGeo.cxx:728
geo::PlaneID::Plane
PlaneID_t Plane
Index of the plane within its TPC.
Definition: geo_types.h:258
recob::Hit::PeakTime
float PeakTime() const
Time of the signal peak, in tick units.
Definition: Hit.h:219
mf::LogWarning
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
Definition: MessageLogger.h:213
geo::TPCID::TPC
TPCID_t TPC
Index of the TPC within its cryostat.
Definition: geo_types.h:203
template<typename PROD , BranchType B>
ProductID art::EDProducer::getProductID ( std::string const &  instanceName = {}) const
inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

124  {
125  return ProducerBase::getProductID<PROD, B>(moduleDescription_,
126  instanceName);
127  }
art::EDProducer::moduleDescription_
ModuleDescription moduleDescription_
Definition: EDProducer.h:115
template<typename PROD , BranchType B>
ProductID art::ProducerBase::getProductID ( ModuleDescription const &  moduleDescription,
std::string const &  instanceName 
) const
inherited

Definition at line 56 of file ProducerBase.h.

References B, and art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

58  {
59  auto const& pd =
60  get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
61  return pd.productID();
62  }
B
Int_t B
Definition: plot.C:25
template<typename T , BranchType = InEvent>
ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it)
inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

191 {
192  if (!moduleContext_)
193  return ProductToken<T>::invalid();
194 
195  consumables_[BT].emplace_back(ConsumableType::Product,
196  TypeID{typeid(T)},
197  it.label(),
198  it.instance(),
199  it.process());
200  return ProductToken<T>{it};
201 }
art::ProductToken::invalid
static ProductToken< T > invalid()
Definition: ProductToken.h:47
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
template<typename T , art::BranchType BT>
void art::Consumer::mayConsumeMany ( )
inherited

Definition at line 205 of file Consumer.h.

206 {
207  if (!moduleContext_)
208  return;
209 
210  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
211 }
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it)
inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

216 {
217  if (!moduleContext_)
218  return ViewToken<T>::invalid();
219 
220  consumables_[BT].emplace_back(ConsumableType::ViewElement,
221  TypeID{typeid(T)},
222  it.label(),
223  it.instance(),
224  it.process());
225  return ViewToken<T>{it};
226 }
art::ViewToken::invalid
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
bool art::ProducerBase::modifiesEvent ( ) const
inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

41  {
42  return true;
43  }
cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )
staticinherited

Definition at line 76 of file Consumer.cc.

Referenced by art::RootOutput::beginSubRun(), art::OutputModule::doBeginRun(), art::OutputModule::doBeginSubRun(), art::OutputModule::doEndRun(), art::OutputModule::doEndSubRun(), art::ProducingService::doPostReadEvent(), art::ProducingService::doPostReadRun(), art::ProducingService::doPostReadSubRun(), art::OutputModule::doWriteEvent(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EventProcessor::readEvent(), art::EventProcessor::readRun(), art::EmptyEvent::readRun_(), art::EventProcessor::readSubRun(), and art::EmptyEvent::readSubRun_().

77 {
78  static Consumer invalid{InvalidTag{}};
79  return &invalid;
80 }
void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset)
protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

90 {
91  if (!moduleContext_)
92  return;
93 
94  pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
95  for (auto& consumablesPerBranch : consumables_) {
96  cet::sort_all(consumablesPerBranch);
97  }
98 }
art::Consumer::requireConsumes_
bool requireConsumes_
Definition: Consumer.h:137
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
void calo::Calorimetry::produce ( art::Event &  evt)
virtual

Implements art::EDProducer.

Definition at line 150 of file Calorimetry_module.cc.

References util::absDiff(), caloAlg, util::CreateAssn(), geo::CryostatID::Cryostat, art::FindManyP< ProdB, Data >::data(), calo::CalorimetryAlg::dEdx_AMP(), calo::CalorimetryAlg::dEdx_AREA(), dir, e, fdEdx, fdQdx, fetime, fFlipTrack_dQdx, art::fill_ptr_vector(), geo::GeometryCore::FindTPCAtPosition(), fMIPs, fnsps, fpitch, fResRng, fSpacePointModuleLabel, fstime, fT0ModuleLabel, ftime, fTrackModuleLabel, fUseArea, fwire, fXYZ, art::DataViewImpl::getByLabel(), GetPitch(), hits(), geo::CryostatID::isValid, util::kBogusD, LOG_DEBUG, geo::GeometryCore::Nplanes(), geo::PlaneID::Plane, geo::GeometryCore::Plane(), geo::GeometryCore::PlaneWireToChannel(), art::Event::put(), geo::TPCID::TPC, lar::util::TrackPitchInView(), geo::PlaneGeo::View(), w, geo::GeometryCore::WireAngleToVertical(), geo::GeometryCore::WirePitch(), x, xx, and zz.

151 {
152  auto const* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
153 
154  art::Handle< std::vector<recob::Track> > trackListHandle;
155  std::vector<art::Ptr<recob::Track> > tracklist;
156  if (evt.getByLabel(fTrackModuleLabel,trackListHandle))
157  art::fill_ptr_vector(tracklist, trackListHandle);
158 
159  // Get Geometry
160  art::ServiceHandle<geo::Geometry> geom;
161 
162  // channel quality
163  lariov::ChannelStatusProvider const& channelStatus
164  = art::ServiceHandle<lariov::ChannelStatusService>()->GetProvider();
165 
166  size_t nplanes = geom->Nplanes();
167 
168  //create anab::Calorimetry objects and make association with recob::Track
169  std::unique_ptr< std::vector<anab::Calorimetry> > calorimetrycol(new std::vector<anab::Calorimetry>);
170  std::unique_ptr< art::Assns<recob::Track, anab::Calorimetry> > assn(new art::Assns<recob::Track, anab::Calorimetry>);
171 
172  //art::FindManyP<recob::SpacePoint> fmsp(trackListHandle, evt, fTrackModuleLabel);
173  art::FindManyP<recob::Hit> fmht(trackListHandle, evt, fTrackModuleLabel);
174  art::FindManyP<recob::Hit, recob::TrackHitMeta> fmthm(trackListHandle, evt, fTrackModuleLabel); //this has more information about hit-track association, only available in PMA for now
175  art::FindManyP<anab::T0> fmt0(trackListHandle, evt, fT0ModuleLabel);
176 
177  for(size_t trkIter = 0; trkIter < tracklist.size(); ++trkIter){
178 
179  decltype(auto) larEnd = tracklist[trkIter]->Trajectory().End();
180  //store track directional cosines
181  double trackCosStart[3]={0.,0.,0.};
182  double trackCosEnd[3]={0.,0.,0.};
183  tracklist[trkIter]->Direction(trackCosStart,trackCosEnd);
184 
185  // Some variables for the hit
186  float time; //hit time at maximum
187  float stime; //hit start time
188  float etime; //hit end time
189  uint32_t channel = 0;//channel number
190  unsigned int cstat = 0; //hit cryostat number
191  unsigned int tpc = 0; //hit tpc number
192  unsigned int wire = 0; //hit wire number
193  unsigned int plane = 0; //hit plane number
194 
195  std::vector< art::Ptr<recob::Hit> > allHits = fmht.at(trkIter);
196  double T0 =0;
197  double TickT0 =0;
198  if ( fmt0.isValid() ) {
199  std::vector< art::Ptr<anab::T0> > allT0 = fmt0.at(trkIter);
200  if ( allT0.size() ) T0 = allT0[0]->Time();
201  TickT0 = T0 / detprop->SamplingRate();
202  }
203 
204  std::vector< std::vector<unsigned int> > hits(nplanes);
205 
206  art::FindManyP<recob::SpacePoint> fmspts(allHits, evt, fSpacePointModuleLabel);
207  for (size_t ah = 0; ah< allHits.size(); ++ah){
208  hits[allHits[ah]->WireID().Plane].push_back(ah);
209  }
210  //get hits in each plane
211  for (size_t ipl = 0; ipl < nplanes; ++ipl){//loop over all wire planes
212 
213  geo::PlaneID planeID;//(cstat,tpc,ipl);
214 
215  fwire.clear();
216  ftime.clear();
217  fstime.clear();
218  fetime.clear();
219  fMIPs.clear();
220  fdQdx.clear();
221  fdEdx.clear();
222  fpitch.clear();
223  fResRng.clear();
224  fXYZ.clear();
225 
226  double Kin_En = 0.;
227  double Trk_Length = 0.;
228  std::vector<double> vdEdx;
229  std::vector<double> vresRange;
230  std::vector<double> vdQdx;
231  std::vector<double> deadwire; //residual range for dead wires
232  std::vector<TVector3> vXYZ;
233 
234  //range of wire signals
235  unsigned int wire0 = 100000;
236  unsigned int wire1 = 0;
237  double PIDA = 0;
238  int nPIDA = 0;
239 
240  // determine track direction. Fill residual range array
241  bool GoingDS = true;
242  // find the track direction by comparing US and DS charge BB
243  double USChg = 0;
244  double DSChg = 0;
245  // temp array holding distance betweeen space points
246  std::vector<double> spdelta;
247  //int nht = 0; //number of hits
248  fnsps = 0; //number of space points
249  std::vector<double> ChargeBeg;
250  std::stack<double> ChargeEnd;
251 
252  // find track pitch
253  double fTrkPitch = 0;
254  for (size_t itp = 0; itp < tracklist[trkIter]->NumberTrajectoryPoints(); ++itp){
255  const TVector3& pos = tracklist[trkIter]->LocationAtPoint(itp);
256  const double Position[3] = { pos.X(), pos.Y(), pos.Z() };
257  geo::TPCID tpcid = geom->FindTPCAtPosition ( Position );
258  if (tpcid.isValid) {
259  try{
260  fTrkPitch = lar::util::TrackPitchInView(*tracklist[trkIter], geom->Plane(ipl).View(), itp);
261  }
262  catch( cet::exception &e){
263  mf::LogWarning("Calorimetry") << "caught exception "
264  << e << "\n setting pitch (C) to "
265  << util::kBogusD;
266  fTrkPitch = 0;
267  }
268  break;
269  }
270  }
271 
272  // find the separation between all space points
273  double xx = 0.,yy = 0.,zz = 0.;
274 
275  //save track 3d points
276  std::vector<double> trkx;
277  std::vector<double> trky;
278  std::vector<double> trkz;
279  std::vector<double> trkw;
280  std::vector<double> trkx0;
281  for (size_t i = 0; i<hits[ipl].size(); ++i){
282  //Get space points associated with the hit
283  std::vector< art::Ptr<recob::SpacePoint> > sptv = fmspts.at(hits[ipl][i]);
284  for (size_t j = 0; j < sptv.size(); ++j){
285 
286  double t = allHits[hits[ipl][i]]->PeakTime() - TickT0; // Want T0 here? Otherwise ticks to x is wrong?
287  double x = detprop->ConvertTicksToX(t, allHits[hits[ipl][i]]->WireID().Plane, allHits[hits[ipl][i]]->WireID().TPC, allHits[hits[ipl][i]]->WireID().Cryostat);
288  double w = allHits[hits[ipl][i]]->WireID().Wire;
289  if (TickT0){
290  trkx.push_back(sptv[j]->XYZ()[0]-detprop->ConvertTicksToX(TickT0, allHits[hits[ipl][i]]->WireID().Plane, allHits[hits[ipl][i]]->WireID().TPC, allHits[hits[ipl][i]]->WireID().Cryostat));
291  }
292  else{
293  trkx.push_back(sptv[j]->XYZ()[0]);
294  }
295  trky.push_back(sptv[j]->XYZ()[1]);
296  trkz.push_back(sptv[j]->XYZ()[2]);
297  trkw.push_back(w);
298  trkx0.push_back(x);
299  }
300  }
301  for (size_t ihit = 0; ihit < hits[ipl].size(); ++ihit){//loop over all hits on each wire plane
302 
303  //std::cout<<ihit<<std::endl;
304 
305  if (!planeID.isValid){
306  plane = allHits[hits[ipl][ihit]]->WireID().Plane;
307  tpc = allHits[hits[ipl][ihit]]->WireID().TPC;
308  cstat = allHits[hits[ipl][ihit]]->WireID().Cryostat;
309  planeID.Cryostat = cstat;
310  planeID.TPC = tpc;
311  planeID.Plane = plane;
312  planeID.isValid = true;
313  }
314 
315  wire = allHits[hits[ipl][ihit]]->WireID().Wire;
316  time = allHits[hits[ipl][ihit]]->PeakTime(); // What about here? T0
317  stime = allHits[hits[ipl][ihit]]->PeakTimeMinusRMS();
318  etime = allHits[hits[ipl][ihit]]->PeakTimePlusRMS();
319 
320  double charge = allHits[hits[ipl][ihit]]->PeakAmplitude();
321  if (fUseArea) charge = allHits[hits[ipl][ihit]]->Integral();
322  //get 3d coordinate and track pitch for the current hit
323  //not all hits are associated with space points, the method uses neighboring spacepts to interpolate
324  double xyz3d[3];
325  double pitch;
326  bool fBadhit = false;
327  if (fmthm.isValid()){
328  auto vhit = fmthm.at(trkIter);
329  auto vmeta = fmthm.data(trkIter);
330  for (size_t ii = 0; ii<vhit.size(); ++ii){
331  if (vhit[ii].key() == allHits[hits[ipl][ihit]].key()){
332  if (vmeta[ii]->Index()>=tracklist[trkIter]->NumberTrajectoryPoints()){
333  throw cet::exception("Calorimetry_module.cc") << "Requested track trajectory index "<<vmeta[ii]->Index()<<" exceeds the total number of trajectory points "<<tracklist[trkIter]->NumberTrajectoryPoints()<<" for track index "<<trkIter<<". Something is wrong with the track reconstruction. Please contact tjyang@fnal.gov";
334  }
335  if (!tracklist[trkIter]->HasValidPoint(vmeta[ii]->Index())){
336  fBadhit = true;
337  continue;
338  }
339  double angleToVert = geom->WireAngleToVertical(vhit[ii]->View(), vhit[ii]->WireID().TPC, vhit[ii]->WireID().Cryostat) - 0.5*::util::pi<>();
340  const TVector3& dir = tracklist[trkIter]->DirectionAtPoint(vmeta[ii]->Index());
341  double cosgamma = std::abs(std::sin(angleToVert)*dir.Y() + std::cos(angleToVert)*dir.Z());
342  if (cosgamma){
343  pitch = geom->WirePitch(0)/cosgamma;
344  }
345  else{
346  pitch = 0;
347  }
348  TVector3 loc = tracklist[trkIter]->LocationAtPoint(vmeta[ii]->Index());
349  xyz3d[0] = loc.X();
350  xyz3d[1] = loc.Y();
351  xyz3d[2] = loc.Z();
352  break;
353  }
354  }
355  }
356  else
357  GetPitch(allHits[hits[ipl][ihit]], trkx, trky, trkz, trkw, trkx0, xyz3d, pitch, TickT0);
358 
359  if (fBadhit) continue;
360  if (xyz3d[2]<-100) continue; //hit not on track
361  if (pitch<=0) pitch = fTrkPitch;
362  if (!pitch) continue;
363 
364  if(fnsps == 0) {
365  xx = xyz3d[0];
366  yy = xyz3d[1];
367  zz = xyz3d[2];
368  spdelta.push_back(0);
369  } else {
370  double dx = xyz3d[0] - xx;
371  double dy = xyz3d[1] - yy;
372  double dz = xyz3d[2] - zz;
373  spdelta.push_back(sqrt(dx*dx + dy*dy + dz*dz));
374  Trk_Length += spdelta.back();
375  xx = xyz3d[0];
376  yy = xyz3d[1];
377  zz = xyz3d[2];
378  }
379 
380  ChargeBeg.push_back(charge);
381  ChargeEnd.push(charge);
382 
383  double MIPs = charge;
384  double dQdx = MIPs/pitch;
385  double dEdx = 0;
386  if (fUseArea) dEdx = caloAlg.dEdx_AREA(allHits[hits[ipl][ihit]], pitch, T0);
387  else dEdx = caloAlg.dEdx_AMP(allHits[hits[ipl][ihit]], pitch, T0);
388 
389  Kin_En = Kin_En + dEdx * pitch;
390 
391  if (allHits[hits[ipl][ihit]]->WireID().Wire < wire0) wire0 = allHits[hits[ipl][ihit]]->WireID().Wire;
392  if (allHits[hits[ipl][ihit]]->WireID().Wire > wire1) wire1 = allHits[hits[ipl][ihit]]->WireID().Wire;
393 
394  fMIPs.push_back(MIPs);
395  fdEdx.push_back(dEdx);
396  fdQdx.push_back(dQdx);
397  fwire.push_back(wire);
398  ftime.push_back(time);
399  fstime.push_back(stime);
400  fetime.push_back(etime);
401  fpitch.push_back(pitch);
402  TVector3 v(xyz3d[0],xyz3d[1],xyz3d[2]);
403  //std::cout << "Adding these positions to v and then fXYZ " << xyz3d[0] << " " << xyz3d[1] << " " << xyz3d[2] << "\n" <<std::endl;
404  fXYZ.push_back(v);
405  ++fnsps;
406  }
407  if (!fnsps){
408  //std::cout << "Adding the aforementioned positions..." << std::endl;
409  calorimetrycol->push_back(anab::Calorimetry(util::kBogusD,
410  vdEdx,
411  vdQdx,
412  vresRange,
413  deadwire,
414  util::kBogusD,
415  fpitch,
416  vXYZ,
417  planeID));
418  util::CreateAssn(*this, evt, *calorimetrycol, tracklist[trkIter], *assn);
419  continue;
420  }
421  for (int isp = 0; isp<fnsps; ++isp){
422  if (isp>3) break;
423  USChg += ChargeBeg[isp];
424  }
425  int countsp = 0;
426  while (!ChargeEnd.empty()){
427  if (countsp>3) break;
428  DSChg += ChargeEnd.top();
429  ChargeEnd.pop();
430  ++countsp;
431  }
432  // Going DS if charge is higher at the end
433  GoingDS = (DSChg > USChg) || (!fFlipTrack_dQdx);
434  // determine the starting residual range and fill the array
435  fResRng.resize(fnsps);
436  if(GoingDS) {
437  fResRng[fnsps - 1] = spdelta[fnsps - 1] / 2;
438  for(int isp = fnsps - 2; isp > -1; isp--) {
439  fResRng[isp] = fResRng[isp+1] + spdelta[isp+1];
440  }
441  } else {
442  fResRng[0] = spdelta[1] / 2;
443  for(int isp = 1; isp < fnsps; isp++) {
444  fResRng[isp] = fResRng[isp-1] + spdelta[isp];
445  }
446  }
447 
448  LOG_DEBUG("CaloPrtHit") << " pt wire time ResRng MIPs pitch dE/dx Ai X Y Z\n";
449 
450  double Ai = -1;
451  for (int i = 0; i < fnsps; ++i){//loop over all 3D points
452  vresRange.push_back(fResRng[i]);
453  vdEdx.push_back(fdEdx[i]);
454  vdQdx.push_back(fdQdx[i]);
455  vXYZ.push_back(fXYZ[i]);
456  if (i!=0 && i!= fnsps-1){//ignore the first and last point
457  // Calculate PIDA
458  Ai = fdEdx[i] * pow(fResRng[i],0.42);
459  nPIDA++;
460  PIDA += Ai;
461  }
462  LOG_DEBUG("CaloPrtHit") <<std::setw(4)<< trkIter
463  //std::cout<<std::setw(4)<< trkIter
464  <<std::setw(4)<< ipl
465  <<std::setw(4) << i
466  <<std::setw(4) << fwire[i]
467  << std::setw(6) << (int)ftime[i]
468  << std::setiosflags(std::ios::fixed | std::ios::showpoint)
469  << std::setprecision(2)
470  << std::setw(8) << fResRng[i]
471  << std::setprecision(1)
472  << std::setw(8) << fMIPs[i]
473  << std::setprecision(2)
474  << std::setw(8) << fpitch[i]
475  << std::setw(8) << fdEdx[i]
476  << std::setw(8) << Ai
477  << std::setw(8) << fXYZ[i].x()
478  << std::setw(8) << fXYZ[i].y()
479  << std::setw(8) << fXYZ[i].z()
480  << "\n";
481  }//end looping over 3D points
482  if(nPIDA > 0) {
483  PIDA = PIDA / (double)nPIDA;
484  }
485  else {
486  PIDA = -1;
487  }
488  LOG_DEBUG("CaloPrtTrk") << "Plane # "<< ipl
489  << "TrkPitch= "
490  << std::setprecision(2) << fTrkPitch
491  << " nhits= " << fnsps
492  << "\n"
493  << std::setiosflags(std::ios::fixed | std::ios::showpoint)
494  << "Trk Length= " << std::setprecision(1)
495  << Trk_Length << " cm,"
496  << " KE calo= " << std::setprecision(1)
497  << Kin_En << " MeV,"
498  << " PIDA= " << PIDA
499  << "\n";
500 
501  // look for dead wires
502  for (unsigned int iw = wire0; iw<wire1+1; ++iw){
503  plane = allHits[hits[ipl][0]]->WireID().Plane;
504  tpc = allHits[hits[ipl][0]]->WireID().TPC;
505  cstat = allHits[hits[ipl][0]]->WireID().Cryostat;
506  channel = geom->PlaneWireToChannel(plane,iw,tpc,cstat);
507  if (channelStatus.IsBad(channel)){
508  LOG_DEBUG("Calorimetry") << "Found dead wire at Plane = " << plane
509  << " Wire =" << iw;
510  unsigned int closestwire = 0;
511  unsigned int endwire = 0;
512  unsigned int dwire = 100000;
513  double mindis = 100000;
514  double goodresrange = 0;
515  //hitCtr = 0;
516  for (size_t ihit = 0; ihit <hits[ipl].size(); ++ihit){
517  // for(art::PtrVector<recob::Hit>::const_iterator hitIter = hitsV.begin();
518  // hitIter != hitsV.end();
519  // ++hitCtr, hitIter++){
520  channel = allHits[hits[ipl][ihit]]->Channel();
521  if (channelStatus.IsBad(channel)) continue;
522  // grab the space points associated with this hit
523  std::vector< art::Ptr<recob::SpacePoint> > sppv = fmspts.at(hits[ipl][ihit]);
524  if(sppv.size() < 1) continue;
525  // only use the first space point in the collection, really each hit should
526  // only map to 1 space point
527  const recob::Track::Point_t xyz{sppv[0]->XYZ()[0],
528  sppv[0]->XYZ()[1],
529  sppv[0]->XYZ()[2]};
530  double dis1 = (larEnd - xyz).Mag2();
531  if (dis1) dis1 = std::sqrt(dis1);
532  if (dis1 < mindis){
533  endwire = allHits[hits[ipl][ihit]]->WireID().Wire;
534  mindis = dis1;
535  }
536  if (util::absDiff(wire, iw) < dwire){
537  closestwire = allHits[hits[ipl][ihit]]->WireID().Wire;
538  dwire = util::absDiff(allHits[hits[ipl][ihit]]->WireID().Wire, iw);
539  goodresrange = dis1;
540  }
541  }
542  if (closestwire){
543  if (iw < endwire){
544  deadwire.push_back(goodresrange+(int(closestwire)-int(iw))*fTrkPitch);
545  }
546  else{
547  deadwire.push_back(goodresrange+(int(iw)-int(closestwire))*fTrkPitch);
548  }
549  }
550  }
551  }
552  //std::cout << "Adding at the end but still same fXYZ" << std::endl;
553  calorimetrycol->push_back(anab::Calorimetry(Kin_En,
554  vdEdx,
555  vdQdx,
556  vresRange,
557  deadwire,
558  Trk_Length,
559  fpitch,
560  vXYZ,
561  planeID));
562  util::CreateAssn(*this, evt, *calorimetrycol, tracklist[trkIter], *assn);
563 
564  }//end looping over planes
565  }//end looping over tracks
566 
567  evt.put(std::move(calorimetrycol));
568  evt.put(std::move(assn));
569 
570  return;
571 }
x
Float_t x
Definition: compare.C:6
calo::Calorimetry::fdEdx
std::vector< double > fdEdx
Definition: Calorimetry_module.cc:110
xx
Double_t xx
Definition: macro.C:12
geo::GeometryCore::Plane
PlaneGeo const & Plane(unsigned int const p, unsigned int const tpc=0, unsigned int const cstat=0) const
Returns the specified wire.
Definition: GeometryCore.h:2938
calo::CalorimetryAlg::dEdx_AREA
double dEdx_AREA(art::Ptr< recob::Hit > hit, double pitch, double T0=0) const
Definition: CalorimetryAlg.cxx:88
geo::PlaneID
The data type to uniquely identify a Plane.
Definition: geo_types.h:250
geo::CryostatID::isValid
bool isValid
Whether this ID points to a valid element.
Definition: geo_types.h:129
calo::Calorimetry::fpitch
std::vector< double > fpitch
Definition: Calorimetry_module.cc:112
geo::CryostatID::Cryostat
CryostatID_t Cryostat
Index of cryostat.
Definition: geo_types.h:130
zz
Double_t zz
Definition: plot.C:279
calo::Calorimetry::fXYZ
std::vector< TVector3 > fXYZ
Definition: Calorimetry_module.cc:113
art::End
typename BeginEndPackage< L >::End End
Definition: PrincipalPackages.h:259
geo::GeometryCore::WirePitch
geo::Length_t WirePitch(geo::PlaneID const &planeid) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:751
calo::Calorimetry::fResRng
std::vector< double > fResRng
Definition: Calorimetry_module.cc:111
art::Event::put
ProductID put(std::unique_ptr< PROD > &&product)
Definition: Event.h:102
geo::GeometryCore::FindTPCAtPosition
geo::TPCID FindTPCAtPosition(double const worldLoc[3]) const
Returns the ID of the TPC at specified location.
Definition: GeometryCore.h:2417
geo::PlaneGeo::View
View_t View() const
Which coordinate does this plane measure.
Definition: PlaneGeo.h:171
geo::GeometryCore::Nplanes
unsigned int Nplanes(unsigned int tpc=0, unsigned int cstat=0) const
Returns the total number of wire planes in the specified TPC.
Definition: GeometryCore.h:2792
hits
void hits()
Definition: readHits.C:15
calo::Calorimetry::fT0ModuleLabel
std::string fT0ModuleLabel
Definition: Calorimetry_module.cc:98
calo::Calorimetry::fSpacePointModuleLabel
std::string fSpacePointModuleLabel
Definition: Calorimetry_module.cc:97
calo::Calorimetry::ftime
std::vector< double > ftime
Definition: Calorimetry_module.cc:105
calo::Calorimetry::fetime
std::vector< double > fetime
Definition: Calorimetry_module.cc:107
calo::Calorimetry::GetPitch
void GetPitch(art::Ptr< recob::Hit > hit, std::vector< double > trkx, std::vector< double > trky, std::vector< double > trkz, std::vector< double > trkw, std::vector< double > trkx0, double *xyz3d, double &pitch, double TickT0)
Definition: Calorimetry_module.cc:573
util::CreateAssn
bool CreateAssn(PRODUCER const &prod, art::Event &evt, std::vector< T > const &a, art::Ptr< U > const &b, art::Assns< U, T > &assn, std::string a_instance, size_t indx=UINT_MAX)
Creates a single one-to-one association.
Definition: AssociationUtil.h:591
geo::TPCID
The data type to uniquely identify a TPC.
Definition: geo_types.h:195
geo::PlaneID::Plane
PlaneID_t Plane
Index of the plane within its TPC.
Definition: geo_types.h:258
calo::Calorimetry::fMIPs
std::vector< double > fMIPs
Definition: Calorimetry_module.cc:108
util::absDiff
constexpr auto absDiff(A const &a, B const &b)
Returns the absolute value of the difference between two values.
Definition: NumericUtils.h:43
calo::Calorimetry::caloAlg
CalorimetryAlg caloAlg
Definition: Calorimetry_module.cc:101
geo::GeometryCore::PlaneWireToChannel
raw::ChannelID_t PlaneWireToChannel(WireID const &wireid) const
Returns the ID of the TPC channel connected to the specified wire.
Definition: GeometryCore.cxx:1167
dir
TDirectory * dir
Definition: macro.C:5
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
Definition: DataViewImpl.h:344
calo::Calorimetry::fdQdx
std::vector< double > fdQdx
Definition: Calorimetry_module.cc:109
mf::LogWarning
MaybeLogger_< ELseverityLevel::ELsev_warning, false > LogWarning
Definition: MessageLogger.h:213
recob::Track::Point_t
tracking::Point_t Point_t
Definition: Track.h:55
calo::Calorimetry::fwire
std::vector< int > fwire
Definition: Calorimetry_module.cc:104
LOG_DEBUG
#define LOG_DEBUG(id)
Definition: MessageLogger.h:263
calo::Calorimetry::fTrackModuleLabel
std::string fTrackModuleLabel
Definition: Calorimetry_module.cc:96
util::kBogusD
constexpr double kBogusD
obviously bogus double value
Definition: PhysicalConstants.h:71
lar::util::TrackPitchInView
double TrackPitchInView(recob::Track const &track, geo::View_t view, size_t trajectory_point=0)
Provides projected wire pitch for the view.
Definition: TrackUtils.cxx:76
calo::Calorimetry::fstime
std::vector< double > fstime
Definition: Calorimetry_module.cc:106
art::fill_ptr_vector
void fill_ptr_vector(std::vector< Ptr< T >> &ptrs, H const &h)
Definition: Ptr.h:464
geo::TPCID::TPC
TPCID_t TPC
Index of the TPC within its cryostat.
Definition: geo_types.h:203
e
Float_t e
Definition: plot.C:34
trkf::Plane
recob::tracking::Plane Plane
Definition: TrackState.h:17
w
Float_t w
Definition: plot.C:23
geo::GeometryCore::WireAngleToVertical
double WireAngleToVertical(geo::View_t view, geo::TPCID const &tpcid) const
Returns the angle of the wires in the specified view from vertical.
Definition: GeometryCore.cxx:770
calo::CalorimetryAlg::dEdx_AMP
double dEdx_AMP(art::Ptr< recob::Hit > hit, double pitch, double T0=0) const
Definition: CalorimetryAlg.cxx:55
fhicl::exception
cet::coded_exception< error, detail::translate > exception
Definition: exception.h:33
void calo::Calorimetry::ReadCaloTree ( )
private
void art::Consumer::showMissingConsumes ( ) const
protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

126 {
127  if (!moduleContext_)
128  return;
129 
130  // If none of the branches have missing consumes statements, exit early.
131  if (std::all_of(cbegin(missingConsumes_),
132  cend(missingConsumes_),
133  [](auto const& perBranch) { return perBranch.empty(); }))
134  return;
135 
136  constexpr cet::HorizontalRule rule{60};
137  mf::LogPrint log{"MTdiagnostics"};
138  log << '\n'
139  << rule('=') << '\n'
140  << "The following consumes (or mayConsume) statements are missing from\n"
141  << module_context(moduleDescription_) << '\n'
142  << rule('-') << '\n';
143 
144  cet::for_all_with_index(
145  missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
146  for (auto const& pi : perBranch) {
147  log << " "
148  << assemble_consumes_statement(static_cast<BranchType>(i), pi)
149  << '\n';
150  }
151  });
152  log << rule('=');
153 }
mf::MaybeLogger_
Definition: MessageLogger.h:114
art::Consumer::moduleDescription_
cet::exempt_ptr< ModuleDescription const > moduleDescription_
Definition: Consumer.h:140
util::pi
constexpr T pi()
Returns the constant pi (up to 35 decimal digits of precision)
Definition: PhysicalConstants.h:80
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
art::Consumer::missingConsumes_
ConsumableProductSets missingConsumes_
Definition: Consumer.h:139
void art::Consumer::validateConsumedProduct ( BranchType const  bt,
ProductInfo const &  pi 
)
protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

103 {
104  // Early exits if consumes tracking has been disabled or if the
105  // consumed product is an allowed consumable.
106  if (!moduleContext_)
107  return;
108 
109  if (cet::binary_search_all(consumables_[bt], pi))
110  return;
111 
112  if (requireConsumes_) {
113  throw Exception(errors::ProductRegistrationFailure,
114  "Consumer: an error occurred during validation of a "
115  "retrieved product\n\n")
116  << "The following consumes (or mayConsume) statement is missing from\n"
117  << module_context(moduleDescription_) << ":\n\n"
118  << " " << assemble_consumes_statement(bt, pi) << "\n\n";
119  }
120 
121  missingConsumes_[bt].insert(pi);
122 }
art::Consumer::moduleDescription_
cet::exempt_ptr< ModuleDescription const > moduleDescription_
Definition: Consumer.h:140
art::Consumer::requireConsumes_
bool requireConsumes_
Definition: Consumer.h:137
art::Exception
cet::coded_exception< errors::ErrorCodes, ExceptionDetail::translate > Exception
Definition: Exception.h:66
util::pi
constexpr T pi()
Returns the constant pi (up to 35 decimal digits of precision)
Definition: PhysicalConstants.h:80
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:138
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:136
art::Consumer::missingConsumes_
ConsumableProductSets missingConsumes_
Definition: Consumer.h:139

Member Data Documentation

CalorimetryAlg calo::Calorimetry::caloAlg
private

Definition at line 101 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<double> calo::Calorimetry::fdEdx
private

Definition at line 110 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<double> calo::Calorimetry::fdQdx
private

Definition at line 109 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<double> calo::Calorimetry::fetime
private

Definition at line 107 of file Calorimetry_module.cc.

Referenced by produce().

bool calo::Calorimetry::fFlipTrack_dQdx
private

Definition at line 100 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<double> calo::Calorimetry::fMIPs
private

Definition at line 108 of file Calorimetry_module.cc.

Referenced by produce().

int calo::Calorimetry::fnsps
private

Definition at line 103 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<double> calo::Calorimetry::fpitch
private

Definition at line 112 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<double> calo::Calorimetry::fResRng
private

Definition at line 111 of file Calorimetry_module.cc.

Referenced by produce().

std::string calo::Calorimetry::fSpacePointModuleLabel
private

Definition at line 97 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<double> calo::Calorimetry::fstime
private

Definition at line 106 of file Calorimetry_module.cc.

Referenced by produce().

std::string calo::Calorimetry::fT0ModuleLabel
private

Definition at line 98 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<double> calo::Calorimetry::ftime
private

Definition at line 105 of file Calorimetry_module.cc.

Referenced by produce().

std::string calo::Calorimetry::fTrackModuleLabel
private

Definition at line 96 of file Calorimetry_module.cc.

Referenced by produce().

bool calo::Calorimetry::fUseArea
private

Definition at line 99 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<int> calo::Calorimetry::fwire
private

Definition at line 104 of file Calorimetry_module.cc.

Referenced by produce().

std::vector<TVector3> calo::Calorimetry::fXYZ
private

Definition at line 113 of file Calorimetry_module.cc.

Referenced by produce().

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