gar Namespace Reference

General GArSoft Utilities. More...

Namespaces
	cheat
	debug
	details
	detinfo
	dump
	Namespace for GArSoft dumping utilities.
	evd
	evd3d
	evgen
	filt
	gar_pandora
	garg4
	geo
	raw
	rec
	rosim
	sdp
	sumdata
	util

Classes
class	anatest
class	anatree
class	CAFHelper
class	CAFMaker
class	HNLAnalysis
class	MatchingPerformance
class	MomentumPerformance
class	NestAlg
class	ParamSim
class	ProviderPack
	Container for a list of pointers to providers. More...
class	ServiceProviderImplementationWrapper
	Service implementation returning a provider. More...
class	ServiceProviderInterfaceWrapper
	Service returning a provider interface. More...
class	SimpleServiceProviderWrapper
	Service returning a provider. More...
class	StructuredTree
struct	UncopiableAndUnmovableClass
	An empty class that can't be copied nor moved. More...

Typedefs
typedef std::pair< float, std::string >	P
template<typename... Services>
using	ProviderPackFromServices = gar::ProviderPack< typename Services::provider_type... >
	Type of provider pack with providers from all specified Services. More...

Functions
garana::GTruth	MakeAnaGTruth (const simb::GTruth &gt, const int &vtxregion)
garana::FSParticle	MakeFSParticle (const simb::MCParticle &mcp)
garana::G4Particle	MakeG4Particle (const simb::MCParticle &mcp, int parentPdg, int progenitorPdg, int progenitorTrackId, const vector< pair< TLorentzVector, TLorentzVector >> &positions, const vector< pair< TLorentzVector, TLorentzVector >> &momenta, const vector< int > &regions, const vector< size_t > &nptsPerRegion)
garana::CaloCluster	MakeAnaCalCluster (const rec::Cluster &clust, const int &region, const std::vector< std::pair< int, float >> &edeps)
garana::Vee	MakeAnaVee (const rec::Vee &vee)
garana::Vertex	MakeAnaVtx (const rec::Vertex &vtx)
int	ProcessNameToCode (std::string const &p)
void	FillGTruth (const simb::GTruth &gt, garana::GTruth &outtruth)
const garana::Track	MakeAnaTrack (const rec::Track &trk, const vector< pair< int, float >> &pidf, const vector< pair< int, float >> &pidb, float ionf, float ionb, const vector< pair< UInt_t, TLorentzVector >> &posBeg, const vector< pair< UInt_t, TLorentzVector >> &posEnd, const vector< pair< UInt_t, TLorentzVector >> &momBeg, const vector< pair< UInt_t, TLorentzVector >> &momEnd, const vector< pair< int, float >> &edeps)
garana::CaloCluster	MakeAnaCalCluster (const rec::Cluster &clust, const int &region, const vector< pair< int, float >> &edeps)
template<typename... Providers>
ProviderPack< Providers... >	makeProviderPack (Providers const *...providers)
	Function to create a ParameterPack from the function arguments. More...
std::set< std::string > const &	IgnorableProviderConfigKeys ()
	Returns a list of configuration keys that providers should ignore. More...
template<class T >
T::provider_type const *	providerFrom ()
	Returns a constant pointer to the provider of specified service. More...
template<typename... Services>
ProviderPackFromServices< Services... >	extractProviders ()
	Returns a provider pack with providers from specified services. More...
G4double	UnivScreenFunc (G4double E, G4double Z, G4double A)

Variables
std::unordered_map< std::string, int >	processMap
G4bool	diffusion = true
G4bool	SinglePhase =false
G4bool	ThomasImelTail =true
G4bool	OutElectrons =true
G4double	biExc = 0.77

Detailed Description

General GArSoft Utilities.

Monte Carlo sdpulation.

Raw data description and utilities.

PseudoLayerPlugin class.

Namespace collecting geometry-related classes utilities.

Geant4 interface.

Event Generation using GENIE, cosmics or single particles.

Event display classes.

GArSoft includes.

Version

Id

Style.cxx,v 1.1.1.1 2010/11/10 19:44:54 p-novaart Exp

Author

messi.nosp@m.er@i.nosp@m.ndian.nosp@m.a.ed.nosp@m.u

bjpjo.nosp@m.nes@.nosp@m.mit.e.nosp@m.du

File with various classes for a generic nearest neighbour type clustering.

Typedef Documentation

typedef pair< float, string > gar::P

Definition at line 75 of file anatree_module.cc.

template<typename... Services>

using gar::ProviderPackFromServices = typedef gar::ProviderPack<typename Services::provider_type...>

Type of provider pack with providers from all specified Services.

Definition at line 26 of file ServicePack.h.

Function Documentation

template<typename... Services>

ProviderPackFromServices<Services...> gar::extractProviders

(

)

Returns a provider pack with providers from specified services.

Template Parameters

Services

the services to extract the providers from

Returns

a ProviderPack containing the current service providers

This convenience function automatically extracts all the service providers from a list of services. This is convenient if an algorithm or service accepts a provider pack for setup:

algo->Setup(extractProviders<
  detinfo::DetectorPropertiesService, detinfo::GArPropertiesService
  >());

Also note that the provider packs can rearrange their elements, so the call above should work just the same as:

algo->Setup(extractProviders<
  detinfo::GArPropertiesService, detinfo::DetectorPropertiesService
  >());

If a provider is needed in the setup argument that is not provided by any of the specified services, a compilation error will occur.

Definition at line 53 of file ServicePack.h.

    { return { gar::providerFrom<Services>()... }; }

void gar::FillGTruth	(	const simb::GTruth &	gt,
		garana::GTruth &	outtruth
	)

Definition at line 151 of file AnaUtils.cxx.

                                                                  {
  
        outtruth.fVertex = gt.fVertex;
        outtruth.fVertex.SetVect(100.*gt.fVertex.Vect()); // GENIE uses meters - convert to cm (same as MCParticle)
        outtruth.fVertex.SetT(1.e9*gt.fVertex.T()); // GENIE uses seconds - convert to ns (same as MCParticle)
        outtruth.fweight = gt.fweight;
        outtruth.fprobability = gt.fprobability;
        outtruth.fXsec = gt.fXsec;
        outtruth.fDiffXsec = gt.fDiffXsec;
        outtruth.fGPhaseSpace = gt.fGPhaseSpace;

        outtruth.fProbePDG = gt.fProbePDG;
        outtruth.fProbeP4 = gt.fProbeP4;
        outtruth.fTgtP4 = gt.fTgtP4;

        outtruth.ftgtZ = gt.ftgtZ;
        outtruth.ftgtA = gt.ftgtA;
        outtruth.ftgtPDG = gt.ftgtPDG;
        outtruth.fHitNucPDG = gt.fHitNucPDG;
        outtruth.fHitQrkPDG = gt.fHitQrkPDG;
        outtruth.fIsSeaQuark = gt.fIsSeaQuark;
        outtruth.fHitNucP4 = gt.fHitNucP4;
        outtruth.fHitNucPos = gt.fHitNucPos;

        outtruth.fGscatter = gt.fGscatter;
        outtruth.fGint = gt.fGint;

        outtruth.fgQ2 = gt.fgQ2;
        outtruth.fgq2 = gt.fgq2;
        outtruth.fgW = gt.fgW ;
        outtruth.fgT = gt.fgT;
        outtruth.fgX = gt.fgX;
        outtruth.fgY = gt.fgY;
        outtruth.fFSleptonP4 = gt.fFSleptonP4;
        outtruth.fFShadSystP4 = gt.fFShadSystP4;
  
        outtruth.fIsCharm = gt.fIsCharm;
        outtruth.fCharmHadronPdg = gt.fCharmHadronPdg;
        outtruth.fIsStrange = gt.fIsStrange;
        outtruth.fStrangeHadronPdg = gt.fStrangeHadronPdg;
        outtruth.fNumProton = gt.fNumProton;
        outtruth.fNumNeutron = gt.fNumNeutron;
        outtruth.fNumPi0 = gt.fNumPi0;
        outtruth.fNumPiPlus = gt.fNumPiPlus;
        outtruth.fNumPiMinus = gt.fNumPiMinus;
        outtruth.fResNum = gt.fResNum;
        outtruth.fDecayMode = gt.fDecayMode;
    }

std::set<std::string> const& gar::IgnorableProviderConfigKeys ( )

inline

Returns a list of configuration keys that providers should ignore.

---

Returns

a reference to a key list

This function may be used for parameter validation, like in:

fhicl::Table<Config> cfg { pset, gar::IgnorableProviderConfigKeys() };

where pset is a fhicl::ParameterSet. This will inform cfg that some keys can be unexpectedly present, or missing.

This implementation includes:

• art framework service keywords

Definition at line 35 of file ProviderUtil.h.

    {
      static std::set<std::string> const ignorable {
        "service_type",     // added by art: service name (possibly interface)
        "service_provider"  // art: service implementation name
        };
      return ignorable;
    } // IgnorableProviderConfigKeys()

garana::CaloCluster gar::MakeAnaCalCluster	(	const rec::Cluster &	clust,
		const int &	region,
		const std::vector< std::pair< int, float >> &	edeps
	)

Definition at line 48 of file AnaUtils.cxx.

                                                                                                                                {

       TLorentzVector pos( clust.Position()[0], clust.Position()[1], clust.Position()[2], clust.Time() );

       //trade in our arrays for vectors/TVector3s
       const float* eigs = clust.EigenVectors(); 
       std::vector<TVector3> vecs;
       for(int i=0; i<3; i++){
           TVector3 eig( eigs[i], eigs[i+1], eigs[i+2] );
           vecs.push_back(eig);
       }

       garana::CaloCluster outclust(pos, region, clust.Energy(), clust.EnergyError(), clust.TimeDiffFirstLast(),
                                    clust.Shape(), clust.ITheta(), clust.IPhi(), vecs, edeps );
        
       return outclust;
    }//

garana::CaloCluster gar::MakeAnaCalCluster	(	const rec::Cluster &	clust,
		const int &	region,
		const vector< pair< int, float >> &	edeps
	)

Definition at line 48 of file AnaUtils.cxx.

                                                                                                                                {

       TLorentzVector pos( clust.Position()[0], clust.Position()[1], clust.Position()[2], clust.Time() );

       //trade in our arrays for vectors/TVector3s
       const float* eigs = clust.EigenVectors(); 
       std::vector<TVector3> vecs;
       for(int i=0; i<3; i++){
           TVector3 eig( eigs[i], eigs[i+1], eigs[i+2] );
           vecs.push_back(eig);
       }

       garana::CaloCluster outclust(pos, region, clust.Energy(), clust.EnergyError(), clust.TimeDiffFirstLast(),
                                    clust.Shape(), clust.ITheta(), clust.IPhi(), vecs, edeps );
        
       return outclust;
    }//

garana::GTruth gar::MakeAnaGTruth	(	const simb::GTruth &	gt,
		const int &	vtxregion
	)

Definition at line 17 of file AnaUtils.cxx.

                                                                           {

        garana::GTruth outtruth;
        FillGTruth(gt, outtruth);

        outtruth.fVertexRegion = vtxregion;

        return outtruth;

    }

const garana::Track gar::MakeAnaTrack	(	const rec::Track &	trk,
		const vector< pair< int, float >> &	pidf,
		const vector< pair< int, float >> &	pidb,
		float	ionf,
		float	ionb,
		const vector< pair< UInt_t, TLorentzVector >> &	posBeg,
		const vector< pair< UInt_t, TLorentzVector >> &	posEnd,
		const vector< pair< UInt_t, TLorentzVector >> &	momBeg,
		const vector< pair< UInt_t, TLorentzVector >> &	momEnd,
		const vector< pair< int, float >> &	edeps
	)

Definition at line 202 of file AnaUtils.cxx.

    {

        const float* tmp = trk.Vertex();
        TLorentzVector vtx(tmp[0], tmp[1], tmp[2], trk.Time());

        tmp = trk.End();
        TLorentzVector end(tmp[0], tmp[1], tmp[2], trk.Time()); // eventually track end times should be different

        tmp = trk.VtxDir();
        TVector3 vtxDir(tmp[0], tmp[1], tmp[2]);

        tmp = trk.EndDir();
        TVector3 endDir(tmp[0], tmp[1], tmp[2]);

        garana::Track anatrk( trk.LengthForward(), trk.LengthBackward(), trk.Momentum_beg(), trk.Momentum_end(),
                        vtx, end, vtxDir, endDir, trk.ChisqForward(), trk.ChisqBackward(), trk.NHits(),
                        trk.TrackParBeg(), trk.TrackParEnd(), trk.CovMatBegPacked(), trk.CovMatEndPacked(),
                        trk.ChargeBeg(), trk.ChargeEnd(), pidf, pidb, ionf, ionb, posBeg, posEnd, momBeg, 
                        momEnd, edeps);
 
        return anatrk;

    }

garana::Vee gar::MakeAnaVee

(

const rec::Vee &

vee

)

Definition at line 67 of file AnaUtils.cxx.

                                                    {

       TLorentzVector vtx( vee.Position()[0], vee.Position()[1], vee.Position()[2], vee.Time() );
       std::vector<TLorentzVector> vecs;
       for(int i=0; i<3; i++) {
           vecs.push_back(vee.FourMomentum(i));
       }

       garana::Vee outvee( vtx, vecs, vee.Chisq(), vee.VertexCov() );
        
       return outvee;
    }//

garana::Vertex gar::MakeAnaVtx

(

const rec::Vertex &

vtx

)

Definition at line 80 of file AnaUtils.cxx.

                                                       {

      TLorentzVector pos( vtx.Position()[0], vtx.Position()[1], vtx.Position()[2], vtx.Time() );
      garana::Vertex outvtx(pos, vtx.CovMat());

      return outvtx;
    }//

garana::FSParticle gar::MakeFSParticle

(

const simb::MCParticle &

mcp

)

Definition at line 29 of file AnaUtils.cxx.

                                                                {
        return garana::FSParticle(mcp.TrackId(), mcp.PdgCode(),mcp.Position(),mcp.Momentum());
    }

garana::G4Particle gar::MakeG4Particle	(	const simb::MCParticle &	mcp,
		int	parentPdg,
		int	progenitorPdg,
		int	progenitorTrackId,
		const vector< pair< TLorentzVector, TLorentzVector >> &	positions,
		const vector< pair< TLorentzVector, TLorentzVector >> &	momenta,
		const vector< int > &	regions,
		const vector< size_t > &	nptsPerRegion
	)

Definition at line 33 of file AnaUtils.cxx.

                                                                            {
        return garana::G4Particle(
            mcp.NumberTrajectoryPoints(), mcp.PdgCode(), parentPdg, progenitorPdg,
            mcp.TrackId(), mcp.Mother(), progenitorTrackId,
            gar::ProcessNameToCode(mcp.Process()),
            gar::ProcessNameToCode(mcp.EndProcess()),
            positions, momenta, regions, nptsPerRegion
          );
    }

template<typename... Providers>

ProviderPack<Providers...> gar::makeProviderPack

(

Providers const *...

providers

)

Function to create a ParameterPack from the function arguments.

Template Parameters

Providers

types of the providers in the parameter pack

Parameters

providers

constant pointers to the providers

Returns

a ParameterPack object containing all the specified providers

This is an convevience function to reduce the typing needed to instantiate a ParameterPack.

Definition at line 182 of file ProviderPack.h.

    { return ProviderPack<Providers...>(providers...); }

int gar::ProcessNameToCode

(

std::string const &

p

)

Definition at line 137 of file AnaUtils.cxx.

    {
        if(processMap.find(p)==processMap.end()){
            std::cout << "WARNING(AnaUtils::ProcessNameToCode): "
                      << "unknown process name, '" << p << "'" << std::endl;
            return -9999;  
        }
        
        else 
            return processMap[p];
  
    }

template<class T >

T::provider_type const* gar::providerFrom

(

)

Returns a constant pointer to the provider of specified service.

---

Template Parameters

T	type of the service

Returns

a constant pointer to the provider of specified service

Exceptions

art::Exception

(category art::errors::NotFound) if pointer is null

This function relies on the following service and provider interface:

• provider is not movable nor copiable
• service contains a type "provider_type" defined as the class of the service provider
• service contains a method "provider()" that returns a non-null pointer to a service provider; the service provider is owned and managed by the service, and the caller is not responsible of regulating the object lifetime, nor it should attempt to

Violations of the protocol yield compilation errors (in case non-compliance can be statically detected), or throw of exceptions.

Example of usage:

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

retrieves the service provider for GArSoft geometry. This requires the inclusion of "Geometry/GeometryGAr.h" header, where the service is declared. Typically, both ServiceUtil.h and the header of the provider class are included in the service header.

Definition at line 63 of file ServiceUtil.h.

    {
      details::ServiceRequirementsChecker<T>(); // instantiate a temporary...
      
      // retrieve the provider
      ::art::ServiceHandle<T> h;
      typename T::provider_type const* pProvider = h->provider();
      if (!pProvider) {
        throw ::art::Exception(::art::errors::NotFound)
          << "Service <" << cet::demangle_symbol(typeid(T).name())
          << "> offered a null provider";
      }
      
      return pProvider;
      
    } // providerFrom()

G4double gar::UnivScreenFunc	(	G4double	E,
		G4double	Z,
		G4double	A
	)

Definition at line 1258 of file NestAlg.cxx.

                                                                 {
    G4double a_0 = 5.29e-11*CLHEP::m; G4double a = 0.626*a_0*pow(Z,(-1./3.));
    G4double epsilon_0 = 8.854e-12*(CLHEP::farad/CLHEP::m);
    G4double epsilon = (a*E*2.*CLHEP::twopi*epsilon_0)/(2*pow(CLHEP::eplus,2.)*pow(Z,2.));
    G4double zeta_0 = pow(Z,(1./6.)); G4double m_N = A*1.66e-27*CLHEP::kg;
    G4double hbar = 6.582e-16*CLHEP::eV*CLHEP::s;
    if ( Z == 54 ) {
      epsilon *= 1.068; //zeta_0 = 1.63;
    } //special case for LXe from Bezrukov et al. 2011
    G4double s_n = log(1+1.1383*epsilon)/(2.*(epsilon +
                                              0.01321*pow(epsilon,0.21226) +
                                              0.19593*sqrt(epsilon)));
    G4double s_e = (a_0*zeta_0/a)*hbar*sqrt(8*epsilon*2.*CLHEP::twopi*epsilon_0/
                                            (a*m_N*pow(CLHEP::eplus,2.)));
    return 1.38e5*0.5*(1+tanh(50*epsilon-0.25))*epsilon*(s_e/s_n);
  }

Variable Documentation

G4double gar::biExc = 0.77

Definition at line 47 of file NestAlg.cxx.

G4bool gar::diffusion = true

Definition at line 43 of file NestAlg.cxx.

G4bool gar::OutElectrons =true

Definition at line 45 of file NestAlg.cxx.

std::unordered_map<std::string,int> gar::processMap

Definition at line 89 of file AnaUtils.cxx.

G4bool gar::SinglePhase =false

Definition at line 45 of file NestAlg.cxx.

G4bool gar::ThomasImelTail =true

Definition at line 45 of file NestAlg.cxx.