EEezClusterQA

class EEezClusterQA : public TDirectory

    private:

      void bookHistograms()

    public:

                      EEezClusterQA(const Char_t* name, const Char_t* title = "EEMC cluster QA histograms")
              virtual ~EEezClusterQA()
                 void bookCut(const Char_t* name, const Char_t* title, Int_t nbin, Float_t min, Float_t max)
                 void bookCut(const Char_t* name, const Char_t* title, Int_t nbinx, Float_t minx, Float_t maxx, Int_t nbiny, Float_t miny, Float_t maxy)
       static TClass* Class()
                 void fail(EEezCluster* cluster, const Char_t* tag)
                 void Fill(EEezClusterPtrVec_t* clusters)
                 void Fill(EEezCluster* cluster)
                 void fillCut(const Char_t* name, Float_t value)
                 void fillCut(const Char_t* name, Float_t valuex, Float_t valuey)
                 void flag(EEezCluster* cluster, const Char_t* tag)
                 void Init()
      virtual TClass* IsA() const
                 void pass(EEezCluster* cluster, const Char_t* tag)
                 void save(EEezCluster* cluster, Int_t iuv, const Char_t* tag)
                 void save(TH1F* histo, Int_t iuv = 0, const Char_t* tag)
                 void setAnalysis(EEezAnalysis* a)
                 void setSmdProfiler(EEsmdProfile* p)
         virtual void ShowMembers(TMemberInspector& insp, char* parent)
         virtual void Streamer(TBuffer& b)
                 void StreamerNVirtual(TBuffer& b)

Data Members

    private:

             TList* m_ListOfCuts      
             TList* m_ListOf2DCuts    
              TH2F* m_XY              Tower-weighted y vs x
              TH2F* m_XY_smd          X,Y from mean SMD response
              TH1F* m_Energy          Tower energy spectrum
              TH2F* m_EcVsEtot        Cluster energy vs sector total energy
              TH1F* m_HitStrips[2]    Frequency plot of hit strips vs index
              TH2F* m_Sum9VsEc[2]     Sum of 9 strips near peak vs eCluster
              TH2F* m_Sum5VsEc[2]     Sum of 5 strips near peak vs eCluster
              TH1F* m_Towers          Tower energy response
              TH2F* m_Sum5            Sum over 5 strips U vs V
              TH2F* m_Sum9            Sum over 9 strips U vs V
              TH1F* m_Frac5to9[2]     Fractional sum of 5 bins to 9 bins
              TH2F* m_Yield           Fit yield U vs V
              TH2F* m_Sigma           Fit sigma U vs V
              TH2F* m_RedChi2         Fit chi2/ndf U vs V
              TH2F* m_YieldVsEc[2]    Fit yield versus Ecluster
              TH1F* m_YieldToSum[2]   Fit yield to sum over tower
              TH1F* m_FreqBigChi2[2]  Frequency plot of strip w/ largest chi2 from fit
              TH1F* m_FreqEmpty[2]    Frequency plot that strip shows no response in fit 
              TH1F* m_MipStrips[2]    Cumulative # of mips
              TH1F* m_SeedAdc[4]      ADC response of TPQR for seeds
              TH1F* m_NeighborAdc[4]  ADC response of TPQR for neighboring towers
              TH2F* m_SumVsEc[2]      Sum of number of mips in seed tower vs Ec
          TProfile* m_SumVsEc_pfx[2]  Sum of number of mips in seed tower vs Ec
              TH2F* m_SumSmd          Simple sum of U plane vs V plane beneath seed tower
              TH1F* m_Parameters[5]   Histogram of fit parameters
      EEezAnalysis* m_EEezAnaly       
      EEsmdProfile* m_EEsmdProf       
       EEmcSmdGeom* m_EEsmdGeom       
        TDirectory* m_UFail           
        TDirectory* m_VFail           
        TDirectory* m_UPass           
        TDirectory* m_VPass           
        TDirectory* m_UFlag           
        TDirectory* m_VFlag           

Class Description

 EEezClusterQA

 Author: Jason C. Webb <jwebb@iucf.indiana.edu>

 EEezClusterQA provides a set of QA histograms which may be filled
 using a single call to EEezClusterQA::Fill().  The arguement of fill
 may either be a pointer to a cluster or a vector of cluster pointers.

 The code assumes that the user has set a pointer to the EEsmdProfile
 object, via setSmdProfiler().  An instance of EEezAnalysis should
 also be set, via setAnalysis(), but is not required.

 Class constructor

 Initialize the class

 Book all QA histograms

 Given a pointer to a cluster, fill the QA histograms

 Loop over all clusters in the vector, filling histograms
 from the clusters contained therein.

 Save the SMD profile and associated histograms in a directory
 as an example of a cluster which passed a cut.  One should not
 call this on every event, just a subset of the data.

 Save the SMD profile and associated histograms in a directory
 as an example of a cluster which failed a cut.  One should not
 call this on every event, just a subset of the data.

 Save the specified cluster in the "flag" directory, just as
 with pass and fail.  Here we do not limit ourselves to the
 first 50 clusters.  Use with care.

 Saves the specified cluster's SMD profile (and other assoicated
 histograms) in the current directory.

 Saves the specified histogram

 Book a histogram for quantities which are cut on.   Histograms are
 stored in a TList, which is searched whenever we fill.  This is
 slow, and should be changed to a hash table.

 Fills the specified cut

 Book a histogram for quantities which are cut on.   Histograms are
 stored in a TList, which is searched whenever we fill.  This is
 slow, and should be changed to a hash table.

 Fills the specified cut
  ((TH1F*)m_ListOf2DCuts -> FindObject( name )) -> Fill(value,valuey);

Inline Functions

               void setAnalysis(EEezAnalysis* a)
               void setSmdProfiler(EEsmdProfile* p)
            TClass* Class()
            TClass* IsA() const
               void ShowMembers(TMemberInspector& insp, char* parent)
               void Streamer(TBuffer& b)
               void StreamerNVirtual(TBuffer& b)
               void ~EEezClusterQA()