Inheritance diagram for larg4::CheckSimEnergyDeposit:

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

Private Member Functions
void	beginJob () override

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

Private Attributes
TH1F *	_hnHits {nullptr}

TH1F *	_hEdep {nullptr}

TH1F *	_hnumPhotons {nullptr}

TH1F *	_hLandauPhotons {nullptr}

TH1F *	_hLandauEdep {nullptr}

TH1F *	_hSteplength {nullptr}

TNtuple *	_ntuple {nullptr}

Detailed Description

Definition at line 27 of file CheckSimEnergyDeposit_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 45 of file CheckSimEnergyDeposit_module.cc.

                                                                             :
   art::EDAnalyzer(p)
 {}

Member Function Documentation

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

overrideprivate

Definition at line 62 of file CheckSimEnergyDeposit_module.cc.

 {
   auto allSims = event.getMany<sim::SimEnergyDepositCollection>();
   for (auto const& sims : allSims) {
     double sumPhotons=0.0;
     double sumE = 0.0;
     _hnHits->Fill(sims->size());
     for (auto const& hit : *sims) {
       // sum up energy deposit in a 1cm slice of liquid Argon.
       if (std::abs(hit.EndZ())<0.5) {
         sumPhotons= sumPhotons + hit.NumPhotons();
         sumE= sumE +hit.Energy();
       }
       _hnumPhotons->Fill( hit.NumPhotons());
       _hEdep->Fill( hit.Energy());   // energy deposit in MeV
       _hSteplength->Fill( hit.StepLength()); // step length in cm
     }
     _hLandauPhotons->Fill(sumPhotons);
     _hLandauEdep->Fill(sumE);
   }
 } // end analyze

void larg4::CheckSimEnergyDeposit::beginJob ( )

overrideprivate

Definition at line 49 of file CheckSimEnergyDeposit_module.cc.

 {
   art::ServiceHandle<art::TFileService const> tfs;
   _hnHits = tfs->make<TH1F>("hnHits", "Number of SimEnergyDeposits", 300, 0, 0);
   _hEdep = tfs->make<TH1F>("hEdep", "Energy deposition in SimEnergyDeposits", 100,0.,0.02);
   _hnumPhotons = tfs->make<TH1F>("hnumPhotons", "number of photons per  SimEnergyDeposit", 100,0.,500.);
   _hLandauPhotons= tfs->make<TH1F>("hLandauPhotons", "number of photons/cm", 100,0.,2000000.);
   _hLandauEdep= tfs->make<TH1F>("hLandauEdep", "Edep/cm", 100,0.,10.);
   _hSteplength= tfs->make<TH1F>("hSteplength", "geant 4 step length", 100,0.,0.05);
   _ntuple = tfs->make<TNtuple>("ntuple","Demo ntuple",
                                "Event:Edep:em_Edep:nonem_Edep:xpos:ypos:zpos:time");
 } // end beginJob