NCDeltaRadiative Class Reference

Inheritance diagram for NCDeltaRadiative:

Public Member Functions
	NCDeltaRadiative (fhicl::ParameterSet const &p)
	NCDeltaRadiative (NCDeltaRadiative const &)=delete
	NCDeltaRadiative (NCDeltaRadiative &&)=delete
NCDeltaRadiative &	operator= (NCDeltaRadiative const &)=delete
NCDeltaRadiative &	operator= (NCDeltaRadiative &&)=delete
void	cout_stuff (art::Event &e, bool passed)
void	FillTree (art::Event &e, size_t const mct_index, size_t const parent_index, bool const is_nc_delta_radiative)
void	Reset ()
bool	filter (art::Event &e) override

Private Attributes
TTree *	ftree
int	frun
int	fsubrun
int	fevent
int	fnu_pdg
int	fccnc
int	fmode
int	finteraction_type
int	fis_nc_delta_radiative
int	fparent_status_code
int	fparent_pdg

Detailed Description

Definition at line 30 of file NCDeltaRadiative_module.cc.

Constructor & Destructor Documentation

NCDeltaRadiative::NCDeltaRadiative ( fhicl::ParameterSet const &  p )

explicit

Definition at line 64 of file NCDeltaRadiative_module.cc.

                                                              :
  art::EDFilter(p),
  ftree(nullptr) {

    art::ServiceHandle<art::TFileService> tfs;
    ftree = tfs->make<TTree>("NCDeltaRadiativeFilter", "");

    ftree->Branch("run", &frun, "run/I");
    ftree->Branch("subrun", &fsubrun, "subrun/I");
    ftree->Branch("event", &fevent, "event/I");
    ftree->Branch("nu_pdg", &fnu_pdg, "nu_pdg/I");
    ftree->Branch("ccnc", &fccnc, "ccnc/I");
    ftree->Branch("mode", &fmode, "mode/I");
    ftree->Branch("is_nc_delta_radiative", &fis_nc_delta_radiative, "is_nc_delta_radiative/I");
    ftree->Branch("parent_status_code", &fparent_status_code, "parent_status_code/I");
    ftree->Branch("parent_pdg", &fparent_pdg, "parent_pdg/I");

}

NCDeltaRadiative::NCDeltaRadiative ( NCDeltaRadiative const &  )

delete

NCDeltaRadiative::NCDeltaRadiative ( NCDeltaRadiative &&  )

delete

Member Function Documentation

void NCDeltaRadiative::cout_stuff	(	art::Event &	e,
		bool	passed = false
	)

Definition at line 84 of file NCDeltaRadiative_module.cc.

                                                                   {

  art::ValidHandle<std::vector<simb::MCTruth>> const & ev_mct =
    e.getValidHandle<std::vector<simb::MCTruth>>("generator");

  std::cout << passed << "\n"
      << "==========================\n";
  for(simb::MCTruth const & mct : *ev_mct) {
    std::cout << "----------------------------\n";
    for(int i = 0; i < mct.NParticles(); ++i) {
      simb::MCParticle const & mcp = mct.GetParticle(i);
      std::cout << mcp.TrackId() << " " << mcp.PdgCode() << " " << mcp.Mother() << " " << mcp.StatusCode() << "\n";
    }
  }

}

void NCDeltaRadiative::FillTree	(	art::Event &	e,
		size_t const	mct_index,
		size_t const	parent_index,
		bool const	is_nc_delta_radiative
	)

Definition at line 118 of file NCDeltaRadiative_module.cc.

                                          {

  Reset();

  frun = e.id().run();
  fsubrun = e.id().subRun();
  fevent = e.id().event();

  art::ValidHandle<std::vector<simb::MCTruth>> const & ev_mct =
    e.getValidHandle<std::vector<simb::MCTruth>>("generator");
  simb::MCNeutrino const & mcn = ev_mct->at(mct_index).GetNeutrino(); 

  fnu_pdg = mcn.Nu().PdgCode();
  fccnc = mcn.CCNC();
  fmode = mcn.Mode();
  finteraction_type = mcn.InteractionType();

  fis_nc_delta_radiative = is_nc_delta_radiative;
  if(parent_index != SIZE_MAX) {
    fparent_status_code = ev_mct->at(mct_index).GetParticle(parent_index).StatusCode();
    fparent_pdg = ev_mct->at(mct_index).GetParticle(parent_index).PdgCode();
  }    

  ftree->Fill();

}

bool NCDeltaRadiative::filter ( art::Event &  e )

override

Definition at line 149 of file NCDeltaRadiative_module.cc.

                                          {

  art::ValidHandle<std::vector<simb::MCTruth>> const & ev_mct =
    e.getValidHandle<std::vector<simb::MCTruth>>("generator");

  for(size_t i = 0; i < ev_mct->size(); ++i) {

    simb::MCTruth const & mct = ev_mct->at(i);
    if(mct.GetNeutrino().CCNC() != 1) continue;

    std::vector<size_t> exiting_photon_parents;
    for(int i = 0; i < mct.NParticles(); ++i) {
      simb::MCParticle const & mcp = mct.GetParticle(i);
      if(mcp.TrackId() != i) {
  std::cout << "ERROR: " << __LINE__ << " " << __PRETTY_FUNCTION__ << "\nTrackId does not match index\n";
  exit(1);
      }
      if(!(mcp.StatusCode() == 1 && mcp.PdgCode() == 22)) continue;
      exiting_photon_parents.push_back(mcp.Mother());
    }

    std::vector<size_t> in_nucleus_photons;
    for(size_t const s : exiting_photon_parents) {
      simb::MCParticle const & mcp = mct.GetParticle(s);
    //CHEKC hardcode, TPC filter:
    if(abs(mcp.Vx())>210 ||  abs(mcp.Vy())>210||mcp.Vz()>510 || mcp.Vz()<-1){
      std::cout<<"OUTSIDE TPC x y z ="<<mcp.Vx()<<" "<<mcp.Vy()<<" "<<mcp.Vz()<<std::endl;
      exit(0);
    }
    if(abs(mcp.PdgCode()) == 2114 || abs(mcp.PdgCode()) == 2214) {
  if(ftree) FillTree(e, i, mcp.PdgCode(), true);
  return true;
      }
      else if(mcp.PdgCode() == 22) {
  in_nucleus_photons.push_back(mcp.Mother());
      }
    }

    for(size_t const s : in_nucleus_photons) {
      simb::MCParticle const & mcp = mct.GetParticle(s);
      if(abs(mcp.PdgCode()) == 2114 || abs(mcp.PdgCode()) == 2214) {
  if(ftree) FillTree(e, i, s, true);
  return true;
      }
    }

  }

  if(ftree) FillTree(e, 0, SIZE_MAX, false);
  return false;

}

NCDeltaRadiative& NCDeltaRadiative::operator= ( NCDeltaRadiative const &  )

delete

NCDeltaRadiative& NCDeltaRadiative::operator= ( NCDeltaRadiative &&  )

delete

void NCDeltaRadiative::Reset

(

)

Definition at line 102 of file NCDeltaRadiative_module.cc.

                             {

  frun = -1;
  fsubrun = -1;
  fevent = -1;
  fnu_pdg = 0;
  fccnc = -1;
  fmode = -2;
  finteraction_type = -2;
  fis_nc_delta_radiative = -1;
  fparent_status_code = -1; 
  fparent_pdg = 0;

}

Member Data Documentation

int NCDeltaRadiative::fccnc

private

Definition at line 38 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::fevent

private

Definition at line 36 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::finteraction_type

private

Definition at line 40 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::fis_nc_delta_radiative

private

Definition at line 41 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::fmode

private

Definition at line 39 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::fnu_pdg

private

Definition at line 37 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::fparent_pdg

private

Definition at line 43 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::fparent_status_code

private

Definition at line 42 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::frun

private

Definition at line 34 of file NCDeltaRadiative_module.cc.

int NCDeltaRadiative::fsubrun

private

Definition at line 35 of file NCDeltaRadiative_module.cc.

TTree* NCDeltaRadiative::ftree

private

Definition at line 32 of file NCDeltaRadiative_module.cc.

---

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

• NCDeltaRadiative_module.cc