module to produce single or multiple specified particles in the detector More...

Inheritance diagram for evgen::FileMuons:

Public Member Functions
	FileMuons (fhicl::ParameterSet const &pset)

void	produce (art::Event &evt)

void	beginJob ()

void	beginRun (art::Run &run)

void	endJob ()

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	ReadEvents (simb::MCTruth &mct)

Private Attributes
int	fEventNumberOffset

std::vector< int >	fPDG

std::vector< double >	fXYZ_Off

std::string	fFileName

std::string	fMuonsFileType

std::string	fTreeName

std::vector< std::string >	fBranchNames

std::ifstream *	fMuonFile

TFile *	fMuonFileR

TTree *	TNtuple

unsigned int	countFile

Float_t	xtmp

Float_t	ytmp

Float_t	ztmp

Float_t	pxtmp

Float_t	pytmp

Float_t	pztmp

Float_t	charge

Float_t	E

Float_t	costheta

Float_t	phi

Float_t	xdet

Float_t	ydet

Float_t	zdet

TBranch *	b_x

TBranch *	b_y

TBranch *	b_z

TBranch *	b_E

TBranch *	b_costheta

TBranch *	b_phi

TBranch *	b_xdet

TBranch *	b_ydet

TBranch *	b_zdet

TBranch *	b_px

TBranch *	b_py

TBranch *	b_pz

TBranch *	b_charge

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

module to produce single or multiple specified particles in the detector

Definition at line 44 of file FileMuons_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 105 of file FileMuons_module.cc.

      : EDProducer{pset}
      , fEventNumberOffset(pset.get<      int                 >("EventNumberOffset"))
      , fPDG              (pset.get< std::vector<int>         >("PDG")              )
      , fXYZ_Off          (pset.get< std::vector<double>      >("InitialXYZOffsets"))
      , fFileName         (pset.get< std::string              >("FileName")         )
      , fMuonsFileType    (pset.get< std::string              >("MuonsFileType")    )
      , fTreeName         (pset.get< std::string              >("TreeName")         )
      , fBranchNames      (pset.get< std::vector<std::string> >("BranchNames")      )
   {
 
     produces< std::vector<simb::MCTruth> >();
     produces< sumdata::RunData, art::InRun >();
 
   }

Member Function Documentation

void evgen::FileMuons::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 122 of file FileMuons_module.cc.

   {
     countFile=fEventNumberOffset;
     mf::LogInfo("FileMuons : starting at event ")  << countFile <<std::endl;
 
     if (fMuonsFileType.compare("source")==0)
       {
         std::cout << "FileMuons: Not yet equipped to walk through muons with TFS mojo."<< std::endl;
       }
     else if (fMuonsFileType.compare("root")==0)
       {
         std::cout << "FileMuons: You have chosen to read muons from Root File " << fFileName << std::endl;
       }
     else if (fMuonsFileType.compare("text")==0)
       {
         std::cout << "FileMuons: You have chosen to read muons from " << fFileName << "." << std::endl;
       }
     else
       {
         std::cout << "FileMuons: You must specify one of source/text/root file to read for muons."<< std::endl;
 
       }
 
     if (fMuonsFileType.compare("text")==0)
       {
         fMuonFile = new std::ifstream(fFileName.c_str());
         long begin = fMuonFile->tellg();
         fMuonFile->seekg (0, std::ios::end);
         long end = fMuonFile->tellg();
         std::cout << "FileMuons: " << fFileName << " size is: " << (end-begin) << " bytes.\n";
         fMuonFile->seekg (0, std::ios::beg);
 
         for (unsigned int header=0; header<3 && fMuonFile->good(); ++header)
           {
             std::string line;
             getline(*fMuonFile,line);
           }
         if (!fMuonFile->good())
           {
             std::cout << "FileMuons: Problem reading muon file header."<< std::endl;
           }
       }  // fMuonsFileType is a text file.
     else if (fMuonsFileType.compare("root")==0)
       {
         fMuonFileR = new TFile(fFileName.c_str(),"READ");
         TNtuple = (TTree*)(fMuonFileR->Get(fTreeName.c_str()));
 
         TNtuple->SetBranchAddress("x", &xtmp, &b_x);
         TNtuple->SetBranchAddress("y", &ytmp, &b_y);
         TNtuple->SetBranchAddress("z", &ztmp, &b_z);
         TNtuple->SetBranchAddress("E", &E, &b_E);
         TNtuple->SetBranchAddress("costheta", &costheta, &b_costheta);
         TNtuple->SetBranchAddress("phi", &phi, &b_phi);
         TNtuple->SetBranchAddress("xdet", &xdet, &b_xdet);
         TNtuple->SetBranchAddress("ydet", &ydet, &b_ydet);
         TNtuple->SetBranchAddress("zdet", &zdet, &b_zdet);
         TNtuple->SetBranchAddress("px", &pxtmp, &b_px);
         TNtuple->SetBranchAddress("py", &pytmp, &b_py);
         TNtuple->SetBranchAddress("pz", &pztmp, &b_pz);
         TNtuple->SetBranchAddress("charge", &charge, &b_charge);
 
       }  // fMuonsFileType is a root file.
 
   }

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

virtual

Reimplemented from art::EDProducer.

Definition at line 197 of file FileMuons_module.cc.

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

void evgen::FileMuons::endJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 188 of file FileMuons_module.cc.

   {
     if (fMuonsFileType.compare("text")==0)
       fMuonFile->close();
     if (fMuonsFileType.compare("root")==0)
       fMuonFileR->Close();
   }

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

virtual

unique_ptr allows ownership to be transferred to the art::Event after the put statement

Implements art::EDProducer.

Definition at line 204 of file FileMuons_module.cc.

   {
 
     ///unique_ptr allows ownership to be transferred to the art::Event after the put statement
     std::unique_ptr< std::vector<simb::MCTruth> > truthcol(new std::vector<simb::MCTruth>);
 
     simb::MCTruth truth;
     truth.SetOrigin(simb::kSingleParticle);
     ReadEvents(truth);
 
 //     std::cout << "put mctruth into the vector" << std::endl;
     truthcol->push_back(truth);
 
 //     std::cout << "add vector to the event " << truthcol->size() << std::endl;
     evt.put(std::move(truthcol));
 
     return;
   }

void evgen::FileMuons::ReadEvents ( simb::MCTruth & mct )

private

every event will have one of each particle species in the fPDG array

Definition at line 224 of file FileMuons_module.cc.

   {
 
 //     std::cout << "size of particle vector is " << fPDG.size() << std::endl;
 
     ///every event will have one of each particle species in the fPDG array
     for (unsigned int i=0; i<fPDG.size(); ++i) {
 
       // Choose momentum
       //double p = 0.0;
       double m(0.108);
 
       TVector3 x;
       TVector3 p;
       Double_t q = 0.;
       Int_t pdgLocal;
 
       if (fMuonsFileType.compare("text")==0)
         {
 
           std::string line;
           getline(*fMuonFile,line);
           if (!fMuonFile->good())
             {
               std::cout << "FileMuons: Problem reading muon file line ...."<< countFile << ". Perhaps you've exhausted the events in " << fFileName << std::endl; exit(0);
             }
           else
             {
               //          std::cout << "FileMuons: getline() gives "<< line << " for event " << countFile << std::endl;
             }
           countFile++;
 
           MF_LOG_DEBUG("FileMuons: countFile is ") << countFile <<std::endl;
           char * cstr, *ptok;
 
       // Split this line into tokens
           cstr = new char [line.size()+1];
           strcpy (cstr, line.c_str());
           // cstr now contains a c-string copy of str
           ptok=strtok (cstr,"*");
           unsigned int fieldCount = 0;
           unsigned int posIndex = 0;
           unsigned int pIndex = 0;
           while (ptok!=NULL)
             {
               // std::cout << ptok << std::endl;
               ptok=strtok(NULL,"*");
               if (fieldCount==9 || fieldCount==10 || fieldCount==11)
                 {
                   p[pIndex] = atof(ptok); pIndex++;
                   //   std::cout << ptok << std::endl;
                 }
               if (fieldCount==6 || fieldCount==7 || fieldCount==8)
                 {
                   x[posIndex] = atof(ptok);
                   // make the z axis point up for x, as with p
                   if (posIndex==2) {x[posIndex] = -1.0*x[posIndex];}
                   posIndex++;
                 }
               if (fieldCount==12)
                 {
                   q = atof(ptok);
                 }
               fieldCount++;
             }
 
           delete[] cstr;
 
         }
       else if (fMuonsFileType.compare("root")==0) // from root file
         {
             /*
               // Don't use this yet. Keep the specific branch-by-branch identification.
               for (unsigned int ii=0;ii<fBranchNames.size();ii++)
               {
                TNtuple->SetBranchAddress(fBranchNames[ii], x+ii);
               }
             */
           //      TNtuple->ResetBranchAddresses();
           TNtuple->GetEntry(countFile);
           //TNtuple->Show(countFile);
 
           x.SetXYZ(xdet,ydet,-zdet); // as with txt file, make z point up.
           // Watch for units change to mm in Modern JdJ files!!
           // This is for pre Spring-2012 JdJ Ntuples.
           p.SetXYZ(pxtmp,pytmp,pztmp);
           q = charge;
 
           countFile++;
 
         } // End read.
 
       static TDatabasePDG  pdgt;
       pdgLocal = -q*fPDG[i];
 
       TParticlePDG* pdgp = pdgt.GetParticle(pdgLocal);
       if (pdgp) m = pdgp->Mass();
 
 
       //       std::cout << "set the position "<<std::endl;
       // This gives coordinates at the center of the 300mx300m plate that is 3m above top of
       // cavern. Got these by histogramming deJong's xdet,ydet,zdet.
       const double cryoGap = 15.0;
       x[0] -= fXYZ_Off[0];
       x[1] -= fXYZ_Off[1];
       x[2] -= fXYZ_Off[2]; // 3 for plate height above top of cryostat.
       // Now, must rotate to TPC coordinates. Let's orient TPC axis along z axis,
       // Cosmics, mostly going along deJong's +z axis must be going along TPC -y axis.
       x.RotateX(-M_PI/2);
       p.RotateX(-M_PI/2);
       //add vector of the position of the center of the point between Cryostats
       // level with top. (To which I've added 3m - in above code - in height.)
       // This is referenced from origin at center-right of first cryostat.
       art::ServiceHandle<geo::Geometry const> geom;
       TVector3 off3(geom->CryostatHalfWidth()*0.01,geom->CryostatHalfHeight()*0.01,geom->CryostatLength()*0.01+cryoGap*0.01/2.0) ;
       x += off3;
 
       TLorentzVector pos(x[0]*100.0, x[1]*100.0, x[2]*100.0, 0.0);
       TLorentzVector pvec(p[0]*1000.0,p[1]*1000.0,p[2]*1000.0,std::sqrt(p.Mag2()*1000.0*1000.0+m*m));
       std::cout << "x[m] and p [TeV] are " << std::endl;
       x.Print();
       p.Print();
 
       int trackid = -1*(i+1); // set track id to -i as these are all primary particles and have id <= 0
       std::string primary("primary");
       simb::MCParticle part(trackid, pdgLocal, primary);
       part.AddTrajectoryPoint(pos, pvec);
 
       //       std::cout << "add the particle to the primary" << std::endl;
 
       mct.Add(part);
 
     }//end loop over particles
 
     return;
   }

Member Data Documentation

TBranch* evgen::FileMuons::b_charge

private

Definition at line 94 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_costheta

private

Definition at line 86 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_E

private

Definition at line 85 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_phi

private

Definition at line 87 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_px

private

Definition at line 91 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_py

private

Definition at line 92 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_pz

private

Definition at line 93 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_x

private

Definition at line 82 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_xdet

private

Definition at line 88 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_y

private

Definition at line 83 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_ydet

private

Definition at line 89 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_z

private

Definition at line 84 of file FileMuons_module.cc.

TBranch* evgen::FileMuons::b_zdet

private

Definition at line 90 of file FileMuons_module.cc.

Float_t evgen::FileMuons::charge

private

Definition at line 74 of file FileMuons_module.cc.

Float_t evgen::FileMuons::costheta

private

Definition at line 76 of file FileMuons_module.cc.

unsigned int evgen::FileMuons::countFile

private

Definition at line 70 of file FileMuons_module.cc.

Float_t evgen::FileMuons::E

private

Definition at line 75 of file FileMuons_module.cc.

std::vector<std::string> evgen::FileMuons::fBranchNames

private

Definition at line 65 of file FileMuons_module.cc.

int evgen::FileMuons::fEventNumberOffset

private

Definition at line 59 of file FileMuons_module.cc.

std::string evgen::FileMuons::fFileName

private

Definition at line 62 of file FileMuons_module.cc.

std::ifstream* evgen::FileMuons::fMuonFile

private

Definition at line 67 of file FileMuons_module.cc.

TFile* evgen::FileMuons::fMuonFileR

private

Definition at line 68 of file FileMuons_module.cc.

std::string evgen::FileMuons::fMuonsFileType

private

Definition at line 63 of file FileMuons_module.cc.

std::vector<int> evgen::FileMuons::fPDG

private

Definition at line 60 of file FileMuons_module.cc.

std::string evgen::FileMuons::fTreeName

private

Definition at line 64 of file FileMuons_module.cc.

std::vector<double> evgen::FileMuons::fXYZ_Off

private

Definition at line 61 of file FileMuons_module.cc.

Float_t evgen::FileMuons::phi

private

Definition at line 77 of file FileMuons_module.cc.

Float_t evgen::FileMuons::pxtmp

private

Definition at line 73 of file FileMuons_module.cc.

Float_t evgen::FileMuons::pytmp

private

Definition at line 73 of file FileMuons_module.cc.

Float_t evgen::FileMuons::pztmp

private

Definition at line 73 of file FileMuons_module.cc.

TTree* evgen::FileMuons::TNtuple

private

Definition at line 69 of file FileMuons_module.cc.

Float_t evgen::FileMuons::xdet

private

Definition at line 78 of file FileMuons_module.cc.

Float_t evgen::FileMuons::xtmp

private

Definition at line 72 of file FileMuons_module.cc.

Float_t evgen::FileMuons::ydet

private

Definition at line 79 of file FileMuons_module.cc.

Float_t evgen::FileMuons::ytmp

private

Definition at line 72 of file FileMuons_module.cc.

Float_t evgen::FileMuons::zdet

private

Definition at line 80 of file FileMuons_module.cc.

Float_t evgen::FileMuons::ztmp

private

Definition at line 72 of file FileMuons_module.cc.

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

larsim/larsim/EventGenerator/FileMuons_module.cc

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation