PndMvdRadDamTask Class Reference

#include <PndMvdRadDamTask.h>

Inheritance diagram for PndMvdRadDamTask:

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

Private Member Functions
void	InitWeightLists ()
	ClassDef (PndMvdRadDamTask, 1)

Private Attributes
Bool_t	fPersistance
TClonesArray *	fMCTracks
TClonesArray *	fMCHits
TClonesArray *	fRadDamHits
PndMvdRadDamList *	fElectronList
PndMvdRadDamList *	fProtonList
PndMvdRadDamList *	fNeutronList
PndMvdRadDamList *	fPionList
PndGeoHandling *	fGeoH
std::map< std::string, TH2 * >	fMapDetHistos
TH1D *	fRadDamHisto
std::map< Int_t, PndMvdRadDamList * >	fWeightListsMap

Detailed Description

Definition at line 21 of file PndMvdRadDamTask.h.

Constructor & Destructor Documentation

PndMvdRadDamTask::PndMvdRadDamTask

(

)

Definition at line 19 of file PndMvdRadDamTask.cxx.

                                   :
    fPersistance(kTRUE),
          fMCTracks(NULL),
          fMCHits(NULL),
          fRadDamHits(NULL),
          fElectronList(NULL),
          fProtonList(NULL),
          fNeutronList(NULL),
          fPionList(NULL),
          fGeoH(PndGeoHandling::Instance()),
          fMapDetHistos(),
          fRadDamHisto(NULL),
          fWeightListsMap()
{
}

PndMvdRadDamTask::~PndMvdRadDamTask

(

)

Definition at line 35 of file PndMvdRadDamTask.cxx.

References fMapDetHistos, and fRadDamHisto.

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

PndMvdRadDamTask::PndMvdRadDamTask ( const PndMvdRadDamTask &  )

delete

Member Function Documentation

PndMvdRadDamTask::ClassDef ( PndMvdRadDamTask  , 1    )

private

void PndMvdRadDamTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 115 of file PndMvdRadDamTask.cxx.

References detname, Double_t, fGeoH, fMapDetHistos, fMCHits, fMCTracks, fRadDamHisto, fRadDamHits, fWeightListsMap, PndGeoHandling::GetPath(), PndMCTrack::GetPdgCode(), PndSdsMCPoint::GetPosition(), PndGeoHandling::GetSensorDimensionsShortId(), PndSdsMCPoint::GetSensorID(), i, PndGeoHandling::MasterToLocalShortId(), mom, pid(), and TString.

{
  // Reset output array
  if ( ! fRadDamHits )
    Fatal("Exec", "No RadDamArray");
  fRadDamHits->Delete();
  
  PndSdsMCPoint* mcPoint;
  PndMCTrack* mcTrack;
  
  for (int i = 0; i < fMCHits->GetEntriesFast(); i++){
          mcPoint = (PndSdsMCPoint*)fMCHits->At(i);
          mcTrack = (PndMCTrack*)(fMCTracks->At(mcPoint->GetTrackID()));
          TVector3 mom;
          mcPoint->Momentum(mom);
          Double_t Ekin = mom.Mag();
          Int_t pid = mcTrack->GetPdgCode();
          Double_t weight = 0;
          if (fWeightListsMap[pid] != NULL){
                  weight = fWeightListsMap[pid]->GetWeight(Ekin);
          }
          //std::cout << "WeightCalc: pid: " << pid << " Energy: " << Ekin << " Weight: " << weight << std::endl;
          new ((*fRadDamHits)[i]) PndMvdRadDamHit(mcPoint->GetTrackID(), i, mcPoint->GetSensorID(), pid, Ekin,
                                            mcPoint->GetPosition(), mom, weight);
    TString detname = fGeoH->GetPath(mcPoint->GetSensorID());
          if (fMapDetHistos[detname.Data()] == 0){
                  PndStringSeparator svec(detname.Data(),"/");
                  TVector3 sensDim = fGeoH->GetSensorDimensionsShortId(mcPoint->GetSensorID());
                  fMapDetHistos[detname.Data()] = new TH2D(svec.Replace("/","o").c_str(),
                                                             fGeoH->GetPath(mcPoint->GetSensorID()),
                                                             (Int_t)(2*sensDim.X()*10),-sensDim.X(),sensDim.X(),        // point resolution mm^2
                                                             (Int_t)(2*sensDim.Y()*10),-sensDim.Y(),sensDim.Y());
          }
          TVector3 localHit = fGeoH->MasterToLocalShortId(mcPoint->GetPosition(), mcPoint->GetSensorID());
          ((TH2D*)(fMapDetHistos[detname.Data()]))->Fill(localHit.X(), localHit.Y(), weight);
          fRadDamHisto->Fill(weight);
  }
}

void PndMvdRadDamTask::Finish ( )

virtual

Definition at line 154 of file PndMvdRadDamTask.cxx.

References fMapDetHistos, and fRadDamHisto.

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

Bool_t PndMvdRadDamTask::GetPersistance ( )

inline

Definition at line 39 of file PndMvdRadDamTask.h.

References fPersistance.

{return fPersistance;};

InitStatus PndMvdRadDamTask::Init ( )

virtual

Definition at line 54 of file PndMvdRadDamTask.cxx.

References fMCHits, fMCTracks, fPersistance, fRadDamHisto, fRadDamHits, and InitWeightLists().

{
        FairRootManager* ioman = FairRootManager::Instance();
  
        if (!ioman)
        {
                std::cout << "-E- PndMvdRadDamTask::Init: "
    << "RootManager not instantiated!" << std::endl;
                return kFATAL;
        }
  
        fMCTracks = (TClonesArray*) ioman->GetObject("MCTrack");
        if (!fMCTracks)
        {
                std::cout << "-W- PndMvdRadDamTask::Init: " << "No MCTrack array!"
    << std::endl;
                return kERROR;
        }
  
        fMCHits = (TClonesArray*) ioman->GetObject("MVDPoint");
        if (!fMCHits)
        {
                std::cout << "-W- PndMvdRadDamTask::Init: " << "No MVDPoint array!"
    << std::endl;
                return kERROR;
        }
  
        // Create and register output array
        fRadDamHits = new TClonesArray("PndMvdRadDamHit");
        ioman->Register("MVDRadDamHit", "MVD", fRadDamHits, fPersistance);
  
        InitWeightLists();
        fRadDamHisto = new TH1D("radDamH","Weight Factors", 1000,0,100);
  
        std::cout << "-I- PndMvdRadDamTask: Initialization successful" << std::endl;
  
        return kSUCCESS;
}

void PndMvdRadDamTask::InitWeightLists ( )

private

Definition at line 93 of file PndMvdRadDamTask.cxx.

References fElectronList, fNeutronList, fPionList, fProtonList, and fWeightListsMap.

Referenced by Init().

{
        fElectronList = new PndMvdRadDamList("$VMCWORKDIR/pandaroot/mvd/MvdTools/MvdRadDamage/electronsWeight.root");
        fProtonList = new PndMvdRadDamList("$VMCWORKDIR/pandaroot/mvd/MvdTools/MvdRadDamage/protonsWeight.root");
        fNeutronList = new PndMvdRadDamList("$VMCWORKDIR/pandaroot/mvd/MvdTools/MvdRadDamage/neutronsWeight.root");
        fPionList = new PndMvdRadDamList("$VMCWORKDIR/pandaroot/mvd/MvdTools/MvdRadDamage/pionsWeight.root");
  
        fWeightListsMap[11] = fElectronList;            //e-
        fWeightListsMap[-11] = fElectronList;           //e+
        fWeightListsMap[2212] = fProtonList;            //p+
        fWeightListsMap[-2212] = fProtonList;           //p-
        fWeightListsMap[2112] = fNeutronList;           //n
        fWeightListsMap[-2112] = fNeutronList;          //n_bar
        fWeightListsMap[111] = fPionList;                       //pi0
        fWeightListsMap[211] = fPionList;                       //pi+
        fWeightListsMap[-211] = fPionList;                      //pi-
        fWeightListsMap[1000010030] = fProtonList;      //Tritium
        fWeightListsMap[1000010020] = fProtonList;      //Deuterium
  
  
}

PndMvdRadDamTask& PndMvdRadDamTask::operator= ( const PndMvdRadDamTask &  )

delete

InitStatus PndMvdRadDamTask::ReInit ( )

virtual

Definition at line 48 of file PndMvdRadDamTask.cxx.

References SetParContainers().

{
  SetParContainers();
  return kSUCCESS;
}

void PndMvdRadDamTask::SetParContainers ( )

virtual

Definition at line 44 of file PndMvdRadDamTask.cxx.

Referenced by ReInit().

{
}

void PndMvdRadDamTask::SetPersistance ( Bool_t  p = kTRUE )

inline

Definition at line 38 of file PndMvdRadDamTask.h.

References fPersistance, and p.

{fPersistance=p;};

Member Data Documentation

PndMvdRadDamList* PndMvdRadDamTask::fElectronList

private

Definition at line 47 of file PndMvdRadDamTask.h.

Referenced by InitWeightLists().

PndGeoHandling* PndMvdRadDamTask::fGeoH

private

Definition at line 52 of file PndMvdRadDamTask.h.

Referenced by Exec().

std::map<std::string, TH2*> PndMvdRadDamTask::fMapDetHistos

private

Definition at line 54 of file PndMvdRadDamTask.h.

Referenced by Exec(), Finish(), and ~PndMvdRadDamTask().

TClonesArray* PndMvdRadDamTask::fMCHits

private

Definition at line 44 of file PndMvdRadDamTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndMvdRadDamTask::fMCTracks

private

Definition at line 43 of file PndMvdRadDamTask.h.

Referenced by Exec(), and Init().

PndMvdRadDamList* PndMvdRadDamTask::fNeutronList

private

Definition at line 49 of file PndMvdRadDamTask.h.

Referenced by InitWeightLists().

Bool_t PndMvdRadDamTask::fPersistance

private

Definition at line 39 of file PndMvdRadDamTask.h.

Referenced by GetPersistance(), Init(), and SetPersistance().

PndMvdRadDamList* PndMvdRadDamTask::fPionList

private

Definition at line 50 of file PndMvdRadDamTask.h.

Referenced by InitWeightLists().

PndMvdRadDamList* PndMvdRadDamTask::fProtonList

private

Definition at line 48 of file PndMvdRadDamTask.h.

Referenced by InitWeightLists().

TH1D* PndMvdRadDamTask::fRadDamHisto

private

Definition at line 55 of file PndMvdRadDamTask.h.

Referenced by Exec(), Finish(), Init(), and ~PndMvdRadDamTask().

TClonesArray* PndMvdRadDamTask::fRadDamHits

private

Definition at line 45 of file PndMvdRadDamTask.h.

Referenced by Exec(), and Init().

std::map<Int_t, PndMvdRadDamList*> PndMvdRadDamTask::fWeightListsMap

private

Definition at line 57 of file PndMvdRadDamTask.h.

Referenced by Exec(), and InitWeightLists().

---

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

• PndMvdRadDamTask.h
• PndMvdRadDamTask.cxx