CSCCrossGap.h

Go to the documentation of this file.

#ifndef MU_END_CROSS_GAP_H
#define MU_END_CROSS_GAP_H

#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include <vector>

class CSCCrossGap {
public:
  CSCCrossGap(double mass, float mom, LocalVector gap);
  ~CSCCrossGap(){};

  std::vector<LocalPoint> ionClusters() const { return clusters; }
  int noOfClusters() const { return clusters.size(); }
  std::vector<int> electrons() const { return electronsInClusters; }
  int noOfElectrons() const { return electronsInClusters.size(); }
  std::vector<double> stepLengths() const { return steps; }
  int noOfSteps() const { return steps.size(); }
  std::vector<float> eLossPerStep() const { return elosses; }
  int noOfElosses() const { return elosses.size(); }

  void addCluster(LocalPoint here) { clusters.push_back(here); }
  void addElectrons(int nelec = 1) { electronsInClusters.push_back(nelec); }
  void addElectronToBack() { ++electronsInClusters.back(); }

  void addStep(double step) { steps.push_back(step); }
  void addEloss(float eloss) { elosses.push_back(eloss); }

  double logGamma(double mass, float momentum);
  double logGamma() { return loggam; }
  double beta2() const { return theBeta2; }
  double gamma() const { return theGamma; }
  LocalVector gapVector() const { return theGap; }
  LocalVector unitVector() const { return theGap.unit(); }
  float length() const { return theGap.mag(); }

private:
  double theBeta2;  // Lorentz beta^2
  double theGamma;  // Lorentz gamma
  double loggam;
  LocalVector theGap;

  std::vector<LocalPoint> clusters;
  std::vector<int> electronsInClusters;
  std::vector<double> steps;
  std::vector<float> elosses;
};

#endif