PreshowerHitMaker.h

Go to the documentation of this file.

#ifndef PreshowerHitMaker_h
#define PreshowerHitMaker_h

#include "FastSimulation/CaloHitMakers/interface/CaloHitMaker.h"
#include "FastSimulation/CaloGeometryTools/interface/Transform3DPJ.h"

class CaloGeometryHelper;
class LandauFluctuationGenerator;
class RandomEngineAndDistribution;

class PreshowerHitMaker : public CaloHitMaker {
public:
  typedef math::XYZVector XYZVector;
  typedef math::XYZVector XYZPoint;
  typedef ROOT::Math::Transform3DPJ Transform3D;

  PreshowerHitMaker(CaloGeometryHelper* calo,
                    const XYZPoint&,
                    const XYZVector&,
                    const XYZPoint&,
                    const XYZVector&,
                    const LandauFluctuationGenerator* aGenerator,
                    const RandomEngineAndDistribution* engine);

  ~PreshowerHitMaker() override { ; }

  inline void setSpotEnergy(double e) override { spotEnergy = e; }
  bool addHit(double r, double phi, unsigned layer = 0) override;

  const std::map<CaloHitID, float>& getHits() override { return hitMap_; };
  // for tuning
  inline void setMipEnergy(double e1, double e2) {
    mip1_ = e1;
    mip2_ = e2;
  }

  float totalLayer1() const { return totalLayer1_; }
  float totalLayer2() const { return totalLayer2_; }
  float layer1Calibrated() const { return 0.024 / 81.1E-6 * totalLayer1_; }
  float layer2Calibrated() const { return 0.024 * 0.7 / 81.1E-6 * totalLayer2_; }
  float totalCalibrated() const { return 0.024 / 81.1E-6 * (totalLayer1_ + 0.7 * totalLayer2_); }

private:
  XYZPoint psLayer1Entrance_;
  XYZVector psLayer1Dir_;
  XYZPoint psLayer2Entrance_;
  XYZVector psLayer2Dir_;
  bool layer1valid_;
  bool layer2valid_;
  Transform3D locToGlobal1_;
  Transform3D locToGlobal2_;
  float anglecorrection1_;
  float anglecorrection2_;
  double mip1_, mip2_;
  float totalLayer1_;
  float totalLayer2_;
  const LandauFluctuationGenerator* theGenerator;
  const RandomEngineAndDistribution* random;
};

#endif