mirror of
https://github.com/gwm17/Mask.git
synced 2024-11-13 14:08:49 -05:00
60 lines
2.4 KiB
C
60 lines
2.4 KiB
C
|
#ifndef ANASEN_EFFICIENCY_H
|
||
|
|
#define ANASEN_EFFICIENCY_H
|
||
|
|
|
||
|
|
#include <string>
|
||
|
|
#include <THashTable.h>
|
||
|
|
|
||
|
|
#include "StripDetector.h"
|
||
|
|
#include "QQQDetector.h"
|
||
|
|
|
||
|
|
class AnasenEfficiency {
|
||
|
|
public:
|
||
|
|
AnasenEfficiency();
|
||
|
|
~AnasenEfficiency();
|
||
|
|
void CalculateEfficiency(const std::string& filename);
|
||
|
|
inline void SetReactionType(int type) { m_rxn_type = type; };
|
||
|
|
void DrawDetectorSystem(const std::string& filename);
|
||
|
|
double RunConsistencyCheck();
|
||
|
|
|
||
|
|
private:
|
||
|
|
void MyFill(THashTable* table, const std::string& name, const std::string& title, int bins, float min, float max, double val);
|
||
|
|
void MyFill(THashTable* table, const std::string& name, const std::string& title, int binsx, float minx, float maxx, double valx, int binsy, float miny, float maxy, double valy);
|
||
|
|
void RunDecay(const std::string& filename);
|
||
|
|
void RunTwoStep(const std::string& filename);
|
||
|
|
void RunThreeStep(const std::string& filename);
|
||
|
|
|
||
|
|
bool IsRing1(double theta, double phi);
|
||
|
|
bool IsRing2(double theta, double phi);
|
||
|
|
bool IsQQQ(double theta, double phi);
|
||
|
|
|
||
|
|
inline bool IsDoubleEqual(double x, double y) { return std::fabs(x-y) < epsilon ? true : false; };
|
||
|
|
|
||
|
|
int m_rxn_type;
|
||
|
|
std::vector<StripDetector> m_Ring1, m_Ring2;
|
||
|
|
std::vector<QQQDetector> m_forwardQQQs;
|
||
|
|
std::vector<QQQDetector> m_backwardQQQs;
|
||
|
|
|
||
|
|
/**** ANASEN geometry constants *****/
|
||
|
|
const int n_sx3_per_ring = 12;
|
||
|
|
const int n_qqq = 4;
|
||
|
|
const double sx3_length = 0.075;
|
||
|
|
const double sx3_width = 0.04;
|
||
|
|
const double barrel_gap = 0.013 + 0.049; //0.049 is base gap due to frames
|
||
|
|
const double ring1_z = sx3_length/2.0 + barrel_gap/2.0;
|
||
|
|
//const double ring2_z = -0.124 + sx3_length/2.0 + 0.0245 - barrel_gap/2.0;
|
||
|
|
const double qqq_nom_z = 0.025 + sx3_length + 0.0245 + barrel_gap/2.0;
|
||
|
|
const double qqq_rinner = 0.0501;
|
||
|
|
const double qqq_router = 0.0990;
|
||
|
|
const double qqq_deltaphi = 1.52119;
|
||
|
|
const double qqq_z[4] = {qqq_nom_z, qqq_nom_z - 0.00828, qqq_nom_z, qqq_nom_z};
|
||
|
|
const double qqq_phi[4] = {5.49779, 0.785398, 2.35619, 3.92699};
|
||
|
|
const double ring_rho[12] = {0.0890601, 0.0889871, 0.0890354, 0.0890247, 0.0890354, 0.0890354, 0.0890247, 0.0890354, 0.0890354, 0.0890247, 0.0890354, 0.0890354};
|
||
|
|
const double ring_phi[12] = {0.785795, 0.262014, 6.02132, 5.49779, 4.97426, 4.45052, 3.92699, 3.40346, 2.87972, 2.35619, 1.83266, 1.30893};
|
||
|
|
/*************************/
|
||
|
|
|
||
|
|
static constexpr double threshold = 0.2; //MeV
|
||
|
|
static constexpr double epsilon = 1.0e-6;
|
||
|
|
static constexpr double deg2rad = M_PI/180.0;
|
||
|
|
};
|
||
|
|
|
||
|
|
#endif
|