#include <PndMvdAllDataEventAna.h>

Inheritance diagram for PndMvdAllDataEventAna:

Public Member Functions
	PndMvdAllDataEventAna ()

	PndMvdAllDataEventAna (TString fileName)

virtual	~PndMvdAllDataEventAna ()

virtual void	Init (TString filename)

virtual void	InitBranch ()

virtual void	InitHistos ()

virtual void	AnaHits ()

virtual void	BeginOfEventAction ()

virtual void	EndOfEventAction ()

virtual void	ClearHistos ()

void	PrintHitArray ()

void	PrintDigiArray ()

void	PrintClusterArray ()

void	PrintRecoArray ()

void	FillHitHistos ()

void	FillDigiHistos ()

void	FillClusterHistos ()

void	FillRecoHistos ()

void	FillHitPerClusterHistos ()

void	FillHitResolutionHistos ()

void	Fill3DHisto ()

void	FillHitProjHistos ()

void	DrawHitHisto (TString detName, TCanvas *extCan=0, Int_t pad=0)

void	DrawDigiHisto (TString detName, TCanvas *extCan=0, Int_t pad=0)

void	DrawClusterHisto (TString detName, TCanvas *extCan=0, Int_t pad=0)

void	DrawRecoHisto (TString detName, TCanvas *extCan=0, Int_t pad=0)

void	DrawAllHistos (TString detName, TCanvas *extCan=0)

void	DrawAllHistos (Int_t i, TCanvas *extCan=0)

void	Draw3D (TString opt="", TCanvas *extCan=0, Int_t pad=0)

void	DrawAllTracks (TCanvas *extCan=0, Int_t pad=0)

void	DrawHitTracks (TCanvas *extCan=0, Int_t pad=0)

void	DrawEvent (bool tracks=false, TCanvas *extCan=0)

void	DrawTopVolume (TCanvas extCan=0, Int_t pad=0, const char opt="")

void	DrawHisxy (TCanvas *extCan=0, Int_t pad=0)

void	DrawHisrz (TCanvas *extCan=0, Int_t pad=0)

void	DrawHistoVec (std::vector< TH1 * > vec, TCanvas extCan=0, Int_t pad=0) const

void	DrawHisto (TH1 histo, TCanvas extCan=0, Int_t pad=0) const

void	DrawHitPerCluster (TCanvas *extCan=0, Int_t pad=0)

void	DrawResolution (TCanvas *extCan=0, Int_t pad=0)

void	Create3DGeoHits ()

void	ClearHistoMaps (std::map< TString, TH1 * > *myHistos) const

void	ClearAllHMaps ()

void	ClearHistoVector (std::vector< TH1 * > *myVectors) const

void	ClearAllVectors ()

TVector3	GetLocalHitPoints (TString detName, TVector3 input)

std::vector< Int_t >	GetHitPerCluster (PndSdsCluster *clusterCand)

std::vector< TString >	GetModulesHit ()

TVector3	CalcMeanHitPos (std::vector< Int_t > points)

void	InitCanvas ()

void	AnaEvents (Int_t count)

void	AnaEventNr (Int_t Nr)

void	NextEvent ()

void	DrawHistos ()

void	DrawHisto (TString path, TString opt="")

void	SetCanvasColumns (Int_t col)

void	SetDrawOption (TString histoName, TString option)

void	SetSelection (TString sel)

void	FillVector ()

void	DrawVector (Int_t index, TString opt)

Int_t	GetEventNr () const

TString	GetSelection () const

TH1 *	GetHisto (TString histoName)

TH1 *	GetHisto (Int_t i)

virtual void	AnaHitNr (Int_t i)

Protected Member Functions
	ClassDef (PndMvdEventAna, 1)

Protected Attributes
TClonesArray *	fDigiArray

TClonesArray *	fClusterArray

TClonesArray *	fRecoArray

TClonesArray *	fGeoTrackArray

TClonesArray *	fTrackFArray

Int_t	fActiveEvent

Int_t	fAnaEvents

TFile *	fFile

TTree *	fTree

TClonesArray *	fHitArray

TString	fClassName

TString	fBranchName

TString	fSelection

TCanvas *	fCan

std::map< TString, TH1 * >	fHistos

std::map< TString, TString >	fDrawOption

std::vector< TH1 * >	fHistoVector

Int_t	fCancol

PndGeoHandling *	fGeoH

Private Member Functions
	ClassDef (PndMvdAllDataEventAna, 1)

Private Attributes
std::map< TString, TH1 * >	fDigiHistos

std::map< TString, TH1 * >	fClusterHistos

std::map< TString, TH1 * >	fRecoHistos

std::map< TString, TH1 * >	fHitPerClusterHistos

std::map< TString, TH1 * >	fHitResolutionHistos

std::vector< TH1 * >	fHitHistoVec

std::vector< TH1 * >	fDigiHistoVec

std::vector< TH1 * >	fClusterHistoVec

std::vector< TH1 * >	fRecoHistoVec

TGeoVolume *	fMvdTopVolume

TGeoVolume *	fRecoTopVolume

TGeoVolume *	fRecoVolume

TH1 *	fRecoErrorHisto

TH1 *	fAllHitPerClusterHistos

TH1 *	fAllHitResolutionHistos

TH3D *	f3DMCHisto

TH3D *	f3DRecoHisto

std::vector< TH1 * >	fSimHisxy

std::vector< TH1 * >	fSimHisrz

std::vector< TH1 * >	fRecoHisxy

std::vector< TH1 * >	fRecoHisrz

TCanvas *	fCan1

TCanvas *	fCan2

PndSdsCalcFePixel *	fPixelCon

PndEventDisplay *	fGeoList

Detailed Description

PndMvdAllDataEventAna.h

Author: t.stockmanns t.sto.nosp@m.ckma.nosp@m.nns@f.nosp@m.z-ju.nosp@m.elich.nosp@m..de Analysis class which loads in a complete simulation chain

Definition at line 39 of file PndMvdAllDataEventAna.h.

Constructor & Destructor Documentation

PndMvdAllDataEventAna::PndMvdAllDataEventAna ( )

inline

Definition at line 41 of file PndMvdAllDataEventAna.h.

41 {};

PndMvdAllDataEventAna::PndMvdAllDataEventAna ( TString fileName )

Definition at line 17 of file PndMvdAllDataEventAna.cxx.

References PndEventDisplay::AddNewGroup(), PndMvdEventAna::fGeoH, fGeoList, fMvdTopVolume, fPixelCon, fRecoVolume, gGeoManager, Init(), PndGeoHandling::Instance(), and PndMvdEventAna::SetCanvasColumns().

                                                             :PndMvdEventAna(fileName)
 {
   Init(fileName);
   SetCanvasColumns(6);
 
 
   fGeoH = PndGeoHandling::Instance();
 
   fRecoVolume = gGeoManager->MakeSphere("RecoHit",gGeoManager->GetMedium("vacuum"),0,0.1);
   fRecoVolume->SetLineColor(kRed);
 
   TGeoVolume* mcHit = gGeoManager->MakeSphere("MCHit",gGeoManager->GetMedium("vacuum"),0,0.1);
   mcHit->SetLineColor(kGreen);
 
   TGeoShape* stripShape = gGeoManager->GetVolume("StripSensor")->GetShape();
   TGeoVolume* stripVol = new TGeoVolume("StripSensorRed",stripShape, gGeoManager->GetMedium("vacuum"));
   stripVol->SetLineColor(kRed);
   gGeoManager->AddVolume(stripVol);
   stripVol = gGeoManager->GetVolume("StripSensorRed");
 
   TGeoShape* pixShape1 = gGeoManager->GetVolume("SensorActiveAreao[NrFE=8ooE]")->GetShape();
   TGeoVolume* pixVol1 = new TGeoVolume("PixelSensor8",pixShape1, gGeoManager->GetMedium("vacuum"));
   pixVol1->SetLineColor(kRed);
   gGeoManager->AddVolume(pixVol1);
   pixVol1 = gGeoManager->GetVolume("PixelSensor8");
 
   TGeoShape* pixShape2 = gGeoManager->GetVolume("SensorActiveAreao[NrFE=4]")->GetShape();
   TGeoVolume* pixVol2 = new TGeoVolume("PixelSensor4",pixShape2, gGeoManager->GetMedium("vacuum"));
   pixVol2->SetLineColor(kRed);
   gGeoManager->AddVolume(pixVol2);
   pixVol2 = gGeoManager->GetVolume("PixelSensor4");
 
   TGeoShape* pixShape3 = gGeoManager->GetVolume("SensorActiveAreao[NrFE=5]")->GetShape();
   TGeoVolume* pixVol3 = new TGeoVolume("PixelSensor5",pixShape3, gGeoManager->GetMedium("vacuum"));
   pixVol3->SetLineColor(kRed);
   gGeoManager->AddVolume(pixVol3);
   pixVol3 = gGeoManager->GetVolume("PixelSensor5");
 
   TGeoVolume* stripVolG = new TGeoVolume("StripSensorGreen",stripShape, gGeoManager->GetMedium("vacuum"));
   stripVolG->SetLineColor(kGreen);
   gGeoManager->AddVolume(stripVolG);
   stripVolG = gGeoManager->GetVolume("StripSensorGreen");
 
   TGeoVolume* pixVolG1 = new TGeoVolume("PixelSensor8Green",pixShape1, gGeoManager->GetMedium("vacuum"));
   pixVolG1->SetLineColor(kGreen);
   gGeoManager->AddVolume(pixVolG1);
   pixVolG1 = gGeoManager->GetVolume("PixelSensor8Green");
 
   TGeoVolume* pixVolG2 = new TGeoVolume("PixelSensor4Green",pixShape2, gGeoManager->GetMedium("vacuum"));
   pixVolG2->SetLineColor(kGreen);
   gGeoManager->AddVolume(pixVolG2);
   pixVolG2 = gGeoManager->GetVolume("PixelSensor4Green");
 
   TGeoVolume* pixVolG3 = new TGeoVolume("PixelSensor5Green",pixShape3, gGeoManager->GetMedium("vacuum"));
   pixVolG3->SetLineColor(kGreen);
   gGeoManager->AddVolume(pixVolG3);
   pixVolG3 = gGeoManager->GetVolume("PixelSensor5Green");
 
   fMvdTopVolume = gGeoManager->GetVolume("MVDoOption1");
   fGeoList = new PndEventDisplay();
 
   fGeoList->AddNewGroup("RecoHits",new PndGeoHitList("RecoHits","MVDoOption1", fRecoVolume));
   fGeoList->AddNewGroup("MCHits",new PndGeoHitList("MCHits","MVDoOption1", mcHit));
   fGeoList->AddNewGroup("StripHits", new PndGeoHitList("StripHits","MVDoOption1",stripVol));
   fGeoList->AddNewGroup("PixHits8", new PndGeoHitList("PixHits8","MVDoOption1",pixVol1));
   fGeoList->AddNewGroup("PixHits4", new PndGeoHitList("PixHits4","MVDoOption1",pixVol2));
   fGeoList->AddNewGroup("PixHits5", new PndGeoHitList("PixHits5","MVDoOption1",pixVol3));
   fGeoList->AddNewGroup("StripHitsGreen", new PndGeoHitList("StripHitsGreen","MVDoOption1",stripVolG));
   fGeoList->AddNewGroup("PixHits8Green", new PndGeoHitList("PixHits8Green","MVDoOption1",pixVolG1));
   fGeoList->AddNewGroup("PixHits4Green", new PndGeoHitList("PixHits4Green","MVDoOption1",pixVolG2));
   fGeoList->AddNewGroup("PixHits5Green", new PndGeoHitList("PixHits5Green","MVDoOption1",pixVolG3));
 
   fPixelCon = new PndSdsCalcFePixel(76,84,10);
 }

virtual PndMvdAllDataEventAna::~PndMvdAllDataEventAna ( )

inlinevirtual

Definition at line 43 of file PndMvdAllDataEventAna.h.

43 {};

Member Function Documentation

void PndMvdEventAna::AnaEventNr ( Int_t Nr )

inherited

Definition at line 28 of file PndMvdEventAna.cxx.

References PndMvdEventAna::AnaHits(), PndMvdEventAna::BeginOfEventAction(), PndMvdEventAna::ClearHistos(), PndMvdEventAna::EndOfEventAction(), PndMvdEventAna::fActiveEvent, PndMvdEventAna::fAnaEvents, PndMvdEventAna::fTree, PndMvdEventAna::InitHistos(), and vecs.

Referenced by PndMvdEventAna::NextEvent().

 {
         TVector3 vecs;
         BeginOfEventAction();
 
         fAnaEvents = 1;
         InitHistos();
         ClearHistos();
         if (Nr > fTree->GetEntriesFast())
                 return;
 
         fActiveEvent = Nr;
         AnaHits();
         EndOfEventAction();
         //InitCanvas();
         //DrawHistos();
 }

void PndMvdEventAna::AnaEvents ( Int_t count )

inherited

Definition at line 79 of file PndMvdEventAna.cxx.

References PndMvdEventAna::AnaHits(), count, PndMvdEventAna::fActiveEvent, PndMvdEventAna::fAnaEvents, PndMvdEventAna::fTree, i, and PndMvdEventAna::InitHistos().

 {
         fAnaEvents = count;
         InitHistos();
         if (count > fTree->GetEntriesFast())
                 count = fTree->GetEntriesFast();
         for (Int_t i = 0; i < count; i++)
         {
                 fActiveEvent = i;
                 AnaHits();
         }
         //InitCanvas();
         //DrawHistos();
 }

virtual void PndMvdEventAna::AnaHitNr ( Int_t i )

inlinevirtualinherited

Reimplemented in PndMvdDigiEventAna, and PndMvdMCEventAna.

Definition at line 59 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::AnaHits().

59 {};

void PndMvdAllDataEventAna::AnaHits ( )

virtual

Reimplemented from PndMvdEventAna.

Definition at line 149 of file PndMvdAllDataEventAna.cxx.

References Create3DGeoHits(), PndMvdEventAna::fActiveEvent, fClusterArray, fDigiArray, PndMvdEventAna::fHitArray, Fill3DHisto(), FillClusterHistos(), FillDigiHistos(), FillHitHistos(), FillHitPerClusterHistos(), FillHitProjHistos(), FillHitResolutionHistos(), FillRecoHistos(), fRecoArray, PndMvdEventAna::fTree, PrintClusterArray(), PrintDigiArray(), PrintHitArray(), and PrintRecoArray().

 {
   fHitArray->Delete();
   fDigiArray->Delete();
   fClusterArray->Delete();
   fRecoArray->Delete();
 
   fTree->GetEntry(fActiveEvent);
 
   std::cout << "hitArray: " << fHitArray->GetEntries() << std::endl;
   std::cout << "digiArray: " << fDigiArray->GetEntries() << std::endl;
   std::cout << "clusterArray: " << fClusterArray->GetEntries() << std::endl;
   std::cout << "recoArray: " << fRecoArray->GetEntries() << std::endl;
 
   PrintHitArray();
   PrintDigiArray();
   PrintClusterArray();
   PrintRecoArray();
 
   FillHitHistos();
   FillDigiHistos();
   FillClusterHistos();
   FillRecoHistos();
   FillHitPerClusterHistos();
   FillHitResolutionHistos();
   Fill3DHisto();
   FillHitProjHistos();
 
   Create3DGeoHits();
 }

virtual void PndMvdAllDataEventAna::BeginOfEventAction ( )

inlinevirtual

Reimplemented from PndMvdEventAna.

Definition at line 49 of file PndMvdAllDataEventAna.h.

References ClearAllVectors().

                                          {
                 ClearAllVectors();
         };

TVector3 PndMvdAllDataEventAna::CalcMeanHitPos ( std::vector< Int_t > points )

Definition at line 426 of file PndMvdAllDataEventAna.cxx.

References Double_t, energy, PndMvdEventAna::fHitArray, PndSdsMCPoint::GetPosition(), PndSdsMCPoint::GetPositionOut(), and i.

Referenced by FillHitResolutionHistos().

 {
   TVector3 result;
   Double_t energy = 0;
   for (UInt_t i = 0; i < points.size(); i++){
     PndSdsMCPoint *myPoint = (PndSdsMCPoint*)fHitArray->At(points[i]);
     TVector3 posIn = myPoint->GetPosition();
     std::cout << "posIn: " << posIn.X() << " " << posIn.Y() << " " << posIn.Z() << std::endl;
     TVector3 posOut = myPoint->GetPositionOut();
     std::cout << "posOut: " << posOut.X() << " " << posOut.Y() << " " << posOut.Z() << std::endl;
 
     TVector3 posHit = posIn + posOut;
     posHit *= 0.5;
 
     std::cout << "posHit: " << posHit.X() << " " << posHit.Y() << " " << posHit.Z() << std::endl;
     //posHit *= myPoint->GetEnergyLoss();
     energy += myPoint->GetEnergyLoss();
     result += posHit;
   }
   //Double_t div = 1/energy;
   Double_t div = 1/(Double_t)points.size();
   result *= div;
   return result;
 }

PndMvdEventAna::ClassDef	(	PndMvdEventAna	,
		1
	)

protectedinherited

PndMvdAllDataEventAna::ClassDef	(	PndMvdAllDataEventAna	,
		1
	)

private

void PndMvdAllDataEventAna::ClearAllHMaps ( )

Definition at line 691 of file PndMvdAllDataEventAna.cxx.

References ClearHistoMaps(), fClusterHistos, fDigiHistos, PndMvdEventAna::fHistos, fHitPerClusterHistos, fHitResolutionHistos, and fRecoHistos.

Referenced by ClearHistos().

 {
   ClearHistoMaps(&fHistos);
   ClearHistoMaps(&fRecoHistos);
   ClearHistoMaps(&fDigiHistos);
   ClearHistoMaps(&fClusterHistos);
   ClearHistoMaps(&fHitPerClusterHistos);
   ClearHistoMaps(&fHitResolutionHistos);
 }

void PndMvdAllDataEventAna::ClearAllVectors ( )

Definition at line 701 of file PndMvdAllDataEventAna.cxx.

References ClearHistoVector(), fClusterHistoVec, fDigiHistoVec, fHitHistoVec, and fRecoHistoVec.

Referenced by BeginOfEventAction(), and ClearHistos().

 {
   ClearHistoVector(&fHitHistoVec);
   ClearHistoVector(&fDigiHistoVec);
   ClearHistoVector(&fClusterHistoVec);
   ClearHistoVector(&fRecoHistoVec);
 }

void PndMvdAllDataEventAna::ClearHistoMaps ( std::map< TString, TH1 * > * myHistos ) const

Definition at line 681 of file PndMvdAllDataEventAna.cxx.

Referenced by ClearAllHMaps().

 {
   for (std::map<TString,TH1*>::const_iterator ki = myMaps->begin(); ki != myMaps->end(); ki++){
     ki->second->Delete();
     (*myMaps)[ki->first] = 0;
   }
   myMaps->clear();
   std::cout << " Maps.size: " << myMaps->size() << std::endl;
 }

void PndMvdAllDataEventAna::ClearHistos ( )

virtual

Reimplemented from PndMvdEventAna.

Definition at line 132 of file PndMvdAllDataEventAna.cxx.

References ClearAllHMaps(), ClearAllVectors(), PndEventDisplay::ClearHits(), and fGeoList.

 {
   ClearAllHMaps();
   ClearAllVectors();
 
   fGeoList->ClearHits("RecoHits");
   fGeoList->ClearHits("MCHits");
   fGeoList->ClearHits("StripHits");
   fGeoList->ClearHits("PixHits4");
   fGeoList->ClearHits("PixHits5");
   fGeoList->ClearHits("PixHits8");
   fGeoList->ClearHits("StripHitsGreen");
   fGeoList->ClearHits("PixHits4Green");
   fGeoList->ClearHits("PixHits5Green");
   fGeoList->ClearHits("PixHits8Green");
 }

void PndMvdAllDataEventAna::ClearHistoVector ( std::vector< TH1 * > * myVectors ) const

Definition at line 708 of file PndMvdAllDataEventAna.cxx.

References i.

Referenced by ClearAllVectors(), and FillHitProjHistos().

 {
   for (UInt_t i = 0; i < myVectors->size(); i++)
   {
     delete ((*myVectors)[i]);
   }
   myVectors->clear();
 }

void PndMvdAllDataEventAna::Create3DGeoHits ( )

Definition at line 652 of file PndMvdAllDataEventAna.cxx.

References PndEventDisplay::AddHit(), PndMvdEventAna::fGeoH, fGeoList, PndMvdEventAna::fHitArray, fRecoArray, PndGeoHandling::GetPath(), gGeoManager, i, and PndEventDisplay::SetHits().

Referenced by AnaHits(), and EndOfEventAction().

 {
   fGeoList->SetHits("RecoHits",fRecoArray);
 
   for (Int_t i = 0; i < fHitArray->GetEntriesFast(); i++){
     PndSdsMCPoint* myPoint = (PndSdsMCPoint*)fHitArray->At(i);
     fGeoList->AddHit("MCHits", myPoint->GetX(), myPoint->GetY(), myPoint->GetZ());
     std::cout << "gGeoManager cd: " << myPoint->GetDetName() << " " << gGeoManager->cd(fGeoH->GetPath(myPoint->GetDetName())) << std::endl;
     TGeoHMatrix* mat = gGeoManager->GetCurrentMatrix();
     mat->Print();
     if (myPoint->GetDetName().Contains("71")){
       std::cout << "StripHits added!" << std::endl;
       fGeoList->AddHit("StripHits", mat);
       fGeoList->AddHit("StripHitsGreen", mat, kFALSE);
     }
     else if (myPoint->GetDetName().Contains("89")){
       std::cout << "PixHits4 added!" << std::endl;
       fGeoList->AddHit("PixHits4", mat);
       fGeoList->AddHit("PixHits4Green", mat, kFALSE);
     }
     else if (myPoint->GetDetName().Contains("56")){
       std::cout << "PixHits8 added!" << std::endl;
       fGeoList->AddHit("PixHits8", mat);
       fGeoList->AddHit("PixHits8Green", mat, kFALSE);
     }
   }
 }

void PndMvdAllDataEventAna::Draw3D	(	TString	opt = `""`,
		TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 517 of file PndMvdAllDataEventAna.cxx.

References f3DMCHisto, and f3DRecoHisto.

 {
   if (extCan == 0)
       extCan = new TCanvas();
   extCan->cd(pad);
   f3DMCHisto->SetMarkerStyle(8);
   f3DMCHisto->SetMarkerSize(0.5);
   f3DMCHisto->SetMarkerColor(2);
   f3DMCHisto->Draw(opt.Data());
   f3DRecoHisto->SetMarkerStyle(8);
   f3DRecoHisto->SetMarkerSize(0.5);
   f3DRecoHisto->SetMarkerColor(1);
   f3DRecoHisto->Draw("same");
 
 }

void PndMvdAllDataEventAna::DrawAllHistos	(	TString	detName,
		TCanvas *	extCan = `0`
	)

Definition at line 485 of file PndMvdAllDataEventAna.cxx.

References DrawClusterHisto(), DrawDigiHisto(), DrawHitHisto(), DrawRecoHisto(), fCan1, fHitPerClusterHistos, and fHitResolutionHistos.

Referenced by DrawAllHistos().

 {
   if (extCan)
       fCan1 = extCan;
   if (fCan1 == 0)
     fCan1 = new TCanvas();
   DrawHitHisto(detName, fCan1, 1);
   DrawDigiHisto(detName, fCan1, 2);
   DrawClusterHisto(detName, fCan1, 3);
   DrawRecoHisto(detName, fCan1, 4);
   fCan1->cd(5);
   fHitPerClusterHistos[detName]->Draw();
   fCan1->cd(6);
   fHitResolutionHistos[detName]->Draw();
   fCan1->Modified(kTRUE);
   fCan1->Update();
 }

void PndMvdAllDataEventAna::DrawAllHistos	(	Int_t	i,
		TCanvas *	extCan = `0`
	)

Definition at line 503 of file PndMvdAllDataEventAna.cxx.

References DrawAllHistos(), PndMvdEventAna::fHistos, i, and TString.

 {
   std::map<TString,TH1*>::const_iterator ki;
   Int_t histoSize = fHistos.size();
 
   if (index >= histoSize) return;
   ki = fHistos.begin();
   for (Int_t i = 0; i < index; i++)
     ki++;
   TString detName = ki->first;
   std::cout << "DetName: " << detName.Data() << std::endl;
   DrawAllHistos(detName, extCan);
 }

void PndMvdAllDataEventAna::DrawAllTracks	(	TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 533 of file PndMvdAllDataEventAna.cxx.

References fGeoTrackArray, and i.

 {
   if (extCan == 0)
     extCan = new TCanvas();
   extCan->cd(pad);
 
   for (Int_t i = 0; i < fGeoTrackArray->GetEntries(); i++){
     TGeoTrack* myTrack = (TGeoTrack*) fGeoTrackArray->At(i);
     myTrack->Print();
     myTrack->Draw();
   }
 }

void PndMvdAllDataEventAna::DrawClusterHisto	(	TString	detName,
		TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 468 of file PndMvdAllDataEventAna.cxx.

References fClusterHistos.

Referenced by DrawAllHistos().

 {
   if (extCan)
       extCan->cd(pad);
 
   if(fClusterHistos[detName] != 0)
     fClusterHistos[detName]->Draw("colz");
 }

void PndMvdAllDataEventAna::DrawDigiHisto	(	TString	detName,
		TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 460 of file PndMvdAllDataEventAna.cxx.

References fDigiHistos.

Referenced by DrawAllHistos().

 {
   if (extCan)
     extCan->cd(pad);
 
   if(fDigiHistos[detName] != 0)
     fDigiHistos[detName]->Draw("colz");
 }

void PndMvdAllDataEventAna::DrawEvent	(	bool	tracks = `false`,
		TCanvas *	extCan = `0`
	)

Definition at line 604 of file PndMvdAllDataEventAna.cxx.

References DrawHistoVec(), DrawHitTracks(), DrawTopVolume(), fAllHitResolutionHistos, fCan2, fRecoHisrz, and fRecoHisxy.

 {
   if (extCan)
         fCan2 = extCan;
   if (fCan2 == 0)
     fCan2 = new TCanvas();
 
   fCan2->Divide(2,2);
 
   fCan2->cd(1);
   DrawTopVolume();
 
   if (tracks)
     DrawHitTracks();
 
   fCan2->cd(2);
   fAllHitResolutionHistos->Draw();
 
   fCan2->cd(3);
   DrawHistoVec(&fRecoHisxy);
 
   fCan2->cd(4);
   DrawHistoVec(&fRecoHisrz);
 }

void PndMvdAllDataEventAna::DrawHisrz	(	TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 577 of file PndMvdAllDataEventAna.cxx.

References DrawHistoVec(), and fRecoHisxy.

 {
   DrawHistoVec(&fRecoHisxy, extCan, pad);
 }

void PndMvdEventAna::DrawHisto	(	TString	path,
		TString	opt = `""`
	)

inherited

Definition at line 74 of file PndMvdEventAna.cxx.

References PndMvdEventAna::fHistos.

 {
         fHistos[path]->Draw(opt.Data());
 }

void PndMvdAllDataEventAna::DrawHisto	(	TH1 *	histo,
		TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)		const

Definition at line 597 of file PndMvdAllDataEventAna.cxx.

Referenced by DrawHitPerCluster(), and DrawResolution().

 {
   if (extCan)
     extCan->cd(pad);
   histo->Draw();
 }

void PndMvdEventAna::DrawHistos ( )

inherited

Definition at line 63 of file PndMvdEventAna.cxx.

References PndMvdEventAna::fCan, PndMvdEventAna::fHistos, and i.

 {
         Int_t i = 0;
         for (std::map<TString, TH1*>::const_iterator ki = fHistos.begin(); ki != fHistos.end(); ki++){
                 i++;
                 fCan->cd(i);
                 std::cout << ki->second << std::endl;
                 ki->second->Draw();//(fDrawOption[ci->second->GetName()]);
         }
 }

void PndMvdAllDataEventAna::DrawHistoVec	(	std::vector< TH1 * > *	vec,
		TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)		const

Definition at line 587 of file PndMvdAllDataEventAna.cxx.

References i.

Referenced by DrawEvent(), DrawHisrz(), and DrawHisxy().

 {
   if (extCan)
     extCan->cd(pad);
   if (vec->size() > 0)
       (*vec)[0]->Draw();
     for (UInt_t i = 1; i < vec->size(); i++)
       (*vec)[i]->Draw("same");
 }

void PndMvdAllDataEventAna::DrawHisxy	(	TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 582 of file PndMvdAllDataEventAna.cxx.

References DrawHistoVec(), and fRecoHisrz.

 {
   DrawHistoVec(&fRecoHisrz, extCan, pad);
 }

void PndMvdAllDataEventAna::DrawHitHisto	(	TString	detName,
		TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 451 of file PndMvdAllDataEventAna.cxx.

References PndMvdEventAna::fHistos.

Referenced by DrawAllHistos().

 {
   if (extCan)
     extCan->cd(pad);
 
   if (fHistos[detName] != 0)
     fHistos[detName]->Draw();
 }

void PndMvdAllDataEventAna::DrawHitPerCluster	(	TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

inline

Definition at line 84 of file PndMvdAllDataEventAna.h.

References DrawHisto(), and fAllHitPerClusterHistos.

84 {DrawHisto(fAllHitPerClusterHistos, extCan, pad);};

PndMvdAllDataEventAna::fAllHitPerClusterHistos

TH1 * fAllHitPerClusterHistos

Definition: PndMvdAllDataEventAna.h:125

PndMvdAllDataEventAna::DrawHisto

void DrawHisto(TH1 *histo, TCanvas *extCan=0, Int_t pad=0) const

Definition: PndMvdAllDataEventAna.cxx:597

void PndMvdAllDataEventAna::DrawHitTracks	(	TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 546 of file PndMvdAllDataEventAna.cxx.

References fGeoTrackArray, PndMvdEventAna::fHitArray, and i.

Referenced by DrawEvent().

 {
   if (extCan == 0)
     extCan = new TCanvas();
   extCan->cd(pad);
 
   std::map<Int_t, Int_t> trackUsed;
   for (Int_t i = 0; i < fHitArray->GetEntries(); i++){
     PndSdsMCPoint* myPoint = (PndSdsMCPoint*)fHitArray->At(i);
     std::cout << "Hit " << i << " TrackID " << myPoint->GetTrackID() << std::endl;
     if (trackUsed[myPoint->GetTrackID()] == 0){
       for (Int_t j = 0; j < fGeoTrackArray->GetEntries(); j++){
         TGeoTrack* myTrack = (TGeoTrack*)fGeoTrackArray->At(j);
         if (myPoint->GetTrackID() == myTrack->GetId()){
           myTrack->Draw();
           myTrack->Print();
         }
       }
       trackUsed[myPoint->GetTrackID()] = 1;
     }
   }
 }

void PndMvdAllDataEventAna::DrawRecoHisto	(	TString	detName,
		TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

Definition at line 476 of file PndMvdAllDataEventAna.cxx.

References fRecoHistos.

Referenced by DrawAllHistos().

 {
   if (extCan)
       extCan->cd(pad);
 
   if (fRecoHistos[detName] != 0)
     fRecoHistos[detName]->Draw();
 }

void PndMvdAllDataEventAna::DrawResolution	(	TCanvas *	extCan = `0`,
		Int_t	pad = `0`
	)

inline

Definition at line 85 of file PndMvdAllDataEventAna.h.

References DrawHisto(), and fAllHitResolutionHistos.

85 {DrawHisto(fAllHitResolutionHistos, extCan, pad);};

PndMvdAllDataEventAna::fAllHitResolutionHistos

TH1 * fAllHitResolutionHistos

Definition: PndMvdAllDataEventAna.h:126

PndMvdAllDataEventAna::DrawHisto

void DrawHisto(TH1 *histo, TCanvas *extCan=0, Int_t pad=0) const

Definition: PndMvdAllDataEventAna.cxx:597

void PndMvdAllDataEventAna::DrawTopVolume	(	TCanvas *	extCan = `0`,
		Int_t	pad = `0`,
		const char *	opt = `""`
	)

Definition at line 569 of file PndMvdAllDataEventAna.cxx.

References fMvdTopVolume.

Referenced by DrawEvent().

 {
   if (extCan)
     extCan->cd(pad);
 //  fMvdTopVolume->Draw(""); //ogl for OpenGL Viewer
   fMvdTopVolume->Draw(opt);
 }

void PndMvdEventAna::DrawVector	(	Int_t	index,
		TString	opt
	)

inherited

Definition at line 110 of file PndMvdEventAna.cxx.

References PndMvdEventAna::fHistoVector.

 {
         fHistoVector[index]->Draw(opt.Data());
 }

virtual void PndMvdAllDataEventAna::EndOfEventAction ( )

inlinevirtual

Reimplemented from PndMvdEventAna.

Definition at line 52 of file PndMvdAllDataEventAna.h.

References Create3DGeoHits().

52 {Create3DGeoHits();};

PndMvdAllDataEventAna::Create3DGeoHits

void Create3DGeoHits()

Definition: PndMvdAllDataEventAna.cxx:652

void PndMvdAllDataEventAna::Fill3DHisto ( )

Definition at line 383 of file PndMvdAllDataEventAna.cxx.

References f3DMCHisto, f3DRecoHisto, PndMvdEventAna::fHitArray, fRecoArray, and i.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fRecoArray->GetEntries(); i++){
     PndSdsHit* myHit = (PndSdsHit*)fRecoArray->At(i);
     f3DRecoHisto->Fill(myHit->GetX(), myHit->GetY(), myHit->GetZ());
   }
   for (Int_t j = 0; j < fHitArray->GetEntries(); j++){
       PndSdsMCPoint* myPoint = (PndSdsMCPoint*)fHitArray->At(j);
       f3DMCHisto->Fill(myPoint->GetX(), myPoint->GetY(), myPoint->GetZ());
   }
 
 }

void PndMvdAllDataEventAna::FillClusterHistos ( )

Definition at line 292 of file PndMvdAllDataEventAna.cxx.

References fClusterArray, fClusterHistos, fDigiArray, fPixelCon, PndSdsDigi::GetCharge(), PndSdsCluster::GetClusterList(), PndSdsDigi::GetFE(), PndSdsDigiPixel::GetPixelColumn(), PndSdsDigiPixel::GetPixelRow(), PndSdsDigi::GetSensorID(), hit(), and i.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fClusterArray->GetEntries(); i++){
     PndSdsCluster *cluster = (PndSdsCluster*)fClusterArray->At(i);
   //if (hit->GetDetName().Contains("119_2")){
     std::vector<Int_t> clusterList = cluster->GetClusterList();
     PndSdsDigiPixel* hit = (PndSdsDigiPixel*)fDigiArray->At(clusterList[0]);
     Int_t sensorID = hit->GetSensorID();
     if (fClusterHistos[sensorID] == 0){
       fClusterHistos[sensorID] = new TH2I("clusterHisto",detName.Data(), 1000,0,1000,201,0,200);
       //fDrawOption[hit->GetDetName()] = "colz";
     }
     TH2* tempHisto = (TH2*)(fClusterHistos[detName]);
     for (Int_t j = 0; j < clusterList.size(); j++){
       hit = (PndSdsDigiPixel*)fDigiArray->At(clusterList[j]);
       PndSdsPixel myFePixel(sensorID, hit->GetFE(), hit->GetPixelColumn(), hit->GetPixelRow(), hit->GetCharge());
       PndSdsPixel mySensorPixel = fPixelCon->CalcSensorHit(myFePixel);
       tempHisto->Fill(mySensorPixel.GetCol(), mySensorPixel.GetRow());
     }
   }
 }

void PndMvdAllDataEventAna::FillDigiHistos ( )

Definition at line 277 of file PndMvdAllDataEventAna.cxx.

References fDigiArray, fDigiHistos, fPixelCon, PndSdsDigi::GetCharge(), PndSdsDigi::GetFE(), PndSdsDigiPixel::GetPixelColumn(), PndSdsDigiPixel::GetPixelRow(), PndSdsDigi::GetSensorID(), hit(), and i.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fDigiArray->GetEntries(); i++){
     PndSdsDigiPixel *hit = (PndSdsDigiPixel*)fDigiArray->At(i);
   //if (hit->GetDetName().Contains("119_2")){
     if (fDigiHistos[hit->GetDetName()] == 0){
       fDigiHistos[hit->GetDetName()] = new TH2I("digiHisto",hit->GetDetName().Data(), 1000,0,1000,201,0,200);
       //fDrawOption[hit->GetDetName()] = "colz";
     }
     TH2* tempHisto = (TH2*)(fDigiHistos[hit->GetDetName()]);
     PndSdsPixel myFePixel(hit->GetSensorID(), hit->GetFE(), hit->GetPixelColumn(), hit->GetPixelRow(), hit->GetCharge());
     PndSdsPixel mySensorPixel = fPixelCon->CalcSensorHit(myFePixel);
     tempHisto->Fill(mySensorPixel.GetCol(), mySensorPixel.GetRow(),mySensorPixel.GetCharge());//, (Double_t)(hit->GetCharge()));
   }
 }

void PndMvdAllDataEventAna::FillHitHistos ( )

Definition at line 247 of file PndMvdAllDataEventAna.cxx.

References PndMvdEventAna::fGeoH, PndMvdEventAna::fHistos, PndMvdEventAna::fHitArray, GetLocalHitPoints(), PndGeoHandling::GetPath(), PndSdsMCPoint::GetXOut(), PndSdsMCPoint::GetYOut(), PndSdsMCPoint::GetZOut(), gGeoManager, i, out, and TString.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fHitArray->GetEntries(); i++){
     PndSdsMCPoint *myPoint = (PndSdsMCPoint*)fHitArray->At(i);
     TString detName = myPoint->GetDetName();
   //if (hit->GetDetName().Contains("119_2")){
     if (fHistos[detName] == 0){
       fHistos[detName] = new TH2I("MCHisto",detName.Data(), 1000,0,1000,201,0,200);
       //fDrawOption[hit->GetDetName()] = "colz";
     }
     //TH2* tempHisto = (TH2*)(fHistos[detName]);
     gGeoManager->cd(fGeoH->GetPath(detName.Data()));
      // TGeoHMatrix* transMat = gGeoManager->GetCurrentMatrix();
 
     TVector3 in(myPoint->GetX(), myPoint->GetY(), myPoint->GetZ());
     TVector3 out(myPoint->GetXOut(), myPoint->GetYOut(), myPoint->GetZOut());
     TVector3 inLocal, outLocal;
     inLocal = GetLocalHitPoints(detName, in);
     outLocal = GetLocalHitPoints(detName, out);
 
     fHistos[detName]->Fill(inLocal.X()*100, inLocal.Y()*100);
     fHistos[detName]->Fill(outLocal.X()*100, outLocal.Y()*100);
 
 //    std::cout << i << ": " << std::endl;
 //     for (Int_t i = 0; i < 3; i++){
 //          std::cout << "posInLocal "<< i << ": " << inLocal(i) << std::endl;
 //            std::cout << "posOutLocal "<< i << ": " << outLocal(i) << std::endl;
 //        }
   }
 }

void PndMvdAllDataEventAna::FillHitPerClusterHistos ( )

Definition at line 339 of file PndMvdAllDataEventAna.cxx.

References fAllHitPerClusterHistos, fClusterArray, PndMvdEventAna::fHitArray, fHitPerClusterHistos, GetHitPerCluster(), hits, i, point, and TString.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fClusterArray->GetEntries(); i++){
     PndSdsCluster* cand = (PndSdsCluster*)fClusterArray->At(i);
     std::vector<Int_t> hits = GetHitPerCluster(cand);
     PndSdsMCPoint* point = (PndSdsMCPoint*)(fHitArray->At(hits[0]));
     TString detName = point->GetDetName();
     if (fHitPerClusterHistos[detName] == 0){
       fHitPerClusterHistos[detName] = new TH1I("HitPerCluster",detName.Data(), 21,0,20);
       //fDrawOption[hit->GetDetName()] = "colz";
     }
     fHitPerClusterHistos[detName]->Fill(hits.size());
     fAllHitPerClusterHistos->Fill(hits.size());
     std::cout << "Hits in Cluster: " << hits.size() << " for Det: " << detName.Data() << std::endl;
   }
 }

void PndMvdAllDataEventAna::FillHitProjHistos ( )

Definition at line 396 of file PndMvdAllDataEventAna.cxx.

References ClearHistoVector(), fRecoArray, fRecoHisrz, fRecoHisxy, fTrackFArray, GFTrackCand::getHit(), GFTrackCand::getNHits(), i, and vec.

Referenced by AnaHits().

 {
   ClearHistoVector(&fRecoHisxy);
   ClearHistoVector(&fRecoHisrz);
   unsigned int detID, hitID;
   TVector3 vec;
 
   for (Int_t i = 0; i < fTrackFArray->GetEntries(); i++){
     GFTrackCand* trackC = (GFTrackCand*)fTrackFArray->At(i);
     TH2D* myHistoXY = new TH2D("recohisxy","MVD Reco Points, xy view",400,-15.,15.,400,-15.,15.);
     myHistoXY->SetMarkerStyle(7);
     myHistoXY->SetMarkerColor(i+1);
 
     TH2D* myHistoRZ = new TH2D("recohisrz","MVD Reco Points, rz view",400,-20.,20.,400,-20.,20.);
     myHistoRZ->SetMarkerStyle(7);
     myHistoRZ->SetMarkerColor(i+1);
 
     fRecoHisxy.push_back(myHistoXY);
     fRecoHisrz.push_back(myHistoRZ);
     for (UInt_t j = 0; j < trackC->getNHits(); j++){
       trackC->getHit(j, detID, hitID);
       PndSdsHit* myHit = (PndSdsHit*)fRecoArray->At(hitID);
       vec.SetXYZ(myHit->GetX(), myHit->GetY(), myHit->GetZ());
       myHistoXY->Fill(vec.x(), vec.y());
       myHistoRZ->Fill(vec.z(), vec.Perp());
     }
   }
 }

void PndMvdAllDataEventAna::FillHitResolutionHistos ( )

Definition at line 356 of file PndMvdAllDataEventAna.cxx.

References CalcMeanHitPos(), fAllHitResolutionHistos, fClusterArray, fHitResolutionHistos, fRecoArray, PndSdsHit::GetClusterIndex(), GetHitPerCluster(), PndSdsHit::GetPosition(), i, and TString.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fRecoArray->GetEntries(); i++){
     PndSdsHit *myHit = (PndSdsHit*)fRecoArray->At(i);
     TString detName = myHit->GetDetName();
     std::cout << "HitResolution for: " << detName << std::endl;
     TVector3 recoPos = myHit->GetPosition();
     std::cout << "RecoPos: " << recoPos.X() << " " << recoPos.Y() << " " << recoPos.Z() << std::endl;
     std::cout << "MC point RefIndex: " << myHit->GetRefIndex() << std::endl;
     std::cout << "Cluster Index: " << myHit->GetClusterIndex() << std::endl;
     PndSdsCluster *myCand = (PndSdsCluster*)(fClusterArray->At(myHit->GetClusterIndex()));
     std::vector<Int_t> points = GetHitPerCluster(myCand);
     for (UInt_t k = 0; k < points.size(); i++)
       std::cout << "ClusterPoints: " << points[k] << std::endl;
     TVector3 hitPos = CalcMeanHitPos(points);
     std::cout << "HitPos: " << hitPos.X() << " " << hitPos.Y() << " " << hitPos.Z() << std::endl;
     if (fHitResolutionHistos[detName] == 0){
       fHitResolutionHistos[detName] = new TH1D("HitRes",detName.Data(), 1000,0,0.01);
       //fDrawOption[hit->GetDetName()] = "colz";
     }
     TVector3 resPos = recoPos - hitPos;
     fHitResolutionHistos[detName]->Fill(resPos.Mag());
     fAllHitResolutionHistos->Fill(resPos.Mag());
     std::cout << "HitRes: " << resPos.Mag() << std::endl;
   }
 }

void PndMvdAllDataEventAna::FillRecoHistos ( )

Definition at line 314 of file PndMvdAllDataEventAna.cxx.

References fRecoArray, fRecoHistos, GetLocalHitPoints(), i, and TString.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fRecoArray->GetEntries(); i++){
     PndSdsHit *myHit = (PndSdsHit*)fRecoArray->At(i);
     TString detName = myHit->GetDetName();
   //if (hit->GetDetName().Contains("119_2")){
     if (fRecoHistos[detName] == 0){
       fRecoHistos[detName] = new TH2I("RecoHisto",detName.Data(), 1000,0,1000,201,0,200);
       //fDrawOption[hit->GetDetName()] = "colz";
     }
     //TH2* tempHisto = (TH2*)(fRecoHistos[detName]);
 
     TVector3 in(myHit->GetX(), myHit->GetY(), myHit->GetZ());
     TVector3 inLocal;
     inLocal = GetLocalHitPoints(detName, in);
 
     fRecoHistos[detName]->Fill(inLocal.X()*100, inLocal.Y()*100);
 
 //    std::cout << i << ": " << std::endl;
 //     for (Int_t i = 0; i < 3; i++){
 //          std::cout << "posInLocal "<< i << ": " << inLocal(i) << std::endl;
 //        }
   }
 }

void PndMvdEventAna::FillVector ( )

inherited

Definition at line 102 of file PndMvdEventAna.cxx.

References PndMvdEventAna::fHistos, and PndMvdEventAna::fHistoVector.

 {
         fHistoVector.clear();
         for (std::map<TString, TH1*>::const_iterator ki = fHistos.begin(); ki != fHistos.end(); ki++){
                 fHistoVector.push_back(ki->second);
         }
 }

Int_t PndMvdEventAna::GetEventNr ( ) const

inlineinherited

Definition at line 52 of file PndMvdEventAna.h.

References PndMvdEventAna::fActiveEvent.

52 {return fActiveEvent;};

PndMvdEventAna::fActiveEvent

Int_t fActiveEvent

Definition: PndMvdEventAna.h:62

TH1* PndMvdEventAna::GetHisto ( TString histoName )

inlineinherited

Definition at line 54 of file PndMvdEventAna.h.

References PndMvdEventAna::fHistos.

54 {return fHistos[histoName];};

PndMvdEventAna::fHistos

std::map< TString, TH1 * > fHistos

Definition: PndMvdEventAna.h:74

TH1* PndMvdEventAna::GetHisto ( Int_t i )

inlineinherited

Definition at line 55 of file PndMvdEventAna.h.

References PndMvdEventAna::fHistoVector, and i.

55 {return fHistoVector[i];};

i

Int_t i

Definition: run_full.C:25

PndMvdEventAna::fHistoVector

std::vector< TH1 * > fHistoVector

Definition: PndMvdEventAna.h:76

std::vector< Int_t > PndMvdAllDataEventAna::GetHitPerCluster ( PndSdsCluster * clusterCand )

Definition at line 629 of file PndMvdAllDataEventAna.cxx.

References fDigiArray, PndSdsCluster::GetClusterList(), PndSdsDigi::GetIndex(), and i.

Referenced by FillHitPerClusterHistos(), and FillHitResolutionHistos().

 {
   std::vector<Int_t> result;
   std::vector<Int_t> digiPos = clusterCand->GetClusterList();
   bool isInResult = false;
   for (UInt_t i = 0; i < digiPos.size(); i++){
     PndSdsDigiPixel* digiHit = (PndSdsDigiPixel*)(fDigiArray->At(digiPos[i]));
     Int_t mcID = digiHit->GetIndex(0);
     std::cout << " -I- GetHitPerCluster: mcID: " << mcID << std::endl;
       for (UInt_t j = 0; j < result.size() && isInResult == false; j++){
         //std::cout << "Result: " << result[j] << std::endl;
         if (mcID == result[j])
           isInResult = true;
       }
     if (isInResult == false)
       result.push_back(mcID);
     isInResult = false;
   }
   for (UInt_t k = 0; k < result.size(); k++)
     std::cout << " Result: " << k << ": " << result[k] << std::endl;
   return result;
 }

TVector3 PndMvdAllDataEventAna::GetLocalHitPoints	(	TString	detName,
		TVector3	input
	)

Definition at line 222 of file PndMvdAllDataEventAna.cxx.

References Double_t, PndMvdEventAna::fGeoH, PndGeoHandling::GetPath(), and gGeoManager.

Referenced by FillHitHistos(), and FillRecoHistos().

 {
   TVector3 result;
   Double_t in[3];
   Double_t local[3];
 
   in[0] = input.X();
   in[1] = input.Y();
   in[2] = input.Z();
 
   gGeoManager->cd(fGeoH->GetPath(detName.Data()));
     TGeoHMatrix* transMat = gGeoManager->GetCurrentMatrix();
 
   transMat->MasterToLocal(in, local);
 
   TGeoVolume* actVolume = gGeoManager->GetCurrentVolume();
   TGeoBBox* actBox = (TGeoBBox*)(actVolume->GetShape());
 
   result.SetX(local[0] + actBox->GetDX());
   result.SetY(local[1] + actBox->GetDY());
   result.SetZ(local[2] + actBox->GetDZ());
 
   return result;
 }

std::vector< TString > PndMvdAllDataEventAna::GetModulesHit ( )

Definition at line 717 of file PndMvdAllDataEventAna.cxx.

References PndMvdEventAna::fGeoH, fRecoHistos, and PndGeoHandling::GetPath().

 {
     std::vector<TString> result;
     for(std::map<TString, TH1*>::const_iterator ki= fRecoHistos.begin(); ki != fRecoHistos.end(); ++ki){
       std::cout << "First: " << ki->first << std::endl;
       std::cout << "Path: " << fGeoH->GetPath(ki->first) << std::endl;
       result.push_back(fGeoH->GetPath(ki->first));
     }
     return result;
 }

TString PndMvdEventAna::GetSelection ( ) const

inlineinherited

Definition at line 53 of file PndMvdEventAna.h.

References PndMvdEventAna::fSelection.

53 {return fSelection;};

PndMvdEventAna::fSelection

TString fSelection

Definition: PndMvdEventAna.h:72

void PndMvdAllDataEventAna::Init ( TString filename )

virtual

Reimplemented from PndMvdEventAna.

Definition at line 92 of file PndMvdAllDataEventAna.cxx.

References f3DMCHisto, f3DRecoHisto, fAllHitPerClusterHistos, fAllHitResolutionHistos, fClusterArray, fDigiArray, PndMvdEventAna::fFile, fGeoTrackArray, PndMvdEventAna::fHitArray, fRecoArray, fTrackFArray, and PndMvdEventAna::fTree.

Referenced by PndMvdAllDataEventAna().

 {
   PndFileNameCreator nameCreator(fileName.Data());
   fFile = new TFile(fileName.Data());
   fTree = (TTree*)fFile->Get("pndsim");
 
   fTree->AddFriend("digi=pndsim",nameCreator.GetDigiFileName().c_str());
   fTree->AddFriend("reco=pndsim",nameCreator.GetRecoFileName().c_str());
   fTree->AddFriend("trackF=pndsim",nameCreator.GetTrackFindingFileName().c_str());
 
   fHitArray=new TClonesArray("PndSdsMCPoint");
   fDigiArray = new TClonesArray("PndSdsDigiPixel");
   fClusterArray = new TClonesArray("PndSdsCluster");
   fRecoArray = new TClonesArray("PndSdsHit");
   fGeoTrackArray = new TClonesArray("TGeoTrack");
   fTrackFArray = new TClonesArray("GFTrackCand");
 
   fTree->SetBranchAddress("MVDPoint",&fHitArray);
   fTree->SetBranchAddress("MVDPixelDigis",&fDigiArray);
   fTree->SetBranchAddress("MVDClusterCand",&fClusterArray);
   fTree->SetBranchAddress("MVDClusterHit",&fRecoArray);
   fTree->SetBranchAddress("MVDTrackCand",&fTrackFArray);
   if (fTree->FindBranch("GeoTracks") != 0)
     fTree->SetBranchAddress("GeoTracks", &fGeoTrackArray);
 
   fAllHitPerClusterHistos = new TH1I("hitPerCluster","hitPerCluster",21,0,20);
   fAllHitResolutionHistos = new TH1D("AllHitResolution","AllHitResolution", 1000,0,0.1);
   f3DMCHisto = new TH3D("h3D","h3D", 50,-15,15, 50, -15,15, 50, -30,30);
   f3DRecoHisto = new TH3D("h3D","h3D", 50,-15,15, 50, -15,15, 50, -30,30);
 
 }

void PndMvdAllDataEventAna::InitBranch ( )

virtual

Implements PndMvdEventAna.

Definition at line 124 of file PndMvdAllDataEventAna.cxx.

125 {

126 }

void PndMvdEventAna::InitCanvas ( )

inherited

Definition at line 55 of file PndMvdEventAna.cxx.

References Double_t, PndMvdEventAna::fCan, PndMvdEventAna::fCancol, and PndMvdEventAna::fHistos.

 {
         if (fCan != 0) delete(fCan);
         fCan = new TCanvas("can1","MCHit view in MVD",0,0,1000,1000);
         std::cout << "InitCanvas fHistos.size() " << fHistos.size() << std::endl;
         fCan->Divide(fCancol, (Int_t)(((Double_t)(fHistos.size())/fCancol) + 0.5));
 }

void PndMvdAllDataEventAna::InitHistos ( )

virtual

Implements PndMvdEventAna.

Definition at line 128 of file PndMvdAllDataEventAna.cxx.

129 {

130 }

void PndMvdEventAna::NextEvent ( )

inherited

Definition at line 46 of file PndMvdEventAna.cxx.

References PndMvdEventAna::AnaEventNr(), PndMvdEventAna::fActiveEvent, PndMvdEventAna::fAnaEvents, and PndMvdEventAna::fTree.

 {
         fAnaEvents = 1;
         if (fActiveEvent + 1 < fTree->GetEntriesFast()){
                 fActiveEvent++;
                 AnaEventNr(fActiveEvent);
         }
 }

void PndMvdAllDataEventAna::PrintClusterArray ( )

Definition at line 197 of file PndMvdAllDataEventAna.cxx.

References fClusterArray, fDigiArray, PndSdsCluster::GetClusterList(), and i.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fClusterArray->GetEntries(); i++){
     std::cout << "Cluster " << i << ": " << std::endl;
     PndSdsCluster* myCluster = (PndSdsCluster*)fClusterArray->At(i);
     //myCluster->Print();
     std::vector<Int_t> myClusterList = myCluster->GetClusterList();
     for (UInt_t j = 0; j < myClusterList.size(); j++){
       PndSdsDigiPixel* myDigi = (PndSdsDigiPixel*)fDigiArray->At(myClusterList[j]);
       std::cout << myClusterList[j] << ": ";
       std::cout << *myDigi;
     }
   }
 
 }

void PndMvdAllDataEventAna::PrintDigiArray ( )

Definition at line 188 of file PndMvdAllDataEventAna.cxx.

References fDigiArray, and i.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fDigiArray->GetEntries(); i++){
     PndSdsDigiPixel* myDigi = (PndSdsDigiPixel*)fDigiArray->At(i);
     std::cout << i << ": ";
     //myDigi->Print();
     std::cout << *myDigi;
   }
 }

void PndMvdAllDataEventAna::PrintHitArray ( )

Definition at line 180 of file PndMvdAllDataEventAna.cxx.

References PndMvdEventAna::fHitArray, and i.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fHitArray->GetEntries(); i++){
     std::cout << i << ": ";
     PndSdsMCPoint* myPoint = (PndSdsMCPoint*)fHitArray->At(i);
     std::cout << *myPoint;
   }
 }

void PndMvdAllDataEventAna::PrintRecoArray ( )

Definition at line 212 of file PndMvdAllDataEventAna.cxx.

References fRecoArray, and i.

Referenced by AnaHits().

 {
   for (Int_t i = 0; i < fRecoArray->GetEntries(); i++){
     std::cout << i << ": ";
     PndSdsHit* myHit = (PndSdsHit*)fRecoArray->At(i);
     std::cout << *myHit;
   }
 }

void PndMvdEventAna::SetCanvasColumns ( Int_t col )

inlineinherited

Definition at line 47 of file PndMvdEventAna.h.

References col, and PndMvdEventAna::fCancol.

Referenced by PndMvdAllDataEventAna(), and PndMvdDigiEventAna::PndMvdDigiEventAna().

47 {fCancol = col;};

col

int col

Definition: anaLmdDigi.C:67

PndMvdEventAna::fCancol

Int_t fCancol

Definition: PndMvdEventAna.h:77

void PndMvdEventAna::SetDrawOption	(	TString	histoName,
		TString	option
	)

inherited

Definition at line 115 of file PndMvdEventAna.cxx.

References PndMvdEventAna::fDrawOption.

 {
         fDrawOption[histoName] = option;
 }

void PndMvdEventAna::SetSelection ( TString sel )

inlineinherited

Definition at line 49 of file PndMvdEventAna.h.

References PndMvdEventAna::fSelection.

49 {fSelection = sel;};

PndMvdEventAna::fSelection

TString fSelection

Definition: PndMvdEventAna.h:72

Member Data Documentation

TH3D* PndMvdAllDataEventAna::f3DMCHisto

private

Definition at line 127 of file PndMvdAllDataEventAna.h.

Referenced by Draw3D(), Fill3DHisto(), and Init().

TH3D* PndMvdAllDataEventAna::f3DRecoHisto

private

Definition at line 128 of file PndMvdAllDataEventAna.h.

Referenced by Draw3D(), Fill3DHisto(), and Init().

Int_t PndMvdEventAna::fActiveEvent

protectedinherited

Definition at line 62 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::AnaEventNr(), PndMvdEventAna::AnaEvents(), PndMvdDigiEventAna::AnaHitNr(), PndMvdMCEventAna::AnaHitNr(), AnaHits(), PndMvdEventAna::GetEventNr(), PndMvdEventAna::Init(), and PndMvdEventAna::NextEvent().

TH1* PndMvdAllDataEventAna::fAllHitPerClusterHistos

private

Definition at line 125 of file PndMvdAllDataEventAna.h.

Referenced by DrawHitPerCluster(), FillHitPerClusterHistos(), and Init().

TH1* PndMvdAllDataEventAna::fAllHitResolutionHistos

private

Definition at line 126 of file PndMvdAllDataEventAna.h.

Referenced by DrawEvent(), DrawResolution(), FillHitResolutionHistos(), and Init().

Int_t PndMvdEventAna::fAnaEvents

protectedinherited

Definition at line 66 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::AnaEventNr(), PndMvdEventAna::AnaEvents(), PndMvdDigiEventAna::CreateMaxHitsHisto(), and PndMvdEventAna::NextEvent().

TString PndMvdEventAna::fBranchName

protectedinherited

Definition at line 71 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::Init(), PndMvdMCEventAna::InitBranch(), and PndMvdDigiEventAna::InitBranch().

TCanvas* PndMvdEventAna::fCan

protectedinherited

Definition at line 73 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::DrawHistos(), and PndMvdEventAna::InitCanvas().

TCanvas* PndMvdAllDataEventAna::fCan1

private

Definition at line 134 of file PndMvdAllDataEventAna.h.

Referenced by DrawAllHistos().

TCanvas* PndMvdAllDataEventAna::fCan2

private

Definition at line 135 of file PndMvdAllDataEventAna.h.

Referenced by DrawEvent().

Int_t PndMvdEventAna::fCancol

protectedinherited

Definition at line 77 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::InitCanvas(), and PndMvdEventAna::SetCanvasColumns().

TString PndMvdEventAna::fClassName

protectedinherited

Definition at line 70 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::Init(), PndMvdMCEventAna::InitBranch(), and PndMvdDigiEventAna::InitBranch().

TClonesArray* PndMvdAllDataEventAna::fClusterArray

protected

Definition at line 104 of file PndMvdAllDataEventAna.h.

Referenced by AnaHits(), FillClusterHistos(), FillHitPerClusterHistos(), FillHitResolutionHistos(), Init(), and PrintClusterArray().

std::map<TString, TH1*> PndMvdAllDataEventAna::fClusterHistos

private

Definition at line 110 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllHMaps(), DrawClusterHisto(), and FillClusterHistos().

std::vector<TH1*> PndMvdAllDataEventAna::fClusterHistoVec

private

Definition at line 117 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllVectors().

TClonesArray* PndMvdAllDataEventAna::fDigiArray

protected

Definition at line 103 of file PndMvdAllDataEventAna.h.

Referenced by AnaHits(), FillClusterHistos(), FillDigiHistos(), GetHitPerCluster(), Init(), PrintClusterArray(), and PrintDigiArray().

std::map<TString, TH1*> PndMvdAllDataEventAna::fDigiHistos

private

Definition at line 109 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllHMaps(), DrawDigiHisto(), and FillDigiHistos().

std::vector<TH1*> PndMvdAllDataEventAna::fDigiHistoVec

private

Definition at line 116 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllVectors().

std::map<TString, TString> PndMvdEventAna::fDrawOption

protectedinherited

Definition at line 75 of file PndMvdEventAna.h.

Referenced by PndMvdDigiEventAna::AnaHitNr(), PndMvdMCEventAna::InitHistos(), and PndMvdEventAna::SetDrawOption().

TFile* PndMvdEventAna::fFile

protectedinherited

Definition at line 67 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::Init(), and Init().

PndGeoHandling* PndMvdEventAna::fGeoH

protectedinherited

Definition at line 78 of file PndMvdEventAna.h.

Referenced by Create3DGeoHits(), PndMvdDigiEventAna::CreateMaxHitsHisto(), FillHitHistos(), GetLocalHitPoints(), GetModulesHit(), PndMvdAllDataEventAna(), and PndMvdDigiEventAna::PndMvdDigiEventAna().

PndEventDisplay* PndMvdAllDataEventAna::fGeoList

private

Definition at line 138 of file PndMvdAllDataEventAna.h.

Referenced by ClearHistos(), Create3DGeoHits(), and PndMvdAllDataEventAna().

TClonesArray* PndMvdAllDataEventAna::fGeoTrackArray

protected

Definition at line 106 of file PndMvdAllDataEventAna.h.

Referenced by DrawAllTracks(), DrawHitTracks(), and Init().

std::map<TString, TH1*> PndMvdEventAna::fHistos

protectedinherited

Definition at line 74 of file PndMvdEventAna.h.

Referenced by PndMvdDigiEventAna::AnaHitNr(), PndMvdMCEventAna::AnaHitNr(), ClearAllHMaps(), DrawAllHistos(), PndMvdEventAna::DrawHisto(), PndMvdEventAna::DrawHistos(), DrawHitHisto(), FillHitHistos(), PndMvdEventAna::FillVector(), PndMvdEventAna::GetHisto(), PndMvdEventAna::InitCanvas(), PndMvdMCEventAna::InitHistos(), and PndMvdDigiEventAna::InitHistos().

std::vector<TH1*> PndMvdEventAna::fHistoVector

protectedinherited

Definition at line 76 of file PndMvdEventAna.h.

Referenced by PndMvdDigiEventAna::CreateMaxHitsHisto(), PndMvdEventAna::DrawVector(), PndMvdEventAna::FillVector(), and PndMvdEventAna::GetHisto().

TClonesArray* PndMvdEventAna::fHitArray

protectedinherited

Definition at line 69 of file PndMvdEventAna.h.

Referenced by PndMvdDigiEventAna::AnaHitNr(), PndMvdMCEventAna::AnaHitNr(), AnaHits(), PndMvdEventAna::AnaHits(), CalcMeanHitPos(), Create3DGeoHits(), DrawHitTracks(), Fill3DHisto(), FillHitHistos(), FillHitPerClusterHistos(), PndMvdEventAna::Init(), Init(), and PrintHitArray().

std::vector<TH1*> PndMvdAllDataEventAna::fHitHistoVec

private

Definition at line 115 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllVectors().

std::map<TString, TH1*> PndMvdAllDataEventAna::fHitPerClusterHistos

private

Definition at line 112 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllHMaps(), DrawAllHistos(), and FillHitPerClusterHistos().

std::map<TString, TH1*> PndMvdAllDataEventAna::fHitResolutionHistos

private

Definition at line 113 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllHMaps(), DrawAllHistos(), and FillHitResolutionHistos().

TGeoVolume* PndMvdAllDataEventAna::fMvdTopVolume

private

Definition at line 120 of file PndMvdAllDataEventAna.h.

Referenced by DrawTopVolume(), and PndMvdAllDataEventAna().

PndSdsCalcFePixel* PndMvdAllDataEventAna::fPixelCon

private

Definition at line 137 of file PndMvdAllDataEventAna.h.

Referenced by FillClusterHistos(), FillDigiHistos(), and PndMvdAllDataEventAna().

TClonesArray* PndMvdAllDataEventAna::fRecoArray

protected

Definition at line 105 of file PndMvdAllDataEventAna.h.

Referenced by AnaHits(), Create3DGeoHits(), Fill3DHisto(), FillHitProjHistos(), FillHitResolutionHistos(), FillRecoHistos(), Init(), and PrintRecoArray().

TH1* PndMvdAllDataEventAna::fRecoErrorHisto

private

Definition at line 124 of file PndMvdAllDataEventAna.h.

std::vector<TH1*> PndMvdAllDataEventAna::fRecoHisrz

private

Definition at line 132 of file PndMvdAllDataEventAna.h.

Referenced by DrawEvent(), DrawHisxy(), and FillHitProjHistos().

std::map<TString, TH1*> PndMvdAllDataEventAna::fRecoHistos

private

Definition at line 111 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllHMaps(), DrawRecoHisto(), FillRecoHistos(), and GetModulesHit().

std::vector<TH1*> PndMvdAllDataEventAna::fRecoHistoVec

private

Definition at line 118 of file PndMvdAllDataEventAna.h.

Referenced by ClearAllVectors().

std::vector<TH1*> PndMvdAllDataEventAna::fRecoHisxy

private

Definition at line 131 of file PndMvdAllDataEventAna.h.

Referenced by DrawEvent(), DrawHisrz(), and FillHitProjHistos().

TGeoVolume* PndMvdAllDataEventAna::fRecoTopVolume

private

Definition at line 121 of file PndMvdAllDataEventAna.h.

TGeoVolume* PndMvdAllDataEventAna::fRecoVolume

private

Definition at line 122 of file PndMvdAllDataEventAna.h.

Referenced by PndMvdAllDataEventAna().

TString PndMvdEventAna::fSelection

protectedinherited

Definition at line 72 of file PndMvdEventAna.h.

Referenced by PndMvdMCEventAna::AnaHitNr(), PndMvdEventAna::GetSelection(), and PndMvdEventAna::SetSelection().

std::vector<TH1*> PndMvdAllDataEventAna::fSimHisrz

private

Definition at line 130 of file PndMvdAllDataEventAna.h.

std::vector<TH1*> PndMvdAllDataEventAna::fSimHisxy

private

Definition at line 129 of file PndMvdAllDataEventAna.h.

TClonesArray* PndMvdAllDataEventAna::fTrackFArray

protected

Definition at line 107 of file PndMvdAllDataEventAna.h.

Referenced by FillHitProjHistos(), and Init().

TTree* PndMvdEventAna::fTree

protectedinherited

Definition at line 68 of file PndMvdEventAna.h.

Referenced by PndMvdEventAna::AnaEventNr(), PndMvdEventAna::AnaEvents(), PndMvdDigiEventAna::AnaHitNr(), PndMvdMCEventAna::AnaHitNr(), AnaHits(), PndMvdEventAna::Init(), Init(), and PndMvdEventAna::NextEvent().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation