larg4::CheckMCParticle Class Reference

Inheritance diagram for larg4::CheckMCParticle:

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

Private Member Functions
void	beginJob () override
void	analyze (const art::Event &event) override

Private Attributes
std::string const	_myName
TH1F *	_hnParts {nullptr}

Detailed Description

Definition at line 34 of file CheckMCParticle_module.cc.

Constructor & Destructor Documentation

larg4::CheckMCParticle::CheckMCParticle ( fhicl::ParameterSet const &  p )

explicit

Definition at line 46 of file CheckMCParticle_module.cc.

                                                                :
  art::EDAnalyzer(p),
  _myName(p.get<std::string>("name", "CheckMCParticle"))
{}

Member Function Documentation

void larg4::CheckMCParticle::analyze ( const art::Event &  event )

overrideprivate

Definition at line 57 of file CheckMCParticle_module.cc.

{
#if defined _verbose_
  
  auto allDropped =  event.getMany<std::map<int,std::set<int>>>();

  for (auto const& maps : allDropped ) {
    for  (auto const& element : *maps ) {
    std::cout << "Parent of dropped Tracks: "<<element.first<< std::endl;
    std::set<int> droppedset= element.second; 
    std::cout<< " droppedid size:  "<< droppedset.size()<<std::endl;
    for(auto const& droppedid :  droppedset  )
      {
        std::cout << droppedid << " ";
      }
    std::cout <<std::endl;   
    }
  }

#endif
  auto allGens = event.getMany<std::vector<simb::MCParticle>>();
  for (auto const& gens : allGens) {
    _hnParts->Fill(gens->size());
#if defined _verbose_
    for (auto const& genpart : *gens) {
      if (genpart.Mother()==0) {
        cout << "Primary momentum:  " <<   genpart.P();
        cout << "  position:  " << genpart.Vx()<< "  "<< genpart.Vy()<<"  "<< genpart.Vz()  << endl;
      }
      cout << "Part id:  " << genpart.TrackId();
      cout << " PDG id:  " << genpart.PdgCode();
      cout << " Status Code:  " << genpart.StatusCode();
      cout << " Mother:  " << genpart.Mother();
      cout << " Creation Process: "<<genpart.Process();
      cout << " End Process: "<<genpart.EndProcess();
      /*
      auto trajectory = genpart.Trajectory();
      cout <<" trajectory size:   " << trajectory.size();
      */
      cout <<" Nr. of Daughters: " << genpart.NumberDaughters();
      cout <<" FirstDaughter:" << genpart.FirstDaughter()<<endl;
      //      cout <<" LastDaughter: " << genpart.LastDaughter() <<endl;
      for (int i=0;i<genpart.NumberDaughters();i++)
        {
          cout << genpart.Daughter(i)<<",";
        }
      cout<<endl;
    }
#endif
  }
} // end analyze

void larg4::CheckMCParticle::beginJob ( )

overrideprivate

Definition at line 51 of file CheckMCParticle_module.cc.

{
  art::ServiceHandle<art::TFileService const> tfs;
  _hnParts = tfs->make<TH1F>("hnParts", "Number of generated Particles", 100, 0., 2000.);
} // end beginJob

Member Data Documentation

TH1F* larg4::CheckMCParticle::_hnParts {nullptr}

private

Definition at line 43 of file CheckMCParticle_module.cc.

std::string const larg4::CheckMCParticle::_myName

private

Definition at line 42 of file CheckMCParticle_module.cc.

---

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

• CheckMCParticle_module.cc