Inheritance diagram for gar::anatest:

Public Member Functions
	anatest (fhicl::ParameterSet const &p)

	anatest (anatest const &)=delete

	anatest (anatest &&)=delete

anatest &	operator= (anatest const &)=delete

anatest &	operator= (anatest &&)=delete

virtual void	beginJob () override

void	analyze (art::Event const &e) override

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 Member Functions
void	ClearVectors ()

void	FillVectors (art::Event const &e)

Private Attributes
double	ItsInTulsa [3]

std::vector< std::string >	fGeneratorLabels

std::string	fGeantLabel

std::string	fHitLabel

std::string	fTPCClusterLabel

std::string	fTrackLabel

std::string	fVertexLabel

std::string	fRawCaloHitLabel

std::string	fCaloHitLabel

std::string	fClusterLabel

std::string	fECALAssnLabel

TTree *	fTree

Int_t	fEvent
	number of the event being processed More...

Int_t	fRun
	number of the run being processed More...

Int_t	fSubRun
	number of the sub-run being processed More...

std::vector< Int_t >	fNeutrinoType

std::vector< Int_t >	fInteractionType

std::vector< Float_t >	fQ2

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 59 of file anatest_module.cc.

Constructor & Destructor Documentation

gar::anatest::anatest ( fhicl::ParameterSet const & p )

explicit

Definition at line 119 of file anatest_module.cc.

                                                : EDAnalyzer(p) {
 
     bool usegenlabels =
         p.get_if_present<std::vector<std::string> >("GeneratorLabels",fGeneratorLabels);
     if (!usegenlabels) fGeneratorLabels.clear();
 
     fGeantLabel       = p.get<std::string>("GEANTLabel","geant");
     fHitLabel         = p.get<std::string>("HitLabel","hit");
     fTPCClusterLabel  = p.get<std::string>("TPCClusterLabel","tpccluster");
     fTrackLabel       = p.get<std::string>("TrackLabel","track");
     fVertexLabel      = p.get<std::string>("VertexLabel","vertex");
     fRawCaloHitLabel  = p.get<std::string>("RawCaloHitLabel","daqecal");
     fCaloHitLabel     = p.get<std::string>("CaloHitLabel","calohit");
     fClusterLabel     = p.get<std::string>("ClusterLabel","calocluster");
     fECALAssnLabel    = p.get<std::string>("ECALAssnLabel","trkecalassn");
 
     fTree = nullptr;
 
     if (usegenlabels) {
         for (size_t i=0; i<fGeneratorLabels.size(); ++i) {
             consumes<std::vector<simb::MCTruth> >(fGeneratorLabels.at(i));
         }
     } else {
         consumesMany<std::vector<simb::MCTruth> >();
     }
 
 
 
     consumes<art::Assns<simb::MCTruth, simb::MCParticle> >(fGeantLabel);
 
     consumes<std::vector<rec::Hit> >(fHitLabel);
     consumes<std::vector<rec::TPCCluster> >(fTPCClusterLabel);
     consumes<std::vector<rec::Track> >(fTrackLabel);
     consumes<art::Assns<rec::Track, rec::TPCCluster> >(fTPCClusterLabel);
     consumes<std::vector<rec::Vertex> >(fVertexLabel);
     consumes<art::Assns<rec::Track, rec::Vertex> >(fVertexLabel);
 
     consumes<std::vector<gar::sdp::CaloDeposit> >(fGeantLabel);
     consumes<std::vector<raw::CaloRawDigit> >(fRawCaloHitLabel);
     consumes<std::vector<rec::CaloHit> >(fCaloHitLabel);
     consumes<std::vector<rec::Cluster> >(fClusterLabel);
     consumes<art::Assns<rec::Cluster, rec::Track>>(fECALAssnLabel);
 
     return;
 } // end constructor

gar::anatest::anatest ( anatest const & )

delete

gar::anatest::anatest ( anatest && )

delete

Member Function Documentation

void gar::anatest::analyze ( art::Event const & e )

overridevirtual

Implements art::EDAnalyzer.

Definition at line 200 of file anatest_module.cc.

                                            {
 
     ClearVectors();
     FillVectors(e);
     fTree->Fill();
     return;
 }

void gar::anatest::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 170 of file anatest_module.cc.

                           {
 
     art::ServiceHandle<geo::GeometryGAr> euclid;
     ItsInTulsa[0] = euclid->TPCXCent();
     ItsInTulsa[1] = euclid->TPCYCent();
     ItsInTulsa[2] = euclid->TPCZCent();
 
 
 
     art::ServiceHandle<art::TFileService> tfs;
     fTree = tfs->make<TTree>("MPDtestTree","MPDtestTree");
 
 
 
     fTree->Branch("Run",           &fRun,         "Run/I");
     fTree->Branch("SubRun",        &fSubRun,      "SubRun/I");
     fTree->Branch("Event",         &fEvent,       "Event/I");
 
     fTree->Branch("NType",         &fNeutrinoType);
     fTree->Branch("InterT",        &fInteractionType);
     fTree->Branch("MC_Q2",         &fQ2);
 
     return;
 }  // End of :anatest::beginJob

void gar::anatest::ClearVectors ( )

private

Definition at line 213 of file anatest_module.cc.

                               {
 
     // clear out all our vectors
     fNeutrinoType.clear();
     fInteractionType.clear();
     fQ2.clear();
 
     return;
 } // end :anatest::ClearVectors

void gar::anatest::FillVectors ( art::Event const & e )

private

Definition at line 228 of file anatest_module.cc.

                                                {
 
 
 
     // =============  Get art handles ==========================================
     // Get handles for MCinfo, also good for MCPTrajectory
     std::vector< art::Handle<std::vector<simb::MCTruth>> > mcthandlelist;
     if (fGeneratorLabels.size()<1) {
         mcthandlelist = e.getMany<std::vector<simb::MCTruth> >();    // get them all (even if there are none)
     } else {
         mcthandlelist.resize(fGeneratorLabels.size());
         for (size_t i=0; i< fGeneratorLabels.size(); ++i) {
           mcthandlelist.at(i) = e.getHandle<std::vector<simb::MCTruth> >(fGeneratorLabels.at(i));
             if (!mcthandlelist.at(i))
                 throw cet::exception("anatest") << " No simb::MCTruth branch."
                      << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
         }
     }
 
 
 
     auto MCPHandle = e.getHandle<std::vector<simb::MCParticle> >(fGeantLabel);
     if (!MCPHandle) {
         throw cet::exception("anatest") << " No simb::MCParticle branch."
             << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
 
 
 
     // Get handles for MCCaloInfo
     auto SimHitHandle = e.getHandle< std::vector<gar::sdp::CaloDeposit> >(fGeantLabel);
     if (!SimHitHandle) {
         throw cet::exception("anatest") << " No gar::sdp::CaloDeposit branch."
             << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
 
 
 
     // Get handle for TPC hit data
     auto HitHandle = e.getHandle< std::vector<rec::Hit> >(fHitLabel);
     if (!HitHandle) {
         throw cet::exception("anatest") << " No rec::Hit branch."
         << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
 
 
 
     // Get handle for TPCClusters
     auto TPCClusterHandle = e.getHandle< std::vector<rec::TPCCluster> >(fTPCClusterLabel);
     if (!TPCClusterHandle) {
         throw cet::exception("anatest") << " No rec::TPCCluster branch."
         << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
 
 
 
     // Get handles for Tracks and also Assn's to TPCClusters, TrackIoniz
     auto TrackHandle = e.getHandle< std::vector<rec::Track> >(fTrackLabel);
     if (!TrackHandle) {
         throw cet::exception("anatest") << " No rec::Track branch."
         << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
     //art::FindManyP<rec::TPCCluster>* findManyTPCClusters = NULL;
     //findManyTPCClusters = new art::FindManyP<rec::TPCCluster>(TrackHandle,e,fTrackLabel);
 
 
 
     // Get handle for Vertices; also Assn's to Tracks
     auto VertexHandle = e.getHandle< std::vector<rec::Vertex> >(fVertexLabel);
     if (!VertexHandle) {
         throw cet::exception("anatest") << " No rec::Vertex branch."
         << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
     //art::FindManyP<rec::Track, rec::TrackEnd>* findManyTrackEnd = NULL;
     //findManyTrackEnd = new art::FindManyP<rec::Track, rec::TrackEnd>(VertexHandle,e,fVertexLabel);
 
 
     // Get handle for CaloDigits
     auto RawHitHandle = e.getHandle< std::vector<gar::raw::CaloRawDigit> >(fRawCaloHitLabel);
     if (!RawHitHandle) {
         throw cet::exception("anatest") << " No :raw::CaloRawDigit branch."
         << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
 
 
     // Get handle for CaloHits
     auto RecoHitHandle = e.getHandle< std::vector<rec::CaloHit> >(fCaloHitLabel);
     if (!RecoHitHandle) {
         throw cet::exception("anatest") << " No rec::CaloHit branch."
         << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
 
 
     // Get handle for CaloClusters; also Assn for matching tracks
     auto RecoClusterHandle = e.getHandle< std::vector<rec::Cluster> >(fClusterLabel);
     if (!RecoClusterHandle) {
         throw cet::exception("anatest") << " No rec::Cluster branch."
         << " Line " << __LINE__ << " in file " << __FILE__ << std::endl;
     }
     //art::FindManyP<rec::Track, rec::TrackEnd>* findManyCALTrackEnd = NULL;
     //findManyCALTrackEnd = new art::FindManyP<rec::Track, rec::TrackEnd>
     //                                   (RecoClusterHandle,e,fECALAssnLabel);
 
 
 
     // =============  Pull art handles =========================================
     fRun    = e.run();
     fSubRun = e.subRun();
     fEvent  = e.id().event();
 
     for (size_t imchl = 0; imchl < mcthandlelist.size(); ++imchl) {
         for ( auto const& mct : (*mcthandlelist.at(imchl)) ) {
             if (mct.NeutrinoSet()) {
                 simb::MCNeutrino nuw = mct.GetNeutrino();
                 fNeutrinoType.push_back(nuw.Nu().PdgCode());
                 fInteractionType.push_back(nuw.InteractionType());
             }
         }
     }
 
     // Need a non-constant backtracker instance, for now.
     cheat::BackTrackerCore const* const_bt = gar::providerFrom<cheat::BackTracker>();
     cheat::BackTrackerCore*             bt = const_cast<cheat::BackTrackerCore*>(const_bt);
 
     bool const dumpMCP = false;
     sim::ParticleList* partList = bt->GetParticleList();
     int nMotherless = 0;
     for ( auto const& mcp : *MCPHandle ) {
         if (mcp.Mother() > 0) {
             simb::MCParticle* TPCeve = bt->FindTPCEve(mcp.TrackId());
             if ( TPCeve->TrackId() == mcp.TrackId() ) ++nMotherless;
         }
         if (dumpMCP) {
             std::cout << "TrackID: " << mcp.TrackId() << " is PDG: " << mcp.PdgCode() <<
                 " with mother " << mcp.Mother() << " produced by process " << mcp.Process()
                 << std::endl;
         }
     }
     if (dumpMCP) {
         std::cout << "Number of children: " << partList->size() << std::endl;
         std::cout << "Number of motherless children: " << nMotherless << std::endl;
     }
 
 
 
 
     std::vector<art::Ptr<rec::Hit>> allhits;
     art::PtrMaker<rec::Hit> makeHitPtr(e,HitHandle.id());
     for (size_t iHit=0; iHit<HitHandle->size(); ++iHit) {
         art::Ptr<rec::Hit> aPtr = makeHitPtr(iHit);
         allhits.push_back(aPtr);
     }
     for (simb::MCParticle mcp : *MCPHandle) {
                 if ( mcp.Mother() > 0 ) continue;
             if ( abs(mcp.PdgCode()) == 13 ) bt->ParticleToHits(&mcp,allhits);
         }
 
 
 
 
 
 /*  for ( rec::Cluster cluster : *RecoClusterHandle ) {
         std::vector<std::pair<simb::MCParticle*,float>> whatMatches;
         whatMatches = bt->ClusterToMCParticles(&cluster);
         std::cout << "\nCluster No. " << cluster.getIDNumber() << " is made of MCParticles: " << std::endl;
         for (auto itr = whatMatches.begin(); itr!=whatMatches.end(); ++itr) {
             std::cout << "G4 track number " << itr->first->TrackId() << "\thas PDG code " <<
                 itr->first->PdgCode() << "\tand its mother is G4 track " << itr->first->Mother()
                 << "\tand energy fraction " << 100*(itr->second) << "%\n";
         }
     }
 
     std::vector<art::Ptr<rec::Cluster>> clusterCol;
     art::PtrMaker<rec::Cluster> makeClusterPtr(e,RecoClusterHandle.id());
     for (size_t iCluster=0; iCluster<RecoClusterHandle->size(); ++iCluster ) {
         art::Ptr<rec::Cluster> aPtr = makeClusterPtr(iCluster);
         clusterCol.push_back(aPtr);
     }
     for ( simb::MCParticle mcp : *MCPHandle ) {
         if (mcp.Mother() == 0) {
             std::cout << "Particle PDG: " << mcp.PdgCode() << " matches to\n";
             std::vector<art::Ptr<rec::Cluster>> clusterList =
                 bt->MCParticleToClusters(&mcp, clusterCol);
             for (art::Ptr<rec::Cluster> iCluster : clusterList)
                 std::cout << iCluster->getIDNumber() << std::endl;
         }
     } */
 
 
 
 
 
     return;
 }

anatest& gar::anatest::operator= ( anatest const & )

delete

anatest& gar::anatest::operator= ( anatest && )

delete

Member Data Documentation

std::string gar::anatest::fCaloHitLabel

private

Definition at line 94 of file anatest_module.cc.

std::string gar::anatest::fClusterLabel

private

Definition at line 95 of file anatest_module.cc.

std::string gar::anatest::fECALAssnLabel

private

Definition at line 96 of file anatest_module.cc.

Int_t gar::anatest::fEvent

private

number of the event being processed

Definition at line 102 of file anatest_module.cc.

std::string gar::anatest::fGeantLabel

private

Definition at line 88 of file anatest_module.cc.

std::vector<std::string> gar::anatest::fGeneratorLabels

private

Definition at line 87 of file anatest_module.cc.

std::string gar::anatest::fHitLabel

private

Definition at line 89 of file anatest_module.cc.

std::vector<Int_t> gar::anatest::fInteractionType

private

Definition at line 108 of file anatest_module.cc.

std::vector<Int_t> gar::anatest::fNeutrinoType

private

Definition at line 107 of file anatest_module.cc.

std::vector<Float_t> gar::anatest::fQ2

private

Definition at line 109 of file anatest_module.cc.

std::string gar::anatest::fRawCaloHitLabel

private

Definition at line 93 of file anatest_module.cc.

Int_t gar::anatest::fRun

private

number of the run being processed

Definition at line 103 of file anatest_module.cc.

Int_t gar::anatest::fSubRun

private

number of the sub-run being processed

Definition at line 104 of file anatest_module.cc.

std::string gar::anatest::fTPCClusterLabel

private

Definition at line 90 of file anatest_module.cc.

std::string gar::anatest::fTrackLabel

private

Definition at line 91 of file anatest_module.cc.

TTree* gar::anatest::fTree

private

Definition at line 99 of file anatest_module.cc.

std::string gar::anatest::fVertexLabel

private

Definition at line 92 of file anatest_module.cc.

double gar::anatest::ItsInTulsa[3]

private

Definition at line 84 of file anatest_module.cc.

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

garsoft/Ana/anatest_module.cc

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation