SingleParticleEvent.h

Go to the documentation of this file.

#ifndef SingleParticleEvent_h
#define SingleParticleEvent_h
//
// SingleParticleEvent by droll (23/DEC/2005)
//

// include files
#include "GeneratorInterface/CosmicMuonGenerator/interface/CosmicMuonParameters.h"
#include "GeneratorInterface/CosmicMuonGenerator/interface/Point5MaterialMap.h"
#include <iostream>

class SingleParticleEvent {
public:
  // constructor
  SingleParticleEvent() {
    ID = 0;
    Px = 0.;
    Py = 0.;
    Pz = 0.;
    E = 0.;
    M = 0.;
    Vx = 0.;
    Vy = 0.;
    Vz = 0.;
    T0 = 0.;
    ID_in = 0;
    Px_in = 0.;
    Py_in = 0.;
    Pz_in = 0.;
    E_in = 0.;
    M_in = 0.;
    Vx_in = 0.;
    Vy_in = 0.;
    Vz_in = 0.;
    T0_in = 0.;
    HitTarget = false;
    PlugVx = PlugOnShaftVx;
    PlugVz = PlugOnShaftVz;
    RhoAir = 0.;
    RhoWall = 0.;
    RhoRock = 0.;
    RhoClay = 0.;
    RhoPlug = 0.;
    ClayWidth = DefaultClayWidth;
  }
  // destructor
  ~SingleParticleEvent() {}

private:
  int ID;
  double Px;
  double Py;
  double Pz;
  double E;
  double M;
  double Vx;
  double Vy;
  double Vz;
  double T0;
  int ID_in;
  double Px_in;
  double Py_in;
  double Pz_in;
  double E_in;
  double M_in;
  double Vx_in;
  double Vy_in;
  double Vz_in;
  double T0_in;
  bool HitTarget;
  bool MTCC;

  // other stuff
  double dX;
  double dY;
  double dZ;
  double tmpVx;
  double tmpVy;
  double tmpVz;
  // update event during propagation
  void update(double stepSize);
  // temporary propagation
  void updateTmp(double stepSize);
  void subtractEloss(double waterEquivalents);  // update 4momentum
  double absVzTmp();                            // |Vz| [mm]
  double rVxyTmp();                             // R_XY [mm]

public:
  // create (initialize) an event with a single particle
  void create(int id, double px, double py, double pz, double e, double m, double vx, double vy, double vz, double t0);
  // propagate particle to target area
  void propagate(double ElossScaleFac,
                 double RadiusTarget,
                 double Z_DistTarget,
                 double Z_CentrTarget,
                 bool TrackerOnly,
                 bool MTCCHalf);
  double Eloss(double waterEquivalents, double Energy);  //return Eloss
  // particle has hit the target volume (during propagation)
  bool hitTarget();
  // event info (direct access)
  //initial state mother particle
  int id_in();                // [HEP particle code]
  double px_in();             // [GeV/c]
  double py_in();             // [GeV/c]
  double pz_in();             // [GeV/c]
  double e_in();              // [GeV]
  double m_in();              // [GeV/c^2]
  double vx_in();             // [mm]
  double vy_in();             // [mm]
  double vz_in();             // [mm]
  double t0_in();             // [mm/c] with c = 299.792458 mm/ns
  double WaterEquivalents();  //[g cm^-2]

  //final state daughter particles
  int id();     // [HEP particle code]
  double px();  // [GeV/c]
  double py();  // [GeV/c]
  double pz();  // [GeV/c]
  double e();   // [GeV]
  double m();   // [GeV/c^2]
  double vx();  // [mm]
  double vy();  // [mm]
  double vz();  // [mm]
  double t0();  // [mm/c] with c = 299.792458 mm/ns
  // event info (calculated)
  double phi();     // in horizontal (x-z) plane [rad]
  double theta();   // off vertical    (y) axis  [rad]
  double absmom();  // |p| [GeV/c]
  double absVz();   // |Vz| [mm]
  double rVxy();    // R_XY [mm]

  void setEug(double Eug);          // [GeV]
  double Eug();                     // [GeV/c]
  double deltaEmin(double Energy);  // [GeV]
  void SurfProj(double Vx_in,
                double Vy_in,
                double Vz_in,
                double Px_in,
                double Py_in,
                double Pz_in,
                double& Vx_up,
                double& Vy_up,
                double& Vz_up);

  double PlugVx;
  double PlugVz;
  double RhoAir;
  double RhoWall;
  double RhoRock;
  double RhoClay;
  double RhoPlug;
  double ClayWidth;
  double waterEquivalents;
  double E_ug;
};
#endif