anal_hit_digi_cluster_fwendcap.C File Reference

Go to the source code of this file.

Functions
gROOT	Macro ("$VMCWORKDIR/gconfig/rootlogon.C")
tsim	SetBranchAddress ("MCTrack",&mctrack_array)
tful	SetBranchAddress ("EmcHit",&hit_array)
tful	SetBranchAddress ("EmcDigi",&digi_array)
tful	SetBranchAddress ("EmcCluster",&cluster_array)
	for (int j=0;j< ncounts;j++) = CAMath::Finite( t.Par(i) )
c1	Divide (2, 3)
gStyle	SetPalette (1)
c1	cd (1)
copies	SetTextSize (0.05)
copies	SetFillColor (0)
copies	SetBorderSize (0)
copies	AddEntry (hz_h,"All copies (G3)","l")
copies	AddEntry (hz_cp1_h,"copy_1","l")
copies	AddEntry (hz_cp2_h,"copy_2","l")
copies	AddEntry (hz_cp3_h,"copy_3","l")
copies	AddEntry (hz_cp4_h,"copy_4","l")
hz_h	Draw ()
hz_cp1_h	SetLineColor (2)
hz_cp1_h	Draw ("same")
hz_cp2_h	SetLineColor (3)
hz_cp3_h	SetLineColor (4)
hz_cp4_h	SetLineColor (6)
hz_h	GetXaxis () -> SetTitle("Z position")
hz_h	GetYaxis () -> SetTitle("Counts")
c1	Modified ()
c1	cd (2)
copies4	AddEntry (hhit_ene,"All copies (G3)","l")
copies4	AddEntry (hhit_ene_cp1,"copy_1","l")
copies4	AddEntry (hhit_ene_cp2,"copy_2","l")
copies4	AddEntry (hhit_ene_cp3,"copy_3","l")
copies4	AddEntry (hhit_ene_cp4,"copy_4","l")
c1	cd (3)
copies6	AddEntry (hz_d,"All copies (G3)","l")
copies6	AddEntry (hz_cp1_d,"copy_1","l")
copies6	AddEntry (hz_cp2_d,"copy_2","l")
copies6	AddEntry (hz_cp3_d,"copy_3","l")
copies6	AddEntry (hz_cp4_d,"copy_4","l")
c1	cd (4)
copies5	AddEntry (hdigi_ene,"All copies (G3)","l")
copies5	AddEntry (hdigi_ene_cp1,"copy_1","l")
copies5	AddEntry (hdigi_ene_cp2,"copy_2","l")
copies5	AddEntry (hdigi_ene_cp3,"copy_3","l")
copies5	AddEntry (hdigi_ene_cp4,"copy_4","l")
c1	cd (5)
copies3	AddEntry (hz_clust,"All copies (G3)","l")
copies3	AddEntry (hz_cp1_c,"copy_1","l")
copies3	AddEntry (hz_cp2_c,"copy_2","l")
copies3	AddEntry (hz_cp3_c,"copy_3","l")
copies3	AddEntry (hz_cp4_c,"copy_4","l")
c1	cd (6)
copies2a	AddEntry (hclus_ene,"All copies (G3)","l")
copies2a	AddEntry (hclus_ene_cp1,"copy_1","l")
copies2a	AddEntry (hclus_ene_cp2,"copy_2","l")
copies2a	AddEntry (hclus_ene_cp3,"copy_3","l")
copies2a	AddEntry (hclus_ene_cp4,"copy_4","l")
c9	Divide (2, 2)

Variables
TFile *	fsim = new TFile("data/fixed/sim_emc_g4_m3_1000eve_fixed.root")
TTree *	tsim =(TTree *) fsim->Get("pndsim")
TClonesArray *	mctrack_array =new TClonesArray("PndMCTrack")
PndEmcMapper *	emcMap =PndEmcMapper::Instance(1)
TFile *	fful = new TFile("data/fixed/full_emc_g4_m3_1000eve_fixed.root")
TTree *	tful =(TTree *) fful->Get("pndsim")
TClonesArray *	hit_array =new TClonesArray("PndEmcHit")
TClonesArray *	digi_array =new TClonesArray("PndEmcDigi")
TClonesArray *	cluster_array =new TClonesArray("PndEmcCluster")
Double_t	theta_mc
Double_t	phi_mc
Double_t	ene_mc
Double_t	cemc
Double_t	cluster_energy
Double_t	cluster_theta
Double_t	cluster_phi
Double_t	hit_ene
Double_t	hit_theta
Double_t	hit_phi
Double_t	x_h
Double_t	y_h
Double_t	z_h
Double_t	module_h
Double_t	crystal_h
Double_t	row_h
Double_t	copy_h
Double_t	crystal_hid
Double_t	row_hid
Double_t	hit_id
Double_t	hit_x
Double_t	hit_y
Double_t	crystal_cid
Double_t	row_cid
Double_t	crystal_c
Double_t	row_c
Double_t	copy_c
Double_t	digi_ene
Double_t	digi_theta
Double_t	digi_phi
Double_t	module_d
Double_t	crystal_d
Double_t	row_d
Double_t	copy_d
Double_t	crystal_did
Double_t	row_did
Double_t	d_id
Int_t	id_h
Int_t	id_d
TVector3	photon_momentum
TLorentzVector	p4mom
int	ndigi
int	nhit
double	max_energy1 =0
TH2F *	hrow_crystal_h = new TH2F("row_crystal","row vs crystal (HIT)",72,-36,36,72,-36,36)
TH2F *	hrow_crystal_d = new TH2F("row_crystal_digi","row vs crystal (DIGI)",72,-36,36,72,-36,36)
TH2F *	hrow_crystal_c = new TH2F("row_crystal_c","row vs crystal (CLUSTER)",72,-36,36,72,-36,36)
TH2F *	hx_y_h = new TH2F("x_y_h","x vs y (HIT)",200,-100,100,200,-100,100)
TH2F *	hx_y_d = new TH2F("x_y_d","GetThetaInt() vs GetPhiInt() (DIGI)",200,200,300,200,200,300)
TH1F *	hid_h = new TH1F("id_h","id (HIT)",140000000,300000000,340000000)
TH1F *	hz_h = new TH1F("z_h","z (HIT)",300,1,300)
TH1F *	hmodule_h = new TH1F("module_h","Module (HIT)",5,1,6)
TH2F *	hxpad_x_h = new TH2F("xpad_x_h","xpad vs x (HIT)",200,-100,100,200,-100,100)
TH2F *	hypad_y_h = new TH2F("ypad_y_h","ypad vs y (HIT)",200,-100,100,200,-100,100)
TH2F *	hxpad_x_d = new TH2F("xpad_x_d","xpad vs GetThetaInt() (DIGI)",200,-100,100,200,200,300)
TH2F *	hypad_y_d = new TH2F("ypad_y_d","ypad vs GetPhiInt() (DIGI)",200,-100,100,200,200,300)
TH2F *	hx_row_h = new TH2F("x_row_h","x vs row (HIT)",200,-100,100,72,-36,36)
TH2F *	hx_row_d = new TH2F("x_row_d","GetThetaInt() vs row (DIGI)",200,200,300,72,-36,36)
TH2F *	hcp_z_h = new TH2F("cp_z_h","Copy vs.Z pos (HITS)",4,1,5,80,200,240)
TH1F *	hz_cp1_h = new TH1F("z_cp1_h","z (HIT)",300,1,300)
TH1F *	hz_cp2_h = new TH1F("z_cp2_h","z (HIT)",300,1,300)
TH1F *	hz_cp3_h = new TH1F("z_cp3_h","z (HIT)",300,1,300)
TH1F *	hz_cp4_h = new TH1F("z_cp4_h","z (HIT)",300,1,300)
TH1F *	hz_d = new TH1F("z_d","z (DIGI)",300,1,300)
TH1F *	hz_cp1_d = new TH1F("z_cp1_d","z (DIGI)",300,1,300)
TH1F *	hz_cp2_d = new TH1F("z_cp2_d","z (DIGI)",300,1,300)
TH1F *	hz_cp3_d = new TH1F("z_cp3_d","z (DIGI)",300,1,300)
TH1F *	hz_cp4_d = new TH1F("z_cp4_d","z (DIGI)",300,1,300)
TH1F *	hz_clust = new TH1F("z_cluster","z (CLUSTER)",300,1,300)
TH1F *	hz_cp1_c = new TH1F("z_cp1_c","z (CLUSTER)",300,1,300)
TH1F *	hz_cp2_c = new TH1F("z_cp2_c","z (CLUSTER)",300,1,300)
TH1F *	hz_cp3_c = new TH1F("z_cp3_c","z (CLUSTER)",300,1,300)
TH1F *	hz_cp4_c = new TH1F("z_cp4_c","z (CLUSTER)",300,1,300)
TH1F *	hene_mc = new TH1F("hene_mc","MC energy HIT (GeV)",100,0.05,1.05)
TH1F *	dene_mc = new TH1F("dene_mc","MC energy DIGI (GeV)",100,0.05,1.05)
TH1F *	hclus_ene_sum_mc = new TH1F("hclus_ene_sum_mc","SUM of MC Cluster energy (GeV)",100,0.05,1.05)
TH1F *	hene_h = new TH1F("hene_h","DATA energy HIT (GeV)",100,0.05,1.05)
TH1F *	hene_d = new TH1F("hene_d","DATA energy DIGI (GeV)",100,0.05,1.05)
TH1F *	hhit_ene = new TH1F("hhit_ene","SUM of Hits energy (GeV)",100,0.05,1.05)
TH1F *	hdigi_ene = new TH1F("hdigi_ene","SUM of Digi energy (GeV)",100,0.05,1.05)
TH1F *	hclus_ene = new TH1F("hclus_ene","SUM of Cluster energy (GeV)",100,0.05,1.05)
TH1F *	hclus_ene_mc = new TH1F("hclus_ene_mc","SUM of MC Cluster energy (GeV)",100,0.05,1.05)
TH1F *	hclus_ene_sum = new TH1F("hclus_ene_sum","SUM of Cluster energy (GeV)",100,0.05,1.05)
TH1F *	hhit_ene_cp1 = new TH1F("hhit_ene_cp1","SUM of Hits energy (GeV) _cp1",100,0.05,1.05)
TH1F *	hhit_ene_cp2 = new TH1F("hhit_ene_cp2","SUM of Hits energy (GeV) _cp2",100,0.05,1.05)
TH1F *	hhit_ene_cp3 = new TH1F("hhit_ene_cp3","SUM of Hits energy (GeV) _cp3",100,0.05,1.05)
TH1F *	hhit_ene_cp4 = new TH1F("hhit_ene_cp4","SUM of Hits energy (GeV) _cp4",100,0.05,1.05)
TH1F *	hclus_ene_cp1 = new TH1F("hclus_ene_cp1","Cluster energy (GeV) _cp1",100,0.05,1.05)
TH1F *	hclus_ene_cp2 = new TH1F("hclus_ene_cp2","Cluster energy (GeV) _cp2",100,0.05,1.05)
TH1F *	hclus_ene_cp3 = new TH1F("hclus_ene_cp3","Cluster energy (GeV) _cp3",100,0.05,1.05)
TH1F *	hclus_ene_cp4 = new TH1F("hclus_ene_cp4","Cluster energy (GeV) _cp4",100,0.05,1.05)
TH1F *	h10 = new TH1F("h10","HIT Theta (MC)",180,0.,180.)
TH1F *	h20 = new TH1F("h20","HIT Phi (MC)",360,-180.,180.)
TH1F *	h10d = new TH1F("h10d","HIT Theta (data)",180,0.,180.)
TH1F *	h20d = new TH1F("h20d","HIT Phi (data)",360,-180.,180.)
TH1F *	h1 = new TH1F("h1","HIT Theta difference",100,-5.,5.)
TH1F *	h2 = new TH1F("h2","HIT Phi difference",100,-5.,5.)
TH1F *	h10digi = new TH1F("h10digi","DIGI Theta (MC)",180,0.,180.)
TH1F *	h20digi = new TH1F("h20digi","DIGI Phi (MC)",360,-180.,180.)
TH1F *	h10ddigi = new TH1F("h10ddigi","DIGI Theta (data)",180,0.,180.)
TH1F *	h20ddigi = new TH1F("h20ddigi","DIGI Phi (data)",360,-180.,180.)
TH1F *	h1d = new TH1F("h1d","DIGI Theta difference",100,-5.,5.)
TH1F *	h2d = new TH1F("h2d","DIGI Phi difference",100,-5.,5.)
TH1F *	hdigi_ene_cp1 = new TH1F("hdigi_ene_cp1","SUM of Digi energy (GeV) _cp1",100,0.05,1.05)
TH1F *	hdigi_ene_cp2 = new TH1F("hdigi_ene_cp2","SUM of Digi energy (GeV) _cp2",100,0.05,1.05)
TH1F *	hdigi_ene_cp3 = new TH1F("hdigi_ene_cp3","SUM of Digi energy (GeV) _cp3",100,0.05,1.05)
TH1F *	hdigi_ene_cp4 = new TH1F("hdigi_ene_cp4","SUM of Digi energy (GeV) _cp4",100,0.05,1.05)
TH1F *	hthc_mc = new TH1F("hthcmc","CLUSTER Theta (MC)",180,0.,180.)
TH1F *	hphc_mc = new TH1F("hphcmc","CLUSTER Phi (MC)",360,-180.,180.)
TH1F *	hthc = new TH1F("hthc","CLUSTER Theta (data)",180,0.,180.)
TH1F *	hphc = new TH1F("hphc","CLUSTER Phi (data)",360,-180.,180.)
TH1F *	hthc_diff = new TH1F("hthc_diff","CLUSTER Theta difference",100,-5.,5.)
TH1F *	hphc_diff = new TH1F("hphc_diff","CLUSTER Phi difference",100,-5.,5.)
TH1F *	hid_d = new TH1F("id_d","id (DIGI)",140000000,300000000,340000000)
TH1F *	hmodule_d = new TH1F("module_d","Module (DIGI)",5,1,6)
TH1F *	hClustNoOrig = new TH1F("hClustNoOrig","hClustNoOrig",10,0,9)
TH1F *	hClustNo = new TH1F("hClustNo","hClustNo",10,0,9)
TH1F *	hClNo_cp1 = new TH1F("hClNo_cp1","hClNo_cp1",10,0,9)
TH1F *	hClNo_cp2 = new TH1F("hClNo_cp2","hClNo_cp2",10,0,9)
TH1F *	hClNo_cp3 = new TH1F("hClNo_cp3","hClNo_cp3",10,0,9)
TH1F *	hClNo_cp4 = new TH1F("hClNo_cp4","hClNo_cp4",10,0,9)
Int_t	ncounts = tsim->GetEntries()
TCanvas *	c1 = new TCanvas("c1", "", 100, 100, 800, 800)
TLegend *	copies = new TLegend(0.65,0.4,0.85,0.6)
TLegend *	copies4 = new TLegend(0.20,0.40,0.5,0.6)
TLegend *	copies6 = new TLegend(0.65,0.4,0.85,0.6)
TLegend *	copies5 = new TLegend(0.20,0.40,0.5,0.6)
TLegend *	copies3 = new TLegend(0.20,0.30,0.5,0.5)
TLegend *	copies2a = new TLegend(0.20,0.4,0.5,0.6)
TCanvas *	c9 = new TCanvas("c9", "", 100, 100, 800, 800)

Function Documentation

copies AddEntry	(	hz_h	,
		"All copies (G3)"	,
		"l"
	)

Referenced by Eff_lE(), and Pic_Compare_FEP_Eff().

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

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

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

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

copies4 AddEntry	(	hhit_ene	,
		"All copies (G3)"	,
		"l"
	)

copies4 AddEntry	(	hhit_ene_cp1	,
		"copy_1"	,
		"l"
	)

copies4 AddEntry	(	hhit_ene_cp2	,
		"copy_2"	,
		"l"
	)

copies4 AddEntry	(	hhit_ene_cp3	,
		"copy_3"	,
		"l"
	)

copies4 AddEntry	(	hhit_ene_cp4	,
		"copy_4"	,
		"l"
	)

copies6 AddEntry	(	hz_d	,
		"All copies (G3)"	,
		"l"
	)

copies6 AddEntry	(	hz_cp1_d	,
		"copy_1"	,
		"l"
	)

copies6 AddEntry	(	hz_cp2_d	,
		"copy_2"	,
		"l"
	)

copies6 AddEntry	(	hz_cp3_d	,
		"copy_3"	,
		"l"
	)

copies6 AddEntry	(	hz_cp4_d	,
		"copy_4"	,
		"l"
	)

copies5 AddEntry	(	hdigi_ene	,
		"All copies (G3)"	,
		"l"
	)

copies5 AddEntry	(	hdigi_ene_cp1	,
		"copy_1"	,
		"l"
	)

copies5 AddEntry	(	hdigi_ene_cp2	,
		"copy_2"	,
		"l"
	)

copies5 AddEntry	(	hdigi_ene_cp3	,
		"copy_3"	,
		"l"
	)

copies5 AddEntry	(	hdigi_ene_cp4	,
		"copy_4"	,
		"l"
	)

copies3 AddEntry	(	hz_clust	,
		"All copies (G3)"	,
		"l"
	)

copies3 AddEntry	(	hz_cp1_c	,
		"copy_1"	,
		"l"
	)

copies3 AddEntry	(	hz_cp2_c	,
		"copy_2"	,
		"l"
	)

copies3 AddEntry	(	hz_cp3_c	,
		"copy_3"	,
		"l"
	)

copies3 AddEntry	(	hz_cp4_c	,
		"copy_4"	,
		"l"
	)

copies2a AddEntry	(	hclus_ene	,
		"All copies (G3)"	,
		"l"
	)

copies2a AddEntry	(	hclus_ene_cp1	,
		"copy_1"	,
		"l"
	)

copies2a AddEntry	(	hclus_ene_cp2	,
		"copy_2"	,
		"l"
	)

copies2a AddEntry	(	hclus_ene_cp3	,
		"copy_3"	,
		"l"
	)

copies2a AddEntry	(	hclus_ene_cp4	,
		"copy_4"	,
		"l"
	)

c1 cd

(

1

)

c1 cd

(

2

)

c1 cd

(

3

)

c1 cd

(

4

)

c cd

(

5

)

c cd

(

6

)

c1 Divide	(	2	,
		3
	)

c9 Divide	(	2	,
		2
	)

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

hMcGemPhi Draw

(

"same"

)

for

(

)

= CAMath::Finite( t.Par(i) )

Definition at line 156 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by PndFilteredPrimaryGenerator::AddFilter(), operator<<(), PndFsmCmpDet::readParameters(), and PndFsmCombiDet::readParameters().

                                   {
   //for (Int_t j=1;j< 3; j++){
     Float_t sum_hit_ene,hit_ene_cp1,hit_ene_cp2,hit_ene_cp3,hit_ene_cp4;
     Float_t sum_digi_ene,digi_ene_cp1,digi_ene_cp2,digi_ene_cp3,digi_ene_cp4;
     Float_t sum_clus_ene=0;
     Float_t cemc_sum=0;

     sum_hit_ene = hit_ene_cp1 = hit_ene_cp2 = hit_ene_cp3 = hit_ene_cp4 = 0;
     sum_digi_ene = digi_ene_cp1 = digi_ene_cp2 = digi_ene_cp3 = digi_ene_cp4 =0;

     tsim->GetEntry(j);
     tful->GetEntry(j);

        PndMCTrack *mctrack=(PndMCTrack *) mctrack_array->At(0);
        photon_momentum=mctrack->GetMomentum();
        theta_mc=photon_momentum.Theta()*(180./TMath::Pi());
        phi_mc=photon_momentum.Phi()*(180./TMath::Pi());

        p4mom=mctrack->Get4Momentum();

        Int_t nhits = hit_array->GetEntries();
        for (Int_t i=0; i<nhits; i++)
          {
            PndEmcHit *hit=(PndEmcHit*)hit_array->At(i);

            module_h = hit->GetModule();

            if (module_h==3){

              hene_mc->Fill(p4mom.E());

              hit_ene=hit->GetEnergy();

              hene_h->Fill(hit_ene);

              sum_hit_ene+=hit_ene;

              crystal_h = hit->GetCrystal();
              row_h = hit->GetRow();
              copy_h = hit->GetCopy();

              if (copy_h==1) hit_ene_cp1+=hit->GetEnergy();
              if (copy_h==2) hit_ene_cp2+=hit->GetEnergy();
              if (copy_h==3) hit_ene_cp3+=hit->GetEnergy();
              if (copy_h==4) hit_ene_cp4+=hit->GetEnergy();

              TVector3 hit_pos(hit->GetX(),hit->GetY(),hit->GetZ());
              hit_theta=hit_pos.Theta()*(180./TMath::Pi());
              hit_phi=hit_pos.Phi()*(180./TMath::Pi());

              //cout << "MC: theta= "<<theta_mc<<", phi= "<<phi_mc<<endl;
              h10->Fill(theta_mc);
              h20->Fill(phi_mc);

              h10d->Fill(hit_theta);
              h20d->Fill(hit_phi);

              h1->Fill(hit_theta-theta_mc);
              h2->Fill(hit_phi-phi_mc);

              if (copy_h==1) hz_cp1_h->Fill(hit->GetZ());
              if (copy_h==2) hz_cp2_h->Fill(hit->GetZ());
              if (copy_h==3) hz_cp3_h->Fill(hit->GetZ());
              if (copy_h==4) hz_cp4_h->Fill(hit->GetZ());

              hz_h->Fill(hit->GetZ());
              hcp_z_h->Fill(copy_h,hit->GetZ());

              hxpad_x_h->Fill(hit->GetXPad(),hit->GetX());
              hypad_y_h->Fill(hit->GetYPad(),hit->GetY());

              hx_row_h->Fill(hit->GetRow(),hit->GetX());

              id_h = hit->GetDetectorID();

              crystal_hid = (id_h%10000);
              row_hid = ((id_h/1000000)%100);

              hmodule_h->Fill(module_h);
              hrow_crystal_h->Fill(crystal_hid,row_hid);
              hx_y_h->Fill(hit->GetX(),hit->GetY());

              hid_h->Fill(id_h);
            }
          }
        hhit_ene->Fill(sum_hit_ene);
        hhit_ene_cp1->Fill(hit_ene_cp1);
        hhit_ene_cp2->Fill(hit_ene_cp2);
        hhit_ene_cp3->Fill(hit_ene_cp3);
        hhit_ene_cp4->Fill(hit_ene_cp4);

        // digits
        Int_t ndigits = digi_array->GetEntries();
        //cout << "DIGI array ==  " << ndigits<< endl;
        for (Int_t i=0; i<ndigits; i++)
          {
            PndEmcDigi *digi=(PndEmcDigi*)digi_array->At(i);

            TVector3 digi_pos=digi->where();

            module_d = digi->GetModule();

            if (module_d==3){

              dene_mc->Fill(p4mom.E());

              digi_ene=digi->GetEnergy();
              hene_d->Fill(digi_ene);

              sum_digi_ene+=digi_ene;

              crystal_d = digi->GetCrystal();
              row_d = digi->GetRow();
              copy_d = digi->GetCopy();

              Double_t digi_z=digi_pos.Z();
              hz_d->Fill(digi_z);

              if (copy_d==1){
                digi_ene_cp1+=digi->GetEnergy();
                hz_cp1_d->Fill(digi_z);
              }
              if (copy_d==2){
                digi_ene_cp2+=digi->GetEnergy();
                hz_cp2_d->Fill(digi_z);
              }
              if (copy_d==3){
                digi_ene_cp3+=digi->GetEnergy();
                hz_cp3_d->Fill(digi_z);
              }
              if (copy_d==4){
                digi_ene_cp4+=digi->GetEnergy();
                hz_cp4_d->Fill(digi_z);
              }

              digi_theta=(digi->GetTheta())*(180./TMath::Pi());
              digi_phi=(digi->GetPhi())*(180./TMath::Pi());

              h10digi->Fill(theta_mc);
              h20digi->Fill(phi_mc);

              h10ddigi->Fill(digi_theta);
              h20ddigi->Fill(digi_phi);

              h1d->Fill(digi_theta-theta_mc);
              h2d->Fill(digi_phi-phi_mc);

              hxpad_x_d->Fill(digi->GetXPad(),digi->GetThetaInt());
              hypad_y_d->Fill(digi->GetYPad(),digi->GetPhiInt());

              hx_row_d->Fill(digi->GetRow(),digi->GetThetaInt());

              id_d = digi->GetDetectorId();

              crystal_did = (id_d%10000);
              row_did = ((id_d/1000000)%100);

              hmodule_d->Fill(module_d);
              hrow_crystal_d->Fill(crystal_did,row_did);
              hx_y_d->Fill(digi->GetThetaInt(),digi->GetPhiInt());

              hid_d->Fill(id_d);
            }
          }
        hdigi_ene->Fill(sum_digi_ene);
        hdigi_ene_cp1->Fill(digi_ene_cp1);
        hdigi_ene_cp2->Fill(digi_ene_cp2);
        hdigi_ene_cp3->Fill(digi_ene_cp3);
        hdigi_ene_cp4->Fill(digi_ene_cp4);

        // reconstructed clusters
        Int_t ClustNoOrig,ClustNo,ClustNo_cp1,ClustNo_cp2,ClustNo_cp3,ClustNo_cp4;
        Int_t detectorID, module, mod3;
        if (cluster_array->GetEntriesFast()>0)
          {
            Int_t    idWithHighestEnergy = 0;
            Double_t highestEnergy = -1.;

            // First find the cluster with the highest energy

            if (module_h==3){

              ClustNoOrig=ClustNo=
                ClustNo_cp1=ClustNo_cp2=ClustNo_cp3=ClustNo_cp4=0;

              ClustNoOrig=cluster_array->GetEntriesFast();
              for (Int_t i=0; i<ClustNoOrig; i++)
                {
                  PndEmcCluster *cluster=(PndEmcCluster*)cluster_array->At(i);
                  cluster_energy=cluster->energy();

                  std::vector<PndEmcDigi*> digiList=cluster->DigiList();
                  ndigi=digiList.size();

                  if(ndigi != 0)
                    {
                      for(int k = 0; k < digiList.size(); k++)
                        {
                          detectorID = digiList[k]->GetDetectorId();
                          module = digiList[k]->GetModule();
                          row_c = digiList[k]->GetRow();
                          crystal_c = digiList[k]->GetCrystal();
                          copy_c = digiList[k]->GetCopy();

                          crystal_cid = (detectorID%10000);
                          row_cid = ((detectorID/1000000)%100);

                          hrow_crystal_c->Fill(crystal_cid,row_cid);
                          //cluster_energy=cluster->energy();
                        }
                    }
                  mod3 = module;

                  if (cluster_energy>highestEnergy)
                    {
                      idWithHighestEnergy = i;
                      highestEnergy = cluster_energy;
                      ClustNo++;
                    }
                } // number of cluster

              // Lets analyze that cluster!
              if (mod3 !=3 )cout << "mod3 = "<< mod3 << endl;
              PndEmcCluster *cluster=(PndEmcCluster*)cluster_array->At(idWithHighestEnergy);

              TVector3 cluster_pos=cluster->where();
              cluster_theta=cluster_pos.Theta()*(180./TMath::Pi());
              cluster_phi=cluster_pos.Phi()*(180./TMath::Pi());//+360;

              cluster_energy=cluster->energy();
              //ClustNo++;

              if (copy_c==1){
                hclus_ene_cp1->Fill(cluster->energy());
                hz_cp1_c->Fill(cluster->z());
                ClustNo_cp1++;
              }
              if (copy_c==2){
                hclus_ene_cp2->Fill(cluster->energy());
                hz_cp2_c->Fill(cluster->z());
                ClustNo_cp2++;
              }
              if (copy_c==3){
                hclus_ene_cp3->Fill(cluster->energy());
                hz_cp3_c->Fill(cluster->z());
                ClustNo_cp3++;
              }
              if (copy_c==4){
                hclus_ene_cp4->Fill(cluster->energy());
                hz_cp4_c->Fill(cluster->z());
                ClustNo_cp4++;
              }

              hz_clust->Fill(cluster->z());

              //cout<<"CLUSTER: th="<<cluster_theta<<", ph="<<cluster_phi<<", energy="<<cluster_energy<<endl;

              sum_clus_ene+=cluster_energy;

              cemc=p4mom.E();
              cemc_sum+=cemc;

              hclus_ene_mc->Fill(cemc);
              hclus_ene->Fill(cluster_energy);

              hthc_mc->Fill(theta_mc); // Monte Carlo
              hphc_mc->Fill(phi_mc);   // Monte Carlo

              hthc->Fill(cluster_theta); // reconstructed
              hphc->Fill(cluster_phi);   // reconstructed

              hthc_diff->Fill(cluster_theta-theta_mc);
              hphc_diff->Fill(cluster_phi-phi_mc);
            }
          }
        hclus_ene_sum_mc->Fill(cemc_sum);

        hclus_ene_sum->Fill(sum_clus_ene);

        hClustNoOrig->Fill(ClustNoOrig);
        hClustNo->Fill(ClustNo);

        hClNo_cp1->Fill(ClustNo_cp1);
        hClNo_cp2->Fill(ClustNo_cp2);
        hClNo_cp3->Fill(ClustNo_cp3);
        hClNo_cp4->Fill(ClustNo_cp4);
   }

hist3 GetXaxis

(

)

-> SetTitle("Z position")

Definition at line 653 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by Pic_FEP_DEGASComparison(), Pic_FEP_Eff_lEE_XiManyTargets(), Pic_OmegaQuadrupoleComparison(), PndFtsHoughSpace::setParametersForHsOption(), PndFtsHoughSpace::WriteHistoOfAllPaths(), and PndFtsHoughSpace::WriteHistoOfAllPathsForEachMcTruthTrack().

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

hresphi_phi GetYaxis

(

)

-> SetTitle("Counts")

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

Referenced by Pic_FEP_DEGASComparison(), Pic_FEP_Eff_lEE_XiManyTargets(), Pic_OmegaQuadrupoleComparison(), PndFtsHoughSpace::setParametersForHsOption(), PndFtsHoughSpace::WriteHistoOfAllPaths(), and PndFtsHoughSpace::WriteHistoOfAllPathsForEachMcTruthTrack().

{reco} - #phi_{truth}, degree");

gROOT Macro

(

"$VMCWORKDIR/gconfig/rootlogon.C"

)

c9 Modified

(

)

copies2a SetBorderSize

(

0

)

tsim SetBranchAddress	(	"MCTrack"	,
		&	mctrack_array
	)

tful SetBranchAddress	(	"EmcHit"	,
		&	hit_array
	)

tful SetBranchAddress	(	"EmcDigi"	,
		&	digi_array
	)

tful SetBranchAddress	(	"EmcCluster"	,
		&	cluster_array
	)

copies2a SetFillColor

(

0

)

ModuleVolume SetLineColor

(

2

)

hclus_ene_cp2 SetLineColor

(

3

)

h20 SetLineColor

(

4

)

hMcGemPhi SetLineColor

(

6

)

gStyle SetPalette

(

1

)

copies2a SetTextSize

(

0.

05

)

Variable Documentation

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

Definition at line 444 of file anal_hit_digi_cluster_fwendcap.C.

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

Definition at line 648 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by if(), and main().

Double_t cemc

Definition at line 31 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by for().

TClonesArray* cluster_array =new TClonesArray("PndEmcCluster")

Definition at line 27 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), emc_correction_hist(), emc_correction_parametrization(), emc_correction_QA(), error_matrix_fit(), fit_resolution(), QAmacro_emc_3(), reco_analys(), and vis_mc_ids().

Double_t cluster_energy

Definition at line 31 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), emc_correction_hist(), emc_correction_parametrization(), emc_correction_QA(), error_matrix_fit(), fit_resolution(), for(), QAmacro_emc_3(), and reco_analys().

Double_t cluster_phi

Definition at line 31 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by emc_correction_hist(), emc_correction_parametrization(), emc_correction_QA(), error_matrix_fit(), QAmacro_emc_3(), and reco_analys().

Double_t cluster_theta

Definition at line 31 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by emc_correction_hist(), emc_correction_parametrization(), emc_correction_QA(), error_matrix_fit(), QAmacro_emc_3(), and reco_analys().

TLegend* copies = new TLegend(0.65,0.4,0.85,0.6)

Definition at line 451 of file anal_hit_digi_cluster_fwendcap.C.

TLegend* copies2a = new TLegend(0.20,0.4,0.5,0.6)

Definition at line 616 of file anal_hit_digi_cluster_fwendcap.C.

TLegend* copies3 = new TLegend(0.20,0.30,0.5,0.5)

Definition at line 583 of file anal_hit_digi_cluster_fwendcap.C.

TLegend* copies4 = new TLegend(0.20,0.40,0.5,0.6)

Definition at line 485 of file anal_hit_digi_cluster_fwendcap.C.

TLegend* copies5 = new TLegend(0.20,0.40,0.5,0.6)

Definition at line 551 of file anal_hit_digi_cluster_fwendcap.C.

TLegend* copies6 = new TLegend(0.65,0.4,0.85,0.6)

Definition at line 518 of file anal_hit_digi_cluster_fwendcap.C.

Double_t copy_c

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t copy_d

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Double_t copy_h

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t crystal_c

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t crystal_cid

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t crystal_d

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Double_t crystal_did

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

Double_t crystal_h

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t crystal_hid

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t d_id

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* dene_mc = new TH1F("dene_mc","MC energy DIGI (GeV)",100,0.05,1.05)

Definition at line 93 of file anal_hit_digi_cluster_fwendcap.C.

TClonesArray* digi_array =new TClonesArray("PndEmcDigi")

Definition at line 24 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), QAmacro_emc_3(), and reco_analys().

Double_t digi_ene

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), and for().

Double_t digi_phi

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

Double_t digi_theta

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

PndEmcMapper* emcMap =PndEmcMapper::Instance(1)

Definition at line 13 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by MyMainFrame::DoRecreate(), PndEmcWaveform::GetTCI(), PndEmcDigi::GetTCI(), PndEmcDigi::PndEmcDigi(), PndEmcDigi::SetDetectorId(), and structure_analysis_emc().

Double_t ene_mc

Definition at line 31 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

TFile* fful = new TFile("data/fixed/full_emc_g4_m3_1000eve_fixed.root")

Definition at line 17 of file anal_hit_digi_cluster_fwendcap.C.

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

Definition at line 8 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), MyMainFrame::DoRecreate(), QAmacro_emc_3(), QAmacro_emc_4(), PndParticleQATask::SetFastSim(), and structure_analysis_emc().

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

Definition at line 121 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), CalculateCrystalMatricesZ(), checkphsp2_2(), PndMvdRiemannVertexFinderTask::CheckTwoCands(), create_trap(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), cutfinderx(), diffQA(), DrawField(), PndFTSCADisplay::DrawGBTrackFast(), drawHistos(), drawMvdHistos(), emc_correction_hist(), EventListing(), FCDrawVariable(), FCQaVar(), PndFTSCATrackParamVector::FilterWithMaterial(), fit_resolution(), PndMvdRiemannVertexFinderTask::FoundCandInMCCands(), KFParticle::GetDistanceFromVertexXY(), KFParticleSIMD::GetDistanceFromVertexXY(), PndRiemannTrackFinder::HitTooClose(), KFParticleBase::Initialize(), KFParticleBaseSIMD::Initialize(), MergeRootfile(), NHitsPerEvent(), PndEmcStructure::PndEmcStructure(), PndLine::PndLine(), QAmacro_emc_3(), QAmacro_emc_4(), simanalys(), testFilter(), TestIt(), and TrackHits().

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

Definition at line 117 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

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

Definition at line 119 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

TH1F* h10ddigi = new TH1F("h10ddigi","DIGI Theta (data)",180,0.,180.)

Definition at line 126 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* h10digi = new TH1F("h10digi","DIGI Theta (MC)",180,0.,180.)

Definition at line 124 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* h1d = new TH1F("h1d","DIGI Theta difference",100,-5.,5.)

Definition at line 128 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by PndCAGBTracker::CATrackFinder(), and PndCAGBTracker::PickUpHits().

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

Definition at line 122 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), PndAnalysisForwardTask::analyzeMCTracks(), PndTrkTrackFinder::Apollonius(), CalculateCrystalMatricesZ(), check_complete(), check_pid(), checkphsp2_2(), PndMvdRiemannVertexFinderTask::CheckTwoCands(), comp_complete(), comp_pid(), comp_recoqa(), PndTrackCollection::copyHit(), PndFtsLineApproximator::copyHitWithNewPosition(), PndLineApproximation::copyHitWithNewPosition(), create_trap(), cutfinder(), cutfinder_fullev(), cutfinder_toyev(), cutfinderx(), diffQA(), PndFTSCADisplay::DrawGBTrackFast(), drawHistos(), PndFTSCAGBTracker::EstimatePV(), FCDrawVariable(), FCQaVar(), PndMvdRiemannVertexFinderTask::FoundCandInMCCands(), Hough(), HypHitAnasigBg_AR(), if(), KFParticleBase::Initialize(), KFParticleBaseSIMD::Initialize(), main(), MergeRootfile(), overall_perf(), overall_perf_full(), overall_perf_toy(), PndEmcStructure::PndEmcStructure(), QAmacro_emc_3(), QAmacro_emc_4(), PndMQTopix4OnlineHisto::Run(), testFilter(), TestIt(), TrackHits(), PndCATrackParamVector::Transport0ToX(), PndCATrackParamVector::TransportJ0ToX0(), PndFTSCATrackParamVector::TransportJ0ToX0(), PndCATrackParam::TransportToX(), PndCATrackParamVector::TransportToX(), PndFTSCATrackParam::TransportToX0(), and PndFTSCATrackParamVector::TransportToX0().

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

Definition at line 118 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

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

Definition at line 120 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

TH1F* h20ddigi = new TH1F("h20ddigi","DIGI Phi (data)",360,-180.,180.)

Definition at line 127 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* h20digi = new TH1F("h20digi","DIGI Phi (MC)",360,-180.,180.)

Definition at line 125 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* h2d = new TH1F("h2d","DIGI Phi difference",100,-5.,5.)

Definition at line 129 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hClNo_cp1 = new TH1F("hClNo_cp1","hClNo_cp1",10,0,9)

Definition at line 149 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hClNo_cp2 = new TH1F("hClNo_cp2","hClNo_cp2",10,0,9)

Definition at line 150 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hClNo_cp3 = new TH1F("hClNo_cp3","hClNo_cp3",10,0,9)

Definition at line 151 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hClNo_cp4 = new TH1F("hClNo_cp4","hClNo_cp4",10,0,9)

Definition at line 152 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hclus_ene = new TH1F("hclus_ene","SUM of Cluster energy (GeV)",100,0.05,1.05)

Definition at line 101 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hclus_ene_cp1 = new TH1F("hclus_ene_cp1","Cluster energy (GeV) _cp1",100,0.05,1.05)

Definition at line 111 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hclus_ene_cp2 = new TH1F("hclus_ene_cp2","Cluster energy (GeV) _cp2",100,0.05,1.05)

Definition at line 112 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hclus_ene_cp3 = new TH1F("hclus_ene_cp3","Cluster energy (GeV) _cp3",100,0.05,1.05)

Definition at line 113 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hclus_ene_cp4 = new TH1F("hclus_ene_cp4","Cluster energy (GeV) _cp4",100,0.05,1.05)

Definition at line 114 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hclus_ene_mc = new TH1F("hclus_ene_mc","SUM of MC Cluster energy (GeV)",100,0.05,1.05)

Definition at line 103 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hclus_ene_sum = new TH1F("hclus_ene_sum","SUM of Cluster energy (GeV)",100,0.05,1.05)

Definition at line 104 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hclus_ene_sum_mc = new TH1F("hclus_ene_sum_mc","SUM of MC Cluster energy (GeV)",100,0.05,1.05)

Definition at line 94 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hClustNo = new TH1F("hClustNo","hClustNo",10,0,9)

Definition at line 147 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hClustNoOrig = new TH1F("hClustNoOrig","hClustNoOrig",10,0,9)

Definition at line 146 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hcp_z_h = new TH2F("cp_z_h","Copy vs.Z pos (HITS)",4,1,5,80,200,240)

Definition at line 72 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hdigi_ene = new TH1F("hdigi_ene","SUM of Digi energy (GeV)",100,0.05,1.05)

Definition at line 100 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hdigi_ene_cp1 = new TH1F("hdigi_ene_cp1","SUM of Digi energy (GeV) _cp1",100,0.05,1.05)

Definition at line 131 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hdigi_ene_cp2 = new TH1F("hdigi_ene_cp2","SUM of Digi energy (GeV) _cp2",100,0.05,1.05)

Definition at line 132 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hdigi_ene_cp3 = new TH1F("hdigi_ene_cp3","SUM of Digi energy (GeV) _cp3",100,0.05,1.05)

Definition at line 133 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hdigi_ene_cp4 = new TH1F("hdigi_ene_cp4","SUM of Digi energy (GeV) _cp4",100,0.05,1.05)

Definition at line 134 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hene_d = new TH1F("hene_d","DATA energy DIGI (GeV)",100,0.05,1.05)

Definition at line 97 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hene_h = new TH1F("hene_h","DATA energy HIT (GeV)",100,0.05,1.05)

Definition at line 96 of file anal_hit_digi_cluster_fwendcap.C.

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

Definition at line 92 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

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

Definition at line 99 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hhit_ene_cp1 = new TH1F("hhit_ene_cp1","SUM of Hits energy (GeV) _cp1",100,0.05,1.05)

Definition at line 106 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hhit_ene_cp2 = new TH1F("hhit_ene_cp2","SUM of Hits energy (GeV) _cp2",100,0.05,1.05)

Definition at line 107 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hhit_ene_cp3 = new TH1F("hhit_ene_cp3","SUM of Hits energy (GeV) _cp3",100,0.05,1.05)

Definition at line 108 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hhit_ene_cp4 = new TH1F("hhit_ene_cp4","SUM of Hits energy (GeV) _cp4",100,0.05,1.05)

Definition at line 109 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hid_d = new TH1F("id_d","id (DIGI)",140000000,300000000,340000000)

Definition at line 143 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hid_h = new TH1F("id_h","id (HIT)",140000000,300000000,340000000)

Definition at line 59 of file anal_hit_digi_cluster_fwendcap.C.

TClonesArray* hit_array =new TClonesArray("PndEmcHit")

Definition at line 21 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by ana_MCOpt(), ana_MCpid(), anaGemPointrate(), anaGemSmearing(), QAmacro_emc_4(), PndDrcAccuDigiPixelDraw::ReadAllHits(), reco_analys(), and vis_mc_ids().

Double_t hit_ene

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by for().

Double_t hit_id

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t hit_phi

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t hit_theta

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t hit_x

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by PndDiscTaskDigitization::Exec().

Double_t hit_y

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hmodule_d = new TH1F("module_d","Module (DIGI)",5,1,6)

Definition at line 144 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hmodule_h = new TH1F("module_h","Module (HIT)",5,1,6)

Definition at line 61 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hphc = new TH1F("hphc","CLUSTER Phi (data)",360,-180.,180.)

Definition at line 139 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hphc_diff = new TH1F("hphc_diff","CLUSTER Phi difference",100,-5.,5.)

Definition at line 141 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hphc_mc = new TH1F("hphcmc","CLUSTER Phi (MC)",360,-180.,180.)

Definition at line 137 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hrow_crystal_c = new TH2F("row_crystal_c","row vs crystal (CLUSTER)",72,-36,36,72,-36,36)

Definition at line 54 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hrow_crystal_d = new TH2F("row_crystal_digi","row vs crystal (DIGI)",72,-36,36,72,-36,36)

Definition at line 53 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hrow_crystal_h = new TH2F("row_crystal","row vs crystal (HIT)",72,-36,36,72,-36,36)

Definition at line 52 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hthc = new TH1F("hthc","CLUSTER Theta (data)",180,0.,180.)

Definition at line 138 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hthc_diff = new TH1F("hthc_diff","CLUSTER Theta difference",100,-5.,5.)

Definition at line 140 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hthc_mc = new TH1F("hthcmc","CLUSTER Theta (MC)",180,0.,180.)

Definition at line 136 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hx_row_d = new TH2F("x_row_d","GetThetaInt() vs row (DIGI)",200,200,300,72,-36,36)

Definition at line 70 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hx_row_h = new TH2F("x_row_h","x vs row (HIT)",200,-100,100,72,-36,36)

Definition at line 69 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hx_y_d = new TH2F("x_y_d","GetThetaInt() vs GetPhiInt() (DIGI)",200,200,300,200,200,300)

Definition at line 57 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hx_y_h = new TH2F("x_y_h","x vs y (HIT)",200,-100,100,200,-100,100)

Definition at line 56 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hxpad_x_d = new TH2F("xpad_x_d","xpad vs GetThetaInt() (DIGI)",200,-100,100,200,200,300)

Definition at line 66 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hxpad_x_h = new TH2F("xpad_x_h","xpad vs x (HIT)",200,-100,100,200,-100,100)

Definition at line 63 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hypad_y_d = new TH2F("ypad_y_d","ypad vs GetPhiInt() (DIGI)",200,-100,100,200,200,300)

Definition at line 67 of file anal_hit_digi_cluster_fwendcap.C.

TH2F* hypad_y_h = new TH2F("ypad_y_h","ypad vs y (HIT)",200,-100,100,200,-100,100)

Definition at line 64 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_clust = new TH1F("z_cluster","z (CLUSTER)",300,1,300)

Definition at line 85 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp1_c = new TH1F("z_cp1_c","z (CLUSTER)",300,1,300)

Definition at line 86 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp1_d = new TH1F("z_cp1_d","z (DIGI)",300,1,300)

Definition at line 80 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp1_h = new TH1F("z_cp1_h","z (HIT)",300,1,300)

Definition at line 74 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp2_c = new TH1F("z_cp2_c","z (CLUSTER)",300,1,300)

Definition at line 87 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp2_d = new TH1F("z_cp2_d","z (DIGI)",300,1,300)

Definition at line 81 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp2_h = new TH1F("z_cp2_h","z (HIT)",300,1,300)

Definition at line 75 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp3_c = new TH1F("z_cp3_c","z (CLUSTER)",300,1,300)

Definition at line 88 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp3_d = new TH1F("z_cp3_d","z (DIGI)",300,1,300)

Definition at line 82 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp3_h = new TH1F("z_cp3_h","z (HIT)",300,1,300)

Definition at line 76 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp4_c = new TH1F("z_cp4_c","z (CLUSTER)",300,1,300)

Definition at line 89 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp4_d = new TH1F("z_cp4_d","z (DIGI)",300,1,300)

Definition at line 83 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_cp4_h = new TH1F("z_cp4_h","z (HIT)",300,1,300)

Definition at line 77 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_d = new TH1F("z_d","z (DIGI)",300,1,300)

Definition at line 79 of file anal_hit_digi_cluster_fwendcap.C.

TH1F* hz_h = new TH1F("z_h","z (HIT)",300,1,300)

Definition at line 60 of file anal_hit_digi_cluster_fwendcap.C.

Int_t id_d

Definition at line 43 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

Int_t id_h

Definition at line 43 of file anal_hit_digi_cluster_fwendcap.C.

double max_energy1 =0

Definition at line 49 of file anal_hit_digi_cluster_fwendcap.C.

TClonesArray* mctrack_array =new TClonesArray("PndMCTrack")

Definition at line 11 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), emc_correction_QA(), error_matrix_fit(), QAmacro_emc_3(), QAmacro_evtgen_3(), QAmacro_evtgen_4(), and vis_mc_ids().

Double_t module_d

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

Double_t module_h

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Int_t ncounts = tsim->GetEntries()

Definition at line 155 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_cluster_energyCorrection_Ntuple_suplement(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), QAmacro_evtgen_3(), QAmacro_evtgen_4(), and vis_mc_ids().

int ndigi

Definition at line 48 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), QAmacro_emc_3(), and reco_analys().

int nhit

Definition at line 48 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by PndTrkLegendreSecTask::Cleanup(), PndTrkLegendreTask::Cleanup(), PndTrkLegendreSecTask2::Cleanup(), PndSciTAnaIdeal::Exec(), PndSciTAnaIdeal::SetSciTEventCorr(), and vis_mc_ids().

TLorentzVector p4mom

Definition at line 46 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), emc_correction_hist(), emc_correction_parametrization(), and reco_analys().

Double_t phi_mc

Definition at line 31 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

TVector3 photon_momentum

Definition at line 45 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), error_matrix_fit(), PndDiscTaskDigitization::Exec(), and QAmacro_emc_3().

Double_t row_c

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t row_cid

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t row_d

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Double_t row_did

Definition at line 39 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

Double_t row_h

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t row_hid

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

TTree* tful =(TTree *) fful->Get("pndsim")

Definition at line 19 of file anal_hit_digi_cluster_fwendcap.C.

Double_t theta_mc

Definition at line 31 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), and analysis_digi_cluster_fwendcap().

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

Definition at line 10 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by analysis_cluster_energyCorrection_Ntuple(), analysis_digi_cluster_7gammaAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_FullEmc(), analysis_digi_cluster_allGammasAnalysis_fwendcap(), analysis_digi_cluster_fwendcap(), MyMainFrame::DoRecreate(), QAmacro_emc_3(), QAmacro_emc_4(), rad_dose_studies(), structure_analysis_emc(), and time_based_studies().

Double_t x_h

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Referenced by PndSdsChargeWeightingAlgorithms::HeadTail().

Double_t y_h

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.

Double_t z_h

Definition at line 33 of file anal_hit_digi_cluster_fwendcap.C.