Functions
gROOT	Macro ("$VMCWORKDIR/gconfig/rootlogon.C")

tsim	SetBranchAddress ("MCTrack",&mctrack_array)

tsim	SetBranchAddress ("EmcPoint",&point_array)

cout<< "POINTs for new FwEndCap == "<< tsim->	GetEntriesFast ()

c1	Divide (3, 3)

gStyle	SetPalette (1)

c1	cd (1)

hrow_crystal_p	Draw ("colz")

hrow_crystal_p	GetXaxis () -> SetTitle("crystal number")

hrow_crystal_p	GetYaxis () -> SetTitle("row number")

c1	cd (2)

c1	cd (3)

copies	SetTextSize (0.05)

copies	SetFillColor (0)

copies	SetBorderSize (0)

copies	AddEntry (hz_p,"All copies","l")

copies	AddEntry (hz_cp1_p,"copy_1","l")

copies	AddEntry (hz_cp2_p,"copy_2","l")

copies	AddEntry (hz_cp3_p,"copy_3","l")

copies	AddEntry (hz_cp4_p,"copy_4","l")

hz_p	Draw ()

hz_cp1_p	SetLineColor (2)

hz_cp1_p	Draw ("same")

hz_cp2_p	SetLineColor (3)

hz_cp3_p	SetLineColor (4)

hz_cp4_p	SetLineColor (6)

c1	Modified ()

c1	cd (4)

h10d	SetLineColor (1)

c1	cd (5)

c1	cd (6)

c1	cd (7)

c1	cd (8)

c1	cd (9)

c2	Update ()

c3a	Divide (2, 2)

Variables
TFile *	fsim = new TFile("data/fixed/sim_emc_g3_m3_1000eve_fixed.root")

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

TClonesArray *	mctrack_array =new TClonesArray("PndMCTrack")

TClonesArray *	point_array =new TClonesArray("PndEmcPoint")

PndEmcMapper *	emcMap =PndEmcMapper::Instance(2)

Double_t	theta_mc

Double_t	phi_mc

Double_t	ene_mc

Double_t	point_ene

Double_t	point_theta

Double_t	point_phi

Double_t	x_p

Double_t	y_p

Double_t	z_p

Double_t	module_p

Double_t	crystal_p

Double_t	row_p

Double_t	copy_p

Double_t	crystal_pid

Double_t	row_pid

Double_t	hit_ene

Double_t	hit_theta

Double_t	hit_phi

Double_t	hit_id

Double_t	hit_x

Double_t	hit_y

Double_t	z_cp1_p

Double_t	z_cp2_p

Double_t	z_cp3_p

Double_t	z_cp4_p

Int_t	id_p

TVector3	photon_momentum

TLorentzVector	phot

int	ndigi

int	npoint

double	max_energy1 =0

TH2F *	hrow_crystal_p = new TH2F("row_crystal","row vs crystal (POINT)",72,-36,36,72,-36,36)

TH2F *	hx_y_p = new TH2F("x_y_p","x vs y (POINT)",200,-100,100,200,-100,100)

TH1F *	hid_p = new TH1F("id_p","id (POINT)",140000000,300000000,340000000)

TH1F *	hx_p = new TH1F("x_p","X (POINT)",100,-100,100)

TH1F *	hy_p = new TH1F("y_p","Y (POINT)",100,-100,100)

TH1F *	hz_p = new TH1F("z_p","Z (POINT)",100,200,300)

TH1F *	hmodule_p = new TH1F("module_p","Module (POINT)",5,1,6)

TH2F *	hcp_z_p = new TH2F("cp_z_p","Copy vs.Z pos (POINT)",4,1,5,80,200,240)

TH1F *	hz_cp1_p = new TH1F("z_cp1_p","z (POINT)",300,1,300)

TH1F *	hz_cp2_p = new TH1F("z_cp2_p","z (POINT)",300,1,300)

TH1F *	hz_cp3_p = new TH1F("z_cp3_p","z (POINT)",300,1,300)

TH1F *	hz_cp4_p = new TH1F("z_cp4_p","z (POINT)",300,1,300)

TH2F *	hxpad_x_p = new TH2F("xpad_x_p","x vs xpad (POINT)",200,-100,100,200,-100,100)

TH2F *	hypad_y_p = new TH2F("ypad_y_p","y vs ypad (POINT)",200,-100,100,200,-100,100)

TH2F *	hxpad_x_cp1_p = new TH2F("xpad_x_cp1_p","x vs xpad _cp1_ (POINT)",200,-100,100,200,-100,100)

TH2F *	hxpad_x_cp2_p = new TH2F("xpad_x_cp2_p","x vs xpad _cp2_ (POINT)",200,-100,100,200,-100,100)

TH2F *	hxpad_x_cp3_p = new TH2F("xpad_x_cp3_p","x vs xpad _cp3_ (POINT)",200,-100,100,200,-100,100)

TH2F *	hxpad_x_cp4_p = new TH2F("xpad_x_cp4_p","x vs xpad _cp4_ (POINT)",200,-100,100,200,-100,100)

TH2F *	hypad_y_cp1_p = new TH2F("ypad_y_cp1_p","y vs ypad _cp1_ (POINT)",200,-100,100,200,-100,100)

TH2F *	hypad_y_cp2_p = new TH2F("ypad_y_cp2_p","y vs ypad _cp2_ (POINT)",200,-100,100,200,-100,100)

TH2F *	hypad_y_cp3_p = new TH2F("ypad_y_cp3_p","y vs ypad _cp3_ (POINT)",200,-100,100,200,-100,100)

TH2F *	hypad_y_cp4_p = new TH2F("ypad_y_cp4_p","y vs ypad _cp4_ (POINT)",200,-100,100,200,-100,100)

TH1F *	hene_mc = new TH1F("hene_mc","MC energy (GeV)",100,0.0,1.05)

TH1F *	hene_p = new TH1F("hene_p","DATA energy (GeV)",1000,0.0,0.1)

TH1F *	hpoi_ene = new TH1F("hpoi_ene","SUM of Points energy (GeV)",100,0.0,1.05)

TH1F *	hhit_ene = new TH1F("hhit_ene","SUM of Hits energy (GeV)",100,0.0,1.05)

TH1F *	h10 = new TH1F("h10","POINT Theta (MC)",180,0.,180.)

TH1F *	h20 = new TH1F("h20","POINT Phi (MC)",360,-180.,180.)

TH1F *	h10d = new TH1F("h10d","POINT Theta (data)",180,0.,180.)

TH1F *	h20d = new TH1F("h20d","POINT Phi (data)",360,-180.,180.)

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

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

TH1F *	h1a = new TH1F("h1a","POINT Theta difference -changed id ",100,-5.,5.)

TH1F *	h2a = new TH1F("h2a","POINT Phi difference - changed id",100,-5.,5.)

TH1F *	h3 = new TH1F("h3","HIT Theta difference",100,-5.,5.)

TH1F *	h4 = new TH1F("h4","HIT Phi difference",100,-5.,5.)

TH1F *	h3a = new TH1F("h3a","HIT Theta difference -changed id",100,-5.,5.)

TH1F *	h4a = new TH1F("h4a","HIT Phi difference -changed id",100,-5.,5.)

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

TLegend *	copies = new TLegend(0.60,0.50,0.85,0.85)

TCanvas *	c2 = new TCanvas("c2", "", 100, 100, 800, 800)

TCanvas *	c3a = new TCanvas("c3a", "", 100, 100, 800, 800)

TCanvas *	c3b = new TCanvas("c3b", "", 100, 100, 800, 800)

Function Documentation

copies AddEntry	(	hz_p	,
		"All copies"	,
		"l"
	)

copies AddEntry	(	hz_cp1_p	,
		"copy_1"	,
		"l"
	)

copies AddEntry	(	hz_cp2_p	,
		"copy_2"	,
		"l"
	)

copies AddEntry	(	hz_cp3_p	,
		"copy_3"	,
		"l"
	)

copies AddEntry	(	hz_cp4_p	,
		"copy_4"	,
		"l"
	)

c1 cd ( 1 )

c1 cd ( 2 )

c1 cd ( 3 )

c1 cd ( 4 )

c1 cd ( 5 )

c1 cd ( 6 )

myCan4 cd ( 7 )

myCan4 cd ( 8 )

c1 cd ( 9 )

c2 Divide	(	3	,
		3
	)

c3a Divide	(	2	,
		2
	)

h5 Draw ( "colz" )

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

hz_cp1_p Draw ( "same" )

cout<< "POINTs for new FwEndCap == " << tsim-> GetEntriesFast ( )

Definition at line 90 of file anal_point_fwendcap.C.

References PndEmcPoint::GetCopy(), PndEmcPoint::GetCrystal(), PndEmcPoint::GetModule(), PndMCTrack::GetMomentum(), PndEmcPoint::GetRow(), PndEmcPoint::GetXPad(), PndEmcPoint::GetYPad(), i, mctrack, npoints, Pi, point, point_ene, x_p, and y_p.

Referenced by anaRadLength(), PndPidCorrelator::ConstructChargedCandidate(), PndSttTrackFinderReal::DoFind(), PndSttTrackFinderIdeal::DoFind(), PndKFParticleFinderQA::Exec(), PndFtsDataAccessor::Exec(), PndSolCorrTask::Exec(), PndFtsCATracking::Exec(), PndFtsTrackerIdeal::Exec(), PndCATracking::Exec(), PndSimpleAnalysis::FillGenericLists(), PndMvdSttGemRiemannTrackFinderTask::FillHitArray(), PndMvdRiemannTrackFinderTask::FillHitArray(), PndAnaTask::FillList(), PndBranchBurstBuilder_timeCut_STB::GetBurstData(), PndTrkNeighboringMap::GetCandseeds(), PndSttMatchTracks::GetHitFromCollections(), PndSttTrackFinderIdeal::GetHitFromCollections(), PndSttHelixTrackFitter::GetHitFromCollections(), PndSttTrackFinderReal::GetHitFromCollections(), PndTrkNeighboringMap::GetIndivisiblesToHit(), PndTrkNeighboringMap::GetNeighboringsToHit(), PndSttMatchTracks::GetPointFromCollections(), PndSttTrackFinderIdeal::GetPointFromCollections(), PndSttTrackFinderReal::GetPointFromCollections(), PndDrcAccuDigiPixelDraw::ReadAllHits(), and GFTrack::setCardinalRep().

                                   {
     //for (Int_t j=1;j< 20; j++) {
     Float_t sum_point_ene=0;
     tsim->GetEntry(j);
 
     PndMCTrack *mctrack=(PndMCTrack *) mctrack_array->At(0); // 1st element in mctrack_array
     photon_momentum=mctrack->GetMomentum();
     theta_mc=photon_momentum.Theta()*(180./TMath::Pi());
     phi_mc=photon_momentum.Phi()*(180./TMath::Pi());//+360.;
     //phot=mctrack->Get4Momentum();
 
     Int_t npoints = point_array->GetEntries();
     cout << "POINT array ==  " << npoints<< endl;
     for (Int_t i=0; i<npoints; i++)
       {
         PndEmcPoint *point=(PndEmcPoint*)point_array->At(i);
 
         module_p = point->GetModule();
 
         if (module_p==3){
           //hene_mc->Fill(phot.E());
 
           point_ene=point->GetEnergyLoss();
           hene_p->Fill(point_ene);
 
           sum_point_ene+=point_ene;
 
           TVector3 point_pos(point->GetX(),point->GetY(),point->GetZ());
           point_theta=point_pos.Theta()*(180./TMath::Pi());
           point_phi=point_pos.Phi()*(180./TMath::Pi());//+360.;
 
 
           h10->Fill(theta_mc);
           h20->Fill(phi_mc);
 
           h10d->Fill(point_theta);
           h20d->Fill(point_phi);
 
           h1->Fill(point_theta-theta_mc);
           h2->Fill(point_phi-phi_mc);
 
           crystal_p = point->GetCrystal();
           row_p     = point->GetRow();
           copy_p    = point->GetCopy();
 
           x_p = point->GetX();
           y_p = point->GetY();
           z_p = point->GetZ();
 
           hz_p->Fill(z_p);
           if (copy_p==1) hz_cp1_p->Fill(z_p);
           if (copy_p==2) hz_cp2_p->Fill(z_p);
           if (copy_p==3) hz_cp3_p->Fill(z_p);
           if (copy_p==4) hz_cp4_p->Fill(z_p);
 
           hxpad_x_p->Fill(point->GetXPad(),x_p);
           hypad_y_p->Fill(point->GetYPad(),y_p);
 
           if (copy_p==1){
             hxpad_x_cp1_p->Fill(point->GetXPad(),x_p);
             hypad_y_cp1_p->Fill(point->GetYPad(),y_p);
           }
           if (copy_p==2){
             hxpad_x_cp2_p->Fill(point->GetXPad(),x_p);
             hypad_y_cp2_p->Fill(point->GetYPad(),y_p);
           }
           if (copy_p==3){
             hxpad_x_cp3_p->Fill(point->GetXPad(),x_p);
             hypad_y_cp3_p->Fill(point->GetYPad(),y_p);
           }
           if (copy_p==4){
             hxpad_x_cp4_p->Fill(point->GetXPad(),x_p);
             hypad_y_cp4_p->Fill(point->GetYPad(),y_p);
           }
 
           id_p = point->GetDetectorID();
 
           crystal_pid = (id_p%10000);
           row_pid     = ((id_p/1000000)%100);
 
           hmodule_p->Fill(module_p);
           hcp_z_p->Fill(copy_p,z_p);
           hrow_crystal_p->Fill(crystal_pid,row_pid);
           hx_p->Fill(x_p);
           hy_p->Fill(y_p);
           hx_y_p->Fill(x_p,y_p);
           hid_p->Fill(id_p);
         }
       }
     hpoi_ene->Fill(sum_point_ene);
   }