anaLmdSim.C File Reference

#include "TCanvas.h"

Go to the source code of this file.

Functions
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/rootlogon.C")
	rootlogon ()
gSystem	Load ("libLmd")
timer	Start ()
tree	SetBranchAddress ("LMDPoint",&point)
h1	SetFillColor (kBlue)
h2	SetTitle ("the MC point; #X cm;#Y cm;# Z cm")
h2	SetFillColor (kRed)
h3	SetTitle ("the EnergyLoss: #X keV;")
	for (int i=0;i< nEvents &&i< tree->GetEntriesFast();i++)
h1	DrawCopy ()

Variables
TStopwatch	timer
Int_t	nEvents = 1e5
TFile *	MCPoint = new TFile("/private/huagen/simdata/Lmd_DPM_elastic_6.2_1.9mrad_5M_1.root","READ")
TTree *	tree = (TTree*)MCPoint->Get("pndsim")
TClonesArray *	point = new TClonesArray("PndSdsMCPoint")
TH2F *	h1 = new TH2F("c1","Point on X versus Y",100, -10, 10, 100, -10, 10)
TH3F *	h2 = new TH3F("c2","Point on X,Y and Z", 100, -10, 10, 100, -10, 10, 100, 1070,1150)
TH1F *	h3 = new TH1F("c3","Energy Loss", 1000, 0, 1000)
double	X =0
double	Y =0
double	Z =0
double	eloss
TCanvas *	c = new TCanvas("c2")

Function Documentation

h1 DrawCopy

(

)

for

(

)

Definition at line 36 of file anaLmdSim.C.

References PndSdsMCPoint::GetPosition(), PndSdsMCPoint::GetPositionOut(), i, Y, and Z.

                                                                {
          tree->GetEntry(i);
          for(int j=0;j<point->GetEntriesFast();j++){
                  PndSdsMCPoint* MyPoint = (PndSdsMCPoint*)point->At(j);
                  X = (MyPoint->GetPosition().X()+MyPoint->GetPositionOut().X())/2.;
                  Y = (MyPoint->GetPosition().Y()+MyPoint->GetPositionOut().Y())/2;
                  Z = (MyPoint->GetPosition().Z()+MyPoint->GetPositionOut().Z())/2;
                  eloss = MyPoint->GetEnergyLoss()*1000000;
                 // cout<<"the eloss is:"<<eloss<<endl;
                  h1->Fill(X,Y);
                  h1->SetFillColor(kBlue);
                  h2->Fill(X,Y,Z);
                  h2->SetFillColor(kRed);
                  h3->Fill(eloss);
          }
  }

gSystem Load

(

"libLmd"

)

gROOT LoadMacro

(

"$VMCWORKDIR/gconfig/rootlogon.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);
}

tree SetBranchAddress	(	"LMDPoint"	,
		&	point
	)

h1 SetFillColor

(

kBlue

)

h2 SetFillColor

(

kRed

)

h2 SetTitle

(

"the MC point; #X cm;#Y cm;# Z cm"

)

h3 SetTitle

(

"the EnergyLoss: #X keV;"

)

timer Start

(

)

Variable Documentation

c = new TCanvas("c2")

Definition at line 52 of file anaLmdSim.C.

double eloss

Definition at line 34 of file anaLmdSim.C.

Referenced by PndSttTrackFinderIdeal::DoFind(), PndSttHitProducerIdeal::Exec(), and PndFtsHitProducerIdeal::Exec().

TH2F* h1 = new TH2F("c1","Point on X versus Y",100, -10, 10, 100, -10, 10)

Definition at line 24 of file anaLmdSim.C.

TH3F* h2 = new TH3F("c2","Point on X,Y and Z", 100, -10, 10, 100, -10, 10, 100, 1070,1150)

Definition at line 26 of file anaLmdSim.C.

TH1F* h3 = new TH1F("c3","Energy Loss", 1000, 0, 1000)

Definition at line 29 of file anaLmdSim.C.

TFile* MCPoint = new TFile("/private/huagen/simdata/Lmd_DPM_elastic_6.2_1.9mrad_5M_1.root","READ")

Definition at line 18 of file anaLmdSim.C.

Int_t nEvents = 1e5

Definition at line 16 of file anaLmdSim.C.

TClonesArray* point = new TClonesArray("PndSdsMCPoint")

Definition at line 21 of file anaLmdSim.C.

TStopwatch timer

Definition at line 13 of file anaLmdSim.C.

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

Definition at line 19 of file anaLmdSim.C.

double X =0

Definition at line 33 of file anaLmdSim.C.

double Y =0

Definition at line 33 of file anaLmdSim.C.

double Z =0

Definition at line 33 of file anaLmdSim.C.