TBPositionCalc.h

Go to the documentation of this file.

#ifndef TBPositionCalc_h
#define TBPositionCalc_h

#include "CLHEP/Vector/ThreeVector.h"
#include "CLHEP/Vector/Rotation.h"
#include "Rtypes.h"

#include <fstream>
#include <vector>
#include <cmath>
#include <map>

#include "DataFormats/EcalRecHit/interface/EcalRecHit.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/Math/interface/Point3D.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "Geometry/CaloGeometry/interface/TruncatedPyramid.h"
#include "Geometry/EcalTestBeam/interface/EcalTBCrystalMap.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CLHEP/Units/GlobalSystemOfUnits.h" 


class TBPositionCalc
{
 public:
  
  TBPositionCalc(const std::map<std::string,double>& providedParameters, const std::string& mapFile, const CaloSubdetectorGeometry *passedGeometry);  

  TBPositionCalc() { };

  ~TBPositionCalc();

  CLHEP::Hep3Vector CalculateTBPos(const std::vector<EBDetId>& passedDetIds, int myCrystal, EcalRecHitCollection const *passedRecHitsMap);
  
  CLHEP::Hep3Vector CalculateCMSPos(const std::vector<EBDetId>& passedDetIds, int myCrystal, EcalRecHitCollection const *passedRecHitsMap);
  
  void computeRotation(int myCrystal, CLHEP::HepRotation & CMStoTB );
    

 private:
  bool        param_LogWeighted_;
  Double32_t  param_X0_;
  Double32_t  param_T0_;
  Double32_t  param_W0_;

  EcalTBCrystalMap * theTestMap;

  const CaloSubdetectorGeometry *theGeometry_;
};
  
#endif