reco_apd.C File Reference

Go to the source code of this file.

Functions
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/basiclibs.C")
	rootlogon ()
	basiclibs ()
t	SetBranchAddress ("EmcApdHit",&hit_array)
tsim	SetBranchAddress ("MCTrack",&mctrack_array)
hEnergy	GetXaxis () -> SetTitle("Energy Loss [MeV]")
	for (Int_t j=0;j< t->GetEntriesFast();j++)
hEnergy	Draw ()

Variables
TFile *	f = new TFile("hit_emc_apd.root")
TTree *	t =(TTree *) f->Get("pndsim")
TClonesArray *	hit_array =new TClonesArray("PndEmcApdHit")
TFile *	fsim = new TFile("sim_emc_apd.root")
TTree *	tsim =(TTree *) fsim->Get("pndsim")
TClonesArray *	mctrack_array =new TClonesArray("PndMCTrack")
Float_t	energy = 0
Int_t	ntrack = 0
TH1F *	hEnergy = new TH1F("hEnergy","Energy Loss",200,0,1)
TH1F *	hMult = new TH1F("hMult","Particle Multiplicity",15,-0.5,14.5)
TH1F *	hMom = new TH1F("hMom","Particle Momentum",150,0,1.5)
TCanvas *	c1 = new TCanvas("c1", "Energy Loss in APDs", 100, 100, 800, 800)
TCanvas *	c2 = new TCanvas("c2", "Particles in APDs", 100, 100, 800, 800)
TCanvas *	c3 = new TCanvas("c3", "Particle Momentum in APDs", 100, 100, 800, 800)

Function Documentation

basiclibs

(

)

hEnergy 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_apd.C.

References PndEmcApdHit::GetEnergy(), PndEmcApdHit::GetMCIndex(), PndMCTrack::GetMomentum(), PndEmcApdHit::GetNPoints(), PndMCTrack::GetPdgCode(), hit_energy, i, mctrack, and ntrack.

{
  t->GetEntry(j);
  tsim->GetEntry(j);

    for (Int_t i=0; i<hit_array->GetEntriesFast(); i++)
      {
        PndEmcApdHit *apd=(PndEmcApdHit*)hit_array->At(i);
        hit_energy=apd->GetEnergy();
        ntrack = apd->GetNPoints();

        hEnergy->Fill(hit_energy*1000.);
        hMult->Fill(ntrack);
        cout << "# Particles in APD: " << ntrack << endl;
        for (Int_t tt=0; tt<ntrack; tt++)
          {
            cout << apd->GetMCIndex(tt) << endl;
            PndMCTrack* mctrack=(PndMCTrack *) mctrack_array->At(apd->GetMCIndex(tt));
            cout << "Particle ID: " << mctrack->GetPdgCode()<< endl;
            hMom->Fill(mctrack->GetMomentum().Mag());
          } // end of MC loop
      } // end of hit loop
}

hEnergy GetXaxis

(

)

-> SetTitle("Energy Loss [MeV]")

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	(	"EmcApdHit"	,
		&	hit_array
	)

tsim SetBranchAddress	(	"MCTrack"	,
		&	mctrack_array
	)

Variable Documentation

TCanvas* c1 = new TCanvas("c1", "Energy Loss in APDs", 100, 100, 800, 800)

Definition at line 55 of file reco_apd.C.

TCanvas* c2 = new TCanvas("c2", "Particles in APDs", 100, 100, 800, 800)

Definition at line 58 of file reco_apd.C.

TCanvas* c3 = new TCanvas("c3", "Particle Momentum in APDs", 100, 100, 800, 800)

Definition at line 61 of file reco_apd.C.

Float_t energy = 0

Definition at line 20 of file reco_apd.C.

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

Definition at line 9 of file reco_apd.C.

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

Definition at line 15 of file reco_apd.C.

TH1F* hEnergy = new TH1F("hEnergy","Energy Loss",200,0,1)

Definition at line 24 of file reco_apd.C.

TClonesArray* hit_array =new TClonesArray("PndEmcApdHit")

Definition at line 12 of file reco_apd.C.

TH1F* hMom = new TH1F("hMom","Particle Momentum",150,0,1.5)

Definition at line 26 of file reco_apd.C.

TH1F* hMult = new TH1F("hMult","Particle Multiplicity",15,-0.5,14.5)

Definition at line 25 of file reco_apd.C.

TClonesArray* mctrack_array =new TClonesArray("PndMCTrack")

Definition at line 17 of file reco_apd.C.

Int_t ntrack = 0

Definition at line 21 of file reco_apd.C.

Referenced by PndSciTAnaIdeal::Exec(), and for().

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

Definition at line 10 of file reco_apd.C.

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

Definition at line 16 of file reco_apd.C.