HMdcLowLevelQA

class HMdcLowLevelQA : public HReconstructor

    public:

                      HMdcLowLevelQA()
                      HMdcLowLevelQA(Text_t* name, Text_t* title)
                      ~HMdcLowLevelQA()
                 void checkNeighboringCells(Int_t s, Int_t m, Int_t c)
       static TClass* Class()
                 void createHistograms()
        virtual Int_t execute()
                 void fillHistogramsCal(Int_t, Int_t, Int_t, Int_t, Float_t, Float_t)
                 void fillHistogramsRaw(Int_t, Int_t, Int_t, Int_t, Float_t, Float_t)
       virtual Bool_t finalize()
               Bool_t getCategories()
                 void getSetup()
       virtual Bool_t init()
                 void initVariables()
      virtual TClass* IsA() const
                 void lookupSlopes(Int_t s, Int_t m)
          TDirectory* Mkdir(TDirectory*, Char_t*, Int_t, Int_t p = 1)
                 void printStatus()
       virtual Bool_t reinit()
                 void resetCounters()
                 void setOutputRoot(Char_t*)
                 void setParContainers()
                 void setStepSize(Int_t size)
         virtual void ShowMembers(TMemberInspector& insp, char* parent)
         virtual void Streamer(TBuffer& b)
                 void StreamerNVirtual(TBuffer& b)
                 void testTimeCuts(Int_t s, Int_t m, Float_t t1, Float_t t2)
                 void writeHistograms(Int_t, Int_t, Int_t, Int_t, Int_t, Int_t)
                 void writeTimeCuts()

Data Members

    protected:

                TH1F* hmdc_time1_everithing_summ[6][4]   Regardless of any time cuts summary for the file
                TH1F* hmdc_time1_good_t1_summ[6][4]      All cuts are OK, summary for the file
                TH1F* hmdc_time1_not_t1_summ[6][4]       Time1 cut is not done, summary for the file
                TH1F* hmdc_trend_meantime1[6][4]         PASS        //trend of mean of time1 per sector and module
                TH1F* hmdc_time1[6][4]                   all t-cuts  //time1 per sector and module
                TH1F* hmdc_trend_meantime1not_t1[6][4]   DOES NOT    //trend of mean of time1 per sector and module
                TH1F* hmdc_time1not_t1[6][4]             PASS t1-cut //time1 per sector and module
                TH1F* hmdc_trend_meantime1not_t2[6][4]   DOES NOT    //trend of mean of time1 per sector and module
                TH1F* hmdc_time1not_t2[6][4]             PASS t2-cut //time1 per sector and module
                TH1F* hmdc_trend_meantime1not_t12[6][4]  DOES NOT    //trend of mean of time1 per sector and module
                TH1F* hmdc_time1not_t12[6][4]            PASS t12cut //time1 per sector and module
                TH1F* hmdc_trend_meantime1_plane[4]      trend of mean of time1 per plane
                TH1F* hmdc_time1_plane[4]                time1 per plane
                TH1F* hmdc_trend_meanSlope[6][4]         slope of CAL1 as a function of file
                TH1F* hmdc_Slope[6][4]                   technical histogramm for calculating slopes mean value for trend histogramm
                TH1F* hmdc_trend_neighbors[6][4]         trend of ratio single to neighboring cells per sector and module
                Int_t previousCellStatus                 buffer for cell number of previous hit
                Int_t double_all[6][4][2]                keeps conter for cells with neighbors and all cells
                TH1F* hmdc_trend_123ratio[6][4]          for 6 modules in 4 planes
                Int_t cnt123[6][4]                       number of "noise" wires (we assume this to be a noise if only 1 or 2 or 3 wires have fired in a module)
                Int_t sect_module_NumberOfHits[6][4]     number of hits per module=(sector,plane)
                Int_t sm_old[2]                          sector => sm_old[0], module =>  sm_old[1]
                TH1F* hmdc_multiplicity[6][4]            hits multiplicity event by event storage
                TH1F* hmdc_trend_multiplicity[6][4]      mean value over N-step events of multiplicity
                TH1F* a_thisFileName                     
                TH2F* hmdc_summ_MBOvsTDC[6][4]           
                TH1F* hmdc_summ_dNdNhits[6][4]           
                TH2F* hmdc_summ_t2mt1_vs_t1[6][4]        
                TH1F* hmdc_summ_tCuts_info[4]            histogram contains 6 bins for each time cut
                TH1F* hmdc_summ_dNdWire[6][4][6][3]      sector//plane//layer//histogramm_kind_(four different histograms: for different time cut conditions)
                TH1F* hmdc_summ_trend_meantime1_plane    
                TH1F* hmdc_summ_trend_meantime1          
                TH1F* hmdc_summ_trend_meantime1not_t1    
                TH1F* hmdc_summ_trend_meantime1not_t2    
                TH1F* hmdc_summ_trend_meantime1not_t12   
                TH1F* hmdc_summ_trend_neighbors          
                TH1F* hmdc_summ_trend_123ratio           
                TH1F* hmdc_summ_trend_multiplicity       
                TH1F* hmdc_trend_NhitsPERmbo[6][4][16]   for 6 modules in 4 planes
                Int_t NhitsPERmbo[6][4][16]              for 6 modules in 4 planes
                Int_t smmb_old[3]                        [sector,module,mbo] combination
              Char_t* fNameRootHists                     file name of root output file
         static Int_t colorsSector[6]                    colors for hists to be compared in one plane
         static Int_t colorsModule[4]                    colors for hists to be compared in one sector
         static Int_t colorsCuts[3]                      colors for hists to be compared for cuts
                Int_t Eventcounter                       counts the number of events
                Int_t step_trend_meantime1_counter       counts the number of bin to be filled in trend hists
                Int_t stepsize_trend_meantime1           defines the number of events to be filled for one data point
                Int_t numberOfSectorsPerModule[4]        stores number of chambers per plane
                Int_t mySetup[6][4]                      actual setup of the detector
         static Int_t cutResults[4]                      contains 0/1 for cuts in t1,t2 and t2-1 after testTimeCuts()
               Bool_t isInitialized                      flag for initialization of histograms
                Int_t nReqEvents                         nuber of Events from Hades::EventLoop
                Int_t nBinsTrendX                        max number of bins for trend histograms
           HCategory* calCat                             ! pointer to the cal data
           HCategory* rawCat                             ! pointer to the raw data
            HLocation loccal1                            ! location for new object.
            HLocation locraw                             ! location for new object.
           HIterator* itercal1                           ! iterator on cal1 data.
           HIterator* iterraw                            ! iterator on raw data.
      HMdcLookupGeom* lookup                             ! lookup table for mapping
       HMdcLookupRaw* lookupRaw                          ! lookup table for mapping
         HMdcTimeCut* timecut                            ! container for time cuts

Class Description

 HMdcLowLevelQA

 Default constructor calls the function setParContainers().

 Constructor calls the constructor of class HReconstructor with the name
 and the title as arguments. It also calls the function setParContainers().

 destructor deletes the iterator

 creates the parameter containers MdcLookupGeom/MdcLookupRaw/MdcTimeCut if they do not
 exist and adds them to the list of parameter containers in the runtime
 database

 Sets rootfile output of HMdcFillTimeHists where all created histograms were written.

 creates the MdcRaw  and MdcCal1 categories and adds them to the current
 event
 creates an iterator which loops over all fired cells
 calls the function setParContainers()

for(Int_t s=0; s<6; s++){ //time cut is the same for the same plane (different sectors)

cout << "s="<<s<<" m="<<m<<" Slope="<<calparraw->getSlope(s,m) << endl;

 Writes hists to file dependend on level (0=raw,1=cal), main (1=summary,2=plane),
 sector,module,layer/mb,cell/channel.-99 means that this ind3ex is not used.

 hmdc_summ_MBOvsTDC (independent from any cuts)

 Function to create subdirectories

cout << "----   ---- event ---- ----" << endl;

Inline Functions

               void resetCounters()
               void checkNeighboringCells(Int_t s, Int_t m, Int_t c)
               void setStepSize(Int_t size)
            TClass* Class()
            TClass* IsA() const
               void ShowMembers(TMemberInspector& insp, char* parent)
               void Streamer(TBuffer& b)
               void StreamerNVirtual(TBuffer& b)