Inheritance diagram for calo::GeneralCalorimetry:

Public Member Functions
	GeneralCalorimetry (fhicl::ParameterSet const &pset)

Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDProducer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept

	Producer (fhicl::ParameterSet const &)

	Producer (Producer const &)=delete

	Producer (Producer &&)=delete

Producer &	operator= (Producer const &)=delete

Producer &	operator= (Producer &&)=delete

void	doBeginJob ()

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::Modifier
	~Modifier () noexcept

	Modifier ()

	Modifier (Modifier const &)=delete

	Modifier (Modifier &&)=delete

Modifier &	operator= (Modifier const &)=delete

Modifier &	operator= (Modifier &&)=delete

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 Member Functions
void	produce (art::Event &evt)

Private Attributes
std::string	fTrackModuleLabel
	module creating the track objects and assns to hits More...

double	fADCToElectrons
	filled using the detinfo::DetectorPropertiesService service More...

geo::View_t	fCollectionView
	view of the collection plane More...

unsigned int	fCollectionPlane
	plane of the collection plane More...

art::ServiceHandle< geo::Geometry const >	fGeo

CalorimetryAlg	caloAlg

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >

Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)

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 30 of file GeneralCalorimetry_module.cc.

Constructor & Destructor Documentation

calo::GeneralCalorimetry::GeneralCalorimetry ( fhicl::ParameterSet const & pset )

explicit

Definition at line 54 of file GeneralCalorimetry_module.cc.

   : EDProducer{pset}
   , fCollectionView(geo::kUnknown)
   , fCollectionPlane(0)
   , caloAlg(pset.get< fhicl::ParameterSet >("CaloAlg"))
 {
   auto const* dp = lar::providerFrom<detinfo::DetectorPropertiesService>();
   fADCToElectrons = 1./dp->ElectronsToADC();
 
   // determine the view of the collection plane
   // just look at one cryostat, the first TPC and loop over those
   // planes
   geo::TPCID FirstTPC(0, 0);
   for(unsigned int p = 0; p < fGeo->Nplanes(FirstTPC); ++p){
     geo::PlaneID planeid{ FirstTPC, p };
     if(fGeo->SignalType(planeid) == geo::kCollection){
       fCollectionView = fGeo->View(planeid);
       fCollectionPlane = p;
     }
   }
 
   fTrackModuleLabel = pset.get< std::string >("TrackModuleLabel");
 
   produces< std::vector<anab::Calorimetry>              >();
   produces< art::Assns<recob::Track, anab::Calorimetry> >();
 }

Member Function Documentation

void calo::GeneralCalorimetry::produce ( art::Event & evt )

privatevirtual

Implements art::EDProducer.

Definition at line 82 of file GeneralCalorimetry_module.cc.

 {
   /*************************************************************/
   /*                          WARNING                          */
   /*************************************************************/
   /* The dQdx information in recob::Track has been deprecated  */
   /* since 2016 and in 11/2018 the recob::Track interface was  */
   /* changed and DQdxAtPoint and NumberdQdx were removed.      */
   /* Therefore the code below is now commented out             */
   /* (note that it was most likely not functional anyways).    */
   /* For any issue please contact: larsoft-team@fnal.gov       */
   /*************************************************************/
   /*
   art::Handle< std::vector<recob::Track> > trackHandle;
   evt.getByLabel(fTrackModuleLabel, trackHandle);
   std::vector< art::Ptr<recob::Track> > tracks;
   art::fill_ptr_vector(tracks, trackHandle);
 
   auto const* dp = lar::providerFrom<detinfo::DetectorPropertiesService>();
 
   //create anab::Calorimetry objects and make association with recob::Track
   std::unique_ptr< std::vector<anab::Calorimetry> > calorimetrycol(new std::vector<anab::Calorimetry>);
   std::unique_ptr< art::Assns<recob::Track, anab::Calorimetry> > assn(new art::Assns<recob::Track, anab::Calorimetry>);
 
   art::FindMany<recob::Hit> fmht(trackHandle, evt, fTrackModuleLabel);
 
   // loop over the tracks
   for(size_t t = 0; t < tracks.size(); ++t){
 
     art::Ptr<recob::Track> trk = tracks.at(t);
 
     double viewPitch = 0.;
     try{
       viewPitch = lar::util::TrackPitchInView(*trk, fCollectionView);
     }
     catch( cet::exception &e){
       mf::LogWarning("GeneralCalorimetry") << "caught exception "
                                            << e
                                            << "\n pitch now set to 0";
     }
 
     // loop over the track trajectory to get the kinetic energy,
     // residual range and the dQdx
     float kineticEnergy = 0.;
     std::vector<float> vdEdx;
     std::vector<float> vresRange;
     std::vector<float> vdQdx;
     std::vector<float> deadwire; //residual range for dead wires
 
     // Check that the number of trajectory points and number of dQdx for the
     // collection view are the same
     if( trk->NumberTrajectoryPoints() != trk->NumberdQdx(fCollectionView) )
       throw cet::exception("GeneralCalorimetry") << "inconsistent number of track trajectory "
                                                  << " and dQdx points\n";
 
     if (trk->NumberTrajectoryPoints()>2){
       for(size_t p = 1; p < trk->NumberTrajectoryPoints()-1; ++p){
         if (!trk->DQdxAtPoint(p, fCollectionView)) continue;
         vresRange.push_back(trk->Length(p));
         vdQdx.push_back(trk->DQdxAtPoint(p, fCollectionView));
         //vdEdx.push_back(vdQdx.back() * fADCToElectrons/util::kGeVToElectrons*1000);
         vdEdx.push_back(caloAlg.dEdx_AMP(vdQdx.back(),
                                          dp->ConvertXToTicks(trk->LocationAtPoint(p).X(),fCollectionPlane,0,0),
                                          fCollectionPlane));
         kineticEnergy += vdEdx.back(); // \todo should this be converted from electrons to energy?
         std::cout<<vresRange.back()<<" "<<vdQdx.back()<<" "<<vdEdx.back()<<std::endl;
       }
 
       geo::PlaneID planeID(0,0,fCollectionPlane);
       calorimetrycol->push_back(anab::Calorimetry(kineticEnergy,
                                                   vdEdx,
                                                   vdQdx,
                                                   vresRange,
                                                   deadwire,
                                                   trk->Length(),
                                                   viewPitch,
                                                   planeID));
 
       util::CreateAssn(*this, evt, *calorimetrycol, trk, *assn);
 
     }
   }// end of loop over all tracks
 
   if (calorimetrycol->size()){
     evt.put(std::move(calorimetrycol));
     evt.put(std::move(assn));
   }
 
   */
   /*************************************************************/
   return;
 }