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

Inheritance diagram for evgen::MUSUN:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

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

virtual	~MUSUN ()

void	produce (art::Event &evt)

void	beginJob ()

void	beginRun (art::Run &run)

void	reconfigure (fhicl::ParameterSet const &p)

void	endRun (art::Run &run)

template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const

template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const

bool	modifiesEvent () const

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)

base_engine_t &	createEngine (seed_t seed)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)

seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const

void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)

void	prepareForJob (fhicl::ParameterSet const &pset)

void	showMissingConsumes () const

Private Member Functions
void	SampleOne (unsigned int i, simb::MCTruth &mct)

void	initialization (double theta1, double theta2, double phi1, double phi2, int figflag, double s_hor, double s_ver1, double s_ver2, double &FI)

void	sampling (double &E, double &theta, double &phi, double &dep)

Private Attributes
int	fPDG
	PDG code of particles to generate. More...

double	fChargeRatio
	Charge ratio of particle / anti-particle. More...

std::string	fInputDir
	Input Directory. More...

std::string	fInputFile1
	Input File 1. More...

std::string	fInputFile2
	Input File 2. More...

std::string	fInputFile3
	Input File 3. More...

double	fCavernAngle
	Angle of the detector from the North to the East. More...

double	fRockDensity

double	fEmin
	Minimum Kinetic Energy (GeV) More...

double	fEmax
	Maximum Kinetic Energy (GeV) More...

double	fThetamin
	Minimum theta. More...

double	fThetamax
	Maximum theta. More...

double	fPhimin
	Minimum phi. More...

double	fPhimax
	Maximum phi. More...

int	figflag
	If want sampled from sphere or parallelepiped. More...

double	fXmin
	Minimum X position. More...

double	fYmin
	Minimum Y position. More...

double	fZmin
	Minimum Z position. More...

double	fXmax
	Maximum X position. More...

double	fYmax
	Maximum Y position. More...

double	fZmax
	Maximum Z position. More...

double	fT0
	Central t position (ns) in world coordinates. More...

double	fSigmaT
	Variation in t position (ns) More...

int	fTDist
	How to distribute t (gaus, or uniform) More...

int	PdgCode

double	Energy

double	phi

double	theta

double	dep

double	Time

double	Momentum

double	px0

double	py0

double	pz0

double	x0

double	y0

double	z0

double	cx

double	cy

double	cz

double	FI = 0.

double	s_hor = 0.

double	s_ver1 = 0.

double	s_ver2 = 0.

double	spmu [121][62][51]

double	fnmu [32401]

double	depth [360][91]

double	fmu [360][91]

double	the1

double	the2

double	ph1

double	ph2

double	se = 0.

double	st = 0.

double	sp = 0.

double	sd = 0.

unsigned int	NEvents = 0

TTree *	fTree

Static Private Attributes
static const int	kGAUS = 1

Detailed Description

module to produce single or multiple specified particles in the detector

Definition at line 210 of file MUSUN_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 333 of file MUSUN_module.cc.

   {
     this->reconfigure(pset);
 
     // create a default random engine; obtain the random seed from NuRandomService,
     // unless overridden in configuration with key "Seed"
     art::ServiceHandle<rndm::NuRandomService>()
       ->createEngine(*this, pset, "Seed");
 
     produces< std::vector<simb::MCTruth> >();
     produces< sumdata::RunData, art::InRun >();
   }

evgen::MUSUN::~MUSUN ( )

virtual

Definition at line 346 of file MUSUN_module.cc.

347 {

348 }

Member Function Documentation

void evgen::MUSUN::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 395 of file MUSUN_module.cc.

References art::TFileDirectory::make().

   {
     // Make the Histograms....
     art::ServiceHandle<art::TFileService> tfs;
     /*
     hPDGCode      = tfs->make<TH1D>("hPDGCode"     ,"PDG Code"      ,30 , -15           , 15             );
     hPositionX    = tfs->make<TH1D>("hPositionX"   ,"Position (cm)" ,500, ( fXmin - 10 ), ( fXmax + 10 ) );
     hPositionY    = tfs->make<TH1D>("hPositionY"   ,"Position (cm)" ,500, ( fYmin - 10 ), ( fYmax + 10 ) );
     hPositionZ    = tfs->make<TH1D>("hPositionZ"   ,"Position (cm)" ,500, ( fZmin - 10 ), ( fZmax + 10 ) );
     hTime         = tfs->make<TH1D>("hTime"        ,"Time (s)"      ,500, 0    , 1e6       );
     hMomentumHigh = tfs->make<TH1D>("hMomentumHigh","Momentum (GeV)",500, fEmin, fEmax     ); 
     hMomentum     = tfs->make<TH1D>("hMomentum"    ,"Momentum (GeV)",500, fEmin, fEmin+1e3 );
     hEnergyHigh   = tfs->make<TH1D>("hEnergyHigh"  ,"Energy (GeV)"  ,500, fEmin, fEmax     );
     hEnergy       = tfs->make<TH1D>("hEnergy"      ,"Energy (GeV)"  ,500, fEmin, fEmin+1e3 );
     hDepth        = tfs->make<TH1D>("hDepth"       ,"Depth (m)"     ,800, 0    , 14000     );
 
     hDirCosineX = tfs->make<TH1D>("hDirCosineX","Normalised Direction cosine",500, -1, 1 ); 
     hDirCosineY = tfs->make<TH1D>("hDirCosineY","Normalised Direction cosine",500, -1, 1 ); 
     hDirCosineZ = tfs->make<TH1D>("hDirCosineZ","Normalised Direction cosine",500, -1, 1 );
 
     hTheta      = tfs->make<TH1D>("hTheta"     ,"Angle (degrees)",500, 0, 90  );
     hPhi        = tfs->make<TH1D>("hPhi"       ,"Angle (degrees)",500, 0, 365 );
     */
     fTree = tfs->make<TTree>("Generator","analysis tree");
     fTree->Branch("particleID",&PdgCode, "particleID/I");
     fTree->Branch("energy" ,&Energy , "energy/D");
     fTree->Branch("time"   ,&Time   , "Time/D"  );
     fTree->Branch("posX"   ,&x0     , "posX/D"  );
     fTree->Branch("posY"   ,&y0     , "posY/D"  );
     fTree->Branch("posZ"   ,&z0     , "posZ/D"  );
     fTree->Branch("cosX"   ,&cx     , "cosX/D"  );
     fTree->Branch("cosY"   ,&cy     , "cosY/D"  );
     fTree->Branch("cosZ"   ,&cz     , "cosZ/D"  );
     fTree->Branch("theta"  ,&theta  , "theta/D" );
     fTree->Branch("phi"    ,&phi    , "phi/D"   );
     fTree->Branch("depth"  ,&dep    , "dep/D"   );
     /*
     fCryos = tfs->make<TTree>("CryoSizes","cryo tree");
     fCryos->Branch("NumTPCs" , &NumTPCs , "NumTPCs/I"       );
     fCryos->Branch("TPCMinX" , &TPCMinX , "TPCMinX[NumTPCs]/D");
     fCryos->Branch("TPCMaxX" , &TPCMaxX , "TPCMaxX[NumTPCs]/D");
     fCryos->Branch("TPCMinY" , &TPCMinY , "TPCMinY[NumTPCs]/D");
     fCryos->Branch("TPCMaxY" , &TPCMaxY , "TPCMaxY[NumTPCs]/D");
     fCryos->Branch("TPCMinZ" , &TPCMinZ , "TPCMinZ[NumTPCs]/D");
     fCryos->Branch("TPCMaxZ" , &TPCMaxZ , "TPCMaxZ[NumTPCs]/D");
     fCryos->Branch("CryoSize", &CryoSize, "CryoSize[6]/D"   );
     fCryos->Branch("DetHall" , &DetHall , "DetHall[6]/D"    );
     */
  }

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

virtual

Reimplemented from art::EDProducer.

Definition at line 448 of file MUSUN_module.cc.

References geo::GeometryCore::DetectorName(), and art::Run::put().

   {
     // grab the geometry object to see what geometry we are using
     art::ServiceHandle<geo::Geometry> geo;
     std::unique_ptr<sumdata::RunData> runcol(new sumdata::RunData(geo->DetectorName()));
     
     // Check fcl parameters were set correctly
     if ( fThetamax > 90.5 ) throw cet::exception("MUSUNGen") << "\nThetamax has to be less than " << M_PI/2 << ", but was entered as " << fThetamax << ", this causes an error so leaving program now...\n\n";  
     if ( fThetamin < 0    ) throw cet::exception("MUSUNGen") << "\nThetamin has to be more than 0, but was entered as " << fThetamin << ", this causes an error so leaving program now...\n\n";      
     if ( fThetamax < fThetamin ) throw cet::exception("MUSUNGen") << "\nMinimum angle is bigger than maximum angle....causes an error so leaving program now....\n\n"; 
     if ( fPhimax > 360.5  ) throw cet::exception("MUSUNGen") << "\nPhimax has to be less than " << 2*M_PI << ", but was entered as " << fPhimax << ", this cause an error so leaving program now...\n\n"; 
     if ( fPhimin < 0      ) throw cet::exception("MUSUNGen") << "\nPhimin has to be more than 0, but was entered as " << fPhimin << ", this causes an error so leaving program now...\n\n"; 
     if ( fPhimax < fPhimin) throw cet::exception("MUSUNGen") << "\nMinimum angle is bigger than maximum angle....causes an error so leaving program now....\n\n";
     if ( fEmax   < fEmin  ) throw cet::exception("MUSUNGen") << "\nMinimum energy is bigger than maximum energy....causes an error so leaving program now....\n\n";
      
     
     run.put(std::move(runcol));
 
     // area of the horizontal plane of the parallelepiped
     s_hor = (fZmax-fZmin)*(fXmax-fXmin);
     // area of the vertical plane of the parallelepiped, perpendicular to z-axis
     s_ver1 = (fXmax-fXmin)*(fYmax-fYmin);
     // area of the vertical plane of the parallelepiped, perpendicular to x-axis
     s_ver2 = (fZmax-fZmin)*(fYmax-fYmin);
     
     //std::cout << s_hor << " " << s_ver1 << " " << s_ver2 << std::endl;
 
     initialization(fThetamin,fThetamax,fPhimin,fPhimax,figflag,s_hor,s_ver1,s_ver2,FI );
     
     std::cout << "Material - SURF rock" << std::endl;
     std::cout << "Density = " << fRockDensity << " g/cm^3" << std::endl;
     std::cout << "Parameters for muon spectrum are from LVD best fit" << std::endl;
     std::cout << "Muon energy range = " << fEmin << " - " << fEmax << " GeV" << std::endl;
     std::cout << "Zenith angle range = " << fThetamin << " - " << fThetamax << " degrees" << std::endl;
     std::cout << "Azimuthal angle range = " << fPhimin << " - " << fPhimax << " degrees" << std::endl;
     std::cout << "Global intensity = " << FI << " (cm^2 s)^(-1) or s^(-1) (for muons on the surface)" << std::endl;
     /*
     NumTPCs = geo->NTPC(0);
     std::cout << "There are " << NumTPCs << " in cryostat 0" << std::endl;
     for (unsigned int c=0; c<geo->Ncryostats(); c++) {
       const geo::CryostatGeo& cryostat=geo->Cryostat(c);
       geo->CryostatBoundaries( CryoSize, 0 );
       std::cout << "Cryo bounds " << CryoSize[0] << " "<< CryoSize[1] << " "<< CryoSize[2] << " "<< CryoSize[3] << " "<< CryoSize[4] << " "<< CryoSize[5] << std::endl;
       for (unsigned int t=0; t<cryostat.NTPC(); t++) {
         geo::TPCID id;
         id.Cryostat=c;
         id.TPC=t;
         id.isValid=true;
         const geo::TPCGeo& tpc=cryostat.TPC(id);
         TPCMinX[t] = tpc.MinX();
         TPCMaxX[t] = tpc.MaxX();
         TPCMinY[t] = tpc.MinY();
         TPCMaxY[t] = tpc.MaxY();
         TPCMinZ[t] = tpc.MinZ();
         TPCMaxZ[t] = tpc.MaxZ();
         std::cout << t << "\t" << TPCMinX[t] << " " << TPCMaxX[t] << " " << TPCMinY[t] << " " << TPCMaxY[t] << " " << TPCMinZ[t] << " " << TPCMaxZ[t] << std::endl;
       }
     }
     fCryos -> Fill();
     */
     return;
   }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const & it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const & it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

 {
   return rng()->createEngine(placeholder_schedule_id(), seed);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make
	)

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make,
		label_t const &	engine_label
	)

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
 }

CurrentProcessingContext const * art::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

   {
     return current_context_.get();
   }

void evgen::MUSUN::endRun ( art::Run & run )

virtual

Reimplemented from art::EDProducer.

Definition at line 514 of file MUSUN_module.cc.

   {
     std::cout << "\n\nNumber of muons = " << NEvents << std::endl;
     std::cout << "Mean muon energy = " << se/NEvents << " GeV" << std::endl;
     std::cout << "Mean zenith angle (theta) = " << st/NEvents << " degrees" << std::endl;
     std::cout << "Mean azimuthal angle (phi)= " << sp/NEvents << " degrees" << std::endl;
     std::cout << "Mean slant depth = " << sd/NEvents << " m w.e." << std::endl;
   }

EngineCreator::seed_t EngineCreator::get_seed_value	(	fhicl::ParameterSet const &	pset,
		char const	key[] = `"seed"`,
		seed_t const	implicit_seed = `-1`
	)

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

 {
   auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
   return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
 }

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const & instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

   {
     return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                                instanceName);
   }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID	(	ModuleDescription const &	moduleDescription,
		std::string const &	instanceName
	)		const

inherited

Definition at line 56 of file ProducerBase.h.

References B, and art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

   {
     auto const& pd =
       get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
     return pd.productID();
   }

void evgen::MUSUN::initialization	(	double	theta1,
		double	theta2,
		double	phi1,
		double	phi2,
		int	figflag,
		double	s_hor,
		double	s_ver1,
		double	s_ver2,
		double &	FI
	)

private

Definition at line 695 of file MUSUN_module.cc.

References e, and in.

   {
     //
     //  Read in the data files
     //
     int lineNumber = 0, index = 0;
     char inputLine[10000];
     std::string fROOTfile; 
 
     for (int a=0;a<121;++a) for (int b=0;b<62;++b) for (int c=0;c<50;++c) spmu[a][b][c]=0;
     for (int a=0;a<23401;++a)                      fnmu[a]    = 0;
     for (int a=0;a<360;++a) for (int b=0;b<91;++b) depth[a][b]= 0;
     for (int a=0;a<360;++a) for (int b=0;b<91;++b) fmu[a][b]  = 0;
 
     std::ostringstream File1LocStream;
     File1LocStream  << fInputDir << fInputFile1;
     std::string File1Loc = File1LocStream. str();
     cet::search_path sp1("FW_SEARCH_PATH");
     if( sp1.find_file(fInputFile1, fROOTfile) ) File1Loc = fROOTfile;
     std::ifstream file1( File1Loc.c_str(), std::ios::in );
     if (!file1.good() ) throw cet::exception("MUSUNGen") << "\nFile1 " << fInputFile1 << " not found in FW_SEARCH_PATH or at " <<  fInputDir <<"\n\n";
     
     while( file1.good() ) {
       //std::cout << "Looking at file 1...." << std::endl;
       file1.getline( inputLine, 9999 );
       char *token;
       token = strtok( inputLine, " " );
       while( token != NULL ) {
         //std::cout << "While loop file 1..." << std::endl;
         fmu[index][lineNumber] = atof( token );
         token = strtok( NULL, " " );
         index++;
         if( index == 360 ) {
           //std::cout << "If statement file 1..." << std::endl;
           index = 0;
           lineNumber++;
         }
       }
     }
     file1.close();
     
     std::ostringstream File2LocStream;
     File2LocStream  << fInputDir << fInputFile2;
     std::string File2Loc = File2LocStream. str();
     cet::search_path sp2("FW_SEARCH_PATH");
     if( sp2.find_file(fInputFile2, fROOTfile) ) File2Loc = fROOTfile;
     std::ifstream file2( File2Loc.c_str(), std::ios::binary|std::ios::in );
     if (!file2.good() ) throw cet::exception("MUSUNGen") << "\nFile2 " << fInputFile2 << " not found in FW_SEARCH_PATH or at " <<  fInputDir <<"\n\n";
   
     int i1 = 0, i2 = 0, i3 = 0;
     float readVal;
     while( file2.good() ) {
       //std::cout << "Looking at file 2...." << std::endl;
       file2.read((char *)(&readVal), sizeof(float));
       spmu[i1][i2][i3] = readVal;
       i1++;
       if( i1 == 121 ) {
         //std::cout << "First if statement file 2..." << std::endl;
         i2++;
         i1 = 0;
       }
       if( i2 == 62 ) {
         //std::cout << "Second if statement file 2..." << std::endl;
         i3++;
         i2 = 0;
       }
     }
     file2.close();
     for( int i=0; i<120; i++ )
       for( int j=0; j<62; j++ )
         for( int k=0; k<51; k++ )
           spmu[i][j][k] = spmu[i+1][j][k];
     spmu[1][1][0] = 0.000853544;
     //std::cout << "Set spmu to some value..." << std::endl;
     
     std::ostringstream File3LocStream;
     File3LocStream  << fInputDir << fInputFile3;
     std::string File3Loc = File3LocStream. str();
     cet::search_path sp3("FW_SEARCH_PATH");
     if( sp3.find_file(fInputFile3, fROOTfile) ) File3Loc = fROOTfile;
     std::ifstream file3( File3Loc.c_str(), std::ios::in );
     if (!file3.good() ) throw cet::exception("MUSUNGen") << "\nFile3 " << fInputFile3 << " not found in FW_SEARCH_PATH or at " <<  fInputDir <<"\n\n";
 
     lineNumber = index = 0;
     while( file3.good() ) {
       //std::cout << "Looking at file 3...." << std::endl;
       file3.getline( inputLine, 9999 );
       char *token;
       token = strtok( inputLine, " " );
       while( token != NULL ) {
         //std::cout << "While loop file 3..." << std::endl;
         depth[index][lineNumber] = atof( token );
         token = strtok( NULL, " " );
         index++;
         if( index == 360 ) {
           //std::cout << "If statement file 3..." << std::endl;
           index = 0;
           lineNumber++;
         }
       }
     }
     file3.close();
     
     //
     //  Set up variables
     //
     
     the1 = theta1;
     the2 = theta2;
     // for c2: c1 and c2 are unused
     //double c1 = cos(M_PI/180.*theta1);
     //double c2 = cos(M_PI/180.*theta2);
     ph1 = M_PI/180.*phi1;
     ph2 = M_PI/180.*phi2;
     // for c2: dph is unused
     //double dph = ph2-ph1;
     
     int ipc = 1;
     double theta = theta1;
     double dc = 1.;
     double sc = 0.;
     int iteration = 0;
     while( theta < theta2-dc/2. ) {
         theta += dc/2.;
         double theta0 = M_PI/180. * theta;
         double cc = cos(theta0);
         double ash = s_hor * cc;
         double asv01 = s_ver1 * sqrt(1. - cc*cc);
         double asv02 = s_ver2 * sqrt(1. - cc*cc);
         int ic1 = (theta + 0.999);
         int ic2 = ic1 + 1;
         if( ic2 > 91 ) ic2 = 91;
         if( ic1 < 1 ) ic1 = 1;
         double phi = phi1;
         double dp = 1.;
         
         while( phi < phi2-dp/2. ) {
             phi += dp/2.;
             //  the long side of the cavern is pointing at 14 deg to the north:
             //  double phi0 = M_PI / 180. * (phi + 90. - 14);
             
             //  Want our co-ord system going from East to South.
             double phi0 = M_PI / 180. * (phi + fCavernAngle);
             
             double asv1 = asv01 * fabs(cos(phi0));
             double asv2 = asv02 * fabs(sin(phi0));
             double asv0 = ash + asv1 + asv2;
             double fl = 1.;
             if( figflag == 1 )
                 fl = asv0;
             int ip1 = (phi + 0.999);
             int ip2 = ip1 + 1;
             if( ip2 > 360 ) ip2 = 1;
             if( ip1 < 1 ) ip1 = 360;
             double sp1 = 0.;
             
             for( int ii=0; ii<4; ii++ ) {
                 int iic = ii/2;
                 int iip = ii%2;
                 if(ip1==360 && (ii==1 || ii==3) ) iip = -359;
                 if( fmu[ip1+iip-1][ic1+iic-1] < 0 ) {
                   if ( pow(10.,fmu[ip1+iip-1][ic1+iic-1]) / 4 > 1e-6 ) {
                     std::cout << "Looking at fmu [ " << ip1 << " + " << iip << " - 1 (" << ip1+iip-1 << ") ] [ " << ic1 << " + " << iic << " - 1 ("<< ic1+iic-1 << ") ] ."
                               << "\nChanging sp1 from " << sp1 << " to " <<  sp1 + pow(10.,fmu[ip1+iip-1][ic1+iic-1]) / 4 << "..........." << sp1 << " + 10 ^ (" << fmu[ip1+iip-1][ic1+iic-1] << ") / 4 "
                               << std::endl;
                   }
                   sp1 = sp1 + pow(10.,fmu[ip1+iip-1][ic1+iic-1]) / 4;
                 }
             }
             /*
               std::cout << iteration<< " time of new sc value! Theta " << theta << ", phi " << phi + dp / 2. << ", sc = " << sc + sp1 * fl * dp * M_PI / 180. * sin(theta0) * dc * M_PI / 180. << " = " 
               << sc << " + " << sp1 << " * " << fl << " * " << dp << " * " << M_PI/180 << " * sin(" << theta0 << ") * " << dc << " * " << M_PI/180 << ".....sin(theta)=" << sin(theta) << "\n" 
               << std::endl; */
             sc = sc + sp1 * fl * dp * M_PI / 180. * sin(theta0) * dc * M_PI / 180.;
             ++iteration;
             ipc = ipc + 1;
             fnmu[ipc-1] = sc;
             phi = phi + dp / 2.;
         }
         
         theta = theta + dc / 2.;
     }
     //std::cout << *FI << " = " << sc << std::endl;
     FI = sc;
     for( int ipc1 = 0; ipc1 < ipc; ipc1++ )
       fnmu[ipc1] = fnmu[ipc1] / fnmu[ipc-1];
   }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::mayConsume ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const & it )

inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const & it )

inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

     {
       return true;
     }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

 {
   static Consumer invalid{InvalidTag{}};
   return &invalid;
 }

void art::Consumer::prepareForJob ( fhicl::ParameterSet const & pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

 {
   if (!moduleContext_)
     return;
 
   pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
   for (auto& consumablesPerBranch : consumables_) {
     cet::sort_all(consumablesPerBranch);
   }
 }

void evgen::MUSUN::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 525 of file MUSUN_module.cc.

References simb::kSingleParticle, LOG_DEBUG, art::Event::put(), and simb::MCTruth::SetOrigin().

   {
     std::unique_ptr< std::vector<simb::MCTruth> > truthcol(new std::vector<simb::MCTruth>);
 
     simb::MCTruth truth;
     truth.SetOrigin(simb::kSingleParticle);
     
     ++NEvents;
     SampleOne(NEvents,truth);
 
     LOG_DEBUG("MUSUN") << truth;
 
     truthcol->push_back(truth);
 
     evt.put(std::move(truthcol));
 
     return;
   }

void evgen::MUSUN::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 353 of file MUSUN_module.cc.

References fhicl::ParameterSet::get().

   {
     // do not put seed in reconfigure because we don't want to reset 
     // the seed midstream
 
     fPDG           = p.get< int                 >("PDG");
     fChargeRatio   = p.get< double              >("ChargeRatio");
 
     fInputDir      = p.get< std::string         >("InputDir");
     fInputFile1    = p.get< std::string         >("InputFile1");
     fInputFile2    = p.get< std::string         >("InputFile2");
     fInputFile3    = p.get< std::string         >("InputFile3");
 
     fCavernAngle   = p.get< double              >("CavernAngle");
     fRockDensity   = p.get< double              >("RockDensity");
 
     fEmin          = p.get< double              >("Emin");
     fEmax          = p.get< double              >("Emax");
 
     fThetamin      = p.get< double              >("Thetamin");
     fThetamax      = p.get< double              >("Thetamax");
 
     fPhimin        = p.get< double              >("Phimin");
     fPhimax        = p.get< double              >("Phimax");
 
     figflag        = p.get<int                  >("igflag");
     fXmin          = p.get<double               >("Xmin");
     fYmin          = p.get<double               >("Ymin");
     fZmin          = p.get<double               >("Zmin");
     fXmax          = p.get<double               >("Xmax");
     fYmax          = p.get<double               >("Ymax");
     fZmax          = p.get<double               >("Zmax");
 
     fT0            = p.get< double              >("T0");
     fSigmaT        = p.get< double              >("SigmaT");
     fTDist         = p.get< int                 >("TDist");
 
     return;
   }

void evgen::MUSUN::SampleOne	(	unsigned int	i,
		simb::MCTruth &	mct
	)

private

Definition at line 549 of file MUSUN_module.cc.

References simb::MCTruth::Add(), simb::MCParticle::AddTrajectoryPoint(), art::RandomNumberGenerator::getEngine(), part, and ss.

                                                        {
 
     // get the random number generator service and make some CLHEP generators
     art::ServiceHandle<art::RandomNumberGenerator> rng;
     CLHEP::HepRandomEngine &engine = rng->getEngine();
     CLHEP::RandFlat   flat(engine);
     CLHEP::RandGaussQ gauss(engine);
 
     Energy = 0;
     theta  = 0;
     phi    = 0; 
     dep    = 0;
     Time   = 0;
     
     sampling( Energy, theta, phi, dep );
 
     theta = theta* M_PI/180;
 
     //  changing the angle phi so z-axis is positioned along the long side
     //  of the cavern pointing at 14 deg from the North to the East.
     //  phi += (90. - 14.0);
     //  Want our co-ord rotation going from East to South.
     phi += fCavernAngle;
     if( phi >= 360. )
       phi -= 360.;
     if( phi < 0 )
       phi += 360.;
     phi *= M_PI / 180.;
     
     // set track id to -i as these are all primary particles and have id <= 0
     int trackid = -1*(i+1);
     std::string primary("primary");
     
     // Work out whether particle/antiparticle, and mass...
     double m = 0.0;
     PdgCode = fPDG;
     double ChargeCheck = 1./ ( 1 + fChargeRatio );
     double pdgfire = flat.fire();
     if ( pdgfire < ChargeCheck ) PdgCode=-PdgCode;
 
     static TDatabasePDG  pdgt;
     TParticlePDG* pdgp = pdgt.GetParticle(PdgCode);
     if (pdgp) m = pdgp->Mass();
 
     //std::cout << pdgfire << " " << ChargeCheck << " " << PdgCode << " " << m << std::endl;
     
     // Work out T0...
     if(fTDist==kGAUS){
       Time = gauss.fire(fT0, fSigmaT);
     }
     else {
       Time = fT0 + fSigmaT*(2.0*flat.fire()-1.0);
     }
         
     //  The minus sign above is for y-axis pointing up, so the y-momentum
     //  is always pointing down
     cx       = -sin(theta)*sin(phi);
     cy       = -cos(theta);
     cz       = +sin(theta)*cos(phi);
     Momentum = std::sqrt(Energy*Energy-m*m); // Get momentum
     px0      = Momentum * cx;
     py0      = Momentum * cy;
     pz0      = Momentum * cz;
     TLorentzVector pvec(px0, py0, pz0, Energy );   
     
     //  Muon coordinates
     double sh1 = s_hor * cos(theta);
     double sv1 = s_ver1 * sin(theta) * fabs(cos(phi));
     double sv2 = s_ver2 * sin(theta) * fabs(sin(phi));
     double ss = sh1 + sv1 + sv2;
     double xfl1 = flat.fire();
     if( xfl1 <= sh1/ss ) {
       x0 = (fXmax - fXmin)*flat.fire() + fXmin;
       y0 = fYmax;
       z0 = (fZmax - fZmin)*flat.fire() + fZmin;
     } else if( xfl1 <= (sh1+sv1)/ss ) {
       x0 = (fXmax - fXmin)*flat.fire() + fXmin;
       y0 = (fYmax - fYmin)*flat.fire() + fYmin;
       if( cz >= 0 ) z0 = fZmax;
       else z0 = fZmin;
     } else {
       if( cx >= 0 ) x0 = fXmin;
       else x0 = fXmax;
       y0 = (fYmax - fYmin)*flat.fire() + fYmin;
       z0 = (fZmax - fZmin)*flat.fire() + fZmin;
     }
     // Make Lorentz vector for x and t....
     TLorentzVector pos(x0, y0, z0, Time); 
    
     //  Parameters written to the file muons_surf_v2_test*.dat
     //      nmu - muon sequential number
     //      id_part - muon charge (10 - positive, 11 - negative )
     //      Energy  - total muon energy in GeV assuming ultrarelativistic muons
     //      x0, y0, z0 - muon coordinates on the surface of parallelepiped
     //          specified above; x-axis and y-axis are pointing in the
     //          directions such that the angle phi (from the slant depth
     //          distribution files) is measured from x to y. z-axis is
     //          pointing upwards.
     //      cx, cy, cz - direction cosines.
   
     simb::MCParticle part(trackid, PdgCode, primary);
     part.AddTrajectoryPoint(pos, pvec);
 
     mct.Add(part);
 
     theta = theta * 180/M_PI;
     phi   = phi   * 180/M_PI;
 
     // Sum energies, angles, depth for average outputs. 
     se += Energy;
     st += theta;
     sp += phi;
     sd += dep;
     
     // Fill Histograms.....
     /*
     hPDGCode      ->Fill (PdgCode);
     hPositionX    ->Fill (x0);
     hPositionY    ->Fill (y0);
     hPositionZ    ->Fill (z0);
     hTime         ->Fill (Time);
     hMomentumHigh ->Fill (Momentum); 
     hMomentum     ->Fill (Momentum); 
     hEnergyHigh   ->Fill (Energy);
     hEnergy       ->Fill (Energy);
     hDepth        ->Fill (dep);
     hDirCosineX   ->Fill (cx); 
     hDirCosineY   ->Fill (cy); 
     hDirCosineZ   ->Fill (cz);
     hTheta        ->Fill (theta);
     hPhi          ->Fill (phi);
     */
     /*
     // Write event by event outsputs.....
     std::cout << "Theta: " << theta << " Phi: " << phi << " Energy: " << Energy << " Momentum: " << Momentum << std::endl; 
     std::cout << "x: " <<  pos.X() << " y: " << pos.Y() << " z: " << pos.Z() << " time: " << pos.T() << std::endl;   
     std::cout << "Px: " <<  pvec.Px() << " Py: " << pvec.Py() << " Pz: " << pvec.Pz() << std::endl;
     std::cout << "Normalised..." << cx << " " << cy << " " << cz << std::endl;
     */
     fTree->Fill();
 
   }

void evgen::MUSUN::sampling	(	double &	E,
		double &	theta,
		double &	phi,
		double &	dep
	)

private

Definition at line 888 of file MUSUN_module.cc.

References DEFINE_ART_MODULE, and art::RandomNumberGenerator::getEngine().

   {
     // get the random number generator service and make some CLHEP generators
     art::ServiceHandle<art::RandomNumberGenerator> rng;
     CLHEP::HepRandomEngine &engine = rng->getEngine();
     CLHEP::RandFlat   flat(engine);
     CLHEP::RandGaussQ gauss(engine);
 
     #if 0 // this code is disabled for good
     double xfl = flat.fire();
     int loIndex = 0, hiIndex = 32400;
     int i = (loIndex+hiIndex)/2;
     bool foundIndex = false;
     if( xfl < fnmu[loIndex] ) {
       i = loIndex;
       foundIndex = true;
     } else if ( xfl > fnmu[hiIndex] ) {
       i = hiIndex;
       foundIndex = true;
     } else if ( xfl > fnmu[i-1] && xfl <= fnmu[i] )
       foundIndex = true;
     while( !foundIndex ) {
       if( xfl < fnmu[i] )
         hiIndex = i;
       else
         loIndex = i;
       i = (loIndex + hiIndex)/2;
       
       if( xfl > fnmu[i-1] && xfl <= fnmu[i] )
         foundIndex = true;
     }
     #else
     double xfl = flat.fire();
     int i = 0;
     while ( xfl > fnmu[i] ) ++i;
     #endif
     int ic = (i-2)/360;
     int ip = i-2-ic*360;
         
     xfl = flat.fire();
     theta = the1 + ((double)ic+xfl);
     xfl = flat.fire();
     phi = ph1 + ((double)ip+xfl);
     if ( phi > 360 ) phi = phi -360;
     dep = depth[ip][ic] * fRockDensity;
     
     int ic1 = cos(M_PI/180.*theta) * 50.;
     if( ic1 < 0 )
       ic1 = 0;
     if( ic1 > 50 )
       ic1 = 50;
     int ip1 = dep / 200. - 16;
     if( ip1 < 0 )
       ip1 = 0;
     if( ip1 > 61 )
       ip1 = 61;
 
     xfl = flat.fire();
     #if 0
     loIndex = 0, hiIndex = 120;
     int j = (loIndex+hiIndex)/2;
     foundIndex = false;
     if( xfl < spmu[loIndex][ip1][ic1] ) {
       j = loIndex;
       foundIndex = true;
     } else if ( xfl > spmu[hiIndex][ip1][ic1] ) {
       j = hiIndex;
       foundIndex = true;
     } else if ( xfl > spmu[j-1][ip1][ic1] && xfl <= spmu[j][ip1][ic1] )
       foundIndex = true;
     while( !foundIndex ) {
       if( xfl < spmu[j][ip1][ic1] )
         hiIndex = j;
       else
         loIndex = j;
       j = (loIndex + hiIndex)/2;
       
       if( xfl > spmu[j-1][ip1][ic1] && xfl <= spmu[j][ip1][ic1] )
         foundIndex = true;
     }
     #else
     int j = 0;
     while ( xfl > spmu[j][ip1][ic1] ) ++j;
     #endif
 
     double En1 = 0.05 * (j-1);
     double En2 = 0.05 * (j);
     E = pow(10.,En1 + (En2 - En1)*flat.fire());
     
     return;
   }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

 {
   if (!moduleContext_)
     return;
 
   // If none of the branches have missing consumes statements, exit early.
   if (std::all_of(cbegin(missingConsumes_),
                   cend(missingConsumes_),
                   [](auto const& perBranch) { return perBranch.empty(); }))
     return;
 
   constexpr cet::HorizontalRule rule{60};
   mf::LogPrint log{"MTdiagnostics"};
   log << '\n'
       << rule('=') << '\n'
       << "The following consumes (or mayConsume) statements are missing from\n"
       << module_context(moduleDescription_) << '\n'
       << rule('-') << '\n';
 
   cet::for_all_with_index(
     missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
       for (auto const& pi : perBranch) {
         log << "  "
             << assemble_consumes_statement(static_cast<BranchType>(i), pi)
             << '\n';
       }
     });
   log << rule('=');
 }

void art::Consumer::validateConsumedProduct	(	BranchType const	bt,
		ProductInfo const &	pi
	)

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

 {
   // Early exits if consumes tracking has been disabled or if the
   // consumed product is an allowed consumable.
   if (!moduleContext_)
     return;
 
   if (cet::binary_search_all(consumables_[bt], pi))
     return;
 
   if (requireConsumes_) {
     throw Exception(errors::ProductRegistrationFailure,
                     "Consumer: an error occurred during validation of a "
                     "retrieved product\n\n")
       << "The following consumes (or mayConsume) statement is missing from\n"
       << module_context(moduleDescription_) << ":\n\n"
       << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
   }
 
   missingConsumes_[bt].insert(pi);
 }

Member Data Documentation

double evgen::MUSUN::cx

private

Definition at line 290 of file MUSUN_module.cc.

double evgen::MUSUN::cy

private

Definition at line 290 of file MUSUN_module.cc.

double evgen::MUSUN::cz

private

Definition at line 290 of file MUSUN_module.cc.

double evgen::MUSUN::dep

private

Definition at line 288 of file MUSUN_module.cc.

double evgen::MUSUN::depth[360][91]

private

Definition at line 300 of file MUSUN_module.cc.

double evgen::MUSUN::Energy

private

Definition at line 288 of file MUSUN_module.cc.

double evgen::MUSUN::fCavernAngle

private

Angle of the detector from the North to the East.

Definition at line 243 of file MUSUN_module.cc.

double evgen::MUSUN::fChargeRatio

private

Charge ratio of particle / anti-particle.

Definition at line 236 of file MUSUN_module.cc.

double evgen::MUSUN::fEmax

private

Maximum Kinetic Energy (GeV)

Definition at line 250 of file MUSUN_module.cc.

double evgen::MUSUN::fEmin

private

Minimum Kinetic Energy (GeV)

Definition at line 249 of file MUSUN_module.cc.

double evgen::MUSUN::FI = 0.

private

Definition at line 293 of file MUSUN_module.cc.

int evgen::MUSUN::figflag

private

If want sampled from sphere or parallelepiped.

Definition at line 257 of file MUSUN_module.cc.

std::string evgen::MUSUN::fInputDir

private

Input Directory.

Definition at line 238 of file MUSUN_module.cc.

std::string evgen::MUSUN::fInputFile1

private

Input File 1.

Definition at line 239 of file MUSUN_module.cc.

std::string evgen::MUSUN::fInputFile2

private

Input File 2.

Definition at line 240 of file MUSUN_module.cc.

std::string evgen::MUSUN::fInputFile3

private

Input File 3.

Definition at line 241 of file MUSUN_module.cc.

double evgen::MUSUN::fmu[360][91]

private

Definition at line 301 of file MUSUN_module.cc.

double evgen::MUSUN::fnmu[32401]

private

Definition at line 299 of file MUSUN_module.cc.

int evgen::MUSUN::fPDG

private

PDG code of particles to generate.

Definition at line 235 of file MUSUN_module.cc.

double evgen::MUSUN::fPhimax

private

Maximum phi.

Definition at line 255 of file MUSUN_module.cc.

double evgen::MUSUN::fPhimin

private

Minimum phi.

Definition at line 254 of file MUSUN_module.cc.

double evgen::MUSUN::fRockDensity

private

Default rock density is 2.70 g cm-3. If this is changed then the three input files need to be remade. If there is a desire for this contact Vitaly Kudryavtsev at V.Kud.nosp@m.ryav.nosp@m.tsev@.nosp@m.shef.nosp@m..ac.u.nosp@m.k

Definition at line 244 of file MUSUN_module.cc.

double evgen::MUSUN::fSigmaT

private

Variation in t position (ns)

Definition at line 266 of file MUSUN_module.cc.

double evgen::MUSUN::fT0

private

Central t position (ns) in world coordinates.

Definition at line 265 of file MUSUN_module.cc.

int evgen::MUSUN::fTDist

private

How to distribute t (gaus, or uniform)

Definition at line 267 of file MUSUN_module.cc.

double evgen::MUSUN::fThetamax

private

Maximum theta.

Definition at line 253 of file MUSUN_module.cc.

double evgen::MUSUN::fThetamin

private

Minimum theta.

Definition at line 252 of file MUSUN_module.cc.

TTree* evgen::MUSUN::fTree

private

Definition at line 313 of file MUSUN_module.cc.

double evgen::MUSUN::fXmax

private

Maximum X position.

Definition at line 261 of file MUSUN_module.cc.

double evgen::MUSUN::fXmin

private

Minimum X position.

Definition at line 258 of file MUSUN_module.cc.

double evgen::MUSUN::fYmax

private

Maximum Y position.

Definition at line 262 of file MUSUN_module.cc.

double evgen::MUSUN::fYmin

private

Minimum Y position.

Definition at line 259 of file MUSUN_module.cc.

double evgen::MUSUN::fZmax

private

Maximum Z position.

Definition at line 263 of file MUSUN_module.cc.

double evgen::MUSUN::fZmin

private

Minimum Z position.

Definition at line 260 of file MUSUN_module.cc.

const int evgen::MUSUN::kGAUS = 1

staticprivate

Definition at line 233 of file MUSUN_module.cc.

double evgen::MUSUN::Momentum

private

Definition at line 289 of file MUSUN_module.cc.

unsigned int evgen::MUSUN::NEvents = 0

private

Definition at line 310 of file MUSUN_module.cc.

int evgen::MUSUN::PdgCode

private

Definition at line 287 of file MUSUN_module.cc.

double evgen::MUSUN::ph1

private

Definition at line 304 of file MUSUN_module.cc.

double evgen::MUSUN::ph2

private

Definition at line 304 of file MUSUN_module.cc.

double evgen::MUSUN::phi

private

Definition at line 288 of file MUSUN_module.cc.

double evgen::MUSUN::px0

private

Definition at line 289 of file MUSUN_module.cc.

double evgen::MUSUN::py0

private

Definition at line 289 of file MUSUN_module.cc.

double evgen::MUSUN::pz0

private

Definition at line 289 of file MUSUN_module.cc.

double evgen::MUSUN::s_hor = 0.

private

Definition at line 294 of file MUSUN_module.cc.

double evgen::MUSUN::s_ver1 = 0.

private

Definition at line 295 of file MUSUN_module.cc.

double evgen::MUSUN::s_ver2 = 0.

private

Definition at line 296 of file MUSUN_module.cc.

double evgen::MUSUN::sd = 0.

private

Definition at line 308 of file MUSUN_module.cc.

double evgen::MUSUN::se = 0.

private

Definition at line 305 of file MUSUN_module.cc.

double evgen::MUSUN::sp = 0.

private

Definition at line 307 of file MUSUN_module.cc.

double evgen::MUSUN::spmu[121][62][51]

private

Definition at line 298 of file MUSUN_module.cc.

double evgen::MUSUN::st = 0.

private

Definition at line 306 of file MUSUN_module.cc.

double evgen::MUSUN::the1

private

Definition at line 304 of file MUSUN_module.cc.

double evgen::MUSUN::the2

private

Definition at line 304 of file MUSUN_module.cc.

double evgen::MUSUN::theta

private

Definition at line 288 of file MUSUN_module.cc.

double evgen::MUSUN::Time

private

Definition at line 288 of file MUSUN_module.cc.

double evgen::MUSUN::x0

private

Definition at line 290 of file MUSUN_module.cc.

double evgen::MUSUN::y0

private

Definition at line 290 of file MUSUN_module.cc.

double evgen::MUSUN::z0

private

Definition at line 290 of file MUSUN_module.cc.

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

larsim/v06_53_00/source/larsim/EventGenerator/MuonPropagation/MUSUN_module.cc

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation