PndMvdRadDamIonizingTask Class Reference

#include <PndMvdRadDamIonizingTask.h>

Inheritance diagram for PndMvdRadDamIonizingTask:

Public Member Functions
	PndMvdRadDamIonizingTask ()
	~PndMvdRadDamIonizingTask ()
	PndMvdRadDamIonizingTask (const PndMvdRadDamIonizingTask &)=delete
PndMvdRadDamIonizingTask &	operator= (const PndMvdRadDamIonizingTask &)=delete
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
virtual void	Exec (Option_t *opt)
virtual void	FinishEvent ()
virtual void	FinishTask ()
void	SetPersistance (Bool_t p=kTRUE)
Bool_t	GetPersistance ()

Private Member Functions
	ClassDef (PndMvdRadDamIonizingTask, 1)

Private Attributes
Bool_t	fPersistance
TClonesArray *	fMCHits
TClonesArray *	fRadDamHits
PndGeoHandling *	fGeoH
std::map< std::string, TProfile2D * >	fMapDetHistos
TH1D *	fRadDamHisto

Detailed Description

Definition at line 20 of file PndMvdRadDamIonizingTask.h.

Constructor & Destructor Documentation

PndMvdRadDamIonizingTask::PndMvdRadDamIonizingTask

(

)

Definition at line 19 of file PndMvdRadDamIonizingTask.cxx.

                                                   :
    fPersistance(kTRUE),
          fMCHits(NULL),
          fRadDamHits(NULL),
          fGeoH(PndGeoHandling::Instance()),
          fMapDetHistos(),
          fRadDamHisto(NULL)
{
}

PndMvdRadDamIonizingTask::~PndMvdRadDamIonizingTask

(

)

Definition at line 29 of file PndMvdRadDamIonizingTask.cxx.

References fMapDetHistos, and fRadDamHisto.

{
        for (std::map<std::string, TProfile2D*>::iterator it = fMapDetHistos.begin(); it != fMapDetHistos.end(); it++)
                delete(it->second);
        fMapDetHistos.clear();
        delete (fRadDamHisto);
}

PndMvdRadDamIonizingTask::PndMvdRadDamIonizingTask ( const PndMvdRadDamIonizingTask &  )

delete

Member Function Documentation

PndMvdRadDamIonizingTask::ClassDef ( PndMvdRadDamIonizingTask  , 1    )

private

void PndMvdRadDamIonizingTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 79 of file PndMvdRadDamIonizingTask.cxx.

References detname, Double_t, fGeoH, fMapDetHistos, fMCHits, fRadDamHisto, PndGeoHandling::GetPath(), PndSdsMCPoint::GetPosition(), PndSdsMCPoint::GetPositionOut(), PndGeoHandling::GetSensorDimensionsShortId(), PndSdsMCPoint::GetSensorID(), i, PndGeoHandling::MasterToLocalShortId(), mom, and TString.

{
  // Reset output array
  
  PndSdsMCPoint* mcPoint;
  

  for (int i = 0; i < fMCHits->GetEntriesFast(); i++){
          mcPoint = (PndSdsMCPoint*)fMCHits->At(i);
          TVector3 length = mcPoint->GetPosition();
          length -= mcPoint->GetPositionOut();
          TVector3 mom;
          mcPoint->Momentum(mom);
          Double_t energyLoss = mcPoint->GetEnergyLoss();
          Double_t dEdx;
          if (length.Mag() > 0)
                  dEdx = energyLoss / length.Mag();
          else
                  dEdx = 0;
          Double_t gray = dEdx * 1E9 * 1.6E-19 * 1/0.01 * 1/2.4 * 1000; // dEdx [GeV/cm] * conversion to electrons * conversion to Joule * sensitive Area mm2 * density silicon * conversion 1/g to 1/kg= Joule / kg

          TString detname = fGeoH->GetPath(mcPoint->GetSensorID());

          if (fMapDetHistos[detname.Data()] == 0){
                  PndStringSeparator svec(detname.Data(),"/");
                  TVector3 sensDim = fGeoH->GetSensorDimensionsShortId(mcPoint->GetSensorID());
                  std::string histoName = "";
                  histoName += svec.Replace("/","o");
                  fMapDetHistos[detname.Data()] = new TProfile2D(histoName.c_str(),
                                                             detname.Data(),
                                                             (Int_t)(2*sensDim.X()*10),-sensDim.X(),sensDim.X(),        // point resolution mm^2
                                                             (Int_t)(2*sensDim.Y()*10),-sensDim.Y(),sensDim.Y());
          }
          //std::cout << "Damage: " << gray << " gray" << std::endl;
          TVector3 localHit = fGeoH->MasterToLocalShortId(mcPoint->GetPosition(), mcPoint->GetSensorID());
          ((TProfile2D*)(fMapDetHistos[detname.Data()]))->Fill(localHit.X(), localHit.Y(), gray);
          fRadDamHisto->Fill(gray);

//        fRadDamHits = FairRootManager::Instance()->GetTClonesArray("MVDRadDamHit");
//        PndMvdRadDamHit* tempHit = new ((*fRadDamHits)[fRadDamHits->GetEntries()]) PndMvdRadDamHit(mcPoint->GetTrackID(), i, mcPoint->GetSensorID(), 0, 0, mcPoint->GetPosition(), mom, dEdx);
        //  std::cout << "TempHit: " << *tempHit << std::endl;
  }
}

void PndMvdRadDamIonizingTask::FinishEvent ( )

virtual

Definition at line 123 of file PndMvdRadDamIonizingTask.cxx.

{
          //fRadDamHits->Delete();
}

void PndMvdRadDamIonizingTask::FinishTask ( )

virtual

Definition at line 128 of file PndMvdRadDamIonizingTask.cxx.

References fMapDetHistos, and fRadDamHisto.

{
        for (std::map<std::string, TProfile2D*>::iterator it = fMapDetHistos.begin(); it != fMapDetHistos.end(); it++)
                it->second->Write();
        fRadDamHisto->Write();

        //fRadDamHisto->Draw();
}

Bool_t PndMvdRadDamIonizingTask::GetPersistance ( )

inline

Definition at line 38 of file PndMvdRadDamIonizingTask.h.

References fPersistance.

{return fPersistance;};

InitStatus PndMvdRadDamIonizingTask::Init ( )

virtual

Definition at line 48 of file PndMvdRadDamIonizingTask.cxx.

References fMCHits, and fRadDamHisto.

{
        FairRootManager* ioman = FairRootManager::Instance();
  
        if (!ioman)
        {
                std::cout << "-E- PndMvdRadDamIonizingTask::Init: "
    << "RootManager not instantiated!" << std::endl;
                return kFATAL;
        }
  
        fMCHits = (TClonesArray*) ioman->GetObject("MVDPoint");
        if (!fMCHits)
        {
                std::cout << "-W- PndMvdRadDamIonizingTask::Init: " << "No MVDPoint array!"
    << std::endl;
                return kERROR;
        }
  
        // Create and register output array
//      fRadDamHits = ioman->Register("MVDRadDamHit", "PndMvdRadDamHit", "MVD", kTRUE);
  

        fRadDamHisto = new TH1D("radDamH","absorbed dose", 10000000,1E-12,1E-4);
  
        std::cout << "-I- PndMvdRadDamIonizingTask: Initialization successful" << std::endl;
  
        return kSUCCESS;
}

PndMvdRadDamIonizingTask& PndMvdRadDamIonizingTask::operator= ( const PndMvdRadDamIonizingTask &  )

delete

InitStatus PndMvdRadDamIonizingTask::ReInit ( )

virtual

Definition at line 42 of file PndMvdRadDamIonizingTask.cxx.

References SetParContainers().

{
  SetParContainers();
  return kSUCCESS;
}

void PndMvdRadDamIonizingTask::SetParContainers ( )

virtual

Definition at line 38 of file PndMvdRadDamIonizingTask.cxx.

Referenced by ReInit().

{
}

void PndMvdRadDamIonizingTask::SetPersistance ( Bool_t  p = kTRUE )

inline

Definition at line 37 of file PndMvdRadDamIonizingTask.h.

References fPersistance, and p.

{fPersistance=p;};

Member Data Documentation

PndGeoHandling* PndMvdRadDamIonizingTask::fGeoH

private

Definition at line 45 of file PndMvdRadDamIonizingTask.h.

Referenced by Exec().

std::map<std::string, TProfile2D*> PndMvdRadDamIonizingTask::fMapDetHistos

private

Definition at line 47 of file PndMvdRadDamIonizingTask.h.

Referenced by Exec(), FinishTask(), and ~PndMvdRadDamIonizingTask().

TClonesArray* PndMvdRadDamIonizingTask::fMCHits

private

Definition at line 42 of file PndMvdRadDamIonizingTask.h.

Referenced by Exec(), and Init().

Bool_t PndMvdRadDamIonizingTask::fPersistance

private

Definition at line 38 of file PndMvdRadDamIonizingTask.h.

Referenced by GetPersistance(), and SetPersistance().

TH1D* PndMvdRadDamIonizingTask::fRadDamHisto

private

Definition at line 48 of file PndMvdRadDamIonizingTask.h.

Referenced by Exec(), FinishTask(), Init(), and ~PndMvdRadDamIonizingTask().

TClonesArray* PndMvdRadDamIonizingTask::fRadDamHits

private

Definition at line 43 of file PndMvdRadDamIonizingTask.h.

---

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

• PndMvdRadDamIonizingTask.h
• PndMvdRadDamIonizingTask.cxx