PndLmdNoiseTrkSuppressionTask Class Reference

#include <PndLmdNoiseTrkSuppressionTask.h>

Inheritance diagram for PndLmdNoiseTrkSuppressionTask:

Public Member Functions
	PndLmdNoiseTrkSuppressionTask ()
	PndLmdNoiseTrkSuppressionTask (Double_t pBeam, TString dir)
	Set up beam momuntum value. More...
	~PndLmdNoiseTrkSuppressionTask ()
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)

Private Member Functions
	ClassDef (PndLmdNoiseTrkSuppressionTask, 1)

Private Attributes
TClonesArray *	fTrkInArray
TClonesArray *	fTrkOutArray
int	fEventNr
Double_t	fPbeam
TString	fdir
TString	fmethodName
Float_t	axrec
Float_t	ayrec
Float_t	athrec
Float_t	aphrec
TMVA::Reader *	reader

Detailed Description

PndLmdNoiseTrkSuppressionTask.h

Author

A.Karavdina karav.nosp@m.din@.nosp@m.kph.u.nosp@m.ni-m.nosp@m.ainz..nosp@m.de

Clean sample with reconstructed tracks before back propagation from background due to noise hits by multivariable cut

Definition at line 22 of file PndLmdNoiseTrkSuppressionTask.h.

Constructor & Destructor Documentation

PndLmdNoiseTrkSuppressionTask::PndLmdNoiseTrkSuppressionTask

(

)

Default constructor

Definition at line 28 of file PndLmdNoiseTrkSuppressionTask.cxx.

    : FairTask("Cleaning Tracks from noise hits Task for PANDA Lmd"),
      fEventNr(0) {
  // tprop = new TNtuple();
  // fEventNr = 0;
}

PndLmdNoiseTrkSuppressionTask::PndLmdNoiseTrkSuppressionTask	(	Double_t	pBeam,
		TString	dir
	)

Set up beam momuntum value.

Definition at line 36 of file PndLmdNoiseTrkSuppressionTask.cxx.

References fdir, and fPbeam.

    : FairTask("Cleaning Tracks from noise hits Task for PANDA Lmd"),
      fEventNr(0) {
  fdir = dir;
  fPbeam = pBeam;
  std::cout << "Beam Momentum in this run is " << fPbeam << std::endl;
}

PndLmdNoiseTrkSuppressionTask::~PndLmdNoiseTrkSuppressionTask

(

)

Destructor

Definition at line 46 of file PndLmdNoiseTrkSuppressionTask.cxx.

{}

Member Function Documentation

PndLmdNoiseTrkSuppressionTask::ClassDef ( PndLmdNoiseTrkSuppressionTask  , 1    )

private

void PndLmdNoiseTrkSuppressionTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 113 of file PndLmdNoiseTrkSuppressionTask.cxx.

References aphrec, athrec, axrec, ayrec, fEventNr, fmethodName, fTrkInArray, fTrkOutArray, fVerbose, PndTrack::GetParamFirst(), and reader.

                                                  {
  fTrkOutArray->Delete();
  if (fVerbose > 2) {
    std::cout << " ---- Info: " << fEventNr << std::endl;
  }

  // go through all tracks
  const int nTrks = fTrkInArray->GetEntries();
  int rec_trk = 0;
  for (Int_t iN = 0; iN < nTrks; iN++) {  // loop over all reconstructed trks
    PndTrack* trkpnd = (PndTrack*)(fTrkInArray->At(iN));
    FairTrackParP fFittedTrkP = trkpnd->GetParamFirst();

    // check the track
    axrec = fFittedTrkP.GetX();
    ayrec = fFittedTrkP.GetY();
    TVector3 MomRecLMD(fFittedTrkP.GetPx(), fFittedTrkP.GetPy(),
                       fFittedTrkP.GetPz());
    MomRecLMD *= 1. / MomRecLMD.Mag();
    athrec = MomRecLMD.Theta();
    aphrec = MomRecLMD.Phi();
    double mva_response = reader->EvaluateMVA(fmethodName);
    if (fVerbose > 2)
      std::cout << "mva_response  = " << mva_response << std::endl;
    bool isSigTrk;
    if (mva_response < 0)
      isSigTrk = false;
    else
      isSigTrk = true;
    if (isSigTrk) {
      // save good trks
      new ((*fTrkOutArray)[rec_trk]) PndTrack(*(trkpnd));  // save Track
      rec_trk++;
    }
  }
  if (fVerbose > 2)
    std::cout << "Ev#" << fEventNr << ": " << rec_trk << " trks saved out of "
              << nTrks << std::endl;
  fEventNr++;
}

InitStatus PndLmdNoiseTrkSuppressionTask::Init ( )

virtual

Definition at line 49 of file PndLmdNoiseTrkSuppressionTask.cxx.

References aphrec, athrec, axrec, ayrec, fdir, fmethodName, fTrkInArray, fTrkOutArray, reader, and TString.

Referenced by runLumiPixel4bNoiseTrkFilter().

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

  fTrkInArray = (TClonesArray*)ioman->GetObject("LMDPndTrackFilt");
  if (!fTrkInArray) {
    std::cout << "-W- PndLmdNoiseTrkSuppressionTask::Init: "
              << "No LMDPndTrack"
              << " array!" << std::endl;
    return kERROR;
  }

  fTrkOutArray = new TClonesArray("PndTrack");
  ioman->Register("LMDPndTrackFilt", "PndLmd", fTrkOutArray, kTRUE);

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

  // TMVA -----------------------------------------------------
  // This loads the library
  // if(fabs(fPbeam-1.5)<1e-1 || fabs(fPbeam-15)<1e-1){
  TMVA::Tools::Instance();

  //   TMVA::Reader *
  reader = new TMVA::Reader("!Color:!Silent");

  // Create a set of variables and declare them to the reader
  // - the variable names MUST corresponds in name and type to those given in
  // the weight file(s) used
  reader->AddVariable("LMDTrackQ.fThetarecLMD", &athrec);
  reader->AddVariable("LMDTrackQ.fPhirecLMD", &aphrec);
  reader->AddVariable("LMDTrackQ.fXrecLMD", &axrec);
  reader->AddVariable("LMDTrackQ.fYrecLMD", &ayrec);

  //  TString dir    = "weights/";
  TString prefix = "TMVAClassification";
  fmethodName = "BDT method";
  TString weightfile =
      fdir + prefix + TString("_BDT") + TString(".weights.xml");
  reader->BookMVA(fmethodName, weightfile);
  // }
  //   else{
  // std::cout << "-W- PndLmdNoiseTrkSuppressionTask::Init: "<< "This method was
  // not trained for momentum"<< fPbeam << std::endl;
  // return kERROR;
  // }
  //------------------------------------------------------
  return kSUCCESS;
}

void PndLmdNoiseTrkSuppressionTask::SetParContainers ( )

virtual

Virtual method Init

FairRun* ana = FairRun::Instance();

Definition at line 106 of file PndLmdNoiseTrkSuppressionTask.cxx.

                                                     {
  // Get Base Container
  // FairRuntimeDb* rtdb=ana->GetRuntimeDb();
}

Member Data Documentation

Float_t PndLmdNoiseTrkSuppressionTask::aphrec

private

Definition at line 53 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec(), and Init().

Float_t PndLmdNoiseTrkSuppressionTask::athrec

private

Definition at line 53 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec(), and Init().

Float_t PndLmdNoiseTrkSuppressionTask::axrec

private

Definition at line 52 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec(), and Init().

Float_t PndLmdNoiseTrkSuppressionTask::ayrec

private

Definition at line 52 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec(), and Init().

TString PndLmdNoiseTrkSuppressionTask::fdir

private

Definition at line 50 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Init(), and PndLmdNoiseTrkSuppressionTask().

int PndLmdNoiseTrkSuppressionTask::fEventNr

private

Definition at line 46 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec().

TString PndLmdNoiseTrkSuppressionTask::fmethodName

private

Definition at line 51 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec(), and Init().

Double_t PndLmdNoiseTrkSuppressionTask::fPbeam

private

Definition at line 47 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by PndLmdNoiseTrkSuppressionTask().

TClonesArray* PndLmdNoiseTrkSuppressionTask::fTrkInArray

private

Definition at line 42 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndLmdNoiseTrkSuppressionTask::fTrkOutArray

private

Definition at line 43 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec(), and Init().

TMVA::Reader* PndLmdNoiseTrkSuppressionTask::reader

private

Definition at line 54 of file PndLmdNoiseTrkSuppressionTask.h.

Referenced by Exec(), and Init().

---

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

• PndLmdNoiseTrkSuppressionTask.h
• PndLmdNoiseTrkSuppressionTask.cxx