reco_analys4.C File Reference

Go to the source code of this file.

Functions
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/basiclibs.C")
	rootlogon ()
	basiclibs ()
t	SetBranchAddress ("EmcBump",&bump_array)
tsim	SetBranchAddress ("MCTrack",&mctrack_array)
	for (Int_t j=0;j< t->GetEntriesFast();j++)
c1	Divide (1, 2)
c1	cd (1)
h1	SetTitle ("Cluster energy of 1 GeV photon")
h1	GetXaxis () -> SetTitle("Energy, GeV")
h1	Draw ()
c1	cd (2)
h2	SetTitle ("Cluster energy of 1 GeV photon (corrected)")

Variables
TFile *	f = new TFile("full_emc.root")
TTree *	t =(TTree *) f->Get("pndsim")
TClonesArray *	bump_array =new TClonesArray("PndEmcBump")
TFile *	fsim = new TFile("sim_emc.root")
TTree *	tsim =(TTree *) fsim->Get("pndsim")
TClonesArray *	mctrack_array =new TClonesArray("PndMCTrack")
double	bump_energy
double	bump_energy_corrected
TH1F *	h1 = new TH1F("h1","Cluster energy",100,0.0,1.10)
TH1F *	h2 = new TH1F("h2","Cluster energy corrected",100,0.0,1.10)
TCanvas *	c1 = new TCanvas("c1", "Cluster Energy", 100, 100, 800, 800)

Function Documentation

basiclibs

(

)

c1 cd

(

1

)

c1 cd

(

2

)

c1 Divide	(	1	,
		2
	)

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 29 of file reco_analys4.C.

References PndEmcCluster::energy(), PndEmcCluster::GetEnergyCorrected(), and i.

        {
                t->GetEntry(j);
                for (Int_t i=0; i<bump_array->GetEntriesFast(); i++)
                {
                        PndEmcBump *bump=(PndEmcBump*)bump_array->At(i);
                        bump_energy=bump->energy();
                        bump_energy_corrected=bump->GetEnergyCorrected();

                        if (bump_energy>0.1){
                                h1->Fill(bump_energy);
                                h2->Fill(bump_energy_corrected);
                        }
                }
        }

h1 GetXaxis

(

)

-> SetTitle("Energy, GeV")

Definition at line 653 of file anal_hit_digi_cluster_fwendcap.C.

{Reconstructed}-#theta_{MC} (#circ)");

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	(	"EmcBump"	,
		&	bump_array
	)

tsim SetBranchAddress	(	"MCTrack"	,
		&	mctrack_array
	)

h1 SetTitle

(

"Cluster energy of 1 GeV photon"

)

h2 SetTitle

(

"Cluster energy of 1 GeV photon (corrected)"

)

Variable Documentation

TClonesArray* bump_array =new TClonesArray("PndEmcBump")

Definition at line 11 of file reco_analys4.C.

double bump_energy

Definition at line 21 of file reco_analys4.C.

double bump_energy_corrected

Definition at line 22 of file reco_analys4.C.

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

Definition at line 45 of file reco_analys4.C.

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

Definition at line 9 of file reco_analys4.C.

TFile* fsim = new TFile("sim_emc.root")

Definition at line 14 of file reco_analys4.C.

TH1F* h1 = new TH1F("h1","Cluster energy",100,0.0,1.10)

Definition at line 24 of file reco_analys4.C.

TH1F* h2 = new TH1F("h2","Cluster energy corrected",100,0.0,1.10)

Definition at line 25 of file reco_analys4.C.

TClonesArray* mctrack_array =new TClonesArray("PndMCTrack")

Definition at line 18 of file reco_analys4.C.

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

Definition at line 10 of file reco_analys4.C.

TTree* tsim =(TTree *) fsim->Get("pndsim")

Definition at line 15 of file reco_analys4.C.