<commit_before>#ifndef HBT_correlation_h
#define HBT_correlation_h

#include "ParameterReader.h"
#include "particleSamples.h"

using namespace std;

struct particle_pair
{
    double K_perp, K_phi, K_rap;
    double q_out, q_side, q_long;
    double cos_qx;
};

class HBT_correlation
{
    private:
        ParameterReader *paraRdr;
        string path;
        particleSamples *particle_list;

        int qnpts;
        double q_min, q_max, delta_q;

        int azimuthal_flag;
        int invariant_radius_flag;
        double psi_ref;
        int n_KT, n_Kphi;
        double dKT, dKphi;
        double KT_min, KT_max;
        double Krap_min, Krap_max;
        double buffer_rapidity;
        double *KT_array, *Kphi_array;

        double *q_out, *q_side, *q_long;
        int number_of_mixed_events;
        int number_of_oversample_events;
        unsigned long long int needed_number_of_pairs;
        unsigned long long int number_pairs_num, number_pairs_denorm;
        unsigned long long int *number_of_pairs_numerator_KTdiff;
        unsigned long long int *number_of_pairs_denormenator_KTdiff;
        unsigned long long int **number_of_pairs_numerator_KTKphidiff;
        unsigned long long int **number_of_pairs_denormenator_KTKphidiff;

        // arrays for the invariant radius
        double **q_inv_mean;
        double **correl_1d_inv_num, **correl_1d_inv_denorm;
        double **correl_1d_inv_num_count;

        double ****q_out_mean, ****q_side_mean, ****q_long_mean;
        double ****correl_3d_num, ****correl_3d_denorm;
        double ****correl_3d_num_count;
        double *****q_out_diff_mean, *****q_side_diff_mean;
        double *****q_long_diff_mean;
        double *****correl_3d_Kphi_diff_num, *****correl_3d_Kphi_diff_denorm;
        double *****correl_3d_Kphi_diff_num_count;

    public:
        HBT_correlation(ParameterReader* paraRdr_in, string path_in, 
                        particleSamples *particle_list_in);
        ~HBT_correlation();

        void calculate_flow_event_plane_angle(int n_order);
        void calculate_HBT_correlation_function();
        void combine_and_bin_particle_pairs(int* event_list);
        void combine_and_bin_particle_pairs_mixed_events(
                                          int event_id, int* mixed_event_list);

        void output_correlation_function_inv();
        void output_correlation_function();
        void output_correlation_function_Kphi_differential();

};

#endif<commit_msg>some clean up<commit_after>#ifndef HBT_correlation_h
#define HBT_correlation_h

#include "ParameterReader.h"
#include "particleSamples.h"

using namespace std;

class HBT_correlation {
 private:
    ParameterReader *paraRdr;
    string path;
    particleSamples *particle_list;

    int qnpts;
    double q_min, q_max, delta_q;

    int azimuthal_flag;
    int invariant_radius_flag;
    double psi_ref;
    int n_KT, n_Kphi;
    double dKT, dKphi;
    double KT_min, KT_max;
    double Krap_min, Krap_max;
    double buffer_rapidity;
    double *KT_array, *Kphi_array;

    double *q_out, *q_side, *q_long;
    int number_of_mixed_events;
    int number_of_oversample_events;
    unsigned long long int needed_number_of_pairs;
    unsigned long long int number_pairs_num, number_pairs_denorm;
    unsigned long long int *number_of_pairs_numerator_KTdiff;
    unsigned long long int *number_of_pairs_denormenator_KTdiff;
    unsigned long long int **number_of_pairs_numerator_KTKphidiff;
    unsigned long long int **number_of_pairs_denormenator_KTKphidiff;

    // arrays for the invariant radius
    double **q_inv_mean;
    double **correl_1d_inv_num, **correl_1d_inv_denorm;
    double **correl_1d_inv_num_count;

    double ****q_out_mean, ****q_side_mean, ****q_long_mean;
    double ****correl_3d_num, ****correl_3d_denorm;
    double ****correl_3d_num_count;
    double *****q_out_diff_mean, *****q_side_diff_mean;
    double *****q_long_diff_mean;
    double *****correl_3d_Kphi_diff_num, *****correl_3d_Kphi_diff_denorm;
    double *****correl_3d_Kphi_diff_num_count;

 public:
    HBT_correlation(ParameterReader* paraRdr_in, string path_in, 
                    particleSamples *particle_list_in);
    ~HBT_correlation();

    void calculate_flow_event_plane_angle(int n_order);
    void calculate_HBT_correlation_function();
    void combine_and_bin_particle_pairs(int* event_list);
    void combine_and_bin_particle_pairs_mixed_events(
                                      int event_id, int* mixed_event_list);

    void output_correlation_function_inv();
    void output_correlation_function();
    void output_correlation_function_Kphi_differential();
};

#endif
{'fraction_non_alphanumeric': 0.0798626019750966, 'fraction_numerical': 0.0038643194504079004, 'mean_word_length': 3.5901477832512314, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'dir': 'git-commits-cleaned', 'id': '2810614', 'n_tokens_mistral': 1650, 'n_tokens_neox': 1544, 'n_words': 309}