BeamMCloss Class Reference

Job module which evaluates track energy loss in GEANT. More...

#include <BeamMCloss.h>

Inheritance diagram for BeamMCloss:

List of all members.

Public Member Functions
	BeamMCloss (const char *version)
	~BeamMCloss ()
JobCResult	Reco (EDMEventHandle &evt)
void	Update (const CfgConfig &)
Private Attributes
TProfile *	fELost [6]
TH2F *	fELost2D [6]

---

Detailed Description

Job module which evaluates track energy loss in GEANT.

Definition at line 18 of file BeamMCloss.h.

---

Constructor & Destructor Documentation

BeamMCloss::BeamMCloss

(

const char *

version

)

Definition at line 28 of file BeamMCloss.cxx.

References name.

                                          :
  JobCModule("BeamMCloss")
{
  SetCfgVersion(version);
  char name[256], title[256];
  for (int i = 0; i < 6; ++i) {
    sprintf(name, "elost%d", i + 1);
    sprintf(title, "Fraction of momentum lost in Chamber %d", i + 1);
    fELost[i] = new TProfile(name, title, 120, 0, 120);

    sprintf(name, "elost2d%d", i + 1);
    fELost2D[i] = new TH2F(name, title, 120, 0, 120, 1000, 0, 0.1);
  }
}

BeamMCloss::~BeamMCloss

(

)

[inline]

Definition at line 21 of file BeamMCloss.h.

{;}

---

Member Function Documentation

JobCResult BeamMCloss::Reco

(

EDMEventHandle &

evt

)

[virtual]

Reimplemented from JobCModule.

Definition at line 45 of file BeamMCloss.cxx.

References cham, GMIPPGeo::Cham(), fELost, fELost2D, geo, EDMEventHandle::Header(), GMIPPGeo::Instance(), EDMEventHeader::IsMC(), kDC1, JobCModule::kFailed, JobCModule::kPassed, kPWC6, MCTruth::Ntrack(), p, MCTrack::Particle(), MCCParticle::Q(), GMIPPGeo::TPC(), MCTruth::Track(), GChamGeo::Z(), and GTPCGeo::ZMin().

{
  if (!evt.Header().IsMC()) return kFailed;
  MCTruth mct(evt);

  for (int i = 0; i < mct.Ntrack(); ++i) {
    const MCTrack& trk = mct.Track(i);

    const MCCParticle& part = trk.Particle();
    // Make sure that track is not empty and came from a charged particle
    if (part.Q() == 0 || trk.empty()) continue;

    // Use only tracks that
    if (trk.back()->Z() < GMIPPGeo::Instance().Cham(kPWC6).Z(1) - 1) 
      continue;

    if (part.Vz() > GMIPPGeo::Instance().TPC().ZMin() + 20) continue;

    for (unsigned int j = 0; j < trk.size(); ++j) {
      const MCCHit* hit = trk[j];
      for (int cham = kPWC6; cham >= kDC1; cham--) {
    GChamGeo& geo = GMIPPGeo::Instance().Cham(cham);
    if (hit->Z() < geo.Z(4) + 1 && 
        hit->Z() > geo.Z(1) - 1) {
      double p = 1 - sqrt(hit->Px() * hit->Px() +
                  hit->Py() * hit->Py() +
                  hit->Pz() * hit->Pz()) / part.P();

      // Do not use outliers
      if (p < 0.1) {
        fELost[cham - 3]->Fill(part.P(), p);
        fELost2D[cham - 3]->Fill(part.P(), p);
      }

      break;
    }
      }
    }
      
  }

  return kPassed;
}

void BeamMCloss::Update

(

const CfgConfig &

)

[inline, virtual]

Implements CfgObserver.

Definition at line 24 of file BeamMCloss.h.

{;}

---

Member Data Documentation

TProfile* BeamMCloss::fELost[6] [private]

Definition at line 27 of file BeamMCloss.h.

Referenced by Reco().

TH2F* BeamMCloss::fELost2D[6] [private]

Definition at line 28 of file BeamMCloss.h.

Referenced by Reco().

---

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

• BeamMCloss.h
• BeamMCloss.cxx

---