de/dff/CgemCosmicRayQA_8h_source.html

#ifndef CGEMCOSMICRAYQA

#define CGEMCOSMICRAYQA


#include "CgemRawEvent/CgemDigi.h"

#include "CgemGeomSvc/CgemGeomSvc.h"

#include "CgemGeomSvc/CgemGeoReadoutPlane.h"

#include "CgemGeomSvc/CgemMidDriftPlane.h"

#include "CgemGeomSvc/CgemGeoAlign.h"

#include "ReadCosmicRayData/CgemLUTReader.h"


/* #include "RealizationSvc/RealizationSvc.h"

#include "RealizationSvc/IRealizationSvc.h"

#include <vector>

#include <string>

#include "TFile.h"

#include "TTree.h"

#include "TString.h"

#include "Identifier/Identifier.h"

#include "TH1F.h"

*/

#include "TFile.h"

#include "TTree.h"

#include "TH1.h"

#include "TH1F.h"

#include "TH1D.h"

#include "TF1.h"

#include "TH2.h"

#include "TH2F.h"

#include "TH2D.h"

#include "TArc.h"

#include "TLine.h"

#include "TCanvas.h"

#include "TGraph2D.h"

#include <Math/Functor.h>

#include <TPolyLine3D.h>

#include <Math/Vector3D.h>

#include <Fit/Fitter.h>

#include "TMath.h"

#include "TCanvas.h"

#include "TGraphErrors.h"


class IDataProviderSvc;


#define MAXNOFTRACKS   5

#define MAXNOFCLUSTERS 500 // HACK THIS 10000


#define MAXNOFHITS     2400

#define MAXNOFFITPOINT 6

#define MAXNOFLAYER    3

#define MAXNOFSHEET    2

// stripID

#define MAXNOFSTRIP_L1_x 856

#define MAXNOFSTRIP_L1_v 1174

#define MAXNOFSTRIP_L2_x 630

#define MAXNOFSTRIP_L2_v 1077

#define MAXNOFSTRIP_L3_x 832

#define MAXNOFSTRIP_L3_v 1395

// strip max length

#define MAXLENGTH_L1_x 532

#define MAXLENGTH_L2_x 690

#define MAXLENGTH_L3_x 847

#define MAXLENGTH_L1_v 532/TMath::Cos(TMath::DegToRad()*46.6877)

#define MAXLENGTH_L2_v 690/TMath::Cos(TMath::DegToRad()*31.0337)

#define MAXLENGTH_L3_v 847/TMath::Cos(TMath::DegToRad()*32.9244)


#define MIN_PHI -180

#define MIN_Z -350

#define MIN_V 0

#define CHI_MIN 0

#define ANG_MIN 0


#define MAX_PHI 180

#define MAX_Z 350

#define MAX_V 856

#define CHI_MAX 50

#define ANG_MAX 90


#define MAX_PHI_BIN 1440//360 //720 //120

#define MAX_Z_BIN  350//700 //100

#define MAX_V_BIN 856

#define CHI_BIN 50

#define ANG_BIN 90 // 27


#define PHI_STEP 0.25//1// 0.5 // 360/MAX_PHI_BIN = 360/120

#define Z_STEP 2//1 //2*MAX_Z/MAX_Z_BIN = 1000/100

#define V_STEP 1

#define MIN_ENTRIES 30


#define maxclusize 10

// -----------------------------------------------------------------------


class CgemCosmicRayQA : public Algorithm{

public:

     // Constructor

     CgemCosmicRayQA(const std::string& name, ISvcLocator* pSvcLocator);

     ~CgemCosmicRayQA();

     // Executable

     StatusCode initialize();

     StatusCode execute   ();

     StatusCode finalize  ();


     bool read_file(TString name);

     void define_output();

     void define_hit_histo();

     void define_cluster1d_histo();

     void define_cluster2d_histo();

     void define_track_histo();

     //

     void define_hit_in_time_histo();

     void define_cluster_selected_histo();

     void define_cluster1d_from_fit_histo();


     void fill_hit_histo();

     void fill_cluster1d_histo();

     void fill_cluster2d_histo();

     void fill_track_histo();

     //

     void fill_hit_in_time_histo();

     void fill_cluster_selected_histo();


     int find_closest_exp_cluster2d(double fphi, double fz);

     int find_closest_exp_cluster1dx(double fphi);

     int find_closest_exp_cluster1dv(double fphi, double fz);

     bool apply_cuts();


     void analyze();

     void histo_cosmetics();


     //

     bool plot_debug_tpc(int id_2D);

     void calibration_tpc(int id_2D);

     void tpc_finalize();


private:


     // A data server

     IDataProviderSvc* m_evtSvc;


     CgemLUTReader *lutreader;

     string lutfile;


     CgemGeoReadoutPlane *anode;

     double anode_mid_gap_L1;

     double anode_mid_gap_L2;


     ICgemGeomSvc *m_SvcCgem;

     CgemMidDriftPlane* middleplane;

     CgemGeoAlign* alignment;


     double minDigiTime ;

     double maxDigiTime;

     double select_size;

     double cut_chi2;

     double cut_ene_L1_x;

     double cut_ene_L1_v;

     double cut_ene_L2_x;

     double cut_ene_L2_v;

     double cut_ene_L3_x;

     double cut_ene_L3_v;

     int cut_size_L1_x;

     int cut_size_L1_v;

     int cut_size_L2_x;

     int cut_size_L2_v;

     int cut_size_L3_x;

     int cut_size_L3_v;

     double cut_phi_min_L1;

     double cut_phi_min_L2;

     double cut_phi_min_L3;

     double cut_phi_max_L1;

     double cut_phi_max_L2;

     double cut_phi_max_L3;

     double cut_z_min_L1;

     double cut_z_min_L2;

     double cut_z_min_L3;

     double cut_z_max_L1;

     double cut_z_max_L2;

     double cut_z_max_L3;


     bool align_flag;

     int NSIGMA_INSIDE;


     double anode_mid_gap_radius[3];


     string filename;

     TFile *input;

     TFile *output;

     TTree *tree;


     // hit

     int event;

     int nhit;

     int hit_strip[MAXNOFHITS];

     int hit_view[MAXNOFHITS];

     int hit_layer[MAXNOFHITS];

     int hit_sheet[MAXNOFHITS];

     double hit_length[MAXNOFHITS];

     int hit_channel[MAXNOFHITS];

     int hit_roc[MAXNOFHITS];

     int hit_feb[MAXNOFHITS];

     int hit_tiger[MAXNOFHITS];

     int hit_chip[MAXNOFHITS];

     double hit_t[MAXNOFHITS];

     double hit_q[MAXNOFHITS];

     double hit_x_tpc[MAXNOFHITS];

     double hit_z_tpc[MAXNOFHITS];

     double hit_ex_tpc[MAXNOFHITS];

     double hit_ez_tpc[MAXNOFHITS];

     double hit_deltax_tpc[MAXNOFHITS];

     double hit_deltaz_tpc[MAXNOFHITS];

     int hit_order[MAXNOFHITS];

     int hit_quality[MAXNOFHITS];

     int hit_selgooddigi[MAXNOFHITS];


     // cluster

     int ncluster;

     int ncluster_1d;

     int ncluster_2d;

     int ncluster_2d_L1_S1;

     int ncluster_2d_L2_S1;

     int ncluster_2d_L2_S2;

     int ncluster_2d_L3_S1;

     int ncluster_2d_L3_S2;

     int ncluster_1d_L1_S1_x;

     int ncluster_1d_L2_S1_x;

     int ncluster_1d_L2_S2_x;

     int ncluster_1d_L3_S1_x;

     int ncluster_1d_L3_S2_x;

     int ncluster_1d_L1_S1_v;

     int ncluster_1d_L2_S1_v;

     int ncluster_1d_L2_S2_v;

     int ncluster_1d_L3_S1_v;

     int ncluster_1d_L3_S2_v;


     // 1D cluster

     double cluster_1d_t[MAXNOFCLUSTERS];

     double cluster_1d_q[MAXNOFCLUSTERS];

     double cluster_1d_r[MAXNOFCLUSTERS];

     double cluster_1d_phi[MAXNOFCLUSTERS];

     double cluster_1d_v[MAXNOFCLUSTERS];

     double cluster_1d_v_cc[MAXNOFCLUSTERS];

     double cluster_1d_phi_cc[MAXNOFCLUSTERS];

     double cluster_1d_v_tpc[MAXNOFCLUSTERS];

     double cluster_1d_phi_tpc[MAXNOFCLUSTERS];

     double cluster_1d_slope_tpc[MAXNOFCLUSTERS];

     double cluster_1d_inter_tpc[MAXNOFCLUSTERS];

     double cluster_1d_size_tpc[MAXNOFCLUSTERS];

     int cluster_1d_layerid[MAXNOFCLUSTERS];

     int cluster_1d_sheetid[MAXNOFCLUSTERS];

     int cluster_1d_view[MAXNOFCLUSTERS];

     int cluster_1d_ID[MAXNOFCLUSTERS];

     int cluster_1d_strip1[MAXNOFCLUSTERS];

     int cluster_1d_strip2[MAXNOFCLUSTERS];

     int cluster_1d_size[MAXNOFCLUSTERS];

     int cluster_1d_hitindex[MAXNOFCLUSTERS][50];


     // 2D cluster

     double cluster_2d_t[MAXNOFCLUSTERS];

     double cluster_2d_q[MAXNOFCLUSTERS];

     double cluster_2d_r[MAXNOFCLUSTERS];

     double cluster_2d_phi[MAXNOFCLUSTERS];

     double cluster_2d_z[MAXNOFCLUSTERS];

     double cluster_2d_z_cc[MAXNOFCLUSTERS];

     double cluster_2d_phi_cc[MAXNOFCLUSTERS];

     double cluster_2d_z_tpc[MAXNOFCLUSTERS];

     double cluster_2d_phi_tpc[MAXNOFCLUSTERS];

     int cluster_2d_layerid[MAXNOFCLUSTERS];

     int cluster_2d_sheetid[MAXNOFCLUSTERS];

     int cluster_2d_view[MAXNOFCLUSTERS];

     int cluster_2d_ID[MAXNOFCLUSTERS];

     int cluster_2d_idx[MAXNOFCLUSTERS];

     int cluster_2d_idv[MAXNOFCLUSTERS];

     int cluster_2d_highest[MAXNOFCLUSTERS];


     // track

     double track_dr;

     double track_phi0;

     double track_dz;

     double track_tanL;

     double track_chi2;

     int track_nfitpoint;

     int track_clusterid[MAXNOFFITPOINT];

     int track_layerid[MAXNOFFITPOINT];

     int track_sheetid[MAXNOFFITPOINT];

     int track_test_layerid;

     int track_test_sheetid;

     double track_xpoca_glo;

     double track_ypoca_glo;

     double track_zpoca_glo;

     double track_x1top_glo;

     double track_y1top_glo;

     double track_z1top_glo;

     double track_phi1top_loc;

     double track_v1top_loc;

     double track_x1bot_glo;

     double track_y1bot_glo;

     double track_z1bot_glo;

     double track_phi1bot_loc;

     double track_v1bot_loc;


     double track_x2top_glo;

     double track_y2top_glo;

     double track_z2top_glo;

     double track_phi2top_loc;

     double track_v2top_loc;

     double track_x2bot_glo;

     double track_y2bot_glo;

     double track_z2bot_glo;

     double track_phi2bot_loc;

     double track_v2bot_loc;


     double track_x3top_glo;

     double track_y3top_glo;

     double track_z3top_glo;

     double track_phi3top_loc;

     double track_v3top_loc;

     double track_x3bot_glo;

     double track_y3bot_glo;

     double track_z3bot_glo;

     double track_phi3bot_loc;

     double track_v3bot_loc;


     // Intersection angles

     // L1

     double ang_xy_L1;

     double ang_yz_L1;

     // L2

     double ang_xy_L2;

     double ang_yz_L2;

     //L3

     double ang_xy_L3;

     double ang_yz_L3;


     // histograms ---------------------------------------

     // HIT

     // nof hits

     // x

     TH1F *h_nofhit_L1_S1_x;

     TH1F *h_nofhit_L1_S2_x;

     TH1F *h_nofhit_L2_S1_x;

     TH1F *h_nofhit_L2_S2_x;

     TH1F *h_nofhit_L3_S1_x;

     TH1F *h_nofhit_L3_S2_x;

     // v

     TH1F *h_nofhit_L1_S1_v;

     TH1F *h_nofhit_L2_S1_v;

     TH1F *h_nofhit_L2_S2_v;

     TH1F *h_nofhit_L3_S1_v;

     TH1F *h_nofhit_L3_S2_v;


     // charge

     // x

     TH1F *h_hit_charge_L1_S1_x;

     TH1F *h_hit_charge_L1_S2_x;

     TH1F *h_hit_charge_L2_S1_x;

     TH1F *h_hit_charge_L2_S2_x;

     TH1F *h_hit_charge_L3_S1_x;

     TH1F *h_hit_charge_L3_S2_x;

     // v

     TH1F *h_hit_charge_L1_S1_v;

     TH1F *h_hit_charge_L2_S1_v;

     TH1F *h_hit_charge_L2_S2_v;

     TH1F *h_hit_charge_L3_S1_v;

     TH1F *h_hit_charge_L3_S2_v;

     // time

     // x

     TH1F *h_hit_time_L1_S1_x;

     TH1F *h_hit_time_L1_S2_x;

     TH1F *h_hit_time_L2_S1_x;

     TH1F *h_hit_time_L2_S2_x;

     TH1F *h_hit_time_L3_S1_x;

     TH1F *h_hit_time_L3_S2_x;

     // v

     TH1F *h_hit_time_L1_S1_v;

     TH1F *h_hit_time_L2_S1_v;

     TH1F *h_hit_time_L2_S2_v;

     TH1F *h_hit_time_L3_S1_v;

     TH1F *h_hit_time_L3_S2_v;

     // charge vs strip ID

     // x

     TH2F *h_hit_charge_vs_strip_L1_S1_x;

     TH2F *h_hit_charge_vs_strip_L1_S2_x;

     TH2F *h_hit_charge_vs_strip_L2_S1_x;

     TH2F *h_hit_charge_vs_strip_L2_S2_x;

     TH2F *h_hit_charge_vs_strip_L3_S1_x;

     TH2F *h_hit_charge_vs_strip_L3_S2_x;

     // v

     TH2F *h_hit_charge_vs_strip_L1_S1_v;

     TH2F *h_hit_charge_vs_strip_L2_S1_v;

     TH2F *h_hit_charge_vs_strip_L2_S2_v;

     TH2F *h_hit_charge_vs_strip_L3_S1_v;

     TH2F *h_hit_charge_vs_strip_L3_S2_v;

     // time vs strip ID

     // x

     TH2F *h_hit_time_vs_strip_L1_S1_x;

     TH2F *h_hit_time_vs_strip_L1_S2_x;

     TH2F *h_hit_time_vs_strip_L2_S1_x;

     TH2F *h_hit_time_vs_strip_L2_S2_x;

     TH2F *h_hit_time_vs_strip_L3_S1_x;

     TH2F *h_hit_time_vs_strip_L3_S2_x;

     // v

     TH2F *h_hit_time_vs_strip_L1_S1_v;

     TH2F *h_hit_time_vs_strip_L2_S1_v;

     TH2F *h_hit_time_vs_strip_L2_S2_v;

     TH2F *h_hit_time_vs_strip_L3_S1_v;

     TH2F *h_hit_time_vs_strip_L3_S2_v;

     // charge vs time

     // x

     TH2F *h_hit_charge_vs_time_L1_S1_x;

     TH2F *h_hit_charge_vs_time_L1_S2_x;

     TH2F *h_hit_charge_vs_time_L2_S1_x;

     TH2F *h_hit_charge_vs_time_L2_S2_x;

     TH2F *h_hit_charge_vs_time_L3_S1_x;

     TH2F *h_hit_charge_vs_time_L3_S2_x;

     // v

     TH2F *h_hit_charge_vs_time_L1_S1_v;

     TH2F *h_hit_charge_vs_time_L2_S1_v;

     TH2F *h_hit_charge_vs_time_L2_S2_v;

     TH2F *h_hit_charge_vs_time_L3_S1_v;

     TH2F *h_hit_charge_vs_time_L3_S2_v;

     // charge vs length

     // v

     TH2F *h_hit_charge_vs_length_L1_S1_v;

     TH2F *h_hit_charge_vs_length_L2_S1_v;

     TH2F *h_hit_charge_vs_length_L2_S2_v;

     TH2F *h_hit_charge_vs_length_L3_S1_v;

     TH2F *h_hit_charge_vs_length_L3_S2_v;

     // HIT IN TIME WINDOW

     TH1F *h_hit_in_time_charge_L1_S1_x;

     TH1F *h_hit_in_time_charge_L1_S2_x;

     TH1F *h_hit_in_time_charge_L2_S1_x;

     TH1F *h_hit_in_time_charge_L2_S2_x;

     TH1F *h_hit_in_time_charge_L3_S1_x;

     TH1F *h_hit_in_time_charge_L3_S2_x;

     // v

     TH1F *h_hit_in_time_charge_L1_S1_v;

     TH1F *h_hit_in_time_charge_L2_S1_v;

     TH1F *h_hit_in_time_charge_L2_S2_v;

     TH1F *h_hit_in_time_charge_L3_S1_v;

     TH1F *h_hit_in_time_charge_L3_S2_v;

     // count on each strip

     // x

     TH1F *h_hit_in_time_strip_L1_S1_x;

     TH1F *h_hit_in_time_strip_L1_S2_x;

     TH1F *h_hit_in_time_strip_L2_S1_x;

     TH1F *h_hit_in_time_strip_L2_S2_x;

     TH1F *h_hit_in_time_strip_L3_S1_x;

     TH1F *h_hit_in_time_strip_L3_S2_x;

     // v

     TH1F *h_hit_in_time_strip_L1_S1_v;

     TH1F *h_hit_in_time_strip_L2_S1_v;

     TH1F *h_hit_in_time_strip_L2_S2_v;

     TH1F *h_hit_in_time_strip_L3_S1_v;

     TH1F *h_hit_in_time_strip_L3_S2_v;

     // charge vs strip ID

     // x

     TH2F *h_hit_in_time_charge_vs_strip_L1_S1_x;

     TH2F *h_hit_in_time_charge_vs_strip_L1_S2_x;

     TH2F *h_hit_in_time_charge_vs_strip_L2_S1_x;

     TH2F *h_hit_in_time_charge_vs_strip_L2_S2_x;

     TH2F *h_hit_in_time_charge_vs_strip_L3_S1_x;

     TH2F *h_hit_in_time_charge_vs_strip_L3_S2_x;

     // v

     TH2F *h_hit_in_time_charge_vs_strip_L1_S1_v;

     TH2F *h_hit_in_time_charge_vs_strip_L2_S1_v;

     TH2F *h_hit_in_time_charge_vs_strip_L2_S2_v;

     TH2F *h_hit_in_time_charge_vs_strip_L3_S1_v;

     TH2F *h_hit_in_time_charge_vs_strip_L3_S2_v;

     // CLUSTER 1d

     // nof cluster 1d

     // x

     TH1F *h_nofcluster1d_L1_S1_x;

     TH1F *h_nofcluster1d_L1_S2_x;

     TH1F *h_nofcluster1d_L2_S1_x;

     TH1F *h_nofcluster1d_L2_S2_x;

     TH1F *h_nofcluster1d_L3_S1_x;

     TH1F *h_nofcluster1d_L3_S2_x;

     // v

     TH1F *h_nofcluster1d_L1_S1_v;

     TH1F *h_nofcluster1d_L2_S1_v;

     TH1F *h_nofcluster1d_L2_S2_v;

     TH1F *h_nofcluster1d_L3_S1_v;

     TH1F *h_nofcluster1d_L3_S2_v;

     // cl. size

     // x

     TH1F *h_cluster1d_size_L1_S1_x;

     TH1F *h_cluster1d_size_L1_S2_x;

     TH1F *h_cluster1d_size_L2_S1_x;

     TH1F *h_cluster1d_size_L2_S2_x;

     TH1F *h_cluster1d_size_L3_S1_x;

     TH1F *h_cluster1d_size_L3_S2_x;

     // v

     TH1F *h_cluster1d_size_L1_S1_v;

     TH1F *h_cluster1d_size_L2_S1_v;

     TH1F *h_cluster1d_size_L2_S2_v;

     TH1F *h_cluster1d_size_L3_S1_v;

     TH1F *h_cluster1d_size_L3_S2_v;

     // charge vs phi (x view)

     TH2F *h_cluster1d_charge_vs_phi_L1_S1_x;

     TH2F *h_cluster1d_charge_vs_phi_L1_S2_x;

     TH2F *h_cluster1d_charge_vs_phi_L2_S1_x;

     TH2F *h_cluster1d_charge_vs_phi_L2_S2_x;

     TH2F *h_cluster1d_charge_vs_phi_L3_S1_x;

     TH2F *h_cluster1d_charge_vs_phi_L3_S2_x;

     // cherge vs v (v view)

     TH2F *h_cluster1d_charge_vs_v_L1_S1_v;

     TH2F *h_cluster1d_charge_vs_v_L2_S1_v;

     TH2F *h_cluster1d_charge_vs_v_L2_S2_v;

     TH2F *h_cluster1d_charge_vs_v_L3_S1_v;

     TH2F *h_cluster1d_charge_vs_v_L3_S2_v;

     // charge of cluster used in the fit (x view)

     TH1F *h_cluster1d_charge_fitted_L1_S1_x;

     TH1F *h_cluster1d_charge_fitted_L1_S2_x;

     TH1F *h_cluster1d_charge_fitted_L2_S1_x;

     TH1F *h_cluster1d_charge_fitted_L2_S2_x;

     TH1F *h_cluster1d_charge_fitted_L3_S1_x;

     TH1F *h_cluster1d_charge_fitted_L3_S2_x;

     // charge of cluster used in the fit (v view)

     TH1F *h_cluster1d_charge_fitted_L1_S1_v;

     TH1F *h_cluster1d_charge_fitted_L1_S2_v;

     TH1F *h_cluster1d_charge_fitted_L2_S1_v;

     TH1F *h_cluster1d_charge_fitted_L2_S2_v;

     TH1F *h_cluster1d_charge_fitted_L3_S1_v;

     TH1F *h_cluster1d_charge_fitted_L3_S2_v;

     // CLUSTER 2d

     // nof cluster 2d

     TH1F *h_nofcluster2d_L1_S1;

     TH1F *h_nofcluster2d_L1_S2;

     TH1F *h_nofcluster2d_L2_S1;

     TH1F *h_nofcluster2d_L2_S2;

     TH1F *h_nofcluster2d_L3_S1;

     TH1F *h_nofcluster2d_L3_S2;

     // total cl. size

     TH1F *h_cluster2d_size_L1_S1;

     TH1F *h_cluster2d_size_L1_S2;

     TH1F *h_cluster2d_size_L2_S1;

     TH1F *h_cluster2d_size_L2_S2;

     TH1F *h_cluster2d_size_L3_S1;

     TH1F *h_cluster2d_size_L3_S2;

     // charge vs phi

     TH2F *h_cluster2d_charge_vs_phi_L1_S1;

     TH2F *h_cluster2d_charge_vs_phi_L1_S2;

     TH2F *h_cluster2d_charge_vs_phi_L2_S1;

     TH2F *h_cluster2d_charge_vs_phi_L2_S2;

     TH2F *h_cluster2d_charge_vs_phi_L3_S1;

     TH2F *h_cluster2d_charge_vs_phi_L3_S2;

     // cherge vs z

     TH2F *h_cluster2d_charge_vs_z_L1_S1;

     TH2F *h_cluster2d_charge_vs_z_L1_S2;

     TH2F *h_cluster2d_charge_vs_z_L2_S1;

     TH2F *h_cluster2d_charge_vs_z_L2_S2;

     TH2F *h_cluster2d_charge_vs_z_L3_S1;

     TH2F *h_cluster2d_charge_vs_z_L3_S2;

     // z vs phi

     TH2F *h_cluster2d_z_vs_phi_L1_S1;

     TH2F *h_cluster2d_z_vs_phi_L1_S2;

     TH2F *h_cluster2d_z_vs_phi_L2_S1;

     TH2F *h_cluster2d_z_vs_phi_L2_S2;

     TH2F *h_cluster2d_z_vs_phi_L3_S1;

     TH2F *h_cluster2d_z_vs_phi_L3_S2;

     // CLUSTER 1d SELECTED

     // charge

     TH1F *h_cluster1d_charge_selected_L1_S1_x;

     TH1F *h_cluster1d_charge_selected_L1_S2_x;

     TH1F *h_cluster1d_charge_selected_L2_S1_x;

     TH1F *h_cluster1d_charge_selected_L2_S2_x;

     TH1F *h_cluster1d_charge_selected_L3_S1_x;

     TH1F *h_cluster1d_charge_selected_L3_S2_x;

     TH1F *h_cluster1d_charge_selected_L1_S1_v;

     TH1F *h_cluster1d_charge_selected_L1_S2_v;

     TH1F *h_cluster1d_charge_selected_L2_S1_v;

     TH1F *h_cluster1d_charge_selected_L2_S2_v;

     TH1F *h_cluster1d_charge_selected_L3_S1_v;

     TH1F *h_cluster1d_charge_selected_L3_S2_v;

     // cl. size

     TH1F *h_cluster1d_size_selected_L1_S1_x;

     TH1F *h_cluster1d_size_selected_L1_S2_x;

     TH1F *h_cluster1d_size_selected_L2_S1_x;

     TH1F *h_cluster1d_size_selected_L2_S2_x;

     TH1F *h_cluster1d_size_selected_L3_S1_x;

     TH1F *h_cluster1d_size_selected_L3_S2_x;

     TH1F *h_cluster1d_size_selected_L1_S1_v;

     TH1F *h_cluster1d_size_selected_L1_S2_v;

     TH1F *h_cluster1d_size_selected_L2_S1_v;

     TH1F *h_cluster1d_size_selected_L2_S2_v;

     TH1F *h_cluster1d_size_selected_L3_S1_v;

     TH1F *h_cluster1d_size_selected_L3_S2_v;

     //Charge maxQ hit in lcuster

     TH1F *h_cl_hitmaxQ_charge_selected_L1_S1_x;

     TH1F *h_cl_hitmaxQ_charge_selected_L1_S2_x;

     TH1F *h_cl_hitmaxQ_charge_selected_L2_S1_x;

     TH1F *h_cl_hitmaxQ_charge_selected_L2_S2_x;

     TH1F *h_cl_hitmaxQ_charge_selected_L3_S1_x;

     TH1F *h_cl_hitmaxQ_charge_selected_L3_S2_x;

     TH1F *h_cl_hitmaxQ_charge_selected_L1_S1_v;

     TH1F *h_cl_hitmaxQ_charge_selected_L1_S2_v;

     TH1F *h_cl_hitmaxQ_charge_selected_L2_S1_v;

     TH1F *h_cl_hitmaxQ_charge_selected_L2_S2_v;

     TH1F *h_cl_hitmaxQ_charge_selected_L3_S1_v;

     TH1F *h_cl_hitmaxQ_charge_selected_L3_S2_v;

     // CLUSTER 2d SELECTED

     TH1F *h_cluster2d_charge_selected_L1_S1;

     TH1F *h_cluster2d_charge_selected_L1_S2;

     TH1F *h_cluster2d_charge_selected_L2_S1;

     TH1F *h_cluster2d_charge_selected_L2_S2;

     TH1F *h_cluster2d_charge_selected_L3_S1;

     TH1F *h_cluster2d_charge_selected_L3_S2;

     // TRACK

     // nof fitted tracks

     TH1F *h_noftrack;


     // chi2 distro

     TH1F *h_track_chi2;


     // xpoca, ypoca, zpoca distro

     TH1F *h_track_xpoca;

     TH1F *h_track_ypoca;

     TH1F *h_track_zpoca;


     // TEST PLANE

     // residual in phi

     TH1F *h_test_residual_rphi;

     TH1F *h_test_residual_rphi_tpc;


     // residual in phi vs chi2

     TH2F *h_test_res_rphi_vs_chi2;

     TH2F *h_test_res_rphi_vs_chi2_tpc;


     // residual in z

     TH1F *h_test_residual_z;

     TH1F *h_test_residual_z_tpc;


     // residual in z vs chi2

     TH2F *h_test_res_z_vs_chi2;

     TH2F *h_test_res_z_vs_chi2_tpc;


     // TRACKER

     // residual in phi

     TH1F *h_residual_rphi_L1_S1;

     TH1F *h_residual_rphi_L1_S2;

     TH1F *h_residual_rphi_L2_S1;

     TH1F *h_residual_rphi_L2_S2;

     TH1F *h_residual_rphi_L3_S1;

     TH1F *h_residual_rphi_L3_S2;


     TH2F *h_res_rphi_vs_chi2_L1_S1;

     TH2F *h_res_rphi_vs_chi2_L1_S2;

     TH2F *h_res_rphi_vs_chi2_L2_S1;

     TH2F *h_res_rphi_vs_chi2_L2_S2;

     TH2F *h_res_rphi_vs_chi2_L3_S1;

     TH2F *h_res_rphi_vs_chi2_L3_S2;


     // TPC residual in phi

     TH1F *h_residual_rphi_L1_S1_tpc;

     TH1F *h_residual_rphi_L1_S2_tpc;

     TH1F *h_residual_rphi_L2_S1_tpc;

     TH1F *h_residual_rphi_L2_S2_tpc;

     TH1F *h_residual_rphi_L3_S1_tpc;

     TH1F *h_residual_rphi_L3_S2_tpc;


     TH2F *h_res_rphi_vs_chi2_L1_S1_tpc;

     TH2F *h_res_rphi_vs_chi2_L1_S2_tpc;

     TH2F *h_res_rphi_vs_chi2_L2_S1_tpc;

     TH2F *h_res_rphi_vs_chi2_L2_S2_tpc;

     TH2F *h_res_rphi_vs_chi2_L3_S1_tpc;

     TH2F *h_res_rphi_vs_chi2_L3_S2_tpc;


     // residual in z

     TH1F *h_residual_z_L1_S1;

     TH1F *h_residual_z_L1_S2;

     TH1F *h_residual_z_L2_S1;

     TH1F *h_residual_z_L2_S2;

     TH1F *h_residual_z_L3_S1;

     TH1F *h_residual_z_L3_S2;


     TH2F *h_res_z_vs_chi2_L1_S1;

     TH2F *h_res_z_vs_chi2_L1_S2;

     TH2F *h_res_z_vs_chi2_L2_S1;

     TH2F *h_res_z_vs_chi2_L2_S2;

     TH2F *h_res_z_vs_chi2_L3_S1;

     TH2F *h_res_z_vs_chi2_L3_S2;


     // TPC residual in z

     TH1F *h_residual_z_L1_S1_tpc;

     TH1F *h_residual_z_L1_S2_tpc;

     TH1F *h_residual_z_L2_S1_tpc;

     TH1F *h_residual_z_L2_S2_tpc;

     TH1F *h_residual_z_L3_S1_tpc;

     TH1F *h_residual_z_L3_S2_tpc;


     TH2F *h_res_z_vs_chi2_L1_S1_tpc;

     TH2F *h_res_z_vs_chi2_L1_S2_tpc;

     TH2F *h_res_z_vs_chi2_L2_S1_tpc;

     TH2F *h_res_z_vs_chi2_L2_S2_tpc;

     TH2F *h_res_z_vs_chi2_L3_S1_tpc;

     TH2F *h_res_z_vs_chi2_L3_S2_tpc;


     // signal

     TH1F *h_signal_charge_2d;

     TH1F *h_signal_charge_1d_x;

     TH1F *h_signal_charge_1d_v;

     TH1F *h_signal_size_1d_x;

     TH1F *h_signal_size_1d_v;


     TH1F *h_signal_charge_2d_chi2cut;

     TH1F *h_signal_charge_1d_x_chi2cut;

     TH1F *h_signal_charge_1d_v_chi2cut;

     TH1F *h_signal_size_1d_x_chi2cut;

     TH1F *h_signal_size_1d_v_chi2cut;


     // background

     TH1F *h_background_charge_2d;

     TH1F *h_background_charge_1d_x;

     TH1F *h_background_charge_1d_v;

     TH1F *h_background_size_1d_x;

     TH1F *h_background_size_1d_v;


     // efficiencies

     TH1F *h_eff_vs_phi;

     TH1F *h_eff_vs_z;


     TH1F *h_eff_vs_phi_1d;

     TH1F *h_eff_vs_v_1d;


     TH2F *h_eff_vs_phi_vs_z;


     TH2F *h_eff_vs_phi_vs_chi2;

     TH2F *h_eff_vs_z_vs_chi2;


     TH1F *h_eff_vs_FEB;


     TH1F *h_track_selection_vs_phi;

     TH1F *h_track_selection_vs_z;

     TH2F *h_track_selection_vs_phi_vs_z;


     // Resolution vs Chi2 Plots

     TH1F *h_resolution_rphi_vs_chi2;

     TH1F *h_resolution_rphi_tpc_vs_chi2;


     TH1F *h_resolution_z_vs_chi2;

     TH1F *h_resolution_z_tpc_vs_chi2;


     // Incident Angles

     TH2F *h_test_residual_rphi_vs_angxy_L1_cc;

     TH2F *h_test_residual_rphi_vs_angxy_L1_tpc;

     TH2F *h_test_residual_rphi_vs_angxy_L2_cc;

     TH2F *h_test_residual_rphi_vs_angxy_L2_tpc;

     TH2F *h_test_residual_rphi_vs_angxy_L3_cc;

     TH2F *h_test_residual_rphi_vs_angxy_L3_tpc;


     TH1F *h_resolution_vs_ang_xy_L1;

     TH1F *h_resolution_vs_ang_xy_L2;

     TH1F *h_resolution_vs_ang_xy_L3;

     TH1F *h_resolution_tpc_vs_ang_xy_L1;

     TH1F *h_resolution_tpc_vs_ang_xy_L2;

     TH1F *h_resolution_tpc_vs_ang_xy_L3;


     //TPC debug

     string output_path_tpc;

     bool bool_plot_tpc;

     TH2F *h_map_xy_event_trk;

     TH2F *h_map_xz_event_trk;

     TH2F *h_map_xy_2d;

     TH2F *h_map_xz_2d;

     TH2F *h_map_xy_1d_tpc;

     TH2F *h_map_xz_1d_tpc;

     TH2F *h_map_xy_tpc_eff;

     TH2F *h_map_xz_tpc_eff;

     TH2F *h_map_xy_tpc_fail;

     TH2F *h_map_xz_tpc_fail;

     TH2F *h_map_xy_cc_eff;

     TH2F *h_map_xz_cc_eff;

     TH2F *h_map_xy_cc_tpc_diff;

     TH2F *h_map_xz_cc_tpc_diff;

     TH2F *h_map_xy_cc_tpc_diff2;

     TH2F *h_map_xz_cc_tpc_diff2;

     TH2F *h_ang_trk_meas;

     TH1F *h_ang_diff;

     TH2F *h_ang_diff_vs_trk;

     TH2F *h_map_xy_q2d;

     TH2F *h_map_xz_q2d;

     TH2F *h_map_xy_dT;

     TH2F *h_map_xz_dT;

     vector <float> vector_time_xy[1000][1000];

     vector <float> vector_time_xz[1000][1000];


     //TPC calibration and correction

     TH1F *h_deltaX_order_size[maxclusize][maxclusize];

     TH1F *h_deltaZ_order_size[maxclusize][maxclusize];

     TGraphErrors *deltaX_order_L1[maxclusize];

     TGraphErrors *deltaZ_order_L1[maxclusize];

     TGraphErrors *deltaX_par0, *deltaX_par1, *deltaZ_par0,*deltaZ_par1, *shift0;

     TF1 *f_deltaX[maxclusize];

     TF1 *f_deltaZ[maxclusize];

     float x_order[maxclusize], xx_order[maxclusize], y_deltaX[maxclusize], ey_deltaX[maxclusize], y_deltaZ[maxclusize], ey_deltaZ[maxclusize];

     float y_deltaX_par0[maxclusize], ey_deltaX_par0[maxclusize], y_deltaX_par1[maxclusize], ey_deltaX_par1[maxclusize];

     float y_deltaZ_par0[maxclusize], ey_deltaZ_par0[maxclusize], y_deltaZ_par1[maxclusize], ey_deltaZ_par1[maxclusize];

     float y_shift0[maxclusize], ey_shift0[maxclusize];

     TH1F *h_ang_tpc;


     // for analysis --------------------------

     int nfittedtrack;

     int nchi2;

     int nstopped;

     int nvalidtrack;

     int ninside;

     int ninside1Dx;

     int ninside1Dz;

     int noutside;

     int nentries_residual_rphi;

     double mean_residual_rphi;

     double sigma_residual_rphi;

     int nentries_residual_z;

     double mean_residual_z;

     double sigma_residual_z;


     // For Efficiency

     int n_validtrack_phi_arr[MAX_PHI_BIN];

     int n_insidetrack_phi_arr[MAX_PHI_BIN];


     int n_validtrack_z_arr[MAX_Z_BIN];

     int n_insidetrack_z_arr[MAX_Z_BIN];


     int n_validtrack_v_arr[MAX_V_BIN];

     int n_insidetrack_v_arr[MAX_V_BIN];


     int n_validtrack_phi_chi2_mat[CHI_BIN][MAX_PHI_BIN];

     int n_insidetrack_phi_chi2_mat[CHI_BIN][MAX_PHI_BIN];


     int n_validtrack_phi_z_mat[MAX_PHI_BIN][MAX_Z_BIN];

     int n_insidetrack_phi_z_mat[MAX_PHI_BIN][MAX_Z_BIN];


     int n_validtrack_z_chi2_mat[CHI_BIN][MAX_Z_BIN];

     int n_insidetrack_z_chi2_mat[CHI_BIN][MAX_Z_BIN];


     int n_validtrack_feb_arr[80];

     int n_insidetrack_feb_arr[80];


     //Efficiency 1D

     int n_validtrack_phi_1d_arr[MAX_PHI_BIN];

     int n_insidetrack_phi_1d_arr[MAX_PHI_BIN];

     int n_validtrack_v_1d_arr[MAX_PHI_BIN];

     int n_insidetrack_v_1d_arr[MAX_PHI_BIN];


     double chi2_cut_arr[CHI_BIN];


     int test_l;

     int test_ml;

     int test_mh;

     int test_h;

     int test_p;


     //


     std::vector< double > vector_test_residual_rphi;

     std::vector< double > vector_test_residual_z;

     std::vector< double > vector_test_charge;

     std::vector< double > vector_test_charge_x;

     std::vector< double > vector_test_charge_v;


     std::vector< double > vector_test_phi;

     std::vector< double > vector_test_z;

     std::vector< double > vector_fit_phi;

     std::vector< double > vector_fit_z;

     std::vector< double > vector_fit_v;

     std::vector< double > vector_chi2;


     std::vector< int > vector_test_size_x;

     std::vector< int > vector_test_size_v;

     std::vector< int > vector_test_clusterid;

     std::vector< int > vector_test_idx;

     std::vector< int > vector_test_idv;


     std::vector< double >  vector_test_phi_1d;

     std::vector< double >  vector_fit_phi_1d;

     std::vector< double >  vector_test_residual_rphi_1d;

     std::vector< double >  vector_test_v_1d;

     std::vector< double >  vector_fit_v_1d;

     std::vector< double >  vector_test_residual_v_1d;


     int test_layer_id_all;

     int test_sheet_id_all;


};


#endif

MAXNOFFITPOINT
#define MAXNOFFITPOINT
Definition CgemCosmicRayQA.h:49

CHI_BIN
#define CHI_BIN
Definition CgemCosmicRayQA.h:82

MAX_V_BIN
#define MAX_V_BIN
Definition CgemCosmicRayQA.h:81

MAXNOFCLUSTERS
#define MAXNOFCLUSTERS
Definition CgemCosmicRayQA.h:46

MAX_PHI_BIN
#define MAX_PHI_BIN
Definition CgemCosmicRayQA.h:79

maxclusize
#define maxclusize
Definition CgemCosmicRayQA.h:90

MAX_Z_BIN
#define MAX_Z_BIN
Definition CgemCosmicRayQA.h:80

MAXNOFHITS
#define MAXNOFHITS
Definition CgemCosmicRayQA.h:48

CgemDigi.h

CgemGeoAlign.h

CgemGeoReadoutPlane.h

CgemGeomSvc.h

CgemLUTReader.h

CgemMidDriftPlane.h

CgemCosmicRayQA
Definition CgemCosmicRayQA.h:92

CgemCosmicRayQA::CgemCosmicRayQA
CgemCosmicRayQA(const std::string &name, ISvcLocator *pSvcLocator)
Definition CgemCosmicRayQA.cxx:49

CgemCosmicRayQA::analyze
void analyze()
Definition CgemCosmicRayQA.cxx:2414

CgemCosmicRayQA::tpc_finalize
void tpc_finalize()
Definition CgemCosmicRayQA.cxx:3415

CgemCosmicRayQA::fill_hit_histo
void fill_hit_histo()
Definition CgemCosmicRayQA.cxx:1165

CgemCosmicRayQA::fill_cluster2d_histo
void fill_cluster2d_histo()
Definition CgemCosmicRayQA.cxx:1716

CgemCosmicRayQA::calibration_tpc
void calibration_tpc(int id_2D)
Definition CgemCosmicRayQA.cxx:3344

CgemCosmicRayQA::define_cluster1d_histo
void define_cluster1d_histo()
Definition CgemCosmicRayQA.cxx:438

CgemCosmicRayQA::define_cluster2d_histo
void define_cluster2d_histo()
Definition CgemCosmicRayQA.cxx:598

CgemCosmicRayQA::define_track_histo
void define_track_histo()
Definition CgemCosmicRayQA.cxx:651

CgemCosmicRayQA::define_cluster_selected_histo
void define_cluster_selected_histo()
Definition CgemCosmicRayQA.cxx:497

CgemCosmicRayQA::apply_cuts
bool apply_cuts()
Definition CgemCosmicRayQA.cxx:2197

CgemCosmicRayQA::histo_cosmetics
void histo_cosmetics()
Definition CgemCosmicRayQA.cxx:2876

CgemCosmicRayQA::execute
StatusCode execute()
Definition CgemCosmicRayQA.cxx:1024

CgemCosmicRayQA::define_output
void define_output()
Definition CgemCosmicRayQA.cxx:255

CgemCosmicRayQA::plot_debug_tpc
bool plot_debug_tpc(int id_2D)
Definition CgemCosmicRayQA.cxx:3079

CgemCosmicRayQA::fill_track_histo
void fill_track_histo()
Definition CgemCosmicRayQA.cxx:1792

CgemCosmicRayQA::find_closest_exp_cluster1dx
int find_closest_exp_cluster1dx(double fphi)
Definition CgemCosmicRayQA.cxx:3016

CgemCosmicRayQA::find_closest_exp_cluster2d
int find_closest_exp_cluster2d(double fphi, double fz)
Definition CgemCosmicRayQA.cxx:2988

CgemCosmicRayQA::fill_hit_in_time_histo
void fill_hit_in_time_histo()
Definition CgemCosmicRayQA.cxx:1313

CgemCosmicRayQA::initialize
StatusCode initialize()
Definition CgemCosmicRayQA.cxx:876

CgemCosmicRayQA::define_cluster1d_from_fit_histo
void define_cluster1d_from_fit_histo()
Definition CgemCosmicRayQA.cxx:568

CgemCosmicRayQA::find_closest_exp_cluster1dv
int find_closest_exp_cluster1dv(double fphi, double fz)
Definition CgemCosmicRayQA.cxx:3045

CgemCosmicRayQA::~CgemCosmicRayQA
~CgemCosmicRayQA()
Definition CgemCosmicRayQA.cxx:91

CgemCosmicRayQA::fill_cluster1d_histo
void fill_cluster1d_histo()
Definition CgemCosmicRayQA.cxx:1411

CgemCosmicRayQA::fill_cluster_selected_histo
void fill_cluster_selected_histo()
Definition CgemCosmicRayQA.cxx:1520

CgemCosmicRayQA::define_hit_histo
void define_hit_histo()
Definition CgemCosmicRayQA.cxx:261

CgemCosmicRayQA::finalize
StatusCode finalize()
Definition CgemCosmicRayQA.cxx:1068

CgemCosmicRayQA::read_file
bool read_file(TString name)
Definition CgemCosmicRayQA.cxx:96

CgemCosmicRayQA::define_hit_in_time_histo
void define_hit_in_time_histo()
Definition CgemCosmicRayQA.cxx:375

CgemGeoAlign
Definition CgemGeoAlign.h:61

CgemGeoReadoutPlane
Definition CgemGeoReadoutPlane.h:14

CgemLUTReader
Definition CgemLUTReader.h:16

CgemMidDriftPlane
Definition CgemMidDriftPlane.h:8

ICgemGeomSvc
Definition ICgemGeomSvc.h:22