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 40 of file MultiEMShowers_module.cc.

Constructor & Destructor Documentation

ems::MCinfo::MCinfo

(

const art::Event &

evt

)

Definition at line 154 of file MultiEMShowers_module.cc.

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

Member Function Documentation

void ems::MCinfo::Findtpcborders

(

const art::Event &

evt

)

Definition at line 161 of file MultiEMShowers_module.cc.

{
  art::ServiceHandle<geo::Geometry const> 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 69 of file MultiEMShowers_module.cc.

  {
    return fCosine;
  }

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

inline

Definition at line 96 of file MultiEMShowers_module.cc.

  {
    return fDirgamma1;
  }

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

inline

Definition at line 101 of file MultiEMShowers_module.cc.

  {
    return fDirgamma2;
  }

double ems::MCinfo::GetMomGamma1 ( ) const

inline

Definition at line 58 of file MultiEMShowers_module.cc.

  {
    return fGammamom1;
  }

double ems::MCinfo::GetMomGamma2 ( ) const

inline

Definition at line 63 of file MultiEMShowers_module.cc.

  {
    return fGammamom2;
  }

double ems::MCinfo::GetMompi0 ( ) const

inline

Definition at line 53 of file MultiEMShowers_module.cc.

  {
    return fMompi0;
  }

int ems::MCinfo::GetNgammas ( ) const

inline

Definition at line 47 of file MultiEMShowers_module.cc.

  {
    return fNgammas;
  }

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

inline

Definition at line 85 of file MultiEMShowers_module.cc.

  {
    return fConvgamma1;
  }

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

inline

Definition at line 90 of file MultiEMShowers_module.cc.

  {
    return fConvgamma2;
  }

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

inline

Definition at line 80 of file MultiEMShowers_module.cc.

  {
    return fPi0pos;
  }

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

inline

Definition at line 75 of file MultiEMShowers_module.cc.

  {
    return fPrimary;
  }

void ems::MCinfo::Info

(

const art::Event &

evt

)

Definition at line 183 of file MultiEMShowers_module.cc.

{
  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 const> 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 262 of file MultiEMShowers_module.cc.

{

  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 118 of file MultiEMShowers_module.cc.

  {
    return fCompton;
  }

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

inline

Definition at line 107 of file MultiEMShowers_module.cc.

  {
    return fInside1;
  }

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

inline

Definition at line 112 of file MultiEMShowers_module.cc.

  {
    return fInside2;
  }

Member Data Documentation

bool ems::MCinfo::fCompton

private

Definition at line 144 of file MultiEMShowers_module.cc.

TVector3 ems::MCinfo::fConvgamma1

private

Definition at line 148 of file MultiEMShowers_module.cc.

TVector3 ems::MCinfo::fConvgamma2

private

Definition at line 149 of file MultiEMShowers_module.cc.

double ems::MCinfo::fCosine

private

Definition at line 142 of file MultiEMShowers_module.cc.

TVector3 ems::MCinfo::fDirgamma1

private

Definition at line 150 of file MultiEMShowers_module.cc.

TVector3 ems::MCinfo::fDirgamma2

private

Definition at line 151 of file MultiEMShowers_module.cc.

double ems::MCinfo::fFidVolCut

private

Definition at line 132 of file MultiEMShowers_module.cc.

double ems::MCinfo::fGammamom1

private

Definition at line 137 of file MultiEMShowers_module.cc.

double ems::MCinfo::fGammamom2

private

Definition at line 139 of file MultiEMShowers_module.cc.

bool ems::MCinfo::fInside1

private

Definition at line 138 of file MultiEMShowers_module.cc.

bool ems::MCinfo::fInside2

private

Definition at line 140 of file MultiEMShowers_module.cc.

double ems::MCinfo::fMaxx

private

Definition at line 126 of file MultiEMShowers_module.cc.

double ems::MCinfo::fMaxy

private

Definition at line 128 of file MultiEMShowers_module.cc.

double ems::MCinfo::fMaxz

private

Definition at line 130 of file MultiEMShowers_module.cc.

double ems::MCinfo::fMinx

private

Definition at line 125 of file MultiEMShowers_module.cc.

double ems::MCinfo::fMiny

private

Definition at line 127 of file MultiEMShowers_module.cc.

double ems::MCinfo::fMinz

private

Definition at line 129 of file MultiEMShowers_module.cc.

double ems::MCinfo::fMompi0

private

Definition at line 136 of file MultiEMShowers_module.cc.

int ems::MCinfo::fNgammas

private

Definition at line 134 of file MultiEMShowers_module.cc.

TVector3 ems::MCinfo::fPi0pos

private

Definition at line 147 of file MultiEMShowers_module.cc.

TVector3 ems::MCinfo::fPrimary

private

Definition at line 146 of file MultiEMShowers_module.cc.

---

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

• MultiEMShowers_module.cc