Functions
t	SetBranchAddress ("EmcCluster",&cluster_array)

t	AddFriend ("pndsim","digi_emc.root")

t	SetBranchAddress ("EmcDigi",&digi_array)

tsim	SetBranchAddress ("MCTrack",&mctrack_array)

	if ((cluster->NumberOfDigis()>1)&&(cluster_energy >0.02))

h3	SetTitle ("Cluster energy of 1 GeV photon")

h3	GetXaxis () -> SetTitle("Energy, GeV")

h3	Draw ()

h1	SetTitle ("Difference between cluster and initial photon theta angle")

h2	SetTitle ("Difference between cluster and initial photon phi angle")

h2theta	GetYaxis () -> SetTitle("#theta_

c6	Divide (2, 2)

c6	cd (1)

c6	cd (2)

c6	cd (3)

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

TClonesArray *	cluster_array =new TClonesArray("PndEmcCluster")

TClonesArray *	digi_array =new TClonesArray("PndEmcDigi")

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

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

TClonesArray *	mctrack_array =new TClonesArray("PndMCTrack")

TVector3	photon_momentum

double	cluster_energy =cluster->energy()

double	cluster_theta

double	cluster_phi

double	theta

double	phi

double	theta_diff

double	phi_diff

int	ndigi

int	npoint

double	max_energy =0

TH1F *	h1 = new TH1F("h1","Theta difference",200,-5.,5.)

TH1F *	h2 = new TH1F("h2","Phi difference",200,-5.,5.)

TH1F *	h3 = new TH1F("h3","Cluster energy",100,0.85,1.05)

TH2F *	h2theta = new TH2F("h2theta","Theta difference",200,0.,180.,200,-5.,5.)

TH2F *	h2phi = new TH2F("h2phi","Phi difference",200,0.,180.,200,-5.,5.)

TH1F *	hE1 = new TH1F("hE1","E1",200,0.,1.05)

TH1F *	hE1E9 = new TH1F("hE1E9","E1 / E9",200,0.,1.05)

TH1F *	hE9E25 = new TH1F("hE9E25","E9 / E25",200,0.,1.05)

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

TCanvas *	c2 = new TCanvas("c2", "#theta_{reco} - #theta_{truth}", 100, 100, 800, 800)

TCanvas *	c3 = new TCanvas("c3", "#phi_{reco} - #phi_{truth}", 100, 100, 800, 800)

TCanvas *	c4 = new TCanvas("c4", "#theta_{reco} - #theta_{truth} vs #theta_{truth}", 100, 100, 800, 800)

	degree

TCanvas *	c5 = new TCanvas("c5", "#phi_{reco} - #phi_{truth} vs #phi_{truth}", 100, 100, 800, 800)

TCanvas *	c6 = new TCanvas("c6", "Cluster Properties", 100, 100, 800, 800)

Function Documentation

t AddFriend	(	"pndsim"	,
		"digi_emc.root"
	)

c6 cd ( 1 )

c6 cd ( 2 )

c6 cd ( 3 )

c6 Divide	(	2	,
		2
	)

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

h3 GetXaxis ( ) -> SetTitle("Energy, GeV")

Definition at line 653 of file anal_hit_digi_cluster_fwendcap.C.

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

h2theta GetYaxis ( ) -> SetTitle("#theta_

Definition at line 131 of file macro/detectors/emc/dedicated/reco_analys.C.

131 {reco} - #theta_{truth}, degree");

degree

Definition: macro/detectors/emc/dedicated/reco_analys.C:130

reco

int reco()

Definition: detectors/disc/reco.C:1

if ( (cluster->NumberOfDigis()>1)&&(cluster_energy >0.02) )

Definition at line 49 of file macro/detectors/emc/dedicated/reco_analys.C.

         {
                 t->GetEntry(j);
                 tsim->GetEntry(j);
 
                 PndMCTrack *mctrack=(PndMCTrack *) mctrack_array->At(0);
                 photon_momentum=mctrack->GetMomentum();
                 theta=photon_momentum.Theta();
                 phi=photon_momentum.Phi();
 
 
                 // Loop over clusters
                 // If we have 1 initial particle and several cluster
                 // we can separate cluster from the first interaction by maximum energy
 
                 max_energy=0;
 
                 for (Int_t i=0; i<cluster_array->GetEntriesFast(); i++)
                 {
                         PndEmcCluster *cluster=(PndEmcCluster*)cluster_array->At(i);
                         cluster_energy=cluster->energy();
                         if (cluster_energy>max_energy)
                         {
                                 max_energy=cluster_energy;
                                 TVector3 cluster_pos=cluster->where();
                                 cluster_theta=cluster_pos.Theta();
                                 cluster_phi=cluster_pos.Phi();
                         }
 
                 }
 
                 if (max_energy>0.6)
                 {
                         theta_diff=(cluster_theta-theta)*180./TMath::Pi();
                         h1->Fill(theta_diff);
                         h2theta->Fill(theta*TMath::RadToDeg(),theta_diff);
 
                         phi_diff=(cluster_phi-phi)*180./TMath::Pi();
                         h2->Fill(phi_diff);
                         h2phi->Fill(phi*TMath::RadToDeg(),phi_diff);
                 }
 
         }