BeamSwimMC Class Reference

Module which retrieves MCTrack's, and swims charged tracks through the spectrometer,. More...

#include <BeamSwimMC.h>

Inheritance diagram for BeamSwimMC:

List of all members.

Public Member Functions
	BeamSwimMC (const char *version)
	~BeamSwimMC ()
void	NewRun (int run, int subrun)
JobCResult	Reco (EDMEventHandle &evt)
void	Update (const CfgConfig &)
Private Attributes
TTree *	fTree
BeamSwimNt *	fNt

---

Detailed Description

Module which retrieves MCTrack's, and swims charged tracks through the spectrometer,.

Definition at line 63 of file BeamSwimMC.h.

---

Constructor & Destructor Documentation

BeamSwimMC::BeamSwimMC

(

const char *

version

)

Definition at line 28 of file BeamSwimMC.cxx.

References fNt, fTree, and JobCModule::SetCfgVersion().

                                          : 
  JobCModule("BeamSwimMC")
{
  SetCfgVersion(version);

  fTree = new TTree("swimmc", "Comparison of swimmer to MC");
  fNt = new BeamSwimNt;
  fTree->Branch("fNt", "BeamSwimNt", &fNt);
}

BeamSwimMC::~BeamSwimMC

(

)

Definition at line 40 of file BeamSwimMC.cxx.

References fNt, and fTree.

{
  fTree->Write();
  delete fTree;
  delete fNt;
}

---

Member Function Documentation

void BeamSwimMC::NewRun	(	int	run,
		int	subrun
	)			[virtual]

Reimplemented from JobCModule.

Definition at line 124 of file BeamSwimMC.cxx.

References config_bfield().

{
  config_bfield(run, subrun);
}

JobCResult BeamSwimMC::Reco

(

EDMEventHandle &

evt

)

[virtual]

Reimplemented from JobCModule.

Definition at line 49 of file BeamSwimMC.cxx.

References fNt, fTree, SwimMIPP::GetXY(), SwimMIPP::kELoss, JobCModule::kPassed, SwimMIPP::kSlice, SwimMIPP::kStepS, BeamSwimNt::nPt, MCTruth::Ntrack(), BeamSwimNt::p0, MCTrack::Particle(), BeamSwimNt::ptot, BeamSwimNt::q, MCCParticle::Q(), SwimMIPP::Swim(), MCTruth::Track(), BeamSwimNt::x, x, BeamSwimNt::x0, BeamSwimNt::xelos, BeamSwimNt::xs, BeamSwimNt::xz, BeamSwimNt::y, y, BeamSwimNt::yelos, BeamSwimNt::ys, BeamSwimNt::yz, BeamSwimNt::z, and Z.

{
  MCTruth tru(evt);
  SwimMIPP swim1, swim2, swim3;

  double x, y;

  for (int itrk = 0; itrk < tru.Ntrack(); ++itrk) {
    const MCTrack& trk = tru.Track(itrk);
    if (trk.empty()) continue;
    const MCCParticle& part = trk.Particle();
    if (part.Q() == 0 || 
    part.P() < 1 ||
    part.Pz() < 0.5 ||
    part.Vz() > 1400) continue;

//     cout << "Tracking " << part.P() << " GeV/c " << part.GetName()
//   << " from " << part.Vx() << ", " << part.Vx() << ", " << part.Vz()
//   << endl;

    fNt->nPt = 0;
    fNt->q = part.Q();
    fNt->x0[0] = part.Vx();
    fNt->x0[1] = part.Vy();
    fNt->x0[2] = part.Vz();
    fNt->p0[0] = part.Px();
    fNt->p0[1] = part.Py();
    fNt->p0[2] = part.Pz();
    fNt->ptot = part.P();

    // Swim with StepS stepper
    swim1.Swim(part.Q(), part.P(), part.Vx(), part.Vy(), part.Vz(),
           part.Px() / part.Pz(), part.Py() / part.Pz(), SwimMIPP::kStepS);

    // Swim with StepZ stepper
    swim2.Swim(part.Q(), part.P(), part.Vx(), part.Vy(), part.Vz(),
           part.Px() / part.Pz(), part.Py() / part.Pz(), SwimMIPP::kSlice);
    
    // Swim with ELoss stepper
    swim3.Swim(part.Q(), part.P(), part.Vx(), part.Vy(), part.Vz(),
           part.Px() / part.Pz(), part.Py() / part.Pz(), SwimMIPP::kELoss);
    
    for (unsigned int i = 0; i < trk.size(); ++i) {
      if (i > 0 && trk[i]->Z() - fNt->z[fNt->nPt - 1] < 1) continue;

      fNt->x[fNt->nPt] = trk[i]->X();
      fNt->y[fNt->nPt] = trk[i]->Y();
      fNt->z[fNt->nPt] = trk[i]->Z();
      swim1.GetXY(trk[i]->Z(), &x, &y);
      fNt->xs[fNt->nPt] = x;
      fNt->ys[fNt->nPt] = y;

      swim2.GetXY(trk[i]->Z(), &x, &y);
      fNt->xz[fNt->nPt] = x;
      fNt->yz[fNt->nPt] = y;

      swim3.GetXY(trk[i]->Z(), &x, &y);
      fNt->xelos[fNt->nPt] = x;
      fNt->yelos[fNt->nPt] = y;

      if (++fNt->nPt >= 256) {
    cout << "Saved " << i << "/" << trk.size() << " points" 
         << endl;
    break;
      }
    }

    fTree->Fill();
  }

  return kPassed;
}

void BeamSwimMC::Update

(

const CfgConfig &

)

[inline, virtual]

Implements CfgObserver.

Definition at line 70 of file BeamSwimMC.h.

{;}

---

Member Data Documentation

BeamSwimNt* BeamSwimMC::fNt [private]

Definition at line 74 of file BeamSwimMC.h.

Referenced by BeamSwimMC(), Reco(), and ~BeamSwimMC().

TTree* BeamSwimMC::fTree [private]

Definition at line 73 of file BeamSwimMC.h.

Referenced by BeamSwimMC(), Reco(), and ~BeamSwimMC().

---

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

• BeamSwimMC.h
• BeamSwimMC.cxx

---