anaMvdPid.C File Reference

Go to the source code of this file.

Functions
gROOT	Macro ("$VMCWORKDIR/gconfig/rootlogon.C")
timer	Start ()
tree	SetBranchAddress ("PndMvdPidCand",&pointlist)
	for (int j=0;j< nEvents &&j< tree->GetEntriesFast();j++)
h1	Draw ()
timer	Stop ()

Variables
TStopwatch	timer
TFile *	inFile = new TFile("MvdPidIdeal.root","READ")
TTree *	tree = (TTree *)inFile->Get("pndsim")
TClonesArray *	pointlist =new TClonesArray("PndMvdPidCand")
TH2F *	h1 =new TH2F("h1", "likelihood count", 100, 0.0, 2.5, 100, 0.0, 1.0)
TH2F *	h2 =new TH2F("h2", "energy loss", 100, 0.0, 2.5, 100, 0.0, 0.05)
TH2F *	h3 =new TH2F("h3", "ideal likelihood count", 100, 0.0, 2.5, 100, 0.0, 1.0)
Int_t	nEvents = 10000
TCanvas *	c =new TCanvas("c2")
Double_t	rtime = timer.RealTime()
Double_t	ctime = timer.CpuTime()

Function Documentation

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

(

)

Definition at line 31 of file anaMvdPid.C.

References PndMvdAdvancedPidAlgo::CalcLikelihood(), dE, dx, PndMvdPidCand::GetLikelihood(), PndMvdPidCand::GetMvdHitdE(), PndMvdPidCand::GetMvdHitdx(), PndMvdPidCand::GetMvdHitMomentum(), PndMvdPidCand::GetMvdHits(), and i.

                                                         {
    tree->GetEntry(j);
    for(int i=0;i<pointlist->GetEntriesFast();i++) {
      PndMvdPidCand* cand=(PndMvdPidCand*)pointlist->At(i);
      PndMvdPidCand* idealcand=new PndMvdPidCand();


      double dE=0;
      double dx=0;
      double momentum=0;

      for (int k=0;k<cand->GetMvdHits(); k++) {
        dE+=cand->GetMvdHitdE(k);
        dx+=cand->GetMvdHitdx(k);
        momentum+=cand->GetMvdHitMomentum(k);
      }

      if (dx>0)
        momentum/=cand->GetMvdHits();

      PndMvdAdvancedPidAlgo::CalcLikelihood(PndMvdAdvancedPidAlgo::proton, momentum, idealcand);
      h2->Fill(momentum, dE/dx);
      h1->Fill(momentum, cand->GetLikelihood(2212));
      h3->Fill(momentum, idealcand->GetLikelihood(2212));
    }
  }

gROOT Macro

(

"$VMCWORKDIR/gconfig/rootlogon.C"

)

tree SetBranchAddress	(	"PndMvdPidCand"	,
		&	pointlist
	)

timer Start

(

)

timer Stop

(

)

Variable Documentation

c =new TCanvas("c2")

Definition at line 57 of file anaMvdPid.C.

Double_t ctime = timer.CpuTime()

Definition at line 67 of file anaMvdPid.C.

TH2F* h1 =new TH2F("h1", "likelihood count", 100, 0.0, 2.5, 100, 0.0, 1.0)

Definition at line 23 of file anaMvdPid.C.

TH2F* h2 =new TH2F("h2", "energy loss", 100, 0.0, 2.5, 100, 0.0, 0.05)

Definition at line 24 of file anaMvdPid.C.

TH2F* h3 =new TH2F("h3", "ideal likelihood count", 100, 0.0, 2.5, 100, 0.0, 1.0)

Definition at line 25 of file anaMvdPid.C.

TFile* inFile = new TFile("MvdPidIdeal.root","READ")

Definition at line 17 of file anaMvdPid.C.

Int_t nEvents = 10000

Definition at line 29 of file anaMvdPid.C.

TClonesArray* pointlist =new TClonesArray("PndMvdPidCand")

Definition at line 19 of file anaMvdPid.C.

Double_t rtime = timer.RealTime()

Definition at line 66 of file anaMvdPid.C.

TStopwatch timer

Definition at line 13 of file anaMvdPid.C.

TTree* tree = (TTree *)inFile->Get("pndsim")

Definition at line 18 of file anaMvdPid.C.