ems::MCinfo Class Reference

Public Member Functions
	MCinfo (const art::Event &evt)
void	Info (const art::Event &evt)
void	Findtpcborders (const art::Event &evt)
int	GetNgammas () const
double	GetMompi0 () const
double	GetMomGamma1 () const
double	GetMomGamma2 () const
double	GetCosine ()
TVector3 const &	GetPrimary () const
TVector3 const &	GetPospi0 () const
TVector3 const &	GetPosgamma1 () const
TVector3 const &	GetPosgamma2 () const
TVector3 const &	GetDirgamma1 () const
TVector3 const &	GetDirgamma2 () const
bool const &	IsInside1 () const
bool const &	IsInside2 () const
bool const &	IsCompton () const

Private Member Functions
bool	insideFidVol (const TLorentzVector &pvtx) const

Private Attributes
double	fMinx
double	fMaxx
double	fMiny
double	fMaxy
double	fMinz
double	fMaxz
double	fFidVolCut
int	fNgammas
double	fMompi0
double	fGammamom1
bool	fInside1
double	fGammamom2
bool	fInside2
double	fCosine
bool	fCompton
TVector3	fPrimary
TVector3	fPi0pos
TVector3	fConvgamma1
TVector3	fConvgamma2
TVector3	fDirgamma1
TVector3	fDirgamma2

Detailed Description

Definition at line 49 of file MultiEMShowers_module.cc.

Constructor & Destructor Documentation

ems::MCinfo::MCinfo

(

const art::Event &

evt

)

Definition at line 105 of file MultiEMShowers_module.cc.

References Findtpcborders(), and Info().

                                     :
fFidVolCut(2.0)
{
        Info(evt);
        Findtpcborders(evt);
}

Member Function Documentation

void ems::MCinfo::Findtpcborders

(

const art::Event &

evt

)

Definition at line 112 of file MultiEMShowers_module.cc.

References fMaxx, fMaxy, fMaxz, fMinx, fMiny, fMinz, and geo::GeometryCore::IterateTPCs().

Referenced by MCinfo().

{
        art::ServiceHandle<geo::Geometry> geom;
        
        fMinx = geom->IterateTPCs().begin()->MinX();
        fMiny = geom->IterateTPCs().begin()->MinY();
        fMinz = geom->IterateTPCs().begin()->MinZ();
        fMaxx = geom->IterateTPCs().begin()->MaxX();
        fMaxy = geom->IterateTPCs().begin()->MaxY();
        fMaxz = geom->IterateTPCs().begin()->MaxZ();

        for (const geo::TPCGeo& tpcg: geom->IterateTPCs())
        {
                if (tpcg.MinX() < fMinx) fMinx = tpcg.MinX();
                if (tpcg.MaxX() > fMaxx) fMaxx = tpcg.MaxX(); 
                if (tpcg.MinY() < fMiny) fMiny = tpcg.MinY();
                if (tpcg.MaxY() > fMaxy) fMaxy = tpcg.MaxY();
                if (tpcg.MinZ() < fMinz) fMinz = tpcg.MinZ();
                if (tpcg.MaxZ() > fMaxz) fMaxz = tpcg.MaxZ();
        }
}

double ems::MCinfo::GetCosine ( )

inline

Definition at line 62 of file MultiEMShowers_module.cc.

References fCosine.

Referenced by ems::MultiEMShowers::analyze().

{ return fCosine; }

TVector3 const& ems::MCinfo::GetDirgamma1 ( ) const

inline

Definition at line 69 of file MultiEMShowers_module.cc.

References fDirgamma1.

Referenced by ems::MultiEMShowers::analyze().

{ return fDirgamma1; }

TVector3 const& ems::MCinfo::GetDirgamma2 ( ) const

inline

Definition at line 70 of file MultiEMShowers_module.cc.

References fDirgamma2.

Referenced by ems::MultiEMShowers::analyze().

{ return fDirgamma2; }

double ems::MCinfo::GetMomGamma1 ( ) const

inline

Definition at line 59 of file MultiEMShowers_module.cc.

References fGammamom1.

Referenced by ems::MultiEMShowers::analyze().

{ return fGammamom1; }

double ems::MCinfo::GetMomGamma2 ( ) const

inline

Definition at line 60 of file MultiEMShowers_module.cc.

References fGammamom2.

Referenced by ems::MultiEMShowers::analyze().

{ return fGammamom2; }

double ems::MCinfo::GetMompi0 ( ) const

inline

Definition at line 58 of file MultiEMShowers_module.cc.

References fMompi0.

Referenced by ems::MultiEMShowers::analyze().

{ return fMompi0; }

int ems::MCinfo::GetNgammas ( ) const

inline

Definition at line 56 of file MultiEMShowers_module.cc.

References fNgammas.

Referenced by ems::MultiEMShowers::analyze().

{ return fNgammas; }

TVector3 const& ems::MCinfo::GetPosgamma1 ( ) const

inline

Definition at line 66 of file MultiEMShowers_module.cc.

References fConvgamma1.

Referenced by ems::MultiEMShowers::analyze(), and ems::MultiEMShowers::convCluster().

{ return fConvgamma1; }

TVector3 const& ems::MCinfo::GetPosgamma2 ( ) const

inline

Definition at line 67 of file MultiEMShowers_module.cc.

References fConvgamma2.

Referenced by ems::MultiEMShowers::analyze(), and ems::MultiEMShowers::convCluster().

{ return fConvgamma2; }

TVector3 const& ems::MCinfo::GetPospi0 ( ) const

inline

Definition at line 65 of file MultiEMShowers_module.cc.

References fPi0pos.

Referenced by ems::MultiEMShowers::analyze().

{ return fPi0pos; }

TVector3 const& ems::MCinfo::GetPrimary ( ) const

inline

Definition at line 64 of file MultiEMShowers_module.cc.

References fPrimary.

Referenced by ems::MultiEMShowers::analyze().

{ return fPrimary; }

void ems::MCinfo::Info

(

const art::Event &

evt

)

Definition at line 134 of file MultiEMShowers_module.cc.

References sim::ParticleList::begin(), simb::MCParticle::Daughter(), sim::ParticleList::end(), simb::MCParticle::EndProcess(), fCompton, fConvgamma1, fConvgamma2, fCosine, fDirgamma1, fDirgamma2, fGammamom1, fGammamom2, fInside1, fInside2, fMompi0, fNgammas, fPi0pos, fPrimary, insideFidVol(), simb::MCParticle::NumberDaughters(), simb::MCParticle::P(), cheat::ParticleInventoryService::ParticleList(), simb::MCParticle::PdgCode(), simb::MCParticle::Position(), simb::MCParticle::Process(), and cheat::ParticleInventoryService::TrackIdToParticle_P().

Referenced by MCinfo().

{
        fMompi0 = 0.0; fPi0pos.SetXYZ(0,0,0); 
        fNgammas = 0;
        fCosine = 0.0;
        fInside1 = false; fInside2 = false;
        fCompton = false;

        fGammamom1 = 0.0; fGammamom2 = 0.0;
        fConvgamma1.SetXYZ(0,0,0); fConvgamma2.SetXYZ(0,0,0); 
        fDirgamma1.SetXYZ(0,0,0); fDirgamma2.SetXYZ(0,0,0);

        art::ServiceHandle< cheat::ParticleInventoryService > pi_serv;
        const sim::ParticleList& plist = pi_serv->ParticleList();
        for (sim::ParticleList::const_iterator ipar = plist.begin(); ipar != plist.end(); ++ipar)
        {
                const simb::MCParticle* particle = ipar->second;

                if (particle->Process() != "primary") continue;

                TLorentzVector posvec = particle->Position();
                TVector3 pose(posvec.X(), posvec.Y(), posvec.Z());
                fPrimary = pose;
                

                if (particle->PdgCode() == 111)
                {
                        fMompi0 = particle->P();
                        
                        TLorentzVector posvec3 = particle->Position();
                        TVector3 pospi0(posvec3.X(), posvec3.Y(), posvec3.Z());
                        fPi0pos =  pospi0;
        
                        if (particle->NumberDaughters() != 2) continue;

                        const simb::MCParticle* daughter1 = pi_serv->TrackIdToParticle_P(particle->Daughter(0));
                        if (daughter1->PdgCode() != 22) continue;

                        const simb::MCParticle* daughter2 = pi_serv->TrackIdToParticle_P(particle->Daughter(1));
                        if (daughter2->PdgCode() != 22) continue; 

                        fNgammas = particle->NumberDaughters();
                        TLorentzVector mom1 = pi_serv->TrackIdToParticle_P(particle->Daughter(0))->Momentum();
                        TLorentzVector mom2 = pi_serv->TrackIdToParticle_P(particle->Daughter(1))->Momentum();

                        // compton process
                        if (daughter1->EndProcess() == "phot") fCompton = true;
                        if (daughter2->EndProcess() == "phot") fCompton = true;

                        TVector3 mom1vec3(mom1.Px(), mom1.Py(), mom1.Pz());
                        fGammamom1 = pi_serv->TrackIdToParticle_P(particle->Daughter(0))->P();
                        TVector3 mom2vec3(mom2.Px(), mom2.Py(), mom2.Pz());
                        fGammamom2 = pi_serv->TrackIdToParticle_P(particle->Daughter(1))->P();

                        TLorentzVector pos1 = pi_serv->TrackIdToParticle_P(particle->Daughter(0))->EndPosition();
                        TLorentzVector pos2 = pi_serv->TrackIdToParticle_P(particle->Daughter(1))->EndPosition();
                                
                        if (insideFidVol(pos1)) fInside1 = true; 
                        if (insideFidVol(pos2)) fInside2 = true;
                                
                        
                        fConvgamma1.SetXYZ(pos1.X(), pos1.Y(), pos1.Z());
                        fConvgamma2.SetXYZ(pos2.X(), pos2.Y(), pos2.Z());

                        TVector3 vecnorm1 = mom1vec3.Unit();
                        fDirgamma1 = vecnorm1;
                        TVector3 vecnorm2 = mom2vec3.Unit();
                        fDirgamma2 = vecnorm2;
                
                        fCosine = fDirgamma1 * fDirgamma2;      
                }
                else
                {
                        fNgammas = particle->NumberDaughters();
                }
        }
}

bool ems::MCinfo::insideFidVol ( const TLorentzVector &  pvtx ) const

private

Definition at line 212 of file MultiEMShowers_module.cc.

References fFidVolCut, fMaxx, fMaxy, fMaxz, fMinx, fMiny, and fMinz.

Referenced by Info(), and IsCompton().

{

        bool inside = false;
        //x
        double dista = fabs(fMinx - pvtx.X());
        double distb = fabs(pvtx.X() - fMaxx); 
        if ((pvtx.X() > fMinx) && (pvtx.X() < fMaxx) &&
                        (dista > fFidVolCut) && (distb > fFidVolCut)) inside = true;
        //y
        dista = fabs(fMaxy - pvtx.Y());
        distb = fabs(pvtx.Y() - fMiny);
        if (inside && (pvtx.Y() > fMiny) && (pvtx.Y() < fMaxy) &&
                (dista > fFidVolCut) && (distb > fFidVolCut)) inside = true;
        else inside = false;

        //z
        dista = fabs(fMaxz - pvtx.Z());
        distb = fabs(pvtx.Z() - fMinz);
        if (inside && (pvtx.Z() > fMinz) && (pvtx.Z() < fMaxz) &&
                (dista > fFidVolCut) && (distb > fFidVolCut)) inside = true;
        else inside = false;
                
        return inside;
}

bool const& ems::MCinfo::IsCompton ( ) const

inline

Definition at line 75 of file MultiEMShowers_module.cc.

References fCompton, and insideFidVol().

Referenced by ems::MultiEMShowers::analyze().

{ return fCompton; }

bool const& ems::MCinfo::IsInside1 ( ) const

inline

Definition at line 72 of file MultiEMShowers_module.cc.

References fInside1.

Referenced by ems::MultiEMShowers::analyze().

{ return fInside1; }

bool const& ems::MCinfo::IsInside2 ( ) const

inline

Definition at line 73 of file MultiEMShowers_module.cc.

References fInside2.

Referenced by ems::MultiEMShowers::analyze().

{ return fInside2; }

Member Data Documentation

bool ems::MCinfo::fCompton

private

Definition at line 95 of file MultiEMShowers_module.cc.

Referenced by Info(), and IsCompton().

TVector3 ems::MCinfo::fConvgamma1

private

Definition at line 99 of file MultiEMShowers_module.cc.

Referenced by GetPosgamma1(), and Info().

TVector3 ems::MCinfo::fConvgamma2

private

Definition at line 100 of file MultiEMShowers_module.cc.

Referenced by GetPosgamma2(), and Info().

double ems::MCinfo::fCosine

private

Definition at line 93 of file MultiEMShowers_module.cc.

Referenced by GetCosine(), and Info().

TVector3 ems::MCinfo::fDirgamma1

private

Definition at line 101 of file MultiEMShowers_module.cc.

Referenced by GetDirgamma1(), and Info().

TVector3 ems::MCinfo::fDirgamma2

private

Definition at line 102 of file MultiEMShowers_module.cc.

Referenced by GetDirgamma2(), and Info().

double ems::MCinfo::fFidVolCut

private

Definition at line 83 of file MultiEMShowers_module.cc.

Referenced by insideFidVol().

double ems::MCinfo::fGammamom1

private

Definition at line 88 of file MultiEMShowers_module.cc.

Referenced by GetMomGamma1(), and Info().

double ems::MCinfo::fGammamom2

private

Definition at line 90 of file MultiEMShowers_module.cc.

Referenced by GetMomGamma2(), and Info().

bool ems::MCinfo::fInside1

private

Definition at line 89 of file MultiEMShowers_module.cc.

Referenced by Info(), and IsInside1().

bool ems::MCinfo::fInside2

private

Definition at line 91 of file MultiEMShowers_module.cc.

Referenced by Info(), and IsInside2().

double ems::MCinfo::fMaxx

private

Definition at line 79 of file MultiEMShowers_module.cc.

Referenced by Findtpcborders(), and insideFidVol().

double ems::MCinfo::fMaxy

private

Definition at line 80 of file MultiEMShowers_module.cc.

Referenced by Findtpcborders(), and insideFidVol().

double ems::MCinfo::fMaxz

private

Definition at line 81 of file MultiEMShowers_module.cc.

Referenced by Findtpcborders(), and insideFidVol().

double ems::MCinfo::fMinx

private

Definition at line 79 of file MultiEMShowers_module.cc.

Referenced by Findtpcborders(), and insideFidVol().

double ems::MCinfo::fMiny

private

Definition at line 80 of file MultiEMShowers_module.cc.

Referenced by Findtpcborders(), and insideFidVol().

double ems::MCinfo::fMinz

private

Definition at line 81 of file MultiEMShowers_module.cc.

Referenced by Findtpcborders(), and insideFidVol().

double ems::MCinfo::fMompi0

private

Definition at line 87 of file MultiEMShowers_module.cc.

Referenced by GetMompi0(), and Info().

int ems::MCinfo::fNgammas

private

Definition at line 85 of file MultiEMShowers_module.cc.

Referenced by GetNgammas(), and Info().

TVector3 ems::MCinfo::fPi0pos

private

Definition at line 98 of file MultiEMShowers_module.cc.

Referenced by GetPospi0(), and Info().

TVector3 ems::MCinfo::fPrimary

private

Definition at line 97 of file MultiEMShowers_module.cc.

Referenced by GetPrimary(), and Info().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/DirOfGamma/MultiEMShowers_module.cc