DrawHits Class Reference

#include <DrawHits.h>

Inheritance diagram for DrawHits:

Public Member Functions
	DrawHits ()
	DrawHits (Int_t verbose)
virtual	~DrawHits ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *option)
virtual void	Finish ()
void	RotBarCoordinate (TVector3 &vector, Int_t barID)
void	SetTreeName (TString str)

Private Member Functions
void	ProcessPhotonHit ()
void	ProcessPhotonMC ()
void	WriteToFile ()
void	CreateHisto ()
void	DrawHisto ()
void	ProcessBarHit ()
void	SetDefaultParameters ()
void	SetPlotStyle ()
Bool_t	InsideBar (Double_t, Double_t, Double_t, Double_t, Double_t, Double_t, Double_t, Double_t, Double_t, Double_t)

Private Attributes
Int_t	fDetectorID
Double_t	fMass
Double_t	fThetaC
Double_t	fErrThetaC
Int_t	fRefIndex
TClonesArray *	fBarPointArray
TClonesArray *	fPDPointArray
TClonesArray *	fHitArray
TClonesArray *	fPDHitArray
TClonesArray *	fMCArray
TList *	fHistoList
PndGeoDrcPar *	fPar
PndGeoDrc *	fGeo
	Basic geometry data of barrel DRC. More...
Int_t	fVerbose
Int_t	nevents
Double_t	xmcp
Double_t	ymcp
Double_t	zmcp
Double_t	pxmc
Double_t	pymc
Double_t	pzmc
Float_t	xhit
Float_t	yhit
Double_t	zhit
Double_t	thit
Double_t	xbar
Double_t	ybar
Double_t	zbar
Double_t	pxMo
Double_t	pyMo
Double_t	pzMo
Double_t	pxPho
Double_t	pyPho
Double_t	pzPho
Double_t	xEnt
Double_t	yEnt
Double_t	zEnt
Double_t	lam
Double_t	fPhiRot
TVector3	fStartVertex
Double_t	fPx
Double_t	fPy
Double_t	fPz
Double_t	fXcross
Double_t	fYcross
TVector3	fPphoB
Double_t	fkxBar
Double_t	fkyBar
Double_t	fkzBar
TString	fTreeName
TTree *	phoTree
std::vector< TH1F * >	histos
TH1F *	currh
TGraph2D *	map
TGraph2D *	pop
TVector3	fPMo
TH1D *	fhThetaC
TH2D *	fhThetaCMass
TH2D *	fhMomAng
TH2D *	fhThetaCMomK
TH2D *	fhThetaCMomP
TH2D *	fhThetaCMomM
TH2D *	fhThetaCMomE
TH1F *	fhPhoTheta
TH2F *	fhPDPlane
TH1D *	fhCHrealMC
TH1D *	fhCHreal
TH2D *	fhCHlam
TH2D *	fhCHlamMC
TH2D *	fhCHlamE
TH2F *	fhHitsD
TH1D *	fhLambda
TH2D *	fhXYPDHit
TH1D *	fhLambdaMC
TH2D *	fhXYPDMCPt
TH1D *	fhPDTime
TH2D *	fhXYPDHitKp
TH2D *	fhXYPDHitKn
TH2D *	fhXYPDHitpip
TH2D *	fhXYPDHitpin
TH2D *	fhXYPDHitmp
TH2D *	fhXYPDHitmn
TH2D *	fhXYPDMCPtKp
TH2D *	fhXYPDMCPtKn
TH2D *	fhXYPDMCPtpip
TH2D *	fhXYPDMCPtpin
TH2D *	fhXYPDMCPtmp
TH2D *	fhXYPDMCPtmn

Detailed Description

DrawHits.h

Class for Analysising DRC Hits

Definition at line 47 of file DrawHits.h.

Constructor & Destructor Documentation

DrawHits::DrawHits

(

)

Default constructor

Definition at line 48 of file DrawHits.cxx.

References fGeo.

:FairTask("DrawHits")
{
  fGeo = new PndGeoDrc();
  
}

DrawHits::DrawHits

(

Int_t

verbose

)

Constructor with verbosity

Definition at line 56 of file DrawHits.cxx.

References fGeo, fVerbose, and verbose.

  :FairTask("DrawHits")
{
  fVerbose = verbose;  
  fGeo = new PndGeoDrc();
}

DrawHits::~DrawHits ( )

virtual

Destructor

Definition at line 63 of file DrawHits.cxx.

References fGeo, and fHistoList.

{
  if (fGeo) delete fGeo;
  fHistoList->Delete();  
  delete fHistoList;    
}

Member Function Documentation

void DrawHits::CreateHisto ( )

private

Definition at line 358 of file DrawHits.cxx.

References fhCHlam, fhCHlamE, fhCHlamMC, fhCHreal, fhCHrealMC, fHistoList, fhLambda, fhLambdaMC, fhMomAng, fhPDPlane, fhPDTime, fhPhoTheta, fhThetaC, fhThetaCMass, fhThetaCMomE, fhThetaCMomK, fhThetaCMomM, fhThetaCMomP, fhXYPDHit, fhXYPDHitKn, fhXYPDHitKp, fhXYPDHitmn, fhXYPDHitmp, fhXYPDHitpin, fhXYPDHitpip, fhXYPDMCPt, fhXYPDMCPtKn, fhXYPDMCPtKp, fhXYPDMCPtmn, fhXYPDMCPtmp, fhXYPDMCPtpin, fhXYPDMCPtpip, fTreeName, phoTree, pxMo, pxPho, pyMo, pyPho, pzMo, pzPho, thit, xbar, xEnt, xhit, ybar, yEnt, yhit, zbar, zEnt, and zhit.

Referenced by Init().

{
  // Creation of the TTree to hold the data for Cherenkov angle 
  phoTree = new TTree(fTreeName, "tree to get cherenkov angle");
  phoTree->Branch("HitX", &xhit, "HitX/D");
  phoTree->Branch("HitY", &yhit, "HitY/D");
  phoTree->Branch("HitZ", &zhit, "HitZ/D");
  phoTree->Branch("HitT", &thit, "HitT/D");
  phoTree->Branch("MotherPx", &pxMo, "MotherPx/D");
  phoTree->Branch("MotherPy", &pyMo, "MotherPy/D");
  phoTree->Branch("MotherPz", &pzMo, "MotherPz/D"); 
  phoTree->Branch("InitPhoPx", &pxPho, "InitPhoPx/D");
  phoTree->Branch("InitPhoPy", &pyPho, "InitPhoPy/D");
  phoTree->Branch("InitPhoPz", &pzPho, "InitPhoPz/D");
  phoTree->Branch("xEntry", &xEnt, "xEntry/D");
  phoTree->Branch("yEntry", &yEnt, "yEntry/D");
  phoTree->Branch("zEntry", &zEnt, "zEntry/D");
  phoTree->Branch("xBar", &xbar, "xBar/D");
  phoTree->Branch("yBar", &ybar, "yBar/D");
  phoTree->Branch("zBar", &zbar, "zBar/D");
  
  // Histogram list
  fHistoList = new TList(); 
  
  //Histogram for Bar Hits

  fhThetaC = new TH1D("fhThetaC", "Cherenkov Angle", 100, 0, 1);
  fhThetaCMass = new TH2D("fhThetaCMass", " ThetaC vs. Mass", 100, 0, 1, 100, 0, 0.6);
  fhThetaCMomK = new TH2D("fhThetaCMomK", " ThetaC vs. Mom", 100,  0, 8.0, 100, 0, 1.);
  fhThetaCMomP = new TH2D("fhThetaCMomP", " ThetaC vs. Mom", 100,  0, 8.0, 100, 0, 1.);
  fhThetaCMomM = new TH2D("fhThetaCMomM", " ThetaC vs. Mom", 100,  0, 8.0, 100, 0, 1.);
  fhThetaCMomE = new TH2D("fhThetaCMomE", " ThetaC vs. Mom", 100,  0, 8.0, 100, 0, 1.);
  fhMomAng = new TH2D("fhMomAng", " Mom vs. Angle", 100, 0, 140.0, 100, 0.0,8.0);
  fhPhoTheta = new TH1F("fhPhoTheta","Number of photons as a function of initial track theta", 150, 0., 150.); 
  fhPDPlane = new TH2F("fhPDPlane", "pixels hit by photons", 338,-109.85,109.85, 338,-109.85,109.85);
 
  //Histogram for Hits in Photon Detector
 
  fhLambda =new TH1D("fhLambda", "No of Photons  vs. Lambda", 200, 100, 1000);
  fhPDTime =new TH1D("fhPDTime", "Time in ns", 100, 0, 500);
  fhXYPDHit = new TH2D("fhXYPDHit", "XY distribution of Photon Hits", 1000, -110, 110, 1000, -110, 110);
  
  fhCHreal = new TH1D("fhCHreal","Generated Cherenkov angle, hits", 100, 0.78, 0.88);
  fhCHlam = new TH2D("fhCHlam","Cherenkov angle as a function of lambda", 1000,0.,1000.,100,0.7,0.9);
  
//Histogram for MCPoints in Photon Detector
  fhXYPDMCPt = new TH2D("fhXYPDMCpt", "XY distribution of Photon MCPt", 1000, -110, 110, 1000, -110, 110);
  fhLambdaMC =new TH1D("fhLambdaMC", "No of Photons  vs. Lambda", 200, 100, 1000);
 
  fhCHrealMC = new TH1D("fhCHrealMC","Generated Cherenkov angle, MC",100, 0.78, 0.88);
  fhCHlamMC = new TH2D("fhCHlamMC","Cherenkov angle as a function of lambda", 1000,0.,1000.,100,0.7,0.9);
  
  fhCHlamE = new TH2D("fhCHlamE","CHerenkov angle as a function of lambda", 1000,0.,1000.,100,0.7,0.9);

//Histogram for visualization
  fhXYPDMCPtKp = new TH2D("fhXYPDMCptKp", "XY distribution of Photon MCPt", 1000, -110, 110, 1000, -110, 110);
  fhXYPDMCPtKn = new TH2D("fhXYPDMCptKn", "XY distribution of Photon MCPt", 1000, -110, 110, 1000, -110, 110);
  fhXYPDMCPtpip = new TH2D("fhXYPDMCptpip", "XY distribution of Photon MCPt", 1000, -110, 110, 1000, -110, 110);
  fhXYPDMCPtpin = new TH2D("fhXYPDMCptpin", "XY distribution of Photon MCPt", 1000, -110, 110, 1000, -110, 110);
  fhXYPDMCPtmp = new TH2D("fhXYPDMCptmp", "XY distribution of Photon MCPt", 1000, -110, 110, 1000, -110, 110);
  fhXYPDMCPtmn = new TH2D("fhXYPDMCptmn", "XY distribution of Photon MCPt", 1000, -110, 110, 1000, -110, 110);
  
  fhXYPDHitKp = new TH2D("fhXYPDHitKp", "XY distribution of Photon Hit", 1000, -110, 110, 1000, -110, 110);
  fhXYPDHitKn = new TH2D("fhXYPDHitKn", "XY distribution of Photon Hit", 1000, -110, 110, 1000, -110, 110);
  fhXYPDHitpip = new TH2D("fhXYPDHitpip", "XY distribution of Photon Hit", 1000, -110, 110, 1000, -110, 110);
  fhXYPDHitpin = new TH2D("fhXYPDHitpin", "XY distribution of Photon Hit", 1000, -110, 110, 1000, -110, 110);
  fhXYPDHitmp = new TH2D("fhXYPDHitmp", "XY distribution of Photon Hit", 1000, -110, 110, 1000, -110, 110);
  fhXYPDHitmn = new TH2D("fhXYPDHitmn", "XY distribution of Photon Hit", 1000, -110, 110, 1000, -110, 110);


  fHistoList->Add(fhThetaC);
  fHistoList->Add(fhThetaCMass);
  fHistoList->Add(fhThetaCMomK);
  fHistoList->Add(fhThetaCMomP);
  fHistoList->Add(fhThetaCMomM);
  fHistoList->Add(fhThetaCMomE);
  fHistoList->Add(fhMomAng);
  fHistoList->Add(fhPhoTheta); 
  fHistoList->Add(fhPDPlane);
 
  fHistoList->Add(fhLambda);
  fHistoList->Add(fhLambdaMC);
  fHistoList->Add(fhXYPDHit);
  fHistoList->Add(fhXYPDMCPt);
  fHistoList->Add(fhPDTime);
  fHistoList->Add(fhCHreal);
  fHistoList->Add(fhCHrealMC);
  fHistoList->Add(fhCHlam);
  fHistoList->Add(fhCHlamMC);
  fHistoList->Add(fhCHlamE);  
}

void DrawHits::DrawHisto ( )

private

Definition at line 465 of file DrawHits.cxx.

References cos(), Double_t, fhCHreal, fhCHrealMC, fhLambda, fhLambdaMC, fhMomAng, fhPDPlane, fhPDTime, fhPhoTheta, fhThetaCMass, fhThetaCMomE, fhThetaCMomK, fhThetaCMomM, fhThetaCMomP, fhXYPDHit, fhXYPDHitKn, fhXYPDHitKp, fhXYPDHitmn, fhXYPDHitmp, fhXYPDHitpin, fhXYPDHitpip, fhXYPDMCPt, fhXYPDMCPtKn, fhXYPDMCPtKp, fhXYPDMCPtmn, fhXYPDMCPtmp, fhXYPDMCPtpin, fhXYPDMCPtpip, i, leg, p1, p2, Pi, SetPlotStyle(), sin(), and theta.

Referenced by Finish().

{
  Double_t deg=180.0/TMath::Pi();
   TCanvas *C1= new TCanvas("C1","Cherenkov Angle vs. Mass ",500,500);
   C1->Divide(1,1);
  //Drawing Histo
   C1->cd(1);
  fhThetaCMass->SetMarkerColor(kBlue);
  fhThetaCMass->SetMarkerStyle(20);
  fhThetaCMass->SetMarkerSize(0.8);
  fhThetaCMass->GetYaxis()->SetTitle("Mass (GeV/c^{2})");
  fhThetaCMass->GetYaxis()->SetTitleOffset(1.0);
  fhThetaCMass->GetXaxis()->SetTitle("#theta_{C} (rad)");
  fhThetaCMass->SetTitle(0);        
  fhThetaCMass->SetStats(0);        
  fhThetaCMass->Draw();
   C1->cd();
   TCanvas *C1A= new TCanvas("C1A"," Mom vs. Cherenkov angle",500,500);
   C1A->Divide(1,2);
   C1A->cd(1);
  fhThetaCMomK->SetMarkerColor(kBlue);
  fhThetaCMomP->SetMarkerColor(kMagenta);
  fhThetaCMomM->SetMarkerColor(kGreen);
  fhThetaCMomE->SetMarkerColor(kBlack);
  fhThetaCMomK->SetMarkerStyle(20);
  fhThetaCMomK->SetMarkerSize(0.6);
  fhThetaCMomP->SetMarkerStyle(20);
  fhThetaCMomP->SetMarkerSize(0.6);
  fhThetaCMomM->SetMarkerStyle(20);
  fhThetaCMomM->SetMarkerSize(0.6);
  fhThetaCMomE->SetMarkerStyle(20);
  fhThetaCMomE->SetMarkerSize(0.6);
  fhThetaCMomK->GetYaxis()->SetTitle("#theta_{C} (rad)");
  fhThetaCMomK->GetYaxis()->SetTitleOffset(1.0);
  fhThetaCMomK->GetXaxis()->SetTitle("p (GeV/c)");
  fhThetaCMomK->SetTitle(0);        
  fhThetaCMomK->SetStats(0);        
  fhThetaCMomK->Draw();
  fhThetaCMomM->Draw("same");
  fhThetaCMomP->Draw("same");
  fhThetaCMomE->Draw("same");
   C1A->cd(2);
  fhPhoTheta->Draw(); 
  
  TCanvas *C1B= new TCanvas("C1B"," Momentum vs. track angle",500,500);
   C1B->Divide(1,1);
   C1B->cd();
  fhMomAng->SetMarkerColor(kBlue);
  fhMomAng->SetMarkerStyle(20);
  fhMomAng->SetMarkerSize(0.3);
  fhMomAng->GetXaxis()->SetTitle("#theta (deg)");
  fhMomAng->GetYaxis()->SetTitleOffset(1.0);
  fhMomAng->GetYaxis()->SetTitle("p (GeV/c)");
  fhMomAng->SetTitle(0);            
//  fhMomAng->SetStats(0);          
  fhMomAng->Draw();
  C1B->cd();
  TCanvas *C2= new TCanvas("C2","Photon Distribution" ,700,700);
  C2->Divide(2,3);
  C2->cd(1);
  fhXYPDHit->SetMarkerColor(kBlue);
  fhXYPDHit->SetMarkerStyle(20);
  fhXYPDHit->SetMarkerSize(0.6);
  fhXYPDHit->GetXaxis()->SetTitle("x_{Hit} (cm)");
  fhXYPDHit->GetXaxis()->SetTitleOffset(1.0);
  fhXYPDHit->GetYaxis()->SetTitle("y_{Hit} (cm)");
  fhXYPDHit->SetTitle(0);           
  fhXYPDHit->Draw();
  C2->cd(2);
  fhXYPDMCPt->SetMarkerColor(kBlue);
  fhXYPDMCPt->SetMarkerStyle(20);
  fhXYPDMCPt->SetMarkerSize(0.6);
  fhXYPDMCPt->GetXaxis()->SetTitle("x_{MC} (cm)");
  fhXYPDMCPt->GetXaxis()->SetTitleOffset(1.0);
  fhXYPDMCPt->GetYaxis()->SetTitle("y_{MC} (cm)");
  fhXYPDMCPt->SetTitle(0);          
  fhXYPDMCPt->Draw();
  C2->cd(3);
  fhLambdaMC->SetLineColor(kBlue);
  fhLambdaMC->SetMarkerStyle(20);
  fhLambdaMC->SetMarkerSize(0.01);
  fhLambdaMC->GetXaxis()->SetTitle("#lambda_{MC} (nm)");
  fhLambdaMC->GetXaxis()->SetTitleOffset(1.0);
  fhLambdaMC->GetYaxis()->SetTitle("Counts");
  fhLambdaMC->SetTitle(0);          
//  fhLambdaMC->SetStats(0);        
  fhLambdaMC->Draw();
  C2->cd(4);
  fhLambda->SetLineColor(kRed);
  fhLambda->SetMarkerStyle(20);
  fhLambda->SetMarkerSize(0.01);
  fhLambda->GetXaxis()->SetTitle("#lambda_{Hit} (nm)");
  fhLambda->GetXaxis()->SetTitleOffset(1.0);
  fhLambda->GetYaxis()->SetTitle("Counts");
  fhLambda->SetTitle(0);            
//  fhLambda->SetStats(0);          
  fhLambda->Draw();
  C2->cd(5);  
  fhCHrealMC->Draw();
  fhCHrealMC->Fit("gaus");
  C2->cd(6);  
  fhCHreal->Draw();
  fhCHreal->Fit("gaus");
  C2->cd();
   TCanvas *C3= new TCanvas("C3","Time ",900,500);
   C3->Divide();
  //Drawing Histo
   C3->cd();
  fhPDTime->SetLineColor(kBlack);
  fhPDTime->SetMarkerStyle(20);
  fhPDTime->SetFillColor(2);
  fhPDTime->GetXaxis()->SetTitle("Time (ns)");
  fhPDTime->GetXaxis()->SetTitleOffset(1.0);
  fhPDTime->GetYaxis()->SetTitle("Counts");
  fhPDTime->SetTitle(0);            
  fhPDTime->Draw();         
  C3->cd();
  TCanvas *C4= new TCanvas("C4"," ",500,500);
   C4->Divide(1,1);
   C4->cd(1);
  fhXYPDHitKp->SetMarkerColor(kBlue);
  fhXYPDHitKn->SetMarkerColor(kBlack);
  fhXYPDHitmp->SetMarkerColor(kYellow);
  fhXYPDHitmn->SetMarkerColor(kGreen);
  fhXYPDHitpip->SetMarkerColor(kMagenta);
  fhXYPDHitpin->SetMarkerColor(kRed);
  fhXYPDHitKp->SetMarkerStyle(20);
  fhXYPDHitKp->SetMarkerSize(0.6);
  fhXYPDHitKn->SetMarkerStyle(20);
  fhXYPDHitKn->SetMarkerSize(0.6);
  fhXYPDHitpip->SetMarkerStyle(20);
  fhXYPDHitpin->SetMarkerSize(0.6);
  fhXYPDHitmp->SetMarkerStyle(20);
  fhXYPDHitmn->SetMarkerSize(0.6);
  fhXYPDHitKp->GetXaxis()->SetTitle("x_{Hit} (cm)");
  fhXYPDHitKp->GetXaxis()->SetTitleOffset(1.0);
  fhXYPDHitKp->GetYaxis()->SetTitle("y_{Hit} (cm)");
  fhXYPDHitKp->SetTitle(0);         
  fhXYPDHitKp->Draw();
  fhXYPDHitKn->Draw("same");
  fhXYPDHitpip->Draw("same");
  fhXYPDHitpip->Draw("same");
  fhXYPDHitmp->Draw("same");
  fhXYPDHitmn->Draw("same");
  Double_t theta=11.25/deg;
   Double_t rad_out=104.0; // PD outer radius
   Double_t rad_in=50.53; // PD outer radius
   Double_t rp_out=rad_out/cos(theta); 
   Double_t rp_in=rad_in/cos(theta);
   Double_t xout[17], yout[17], xin[17], yin[17]; 
  for(Int_t i=0; i<16; i++){
     xout[i]=rp_out*cos(2.0*i*theta); 
     yout[i]=rp_out*sin(2.0*i*theta); 
     xin[i]=rp_in*cos(2.0*i*theta); 
     yin[i]=rp_in*sin(2.0*i*theta);
  TLine *l1 =new TLine(xout[i],yout[i],xin[i],yin[i]);
  l1->Draw("same");
  }
     xout[16]=xout[0]; 
     yout[16]=yout[0]; 
     xin[16]=xin[0]; 
     yin[16]=yin[0];
  TPolyLine *p1 =new TPolyLine(17,xout,yout);
  p1->SetFillColor(kYellow);
  p1->Draw("same");
  TPolyLine *p2 =new TPolyLine(17,xin,yin);
  p2->SetFillColor(0);
  p2->Draw("same");
   C4->cd();
  TCanvas *C5= new TCanvas("C5"," ",500,500);
   C5->Divide(1,1);
   C5->cd(1);
  fhXYPDMCPtKp->SetMarkerColor(kBlue);
  fhXYPDMCPtKn->SetMarkerColor(kBlack);
  fhXYPDMCPtmp->SetMarkerColor(kYellow);
  fhXYPDMCPtmn->SetMarkerColor(kGreen);
  fhXYPDMCPtpip->SetMarkerColor(kMagenta);
  fhXYPDMCPtpin->SetMarkerColor(kRed);
  fhXYPDMCPtKp->SetMarkerStyle(20);
  fhXYPDMCPtKp->SetMarkerSize(0.6);
  fhXYPDMCPtKn->SetMarkerStyle(20);
  fhXYPDMCPtKn->SetMarkerSize(0.6);
  fhXYPDMCPtpip->SetMarkerStyle(20);
  fhXYPDMCPtpip->SetMarkerSize(0.6);
  fhXYPDMCPtpin->SetMarkerSize(0.6);
  fhXYPDMCPtpin->SetMarkerSize(0.6);
  fhXYPDMCPtmp->SetMarkerStyle(20);
  fhXYPDMCPtmp->SetMarkerSize(0.6);
  fhXYPDMCPtmn->SetMarkerSize(0.6);
  fhXYPDMCPtmn->SetMarkerSize(0.6);
  fhXYPDMCPtKp->GetXaxis()->SetTitle("x_{MC} (cm)");
  fhXYPDMCPtKp->GetXaxis()->SetTitleOffset(1.0);
  fhXYPDMCPtKp->GetYaxis()->SetTitle("y_{MC} (cm)");
  fhXYPDMCPtKp->SetTitle(0);        
  fhXYPDMCPtKp->Draw();
  fhXYPDMCPtKn->Draw("same");
  fhXYPDMCPtpip->Draw("same");
  fhXYPDMCPtpip->Draw("same");
  fhXYPDMCPtmp->Draw("same");
  fhXYPDMCPtmn->Draw("same");
  for(Int_t i=0; i<16; i++){
  TLine *l2 =new TLine(xout[i],yout[i],xin[i],yin[i]);
  l2->Draw("same");
  }
  p1->Draw("same");
  p2->Draw("same");
 
  SetPlotStyle();
  //TCanvas *C7=  //[R.K.03/2017] unused variable
  new TCanvas("C7","Full sim study",500,500);  
  fhPDPlane->SetMarkerStyle(20);
  fhPDPlane->SetMarkerSize(0.3);
  fhPDPlane->SetMarkerColor(8);
  fhPDPlane->Draw();
  //fhPDPlane->Draw("COL2Z"); 
  for(Int_t i=0; i<16; i++){
  TLine *l2 =new TLine(xout[i],yout[i],xin[i],yin[i]);
  l2->Draw("same");
  }
  p1->Draw("same");
  p2->Draw("same");
  TLegend* leg = new TLegend(0.5,0.5,0.75,0.6); 
  leg->AddEntry(fhPDPlane,"other photons","p");
  leg->SetTextSize(0.03); 
  leg->Draw("same");
  
}

void DrawHits::Exec ( Option_t *  option )

virtual

Executed task

Definition at line 124 of file DrawHits.cxx.

References fDetectorID, nevents, ProcessBarHit(), ProcessPhotonHit(), and ProcessPhotonMC().

{

  nevents++;
  fDetectorID = 0;
   

  // fNHits = 0;
  ProcessBarHit();
  ProcessPhotonHit();
  ProcessPhotonMC();
}

void DrawHits::Finish ( )

virtual

Finish task

Definition at line 457 of file DrawHits.cxx.

References DrawHisto(), and WriteToFile().

{
  cout << "-I- DrawHits: Finish" << endl;

   WriteToFile();
   DrawHisto();
}

InitStatus DrawHits::Init ( )

virtual

Definition at line 70 of file DrawHits.cxx.

References CreateHisto(), fBarPointArray, fHitArray, fMCArray, fPDHitArray, fPDPointArray, and nevents.

{
  cout << " ---------- INITIALIZATION ------------" << endl;
  nevents = 0;
  // Get RootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) {
    cout << "-E- DrawHits::Init: "
         << "RootManager not instantiated!" << endl;
    return kFATAL;
  }

  // Get input array
  fMCArray = (TClonesArray*) ioman->GetObject("MCTrack");
  if ( ! fMCArray ) {
    cout << "-W- DrawHits::Init: "
         << "No MCTrack array!" << endl;
    return kERROR;
  }
  // Get input array
  fBarPointArray = (TClonesArray*) ioman->GetObject("DrcBarPoint");
  if ( ! fBarPointArray ) {
    cout << "-W- DrawHits::Init: "
         << "No DrcBarPoint array!" << endl;
    return kERROR;
  }
  // Get Photon point array
  fPDPointArray = (TClonesArray*) ioman->GetObject("DrcPDPoint");
  if ( ! fPDPointArray ) {
    cout << "-W- DrawHits::Init: "
         << "No DrcPDPoint array!" << endl;
    return kERROR;
  }
  // Get input array
  fHitArray = (TClonesArray*) ioman->GetObject("DrcHit");
  if ( ! fHitArray ) {
    cout << "-W- DrawHits::Init: "
         << "No DrcHit array!" << endl;
    return kERROR;
  }
  // Get input array
  fPDHitArray = (TClonesArray*) ioman->GetObject("DrcPDHit");
  if ( ! fPDHitArray ) {
    cout << "-W- DrawHits::Init: "
         << "No DrcPDHit array!" << endl;
    return kERROR;
  }
    
   cout << "-I- DrawHits: Intialization successfull" << endl;
   CreateHisto();
   return kSUCCESS;
   
}

Bool_t DrawHits::InsideBar ( Double_t  , Double_t  , Double_t  , Double_t  , Double_t  , Double_t  , Double_t  , Double_t  , Double_t  , Double_t    )

private

void DrawHits::ProcessBarHit ( )

private

Definition at line 137 of file DrawHits.cxx.

References Double_t, fBarPointArray, fErrThetaC, fHitArray, fhMomAng, fhThetaC, fhThetaCMass, fhThetaCMomE, fhThetaCMomK, fhThetaCMomM, fhThetaCMomP, fMCArray, fThetaC, fVerbose, PndDrcBarPoint::GetMass(), PndMCTrack::GetMotherID(), PndMCTrack::GetPdgCode(), PndDrcHit::GetRefIndex(), PndDrcBarPoint::GetThetaC(), hit(), nevents, Pi, printf(), pt(), sqrt(), tr, xbar, xEnt, ybar, yEnt, zbar, and zEnt.

Referenced by Exec().

{
  PndDrcBarPoint* pt=NULL;
  PndDrcHit* hit=NULL;
  PndMCTrack* tr = NULL;
  Double_t deg=180.0/TMath::Pi();

 // Loop over PndDrcHit
  for(Int_t j=0; j<fHitArray->GetEntriesFast(); j++) {
 //for(Int_t j=0; j<fBarPointArray->GetEntriesFast(); j++) {
       if (fVerbose > 1) printf("\n\n=====> Event No. %d\n", nevents); 

       hit = (PndDrcHit*)fHitArray->At(j);

       Double_t xH = hit->GetX();
       Double_t yH = hit->GetY();
       Double_t zH = hit->GetZ();
       
       Int_t mcRef= hit->GetRefIndex();
       pt= (PndDrcBarPoint*)fBarPointArray->At(mcRef);
  //     pt= (PndDrcBarPoint*)fBarPointArray->At(j);
       Int_t chtrID= pt->GetTrackID();
       tr = (PndMCTrack*)fMCArray->At(chtrID);
       Int_t chtrPdg= tr->GetPdgCode();
       Int_t chtrMid=tr->GetMotherID();
       
       Double_t trPx = pt->GetPx();
       Double_t trPy = pt->GetPy();
       Double_t trPz = pt->GetPz();
       Double_t trP = sqrt(trPx*trPx + trPy*trPy +trPz*trPz);
       Double_t mass = pt->GetMass();
       //Double_t energy = sqrt(trP*trP + mass*mass); //[R.K. 01/2017] unused variable?
       
       // for my tree phoTree:
       xEnt = pt->GetX();
       yEnt = pt->GetY();
       zEnt = pt->GetZ();
       
       xbar = xH;
       ybar = yH;
       zbar = zH;
              
       //phoTree->Fill();
       TVector3 barmom;
       pt->Momentum(barmom);
       Double_t angIn= barmom.Theta();
       Double_t angDeg=angIn*deg;
   
       fThetaC = gRandom->Gaus(pt->GetThetaC(),0.008);
       fErrThetaC = 0.; //rad
       
         if(abs(chtrPdg) == 321)fhThetaCMomK->Fill(trP,fThetaC);
         if(abs(chtrPdg) == 211)fhThetaCMomP->Fill(trP,fThetaC);
         if(abs(chtrPdg) == 13)fhThetaCMomM->Fill(trP,fThetaC);
         if(abs(chtrPdg) == 11)fhThetaCMomE->Fill(trP,fThetaC);
       if(chtrMid == -1){
         fhMomAng->Fill(angDeg,trP);
         cout<<"ang = "<<angDeg<<", mom = "<<trP<<endl;

       }//prim particle
       fhThetaC->Fill(fThetaC);
       fhThetaCMass->Fill(fThetaC,mass);

    }
}

void DrawHits::ProcessPhotonHit ( )

private

Definition at line 257 of file DrawHits.cxx.

References Double_t, fhCHlam, fhCHreal, fhLambda, fhPDTime, fhPhoTheta, fhXYPDHit, fhXYPDHitKn, fhXYPDHitKp, fhXYPDHitmn, fhXYPDHitmp, fhXYPDHitpin, fhXYPDHitpip, fMCArray, fPDHitArray, fPDPointArray, fPMo, fPx, fPy, fPz, fStartVertex, fVerbose, fXcross, fYcross, PndMCTrack::GetMomentum(), PndMCTrack::GetMotherID(), PndMCTrack::GetPdgCode(), PndDrcPDHit::GetRefIndex(), PndMCTrack::GetStartVertex(), PndDrcPDHit::GetTime(), nevents, Pi, pxMo, pxPho, pyMo, pyPho, pzMo, pzPho, sqrt(), thit, tr, xhit, yhit, and zhit.

Referenced by Exec().

{
  PndDrcPDPoint* Ppt=NULL;
  PndDrcPDHit* pdhit=NULL;
  PndMCTrack* tr = NULL;
  PndMCTrack* trMr;
  if (fVerbose > 0) {
    cout <<" Number of Detected Photon Hits : "<<fPDHitArray->GetEntries()<<endl;
  }

// Loop over PndDrcPDHits
  for(Int_t k=0; k<fPDHitArray->GetEntriesFast(); k++) {
   
    pdhit = (PndDrcPDHit*)fPDHitArray->At(k);
    
    Int_t mcPDRef= pdhit->GetRefIndex();
    Ppt = (PndDrcPDPoint*)fPDPointArray->At(mcPDRef);
    
    Int_t trID= Ppt->GetTrackID();
    tr = (PndMCTrack*)fMCArray->At(trID);
   
    Int_t trMID=tr->GetMotherID();

    trMr = (PndMCTrack*)fMCArray->At(trMID);
    Int_t trMpdg=trMr->GetPdgCode();
  
    if(trMr->GetMotherID()==-1){
    

      Double_t xPHit= pdhit->GetX();
      Double_t yPHit= pdhit->GetY();
      Double_t zPHit= pdhit->GetZ();
      Double_t time = pdhit->GetTime();
      
      // for my tree:
      xhit = xPHit;
      yhit = yPHit;
      zhit = zPHit;
      thit = time;
      
      pxMo = trMr->GetMomentum().X();
      pyMo = trMr->GetMomentum().Y();
      pzMo = trMr->GetMomentum().Z();
      fPMo.SetXYZ(pxMo, pyMo, pzMo);
      
      fhPhoTheta->Fill(fPMo.Theta()/3.1415*180.);
     
      pxPho = tr->GetMomentum().X(); // initial momentum of a photon
      pyPho = tr->GetMomentum().Y();
      pzPho = tr->GetMomentum().Z();
      
      fStartVertex = tr->GetStartVertex();
      
      // momentum of a photon on the PD Plane
      fPx= Ppt->GetPx();
      fPy= Ppt->GetPy();
      fPz= Ppt->GetPz();
            
      // extrapolation of momentum on the PD Plane to the bar:
      fXcross = xPHit - fPx/fPz*(120. + zPHit);
      fYcross = yPHit - fPy/fPz*(120. + zPHit);
     
      //phoTree->Fill();
      // ------------
      
      //Double_t xP= Ppt->GetX(); //[R.K. 01/2017] unused variable?
      //Double_t yP= Ppt->GetY(); //[R.K. 01/2017] unused variable?
      //Double_t zP= Ppt->GetZ(); //[R.K. 01/2017] unused variable?
      
      Double_t PPx= Ppt->GetPx();
      Double_t PPy= Ppt->GetPy();
      Double_t PPz= Ppt->GetPz();
      Double_t etot = sqrt(PPx*PPx + PPy*PPy +PPz*PPz);
      //Double_t nRefrac=1.467; //[R.K. 01/2017] unused variable?
      //Double_t lambdah=197.0*2.0*TMath::Pi()/nRefrac/(etot*1.0E9);//wavelength of photon in nm
      Double_t lambdah=197.0*2.0*TMath::Pi()/(etot*1.0E9);
      
      fhPDTime->Fill(time);
      //if(lambdah > 350. && lambdah < 600.){ 
      fhLambda->Fill(lambdah);
      TVector3 photon1;
      TVector3 mother1;
      photon1.SetXYZ(tr->GetMomentum().X(), tr->GetMomentum().Y(), tr->GetMomentum().Z());
      mother1.SetXYZ(trMr->GetMomentum().X(), trMr->GetMomentum().Y(), trMr->GetMomentum().Z());
      fhCHreal->Fill(photon1.Angle(mother1));
      fhCHlam->Fill(lambdah, photon1.Angle(mother1));
      //}
    
      fhXYPDHit->Fill(xPHit,yPHit);
            
      if( nevents == 1 && trMpdg ==321)fhXYPDHitKp->Fill(xPHit,yPHit);
      if( nevents == 1 && trMpdg ==-321)fhXYPDHitKn->Fill(xPHit,yPHit);
      if( nevents == 1 && trMpdg ==211)fhXYPDHitpip->Fill(xPHit,yPHit);
      if( nevents == 1 && trMpdg ==-211)fhXYPDHitpin->Fill(xPHit,yPHit);
      if( nevents == 1 && trMpdg ==13)fhXYPDHitmp->Fill(xPHit,yPHit);
      if( nevents == 1 && trMpdg ==-13)fhXYPDHitmn->Fill(xPHit,yPHit);
   }// photon from primary particle 
 }// photon hits
}

void DrawHits::ProcessPhotonMC ( )

private

Definition at line 203 of file DrawHits.cxx.

References Double_t, fhCHlamMC, fhCHrealMC, fhLambdaMC, fhXYPDMCPt, fhXYPDMCPtKn, fhXYPDMCPtKp, fhXYPDMCPtmn, fhXYPDMCPtmp, fhXYPDMCPtpin, fhXYPDMCPtpip, fMCArray, fPDPointArray, PndMCTrack::GetMomentum(), PndMCTrack::GetMotherID(), PndMCTrack::GetPdgCode(), lambda(), nevents, Pi, and sqrt().

Referenced by Exec().

{
  //Double_t deg=180.0/TMath::Pi(); //[R.K. 01/2017] unused variable?


   PndMCTrack* ptr = NULL;
   PndMCTrack* trMr=NULL;
   for(Int_t l=0; l<fPDPointArray->GetEntriesFast(); l++) {
      PndDrcPDPoint *Ppt = (PndDrcPDPoint*)fPDPointArray->At(l);
      Int_t trID= Ppt->GetTrackID();
      ptr = (PndMCTrack*)fMCArray->At(trID);
    
      Int_t trMID=ptr->GetMotherID();
      trMr = (PndMCTrack*)fMCArray->At(trMID);
      Int_t trMpdg=trMr->GetPdgCode();

      if(trMr->GetMotherID()==-1){
         
         Double_t xP= Ppt->GetX();
         Double_t yP= Ppt->GetY();
         //Double_t zP= Ppt->GetZ(); //[R.K. 01/2017] unused variable?

         fhXYPDMCPt->Fill(xP,yP);
        if( nevents==1 && trMpdg ==321)fhXYPDMCPtKp->Fill(xP,yP);
        if( nevents==1 && trMpdg ==-321)fhXYPDMCPtKn->Fill(xP,yP);
        if( nevents==1 && trMpdg ==211)fhXYPDMCPtpip->Fill(xP,yP);
        if( nevents==1 && trMpdg ==-211)fhXYPDMCPtpin->Fill(xP,yP);
        if( nevents==1 && trMpdg ==13)fhXYPDMCPtmp->Fill(xP,yP);
        if( nevents ==1 && trMpdg ==-13)fhXYPDMCPtmn->Fill(xP,yP);
            
        // Momentum of photon at PD   point
         Double_t Pfx= Ppt->GetPx();
         Double_t Pfy= Ppt->GetPy();
         Double_t Pfz= Ppt->GetPz();
         //Double_t Pf = sqrt(Pfx*Pfx + Pfy*Pfy +Pfz*Pfz); //[R.K. 01/2017] unused variable?
         Double_t etot = sqrt(Pfx*Pfx + Pfy*Pfy +Pfz*Pfz);
         //Double_t nRefrac=1.467; //[R.K. 01/2017] unused variable?
         //Double_t lambda=197.0*2.0*TMath::Pi()/nRefrac/(etot*1.0E9);//wavelength of photon in nm
         Double_t lambda=197.0*2.0*TMath::Pi()/(etot*1.0E9);

       //if(lambda > 350. && lambda < 600.){ 
       fhLambdaMC->Fill(lambda);
       TVector3 photon;
       TVector3 mother;
       photon.SetXYZ(ptr->GetMomentum().X(), ptr->GetMomentum().Y(), ptr->GetMomentum().Z());
       mother.SetXYZ(trMr->GetMomentum().X(), trMr->GetMomentum().Y(), trMr->GetMomentum().Z());              
       fhCHrealMC->Fill(photon.Angle(mother));
       fhCHlamMC->Fill(lambda,mother.Angle(photon));
       //}
     }// Prim particle
  }//Photon points

}  

void DrawHits::RotBarCoordinate	(	TVector3 &	vector,
		Int_t	barID
	)

Rotation in Bar Co-ordinate System

void DrawHits::SetDefaultParameters ( )

private

Set the parameters to the default values.

void DrawHits::SetPlotStyle ( )

inlineprivate

Definition at line 107 of file DrawHits.h.

References Double_t.

Referenced by DrawHisto().

 {
    const Int_t NRGBs = 5;
    const Int_t NCont = 255;

    Double_t stops[NRGBs] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
    Double_t red[NRGBs]   = { 0.00, 0.00, 0.87, 1.00, 0.51 };
    Double_t green[NRGBs] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
    Double_t blue[NRGBs]  = { 0.51, 1.00, 0.12, 0.00, 0.00 };
    TColor::CreateGradientColorTable(NRGBs, stops, red, green, blue, NCont);
    gStyle->SetNumberContours(NCont);
 }

void DrawHits::SetTreeName ( TString  str )

inline

Definition at line 71 of file DrawHits.h.

References fTreeName.

{fTreeName = str;}

void DrawHits::WriteToFile ( )

private

Definition at line 450 of file DrawHits.cxx.

References fHistoList, next, and phoTree.

Referenced by Finish().

{
              phoTree->Write();
              TIter next(fHistoList);
                           while ( TH1* histo = ((TH1*)next()) ) histo->Write();
}

Member Data Documentation

TH1F* DrawHits::currh

private

Definition at line 164 of file DrawHits.h.

TClonesArray* DrawHits::fBarPointArray

private

Definition at line 89 of file DrawHits.h.

Referenced by Init(), and ProcessBarHit().

Int_t DrawHits::fDetectorID

private

Definition at line 84 of file DrawHits.h.

Referenced by Exec().

Double_t DrawHits::fErrThetaC

private

Definition at line 86 of file DrawHits.h.

Referenced by ProcessBarHit().

PndGeoDrc* DrawHits::fGeo

private

Basic geometry data of barrel DRC.

Definition at line 98 of file DrawHits.h.

Referenced by DrawHits(), and ~DrawHits().

TH2D* DrawHits::fhCHlam

private

Definition at line 181 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

TH2D* DrawHits::fhCHlamE

private

Definition at line 183 of file DrawHits.h.

Referenced by CreateHisto().

TH2D* DrawHits::fhCHlamMC

private

Definition at line 182 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonMC().

TH1D* DrawHits::fhCHreal

private

Definition at line 180 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH1D* DrawHits::fhCHrealMC

private

Definition at line 179 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

TH2F* DrawHits::fhHitsD

private

Definition at line 186 of file DrawHits.h.

TList* DrawHits::fHistoList

private

Definition at line 94 of file DrawHits.h.

Referenced by CreateHisto(), WriteToFile(), and ~DrawHits().

TClonesArray* DrawHits::fHitArray

private

Definition at line 91 of file DrawHits.h.

Referenced by Init(), and ProcessBarHit().

TH1D* DrawHits::fhLambda

private

Definition at line 188 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH1D* DrawHits::fhLambdaMC

private

Definition at line 190 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

TH2D* DrawHits::fhMomAng

private

Definition at line 171 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessBarHit().

TH2F* DrawHits::fhPDPlane

private

Definition at line 178 of file DrawHits.h.

Referenced by CreateHisto(), and DrawHisto().

TH1D* DrawHits::fhPDTime

private

Definition at line 192 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH1F* DrawHits::fhPhoTheta

private

Definition at line 177 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH1D* DrawHits::fhThetaC

private

Definition at line 169 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessBarHit().

TH2D* DrawHits::fhThetaCMass

private

Definition at line 170 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessBarHit().

TH2D* DrawHits::fhThetaCMomE

private

Definition at line 175 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessBarHit().

TH2D* DrawHits::fhThetaCMomK

private

Definition at line 172 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessBarHit().

TH2D* DrawHits::fhThetaCMomM

private

Definition at line 174 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessBarHit().

TH2D* DrawHits::fhThetaCMomP

private

Definition at line 173 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessBarHit().

TH2D* DrawHits::fhXYPDHit

private

Definition at line 189 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH2D* DrawHits::fhXYPDHitKn

private

Definition at line 195 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH2D* DrawHits::fhXYPDHitKp

private

Definition at line 194 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH2D* DrawHits::fhXYPDHitmn

private

Definition at line 199 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH2D* DrawHits::fhXYPDHitmp

private

Definition at line 198 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH2D* DrawHits::fhXYPDHitpin

private

Definition at line 197 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH2D* DrawHits::fhXYPDHitpip

private

Definition at line 196 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonHit().

TH2D* DrawHits::fhXYPDMCPt

private

Definition at line 191 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

TH2D* DrawHits::fhXYPDMCPtKn

private

Definition at line 202 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

TH2D* DrawHits::fhXYPDMCPtKp

private

Definition at line 201 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

TH2D* DrawHits::fhXYPDMCPtmn

private

Definition at line 206 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

TH2D* DrawHits::fhXYPDMCPtmp

private

Definition at line 205 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

TH2D* DrawHits::fhXYPDMCPtpin

private

Definition at line 204 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

TH2D* DrawHits::fhXYPDMCPtpip

private

Definition at line 203 of file DrawHits.h.

Referenced by CreateHisto(), DrawHisto(), and ProcessPhotonMC().

Double_t DrawHits::fkxBar

private

Definition at line 157 of file DrawHits.h.

Double_t DrawHits::fkyBar

private

Definition at line 158 of file DrawHits.h.

Double_t DrawHits::fkzBar

private

Definition at line 159 of file DrawHits.h.

Double_t DrawHits::fMass

private

Definition at line 85 of file DrawHits.h.

TClonesArray* DrawHits::fMCArray

private

Definition at line 93 of file DrawHits.h.

Referenced by Init(), ProcessBarHit(), ProcessPhotonHit(), and ProcessPhotonMC().

PndGeoDrcPar* DrawHits::fPar

private

Definition at line 96 of file DrawHits.h.

TClonesArray* DrawHits::fPDHitArray

private

Definition at line 92 of file DrawHits.h.

Referenced by Init(), and ProcessPhotonHit().

TClonesArray* DrawHits::fPDPointArray

private

Definition at line 90 of file DrawHits.h.

Referenced by Init(), ProcessPhotonHit(), and ProcessPhotonMC().

Double_t DrawHits::fPhiRot

private

Definition at line 149 of file DrawHits.h.

TVector3 DrawHits::fPMo

private

Definition at line 167 of file DrawHits.h.

Referenced by ProcessPhotonHit().

TVector3 DrawHits::fPphoB

private

Definition at line 156 of file DrawHits.h.

Double_t DrawHits::fPx

private

Definition at line 151 of file DrawHits.h.

Referenced by ProcessPhotonHit().

Double_t DrawHits::fPy

private

Definition at line 152 of file DrawHits.h.

Referenced by ProcessPhotonHit().

Double_t DrawHits::fPz

private

Definition at line 153 of file DrawHits.h.

Referenced by ProcessPhotonHit().

Int_t DrawHits::fRefIndex

private

Definition at line 87 of file DrawHits.h.

TVector3 DrawHits::fStartVertex

private

Definition at line 150 of file DrawHits.h.

Referenced by ProcessPhotonHit().

Double_t DrawHits::fThetaC

private

Definition at line 86 of file DrawHits.h.

Referenced by ProcessBarHit().

TString DrawHits::fTreeName

private

Definition at line 161 of file DrawHits.h.

Referenced by CreateHisto(), and SetTreeName().

Int_t DrawHits::fVerbose

private

Verbosity level

Definition at line 104 of file DrawHits.h.

Referenced by DrawHits(), ProcessBarHit(), and ProcessPhotonHit().

Double_t DrawHits::fXcross

private

Definition at line 154 of file DrawHits.h.

Referenced by ProcessPhotonHit().

Double_t DrawHits::fYcross

private

Definition at line 155 of file DrawHits.h.

Referenced by ProcessPhotonHit().

std::vector<TH1F*> DrawHits::histos

private

Definition at line 163 of file DrawHits.h.

Double_t DrawHits::lam

private

Definition at line 147 of file DrawHits.h.

TGraph2D* DrawHits::map

private

Definition at line 165 of file DrawHits.h.

Int_t DrawHits::nevents

private

Definition at line 122 of file DrawHits.h.

Referenced by Exec(), Init(), ProcessBarHit(), ProcessPhotonHit(), and ProcessPhotonMC().

TTree* DrawHits::phoTree

private

Definition at line 162 of file DrawHits.h.

Referenced by CreateHisto(), and WriteToFile().

TGraph2D* DrawHits::pop

private

Definition at line 166 of file DrawHits.h.

Double_t DrawHits::pxmc

private

Definition at line 128 of file DrawHits.h.

Double_t DrawHits::pxMo

private

Definition at line 138 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::pxPho

private

Definition at line 141 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::pymc

private

Definition at line 129 of file DrawHits.h.

Double_t DrawHits::pyMo

private

Definition at line 139 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::pyPho

private

Definition at line 142 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::pzmc

private

Definition at line 130 of file DrawHits.h.

Double_t DrawHits::pzMo

private

Definition at line 140 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::pzPho

private

Definition at line 143 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::thit

private

Definition at line 134 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::xbar

private

Definition at line 135 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessBarHit().

Double_t DrawHits::xEnt

private

Definition at line 144 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessBarHit().

Float_t DrawHits::xhit

private

Definition at line 131 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::xmcp

private

Definition at line 125 of file DrawHits.h.

Double_t DrawHits::ybar

private

Definition at line 136 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessBarHit().

Double_t DrawHits::yEnt

private

Definition at line 145 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessBarHit().

Float_t DrawHits::yhit

private

Definition at line 132 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::ymcp

private

Definition at line 126 of file DrawHits.h.

Double_t DrawHits::zbar

private

Definition at line 137 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessBarHit().

Double_t DrawHits::zEnt

private

Definition at line 146 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessBarHit().

Double_t DrawHits::zhit

private

Definition at line 133 of file DrawHits.h.

Referenced by CreateHisto(), and ProcessPhotonHit().

Double_t DrawHits::zmcp

private

Definition at line 127 of file DrawHits.h.

---

The documentation for this class was generated from the following files:

• DrawHits.h
• DrawHits.cxx