Inheritance diagram for AnalysisExample::dEdx:

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

virtual	~dEdx ()

void	beginJob ()

void	beginRun (const art::Run &run)

void	reconfigure (fhicl::ParameterSet const &pset)

void	analyze (const art::Event &evt)

void	endJob ()

void	findDepositionVertex (std::vector< sim::SimChannel > SCHandle)

void	checkDepositionVertex (std::vector< sim::SimChannel > SCHandle, double x, double y, double z)

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

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

std::string	workerType () const

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

	Analyzer (fhicl::ParameterSet const &pset)

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

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Observer
	~Observer () noexcept

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

fhicl::ParameterSetID	selectorConfig () const

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

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Attributes
art::ServiceHandle< art::TFileService >	tfs

art::ServiceHandle< geo::Geometry >	fGeom

std::string	fSimulationProducerLabel

TH1D *	fMCdEdxHist

std::vector< double >	fMCdEdxVec

double	fMaxMCdEdx = 0.

double	dE

double	trackLength

TH3D *	fStartEdepHist

TH3D *	fTrackEdepHist

TH3D *	fShowerEdepHist

std::vector< unsigned int >	fNullList

double	xi = 0.

double	yi = 0.

double	zi = 1000.

double	ziFalse = -500.

unsigned int	nearDeps = 0

bool	goodDepVertex = false

double	fConvDist

TH1D *	fConvDistHist

double	Vxi

double	Vyi

double	Vzi

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 53 of file dEdx_module.cc.

Constructor & Destructor Documentation

AnalysisExample::dEdx::dEdx ( fhicl::ParameterSet const & pset )

explicit

Definition at line 112 of file dEdx_module.cc.

     : EDAnalyzer(parameterSet)
   {
     this->reconfigure(parameterSet);
   }

AnalysisExample::dEdx::~dEdx ( )

virtual

Definition at line 120 of file dEdx_module.cc.

121 {

122 }

Member Function Documentation

void AnalysisExample::dEdx::analyze ( const art::Event & evt )

Definition at line 165 of file dEdx_module.cc.

   {
 
     auto simChanHandle = event.getHandle<std::vector<sim::SimChannel>>(fSimulationProducerLabel);
 
     // initialize at the start of every event
     dE = 0.;
     trackLength = 0.;
     goodDepVertex = false;
     unsigned int nTries = 0;
     ziFalse = -500.;
 
     // find conversion point
     while(!goodDepVertex && nTries < 300)
     {
       this->findDepositionVertex( (*simChanHandle) );
       this->checkDepositionVertex( (*simChanHandle), xi, yi, zi );
       nTries++;
     }
 
     std::cout << "Number of attemps at finding deposition vertex = " << nTries << "." << std::endl;
 
     // find energy deposited into first 2.5 cm of "track"
     if(nTries < 300 && goodDepVertex)
     {
       fStartEdepHist->Fill(xi,yi,zi);
       for(auto const& channel : (*simChanHandle))
       {
         if(fGeom->SignalType(channel.Channel()) == geo::kCollection)
         {
           auto const& timeSlices = channel.TDCIDEMap();
           for(auto const& t : timeSlices)
           {
             auto const& eDeps = t.second;
             for(auto const& eDep : eDeps)
             {
               fShowerEdepHist->Fill(eDep.x,eDep.y,eDep.z);
               trackLength = std::sqrt( (eDep.x-xi)*(eDep.x-xi) + (eDep.y-yi)*(eDep.y-yi) + (eDep.z-zi)*(eDep.z-zi) );
               if(trackLength < 2.5 /* && eDep.z >= zi*/)
               {
                 dE += eDep.energy;
                 fTrackEdepHist->Fill(eDep.x,eDep.y,eDep.z);
               }
             } // every energy deposition
           } // in every time slice
         } // for the collection plane
       } // and for every wire
       dE /= 2.5;
       fMCdEdxVec.push_back(dE);
       if(dE > fMaxMCdEdx) fMaxMCdEdx = dE;
       if(dE < 0.5) fNullList.push_back(event.id().event());
     }
 
     auto particleHandle = event.getHandle< std::vector<simb::MCParticle> >(fSimulationProducerLabel);
 
     for(auto const& particle : (*particleHandle))
     {
       if(particle.Process() == "primary")
       {
         Vxi = particle.Vx();
         Vyi = particle.Vy();
         Vzi = particle.Vz();
       }
     }
       fConvDist =  std::sqrt( (Vxi-xi)*(Vxi-xi) + (Vyi-yi)*(Vyi-yi) + (Vzi-zi)*(Vzi-zi) );
       fConvDistHist->Fill(fConvDist);
   }

void AnalysisExample::dEdx::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 126 of file dEdx_module.cc.

   {
 
     // hard-coded fd4apa cryostat dimensions
     fStartEdepHist = tfs->make<TH3D>("fStartEdepHist","fStartEdepHist",50,-240.,240.,50,-750.,800.,50,-50.,560.);
     fTrackEdepHist = tfs->make<TH3D>("fTrackEdepHist","fTrackEdepHist",50,-240.,240.,50,-750.,800.,50,-50.,560.);
     fShowerEdepHist = tfs->make<TH3D>("fShowerEdepHist","fShowerEdepHist",50,-240.,240.,50,-750.,800.,50,-50.,560.);
 
     fStartEdepHist->GetXaxis()->SetTitle("x");
     fStartEdepHist->GetYaxis()->SetTitle("y");
     fStartEdepHist->GetZaxis()->SetTitle("z");
 
     fShowerEdepHist->GetXaxis()->SetTitle("x");
     fShowerEdepHist->GetYaxis()->SetTitle("y");
     fShowerEdepHist->GetZaxis()->SetTitle("z");
 
     fTrackEdepHist->GetXaxis()->SetTitle("x");
     fTrackEdepHist->GetYaxis()->SetTitle("y");
     fTrackEdepHist->GetZaxis()->SetTitle("z");
 
     fConvDistHist = tfs->make<TH1D>("fConvDistHist","fConvDistHist",100,0.,50.);
 
   }

void AnalysisExample::dEdx::beginRun ( const art::Run & run )

Definition at line 152 of file dEdx_module.cc.

153 {

154 }

void AnalysisExample::dEdx::checkDepositionVertex	(	std::vector< sim::SimChannel >	SCHandle,
		double	x,
		double	y,
		double	z
	)

Definition at line 275 of file dEdx_module.cc.

   {
     std::cout << "Checking for Deposition Vertex" << std::endl;
     nearDeps = 0;
     double track = 0;
     for(auto const& channel : SCHandle)
     {
       if(fGeom->SignalType(channel.Channel()) == geo::kCollection)
       {
         auto const& timeSlices = channel.TDCIDEMap();
         for(auto const& t : timeSlices)
         {
           auto const& eDeps = t.second;
           for(auto const& eDep : eDeps)
           {
             track = std::sqrt( (eDep.x-x)*(eDep.x-x) + (eDep.y-y)*(eDep.y-y) + (eDep.z-z)*(eDep.z-z) );
             if(track < 2.5 && eDep.energy > 0.)
             {
               nearDeps++;
             } // if z position is minimum
           } // energy deposition loop
         } // time slice loop 
       } // if collection wire
     } // sim channel loop
     if( nearDeps > 30) goodDepVertex = true;
     else ziFalse = z + 0.05;
   } // checkDepositionVertex

void AnalysisExample::dEdx::endJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 235 of file dEdx_module.cc.

   {
     fMCdEdxHist = tfs->make<TH1D>("fMCdEdxHist",";dE/dx (MeV/cm)",200,0.0,10.);
 
     for(unsigned int e = 0; e < fMCdEdxVec.size(); e++)
     {
       fMCdEdxHist->Fill(fMCdEdxVec[e]);
     }
 
     std::cout << "Events with < 0.5 dE/dx entry: ";
     for(unsigned int a = 0; a < fNullList.size(); a++) std::cout << fNullList[a] << ", ";
 
   }

void AnalysisExample::dEdx::findDepositionVertex ( std::vector< sim::SimChannel > SCHandle )

Definition at line 249 of file dEdx_module.cc.

   {
     zi = 1000.;
     std::cout << "Finding Deposition Vertex... " << std::endl;
     for(auto const& channel : SCHandle)
     {
       if(fGeom->SignalType(channel.Channel()) == geo::kCollection)
       {
         auto const& timeSlices = channel.TDCIDEMap();
         for(auto const& t : timeSlices)
         {
           auto const& eDeps = t.second;
           for(auto const& eDep : eDeps)
           {
             if( (eDep.z < zi) && (eDep.energy > 0.) && (eDep.z > ziFalse) )
             {
               xi = eDep.x;
               yi = eDep.y;
               zi = eDep.z;
             } // if z position is a "valid" minimum
           } // energy deposition loop
         } // time slice loop 
       } // if collection wire
     } // sim channel loop
   } // findDepositionVertex

void AnalysisExample::dEdx::reconfigure ( fhicl::ParameterSet const & pset )

Definition at line 158 of file dEdx_module.cc.

   {
     fSimulationProducerLabel = p.get<std::string>("SimulationLabel");
   }

Member Data Documentation

double AnalysisExample::dEdx::dE

private

Definition at line 77 of file dEdx_module.cc.

double AnalysisExample::dEdx::fConvDist

private

Definition at line 102 of file dEdx_module.cc.

TH1D* AnalysisExample::dEdx::fConvDistHist

private

Definition at line 103 of file dEdx_module.cc.

art::ServiceHandle<geo::Geometry> AnalysisExample::dEdx::fGeom

private

Definition at line 70 of file dEdx_module.cc.

double AnalysisExample::dEdx::fMaxMCdEdx = 0.

private

Definition at line 76 of file dEdx_module.cc.

TH1D* AnalysisExample::dEdx::fMCdEdxHist

private

Definition at line 74 of file dEdx_module.cc.

std::vector<double> AnalysisExample::dEdx::fMCdEdxVec

private

Definition at line 75 of file dEdx_module.cc.

std::vector<unsigned int> AnalysisExample::dEdx::fNullList

private

Definition at line 89 of file dEdx_module.cc.

TH3D* AnalysisExample::dEdx::fShowerEdepHist

private

Definition at line 86 of file dEdx_module.cc.

std::string AnalysisExample::dEdx::fSimulationProducerLabel

private

Definition at line 71 of file dEdx_module.cc.

TH3D* AnalysisExample::dEdx::fStartEdepHist

private

Definition at line 84 of file dEdx_module.cc.

TH3D* AnalysisExample::dEdx::fTrackEdepHist

private

Definition at line 85 of file dEdx_module.cc.

bool AnalysisExample::dEdx::goodDepVertex = false

private

Definition at line 99 of file dEdx_module.cc.

unsigned int AnalysisExample::dEdx::nearDeps = 0

private

Definition at line 98 of file dEdx_module.cc.

art::ServiceHandle<art::TFileService> AnalysisExample::dEdx::tfs

private

Definition at line 69 of file dEdx_module.cc.

double AnalysisExample::dEdx::trackLength

private

Definition at line 80 of file dEdx_module.cc.

double AnalysisExample::dEdx::Vxi

private

Definition at line 104 of file dEdx_module.cc.

double AnalysisExample::dEdx::Vyi

private

Definition at line 105 of file dEdx_module.cc.

double AnalysisExample::dEdx::Vzi

private

Definition at line 106 of file dEdx_module.cc.

double AnalysisExample::dEdx::xi = 0.

private

Definition at line 92 of file dEdx_module.cc.

double AnalysisExample::dEdx::yi = 0.

private

Definition at line 93 of file dEdx_module.cc.

double AnalysisExample::dEdx::zi = 1000.

private

Definition at line 94 of file dEdx_module.cc.

double AnalysisExample::dEdx::ziFalse = -500.

private

Definition at line 97 of file dEdx_module.cc.

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

duneana/duneana/EnergyStudies/dEdx_module.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation