PndMdtPointsToWaveform Class Reference

#include <PndMdtPointsToWaveform.h>

Inheritance diagram for PndMdtPointsToWaveform:

Classes
struct	key

Public Member Functions
	PndMdtPointsToWaveform (Int_t verbose=0, Bool_t storewaves=kTRUE)
virtual	~PndMdtPointsToWaveform ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
void	SetStorageOfData (Bool_t val)
void	RunTimeBased ()
void	FinishTask ()
void	SetPersistency (Bool_t val=kTRUE)
Bool_t	GetPersistency ()

Private Member Functions
virtual void	SetParContainers ()
virtual void	exec_t ()
virtual void	exec_e ()
Int_t	PdgToIndex (Int_t pdg)
	PndMdtPointsToWaveform (const PndMdtPointsToWaveform &L)
PndMdtPointsToWaveform &	operator= (const PndMdtPointsToWaveform &)
Bool_t	Digitize (PndMdtWaveform *theWf, Double_t &time, Double_t &amp, Bool_t isWire)
	ClassDef (PndMdtPointsToWaveform, 2)

Private Attributes
TClonesArray *	fMcTrackArray
TClonesArray *	fPointArray
TClonesArray *	fWaveformArray
Bool_t	fTimeOrderedWaveform
PndMdtWaveformWriteoutBuffer *	fDataBuffer
PndMdtParamDigi *	fParamDigiModel
PndMdtIGeometry *	fGeoIF
Int_t	fVerbose
Int_t	HowManyPoint
Int_t	nWaveformProduced
TFile *	fFile
TTree *	tTree
Double_t	fWirpT
Double_t	fStripT
Double_t	fEvtT
Double_t	fLength
Double_t	fDis
Int_t	fMod
Int_t	fPid

Detailed Description

Definition at line 30 of file PndMdtPointsToWaveform.h.

Constructor & Destructor Documentation

PndMdtPointsToWaveform::PndMdtPointsToWaveform	(	Int_t	verbose = 0,
		Bool_t	storewaves = kTRUE
	)

Definition at line 36 of file PndMdtPointsToWaveform.cxx.

References PndPersistencyTask::SetPersistency().

                                                                              :
  //fDigiPar(new PndMdtDigiPar())
  //, fGeoPar(new PndMdtGeoPar())
    fWaveformArray(0)
  , fTimeOrderedWaveform(kFALSE)
  , fDataBuffer(0)
  , fVerbose(verbose)
{
        SetPersistency(storewaves);
}

PndMdtPointsToWaveform::~PndMdtPointsToWaveform ( )

virtual

Definition at line 50 of file PndMdtPointsToWaveform.cxx.

{
}

PndMdtPointsToWaveform::PndMdtPointsToWaveform ( const PndMdtPointsToWaveform &  L )

private

Member Function Documentation

PndMdtPointsToWaveform::ClassDef ( PndMdtPointsToWaveform  , 2    )

private

Bool_t PndMdtPointsToWaveform::Digitize ( PndMdtWaveform *  theWf, Double_t &  time, Double_t &  amp, Bool_t  isWire  )

private

Definition at line 467 of file PndMdtPointsToWaveform.cxx.

References CAMath::Abs(), Bool_t, Double_t, PndMdtWaveform::GetSignal(), and is.

Referenced by exec_e().

{
  const Double_t fWireNoiseSigma = 0.05;
  const Double_t fStripNoiseSigma = 0.01;
  const Double_t fSamplingInterval = 10.;//nano seconds

  time = -1;
  amp = -9999;
  Double_t fPeak = 0.;
  Bool_t NotFound = kTRUE;
  //Double_t fNoiseSigma; //[R.K. 01/2017] unused variable?
  Double_t fThreshold(1e9);
  if(isWire) fThreshold = 5.*fWireNoiseSigma;
  else fThreshold = 5.*fStripNoiseSigma;

  const std::vector<Double_t>& fSignalData = theWf->GetSignal();
  for(size_t is=0; is < fSignalData.size(); ++ is){
    if(TMath::Abs(fSignalData[is]) > TMath::Abs(fPeak)) fPeak = fSignalData[is];
    if(TMath::Abs(fSignalData[is]) > fThreshold && NotFound) {
      time = is*fSamplingInterval;//nano seconds
      NotFound = kFALSE;
    }
  }
  amp = fPeak;
  return time > 0;
}

void PndMdtPointsToWaveform::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 134 of file PndMdtPointsToWaveform.cxx.

References exec_e(), exec_t(), and fTimeOrderedWaveform.

{
  if(fTimeOrderedWaveform)
    exec_t();
  else
    exec_e();
}

void PndMdtPointsToWaveform::exec_e ( )

privatevirtual

Definition at line 141 of file PndMdtPointsToWaveform.cxx.

References CAMath::Abs(), Bool_t, PndMdtParamDigi::Compute(), ctime, Digitize(), Double_t, fDis, fEvtT, fGeoIF, fLength, fMcTrackArray, fMod, fParamDigiModel, fPid, fPointArray, fStripT, fVerbose, fWaveformArray, fWirpT, PndMdtPoint::GetModule(), PndMdtPoint::GetMomIn(), PndMCTrack::GetPdgCode(), PndMdtPoint::GetPosIn(), PndMdtPoint::GetPosOut(), PndMdtPoint::GetSector(), PndMdtParamDigi::GetStripSignals(), PndMdtParamDigi::GetWireSignal(), HowManyPoint, PndMdtIGeometry::MapWireToStrip(), PndMdtIGeometry::MasterToLocal(), nHits, nWaveformProduced, PdgToIndex(), rtime, PndMdtParamDigi::SetParams(), PndMdtWaveform::SetSignal(), CAMath::Sin(), timer, and tTree.

Referenced by Exec().

{
  TStopwatch timer;
  if (fVerbose>2){
    timer.Start();
  }

  fWaveformArray->Clear();

  Double_t fEventTime = FairRootManager::Instance()->GetEventTime();//nano seconds
  Int_t nHits = fPointArray->GetEntriesFast();
  Int_t evtNo = FairRunAna::Instance()->GetEventHeader()->GetMCEntryNumber();
  if (fVerbose>1){
    cout<<"================PndMdtPointsToWaveform==============================================="<<endl;
    cout<<"Event No.: "<<evtNo<<", Event Time: "<<fEventTime<<", Load in "<<nHits<< " PndMdtPoint."<<endl;
  }
  std::map<key, Int_t> fMcTrackingInfo;
  PndMdtPoint* thePoint(0);
  Int_t TrackID;
  Int_t ParticleType;
  Int_t DetectorID;
  TVector3 EntrancePosition;
  TVector3 ExitingPosition;
  TVector3 Momentum;
  Double_t fWaveformTimeStamp;

  Int_t nPoints = fPointArray->GetEntriesFast();
  HowManyPoint +=  nPoints;
  std::map<Int_t, PndMdtWaveform* > listofWaveforms;

  Double_t fAmplitude;
  Double_t fTimeStamp;

  TVector3 fTubeCenter;
  for (Int_t ip=0; ip<nPoints; ip++) 
  {
    thePoint = (PndMdtPoint*) fPointArray->At(ip);
    if(thePoint->GetEnergyLoss() == 0) continue;
    key _TrkandDetID(thePoint->GetTrackID(), thePoint->GetDetectorID());
    //cout<<"TrackID = "<<thePoint->GetTrackID()<<", DetID = "<<thePoint->GetDetectorID()<<endl;
    std::map<key, Int_t>::iterator it = fMcTrackingInfo.find(_TrkandDetID);
    if( fMcTrackingInfo.end() == it )
    {
      fMcTrackingInfo.insert(std::pair<key,Int_t>(_TrkandDetID, 1));
      TrackID = thePoint->GetTrackID();
      DetectorID = thePoint->GetDetectorID();
      //ParticleType = thePoint->GetParticleType();
      EntrancePosition = thePoint->GetPosIn();
      ExitingPosition = thePoint->GetPosOut();
      Momentum = thePoint->GetMomIn();
      fWaveformTimeStamp = thePoint->GetTime() + fEventTime;

      //fxs = EntrancePosition.X();
      //fys = EntrancePosition.Y();
      //fzs = EntrancePosition.Z();
      fLength = thePoint->GetLength();
      fMod = thePoint->GetModule();

      //
      TVector3 fNewEntrance(0,0,0);
      TVector3 fNewExit(0,0,0);
      Bool_t fOK1 = fGeoIF->MasterToLocal(DetectorID, EntrancePosition, fNewEntrance);
      Bool_t fOK2 = fGeoIF->MasterToLocal(DetectorID, ExitingPosition, fNewExit);
      if( ! (fOK1 && fOK2) ) {
        cout <<"PndMdtPointsToWaveform::exec_e, invliad detector index #"<<DetectorID<<endl;
        continue; 
      }
      TVector2 fexit = fNewExit.XYvector();
      TVector2 fentr = fNewEntrance.XYvector();
      TVector2 fvd = fexit - fentr;

      fDis = TMath::Abs(fexit.Mod()*fentr.Mod()*TMath::Sin(fexit.DeltaPhi(fentr)))/fvd.Mod();
      if(fVerbose >1 ){
        cout<<"================================================"<<endl;
        cout<<"Mod = "<<thePoint->GetModule()<<endl;
        cout<<"Oct = "<<thePoint->GetSector()<<endl
          <<"Entrance "<<EntrancePosition<<endl
          <<"Exiting "<<ExitingPosition<<endl
          <<"new Entrance "<<fNewEntrance<<endl
          <<"new Exiting "<<fNewExit<<endl;
      }
      PndMCTrack *mcTrk = (PndMCTrack*)fMcTrackArray->At(thePoint->GetTrackID());
      ParticleType= PdgToIndex(mcTrk->GetPdgCode());
      fPid = ParticleType;

      fParamDigiModel->SetParams(ParticleType
          , Momentum
          , fNewEntrance 
          , fNewExit);
      //cout<<"================================================"<<endl;
      fParamDigiModel->Compute();

      PndMdtWaveform* aWf= new PndMdtWaveform(TrackID, DetectorID, fWaveformTimeStamp, kTRUE);
      //++nWaveformProduced;
      aWf->SetSignal(fParamDigiModel->GetWireSignal());

      std::map<Int_t, PndMdtWaveform*>::iterator mit3 = listofWaveforms.find(DetectorID);
      if( mit3 == listofWaveforms.end())
      {//the first waveform
        listofWaveforms.insert(std::pair<Int_t, PndMdtWaveform*>(DetectorID, aWf));
      }else{
        *(mit3->second) += (*aWf);
        //delete aWf;
      }

      Digitize(aWf, fTimeStamp, fAmplitude, kTRUE);
      fWirpT = fTimeStamp + fEventTime;

      const std::map<Int_t, std::vector<Double_t> >&  fStripSignals = fParamDigiModel->GetStripSignals();

      std::map<Int_t, std::vector<Double_t> >::const_iterator mit = fStripSignals.begin();
      std::map<Int_t, std::vector<Double_t> >::const_iterator mend = fStripSignals.end();
      Int_t newDetID;
      for(; mit != mend; ++mit){
        //re-make detector id;
        fGeoIF->MapWireToStrip(DetectorID, fNewEntrance, newDetID);//waiting for updates
        Int_t lid = mit->first;
        if(lid%2 == 0) newDetID += lid/2;
        else newDetID -= lid/2;
        PndMdtWaveform* aSf= new PndMdtWaveform(TrackID, newDetID, fWaveformTimeStamp, kFALSE);
        //++nWaveformProduced;
        aSf->SetSignal(mit->second);

        Digitize(aSf, fTimeStamp, fAmplitude, kFALSE);
        fStripT = fTimeStamp + fEventTime;
        fEvtT = fEventTime;

        std::map<Int_t, PndMdtWaveform*>::iterator mit1 = listofWaveforms.find(newDetID);
        if( mit1 == listofWaveforms.end())
        {//the first waveform
          listofWaveforms.insert(std::pair<Int_t, PndMdtWaveform*>(newDetID, aSf));
        }else{
          *(mit1->second) += (*aSf);
          //delete aSf;
        }
      }
      tTree->Fill();
    }else{
      ++ it->second;
    }
  }

  std::map<Int_t, PndMdtWaveform*>::iterator mit2 = listofWaveforms.begin();
  for(; mit2 != listofWaveforms.end(); ++ mit2){
    (*fWaveformArray)[fWaveformArray->GetEntriesFast()] = mit2->second;
    ++ nWaveformProduced;
  }

  if(fVerbose > 1){
    std::map<key, Int_t>::iterator it = fMcTrackingInfo.begin();
    std::map<key, Int_t>::iterator end = fMcTrackingInfo.end();
    for(; it != end; ++it){
      cout<<"Track #"<<it->first.TrkID<<" contributes "<<it->second<<" points in tube "<<it->first.DetID<<endl;
    }
  }
  if(fVerbose > 1){
    cout<<"Number of produced waveforms: "<< fWaveformArray->GetEntriesFast()<<endl;
    cout<<"================PndMdtPointsToWaveform==============================================="<<endl;
  }

  if (fVerbose>2){
    timer.Stop();
    Double_t rtime = timer.RealTime();
    Double_t ctime = timer.CpuTime();
    cout << "PndMdtPointsToWaveform, Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
  }
}

void PndMdtPointsToWaveform::exec_t ( )

privatevirtual

Definition at line 310 of file PndMdtPointsToWaveform.cxx.

References CAMath::Abs(), Bool_t, PndMdtParamDigi::Compute(), ctime, Double_t, fDataBuffer, fDis, fGeoIF, fMcTrackArray, fParamDigiModel, fPid, fPointArray, fVerbose, PndMdtWaveform::GetActiveTime(), PndMdtPoint::GetModule(), PndMdtPoint::GetMomIn(), PndMCTrack::GetPdgCode(), PndMdtPoint::GetPosIn(), PndMdtPoint::GetPosOut(), PndMdtPoint::GetSector(), PndMdtParamDigi::GetStripSignals(), PndMdtParamDigi::GetWireSignal(), HowManyPoint, PndMdtIGeometry::MapWireToStrip(), PndMdtIGeometry::MasterToLocal(), nHits, nWaveformProduced, PdgToIndex(), rtime, PndMdtParamDigi::SetParams(), PndMdtWaveform::SetSignal(), CAMath::Sin(), and timer.

Referenced by Exec().

{
  TStopwatch timer;
  if (fVerbose>2){
    timer.Start();
  }

  Double_t fEventTime = FairRootManager::Instance()->GetEventTime();//nano seconds
  Int_t nHits = fPointArray->GetEntriesFast();
  Int_t evtNo = FairRunAna::Instance()->GetEventHeader()->GetMCEntryNumber();
  if (fVerbose>1){
    cout<<"================PndMdtPointsToWaveform==============================================="<<endl;
    cout<<"Event No.: "<<evtNo<<", Event Time: "<<fEventTime<<", Load in "<<nHits<< " PndMdtHit."<<endl;
    cout<<"Size of PndMdtWaveformWriteoutBuffer: "<< fDataBuffer->GetNData()<<endl;
    cout<<"================PndMdtPointsToWaveform==============================================="<<endl;
  }
  std::map<key, Int_t> fMcTrackingInfo;
  PndMdtPoint* thePoint(0);
  Int_t TrackID;
  Int_t ParticleType;
  Long_t DetectorID;
  TVector3 EntrancePosition;
  TVector3 ExitingPosition;
  TVector3 Momentum;
  Double_t TimeStamp;

  Int_t nPoints = fPointArray->GetEntriesFast();
  HowManyPoint +=  nPoints;
  for (Int_t ip=0; ip<nPoints; ip++) 
  {
    thePoint = (PndMdtPoint*) fPointArray->At(ip);
    if(thePoint->GetEnergyLoss() == 0) continue;
    key _TrkandDetID(thePoint->GetTrackID(), thePoint->GetDetectorID());
    std::map<key, Int_t>::iterator it = fMcTrackingInfo.find(_TrkandDetID);
    if( fMcTrackingInfo.end() == it )
    {
      fMcTrackingInfo.insert(std::pair<key,Int_t>(_TrkandDetID, 1));
      TrackID = thePoint->GetTrackID();
      DetectorID = thePoint->GetDetectorID();
      //ParticleType = thePoint->GetParticleType();
      EntrancePosition = thePoint->GetPosIn();
      ExitingPosition = thePoint->GetPosOut();
      Momentum = thePoint->GetMomIn();
      TimeStamp = thePoint->GetTime() + fEventTime;

      TVector3 fNewEntrance(0,0,0);
      TVector3 fNewExit(0,0,0);
      Bool_t fOK1 = fGeoIF->MasterToLocal(DetectorID, EntrancePosition, fNewEntrance);
      Bool_t fOK2 = fGeoIF->MasterToLocal(DetectorID, ExitingPosition, fNewExit);
      if( ! (fOK1 && fOK2) ) {
        cout <<"PndMdtPointsToWaveform::exec_t, invalid detector index #"<<DetectorID<<endl;
        continue; 
      }
      TVector2 fexit = fNewExit.XYvector();
      TVector2 fentr = fNewEntrance.XYvector();
      TVector2 fvd = fexit - fentr;

      fDis = TMath::Abs(fexit.Mod()*fentr.Mod()*TMath::Sin(fexit.DeltaPhi(fentr)))/fvd.Mod();
      if(fVerbose > 1){
        cout<<"================================================"<<endl;
        cout<<"Mod = "<<thePoint->GetModule()<<endl;
        cout<<"Oct = "<<thePoint->GetSector()<<endl
          <<"Entrance "<<EntrancePosition<<endl
          <<"Exiting "<<ExitingPosition<<endl
          <<"new Entrance "<<fNewEntrance<<endl
          <<"new Exiting "<<fNewExit<<endl;
      }
      PndMCTrack *mcTrk = (PndMCTrack*)fMcTrackArray->At(thePoint->GetTrackID());
      ParticleType= PdgToIndex(mcTrk->GetPdgCode());
      fPid = ParticleType;

      fParamDigiModel->SetParams(ParticleType
          , Momentum
          , fNewEntrance 
          , fNewExit);
      fParamDigiModel->Compute();

      PndMdtWaveform* aWf = new PndMdtWaveform(TrackID, DetectorID, TimeStamp, kTRUE);//wire signal
      ++nWaveformProduced;
      aWf->SetSignal(fParamDigiModel->GetWireSignal());
      fDataBuffer->FillNewData(aWf, aWf->GetTimeStamp(), aWf->GetActiveTime());
      //aWf->Print();
      delete aWf;

      const std::map<Int_t, std::vector<Double_t> >&  fStripSignals = fParamDigiModel->GetStripSignals();
      std::map<Int_t, std::vector<Double_t> >::const_iterator mit = fStripSignals.begin();
      std::map<Int_t, std::vector<Double_t> >::const_iterator mend = fStripSignals.end();
      Int_t detID;
      for(; mit != mend; ++mit, ++nWaveformProduced){
        //re-make detector id;
        fGeoIF->MapWireToStrip(DetectorID, EntrancePosition, detID);//waiting for updates for new geometry versions
        Int_t lid = mit->first;
        if(lid%2 == 0) detID += lid/2;
        else detID -= lid/2;
        PndMdtWaveform* aSf = new PndMdtWaveform(TrackID, detID, TimeStamp, kFALSE);//strip signal
        ++nWaveformProduced;
        aSf->SetSignal(mit->second);
        fDataBuffer->FillNewData(aSf, aSf->GetTimeStamp(), aSf->GetActiveTime());
        delete aSf;
      }
    }else{
      ++ it->second;
    }
  }

  if(fVerbose > 1){
    std::map<key, Int_t>::iterator it = fMcTrackingInfo.begin();
    std::map<key, Int_t>::iterator end = fMcTrackingInfo.end();
    for(; it != end; ++it){
      cout<<"Track #"<<it->first.TrkID<<" contributes "<<it->second<<" points in tube "<<it->first.DetID<<endl;
    }
  }

  if (fVerbose>2){
    timer.Stop();
    Double_t rtime = timer.RealTime();
    Double_t ctime = timer.CpuTime();
    cout << "PndMdtPointsToWaveform, Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
  }
}

void PndMdtPointsToWaveform::FinishTask

(

)

Definition at line 452 of file PndMdtPointsToWaveform.cxx.

References fFile, fVerbose, HowManyPoint, nWaveformProduced, and tTree.

{
  if (fVerbose>0)
    {
      std::cout<<"==================================================="<<std::endl;
      std::cout<<"PndMdtPointsToWaveform::FinishTask"<<std::endl;
      std::cout<<"***************************************************"<<std::endl;
      std::cout<<"Read points# "<<HowManyPoint<<std::endl;
      std::cout<<"Produce waveforms# "<<nWaveformProduced<<std::endl;
      std::cout<<"***************************************************"<<std::endl;
    }
  fFile->cd();
  tTree->Write();
  fFile->Close();
}

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

Referenced by PndLmdPixelHitProducerFast::GetPersistance(), PndMdtDigitization::Init(), PndMdtHitProducerIdeal::Init(), PndMdtClusterTask::Init(), PndFtsHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndRichHitProducer::Init(), PndSttHitProducerRealFast::Init(), PndSttHelixHitProducer::Init(), PndDiscTaskPID::Init(), PndIdealTrackFinder::Init(), PndSttMvdGemTracking::Init(), PndMdtTrkProducer::Init(), PndFtsHitProducerRealFull::Init(), PndLmdPixelClusterTask::Init(), PndSttHitProducerRealFull::Init(), PndLmdStripClusterTask::Init(), PndEmcApdHitProducer::Init(), PndMissingPzCleanerTask::Init(), PndEmcMakeRecoHit::Init(), PndEmcMakeClusterOnline::Init(), PndTrackSmearTask::Init(), PndEmcFWEndcapTimebasedWaveforms::Init(), PndSttHitProducerIdeal::Init(), PndEmcFWEndcapDigi::Init(), PndFtsHitProducerIdeal::Init(), PndEmcMakeCluster::Init(), Init(), PndDiscTaskDigitization::Init(), PndEmcMakeDigi::Init(), PndSdsTimeWalkCorrTask::Init(), PndLmdPixelHitProducerFast::Init(), PndDrcHitFinder::Init(), PndRichHitFinder::Init(), PndEmcMakeCorr::Init(), PndFtofHitProducerIdeal::Init(), PndEmcHitsToWaveform::Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), PndRecoMultiKalmanTask2::Init(), PndEmcHitProducer::Init(), PndDrcHitProducerReal::Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndRecoKalmanTask2::Init(), PndDrcTimeDigiTask::Init(), PndEmcExpClusterSplitter::Init(), PndFtsHoughTrackerTask::Init(), PndSdsNoiseProducer::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsIdealRecoTask::Init(), PndSdsHybridHitProducer::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndGemDigitize::Init(), PndSdsStripHitProducer::Init(), PndGemFindHits::Init(), PndSdsPixelClusterTask::Init(), PndSdsStripClusterTask::Init(), PndMvdGemTrackFinderOnHits::Init(), PndBarrelTrackFinder::Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

{ return fPersistency; }

InitStatus PndMdtPointsToWaveform::Init ( )

virtual

Virtual method Init

Definition at line 55 of file PndMdtPointsToWaveform.cxx.

References PndMdtIGeometry::AddSensor(), Bool_t, fDataBuffer, fDis, fEvtT, fFile, fGeoIF, fLength, fMcTrackArray, fMod, fParamDigiModel, fPid, fPointArray, fStripT, fTimeOrderedWaveform, fVerbose, fWaveformArray, fWirpT, PndPersistencyTask::GetPersistency(), HowManyPoint, PndMdtIGeometry::Init(), PndMdtParamDigi::Init(), PndMdtIGeometry::Instance(), nWaveformProduced, PndMdtIGeometry::Print(), PndMdtParamDigi::SetNoiseWidth(), PndMdtParamDigi::SetOptimization(), PndMdtIGeometry::SetVerbose(), tTree, PndMdtParamDigi::UseDetailedSim(), PndMdtParamDigi::UseGaussianAmp(), PndMdtParamDigi::UseNoise(), and PndMdtParamDigi::UsePlot().

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

  // Get input array
  fMcTrackArray= (TClonesArray*) ioman->GetObject("MCTrack");
  if ( ! fMcTrackArray) {
    cout << "-W- PndMdtPointsToWaveform::Init: "
      << "No MCTrack array!" << endl;
    return kERROR;
  }

  // Get input array
  fPointArray = (TClonesArray*) ioman->GetObject("MdtPoint");
  if ( ! fPointArray) {
    cout << "-W- PndMdtPointsToWaveform::Init: "
      << "No MdtHit array!" << endl;
    return kERROR;
  }
  if(fTimeOrderedWaveform){
    fDataBuffer = new PndMdtWaveformWriteoutBuffer("MdtWaveform", "Mdt", GetPersistency());
    fDataBuffer->ActivateBuffering(kTRUE);
    fDataBuffer->SetVerbose(fVerbose);
    ioman->RegisterWriteoutBuffer("MdtWaveform", fDataBuffer);
  }else{
    fWaveformArray = ioman->Register("MdtWaveform", "PndMdtWaveform", "Mdt",GetPersistency());
  }


  //fDigiPar->printParams();

  nWaveformProduced = 0;
  HowManyPoint = 0;

  fParamDigiModel = new PndMdtParamDigi();
  fParamDigiModel->UseNoise(kTRUE);
  fParamDigiModel->SetNoiseWidth(0.05, 0.02);
  fParamDigiModel->UseDetailedSim(kFALSE);
  fParamDigiModel->UsePlot(kFALSE);
  fParamDigiModel->UseGaussianAmp(kFALSE);
  fParamDigiModel->SetOptimization(10);
  //fParamDigiModel->SetVerbose(2);
  fParamDigiModel->Init();

  fGeoIF = PndMdtIGeometry::Instance();
  fGeoIF->SetVerbose(0);
  fGeoIF->AddSensor("MDT");
  fGeoIF->AddSensor("GasCell");
  Bool_t fGoodGeo = fGeoIF->Init();
  if(fGoodGeo) fGeoIF->Print();


  fFile = new TFile("PndMdtPointsToWaveform.root", "RECREATE");
  tTree = new TTree("pt", "time-stamp diff");
  tTree->Branch("wt", &fWirpT, "wt/D");
  tTree->Branch("st", &fStripT, "st/D");
  tTree->Branch("et", &fEvtT, "et/D");
  tTree->Branch("len", &fLength, "len/D");
  tTree->Branch("dis", &fDis, "dis/D");
  tTree->Branch("mod", &fMod, "mod/I");
  tTree->Branch("pid", &fPid, "pid/I");
  //tTree->Branch("xs", &fxs, "xs/D");
  //tTree->Branch("ys", &fys, "ys/D");
  //tTree->Branch("zs", &fzs, "zs/D");
  //tTree->Branch("xe", &fxe, "xe/D");
  //tTree->Branch("ye", &fye, "ye/D");
  //tTree->Branch("ze", &fze, "ze/D");

  cout << "-I- PndMdtPointsToWaveform: Intialization "<<(fGoodGeo ? "successful." : "failed.") << endl;

  return fGoodGeo ? kSUCCESS : kERROR;
}

PndMdtPointsToWaveform& PndMdtPointsToWaveform::operator= ( const PndMdtPointsToWaveform &  )

private

Int_t PndMdtPointsToWaveform::PdgToIndex ( Int_t  pdg )

private

Definition at line 493 of file PndMdtPointsToWaveform.cxx.

Referenced by exec_e(), and exec_t().

{
  Int_t abspdg = abs(pdg);
  if(abspdg == 11) return 0;
  else if(abspdg == 211) return 2;
  else if(abspdg == 13) return 1;
  else if(abspdg == 2212) return 4;
  else if(abspdg == 321) return 3;
  else return 1;
}

void PndMdtPointsToWaveform::RunTimeBased ( )

inline

Definition at line 52 of file PndMdtPointsToWaveform.h.

References fTimeOrderedWaveform.

Referenced by digi_complete_tb().

{ fTimeOrderedWaveform = kTRUE;}

void PndMdtPointsToWaveform::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 431 of file PndMdtPointsToWaveform.cxx.

References run.

                                              {

  // Get run and runtime database
  FairRun* run = FairRun::Instance();
  if ( ! run ) Fatal("SetParContainers", "No analysis run");

  FairRuntimeDb* db = run->GetRuntimeDb();
  if ( ! db ) Fatal("SetParContainers", "No runtime database");

  // Get Mdt geometry parameter container
  //  fGeoPar = (PndMdtGeoPar*) db->getContainer("PndMdtGeoPar");
  // Get Mdt digitisation parameter container
  //  fDigiPar = (PndMdtDigiPar*) db->getContainer("PndMdtDigiPar");

}

void PndPersistencyTask::SetPersistency ( Bool_t  val = kTRUE )

inlineinherited

Definition at line 31 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency, and val.

Referenced by barrelTrackFinder(), digi_complete(), digi_complete_newSTT(), digiOnly_complete(), PndBarrelTrackFinder::PndBarrelTrackFinder(), PndCATracking::PndCATracking(), PndDrcHitFinder::PndDrcHitFinder(), PndEmc2DLocMaxFinder::PndEmc2DLocMaxFinder(), PndEmcExpClusterSplitter::PndEmcExpClusterSplitter(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcFWEndcapDigi::PndEmcFWEndcapDigi(), PndEmcFWEndcapTimebasedWaveforms::PndEmcFWEndcapTimebasedWaveforms(), PndEmcHitProducer::PndEmcHitProducer(), PndEmcHitsToWaveform::PndEmcHitsToWaveform(), PndEmcMakeBump::PndEmcMakeBump(), PndEmcMakeCluster::PndEmcMakeCluster(), PndEmcMakeClusterOnline::PndEmcMakeClusterOnline(), PndEmcMakeDigi::PndEmcMakeDigi(), PndEmcMakeRecoHit::PndEmcMakeRecoHit(), PndEmcMultiWaveformToCalibratedDigi::PndEmcMultiWaveformToCalibratedDigi(), PndEmcPhiBumpSplitter::PndEmcPhiBumpSplitter(), PndEmcTmpWaveformToDigi::PndEmcTmpWaveformToDigi(), PndEmcWaveformToCalibratedDigi::PndEmcWaveformToCalibratedDigi(), PndEmcWaveformToDigi::PndEmcWaveformToDigi(), PndFtofHitProducerIdeal::PndFtofHitProducerIdeal(), PndFtsCATracking::PndFtsCATracking(), PndFtsHitProducerIdeal::PndFtsHitProducerIdeal(), PndFtsHitProducerRealFast::PndFtsHitProducerRealFast(), PndFtsHitProducerRealFull::PndFtsHitProducerRealFull(), PndFtsHoughTrackerTask::PndFtsHoughTrackerTask(), PndGemDigitize::PndGemDigitize(), PndGemFindHits::PndGemFindHits(), PndIdealTrackFinder::PndIdealTrackFinder(), PndLmdPixelClusterTask::PndLmdPixelClusterTask(), PndLmdPixelHitProducerFast::PndLmdPixelHitProducerFast(), PndMdtClusterTask::PndMdtClusterTask(), PndMdtDigitization::PndMdtDigitization(), PndMdtHitProducerIdeal::PndMdtHitProducerIdeal(), PndMdtPointsToWaveform(), PndMdtTrkProducer::PndMdtTrkProducer(), PndMissingPzCleanerTask::PndMissingPzCleanerTask(), PndMvdGemTrackFinderOnHits::PndMvdGemTrackFinderOnHits(), PndMvdHitProducerIdeal::PndMvdHitProducerIdeal(), PndMvdPixelClusterTask::PndMvdPixelClusterTask(), PndMvdTimeWalkCorrTask::PndMvdTimeWalkCorrTask(), PndMvdToPix4ClusterTask::PndMvdToPix4ClusterTask(), PndRecoKalmanTask::PndRecoKalmanTask(), PndRecoKalmanTask2::PndRecoKalmanTask2(), PndRecoMultiKalmanTask::PndRecoMultiKalmanTask(), PndRecoMultiKalmanTask2::PndRecoMultiKalmanTask2(), PndRichHitFinder::PndRichHitFinder(), PndRichHitProducer::PndRichHitProducer(), PndSciTDigiTask::PndSciTDigiTask(), PndSciTHitProducerIdeal::PndSciTHitProducerIdeal(), PndSdsHitProducerIdeal::PndSdsHitProducerIdeal(), PndSdsHybridHitProducer::PndSdsHybridHitProducer(), PndSdsIdealClusterTask::PndSdsIdealClusterTask(), PndSdsIdealRecoTask::PndSdsIdealRecoTask(), PndSdsNoiseProducer::PndSdsNoiseProducer(), PndSdsPixelClusterTask::PndSdsPixelClusterTask(), PndSdsStripClusterTask::PndSdsStripClusterTask(), PndSdsStripHitProducer::PndSdsStripHitProducer(), PndSdsTimeWalkCorrTask::PndSdsTimeWalkCorrTask(), PndSttFindTracks::PndSttFindTracks(), PndSttHelixHitProducer::PndSttHelixHitProducer(), PndSttHitProducerIdeal::PndSttHitProducerIdeal(), PndSttHitProducerRealFast::PndSttHitProducerRealFast(), PndSttHitProducerRealFull::PndSttHitProducerRealFull(), PndSttMatchTracks::PndSttMatchTracks(), PndSttMvdGemTracking::PndSttMvdGemTracking(), PndTrackSmearTask::PndTrackSmearTask(), PndTrkTracking2::PndTrkTracking2(), reco(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), recoideal_complete(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndLmdPixelHitProducerFast::SetPersistance(), PndSdsHitProducerIdeal::SetPersistance(), PndSttMvdGemTracking::SetPersistenc(), PndMdtClusterTask::SetPersistence(), PndSttHelixHitProducer::SetPersistence(), PndMissingPzCleanerTask::SetPersistence(), PndFtsHitProducerRealFast::SetPersistence(), PndFtsHitProducerRealFull::SetPersistence(), PndSttHitProducerIdeal::SetPersistence(), PndSttHitProducerRealFull::SetPersistence(), PndSttHitProducerRealFast::SetPersistence(), PndFtsHitProducerIdeal::SetPersistence(), PndTrackSmearTask::SetPersistence(), PndSciTHitProducerIdeal::SetPersistence(), PndIdealTrackFinder::SetPersistence(), PndSttMatchTracks::SetPersistence(), PndSttFindTracks::SetPersistence(), PndFtsHoughTrackerTask::SetPersistence(), PndTrkTracking2::SetPersistence(), PndEmcMakeRecoHit::SetStorageOfData(), PndEmcMakeClusterOnline::SetStorageOfData(), PndEmcFWEndcapDigi::SetStorageOfData(), PndEmcFWEndcapTimebasedWaveforms::SetStorageOfData(), PndEmcMakeDigi::SetStorageOfData(), SetStorageOfData(), PndEmc2DLocMaxFinder::SetStorageOfData(), PndEmcMakeCluster::SetStorageOfData(), PndEmcHitsToWaveform::SetStorageOfData(), PndEmcMakeBump::SetStorageOfData(), PndEmcTmpWaveformToDigi::SetStorageOfData(), PndEmcWaveformToDigi::SetStorageOfData(), PndEmcWaveformToCalibratedDigi::SetStorageOfData(), PndEmcMultiWaveformToCalibratedDigi::SetStorageOfData(), PndEmcExpClusterSplitter::SetStorageOfData(), PndEmcPhiBumpSplitter::SetStorageOfData(), standard_tracking(), and PndEmcFullDigiTask::StoreDigi().

{ fPersistency = val; }

void PndMdtPointsToWaveform::SetStorageOfData

(

Bool_t

val

)

Definition at line 447 of file PndMdtPointsToWaveform.cxx.

References PndPersistencyTask::SetPersistency().

{
  SetPersistency(val);
  return;
}

Member Data Documentation

PndMdtWaveformWriteoutBuffer* PndMdtPointsToWaveform::fDataBuffer

private

Definition at line 66 of file PndMdtPointsToWaveform.h.

Referenced by exec_t(), and Init().

Double_t PndMdtPointsToWaveform::fDis

private

Definition at line 112 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), and Init().

Double_t PndMdtPointsToWaveform::fEvtT

private

Definition at line 110 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), and Init().

TFile* PndMdtPointsToWaveform::fFile

private

Definition at line 106 of file PndMdtPointsToWaveform.h.

Referenced by FinishTask(), and Init().

PndMdtIGeometry* PndMdtPointsToWaveform::fGeoIF

private

Definition at line 69 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), and Init().

Double_t PndMdtPointsToWaveform::fLength

private

Definition at line 111 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), and Init().

TClonesArray* PndMdtPointsToWaveform::fMcTrackArray

private

Definition at line 56 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), and Init().

Int_t PndMdtPointsToWaveform::fMod

private

Definition at line 113 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), and Init().

PndMdtParamDigi* PndMdtPointsToWaveform::fParamDigiModel

private

Definition at line 68 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), and Init().

Int_t PndMdtPointsToWaveform::fPid

private

Definition at line 114 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), and Init().

TClonesArray* PndMdtPointsToWaveform::fPointArray

private

Input array of PndMdtHits

Definition at line 59 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), and Init().

Double_t PndMdtPointsToWaveform::fStripT

private

Definition at line 109 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), and Init().

Bool_t PndMdtPointsToWaveform::fTimeOrderedWaveform

private

Definition at line 65 of file PndMdtPointsToWaveform.h.

Referenced by Exec(), Init(), and RunTimeBased().

Int_t PndMdtPointsToWaveform::fVerbose

private

Verbosity level

Definition at line 80 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), FinishTask(), and Init().

TClonesArray* PndMdtPointsToWaveform::fWaveformArray

private

Definition at line 62 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), and Init().

Double_t PndMdtPointsToWaveform::fWirpT

private

Definition at line 108 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), and Init().

Int_t PndMdtPointsToWaveform::HowManyPoint

private

Definition at line 82 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), FinishTask(), and Init().

Int_t PndMdtPointsToWaveform::nWaveformProduced

private

Definition at line 83 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), exec_t(), FinishTask(), and Init().

TTree* PndMdtPointsToWaveform::tTree

private

Definition at line 107 of file PndMdtPointsToWaveform.h.

Referenced by exec_e(), FinishTask(), and Init().

---

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

• PndMdtPointsToWaveform.h
• PndMdtPointsToWaveform.cxx