PndLmdDigiQTask Class Reference

#include <PndLmdDigiQTask.h>

Inheritance diagram for PndLmdDigiQTask:

Public Member Functions
	PndLmdDigiQTask (const PndLmdDigiQTask &)=delete
	PndLmdDigiQTask (TString digiBranch="LMDPixelDigis")
	~PndLmdDigiQTask ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()
void	SetWriteMC (bool wr)

Private Member Functions
	ClassDef (PndLmdDigiQTask, 2)

Private Attributes
TString	fDigiName
TClonesArray *	fMCHits
TClonesArray *	fMCTracks
TClonesArray *	fDigis
TClonesArray *	fDigiQ
int	fEventNr
bool	fWriteAllMC
PndLmdGeometryHelper &	lmd_geo_helper

Detailed Description

Definition at line 32 of file PndLmdDigiQTask.h.

Constructor & Destructor Documentation

PndLmdDigiQTask::PndLmdDigiQTask ( const PndLmdDigiQTask &  )

delete

Default constructor

PndLmdDigiQTask::PndLmdDigiQTask

(

TString

digiBranch = "LMDPixelDigis"

)

Definition at line 47 of file PndLmdDigiQTask.cxx.

References fDigiName, and fWriteAllMC.

                                                   :
                FairTask("Digis Quality Task for PANDA Lmd"), fEventNr(0), lmd_geo_helper(
                                PndLmdGeometryHelper::getInstance()) {
        fWriteAllMC = false;
        fDigiName = digiBranch;
}

PndLmdDigiQTask::~PndLmdDigiQTask

(

)

Destructor

Definition at line 55 of file PndLmdDigiQTask.cxx.

                                  {
}

Member Function Documentation

PndLmdDigiQTask::ClassDef ( PndLmdDigiQTask  , 2    )

private

void PndLmdDigiQTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Set elastic/inelastic flag ----------------------------------------------——

---

Definition at line 110 of file PndLmdDigiQTask.cxx.

References fabs(), fDigiQ, fDigis, fEventNr, fMCHits, fMCTracks, fVerbose, PndLmdGeometryHelper::getHitLocationInfo(), PndSdsDigi::GetIndex(), PndMCTrack::GetMomentum(), PndMCTrack::GetMotherID(), PndSdsDigi::GetNIndices(), PndMCTrack::GetPdgCode(), PndSdsDigi::GetSensorID(), i, lmd_geo_helper, MCPoint, PndLmdDigiQ::SetEvT(), PndLmdDigiQ::SetFlBkg(), PndLmdDigiQ::SetFlSig(), PndLmdDigiQ::SetHalf(), PndLmdDigiQ::SetMCtrkID(), PndLmdDigiQ::SetModule(), PndLmdDigiQ::SetPDG(), PndLmdDigiQ::SetPlane(), PndLmdDigiQ::SetSide(), and PndLmdDigiQ::SetThMC().

                                    {
        FairRootManager* ioman = FairRootManager::Instance();
        double glEvTime = ioman->GetEventTime();
        fDigiQ->Delete();
        //const int nMCHits = fMCHits->GetEntriesFast(); //[R.K. 01/2017] unused variable
        const int nDigis = fDigis->GetEntriesFast();
        const int nParticles = fMCTracks->GetEntriesFast();

        bool elfl = false;
        int sumID = 0;
        //int TotCharge=0; //[R.K. 01/2017] unused variable
        for (Int_t iN = 0; iN < nParticles; iN++) {

                PndMCTrack *mctrk = (PndMCTrack*) fMCTracks->At(iN);
                //int glPDG = mctrk->GetPdgCode(); //[R.K. 01/2017] unused variable
                Int_t mcID = mctrk->GetPdgCode();
                int motherid = mctrk->GetMotherID();
                if (motherid < 0) {
                        sumID += fabs(mcID);
                }
        }
        if (sumID == 4424)
                elfl = true;

        for (Int_t i = 0; i < nDigis; i++) {
                PndSdsDigiPixel* DigiPoint = (PndSdsDigiPixel*) (fDigis->At(i)); // read digi hit
                TClonesArray& clref = *fDigiQ;
                Int_t size = clref.GetEntriesFast();
                PndLmdDigiQ *digiq = new (clref[size]) PndLmdDigiQ(*DigiPoint, elfl);
                int nmchits = DigiPoint->GetNIndices(); // number of PndSdsMCPoint(s) contributing to digi
                int sensorID = DigiPoint->GetSensorID();
                auto const& digi_info(lmd_geo_helper.getHitLocationInfo(sensorID));

                for (Int_t j = 0; j < nmchits; j++) {
                        int mcind = DigiPoint->GetIndex(j);
                        PndSdsMCPoint* MCPoint = (PndSdsMCPoint*) (fMCHits->At(mcind));
                        int MCtrk = MCPoint->GetTrackID();
                        PndMCTrack *mctrk = (PndMCTrack*) fMCTracks->At(MCtrk);
                        int glPDG = mctrk->GetPdgCode(); // barp: -2212
                        digiq->SetMCtrkID(MCtrk);
                        digiq->SetPDG(glPDG);
                        digiq->SetPlane(digi_info.plane);
                        digiq->SetHalf(digi_info.detector_half);
                        digiq->SetModule(digi_info.module);
                        digiq->SetSide(digi_info.module_side);
                        digiq->SetEvT(glEvTime);
                        int movID = mctrk->GetMotherID(); // primary: <0
                        TVector3 MomMC = mctrk->GetMomentum();
                        double thMC = MomMC.Theta();
                        digiq->SetThMC(thMC);
                        if (glPDG == -2212 && movID < 0 && elfl)
                                digiq->SetFlSig();
                        else
                                digiq->SetFlBkg();
                }
        }
        if (fVerbose > 2)
                cout << "PndLmdDigiQTask::Exec END!" << endl;
        fEventNr++;
}

void PndLmdDigiQTask::Finish ( )

virtual

Definition at line 173 of file PndLmdDigiQTask.cxx.

                             {
}

InitStatus PndLmdDigiQTask::Init ( )

virtual

Virtual method Init

Definition at line 59 of file PndLmdDigiQTask.cxx.

References fDigiQ, fDigis, fMCHits, and fMCTracks.

Referenced by runLumiPixel8DigiQA().

                                 {

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

        // fMCHeader = (TClonesArray*) ioman->GetObject("MCEventHeader.");
        // if ( !fMCHeader)     {
        //   std::cout << "-W- PndLmdDigiQTask::Init: "<< "No MCEventHeader "<<" object!" << std::endl;
        //   return kERROR;
        // }

        //Get MC points
        fMCHits = (TClonesArray*) ioman->GetObject("LMDPoint");
        if (!fMCHits) {
                std::cout << "-W- PndLmdDigiQTask::Init: " << "No LMDPoint" << " array!"
                                << std::endl;
                return kERROR;
        }

        //Get MC tracks
        fMCTracks = (TClonesArray*) ioman->GetObject("MCTrack");
        if (!fMCTracks) {
                std::cout << "-W- PndLmdDigiQTask::Init: " << "No MCTrack" << " array!"
                                << std::endl;
                return kERROR;
        }

        //Get digis
        fDigis = (TClonesArray*) ioman->GetObject("LMDPixelDigis");
        if (!fDigis) {
                std::cout << "-W- PndLmdDigiQTask::Init: " << "No LMDPixelDigis"
                                << " array!" << std::endl;
                return kERROR;
        }

        fDigiQ = new TClonesArray("PndLmdDigiQ");
        ioman->Register("LMDPixelDigisQ", "DigiQ", fDigiQ, kTRUE);

        //FairRun* fRun = FairRun::Instance(); //[R.K. 01/2017] unused variable?
        //FairRuntimeDb* rtdb = fRun->GetRuntimeDb(); //[R.K. 01/2017] unused variable

        return kSUCCESS;
}

void PndLmdDigiQTask::SetWriteMC ( bool  wr )

inline

Definition at line 56 of file PndLmdDigiQTask.h.

References fWriteAllMC.

{fWriteAllMC = wr;}

Member Data Documentation

TString PndLmdDigiQTask::fDigiName

private

Definition at line 60 of file PndLmdDigiQTask.h.

Referenced by PndLmdDigiQTask().

TClonesArray* PndLmdDigiQTask::fDigiQ

private

Definition at line 71 of file PndLmdDigiQTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndLmdDigiQTask::fDigis

private

Definition at line 63 of file PndLmdDigiQTask.h.

Referenced by Exec(), and Init().

int PndLmdDigiQTask::fEventNr

private

Definition at line 73 of file PndLmdDigiQTask.h.

Referenced by Exec().

TClonesArray* PndLmdDigiQTask::fMCHits

private

Definition at line 61 of file PndLmdDigiQTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndLmdDigiQTask::fMCTracks

private

Definition at line 62 of file PndLmdDigiQTask.h.

Referenced by Exec(), and Init().

bool PndLmdDigiQTask::fWriteAllMC

private

Definition at line 78 of file PndLmdDigiQTask.h.

Referenced by PndLmdDigiQTask(), and SetWriteMC().

PndLmdGeometryHelper& PndLmdDigiQTask::lmd_geo_helper

private

Definition at line 80 of file PndLmdDigiQTask.h.

Referenced by Exec().

---

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

• PndLmdDigiQTask.h
• PndLmdDigiQTask.cxx