larg4::LArG4Ana Class Reference

Inheritance diagram for larg4::LArG4Ana:

Public Member Functions
	LArG4Ana (fhicl::ParameterSet const &pset)
	Standard constructor and destructor for an FMWK module. More...
void	analyze (const art::Event &evt)
void	beginJob ()
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
std::string	fG4ModuleLabel
	module label for the Geant More...
std::string	fTruthModuleLabel
	module label for the Geant More...
TH1D *	fPDGCodes
TH1D *	fPi0Momentum
TH1D *	fnEnergy
TH1D *	fnDist
TH1D *	fnumChannels
	The number of channels recieving charge per event. More...
TProfile *	fnumIDEs
	Number of drift electrons per channel. More...
TH1D *	fEventCharge
	Charge collected per event. More...
TH1D *	fEventEnergy
	Energy collected per event. More...
TProfile *	fChannelCharge
	Charge per channel. More...
TProfile *	fChannelEnergy
	Energy per channel. More...
TTree *	fTree
Int_t	fTEvt
Int_t	fTSub
Int_t	fTRun
Int_t	fTPdg
Int_t	fTID
Int_t	fTNdsOriginal
Int_t	fTNds
Int_t	fTNds4
Int_t *	fTDID
Int_t *	fTDPdg
Float_t *	fTDWt
Char_t	fTProcess [35]
Char_t	fTVolume [35]
Char_t	fTTVolume [35]
Char_t	fTMaterial [35]
Char_t	fTDProcess [200][35]
Int_t	fTParentID
Int_t	fTStatus
Float_t	fTWeight
Float_t *	fT4Origin
Float_t *	fT4DOrigin
Float_t *	fT4Termination
Float_t *	fT4Momentum
Float_t *	fT4DMomentum

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 38 of file LArG4Ana_module.cc.

Constructor & Destructor Documentation

larg4::LArG4Ana::LArG4Ana ( fhicl::ParameterSet const &  pset )

explicit

Standard constructor and destructor for an FMWK module.

Definition at line 99 of file LArG4Ana_module.cc.

    : EDAnalyzer(pset)
    , fG4ModuleLabel{pset.get< std::string >("GeantModuleLabel")}
    , fTruthModuleLabel{pset.get< std::string >("TruthModuleLabel")}
    , fTNdsOriginal{pset.get< int         >("Ndaughters"      )}
    , fTNds{fTNdsOriginal}
  {}

Member Function Documentation

void larg4::LArG4Ana::analyze

(

const art::Event &

evt

)

Definition at line 183 of file LArG4Ana_module.cc.

  {

    //get the list of particles from this event
    art::ServiceHandle<cheat::ParticleInventoryService const> pi_serv;
    const sim::ParticleList& plist = pi_serv->ParticleList();
    art::ServiceHandle<geo::Geometry const> geom;

    // loop over all sim::SimChannels in the event and make sure there are no
    // sim::IDEs with trackID values that are not in the sim::ParticleList
    std::vector<const sim::SimChannel*> sccol;
    evt.getView(fG4ModuleLabel, sccol);

    double totalCharge=0.0;
    double totalEnergy=0.0;
    fnumChannels->Fill(sccol.size());
    for(size_t sc = 0; sc < sccol.size(); ++sc){
      double numIDEs=0.0;
      double scCharge=0.0;
      double scEnergy=0.0;
      const auto & tdcidemap = sccol[sc]->TDCIDEMap();
      for(auto mapitr = tdcidemap.begin(); mapitr != tdcidemap.end(); mapitr++){
        const std::vector<sim::IDE> idevec = (*mapitr).second;
        numIDEs += idevec.size();
        for(size_t iv = 0; iv < idevec.size(); ++iv){
          if(plist.find( idevec[iv].trackID ) == plist.end()
             && idevec[iv].trackID != sim::NoParticleId)
          mf::LogWarning("LArG4Ana") << idevec[iv].trackID << " is not in particle list";
          totalCharge +=idevec[iv].numElectrons;
          scCharge += idevec[iv].numElectrons;
          totalEnergy +=idevec[iv].energy;
          scEnergy += idevec[iv].energy;
        }
      }
      fnumIDEs->Fill(sc,numIDEs);
      fChannelCharge->Fill(sc,scCharge);
      fChannelEnergy->Fill(sc,scEnergy);
    }
    fEventCharge->Fill(totalCharge);
    fEventEnergy->Fill(totalEnergy);

    // get the particles from the back tracker
    const sim::ParticleList& Particles = pi_serv->ParticleList();
    std::vector<const simb::MCParticle*> pvec;
    pvec.reserve(Particles.size());
    for (const auto& PartPair: Particles) {
      pvec.push_back(PartPair.second);
      fPDGCodes->Fill(PartPair.second->PdgCode());
    }

    // now look for pi0's that decay to 2 gammas
    int pi0loc = -1;
    int numpi0gamma = 0;
    for(unsigned int i = 0; i < pvec.size(); ++i){
      if(pvec[i]->PdgCode() == 111) pi0loc = i;
      if(pvec[i]->Mother() == pi0loc+1 &&
         pi0loc > 0 &&
         pvec[i]->PdgCode() == 22){
        mf::LogInfo("LArG4Ana") << pvec[i]->E() << " gamma energy ";
        ++numpi0gamma;
      }

      // n,Lambda,K0s,K0L,K0
      if (pvec[i]->PdgCode() == 2112 ||
          pvec[i]->PdgCode() == 3122 ||
          pvec[i]->PdgCode() == 130  ||
          pvec[i]->PdgCode() == 310  ||
          pvec[i]->PdgCode() == 311 ) {
        fnEnergy->Fill(pvec[i]->E(),pvec[i]->Weight());
        fnDist->Fill(pvec[i]->Vx(),pvec[i]->Weight());
      }

      fTPdg = pvec[i]->PdgCode();
      fTID = pvec[i]->TrackId();
      // 0 out strings, else there may be cruft in here from prev evt.
      for (unsigned int s = 0; s < 35; ++s){
          *(fTProcess+s)  = 0;
          *(fTProcess+s)  = 0;
          *(fTMaterial+s) = 0;
          *(fTMaterial+s) = 0;
          *(fTVolume+s)   = 0;
          *(fTVolume+s)   = 0;
          *(fTTVolume+s)  = 0;
          *(fTTVolume+s)  = 0;
      }

      for(unsigned int s = 0; s < pvec[i]->Process().length(); ++s) *(fTProcess+s) = pvec[i]->Process()[s];

      TVector3 dum = pvec[i]->Position().Vect();

      for (unsigned int s = 0; s < geom->MaterialName(pvec[i]->Position().Vect()).length(); ++s)
        *(fTMaterial+s) = geom->MaterialName(pvec[i]->Position().Vect())[s];

      for (unsigned int s = 0; s < geom->VolumeName(pvec[i]->Position().Vect()).length(); ++s)
        *(fTVolume+s) = geom->VolumeName(pvec[i]->Position().Vect())[s];

      for (unsigned int s = 0; s < geom->VolumeName(pvec[i]->EndPosition().Vect()).length(); ++s)
        *(fTTVolume+s) = geom->VolumeName(pvec[i]->EndPosition().Vect())[s];

      fTEvt = evt.id().event();
      fTSub = evt.subRun();
      fTRun = evt.run();
      fTParentID = pvec[i]->Mother();
      fTStatus = pvec[i]->StatusCode();
      int daughter = 9999;
      fTNds = TMath::Min(pvec[i]->NumberDaughters(),fTNdsOriginal);
      for( int d = 0; d < fTNds; d++ ){
        daughter = pvec[i]->Daughter(d);
        fTDID[d] = daughter;
        // zero it out.
        for (unsigned int s = 0; s < 35; ++s) *(fTDProcess[d]+s) = 0;

        for(unsigned int jj = i; jj < pvec.size(); ++jj){ // Don't look below i.

          if (fTDID[d] == pvec[jj]->TrackId()){
            fTDPdg[d] = pvec[jj]->PdgCode(); // get the pointer,
            fTDWt[d]  = pvec[jj]->Weight();

            for (unsigned int s = 0; s < pvec[jj]->Process().length(); ++s)
              *(fTDProcess[d]+s) = pvec[jj]->Process()[s];

            for (unsigned int kk = 0; kk < 4; ++kk){
              fT4DOrigin[d*4+kk] = pvec[jj]->Position()[kk];
              fT4DMomentum[d*4+kk] = pvec[jj]->Momentum()[kk];
            }
            break;
          }
        }
      }//end loop over d

      for (unsigned int ii = 0; ii < 4; ++ii){
        fT4Termination[ii] = 1e9;
        fT4Origin[ii] = pvec[i]->Position()[ii];
        if (ii!=3) fT4Termination[ii] = pvec[i]->EndPosition()[ii];
        if (ii==4) fT4Termination[ii] = pvec[i]->Momentum()[ii]; // yes, odd
        fT4Momentum[ii] = pvec[i]->Momentum()[ii];
      }

      fTWeight = pvec[i]->Weight();
      fTree->Fill();

    } // end loop on particles in list
    if(numpi0gamma == 2 && pi0loc > 0){
      mf::LogInfo("LArG4Ana") << pvec[pi0loc]->E();
      fPi0Momentum->Fill(pvec[pi0loc]->E());
    }

    return;
  }

void larg4::LArG4Ana::beginJob ( )

virtual

Reimplemented from art::EDAnalyzer.

Definition at line 108 of file LArG4Ana_module.cc.

  {
    art::ServiceHandle<art::TFileService const> tfs;
    art::ServiceHandle<geo::Geometry const> geo;

    fPDGCodes    = tfs->make<TH1D>("pdgcodes", ";PDG Code;",               5000, -2500, 2500);
    fPi0Momentum = tfs->make<TH1D>("pi0mom",   ";#pi^{0} Momentum (GeV);", 1000, 0.,    1000.);

    fTree = tfs->make<TTree>("MCTTree","MCTTree");
    fnEnergy = tfs->make<TH1D>("nEnergy", ";n,#Lambda^{0},K^{0} Momentum (GeV);", 100, 0., 10.);
    fnDist = tfs->make<TH1D>("nDistance", ";n,#Lambda^{0},K^{0} Distance (m);", 200, -30000.0, +30000.);


    // Some histograms relating to drift electrons, active detector
    // channels and charge/energy on channels
    fnumChannels = tfs->make<TH1D>("fnumChannels",
                                   "Active channels;Active channels;# events",
                                   256, 0, geo->Nchannels());
    fnumIDEs = tfs->make<TProfile>("fnumIDEs",
                                   "Drift Electrons per channel;Channel;Drift electrons",
                                   geo->Nchannels()+1, 0, geo->Nchannels(),
                                   0, 1e4);
    fEventCharge = tfs->make<TH1D>("fEventCharge",
                                   "Charge in event;Total charge per event;# events",
                                   100, 0, 2.5e8);
    fEventEnergy = tfs->make<TH1D>("fEventEnergy",
                                   "Energy in event;Total energy per event;# events",
                                   100, 0, 1e4);
    fChannelCharge = tfs->make<TProfile>("fChannelCharge",
                                         "Charge on channel;Channel;Total charge per channel",
                                         geo->Nchannels()+1,0,geo->Nchannels(),
                                         0,1e5);
    fChannelEnergy = tfs->make<TProfile>("fChannelEnergy",
                                         "Energy on channel;Channel;Total energy per channel",
                                         geo->Nchannels()+1,0,geo->Nchannels(),
                                         0, 1e3);


    fT4Origin = new Float_t[4];
    fT4DOrigin = new Float_t[fTNds*4];
    fT4Termination = new Float_t[4];
    fT4Momentum = new Float_t[4];
    fT4DMomentum = new Float_t[fTNds*4];
    fTDID = new Int_t[fTNds];
    fTDPdg = new Int_t[fTNds];
    fTDWt  = new Float_t[fTNds];
    fTNds4 = fTNds*4; //  TTree/Branch requirement to store this.

    fTree->Branch("MCEvt", &fTEvt, "MCEvt/I");
    fTree->Branch("MCSub", &fTSub, "MCSub/I");
    fTree->Branch("MCRun", &fTRun, "MCRun/I");
    fTree->Branch("MCWt", &fTWeight, "MCWt/F");
    fTree->Branch("MCPdg", &fTPdg, "MCPdg/I");
    fTree->Branch("MCID", &fTID, "MCID/I");
    fTree->Branch("MCParentID", &fTParentID, "MCParentID/I");
    fTree->Branch("MCNumDs", &fTNds, "MCNumDs/I");
    fTree->Branch("MCNumDs4", &fTNds4, "MCNumDs4/I");
    fTree->Branch("MCDID", fTDID, "MCDID[MCNumDs]/I");
    fTree->Branch("MCDPdg", fTDPdg, "MCDPdg[MCNumDs]/I");
    fTree->Branch("MCDWt", fTDWt, "MCDWt[MCNumDs]/I");
    fTree->Branch("MCProcess", fTProcess, "MCProcess/C");
    fTree->Branch("MCVolume", fTVolume, "MCVolume/C");
    fTree->Branch("MCTVolume", fTTVolume, "MCTVolume/C");
    fTree->Branch("MCMaterial", fTMaterial, "MCMaterial/C");
    fTree->Branch("MCDProcess", fTDProcess, "MCDProcess[MCNumDs]/C");
    fTree->Branch("MCStatus", &fTStatus, "MCStatus/I");
    fTree->Branch("MCOrigin", fT4Origin, "MCOrigin[4]/F");
    fTree->Branch("MCDOrigin", fT4DOrigin, "MCDOrigin[MCNumDs4]/F");
    fTree->Branch("MCTermination", fT4Termination, "MCTermination[4]/F");
    fTree->Branch("MCMomentum", fT4Momentum, "MCMomentum[4]/F");
    fTree->Branch("MCDMomentum", fT4DMomentum, "MCDMomentum[MCNumDs4]/F");

  }

Member Data Documentation

TProfile* larg4::LArG4Ana::fChannelCharge

private

Charge per channel.

Definition at line 61 of file LArG4Ana_module.cc.

TProfile* larg4::LArG4Ana::fChannelEnergy

private

Energy per channel.

Definition at line 62 of file LArG4Ana_module.cc.

TH1D* larg4::LArG4Ana::fEventCharge

private

Charge collected per event.

Definition at line 59 of file LArG4Ana_module.cc.

TH1D* larg4::LArG4Ana::fEventEnergy

private

Energy collected per event.

Definition at line 60 of file LArG4Ana_module.cc.

std::string larg4::LArG4Ana::fG4ModuleLabel

private

module label for the Geant

Definition at line 49 of file LArG4Ana_module.cc.

TH1D* larg4::LArG4Ana::fnDist

private

Definition at line 56 of file LArG4Ana_module.cc.

TH1D* larg4::LArG4Ana::fnEnergy

private

Definition at line 55 of file LArG4Ana_module.cc.

TH1D* larg4::LArG4Ana::fnumChannels

private

The number of channels recieving charge per event.

Definition at line 57 of file LArG4Ana_module.cc.

TProfile* larg4::LArG4Ana::fnumIDEs

private

Number of drift electrons per channel.

Definition at line 58 of file LArG4Ana_module.cc.

TH1D* larg4::LArG4Ana::fPDGCodes

private

Definition at line 53 of file LArG4Ana_module.cc.

TH1D* larg4::LArG4Ana::fPi0Momentum

private

Definition at line 54 of file LArG4Ana_module.cc.

Float_t* larg4::LArG4Ana::fT4DMomentum

private

Definition at line 90 of file LArG4Ana_module.cc.

Float_t* larg4::LArG4Ana::fT4DOrigin

private

Definition at line 87 of file LArG4Ana_module.cc.

Float_t* larg4::LArG4Ana::fT4Momentum

private

Definition at line 89 of file LArG4Ana_module.cc.

Float_t* larg4::LArG4Ana::fT4Origin

private

Definition at line 86 of file LArG4Ana_module.cc.

Float_t* larg4::LArG4Ana::fT4Termination

private

Definition at line 88 of file LArG4Ana_module.cc.

Int_t* larg4::LArG4Ana::fTDID

private

Definition at line 75 of file LArG4Ana_module.cc.

Int_t* larg4::LArG4Ana::fTDPdg

private

Definition at line 76 of file LArG4Ana_module.cc.

Char_t larg4::LArG4Ana::fTDProcess[200][35]

private

Definition at line 82 of file LArG4Ana_module.cc.

Float_t* larg4::LArG4Ana::fTDWt

private

Definition at line 77 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTEvt

private

Definition at line 67 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTID

private

Definition at line 71 of file LArG4Ana_module.cc.

Char_t larg4::LArG4Ana::fTMaterial[35]

private

Definition at line 81 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTNds

private

Definition at line 73 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTNds4

private

Definition at line 74 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTNdsOriginal

private

Definition at line 72 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTParentID

private

Definition at line 83 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTPdg

private

Definition at line 70 of file LArG4Ana_module.cc.

Char_t larg4::LArG4Ana::fTProcess[35]

private

Definition at line 78 of file LArG4Ana_module.cc.

TTree* larg4::LArG4Ana::fTree

private

Definition at line 66 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTRun

private

Definition at line 69 of file LArG4Ana_module.cc.

std::string larg4::LArG4Ana::fTruthModuleLabel

private

module label for the Geant

Definition at line 50 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTStatus

private

Definition at line 84 of file LArG4Ana_module.cc.

Int_t larg4::LArG4Ana::fTSub

private

Definition at line 68 of file LArG4Ana_module.cc.

Char_t larg4::LArG4Ana::fTTVolume[35]

private

Definition at line 80 of file LArG4Ana_module.cc.

Char_t larg4::LArG4Ana::fTVolume[35]

private

Definition at line 79 of file LArG4Ana_module.cc.

Float_t larg4::LArG4Ana::fTWeight

private

Definition at line 85 of file LArG4Ana_module.cc.

---

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

• larsim/larsim/LegacyLArG4/LArG4Ana_module.cc