evgen::CosmicsGen Class Reference

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

Inheritance diagram for evgen::CosmicsGen:

Public Member Functions
	CosmicsGen (fhicl::ParameterSet const &pset)
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 Member Functions
void	produce (art::Event &evt) override
void	beginJob () override
void	beginRun (art::Run &run) override

Private Attributes
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 *	fPhotonsInCStat
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 *	fElectronsInCStat
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 *	fMuonsInCStat
TH1F *	fMuonsInTPC
CLHEP::HepRandomEngine &	fEngine
	art-managed random-number engine More...
evgb::CRYHelper	fCRYHelp
	CRY generator object. 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 36 of file CosmicsGen_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 91 of file CosmicsGen_module.cc.

    : art::EDProducer{pset}
    , fbuffbox{pset.get<std::vector<double>>("BufferBox",{0.0, 0.0, 0.0, 0.0, 0.0, 0.0})}
    // create a default random engine; obtain the random seed from NuRandomService,
    // unless overridden in configuration with key "Seed"
    , fEngine(art::ServiceHandle<rndm::NuRandomService>()->createEngine(*this, pset, "Seed"))
    , fCRYHelp{pset,
               fEngine,
               art::ServiceHandle<geo::Geometry const>{}->GetWorldVolumeName()}
  {
    produces< std::vector<simb::MCTruth> >();
    produces< sumdata::RunData, art::InRun >();
  }

Member Function Documentation

void evgen::CosmicsGen::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDProducer.

Definition at line 106 of file CosmicsGen_module.cc.

  {
    art::ServiceHandle<art::TFileService const> 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);
    fPhotonsInCStat   = tfs->make<TH1F>("fPhotonsInCryostat", ";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);
    fElectronsInCStat   = tfs->make<TH1F>("fElectronsInCryotat",  ";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);
    fMuonsInCStat   = tfs->make<TH1F>("fMuonsInCryostat", ";Number Muons;Samples", 100, 0, 1000);
    fMuonsInTPC     = tfs->make<TH1F>("fMuonsInTPC",      ";Number Muons;Samples", 100, 0, 1000);
  }

void evgen::CosmicsGen::beginRun ( art::Run &  run )

overrideprivatevirtual

Reimplemented from art::EDProducer.

Definition at line 142 of file CosmicsGen_module.cc.

  {
    art::ServiceHandle<geo::Geometry const> geo;
    run.put(std::make_unique<sumdata::RunData>(geo->DetectorName()));
  }

void evgen::CosmicsGen::produce ( art::Event &  evt )

overrideprivatevirtual

Implements art::EDProducer.

Definition at line 149 of file CosmicsGen_module.cc.

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

    // fill some histograms about this event
    art::ServiceHandle<geo::Geometry const> geom;

    int nCrossCryostat = 0;

    simb::MCTruth truth;

    while(nCrossCryostat < 1){

      simb::MCTruth pretruth;
      truth.SetOrigin(simb::kCosmicRay);
      fCRYHelp.Sample(pretruth,
                       geom->SurfaceY(),
                       geom->DetLength(),
                       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& v4 = particle.Position();
        const TLorentzVector& p4 = particle.Momentum();
        double x0[3] = {v4.X(),  v4.Y(),  v4.Z() };
        double dx[3] = {p4.Px(), p4.Py(), p4.Pz()};

        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;
        }

        // now check if the particle goes through any cryostat in the detector
        // if so, add it to the truth object.
        for(unsigned int c = 0; c < geom->Ncryostats(); ++c){

          double bounds[6] = {0.};
          geom->CryostatBoundaries(bounds, c);

          //add a buffer box around the cryostat bounds to increase the acceptance
          //(geometrically) at the CRY level to make up for particles we will loose
          //due to multiple scattering effects that pitch in during GEANT4 tracking
          //By default, the buffer box has zero size
          for (unsigned int cb=0; cb<6; cb++)
             bounds[cb] = bounds[cb]+fbuffbox[cb];

          //calculate the intersection point with each cryostat surface
          bool intersects_cryo = false;
          for (int bnd=0; bnd!=6; ++bnd) {
            if (bnd<2) {
              double p2[3] = {bounds[bnd],  x0[1] + (dx[1]/dx[0])*(bounds[bnd] - x0[0]), x0[2] + (dx[2]/dx[0])*(bounds[bnd] - x0[0])};
              if ( p2[1] >= bounds[2] && p2[1] <= bounds[3] &&
                   p2[2] >= bounds[4] && p2[2] <= bounds[5] ) {
                intersects_cryo = true;
                break;
              }
            }
            else if (bnd>=2 && bnd<4) {
              double p2[3] = {x0[0] + (dx[0]/dx[1])*(bounds[bnd] - x0[1]), bounds[bnd], x0[2] + (dx[2]/dx[1])*(bounds[bnd] - x0[1])};
              if ( p2[0] >= bounds[0] && p2[0] <= bounds[1] &&
                   p2[2] >= bounds[4] && p2[2] <= bounds[5] ) {
                intersects_cryo = true;
                break;
              }
            }
            else if (bnd>=4) {
              double p2[3] = {x0[0] + (dx[0]/dx[2])*(bounds[bnd] - x0[2]), x0[1] + (dx[1]/dx[2])*(bounds[bnd] - x0[2]), bounds[bnd]};
              if ( p2[0] >= bounds[0] && p2[0] <= bounds[1] &&
                   p2[1] >= bounds[2] && p2[1] <= bounds[3] ) {
                intersects_cryo = true;
                break;
              }
            }
          }


          if (intersects_cryo) {
            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;

            //The following code no longer works now that we require intersection with the cryostat boundary
            //For example, the particle could intersect this cryostat but miss its TPC, but intersect a TPC
            //in another cryostat
            /*try{
              unsigned int tpc   = 0;
              unsigned int cstat = 0;
              geom->PositionToTPC(x2, tpc, cstat);
              if      (std::abs(particle.PdgCode())==13) ++tpcMuons;
              else if (std::abs(particle.PdgCode())==22) ++tpcPhotons;
              else if (std::abs(particle.PdgCode())==11) ++tpcElectrons;
            }
            catch(cet::exception &e){
              MF_LOG_DEBUG("CosmicsGen") << "current particle does not go through any tpc";
            }*///

            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
            break; //leave loop over cryostats to avoid adding particle multiple times
          }// end if particle goes into a cryostat
        }// end loop over cryostats in the detector

      }// loop on particles

      nCrossCryostat = truth.NParticles();

      fPhotonsPerSample  ->Fill(allPhotons);
      fElectronsPerSample->Fill(allElectrons);
      fMuonsPerSample    ->Fill(allMuons);

      fPhotonsInCStat  ->Fill(numPhotons);
      fElectronsInCStat->Fill(numElectrons);
      fMuonsInCStat    ->Fill(numMuons);

      /*fPhotonsInTPC  ->Fill(tpcPhotons);
      fElectronsInTPC->Fill(tpcElectrons);
      fMuonsInTPC    ->Fill(tpcMuons);*/
    }

    truthcol->push_back(truth);
    evt.put(std::move(truthcol));
  }// end produce

Member Data Documentation

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

private

Definition at line 46 of file CosmicsGen_module.cc.

evgb::CRYHelper evgen::CosmicsGen::fCRYHelp

private

CRY generator object.

Definition at line 84 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fElectronAngles

private

Electron rate vs angle.

Definition at line 60 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fElectronAnglesHi

private

Electron rate vs angle, high momenta.

Definition at line 63 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fElectronAnglesLo

private

Electron rate vs angle, low momenta.

Definition at line 61 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fElectronAnglesMi

private

Electron rate vs angle, middle momenta.

Definition at line 62 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fElectronCosQ

private

Electron rate vs cos(Q)

Definition at line 64 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fElectronEnergy

private

Electron energy (GeV)

Definition at line 65 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fElectronsInCStat

private

number of electrons in the cryostat during the sampled time window

Definition at line 67 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fElectronsInTPC

private

number of electrons in the tpc during the sampled time window

Definition at line 69 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fElectronsPerSample

private

number of electrons in the sampled time window

Definition at line 66 of file CosmicsGen_module.cc.

CLHEP::HepRandomEngine& evgen::CosmicsGen::fEngine

private

art-managed random-number engine

Definition at line 83 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fMuonAngles

private

Muon rate vs angle.

Definition at line 72 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fMuonAnglesHi

private

Muon rate vs angle, high momenta.

Definition at line 75 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fMuonAnglesLo

private

Muon rate vs angle, low momenta.

Definition at line 73 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fMuonAnglesMi

private

Muon rate vs angle, middle momenta.

Definition at line 74 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fMuonCosQ

private

Muon rate vs cos(Q)

Definition at line 76 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fMuonEnergy

private

Muon energy (GeV)

Definition at line 77 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fMuonsInCStat

private

number of muons in the cryostat during the sampled time window

Definition at line 79 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fMuonsInTPC

private

number of muons in the tpc during the sampled time window

Definition at line 81 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fMuonsPerSample

private

number of muons in the sampled time window

Definition at line 78 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fPhotonAngles

private

Photon rate vs angle.

Definition at line 48 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fPhotonAnglesHi

private

Photon rate vs angle, high momenta.

Definition at line 51 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fPhotonAnglesLo

private

Photon rate vs angle, low momenta.

Definition at line 49 of file CosmicsGen_module.cc.

TH2F* evgen::CosmicsGen::fPhotonAnglesMi

private

Photon rate vs angle, middle momenta.

Definition at line 50 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fPhotonCosQ

private

Photon rate vs cos(Q)

Definition at line 52 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fPhotonEnergy

private

Photon energy (GeV)

Definition at line 53 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fPhotonsInCStat

private

number of photons in the cryostat during the sampled time window

Definition at line 55 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fPhotonsInTPC

private

number of photons in the tpc during the sampled time window

Definition at line 57 of file CosmicsGen_module.cc.

TH1F* evgen::CosmicsGen::fPhotonsPerSample

private

number of photons in the sampled time window

Definition at line 54 of file CosmicsGen_module.cc.

---

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

• larsim/larsim/EventGenerator/CRY/CosmicsGen_module.cc