prova.h

//////////////////////////////////////////////////////////
// This class has been automatically generated on
// Fri Mar 30 12:27:06 2018 by ROOT version 6.12/06
// from TTree apv_raw/APVRawData
// found on file: ../TBReco/data/run177.root
//////////////////////////////////////////////////////////

#ifndef prova_h
#define prova_h

#include <TROOT.h>
#include <TChain.h>
#include <TFile.h>
#include <TH1F.h>
#include <TH2F.h>
#include <string.h>
#include <iostream>
#include <math.h>
#include <TCanvas.h>

// Header file for the classes stored in the TTree if any.
#include <vector>

using namespace std;

class prova {
public:
  TTree            *fChain;   //!pointer to the analyzed TTree or TChain
  Int_t            fCurrent; //!current Tree number in a TChain
  // Fixed size dimensions of array or collections stored in the TTree if any.

  // Declaration of leaf types
  //IMPOSTO ANGOLO E LARGHEZZA IN STRIP DELLA ROAD
  //Double_t      cos(angle_2) = 0.982416299; 
  map<int,int>     infncard;
  string           _ch_p1 = "P1_M01", _ch_p2 = "P2_M01", _ch_p3 = "P3_M01", _ch_p4 = "P4_M01";
  Double_t         _super_eta;
  Int_t            _N_lay = 4;
  Int_t            _N_pcb = 5;
  Int_t            eff_num = 0;
  Int_t            make_fit = 0;
  Double_t         pitch = 0.425; //mm
  //double      temp_centroid[4];
  int              n_volte = 0;
  int              eff_evt = 0;
  int              n_event = 0;
  double           temp_centroid_pcb[20];
  int              n_clu_1, n_clu_2, n_clu_3, n_clu_4;
  int              pcb_n_cluster[20];
  double           num_eff[4];
  double           den_eff[4];
  double           num_fit_eff[4];
  double           den_fit_eff[4];
  double           num_hist_eff[4];
  double           den_hist_eff[4];
  int              n_holes_max = 2;
  Double_t         n_sigma = 4;
  Double_t         angle = 0.0261799, phi = 0.183;
  Double_t         major_base = 1319.2;
  Double_t         road_eta = 2;
  double           sigma_road = 1.;
  //Double_t       road_stereo = 2*(()*tan(angle) + sigma_road);
  map<int,int>     ch1, ch2, ch3, ch4;
  vector<double>   res_1, res_2, res_3, res_4;
  vector<TH1F*>    h_residues_pos;
  vector<TH1F*>    h_residues_pos_fit;
  vector<TH2F*>    h2_extrapolate_pos;
  vector<TH2F*>    h2_extrapolate_pos_fit;
  vector<TH1F*>    num_eff_pos; 
  vector<TH1F*>    den_eff_pos;
  vector<TH1F*>    eff_pos;
  ULong64_t        evt;
  UInt_t           error;
  Int_t            daqTimeSec;
  Int_t            daqTimeMicroSec;
  Int_t            srsTimeStamp;
  UInt_t           srsTrigger;
  vector<double>   *ped_stdev;
  vector<double>   *ped_mean;
  vector<char>     *ped_readout;
  vector<string>   *ped_chamb;
  vector<int>      *ped_Strip;
  vector<int>      *ped_layer;
  vector<unsigned int>   *srsFec;
  vector<unsigned int>   *srsChip;
  vector<unsigned int>   *srsChan;
  vector<string>    *mmChamber;
  vector<int>       *mmLayer;
  vector<char>      *mmReadout;
  vector<int>       *mmStrip;
  vector<vector<short> > *raw_q;
  vector<short>     *max_q;
  vector<int>       *t_max_q;
  vector<int>       ped_strip_tmm;
  vector<string>    ped_chamb_tmm;
  vector<int>       strip;
  vector<int>       cluster_max_q;
  vector<string>    chamb;
  vector<int>       pcb;
  vector<double>    centroid;
  vector<int>       q_sum;
  vector<int>       n_cluster;

  
  // List of branches
  TBranch          *b_evt;               //!
  TBranch          *b_error;             //!
  TBranch          *b_daqTimeSec;        //!
  TBranch          *b_daqTimeMicroSec;   //!
  TBranch          *b_srsTimeStamp;      //!
  TBranch          *b_srsTrigger;        //!
  TBranch          *b_srsFec;            //!
  TBranch          *b_srsChip;           //!
  TBranch          *b_srsChan;           //!
  TBranch          *b_mmChamber;         //!
  TBranch          *b_mmLayer;           //!
  TBranch          *b_mmReadout;         //!
  TBranch          *b_mmStrip;           //!
  TBranch          *b_raw_q;             //!
  TBranch          *b_max_q;             //!
  TBranch          *b_t_max_q;           //!

  prova(TTree *tree=nullptr);
  virtual ~prova();
  virtual Int_t    Cut(Long64_t entry);
  virtual Int_t    GetEntry(Long64_t entry);
  virtual Long64_t LoadTree(Long64_t entry);
  virtual void     Init(TTree *tree);
  virtual void     Loop();
  virtual Bool_t   Notify();
  virtual void     Show(Long64_t entry = -1);
  //virtual void   Cluster(map<int,int>&, string&, int&);
  virtual void     Mapping(map<int,int>&,string&);  
  virtual void     Output_ped(ofstream&);
  virtual void     Output_strip(ofstream&);
  //virtual void    Histo_PCB(vector<TH1F*>&, vector<TH1F*>&, vector<TH1F*>&, map<int,int>&, int strip_pcb[], int );
  virtual void     Inizializ();
  virtual bool     DAQCount(string&);
  virtual void     Efficiency_Fit_1();
  virtual void     Efficiency_Fit_2();
  virtual void     Efficiency_Fit_3();
  virtual void     Efficiency_Fit_4();
  virtual void     FitParameters(double&, double&, double&, double&, double&, double&, double par[]);
  virtual void     Efficiency_1();
  virtual void     Efficiency_2();
  virtual void     Efficiency_stereo();
  virtual void     parameters(double&, double&, double&, double&, double par[]);
  virtual void     Efficiency_Position_1();
  virtual void     Efficiency_Position_2();
  virtual void     Efficiency_Position_Stereo();
  virtual int      Stereo_Road(double&);
  virtual void     Histo_Layer(vector<TH1F*>&, vector<TH1F*>&, vector<TProfile*>&, vector<TH2F*>&, vector<TH1F*>&,  map<int,int>&, int );
  virtual void     mapINFNcards();
  //virtual void   Histo_ped(map<int,int>&,vector<TH1F*>&, string&);

private:
  TString          _outdirname;
  TString          _outdirname_pcb;
  TString          _outdirname_ped, _outdirname_lay;

};

#endif