Inheritance diagram for shwf::ShowerReco:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

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

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

virtual	~ShowerReco ()

void	beginJob ()

void	beginRun (art::Run &run)

void	reconfigure (fhicl::ParameterSet const &pset)

void	produce (art::Event &evt)

void	GetVertexAndAnglesFromCluster (art::Ptr< recob::Cluster > clust, unsigned int plane)

void	Get2DVariables (art::PtrVector< recob::Hit > hitlist)

void	LongTransEnergy (unsigned int set, std::vector< art::Ptr< recob::Hit > > hitlist, bool isData=false)

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	ClearandResizeVectors (unsigned int nPlanes)

Private Attributes
int	fRun

int	fEvent

int	fSubRun

float	slope [3]

float	angle [3]

std::string	fClusterModuleLabel

float	ftimetick

double	xyz_vertex [3]

double	totCnrg

double	totCnrg_corr

double	fMean_wire_pitch

fhicl::ParameterSet	fCaloPSet

std::vector< double >	fRMS_2cm

std::vector< int >	fNpoints_2cm

std::vector< double >	fCorr_MeV_2cm

std::vector< double >	fCorr_Charge_2cm

std::vector< int >	fNpoints_corr_ADC_2cm

std::vector< int >	fNpoints_corr_MeV_2cm

std::vector< double >	fTotChargeADC

std::vector< double >	fTotChargeMeV

std::vector< double >	fTotChargeMeV_MIPs

std::vector< double >	fChargeADC_2cm

std::vector< double >	fChargeMeV_2cm

std::vector< double >	fChargeMeV_2cm_refined

std::vector< double >	fChargeMeV_2cm_axsum

std::vector< std::vector< double > >	fDistribChargeADC

std::vector< std::vector< double > >	fDistribChargeMeV

std::vector< std::vector< double > >	fDistribHalfChargeMeV

std::vector< std::vector< double > >	fDistribChargeposition

std::vector< std::vector< double > >	fSingleEvtAngle

std::vector< std::vector< double > >	fSingleEvtAngleVal

std::vector< unsigned int >	fWire_vertex

std::vector< double >	fTime_vertex

std::vector< double >	fWire_vertexError

std::vector< double >	fTime_vertexError

std::vector< unsigned int >	fWire_last

std::vector< double >	fTime_last

std::vector< double >	xyz_vertex_fit

std::vector< std::vector< double > >	fNPitch

float	Kin_En

std::vector< float >	vdEdx

std::vector< float >	vresRange

std::vector< float >	vdQdx

std::vector< float >	deadwire

float	Trk_Length

float	fTrkPitchC

float	fdEdxlength

float	fcalodEdxlength

bool	fUseArea

double	xphi

double	xtheta

unsigned int	fTPC

unsigned int	fNPlanes

unsigned int	fNAngles

TTree *	ftree_shwf

double	fWirePitch

double	fTimeTick

double	fDriftVelocity

double	fWireTimetoCmCm

std::vector< int >	fNhitsperplane

std::vector< double >	fTotADCperplane

art::ServiceHandle< geo::Geometry >	geom

const detinfo::DetectorProperties *	detprop = lar::providerFrom<detinfo::DetectorPropertiesService>()

Detailed Description

Definition at line 79 of file ShowerReco_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

shwf::ShowerReco::ShowerReco ( fhicl::ParameterSet const & pset )

explicit

METHODS global

Definition at line 212 of file ShowerReco_module.cc.

References reconfigure().

 {
   this->reconfigure(pset);
   produces< std::vector<recob::Shower>                >();
   produces< art::Assns<recob::Shower, recob::Cluster> >();
   produces< art::Assns<recob::Shower, recob::Hit>     >();
   produces< std::vector<anab::Calorimetry>              >();
   produces< art::Assns<recob::Shower, anab::Calorimetry> >();
 }

shwf::ShowerReco::~ShowerReco ( )

virtual

Constructor

Definition at line 237 of file ShowerReco_module.cc.

238 {

239 }

Member Function Documentation

void shwf::ShowerReco::beginJob ( )

virtual

Destructor

Get Geometry

Todo:: the call to geo->Nplanes() assumes this is a single cryostat and single TPC detector

Todo:: need to generalize to multiple cryostats and TPCs

Get TFileService and define output Histograms

All-knowing tree with reconstruction information

Reimplemented from art::EDProducer.

Definition at line 250 of file ShowerReco_module.cc.

References detprop, fChargeADC_2cm, fChargeMeV_2cm, fChargeMeV_2cm_axsum, fChargeMeV_2cm_refined, fDistribChargeADC, fDistribChargeMeV, fDistribChargeposition, fDistribHalfChargeMeV, fEvent, fMean_wire_pitch, fNAngles, fNhitsperplane, fNPitch, fNPlanes, fNpoints_2cm, fRMS_2cm, fRun, fSubRun, ftimetick, fTotADCperplane, fTotChargeADC, fTotChargeMeV, fTotChargeMeV_MIPs, ftree_shwf, art::TFileDirectory::make(), geo::GeometryCore::Nplanes(), detinfo::DetectorProperties::SamplingRate(), geo::GeometryCore::WirePitch(), xphi, xtheta, and xyz_vertex_fit.

 {
 
 
 
   art::ServiceHandle<geo::Geometry> geo;
 
   detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
 
   fNPlanes = geo->Nplanes();
   fMean_wire_pitch = geo->WirePitch(); //wire pitch in cm
 
   art::ServiceHandle<art::TFileService> tfs;
 
  // auto const* detp = lar::providerFrom<detinfo::DetectorPropertiesService>();
   ftimetick=detprop->SamplingRate()/1000.;
   
 
   
   
   
 
   ftree_shwf =tfs->make<TTree>("ShowerReco","Results");
   // ftree_shwf->Branch("ftheta_Mean","std::vector<double>",&fTheta_Mean);
   // ftree_shwf->Branch("ftheta_RMS","std::vector<double>",&fTheta_RMS );
   
   ftree_shwf->Branch("run",&fRun,"run/I");
   ftree_shwf->Branch("subrun",&fSubRun,"subrun/I");
   ftree_shwf->Branch("event",&fEvent,"event/I");
   ftree_shwf->Branch("nplanes",&fNPlanes,"nplanes/I");
   ftree_shwf->Branch("nangles",&fNAngles,"nangles/I");
   
 //   ftree_shwf->Branch("fthetaN","std::vector<double>",&fThetaN_ang);
 //   ftree_shwf->Branch("fphiN","std::vector<double>",&fPhiN_ang);
   
   ftree_shwf->Branch("xtheta",&xtheta,"xtheta/D");
   ftree_shwf->Branch("xphi",&xphi,"xphi/D");
   
   ftree_shwf->Branch("ftotChargeADC","std::vector<double>",&fTotChargeADC);
   ftree_shwf->Branch("ftotChargeMeV","std::vector<double>",&fTotChargeMeV);
   ftree_shwf->Branch("fTotChargeMeV_MIPs","std::vector<double>",&fTotChargeMeV_MIPs);
   
   ftree_shwf->Branch("NPitch","std::vector< std::vector<double> >", &fNPitch);
   //  ftree_shwf->Branch("Pitch","std::vector<double>", &fPitch);
   
   // this should be temporary - until the omega is sorted out.
   ftree_shwf->Branch("RMS_2cm","std::vector<double>",&fRMS_2cm);
   ftree_shwf->Branch("Npoints_2cm","std::vector<int>",&fNpoints_2cm);
 //    ftree_shwf->Branch("RMS_4cm","std::vector<double>",&fRMS_4cm);
 //    ftree_shwf->Branch("Npoints_4cm","std::vector<int>",&fNpoints_4cm);
 
   ftree_shwf->Branch("ChargeADC_2cm","std::vector<double>",&fChargeADC_2cm);
   ftree_shwf->Branch("ChargeMeV_2cm","std::vector<double>",&fChargeMeV_2cm);
 //    ftree_shwf->Branch("ChargeADC_4cm","std::vector<double>",&fChargeADC_4cm);
 //    ftree_shwf->Branch("ChargeMeV_4cm","std::vector<double>",&fChargeMeV_4cm);
    
   ftree_shwf->Branch("ChargeMeV_2cm_refined","std::vector<double>",&fChargeMeV_2cm_refined);
 //    ftree_shwf->Branch("ChargeMeV_4cm_refined","std::vector<double>",&fChargeMeV_4cm_refined);
    
   ftree_shwf->Branch("ChargeMeV_2cm_axsum","std::vector<double>",&fChargeMeV_2cm_axsum);
 //    ftree_shwf->Branch("ChargeMeV_4cm_axsum","std::vector<double>",&fChargeMeV_4cm_axsum);
 
 
   ftree_shwf->Branch("fNhitsperplane","std::vector<int>",&fNhitsperplane);
   ftree_shwf->Branch("fTotADCperplane","std::vector<double>",&fTotADCperplane);
 
   
 
   ftree_shwf->Branch("ChargedistributionADC","std::vector<std::vector<double>>",&fDistribChargeADC);
   
   ftree_shwf->Branch("ChargedistributionMeV","std::vector<std::vector<double>>",&fDistribChargeMeV);
   ftree_shwf->Branch("DistribHalfChargeMeV","std::vector<std::vector<double>>",&fDistribHalfChargeMeV);
   ftree_shwf->Branch("ChargedistributionPosition","std::vector<std::vector<double>>",&fDistribChargeposition);
   ftree_shwf->Branch("xyz_vertex_fit","std::vector<double>", &xyz_vertex_fit);
  
 }

void shwf::ShowerReco::beginRun ( art::Run & run )

virtual

Reimplemented from art::EDProducer.

Definition at line 333 of file ShowerReco_module.cc.

References detprop, detinfo::DetectorProperties::DriftVelocity(), detinfo::DetectorProperties::Efield(), fChargeADC_2cm, fChargeMeV_2cm, fChargeMeV_2cm_axsum, fChargeMeV_2cm_refined, fCorr_Charge_2cm, fCorr_MeV_2cm, fDistribChargeADC, fDistribChargeMeV, fDistribChargeposition, fDistribHalfChargeMeV, fDriftVelocity, fNAngles, fNhitsperplane, fNPitch, fNPlanes, fNpoints_2cm, fNpoints_corr_ADC_2cm, fNpoints_corr_MeV_2cm, fRMS_2cm, fTime_last, fTime_vertex, fTime_vertexError, fTimeTick, fTotADCperplane, fTotChargeADC, fTotChargeMeV, fTotChargeMeV_MIPs, fWire_last, fWire_vertex, fWire_vertexError, fWirePitch, fWireTimetoCmCm, geom, detinfo::DetectorProperties::SamplingRate(), detinfo::DetectorProperties::Temperature(), vdEdx, vdQdx, vresRange, geo::GeometryCore::WirePitch(), and xyz_vertex_fit.

 {
  detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
 
  fWirePitch = geom->WirePitch(); //wire pitch in cm
  fTimeTick=detprop->SamplingRate()/1000.;  
  fDriftVelocity=detprop->DriftVelocity(detprop->Efield(),detprop->Temperature());
  fWireTimetoCmCm=(fTimeTick*fDriftVelocity)/fWirePitch; 
    
 }

void shwf::ShowerReco::ClearandResizeVectors ( unsigned int nPlanes )

private

Referenced by produce().

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().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

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().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

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().

 {
   return rng()->createEngine(placeholder_schedule_id(), seed);
 }

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().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make);
 }

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().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
 }

CurrentProcessingContext const * art::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

   {
     return current_context_.get();
   }

void shwf::ShowerReco::Get2DVariables ( art::PtrVector< recob::Hit > hitlist )

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_().

 {
   auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
   return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
 }

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_.

   {
     return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                                instanceName);
   }

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().

   {
     auto const& pd =
       get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
     return pd.productID();
   }

void shwf::ShowerReco::GetVertexAndAnglesFromCluster	(	art::Ptr< recob::Cluster >	clust,
		unsigned int	plane
	)

Actual routine that reconstruct the shower

Definition at line 1089 of file ShowerReco_module.cc.

References angle, DEFINE_ART_MODULE, recob::Cluster::EndTick(), recob::Cluster::EndWire(), fTime_last, fTime_vertex, fTime_vertexError, fWire_last, fWire_vertex, fWire_vertexError, recob::Cluster::SigmaStartTick(), recob::Cluster::SigmaStartWire(), slope, recob::Cluster::StartAngle(), recob::Cluster::StartTick(), and recob::Cluster::StartWire().

Referenced by produce().

 {
   //convert to cm/cm units needed in the calculation
 //  slope[plane]=clust->dTdW();//*(ftimetick*fdriftvelocity)/fMean_wire_pitch;  
   angle[plane]=clust->StartAngle();
   slope[plane]=std::tan(clust->StartAngle());
   fWire_vertex[plane]=clust->StartWire();
   fTime_vertex[plane]=clust->StartTick();
 
   fWire_vertexError[plane]=clust->SigmaStartWire();  // wire coordinate of vertex for each plane
   fTime_vertexError[plane]=clust->SigmaStartTick();  // time coordinate of vertex for each plane
   
   fWire_last[plane]=clust->EndWire();  // wire coordinate of last point for each plane
   fTime_last[plane]=clust->EndTick();
   
 
  // std::cout << "======= setting slope for view: " << plane 
 //                          << " " << slope[plane] << " " << fWire_vertex[plane] 
 //                          << " " << fTime_vertex[plane] << " " << std::endl;
         //                  <<  fWire_vertex[plane]+50<< " "
         //                  << fTime_vertex[plane] + slope[plane]*(fWire_vertex[plane]+50)<< std::endl;
   
 }

void shwf::ShowerReco::LongTransEnergy	(	unsigned int	set,
		std::vector< art::Ptr< recob::Hit > >	hitlist,
		bool	isData = `false`
	)

third loop to get only points inside of 1RMS of value.

Definition at line 847 of file ShowerReco_module.cc.

References calo::CalorimetryAlg::dEdx_AMP(), calo::CalorimetryAlg::dEdx_AREA(), fcalodEdxlength, fCaloPSet, fChargeMeV_2cm, fChargeMeV_2cm_refined, fCorr_MeV_2cm, fdEdxlength, fDistribChargeMeV, fDistribChargeposition, fNPitch, fNPlanes, fNpoints_2cm, fNpoints_corr_MeV_2cm, fRMS_2cm, fTime_vertex, fTotChargeADC, fTotChargeMeV, fTrkPitchC, fUseArea, fWire_vertex, fWireTimetoCmCm, geom, util::GeometryUtilities::Get2DDistance(), util::GeometryUtilities::GetPointOnLine(), geo::kCollection, Kin_En, geo::kMysteryType, recob::Hit::PeakTime(), util::GeometryUtilities::PitchInView(), geo::PlaneID::Plane, geo::GeometryCore::SignalType(), slope, totCnrg, totCnrg_corr, Trk_Length, vdEdx, vdQdx, vresRange, geo::WireID::Wire, recob::Hit::WireID(), xphi, and xtheta.

Referenced by produce().

 {
   // alogorithm for energy vs dx of the shower (roto-translation) COLLECTION VIEW
   // double  wire_cm, time_cm;
   // int loop_nrg = 0;
 
   calo::CalorimetryAlg calalg(fCaloPSet);
   util::GeometryUtilities gser;
   
   
   
   double totCnrg = 0,totCnrg_corr =0;//, totNewCnrg=0 ; // tot enegry of the shower in collection
 //   art::ServiceHandle<geo::Geometry> geom;
 //   auto const* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
 
   double time;
   unsigned int wire=0,plane=fNPlanes-1;
 
   
   double mevav2cm=0.;  
   double sum=0.;
   double npoints_calo=0;
 
   int direction=-1;
   
   //override direction if phi (XZ angle) is less than 90 degrees
   if(fabs(xphi)<90)
        direction=1;
  
   //variables to check whether a hit is close to the shower axis.
   double ortdist,linedist;
   double wire_on_line,time_on_line;
   
   //get effective pitch using 3D angles
   double newpitch=gser.PitchInView(plane,xphi,xtheta);
  
   
   
   for(art::PtrVector<recob::Hit>::const_iterator hitIter = hitlist.begin(); hitIter != hitlist.end();  hitIter++){
     art::Ptr<recob::Hit> theHit = (*hitIter);
     time = theHit->PeakTime() ;  
   
     wire=  theHit->WireID().Wire;
     plane=  theHit->WireID().Plane;
     
     double dEdx_new;
    // double dEdx_MIP;
     
     if(fUseArea)
     { dEdx_new = calalg.dEdx_AREA((*hitIter), newpitch ); 
      // dEdx_MIP = calalg.dEdx_AREA_forceMIP((*hitIter), newpitch ); 
     }
     else  //this will hopefully go away, once all of the calibration factors are calculated.
     {
       dEdx_new = calalg.dEdx_AMP((*hitIter), newpitch ); 
      // dEdx_MIP = calalg.dEdx_AMP_forceMIP((*hitIter), newpitch );
     }
    
     //calculate total energy.
     totCnrg_corr += dEdx_new; 
    //totNewCnrg+=dEdx_MIP;
 
     // calculate the wire,time coordinates of the hit projection on to the 2D shower axis
     gser.GetPointOnLine(slope[plane]/fWireTimetoCmCm,fWire_vertex[plane],fTime_vertex[plane],wire,time,wire_on_line,time_on_line);
     linedist=gser.Get2DDistance(wire_on_line,time_on_line,fWire_vertex[plane],fTime_vertex[plane]);
     ortdist=gser.Get2DDistance(wire_on_line,time_on_line,wire,time);
     
  
     
    //calculate the distance from the vertex using the effective pitch metric 
     double wdist=(((double)wire-(double)fWire_vertex[plane])*newpitch)*direction;  //wdist is always positive
 
     
 //   if( (fabs(wdist)<fcalodEdxlength)&&(fabs(wdist)>0.2)){  
   if( (wdist<fcalodEdxlength)&&(wdist>0.2)){    
   
     vdEdx.push_back(dEdx_new);
     vresRange.push_back(fabs(wdist));
     vdQdx.push_back((*hitIter)->PeakAmplitude()/newpitch);
     Trk_Length=wdist;
     fTrkPitchC=fNPitch[set][plane];
     Kin_En+=dEdx_new*newpitch;
     npoints_calo++;
     sum+=dEdx_new;
     
     
     
 //std::cout << " CALORIMETRY:" << " Pitch " <<newpitch << " dist: " << wdist <<  " dE/dx: " << dEdx_new << "MeV/cm "  << " average: " <<  sum/npoints_calo  << "hit: wire, time " << wire << " " << time << " line,ort " << linedist << " " << ortdist<< " direction " << direction << std::endl;
     
    
      if(wdist<fdEdxlength 
       && ((direction==1 && wire>fWire_vertex[plane])     //take no hits before vertex (depending on direction)
       || (direction==-1 && wire<fWire_vertex[plane])  ) 
              && ortdist<4.5 && linedist < fdEdxlength ){
               fChargeMeV_2cm[set]+= dEdx_new ; 
               fNpoints_2cm[set]++;
              // std::cout << " CALORIMETRY:" << " Pitch " <<newpitch << " dist: " << wdist <<  " dE/dx: " << dEdx_new << "MeV/cm "  << " average: " <<  sum/npoints_calo  << "hit: wire, time " << wire << " " << time << " line,ort " << linedist << " " << ortdist<< " direction " << direction << std::endl;
              }
 
         // fill out for 4cm preshower
                  
         //fDistribChargeADC[set].push_back(ch_adc);  //vector with the first De/Dx points
         fDistribChargeMeV[set].push_back(dEdx_new);  //vector with the first De/Dx points
         //fDistribHalfChargeMeV[set].push_back(Bcorr_half);
         fDistribChargeposition[set].push_back(wdist);  //vector with the first De/Dx points' positions 
       
      }//end inside range if statement
         
   }// end first loop on hits.
   
   auto const signalType
     = hitlist.empty()? geo::kMysteryType: geom->SignalType(hitlist.front()->WireID());
  
   if(signalType == geo::kCollection)
         {
       fTotChargeADC[set]=totCnrg*newpitch; 
       fTotChargeMeV[set]=totCnrg_corr*newpitch;  
       //fTotChargeMeV_MIPs[set]=totNewCnrg*newpitch;
         }
       
 
   //calculate average dE/dx
   if(fNpoints_2cm[set]>0)       {
     mevav2cm=fChargeMeV_2cm[set]/fNpoints_2cm[set];
    }
    //double RMS_ADC_2cm=0.;
   
   
   
   //second loop to calculate RMS      
   for(art::PtrVector<recob::Hit>::const_iterator hitIter = hitlist.begin(); hitIter != hitlist.end(); hitIter++){
     art::Ptr<recob::Hit> theHit = (*hitIter);
     time = theHit->PeakTime() ;  
     wire=  theHit->WireID().Wire;
     plane=  theHit->WireID().Plane;
     double dEdx=0;
      
    if(fUseArea)
     { dEdx = calalg.dEdx_AREA((*hitIter), newpitch ); 
      }
     else  //this will hopefully go away, once all of the calibration factors are calculated.
     {
       dEdx = calalg.dEdx_AMP((*hitIter), newpitch ); 
      }
     
     
     
     gser.GetPointOnLine(slope[plane]/fWireTimetoCmCm,fWire_vertex[plane],fTime_vertex[plane],wire,time,wire_on_line,time_on_line);
     linedist=gser.Get2DDistance(wire_on_line,time_on_line,fWire_vertex[plane],fTime_vertex[plane]);
     ortdist=gser.Get2DDistance(wire_on_line,time_on_line,wire,time);
     
     
     
     double wdist=(((double)wire-(double)fWire_vertex[plane])*newpitch)*direction;
     
  //    //std::cout << dEdx << " MeV, outside of if;; wd " << wdist << " ld " << linedist << " od " << ortdist << std::endl;
     
     if( (wdist<fcalodEdxlength)&&(wdist>0.2)){ 
       if(wdist<fdEdxlength
           && ((direction==1 && wire>fWire_vertex[plane]) || 
           (direction==-1 && wire<fWire_vertex[plane])  ) 
              && ortdist<4.5 && linedist < fdEdxlength)
         {
 //        //std::cout << dEdx << " MeV " << std::endl;
         fRMS_2cm[set]+= (dEdx-mevav2cm)*(dEdx-mevav2cm); 
         }
         
     
       }  //end if on correct hits. 
   }  //end RMS_calculating loop.    
   
   if(fNpoints_2cm[set]>0)       
     {
     fRMS_2cm[set]=TMath::Sqrt(fRMS_2cm[set]/fNpoints_2cm[set]);
      }
   
   //std::cout << " average dE/dx: " << mevav2cm << " RMS::  " << fRMS_2cm[set] << " " << fNpoints_2cm[set] <<  std::endl;
   
     
   for(art::PtrVector<recob::Hit>::const_iterator hitIter = hitlist.begin(); hitIter != hitlist.end(); hitIter++){
     art::Ptr<recob::Hit> theHit = (*hitIter);
     time = theHit->PeakTime() ;  
     wire=  theHit->WireID().Wire;
     plane=  theHit->WireID().Plane;
    
     double dEdx=0;
      if(fUseArea)
     { dEdx = calalg.dEdx_AREA((*hitIter), newpitch ); 
      }
     else  //this will hopefully go away, once all of the calibration factors are calculated.
     {
       dEdx = calalg.dEdx_AMP((*hitIter), newpitch ); 
      }
     
     gser.GetPointOnLine(slope[plane]/fWireTimetoCmCm,fWire_vertex[plane],fTime_vertex[plane],wire,time,wire_on_line,time_on_line);
     linedist=gser.Get2DDistance(wire_on_line,time_on_line,fWire_vertex[plane],fTime_vertex[plane]);
     ortdist=gser.Get2DDistance(wire_on_line,time_on_line,wire,time);
     
     double wdist=(((double)wire-(double)fWire_vertex[plane])*newpitch)*direction;
     
     
     if( (wdist < fcalodEdxlength) && (wdist > 0.2
         && ((direction==1 && wire>fWire_vertex[plane]) ||
         (direction==-1 && wire<fWire_vertex[plane])  ) 
              && ortdist<4.5 && linedist < fdEdxlength ))
       { 
       if(wdist < fdEdxlength)
           {
           if( ((dEdx > (mevav2cm-fRMS_2cm[set]) )
             && (dEdx < (mevav2cm+fRMS_2cm[set]) )) 
                 || (newpitch > 0.3*fdEdxlength ) ) {
             fCorr_MeV_2cm[set]+= dEdx; 
           fNpoints_corr_MeV_2cm[set]++;
           }
         
       } // end if on good hits
       
       
     }
   } //end of third loop on hits  
   
   if(fNpoints_corr_MeV_2cm[set]>0){
         //std::cout << " ++ NPoints 2cm, ADC and MeV " 
         //                        << fNpoints_corr_MeV_2cm[set] << " " 
         //                        << fNpoints_corr_ADC_2cm[set] 
         //                        << " corrected average De/Dx, charge, MeV:  " 
         //                        << fCorr_Charge_2cm[set]/fNpoints_corr_ADC_2cm[set] 
         //                        << " " << fCorr_MeV_2cm[set]/fNpoints_corr_MeV_2cm[set] << std::endl;
         //fCorr_Charge_2cm[set]/=fNpoints_corr_ADC_2cm[set];
         fCorr_MeV_2cm[set]/=fNpoints_corr_MeV_2cm[set];
         fChargeMeV_2cm_refined[set]=fCorr_MeV_2cm[set];
   }
   
 
 // std::cout << " total ENERGY, birks: " << fTotChargeMeV[set] << " MeV "  << " |average:  " << fChargeMeV_2cm_refined[set] <<   std::endl;
 }

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().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

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().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

     {
       return true;
     }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

 {
   static Consumer invalid{InvalidTag{}};
   return &invalid;
 }

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().

 {
   if (!moduleContext_)
     return;
 
   pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
   for (auto& consumablesPerBranch : consumables_) {
     cet::sort_all(consumablesPerBranch);
   }
 }

void shwf::ShowerReco::produce ( art::Event & evt )

virtual

Get Clusters

Todo:: really need to determine the values of the arguments of the recob::Shower ctor

Fill the output tree with all information

Implements art::EDProducer.

Definition at line 464 of file ShowerReco_module.cc.

References angle, geo::GeometryCore::ChannelsIntersect(), ClearandResizeVectors(), util::CreateAssn(), deadwire, detprop, detinfo::DetectorProperties::DriftVelocity(), e, detinfo::DetectorProperties::Efield(), art::EventID::event(), fChargeADC_2cm, fChargeMeV_2cm, fClusterModuleLabel, fCorr_Charge_2cm, fCorr_MeV_2cm, fDriftVelocity, fEvent, fNAngles, fNhitsperplane, fNPlanes, fNpoints_2cm, fNpoints_corr_ADC_2cm, fNpoints_corr_MeV_2cm, fRMS_2cm, fRun, fSubRun, fTime_vertex, fTime_vertexError, fTimeTick, fTotADCperplane, fTotChargeADC, fTotChargeMeV, fTotChargeMeV_MIPs, ftree_shwf, fTrkPitchC, fWire_last, fWire_vertex, fWire_vertexError, geom, art::DataViewImpl::getByLabel(), GetVertexAndAnglesFromCluster(), hits(), art::Event::id(), util::kBogusD, art::Ptr< T >::key(), Kin_En, geo::PlaneGeo::LocalToWorld(), LongTransEnergy(), geo::GeometryCore::NearestWire(), geo::GeometryCore::Nplanes(), geo::origin(), geo::GeometryCore::Plane(), geo::GeometryCore::PlaneWireToChannel(), art::Event::productGetter(), art::PtrVector< T >::push_back(), art::Event::put(), art::EventID::run(), detinfo::DetectorProperties::SamplingRate(), recob::Shower::set_direction(), recob::Shower::set_direction_err(), art::PtrVector< T >::size(), art::EventID::subRun(), detinfo::DetectorProperties::Temperature(), detinfo::DetectorProperties::TriggerOffset(), Trk_Length, vdEdx, vdQdx, vresRange, xphi, xtheta, xx, xyz_vertex_fit, y, and z.

 { 
 
   detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
 
   util::GeometryUtilities gser;
   fNPlanes = geom->Nplanes();
   //fdriftvelocity = detprop->DriftVelocity(Efield_SI,Temperature);
    std::unique_ptr<std::vector<recob::Shower> > Shower3DVector(new std::vector<recob::Shower>);
    std::unique_ptr< art::Assns<recob::Shower, recob::Cluster> > cassn(new art::Assns<recob::Shower, recob::Cluster>);
    std::unique_ptr< art::Assns<recob::Shower, recob::Hit>     > hassn(new art::Assns<recob::Shower, recob::Hit>);
    std::unique_ptr< std::vector<anab::Calorimetry> > calorimetrycol(new std::vector<anab::Calorimetry>);
    std::unique_ptr< art::Assns< anab::Calorimetry,recob::Shower> > calassn(new art::Assns<anab::Calorimetry,recob::Shower>);
   
  
   art::Handle< std::vector<recob::Cluster> > clusterListHandle;
   evt.getByLabel(fClusterModuleLabel,clusterListHandle);
 
   art::Handle< std::vector<art::PtrVector < recob::Cluster> > > clusterAssociationHandle;
   evt.getByLabel(fClusterModuleLabel,clusterAssociationHandle);
   
   
   art::FindManyP<recob::Hit> fmh(clusterListHandle, evt,fClusterModuleLabel);
   
 //   std::cout << " ---------- ShowerReco!!! -------------- " << std::endl;
   fRun = evt.id().run();
   fSubRun = evt.id().subRun();
   fEvent = evt.id().event();
 //   unsigned int nCollections= clusterAssociationHandle->size();
 //   std::vector < art::PtrVector<recob::Cluster> >::const_iterator clusterSet = clusterAssociationHandle->begin();
 //    for(unsigned int iCol=0;iCol<nCollections;iCol++)
 //   {
 //    const art::PtrVector<recob::Cluster> pvcluster(*(clusterSet++));
 //     auto it = pvcluster.begin();
 //     int nClusts = pvcluster.size();
 //     ClearandResizeVectors(nClusts);   // need to do this smarter (coutn at start and have an overall index?)
 //     for(int iClust = 0; iClust < nClusts; ++iClust) {
 //       const art::Ptr<recob::Cluster> pclust(*(it++)); 
 //       auto pcoll { pclust };
 //       art::FindManyP<recob::Hit> fs( pcoll, evt, fClusterModuleLabel);
 //       std::vector< art::Ptr<recob::Hit> > hitlist = fs.at(0);
 //       //std::cout << " hitlist size for coll " << iCol << " clust " << iClust << " " << hitlist.size() << std::endl;
 //     }
 //   } 
   
   
   // find all the hits associated to all the clusters (once and for all);
   // the index of the query matches the index of the cluster in the collection
   // (conveniently carried around in its art pointer)
   art::FindManyP<recob::Hit> ClusterHits(clusterListHandle, evt, fClusterModuleLabel);
   
   std::vector < art::PtrVector<recob::Cluster> >::const_iterator clusterSet = clusterAssociationHandle->begin();
   // loop over vector of vectors (each size of NPlanes) and reconstruct showers from each of those
   for(size_t iClustSet = 0;iClustSet < clusterAssociationHandle->size(); iClustSet++){
  
     const art::PtrVector<recob::Cluster>  CurrentClusters=(*(clusterSet++));
     
     // do some error checking - i.e. are the clusters themselves present.
     if(clusterListHandle->size() < 2 || CurrentClusters.size() < 2) {
       //std::cout << "not enough clusters to reconstruct" << std::endl;
       //std::cout << "emergency filling tree @ run, evt," <<  fRun << " " << fEvent << std::endl;
       ftree_shwf->Fill();
       return;
     }
     //std::cout << " Cluster Set: " << iClustSet << " " << std::endl; 
     
     ClearandResizeVectors( fNPlanes);
     
     std::vector< std::vector< art::Ptr<recob::Hit> > > hitlist_all;
     hitlist_all.resize(fNPlanes);
   
    // int nClusts = CurrentClusters.size();
     for(size_t iClust = 0; iClust < CurrentClusters.size(); iClust++){
       art::Ptr<recob::Cluster> const& pclust = CurrentClusters[iClust];
         //size_t ii=0; 
         //std::cout << " clusterListHandle  " << clusterListHandle->size() << " fNPlanes " << fNPlanes << " "<< CurrentClusters.size() << std::endl;
       
       // get all the hits for this cluster;
       // pclust is a art::Ptr to the original cluster collection stored in the event;
       // its key corresponds to its index in the collection
       // (and therefore to the index in the query)
       std::vector< art::Ptr<recob::Hit> > const& hitlist = ClusterHits.at(pclust.key());
       //std::cout << " hitlist size for coll " << iClustSet << " clust " << iClust << " " << hitlist.size() << std::endl;
     
       unsigned int p(0); //c=channel, p=plane, w=wire
     
       //std::cout << " hitlist size " << hitlist.size() << std::endl;
       if(hitlist.size() == 0) continue;
       
       p=  (*hitlist.begin())->WireID().Plane;
       //art::Ptr<recob::Cluster> cl(clusterListHandle, iClust);
       //get vertex position and slope information to start with - ii is the posistion of the correct cluster:
       GetVertexAndAnglesFromCluster( pclust,p);
     
       
       double ADCcharge=0;
       //loop over cluster hits
       for(art::Ptr<recob::Hit> const& hit: hitlist){
         p=  hit->WireID().Plane;
         hitlist_all[p].push_back(hit);
         ADCcharge+= hit->PeakAmplitude();
       }
       fNhitsperplane[p]=hitlist_all[p].size();
       fTotADCperplane[p]=ADCcharge;
     //   unsigned int nCollections= clusterAssociationHandle->size();
 //   std::vector < art::PtrVector<recob::Cluster> >::const_iterator clusterSet = clusterAssociationHandle->begin();
 //    for(unsigned int iCol=0;iCol<nCollections;iCol++)
 //   {
 //    const art::PtrVector<recob::Cluster> pvcluster(*(clusterSet++));
 //     auto it = pvcluster.begin();
 //     int nClusts = pvcluster.size();
 //     ClearandResizeVectors(nClusts);   // need to do this smarter (coutn at start and have an overall index?)
 //     for(int iClust = 0; iClust < nClusts; ++iClust) {
 //       const art::Ptr<recob::Cluster> pclust(*(it++)); 
 //       auto pcoll { pclust };
 //       art::FindManyP<recob::Hit> fs( pcoll, evt, fClusterModuleLabel);
 //       std::vector< art::Ptr<recob::Hit> > hitlist = fs.at(0);
 //       //std::cout << " hitlist size for coll " << iCol << " clust " << iClust << " " << hitlist.size() << std::endl;
 //     }
 //   } 
     
     
  
     
     
    // Find the right cluster from the standard cluster list -> to get the hitlist associations. 
 //    for( ii = 0; ii < clusterListHandle->size(); ii++){
 //      art::Ptr<recob::Cluster> cl(clusterListHandle, ii);
 //      if((*cl).ID() == (*CurrentClusters[iClust]).ID() )  //find the right cluster out of the list of associated clusters
 //      break;
 //        }
 //     
     //get vertex position and slope information to start with - ii is the posistion of the correct cluster:
    // std::vector< art::Ptr<recob::Hit> > hitlist = fmh.at(ii);
   //  std::sort(hitlist.begin(), hitlist.end(), SortByWire());
   
 
   } // End loop on clusters.
   // Now I have the Hitlists and the relevent clusters parameters saved.
 
   
 //   for(unsigned int i = 0; i < fNPlanes; ++i){
 //     std::sort(hitlist_all[i].begin(), hitlist_all[i].end(),SortByWire());
 //    
 //   }
  
   //find best set:
   unsigned int bp1 = 0,bp2 = 0;
   double minerror1=99999999,minerror2=9999999;
   for(unsigned int ii = 0; ii < fNPlanes; ++ii)
   {
     double locerror=fWire_vertexError[ii]*fWire_vertexError[ii]+ fTime_vertexError[ii]*fTime_vertexError[ii];  // time coordinate of vertex for each plane
 
     if(minerror1 >= locerror )  // the >= sign is to favorize collection
     {
      minerror1=locerror;
      bp1=ii;
     }
   }
   for(unsigned int ij = 0; ij < fNPlanes; ++ij) 
   {
     double locerror=fWire_vertexError[ij]*fWire_vertexError[ij]+ fTime_vertexError[ij]*fTime_vertexError[ij];  // time coordinate of vertex for each plane
     
     if(minerror2 >= locerror && ij!=bp1 )
     {
      minerror2=locerror;
      bp2=ij;
     }
   }
   //bp1=0;
   //bp1=2;
   //std::cout << " best planes " << bp1 << " " << bp2 << std::endl;
 for(unsigned int ij = 0; ij < fNPlanes; ++ij) 
   {
    //std::cout << " wire distances: " << ij << " " << fabs((double)fWire_vertex[ij]-(double)fWire_last[ij]);
     
   }
     
     
  //std::cout << " angles in:  " << bp1 << " " << bp2 << " "   << slope[bp1]*TMath::Pi()/180. << " " << slope[bp2]*TMath::Pi()/180. << " " << slope[bp1] << " " << slope[bp2] << std::endl;
 // gser.Get3DaxisN(bp1,bp2,slope[bp1]*TMath::Pi()/180.,slope[bp2]*TMath::Pi()/180.,xphi,xtheta);
  gser.Get3DaxisN(bp1,bp2,angle[bp1],angle[bp2],xphi,xtheta);
      
  const double origin[3] = {0.};
  std::vector <std::vector <  double > > position;
  // auto const* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();
  double fTimeTick=detprop->SamplingRate()/1000.; 
  double fDriftVelocity=detprop->DriftVelocity(detprop->Efield(),detprop->Temperature());
  // get starting positions for all planes
  for(unsigned int xx=0;xx<fNPlanes;xx++){
       double pos1[3];
       geom->Plane(xx).LocalToWorld(origin, pos1);
       std::vector <double > pos2;
       pos2.push_back(pos1[0]);
       pos2.push_back(pos1[1]);
       pos2.push_back(pos1[2]);
       position.push_back(pos2);
     }
   // Assuming there is no problem ( and we found the best pair that comes close in time )
   // we try to get the Y and Z coordinates for the start of the shower. 
     try{
         int chan1=geom->PlaneWireToChannel(bp1,fWire_vertex[bp1], 0);
         int chan2=geom->PlaneWireToChannel(bp2,fWire_vertex[bp2], 0);
 
         double y,z;
 //      bool wires_cross = geom->ChannelsIntersect(chan1,chan2,y,z);
         geom->ChannelsIntersect(chan1,chan2,y,z);
  //       geom->ChannelsIntersect(chan1,chan2,y,z);
         
         xyz_vertex_fit[1]=y;
         xyz_vertex_fit[2]=z;
         xyz_vertex_fit[0]=(fTime_vertex[bp1]-detprop->TriggerOffset()) *fDriftVelocity*fTimeTick+position[0][0];
 
 
         //std::cout << ":::::: found x,y,z vertex " << wires_cross << " " << xyz_vertex_fit[0] << " " << y << " " << z << " " << wires_cross << std::endl;
     }
     catch(cet::exception e) {
       mf::LogWarning("ShowerReco") << "caught exception \n" << e;
       xyz_vertex_fit[1]=0;
       xyz_vertex_fit[2]=0;
       xyz_vertex_fit[0]=0;
      }
   
         
    // if collection is not best plane, project starting point from that   
       if(bp1!=fNPlanes-1 && bp2!=fNPlanes-1)
       {
         double pos[3];
         unsigned int  wirevertex;
 
         geom->Plane(fNPlanes-1).LocalToWorld(origin, pos);
         //planex[p] = pos[0];
         //std::cout << "plane X positionp " << 2 << " " << pos[0] << std::endl;
 
         pos[1]=xyz_vertex_fit[1];
         pos[2]=xyz_vertex_fit[2];
         wirevertex = geom->NearestWire(pos,fNPlanes-1);
         //geo->ChannelToWire(channel2,cs,t,p,wirevertex);   //\fixme!
         //wirevertex=  (*a)->WireID().Wire;
         
         
         
         
      double drifttick=(xyz_vertex_fit[0]/detprop->DriftVelocity(detprop->Efield(),detprop->Temperature()))*(1./fTimeTick);
      fWire_vertex[fNPlanes-1]= wirevertex;  // wire coordinate of vertex for each plane
      fTime_vertex[fNPlanes-1] = drifttick-(pos[0]/detprop->DriftVelocity(detprop->Efield(),detprop->Temperature()))*(1./fTimeTick)+detprop->TriggerOffset();
         
         
       }
       
       
    //  std::cout << "^^^^^^cross-check xphi and xtheta: " << xphi << " " << xtheta << std::endl;
       
      
       if(fabs(xphi) < 5. )
       { xtheta= gser.Get3DSpecialCaseTheta(bp1,bp2,fWire_last[bp1]-fWire_vertex[bp1], fWire_last[bp2]-fWire_vertex[bp2]);
         
         //std::cout << "xphi, xtheta1:,alt" << xphi << " " << xtheta  <<std::endl;
       }
     //}
     
 // zero the arrays just to make sure    
   for(unsigned int i = 0; i < fNAngles; ++i){
     fTotChargeADC[i]=0;
     fTotChargeMeV[i]=0;
     fTotChargeMeV_MIPs[i]=0;
     fNpoints_corr_ADC_2cm[i]=0;
     fNpoints_corr_MeV_2cm[i]=0;
         
     fRMS_2cm[i]=0; 
     fNpoints_2cm[i]=0; 
 
     fCorr_MeV_2cm[i]=0; 
     fCorr_Charge_2cm[i]=0;
     
     fChargeADC_2cm[i]=0;  //Initial charge in ADC/cm for each plane angle calculation 4cm
     fChargeMeV_2cm[i]=0;  //initial charge in MeV/cm for each angle calculation first 4cm
     
   }
 
 
   // do loop and choose Collection. With ful calorimetry can do all.
   //for(unsigned int set=0;set<fNAngles;set++)
    if(!(fabs(xphi) >89 && fabs(xphi)<91)) // do not calculate pitch for extreme angles
     LongTransEnergy(0,hitlist_all[fNPlanes-1]); //temporary only plane 2. 
   
     
   
    //std::vector< recob::SpacePoint >   spacepoints = std::vector<recob::SpacePoint>()  
   
 
          
   // make an art::PtrVector of the clusters
   art::PtrVector<recob::Cluster> prodvec;
   for(unsigned int i = 0; i < clusterListHandle->size(); ++i){
     art::Ptr<recob::Cluster> prod(clusterListHandle, i);
     prodvec.push_back(prod);
   }
 
   //create a singleSpacePoint at vertex.
   std::vector< recob::SpacePoint > spcpts;
   
   //get direction cosines and set them for the shower 
   // TBD determine which angle to use for the actual shower
   double fPhi=xphi;
   double fTheta=xtheta;
   
   TVector3 dcosVtx(TMath::Cos(fPhi*TMath::Pi()/180)*TMath::Sin(fTheta*TMath::Pi()/180),
                    TMath::Cos(fTheta*TMath::Pi()/180),
                    TMath::Sin(fPhi*TMath::Pi()/180)*TMath::Sin(fTheta*TMath::Pi()/180));
   // fill with bogus values for now
   TVector3 dcosVtxErr(util::kBogusD, util::kBogusD, util::kBogusD);
   //double maxTransWidth[2] = { util::kBogusD };
   //double distMaxWidth = util::kBogusD;
   //recob::Shower  singShower(dcosVtx, dcosVtxErr, maxTransWidth, distMaxWidth, 1);
   recob::Shower singShower;
   singShower.set_direction(dcosVtx);
   singShower.set_direction_err(dcosVtxErr);
   
   Shower3DVector->push_back(singShower);
   // associate the shower with its clusters
   util::CreateAssn(*this, evt, *Shower3DVector, prodvec, *cassn);
   
    // get the hits associated with each cluster and associate those with the shower
    for(size_t p = 0; p < prodvec.size(); ++p){
       std::vector< art::Ptr<recob::Hit> > hits = fmh.at(p);
       util::CreateAssn(*this, evt, *Shower3DVector, hits, *hassn);
    }
   
  
    geo::PlaneID planeID(0,0,fNPlanes-1);
    calorimetrycol->push_back(anab::Calorimetry(Kin_En,
                                                vdEdx,
                                                vdQdx,
                                                vresRange,
                                                deadwire,
                                                Trk_Length,
                                                fTrkPitchC,
                                                planeID));
 
   art::PtrVector < recob::Shower >  ssvec;
         
     //for(unsigned int ip=0;ip<1;ip++)  {
         art::ProductID aid = this->getProductID< std::vector < recob::Shower > >();
         art::Ptr< recob::Shower > aptr(aid, 0, evt.productGetter(aid));
         ssvec.push_back(aptr);
       //}
     
   
   //util::CreateAssn(*this, evt, *Shower3DVector, calorimetrycol, *calassn);
    util::CreateAssn(*this, evt, *calorimetrycol,ssvec,*calassn);
   //std::cout << " filling tree @ run, evt," <<  fRun << " " << fEvent << std::endl;
   ftree_shwf->Fill();
   
   //for(unsigned int iplane = 0; iplane < fNPlanes; ++iplane)
   //fh_theta[iplane]->Write(Form("fh_theta_%d_%d",iplane,evt.id().event()));
   // This needs work, clearly.  
   //for(int p=0;p<2;p++)Shower3DVector->push_back(shower);
  
   } // end loop on Vectors of "Associated clusters"
 
   evt.put(std::move(Shower3DVector));
   evt.put(std::move(cassn));  
   evt.put(std::move(hassn));
   evt.put(std::move(calorimetrycol));
   evt.put(std::move(calassn));
   
   
   
 }

void shwf::ShowerReco::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 223 of file ShowerReco_module.cc.

References fcalodEdxlength, fCaloPSet, fClusterModuleLabel, fdEdxlength, fUseArea, and fhicl::ParameterSet::get().

Referenced by ShowerReco().

 {
   fClusterModuleLabel = pset.get< std::string >("ClusterModuleLabel");
 //  fVertexCLusterModuleLabel=pset.get<std::string > ("VertexClusterModuleLabel");
   fCaloPSet=pset.get< fhicl::ParameterSet >("CaloAlg");
   
   fdEdxlength= pset.get< double >("dEdxlength");   //distance that gets used to determine e/gamma separation
   fcalodEdxlength= pset.get< double >("calodEdxlength");  // cutoff distance for hits saved to the calo object.  
   fUseArea= pset.get< bool >("UseArea"); 
   
   return;
 }

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().

 {
   if (!moduleContext_)
     return;
 
   // If none of the branches have missing consumes statements, exit early.
   if (std::all_of(cbegin(missingConsumes_),
                   cend(missingConsumes_),
                   [](auto const& perBranch) { return perBranch.empty(); }))
     return;
 
   constexpr cet::HorizontalRule rule{60};
   mf::LogPrint log{"MTdiagnostics"};
   log << '\n'
       << rule('=') << '\n'
       << "The following consumes (or mayConsume) statements are missing from\n"
       << module_context(moduleDescription_) << '\n'
       << rule('-') << '\n';
 
   cet::for_all_with_index(
     missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
       for (auto const& pi : perBranch) {
         log << "  "
             << assemble_consumes_statement(static_cast<BranchType>(i), pi)
             << '\n';
       }
     });
   log << rule('=');
 }

void art::Consumer::validateConsumedProduct	(	BranchType const	bt,
		ProductInfo const &	pi
	)

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

 {
   // Early exits if consumes tracking has been disabled or if the
   // consumed product is an allowed consumable.
   if (!moduleContext_)
     return;
 
   if (cet::binary_search_all(consumables_[bt], pi))
     return;
 
   if (requireConsumes_) {
     throw Exception(errors::ProductRegistrationFailure,
                     "Consumer: an error occurred during validation of a "
                     "retrieved product\n\n")
       << "The following consumes (or mayConsume) statement is missing from\n"
       << module_context(moduleDescription_) << ":\n\n"
       << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
   }
 
   missingConsumes_[bt].insert(pi);
 }

Member Data Documentation

float shwf::ShowerReco::angle[3]

private

Definition at line 108 of file ShowerReco_module.cc.

Referenced by GetVertexAndAnglesFromCluster(), and produce().

std::vector<float> shwf::ShowerReco::deadwire

private

Definition at line 170 of file ShowerReco_module.cc.

Referenced by produce().

const detinfo::DetectorProperties* shwf::ShowerReco::detprop = lar::providerFrom<detinfo::DetectorPropertiesService>()

private

Definition at line 194 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and produce().

float shwf::ShowerReco::fcalodEdxlength

private

Definition at line 174 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and reconfigure().

fhicl::ParameterSet shwf::ShowerReco::fCaloPSet

private

Definition at line 119 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and reconfigure().

std::vector<double> shwf::ShowerReco::fChargeADC_2cm

private

Definition at line 133 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and produce().

std::vector<double> shwf::ShowerReco::fChargeMeV_2cm

private

Definition at line 134 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fChargeMeV_2cm_axsum

private

Definition at line 137 of file ShowerReco_module.cc.

Referenced by beginJob(), and beginRun().

std::vector<double> shwf::ShowerReco::fChargeMeV_2cm_refined

private

Definition at line 136 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and LongTransEnergy().

std::string shwf::ShowerReco::fClusterModuleLabel

private

Definition at line 111 of file ShowerReco_module.cc.

Referenced by produce(), and reconfigure().

std::vector<double> shwf::ShowerReco::fCorr_Charge_2cm

private

Definition at line 124 of file ShowerReco_module.cc.

Referenced by beginRun(), and produce().

std::vector<double> shwf::ShowerReco::fCorr_MeV_2cm

private

Definition at line 123 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

float shwf::ShowerReco::fdEdxlength

private

Definition at line 173 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and reconfigure().

std::vector<std::vector<double> > shwf::ShowerReco::fDistribChargeADC

private

Definition at line 139 of file ShowerReco_module.cc.

Referenced by beginJob(), and beginRun().

std::vector<std::vector<double> > shwf::ShowerReco::fDistribChargeMeV

private

Definition at line 140 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and LongTransEnergy().

std::vector<std::vector<double> > shwf::ShowerReco::fDistribChargeposition

private

Definition at line 142 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and LongTransEnergy().

std::vector<std::vector<double> > shwf::ShowerReco::fDistribHalfChargeMeV

private

Definition at line 141 of file ShowerReco_module.cc.

Referenced by beginJob(), and beginRun().

double shwf::ShowerReco::fDriftVelocity

private

Definition at line 186 of file ShowerReco_module.cc.

Referenced by beginRun(), and produce().

int shwf::ShowerReco::fEvent

private

Definition at line 104 of file ShowerReco_module.cc.

Referenced by beginJob(), and produce().

double shwf::ShowerReco::fMean_wire_pitch

private

Definition at line 118 of file ShowerReco_module.cc.

Referenced by beginJob().

unsigned int shwf::ShowerReco::fNAngles

private

Definition at line 180 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and produce().

std::vector< int > shwf::ShowerReco::fNhitsperplane

private

Definition at line 189 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and produce().

std::vector< std::vector<double> > shwf::ShowerReco::fNPitch

private

Definition at line 163 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and LongTransEnergy().

unsigned int shwf::ShowerReco::fNPlanes

private

Definition at line 179 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<int> shwf::ShowerReco::fNpoints_2cm

private

Definition at line 122 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<int> shwf::ShowerReco::fNpoints_corr_ADC_2cm

private

Definition at line 126 of file ShowerReco_module.cc.

Referenced by beginRun(), and produce().

std::vector<int> shwf::ShowerReco::fNpoints_corr_MeV_2cm

private

Definition at line 127 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fRMS_2cm

private

Definition at line 121 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

int shwf::ShowerReco::fRun

private

Definition at line 104 of file ShowerReco_module.cc.

Referenced by beginJob(), and produce().

std::vector<std::vector<double> > shwf::ShowerReco::fSingleEvtAngle

private

Definition at line 144 of file ShowerReco_module.cc.

std::vector<std::vector<double> > shwf::ShowerReco::fSingleEvtAngleVal

private

Definition at line 145 of file ShowerReco_module.cc.

int shwf::ShowerReco::fSubRun

private

Definition at line 104 of file ShowerReco_module.cc.

Referenced by beginJob(), and produce().

std::vector<double> shwf::ShowerReco::fTime_last

private

Definition at line 155 of file ShowerReco_module.cc.

Referenced by beginRun(), and GetVertexAndAnglesFromCluster().

std::vector<double> shwf::ShowerReco::fTime_vertex

private

Definition at line 149 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fTime_vertexError

private

Definition at line 152 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), and produce().

float shwf::ShowerReco::ftimetick

private

Definition at line 113 of file ShowerReco_module.cc.

Referenced by beginJob().

double shwf::ShowerReco::fTimeTick

private

Definition at line 185 of file ShowerReco_module.cc.

Referenced by beginRun(), and produce().

std::vector< double > shwf::ShowerReco::fTotADCperplane

private

Definition at line 190 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and produce().

std::vector<double> shwf::ShowerReco::fTotChargeADC

private

Definition at line 129 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fTotChargeMeV

private

Definition at line 130 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fTotChargeMeV_MIPs

private

Definition at line 131 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and produce().

unsigned int shwf::ShowerReco::fTPC

private

Definition at line 178 of file ShowerReco_module.cc.

TTree* shwf::ShowerReco::ftree_shwf

private

Definition at line 181 of file ShowerReco_module.cc.

Referenced by beginJob(), and produce().

float shwf::ShowerReco::fTrkPitchC

private

Definition at line 172 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and produce().

bool shwf::ShowerReco::fUseArea

private

Definition at line 175 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and reconfigure().

std::vector<unsigned int> shwf::ShowerReco::fWire_last

private

Definition at line 154 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), and produce().

std::vector<unsigned int> shwf::ShowerReco::fWire_vertex

private

Definition at line 148 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), LongTransEnergy(), and produce().

std::vector<double> shwf::ShowerReco::fWire_vertexError

private

Definition at line 151 of file ShowerReco_module.cc.

Referenced by beginRun(), GetVertexAndAnglesFromCluster(), and produce().

double shwf::ShowerReco::fWirePitch

private

Definition at line 184 of file ShowerReco_module.cc.

Referenced by beginRun().

double shwf::ShowerReco::fWireTimetoCmCm

private

Definition at line 187 of file ShowerReco_module.cc.

Referenced by beginRun(), and LongTransEnergy().

art::ServiceHandle<geo::Geometry> shwf::ShowerReco::geom

private

Definition at line 193 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

float shwf::ShowerReco::Kin_En

private

Definition at line 166 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and produce().

float shwf::ShowerReco::slope[3]

private

Definition at line 107 of file ShowerReco_module.cc.

Referenced by GetVertexAndAnglesFromCluster(), and LongTransEnergy().

double shwf::ShowerReco::totCnrg

private

Definition at line 117 of file ShowerReco_module.cc.

Referenced by LongTransEnergy().

double shwf::ShowerReco::totCnrg_corr

private

Definition at line 117 of file ShowerReco_module.cc.

Referenced by LongTransEnergy().

float shwf::ShowerReco::Trk_Length

private

Definition at line 171 of file ShowerReco_module.cc.

Referenced by LongTransEnergy(), and produce().

std::vector<float> shwf::ShowerReco::vdEdx

private

Definition at line 167 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

std::vector<float> shwf::ShowerReco::vdQdx

private

Definition at line 169 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

std::vector<float> shwf::ShowerReco::vresRange

private

Definition at line 168 of file ShowerReco_module.cc.

Referenced by beginRun(), LongTransEnergy(), and produce().

double shwf::ShowerReco::xphi

private

Definition at line 177 of file ShowerReco_module.cc.

Referenced by beginJob(), LongTransEnergy(), and produce().

double shwf::ShowerReco::xtheta

private

Definition at line 177 of file ShowerReco_module.cc.

Referenced by beginJob(), LongTransEnergy(), and produce().

double shwf::ShowerReco::xyz_vertex[3]

private

Definition at line 114 of file ShowerReco_module.cc.

std::vector<double> shwf::ShowerReco::xyz_vertex_fit

private

Definition at line 161 of file ShowerReco_module.cc.

Referenced by beginJob(), beginRun(), and produce().

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

/cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/ShowerFinder/ShowerReco_module.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation