ROOT logo
#include "hshowerrawhist.h"
#include "hshowercalhist.h"
#include "hshowerhithist.h"
#include "hruntimedb.h"
#include "hevent.h"
#include "hspectrometer.h"
#include "hdetector.h"
#include "hshowerdetector.h"
#include "hcategory.h"
#include "hmatrixcatiter.h"
#include "hlocation.h"
#include "hshowerraw.h"
#include "hshowercal.h"
#include "hdebug.h"
#include "hades.h"
#include "hiterator.h"
#include "showerdef.h"

#include "TH1.h"
#include "TH2.h"

ClassImp(HShowerHitHist)

    HShowerHitHist::~HShowerHitHist(void) {
	if (fIter) delete fIter;
    }

Bool_t HShowerHitHist::init() {
    printf("initialization of shower hist\n");
    m_nEvents = 0;

    m_pHitCat=gHades->getCurrentEvent()->getCategory(catShowerHit);
    if (m_pHitCat) {
	fIter = (HIterator*)m_pHitCat->MakeIterator();
    }
    bookHist();
    return kTRUE;
}

Bool_t HShowerHitHist::finalize(void) {
    finalizeHist();
    return kTRUE;
}

Int_t HShowerHitHist::execute()
{
    if(fIter == 0) return 0; // noting todo
    HShowerHit *pHit;

    Int_t n = 0;

    fIter->Reset();
    while((pHit = (HShowerHit *)fIter->Next()))
    {
	fillHist(pHit);
	n++;
    }

    m_nEvents++;
    return 0;
}

/***********************************************************/
/* definition of private function                          */
/***********************************************************/

Bool_t HShowerHitHist::bookHist() {
    Char_t name[80];
    Char_t title[80];

    m_pChargeHitHist = new TH1F("ShowerHitCharge", "Shower Hit Level - Charge",
				128, 0, 256);

    for(Int_t i = 0; i < 3; i++) {
	sprintf(name, "ShowerHitCharge%d", i);
	sprintf(title, "Shower Hit Level - Charge in Module %d", i);
	m_pChargeHitModHist[i] = new TH1F(name, title, 128, 0, 256);

	sprintf(name, "ShowerHitAvg%d", i);
	sprintf(title, "Shower Hit Level - Average Charge in Module %d", i);
	m_pChargeHitAvgHist[i] =new TH2F(name, title, 32, 0, 32, 32, 0, 32);
	//  m_pChargeHitAvgHist[i]->SetDrawOption(colz);
	sprintf(name, "ShowerHitFreq%d", i);
	sprintf(title, "Shower Hit Level - Freq in Module %d", i);
	m_pChargeHitFreqHist[i] =new TH2F(name, title, 32, 0, 32, 32, 0, 32);
    }
    return kTRUE;
}


Bool_t HShowerHitHist::fillHist(Int_t nModule, Int_t nRow,
				Int_t nColumn, Float_t fCharge) {

    m_pChargeHitHist->Fill(fCharge);
    m_pChargeHitModHist[nModule]->Fill(fCharge);

    m_pChargeHitFreqHist[nModule]->Fill(nColumn, nRow, 1 );
    m_pChargeHitAvgHist[nModule]->Fill(nColumn, nRow, fCharge);

    return 1;
}

Bool_t HShowerHitHist::fillHist(HShowerHit* pHit) {
    Int_t mod = pHit->getModule();
    Int_t row = pHit->getRow();
    Int_t col = pHit->getCol();
    Float_t charge = pHit->getCharge();

    fillHist(mod, row, col, charge);

    return kTRUE;
}


Bool_t HShowerHitHist::finalizeHist() {
    for(Int_t i = 0; i < 3; i++) {
	m_pChargeHitAvgHist[i]->Divide(m_pChargeHitFreqHist[i]);
	m_pChargeHitFreqHist[i]->Scale(1.0/m_nEvents);
    }
    writeHist();
    return kTRUE;
}

Bool_t HShowerHitHist::writeHist() {
    printf("writing histograms ...\n");
    m_pChargeHitHist->Write();
    for(Int_t i = 0; i < 3; i++) {
	m_pChargeHitModHist[i]->Write();
	m_pChargeHitFreqHist[i]->Write();
	m_pChargeHitAvgHist[i]->Write();
    }

    return kTRUE;
}

 hshowerhithist.cc:1
 hshowerhithist.cc:2
 hshowerhithist.cc:3
 hshowerhithist.cc:4
 hshowerhithist.cc:5
 hshowerhithist.cc:6
 hshowerhithist.cc:7
 hshowerhithist.cc:8
 hshowerhithist.cc:9
 hshowerhithist.cc:10
 hshowerhithist.cc:11
 hshowerhithist.cc:12
 hshowerhithist.cc:13
 hshowerhithist.cc:14
 hshowerhithist.cc:15
 hshowerhithist.cc:16
 hshowerhithist.cc:17
 hshowerhithist.cc:18
 hshowerhithist.cc:19
 hshowerhithist.cc:20
 hshowerhithist.cc:21
 hshowerhithist.cc:22
 hshowerhithist.cc:23
 hshowerhithist.cc:24
 hshowerhithist.cc:25
 hshowerhithist.cc:26
 hshowerhithist.cc:27
 hshowerhithist.cc:28
 hshowerhithist.cc:29
 hshowerhithist.cc:30
 hshowerhithist.cc:31
 hshowerhithist.cc:32
 hshowerhithist.cc:33
 hshowerhithist.cc:34
 hshowerhithist.cc:35
 hshowerhithist.cc:36
 hshowerhithist.cc:37
 hshowerhithist.cc:38
 hshowerhithist.cc:39
 hshowerhithist.cc:40
 hshowerhithist.cc:41
 hshowerhithist.cc:42
 hshowerhithist.cc:43
 hshowerhithist.cc:44
 hshowerhithist.cc:45
 hshowerhithist.cc:46
 hshowerhithist.cc:47
 hshowerhithist.cc:48
 hshowerhithist.cc:49
 hshowerhithist.cc:50
 hshowerhithist.cc:51
 hshowerhithist.cc:52
 hshowerhithist.cc:53
 hshowerhithist.cc:54
 hshowerhithist.cc:55
 hshowerhithist.cc:56
 hshowerhithist.cc:57
 hshowerhithist.cc:58
 hshowerhithist.cc:59
 hshowerhithist.cc:60
 hshowerhithist.cc:61
 hshowerhithist.cc:62
 hshowerhithist.cc:63
 hshowerhithist.cc:64
 hshowerhithist.cc:65
 hshowerhithist.cc:66
 hshowerhithist.cc:67
 hshowerhithist.cc:68
 hshowerhithist.cc:69
 hshowerhithist.cc:70
 hshowerhithist.cc:71
 hshowerhithist.cc:72
 hshowerhithist.cc:73
 hshowerhithist.cc:74
 hshowerhithist.cc:75
 hshowerhithist.cc:76
 hshowerhithist.cc:77
 hshowerhithist.cc:78
 hshowerhithist.cc:79
 hshowerhithist.cc:80
 hshowerhithist.cc:81
 hshowerhithist.cc:82
 hshowerhithist.cc:83
 hshowerhithist.cc:84
 hshowerhithist.cc:85
 hshowerhithist.cc:86
 hshowerhithist.cc:87
 hshowerhithist.cc:88
 hshowerhithist.cc:89
 hshowerhithist.cc:90
 hshowerhithist.cc:91
 hshowerhithist.cc:92
 hshowerhithist.cc:93
 hshowerhithist.cc:94
 hshowerhithist.cc:95
 hshowerhithist.cc:96
 hshowerhithist.cc:97
 hshowerhithist.cc:98
 hshowerhithist.cc:99
 hshowerhithist.cc:100
 hshowerhithist.cc:101
 hshowerhithist.cc:102
 hshowerhithist.cc:103
 hshowerhithist.cc:104
 hshowerhithist.cc:105
 hshowerhithist.cc:106
 hshowerhithist.cc:107
 hshowerhithist.cc:108
 hshowerhithist.cc:109
 hshowerhithist.cc:110
 hshowerhithist.cc:111
 hshowerhithist.cc:112
 hshowerhithist.cc:113
 hshowerhithist.cc:114
 hshowerhithist.cc:115
 hshowerhithist.cc:116
 hshowerhithist.cc:117
 hshowerhithist.cc:118
 hshowerhithist.cc:119
 hshowerhithist.cc:120
 hshowerhithist.cc:121
 hshowerhithist.cc:122
 hshowerhithist.cc:123
 hshowerhithist.cc:124
 hshowerhithist.cc:125
 hshowerhithist.cc:126
 hshowerhithist.cc:127
 hshowerhithist.cc:128
 hshowerhithist.cc:129
 hshowerhithist.cc:130
 hshowerhithist.cc:131
 hshowerhithist.cc:132
 hshowerhithist.cc:133
 hshowerhithist.cc:134
 hshowerhithist.cc:135