Mask/include/StripDetector.h

#ifndef __STRIP_DETECTOR_H
#define __STRIP_DETECTOR_H

// +z is along beam axis
// +y is vertically "downward" in the lab frame

//angles must be in radians, but distances can be whatever
//PROVIDED all input distances are the same

//Front strips from largest y to smallest y

//Back strips from lowest z to highest z

#include <cmath>
#include <vector>
#include <random>

#include "Vec3.h"
#include "Rotation.h"
#include "RandomGenerator.h"

struct StripHit
{
	int front_strip_index=-1;
	int back_strip_index=-1;
	double front_ratio=0.0;
};

class StripDetector {
public:
  
	StripDetector(int ns, double len, double wid, double cphi, double cz, double crho);
	~StripDetector();
	inline Mask::Vec3 GetFrontStripCoordinates(int stripch, int corner) { return front_strip_coords[stripch][corner]; }
	inline Mask::Vec3 GetBackStripCoordinates(int stripch, int corner) { return back_strip_coords[stripch][corner]; }
	inline Mask::Vec3 GetRotatedFrontStripCoordinates(int stripch, int corner) { return rotated_front_strip_coords[stripch][corner]; }
	inline Mask::Vec3 GetRotatedBackStripCoordinates(int stripch, int corner) { return rotated_back_strip_coords[stripch][corner]; }
	inline Mask::Vec3 GetNormRotated() { return zRot*m_norm_unrot; }

	inline void TurnOnRandomizedCoordinates() { rndmFlag = true; }
	inline void TurnOffRandomizedCoordinates() { rndmFlag = false; }

	Mask::Vec3 GetHitCoordinates(int front_stripch, double front_strip_ratio);
	StripHit GetChannelRatio(double theta, double phi);

private:
	inline bool ValidChannel(int f) { return ((f >= 0 && f < num_strips) ? true : false); };
	inline bool ValidRatio(double r) { return ((r >= -1 && r <= 1) ? true : false); };
	void CalculateCorners();

	int num_strips;
	static constexpr int num_corners = 4;

	double length; //common to all strips, hence total
	double total_width;
	double front_strip_width; //assuming equal widths
	double back_strip_length; //assuming equal widths

	double center_phi; //assuming det centered above x-axis (corresponds to zero phi)
	double center_z;
	double center_rho; //perpendicular radius from axis

	std::vector<std::vector<Mask::Vec3>> front_strip_coords, back_strip_coords;
	std::vector<std::vector<Mask::Vec3>> rotated_front_strip_coords, rotated_back_strip_coords;

	Mask::Vec3 m_norm_unrot;

	Mask::ZRotation zRot;

	std::uniform_real_distribution<double> m_uniform_fraction;

	bool rndmFlag;

};

#endif