CSCCrossGap.h

Go to the documentation of this file.

#ifndef MU_END_CROSS_GAP_H
#define MU_END_CROSS_GAP_H

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

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