NHitsPerEvent.C File Reference

#include <map>

Go to the source code of this file.

Functions
void	Init ()
void	SetEntryNr (Int_t nr)
void	Next ()
void	Draw ()
void	ExtractData (TClonesArray *array, TGraph *g)
int	NHitsPerEvent ()

Variables
Int_t	gEntryNr = 0
TTree *	t
TClonesArray *	mvdHitArray = new TClonesArray("PndSdsHit")
TClonesArray *	mvdProjectedArray = new TClonesArray("PndSdsHit")
TClonesArray *	trackArray = new TClonesArray("PndTrack")
TGraph *	gHitPoints = new TGraph(4)
TGraph *	gProjections = new TGraph(4)
TMultiGraph *	mg = new TMultiGraph()
TH1I *	h1 = new TH1I("h1", "Hits per Event", 101, -0.5, 100.5)

Function Documentation

void Draw

(

)

Definition at line 64 of file NHitsPerEvent.C.

           {

        std::cout << "mvdHitArray: " << mvdHitArray->GetEntries() << " projectedHitArray " << mvdProjectedArray->GetEntries() << std::endl;

//      gHitPoints->Clear();
        ExtractData(mvdHitArray, gHitPoints);
        ExtractData(mvdProjectedArray, gProjections);
        //if (gHitPoints->GetN() != x.size()) gHitPoints->Expand(x.size());

        std::cout << "NHits: " << gHitPoints->GetN() << std::endl;
        gHitPoints->Draw("apl");
        gProjections->Draw("pl");
}

void ExtractData	(	TClonesArray *	array,
		TGraph *	g
	)

Definition at line 78 of file NHitsPerEvent.C.

References i, x, and y.

Referenced by Draw().

                                                {
        std::map<Double_t, Double_t>  x;
        std::map<Double_t, Double_t> y;
        for (int i = 0; i < array->GetEntriesFast(); i++){
                PndSdsHit* myHit = (PndSdsHit*)array->At(i);
                x[myHit->GetZ()] = myHit->GetX();
                y[myHit->GetZ()] = myHit->GetY();
                std::cout << i  << " : (" << myHit->GetX() << "/" << myHit->GetY() << "/" << myHit->GetZ() << ")" << std::endl;
        }
        for (int j = 0; j < g->GetN(); j++) g->RemovePoint(j);

        std::map<Double_t, Double_t >::iterator yit = y.begin();
        int i = 0;
        for (std::map<Double_t, Double_t >::iterator xit = x.begin() ; xit != x.end(); xit++, yit++, i++){
                g->SetPoint(i, xit->second, yit->second);
        }
}

void Init

(

)

Initialize the session

Definition at line 20 of file NHitsPerEvent.C.

{
        TFile* f = new TFile("2014-10-27-00-06-02_50_MHz_Filtered_PosCorrected_Event.root", "READ");
        t = (TTree*)f->Get("pndsim");
//      t->AddFriend("pndsim", "2014-10-27-00-06-02_50_MHz_Filtered_PosCorrected_Residuals.root");
//      t->AddFriend("pndsim", "2014-10-27-00-06-02_50_MHz_Filtered_PosCorrected_Track.root");
//      t->AddFriend("pndsim", "2014-10-27-00-06-02_50_MHz_Filtered_TimeCorrectedSorted_Event.root");
        t->StartViewer();

        t->SetBranchAddress("MVDHitsPixelPosCorrected", &mvdHitArray);
//      t->SetBranchAddress("MvdProjectedHit", &mvdProjectedArray);
//      t->SetBranchAddress("MvdTrack", &trackArray);

        gHitPoints->SetLineColor(kRed);
        gHitPoints->SetMarkerStyle(2);
        gHitPoints->SetMarkerColor(kRed);

        gProjections->SetLineColor(kBlue);
        gProjections->SetMarkerStyle(5);
        gProjections->SetMarkerColor(kBlue);


//      mg->Add(gProjections, "Projected");
        mg->Add(gHitPoints, "Hits");
        t->GetEntry(0);
}

void Next

(

)

Definition at line 52 of file NHitsPerEvent.C.

References Bool_t, Draw(), gEntryNr, mvdProjectedArray, run, and t.

           {
        Bool_t run = kTRUE;
        while(run){
                t->GetEntry(++gEntryNr);
                if (mvdProjectedArray->GetEntriesFast() > 0){
                        run = kFALSE;
                }
        }
        std::cout << "NextEntry: " << gEntryNr << std::endl;
        Draw();
}

int NHitsPerEvent

(

)

Definition at line 96 of file NHitsPerEvent.C.

References h1, i, Init(), mvdHitArray, and t.

{
        Init();
        for (int i = 0; i < t->GetEntries(); i++){
                t->GetEntry(i);
                h1->Fill(mvdHitArray->GetEntriesFast());
        }
        h1->Draw();
  return 0;
}

void SetEntryNr

(

Int_t

nr

)

Definition at line 47 of file NHitsPerEvent.C.

References gEntryNr, and t.

                         {
        gEntryNr = nr;
        t->GetEntry(gEntryNr);
}

Variable Documentation

Int_t gEntryNr = 0

Definition at line 10 of file NHitsPerEvent.C.

Referenced by Next(), and SetEntryNr().

TGraph* gHitPoints = new TGraph(4)

Definition at line 15 of file NHitsPerEvent.C.

Referenced by Draw(), and Init().

TGraph* gProjections = new TGraph(4)

Definition at line 16 of file NHitsPerEvent.C.

Referenced by Draw(), and Init().

TH1I* h1 = new TH1I("h1", "Hits per Event", 101, -0.5, 100.5)

Definition at line 18 of file NHitsPerEvent.C.

TMultiGraph* mg = new TMultiGraph()

Definition at line 17 of file NHitsPerEvent.C.

Referenced by Init(), and plot_hitpattern().

TClonesArray* mvdHitArray = new TClonesArray("PndSdsHit")

Definition at line 12 of file NHitsPerEvent.C.

Referenced by Draw(), Init(), and NHitsPerEvent().

TClonesArray* mvdProjectedArray = new TClonesArray("PndSdsHit")

Definition at line 13 of file NHitsPerEvent.C.

Referenced by Draw(), Init(), and Next().

TTree* t

Definition at line 11 of file NHitsPerEvent.C.

TClonesArray* trackArray = new TClonesArray("PndTrack")

Definition at line 14 of file NHitsPerEvent.C.

Referenced by OnlineDBuilderTask::Exec(), and Init().