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

tsim	SetBranchAddress ("MCTrack",&mctrack_array)

tsim	SetBranchAddress ("EmcPoint",&point_array)

tsim	SetBranchAddress ("EmcHit",&hit_array)

	if (nhits!=0)

HISThit_ene	Fill (sum_hit_ene)

Variables
TFile *	fsim = new TFile("../../emc_complete_1GeVphoton_FwEndCap.root")

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

TClonesArray *	mctrack_array =new TClonesArray("PndMCTrack")

TClonesArray *	point_array =new TClonesArray("PndEmcPoint")

TClonesArray *	hit_array =new TClonesArray("PndEmcHit")

Double_t	theta_mc =photon_momentum.Theta()*(180./TMath::Pi())

Double_t	phi_mc =photon_momentum.Phi()*(180./TMath::Pi())

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	hit_ene

Double_t	hit_theta

Double_t	hit_phi

Double_t	hit_id

Double_t	hit_x

Double_t	hit_y

Int_t	id_p

Int_t	crystal_pid

Int_t	row_pid

TVector3	photon_momentum =mctrack->GetMomentum()

TLorentzVector	p4mom

TH2F *	hypad_xpad_p = new TH2F("ypad_xpad","crystal row vs col number (derived from POINT)",74,-37,37,74,-37,37)

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

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

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

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

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

TH1F *	hmodule_p = new TH1F("module_p","Module (derived from POINT)",5,0.5,5.5)

TH2F *	hxpad_x_p = new TH2F("xpad_x_p","x vs xpad (derived from POINT)",2*37.5,-37.5,37.5,211,-105,105)

TH2F *	hypad_y_p = new TH2F("ypad_y_p","y vs ypad (derived from POINT)",2*37.5,-37.5,37.5,211,-105,105)

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 Point energies (GeV)",100,0.0,1.05)

TH1F *	hpoi_ene11_12 = new TH1F("hpoi_ene11_12","crystal [11,12] dep. energy (GeV)",100,0.0001,0.2)

TH1F *	hpoi_ene12_12 = new TH1F("hpoi_ene12_12","crystal [12,12] dep. energy (GeV)",100,0.0001,0.2)

TH1F *	hpoi_ene13_12 = new TH1F("hpoi_ene13_12","crystal [13,12] dep. energy (GeV)",100,0.0001,0.2)

TH1F *	hpoi_ene14_12 = new TH1F("hpoi_ene14_12","crystal [14,12] dep. energy (GeV)",100,0.0001,0.2)

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

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

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

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

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

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

Int_t	ncounts = tsim->GetEntriesFast()

cout<< "Number of events = " << ncounts<< endl;Int_t module_h, id_h, crystal_hid, row_hid;TH1F HISTene_h=new TH1F("HISTene_h","energy (hit)", 1000, 0.0001, 1.01);TH1F HISThit_ene=new TH1F("HISThit_ene","FwEndCap energy (hit)", 1000, 0.0, 1.05);TH1F HISTz_h=new TH1F("HISTz_h","z (hit)", 100, 200, 300);TH2F HISTrow_x_h=new TH2F("HISTrow_x_h","row vs x (hit)", 200,-100, 100, 2 38.5,-38.5, 38.5);TH2F HISTrow_crystal_h=new TH2F("HISTrow_crystal_h","row vs column (hit)", 2 38.5,-38.5, 38.5, 2 38.5,-38.5, 38.5);TH2F HISTy_x_h=new TH2F("HISTy_x_h","y vs x (hit)", 211,-105, 105, 211,-105, 105);TH1F HISTtheta_h=new TH1F("HISTtheta_h","theta (hit)", 180, 0, 180);TH2I HISTxpad_x_h=new TH2I("HISTxpad_x_h","", 10500,-105, 105, 100 38.5,-38.5, 38.5);TH2I HISTypad_y_h=new TH2I("HISTypad_y_h","", 10500,-105, 105, 100 38.5,-38.5, 38.5);TH1F HISTOhit_ene11_12=new TH1F("HISTOhit_ene11_12","crystal [11,12] dep. energy (GeV)", 100, 0.0001, 1.01);TH1F HISTOhit_ene12_12=new TH1F("HISTOhit_ene12_12","crystal [12,12] dep. energy (GeV)", 100, 0.0001, 1.01);TH1F HISTOhit_ene13_12=new TH1F("HISTOhit_ene13_12","crystal [13,12] dep. energy (GeV)", 100, 0.0001, 1.01);TH1F HISTOhit_ene14_12=new TH1F("HISTOhit_ene14_12","crystal [14,12] dep. energy (GeV)", 100, 0.0001, 1.01);Int_t totalhits=0;Int_t efficiencyCounter=0;Int_t atLeastOneHitPerEventCounter=0;for(int j=0;j< ncounts;j++){Int_t efficiencyFlag=0;Int_t atLeastOneHitPerEventFlag=0;Double_t sum_hit_ene=0;tsim->	GetEntry (j)

PndMCTrack *	mctrack =(PndMCTrack *) mctrack_array->At(0)

Int_t	nhits = hit_array->GetEntries()

	totalhits =nhits

Function Documentation

HISThit_ene Fill ( sum_hit_ene )

Referenced by PndGemFindHitsAna::AnaHistos(), PndHypGeCOSYBackgroundAna::Exec(), PndGemFindHitsAna::Exec(), PndFtsHoughSpace::FillHoughSpace(), PndFtsHoughSpace::FindAllPeaksWithTSpectrum2(), gem_material_ana1(), and materialana().

if ( nhits! = 0 )

Definition at line 258 of file analysis_point_hit_fwendcap.C.

References PndEmcHit::GetEnergy(), PndEmcHit::GetModule(), PndEmcHit::GetXPad(), PndEmcHit::GetYPad(), hit(), hit_ene, i, and Pi.

           {
             PndEmcHit *hit=(PndEmcHit*)hit_array->At(i); // get the i-th hit (detector) of the event
             module_h = hit->GetModule();
 
             if (module_h==3){
               //hene_mc->Fill(p4mom.E());
               hit_ene=hit->GetEnergy();
               HISTene_h->Fill(hit_ene);
               if (hit_ene >= 0.010) // 10 MeV threshold used (here only) to calculate the efficiency of the FwEndCap geometry
                 efficiencyFlag = 1;
 
 
               //filling some histograms for a few crystals
               if (hit->GetXPad()==11 && hit->GetYPad()==12)
                 HISTOhit_ene11_12->Fill(hit_ene);
               else if (hit->GetXPad()==12 && hit->GetYPad()==12)
                 HISTOhit_ene12_12->Fill(hit_ene);
               else if (hit->GetXPad()==13 && hit->GetYPad()==12)
                 HISTOhit_ene13_12->Fill(hit_ene);
               else if (hit->GetXPad()==14 && hit->GetYPad()==12)
                 HISTOhit_ene14_12->Fill(hit_ene);
 
               sum_hit_ene+=hit_ene;
 
               TVector3 hit_pos(hit->GetX(),hit->GetY(),hit->GetZ()); // position of the center of the hit crystal
               hit_theta=hit_pos.Theta()*(180./TMath::Pi());
               hit_phi=hit_pos.Phi()*(180./TMath::Pi());
 
               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);
 
               HISTz_h->Fill(hit->GetZ());
               if (hit->GetXPad() >= 1)
                 HISTxpad_x_h->Fill(hit->GetX(),hit->GetXPad()-1);
               else
                 HISTxpad_x_h->Fill(hit->GetX(),hit->GetXPad());
 
               if (hit->GetYPad() >= 1)
                 HISTypad_y_h->Fill(hit->GetY(),hit->GetYPad()-1);
               else
                 HISTypad_y_h->Fill(hit->GetY(),hit->GetYPad());
 
               if (hit->GetYPad() >= 1)
                 HISTrow_x_h->Fill(hit->GetX(),hit->GetYPad()-1);
               else
                 HISTrow_x_h->Fill(hit->GetX(),hit->GetYPad());
 
               id_h = hit->GetDetectorID();
               crystal_hid = -(id_h%10000 -36);     //crystal column number (the same as GetXPad())
               row_hid = (id_h/1000000)%100 - 37;   //crystal row number (the same as GetYPad())
 
               if (hit->GetXPad() <= -1 && hit->GetYPad() <= -1)
                 HISTrow_crystal_h->Fill(hit->GetXPad(),hit->GetYPad());
               else if (hit->GetXPad() <= -1 && hit->GetYPad() >= 1)
                 HISTrow_crystal_h->Fill(hit->GetXPad(),hit->GetYPad()-1);
               else if (hit->GetXPad() >= 1 && hit->GetYPad() <= -1)
                 HISTrow_crystal_h->Fill(hit->GetXPad()-1,hit->GetYPad());
               else //if (hit->GetXPad() >= 1 && hit->GetYPad() >= 1)
                 HISTrow_crystal_h->Fill(hit->GetXPad()-1,hit->GetYPad()-1);
 
               HISTy_x_h->Fill(hit->GetX(),hit->GetY());
             }
           }