hit_analys.C File Reference

Go to the source code of this file.

Functions
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/basiclibs.C")
	rootlogon ()
	basiclibs ()
t	SetBranchAddress ("EmcHit",&hit_array)
	for (Int_t j=0;j< t->GetEntriesFast();j++)
h1	Draw ()
c1	Update ()

Variables
TFile *	f = new TFile("hit_emc.root")
TTree *	t =(TTree *) f->Get("pndsim")
TClonesArray *	hit_array =new TClonesArray("PndEmcHit")
double	hit_theta
double	point_theta
TH1F *	h1 = new TH1F("h1","Theta",100,-0.1,0.1)
TCanvas *	c1 = new TCanvas("c1", "", 100, 100, 800, 800)

Function Documentation

basiclibs

(

)

h1 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");
}

for

(

Int_t

j = 0; j< t->GetEntriesFast(); j++

)

Definition at line 20 of file hit_analys.C.

References PndEmcHit::GetPointList(), hit(), and i.

        {
                t->GetEntry(j);
                for (Int_t i=0; i<hit_array->GetEntriesFast(); i++)
                {
                        PndEmcHit *hit=(PndEmcHit*)hit_array->At(i);
                        TVector3 hit_pos(hit->GetX(),hit->GetY(),hit->GetZ());
                        hit_theta=hit_pos.Theta();

                        std::vector<PndEmcPoint*> pointList=hit->GetPointList();
                        if (pointList.size()!=0)
                        {
                                TVector3 point_pos(pointList[0]->GetX(),pointList[0]->GetY(),pointList[0]->GetZ());
                                point_theta=point_pos.Theta();
                                h1->Fill(hit_theta-point_theta);
                        }
                }
        }

gROOT LoadMacro

(

"$VMCWORKDIR/gconfig/basiclibs.C"

)

rootlogon

(

)

Definition at line 1 of file outdated/mpiTools/macros/emc/rootlogon.C.

{
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcDigi");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcMC");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcTools");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcReco");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcData");
  gSystem->AddIncludePath("-I$VMCWORKDIR/pnddata");
  gSystem->AddIncludePath("-I$VMCWORKDIR/base");
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();

  gSystem->Load("libDpmEvtGen.so");
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libTrkBase");
  gSystem->Load("libPndData");
  gSystem->Load("libField");
  gSystem->Load("libPassive");
  gSystem->Load("libGen");
  gSystem->Load("libPGen");
  gSystem->Load("libEmc"); 

  delete gRandom;
  gRandom=new TRandom3(0);
}

t SetBranchAddress	(	"EmcHit"	,
		&	hit_array
	)

c1 Update

(

)

Variable Documentation

TCanvas* c1 = new TCanvas("c1", "", 100, 100, 800, 800)

Definition at line 39 of file hit_analys.C.

TFile* f = new TFile("hit_emc.root")

Definition at line 9 of file hit_analys.C.

TH1F* h1 = new TH1F("h1","Theta",100,-0.1,0.1)

Definition at line 18 of file hit_analys.C.

TClonesArray* hit_array =new TClonesArray("PndEmcHit")

Definition at line 12 of file hit_analys.C.

double hit_theta

Definition at line 16 of file hit_analys.C.

double point_theta

Definition at line 16 of file hit_analys.C.

TTree* t =(TTree *) f->Get("pndsim")

Definition at line 10 of file hit_analys.C.