gar::evgen::CosmicsGen Class Reference

A module to check the results from the Monte Carlo generator. More...

Inheritance diagram for gar::evgen::CosmicsGen:

Public Member Functions
	CosmicsGen (fhicl::ParameterSet const &pset)
virtual	~CosmicsGen ()
void	produce (::art::Event &evt)
void	beginJob ()
void	beginRun (::art::Run &run)
void	reconfigure (fhicl::ParameterSet const &p)
Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDProducer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept
	Producer (fhicl::ParameterSet const &)
	Producer (Producer const &)=delete
	Producer (Producer &&)=delete
Producer &	operator= (Producer const &)=delete
Producer &	operator= (Producer &&)=delete
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::Modifier
	~Modifier () noexcept
	Modifier ()
	Modifier (Modifier const &)=delete
	Modifier (Modifier &&)=delete
Modifier &	operator= (Modifier const &)=delete
Modifier &	operator= (Modifier &&)=delete
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
evgb::CRYHelper *	fCRYHelp
	CRY generator object. More...
std::vector< double >	fbuffbox
TH2F *	fPhotonAngles
	Photon rate vs angle. More...
TH2F *	fPhotonAnglesLo
	Photon rate vs angle, low momenta. More...
TH2F *	fPhotonAnglesMi
	Photon rate vs angle, middle momenta. More...
TH2F *	fPhotonAnglesHi
	Photon rate vs angle, high momenta. More...
TH1F *	fPhotonCosQ
	Photon rate vs cos(Q) More...
TH1F *	fPhotonEnergy
	Photon energy (GeV) More...
TH1F *	fPhotonsPerSample
	number of photons in the sampled time window More...
TH1F *	fPhotonsInTPC
TH2F *	fElectronAngles
	Electron rate vs angle. More...
TH2F *	fElectronAnglesLo
	Electron rate vs angle, low momenta. More...
TH2F *	fElectronAnglesMi
	Electron rate vs angle, middle momenta. More...
TH2F *	fElectronAnglesHi
	Electron rate vs angle, high momenta. More...
TH1F *	fElectronCosQ
	Electron rate vs cos(Q) More...
TH1F *	fElectronEnergy
	Electron energy (GeV) More...
TH1F *	fElectronsPerSample
	number of electrons in the sampled time window More...
TH1F *	fElectronsInTPC
TH2F *	fMuonAngles
	Muon rate vs angle. More...
TH2F *	fMuonAnglesLo
	Muon rate vs angle, low momenta. More...
TH2F *	fMuonAnglesMi
	Muon rate vs angle, middle momenta. More...
TH2F *	fMuonAnglesHi
	Muon rate vs angle, high momenta. More...
TH1F *	fMuonCosQ
	Muon rate vs cos(Q) More...
TH1F *	fMuonEnergy
	Muon energy (GeV) More...
TH1F *	fMuonsPerSample
	number of muons in the sampled time window More...
TH1F *	fMuonsInTPC
rndm::NuRandomService::seed_t	fSeed
	override seed with a fcl parameter not equal to zero More...

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >
Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >
Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)
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

A module to check the results from the Monte Carlo generator.

Definition at line 46 of file CosmicsGen_module.cc.

Constructor & Destructor Documentation

gar::evgen::CosmicsGen::CosmicsGen ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 101 of file CosmicsGen_module.cc.

                                                        :
      art::EDProducer{pset}, fCRYHelp(0)
    {

      //the buffer box bounds specified here will extend on the cryostat boundaries
      fbuffbox = pset.get< std::vector<double> >("BufferBox",{0.0, 0.0, 0.0, 0.0, 0.0, 0.0});
      fSeed          = pset.get< rndm::NuRandomService::seed_t >("Seed",0);
      this->reconfigure(pset);

      produces< std::vector<simb::MCTruth> >();
      produces< gar::sumdata::RunData, ::art::InRun >();
    }

gar::evgen::CosmicsGen::~CosmicsGen ( )

virtual

Definition at line 115 of file CosmicsGen_module.cc.

    {
      if(fCRYHelp) delete fCRYHelp;
    }

Member Function Documentation

void gar::evgen::CosmicsGen::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 138 of file CosmicsGen_module.cc.

    {
      ::art::ServiceHandle<::art::TFileService> tfs;

      fPhotonAngles       = tfs->make<TH2F>("fPhotonAngles",        ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fPhotonAnglesLo     = tfs->make<TH2F>("fPhotonAnglesLo",      ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fPhotonAnglesMi     = tfs->make<TH2F>("fPhotonAnglesMi",      ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fPhotonAnglesHi     = tfs->make<TH2F>("fPhotonAnglesHi",      ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fPhotonCosQ         = tfs->make<TH1F>("fPhotonCosQ",          ";cos#theta;tracks",         50,-1.0,1.0);
      fPhotonEnergy       = tfs->make<TH1F>("fPhotonEnergy",        ";E (GeV)",                  5000,0.0,1000.0);
      fPhotonsPerSample   = tfs->make<TH1F>("fPhotonsPerSample",    ";Number Photons;Samples",   100, 0, 1000);
      fPhotonsInTPC       = tfs->make<TH1F>("fPhotonsInTPC",        ";Number Photons;Samples",   100, 0, 1000);

      fElectronAngles     = tfs->make<TH2F>("fElectronAngles",      ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fElectronAnglesLo   = tfs->make<TH2F>("fElectronAnglesLo",    ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fElectronAnglesMi   = tfs->make<TH2F>("fElectronAnglesMi",    ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fElectronAnglesHi   = tfs->make<TH2F>("fElectronAnglesHi",    ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fElectronCosQ       = tfs->make<TH1F>("fElectronCosQ",        ";cos#theta;tracks",         50,-1.0,1.0);
      fElectronEnergy     = tfs->make<TH1F>("fElectronEnergy",      ";E (GeV)",                  5000,0.0,1000.0);
      fElectronsPerSample = tfs->make<TH1F>("fElectronsPerSample",  ";Number Electrons;Samples", 100, 0, 1000);
      fElectronsInTPC     = tfs->make<TH1F>("fElectronsInTPC",      ";Number Electrons;Samples", 100, 0, 1000);

      fMuonAngles         = tfs->make<TH2F>("fMuonAngles",          ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fMuonAnglesLo       = tfs->make<TH2F>("fMuonAnglesLo",        ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fMuonAnglesMi       = tfs->make<TH2F>("fMuonAnglesMi",        ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fMuonAnglesHi       = tfs->make<TH2F>("fMuonAnglesHi",        ";#phi;cos#theta",           36,-180.0,180.0,50,-1.0,1.0);
      fMuonCosQ           = tfs->make<TH1F>("fMuonCosQ",            ";cos#theta;tracks",         50,-1.0,1.0);
      fMuonEnergy         = tfs->make<TH1F>("fMuonEnergy",          ";E (GeV)",                  5000,0.0,1000.0);
      fMuonsPerSample     = tfs->make<TH1F>("fMuonsPerSample",      ";Number Muons;Samples",     100, 0, 1000);
      fMuonsInTPC         = tfs->make<TH1F>("fMuonsInTPC",          ";Number Muons;Samples",     100, 0, 1000);

    }

void gar::evgen::CosmicsGen::beginRun

(

::art::Run &

run

)

Definition at line 172 of file CosmicsGen_module.cc.

    {
      // grab the geometry object to see what geometry we are using
      auto geo = gar::providerFrom<geo::GeometryGAr>();

      std::unique_ptr<gar::sumdata::RunData> runcol(new gar::sumdata::RunData(geo->DetectorName()));

      run.put(std::move(runcol));

      return;
    }

void gar::evgen::CosmicsGen::produce

(

::art::Event &

evt

)

Definition at line 185 of file CosmicsGen_module.cc.

    {
      std::unique_ptr< std::vector<simb::MCTruth> > truthcol(new std::vector<simb::MCTruth>);

      // fill some histograms about this event
      auto geom = gar::providerFrom<geo::GeometryGAr>();

      // int nCrossCryostat = 0;   // no cryostat here.

      simb::MCTruth truth;

      // don't need to have a particle on every event

      //while(nCrossCryostat < 1){

        simb::MCTruth pretruth;
        truth.SetOrigin(simb::kCosmicRay);
        //std::cout << " calling CRY Helper: " << geom->SurfaceY() << " " << geom->TPCLength() << std::endl;
        fCRYHelp->Sample(pretruth,
                         geom->SurfaceY(),
                         geom->TPCLength(),
                         0);

        int numPhotons   = 0;
        int numElectrons = 0;
        int numMuons     = 0;
        int allPhotons   = 0;
        int allElectrons = 0;
        int allMuons     = 0;

          // loop over particles in the truth object
        for(int i = 0; i < pretruth.NParticles(); ++i){
          simb::MCParticle particle = pretruth.GetParticle(i);
          const TLorentzVector& p4 = particle.Momentum();

          if      (std::abs(particle.PdgCode())==13) ++allMuons;
          else if (std::abs(particle.PdgCode())==22) ++allPhotons;
          else if (std::abs(particle.PdgCode())==11) ++allElectrons;

          TH1F* hCosQ     = 0;
          TH2F* hAngles   = 0;
          TH2F* hAnglesLo = 0;
          TH2F* hAnglesMi = 0;
          TH2F* hAnglesHi = 0;
          TH1F* hEnergy   = 0;
          if (std::abs(particle.PdgCode())==13) {
            hCosQ     = fMuonCosQ;
            hAngles   = fMuonAngles;
            hAnglesLo = fMuonAnglesLo;
            hAnglesMi = fMuonAnglesMi;
            hAnglesHi = fMuonAnglesHi;
            hEnergy   = fMuonEnergy;
          }
          else if (std::abs(particle.PdgCode())==22) {
            hCosQ     = fPhotonCosQ;
            hAngles   = fPhotonAngles;
            hAnglesLo = fPhotonAnglesLo;
            hAnglesMi = fPhotonAnglesMi;
            hAnglesHi = fPhotonAnglesHi;
            hEnergy   = fPhotonEnergy;
          }
          else if (std::abs(particle.PdgCode())==11) {
            hCosQ     = fElectronCosQ;
            hAngles   = fElectronAngles;
            hAnglesLo = fElectronAnglesLo;
            hAnglesMi = fElectronAnglesMi;
            hAnglesHi = fElectronAnglesHi;
            hEnergy   = fElectronEnergy;
          }

          // TODO: check if the particle goes through the detector, if it does
          // add it to the list

          truth.Add(particle);

          if      (std::abs(particle.PdgCode())==13) ++numMuons;
          else if (std::abs(particle.PdgCode())==22) ++numPhotons;
          else if (std::abs(particle.PdgCode())==11) ++numElectrons;

          if (hCosQ!=0) {
            double cosq = -p4.Py()/p4.P();
            double phi  = std::atan2(p4.Pz(),p4.Px());
            phi *= 180/M_PI;
            hCosQ->Fill(cosq);
            hAngles->Fill(phi,cosq);
            if      (p4.E()<1.0)  hAnglesLo->Fill(phi,cosq);
            else if (p4.E()<10.0) hAnglesMi->Fill(phi,cosq);
            else                  hAnglesHi->Fill(phi,cosq);
            hEnergy->Fill(p4.E());
          }//end if there is a cos(theta) histogram

        }// loop on particles

        //}

      truthcol->push_back(truth);
      evt.put(std::move(truthcol));

      return;
    }// end produce

void gar::evgen::CosmicsGen::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 121 of file CosmicsGen_module.cc.

    {
      if(fCRYHelp){
        delete fCRYHelp;
        fCRYHelp = 0;
      }

      CLHEP::HepRandomEngine &engine = art::ServiceHandle<rndm::NuRandomService>()->createEngine(*this,p,"Seed");

      auto geo = gar::providerFrom<gar::geo::GeometryGAr>();

      fCRYHelp = new evgb::CRYHelper(p, engine, geo->GetWorldVolumeName());

      return;
    }

Member Data Documentation

std::vector<double> gar::evgen::CosmicsGen::fbuffbox

private

Definition at line 61 of file CosmicsGen_module.cc.

evgb::CRYHelper* gar::evgen::CosmicsGen::fCRYHelp

private

CRY generator object.

Definition at line 59 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fElectronAngles

private

Electron rate vs angle.

Definition at line 73 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fElectronAnglesHi

private

Electron rate vs angle, high momenta.

Definition at line 76 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fElectronAnglesLo

private

Electron rate vs angle, low momenta.

Definition at line 74 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fElectronAnglesMi

private

Electron rate vs angle, middle momenta.

Definition at line 75 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fElectronCosQ

private

Electron rate vs cos(Q)

Definition at line 77 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fElectronEnergy

private

Electron energy (GeV)

Definition at line 78 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fElectronsInTPC

private

number of electrons in the tpc during the sampled time window

Definition at line 80 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fElectronsPerSample

private

number of electrons in the sampled time window

Definition at line 79 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fMuonAngles

private

Muon rate vs angle.

Definition at line 83 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fMuonAnglesHi

private

Muon rate vs angle, high momenta.

Definition at line 86 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fMuonAnglesLo

private

Muon rate vs angle, low momenta.

Definition at line 84 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fMuonAnglesMi

private

Muon rate vs angle, middle momenta.

Definition at line 85 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fMuonCosQ

private

Muon rate vs cos(Q)

Definition at line 87 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fMuonEnergy

private

Muon energy (GeV)

Definition at line 88 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fMuonsInTPC

private

number of muons in the tpc during the sampled time window

Definition at line 90 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fMuonsPerSample

private

number of muons in the sampled time window

Definition at line 89 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fPhotonAngles

private

Photon rate vs angle.

Definition at line 63 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fPhotonAnglesHi

private

Photon rate vs angle, high momenta.

Definition at line 66 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fPhotonAnglesLo

private

Photon rate vs angle, low momenta.

Definition at line 64 of file CosmicsGen_module.cc.

TH2F* gar::evgen::CosmicsGen::fPhotonAnglesMi

private

Photon rate vs angle, middle momenta.

Definition at line 65 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fPhotonCosQ

private

Photon rate vs cos(Q)

Definition at line 67 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fPhotonEnergy

private

Photon energy (GeV)

Definition at line 68 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fPhotonsInTPC

private

number of photons in the tpc during the sampled time window

Definition at line 70 of file CosmicsGen_module.cc.

TH1F* gar::evgen::CosmicsGen::fPhotonsPerSample

private

number of photons in the sampled time window

Definition at line 69 of file CosmicsGen_module.cc.

rndm::NuRandomService::seed_t gar::evgen::CosmicsGen::fSeed

private

override seed with a fcl parameter not equal to zero

Definition at line 93 of file CosmicsGen_module.cc.

---

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

• garsoft/EventGenerator/CRY/CosmicsGen_module.cc