PndGemDigitize Class Reference

#include <PndGemDigitize.h>

Inheritance diagram for PndGemDigitize:

Public Member Functions
	PndGemDigitize ()
	PndGemDigitize (Int_t iVerbose)
	PndGemDigitize (const char *name, Int_t iVerbose)
virtual	~PndGemDigitize ()
virtual void	Exec (Option_t *opt)
void	SaveOutsideHits (Bool_t bt=kTRUE)
void	SetRealisticResponse (Bool_t bt=kTRUE)
void	RunTimeBased (Bool_t bt=kTRUE)
Bool_t	GetTimeBased () const
void	SetShowerSigma (Double_t dd)
void	SetPersistency (Bool_t val=kTRUE)
Bool_t	GetPersistency ()

Private Member Functions
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
void	Reset ()
virtual void	Finish ()
void	DigitizeEvent ()
void	DigitizeRealisticEvent ()
void	SimulateRectangularResponse (Int_t sensorDetId, Int_t side, Double_t channelInd, Double_t stripWidth, Double_t showerSigma, Double_t showerStrength, Int_t iPoint)
void	SimulateGaussianResponse (PndGemSensor *sensor, Int_t side, PndGemMCPoint *gemPoint, Double_t showerSigma, Double_t showerStrength, Int_t iPoint)
void	ActivateChannel (Int_t sensorDetId, Int_t sensorSide, Int_t channelNumber, Double_t signalHeight, Double_t signalTime, Int_t iPoint)
void	PrintDigis ()
	ClassDef (PndGemDigitize, 1)

Private Attributes
PndGemDigiPar *	fDigiPar
TClonesArray *	fPoints
TClonesArray *	fDigis
TClonesArray *	fDigiMatches
Double_t	fSSigma
Bool_t	fSaveOutsideHits
TClonesArray *	fHitOutsideArray
Bool_t	fRealisticResponse
Int_t	fTNofEvents
Int_t	fTNofPoints
Int_t	fTNofDigis
TRandom2 *	fRand
PndGemDigiWriteoutBuffer *	fDataBuffer
Bool_t	fTimeOrderedDigi
Int_t	fNPoints
Int_t	fNFailed
Int_t	fNOutside
Int_t	fNMulti
Int_t	fNDigis
TStopwatch	fTimer
std::map< std::pair< Int_t, Int_t >, Int_t >	fChannelMap

Detailed Description

Definition at line 37 of file PndGemDigitize.h.

Constructor & Destructor Documentation

PndGemDigitize::PndGemDigitize

(

)

Default constructor

Definition at line 50 of file PndGemDigitize.cxx.

References Reset(), and PndPersistencyTask::SetPersistency().

  : PndPersistencyTask("GEM Digitizer", 0),
    fDigiPar(NULL),
    fPoints(NULL),
    fDigis(NULL),
    fDigiMatches(NULL),
    fSSigma(0.03),
    fSaveOutsideHits(kFALSE),
    fHitOutsideArray(NULL),
    fRealisticResponse(kFALSE),
    fTNofEvents(0),
    fTNofPoints(0),
    fTNofDigis(0),
    fRand(new TRandom2()),
    fDataBuffer(0),
    fTimeOrderedDigi(kFALSE)
{
  SetPersistency(kTRUE);
  Reset();
}

PndGemDigitize::PndGemDigitize

(

Int_t

iVerbose

)

Standard constructor

Definition at line 75 of file PndGemDigitize.cxx.

References Reset(), and PndPersistencyTask::SetPersistency().

  : PndPersistencyTask("GEM Digitizer", iVerbose),
    fDigiPar(NULL),
    fPoints(NULL),
    fDigis(NULL),
    fDigiMatches(NULL),
    fSSigma(0.03),
    fSaveOutsideHits(kFALSE),
    fHitOutsideArray(NULL),
    fRealisticResponse(kFALSE),
    fTNofEvents(0),
    fTNofPoints(0),
    fTNofDigis(0),
    fRand(new TRandom2()),
    fDataBuffer(0),
    fTimeOrderedDigi(kFALSE)
{
  SetPersistency(kTRUE);
  Reset();
}

PndGemDigitize::PndGemDigitize	(	const char *	name,
		Int_t	iVerbose
	)

Constructor with name

Definition at line 100 of file PndGemDigitize.cxx.

References Reset(), and PndPersistencyTask::SetPersistency().

  : PndPersistencyTask(name, iVerbose),
    fDigiPar(NULL),
    fPoints(NULL),
    fDigis(NULL),
    fDigiMatches(NULL),
    fSSigma(0.03),
    fSaveOutsideHits(kFALSE),
    fHitOutsideArray(NULL),
    fRealisticResponse(kFALSE),
    fTNofEvents(0),
    fTNofPoints(0),
    fTNofDigis(0),
    fRand(new TRandom2()),
    fDataBuffer(0),
    fTimeOrderedDigi(kFALSE)
{
        SetPersistency(kTRUE);
        Reset();
}

PndGemDigitize::~PndGemDigitize ( )

virtual

Destructor

Definition at line 125 of file PndGemDigitize.cxx.

References fDigiMatches, fDigiPar, fDigis, and Reset().

                                { 
  Reset();
  if ( fDigiPar)   delete fDigiPar;
  if ( fDigis ) {
    fDigis->Delete();
    delete fDigis;
  }
  if ( fDigiMatches ) {
    fDigiMatches->Delete();
    delete fDigiMatches;
  }

}

Member Function Documentation

void PndGemDigitize::ActivateChannel ( Int_t  sensorDetId, Int_t  sensorSide, Int_t  channelNumber, Double_t  signalHeight, Double_t  signalTime, Int_t  iPoint  )

private

Activate channel

Definition at line 504 of file PndGemDigitize.cxx.

References a, PndGemDigi::AddCharge(), PndGemDigi::AddIndex(), fChannelMap, fDigis, fNDigis, fTNofDigis, and kGemDigi.

Referenced by DigitizeEvent(), SimulateGaussianResponse(), and SimulateRectangularResponse().

                                                                                                                                                       {
  Int_t digiDetId = sensorDetId | kGemDigi << 21 | sensorSide << 5;

  pair<Int_t, Int_t> a (digiDetId, channelNumber);
  if ( fChannelMap.find(a) == fChannelMap.end() ) {
    // Channel not yet active, create new digi
    
    //  cout << fNDigis << " : " << sensorDetId << " " << channelNumber << endl;
    
    new ((*fDigis)[fNDigis]) PndGemDigi(digiDetId, channelNumber, iPoint, signalHeight, signalTime);

    fChannelMap[a] = fNDigis;
    fNDigis++;
    fTNofDigis++;
  }
  else {
    Int_t iDigi = fChannelMap[a];
    PndGemDigi* ddigi = dynamic_cast<PndGemDigi*>(fDigis->At(iDigi));

    ddigi->AddCharge(signalHeight);

    //  cout << "another mc point to this digi, last one " << ddigi->GetIndex(ddigi->GetNIndices()-1) << ", now adding " << iPoint << endl;
    ddigi->AddIndex(iPoint);
  }
}

PndGemDigitize::ClassDef ( PndGemDigitize  , 1    )

private

void PndGemDigitize::DigitizeEvent ( )

private

Digitize MC points in one event

Definition at line 186 of file PndGemDigitize.cxx.

References ActivateChannel(), Double_t, fDigiPar, fPoints, fSaveOutsideHits, fTNofPoints, fVerbose, PndGemSensor::GetChannel(), PndGemSensor::GetDetectorId(), PndGemDigiPar::GetNodeName(), PndGemDigiPar::GetSensor(), PndGemMCPoint::GetSensorId(), PndGemDigiPar::GetStationNr(), PndGemMCPoint::GetXOut(), PndGemMCPoint::GetYOut(), gGeoManager, pos, sensor, and TString.

Referenced by Exec().

                                   {
  if ( fVerbose > 0 ) cout << "-I- PndGemDigitize::DigitizeEvent()" << endl;
  
  PndGemSensor* sensor;
  
  Int_t nofHitsOutside = 0;
  Int_t nofPoints = fPoints->GetEntriesFast();

//   cout << "GEM digitize points:" << endl;
//   for ( Int_t iPoint = 0 ; iPoint < nofPoints ; iPoint++ ) {
//     PndGemMCPoint* currentPndGemMCPoint = (PndGemMCPoint*)fPoints->At(iPoint);
//     cout << iPoint << " : " 
//       << currentPndGemMCPoint->GetX() << " "
//       << currentPndGemMCPoint->GetY() << " "
//       << currentPndGemMCPoint->GetZ() << endl;
//   }
//   cout << "GEM digis: " << endl;

  for ( Int_t iPoint = 0 ; iPoint < nofPoints ; iPoint++ ) {
    PndGemMCPoint* currentPndGemMCPoint = (PndGemMCPoint*)fPoints->At(iPoint);

    Double_t posIn[3] = {0.5*(currentPndGemMCPoint->GetX()+currentPndGemMCPoint->GetXOut()),
                         0.5*(currentPndGemMCPoint->GetY()+currentPndGemMCPoint->GetYOut()),
                         currentPndGemMCPoint->GetZ()};
    
    Int_t sensorId = currentPndGemMCPoint->GetSensorId();

    TString nodeName = fDigiPar->GetNodeName(sensorId);

    gGeoManager->cd(nodeName.Data());
    TGeoNode* curNode = gGeoManager->GetCurrentNode();

//     cout << "id " << sensorId << " -> " 
//       << fDigiPar->GetStationNr(sensorId) << "."
//       << fDigiPar->GetSensorNr(sensorId) << "."
//       << fDigiPar->GetSegmentNr(sensorId) << " > " 
//       << nodeName.Data() << endl;

    sensor = (PndGemSensor*)fDigiPar->GetSensor(sensorId);
    if ( !sensor ) {
      cout << " -E- " << GetName() << ":Exec() There is no sensor: \"" 
           << nodeName.Data() << "\"." << endl;
      continue;
    }
    Double_t locPosIn[4];

    fTNofPoints++;
    
    curNode->MasterToLocal(posIn,locPosIn);

    //    cout << "position:  " << posIn[0] << " " << posIn[1] << " " << posIn[2] << endl;
    //   << " transf to " << locPosIn[0] << " " << locPosIn[1] << " " << locPosIn[2] << endl;
    
    //    if ( sensor->GetType()!=1 ) { locPosIn[3] = locPosIn[2]; locPosIn[2] = locPosIn[1]; locPosIn[1] = locPosIn[3]; locPosIn[0] = -locPosIn[0]; }
    
    Int_t sensorDetId = sensor->GetDetectorId();

    Int_t channelNumber = sensor->GetChannel(locPosIn[0],locPosIn[1],0);
    if ( fDigiPar->GetStationNr(sensorId) == 1 ) {
      //      cout << sensor->GetType() << ": for point " << locPosIn[0] << " " << locPosIn[1] << " ---> channel " << channelNumber << endl;
    }
    if ( channelNumber == -1 ) {
      if ( fSaveOutsideHits ) {
        TVector3 pos;
        currentPndGemMCPoint->Position(pos);
        TVector3 dposLocal(0.,0.,0.);


        //Int_t hitDetId = sensorDetId | kGemHit << 21; //[R.K. 01/2017] unused variable?
      
        new ((*fHitOutsideArray)[nofHitsOutside]) PndGemHit(sensorDetId,
                                                            pos,dposLocal,iPoint,currentPndGemMCPoint->GetEnergyLoss(),1);
      }
      nofHitsOutside++;
    }
    else {
      ActivateChannel(sensorDetId,0,channelNumber,1,0.,iPoint);
    }

    channelNumber = sensor->GetChannel(locPosIn[0],locPosIn[1],1);
    if ( channelNumber == -1 ) continue;
    
    ActivateChannel(sensorDetId,1,channelNumber,1.,0.,iPoint);
    
  }
}

void PndGemDigitize::DigitizeRealisticEvent ( )

private

Digitize MC points in one event

Definition at line 275 of file PndGemDigitize.cxx.

References Double_t, fDigiPar, fPoints, fSSigma, fTNofPoints, fVerbose, PndGemDigiPar::GetNChannels(), PndGemDigiPar::GetNodeName(), PndGemDigiPar::GetNSensors(), PndGemDigiPar::GetNStations(), PndGemDigiPar::GetSensor(), PndGemMCPoint::GetSensorId(), gGeoManager, sensor, SimulateGaussianResponse(), and TString.

Referenced by Exec().

                                            {
  if ( fVerbose > 0 ) cout << "-I- PndGemDigitize::DigitizeRealisticEvent()" << endl; 
  //  cout << "DRE " << fTNofEvents << endl;
  
  //Double_t showerSigma =   0.1; // radius of smearing, in cm  
  Double_t showerSigma =   0.03; // radius of smearing, in cm (based on Andrii's sim)
  showerSigma = fSSigma;
  //Double_t totalSignal = 100.;  // total signal strength //[R.K. 01/2017] unused variable?
  Double_t sigMult = 1.e6;
  
  if ( fVerbose > 1 ) {
    cout << " showerSigma=" << showerSigma  
         << " sigMult=" << sigMult
         << endl;
  }
  
  PndGemSensor* sensor;
  
  //Int_t nofHitsOutside = 0; //[R.K. 01/2017] unused variable?
  Int_t nofPoints = fPoints->GetEntriesFast();

  // if ( fVerbose > 0 ) {
  //   cout << "GEM digitize points:" << endl;
  //   for ( Int_t iPoint = 0 ; iPoint < nofPoints ; iPoint++ ) {
  //     PndGemMCPoint* currentPndGemMCPoint = (PndGemMCPoint*)fPoints->At(iPoint);    }
  //   cout << "GEM digis: " << endl;
  // }

  if ( fVerbose ) {
    cout << "GemDigiPar has "
         << fDigiPar->GetNStations() << " stations, "
         << fDigiPar->GetNSensors() << " sensors, "
         << fDigiPar->GetNChannels() << " channels. "
         << endl;
  }

  for ( Int_t iPoint = 0 ; iPoint < nofPoints ; iPoint++ ) {
    PndGemMCPoint* currentPndGemMCPoint = (PndGemMCPoint*)fPoints->At(iPoint);

    Double_t posIn[3] = {currentPndGemMCPoint->GetX(),
                         currentPndGemMCPoint->GetY(),
                         currentPndGemMCPoint->GetZ()};
    
    Int_t sensorId = currentPndGemMCPoint->GetSensorId();

    if ( fVerbose > 0 ) {
      cout << iPoint << " : " 
           << currentPndGemMCPoint->GetX() << " "
           << currentPndGemMCPoint->GetY() << " "
           << currentPndGemMCPoint->GetZ() << " in sensor "
           << sensorId << endl;
    }

    //Double_t stripWidth; //[R.K. 01/2017] unused variable?

    TString nodeName = fDigiPar->GetNodeName(sensorId);

    if ( fVerbose > 0 ) {
      cout << "got nodeName = \"" << nodeName.Data() << "\"" << endl;
    }

    gGeoManager->cd(nodeName.Data());
    TGeoNode* curNode = gGeoManager->GetCurrentNode();

//     cout << "id " << sensorId << " -> " 
//       << fDigiPar->GetStationNr(sensorId) << "."
//       << fDigiPar->GetSensorNr(sensorId) << "."
//       << fDigiPar->GetSegmentNr(sensorId) << " > " 
//       << nodeName.Data() << endl;
    sensor = (PndGemSensor*)fDigiPar->GetSensor(sensorId);

    if ( !sensor ) {
      cout << " -E- " << GetName() << ":Exec() There is no sensor: \"" 
           << nodeName.Data() << "\"." << endl;
      continue;
    }
    Double_t locPosIn[4];

    fTNofPoints++;
    
    curNode->MasterToLocal(posIn,locPosIn);

    if ( fVerbose ) {
      cout << "position:  " << posIn[0] << " " << posIn[1] << " " << posIn[2] << endl;
      cout << " transf to " << locPosIn[0] << " " << locPosIn[1] << " " << locPosIn[2] << endl;
    }

    //    if ( sensor->GetType()!=1 ) { locPosIn[3] = locPosIn[2]; locPosIn[2] = locPosIn[1]; locPosIn[1] = locPosIn[3]; locPosIn[0] = -locPosIn[0]; }
    
    //Int_t sensorDetId = sensor->GetDetectorId(); //[R.K. 01/2017] unused variable?
    //Double_t rectSig  = 100.; //[R.K. 01/2017] unused variable?
    
    //Double_t channelInd = sensor->GetChannel(locPosIn[0],locPosIn[1],0,stripWidth);
    //Int_t channelNumber = (Int_t)channelInd;

    //    cout << "channelNumber = " << channelNumber << " / channelInd = " << channelInd << " / stripWidth = " << stripWidth << endl;

    //    SimulateRectangularResponse(sensorDetId,0,channelInd,stripWidth,showerSigma,rectSig,iPoint);
    //    SimulateGaussianResponse(sensor,0,currentPndGemMCPoint,showerSigma,rectSig,iPoint);    
    SimulateGaussianResponse(sensor,0,currentPndGemMCPoint,showerSigma,sigMult*currentPndGemMCPoint->GetEnergyLoss(),iPoint);

    //channelInd = sensor->GetChannel(locPosIn[0],locPosIn[1],1,stripWidth);
    //channelNumber = (Int_t)channelInd;

    //    cout << "channelNumber = " << channelNumber << " / channelInd = " << channelInd << " / stripWidth = " << stripWidth << endl;

    //    SimulateRectangularResponse(sensorDetId,1,channelInd,stripWidth,showerSigma,rectSig,iPoint);
    //SimulateGaussianResponse(sensor,1,currentPndGemMCPoint,showerSigma,rectSig,iPoint);
    SimulateGaussianResponse(sensor,1,currentPndGemMCPoint,showerSigma,sigMult*currentPndGemMCPoint->GetEnergyLoss(),iPoint);
    
  }
}

void PndGemDigitize::Exec ( Option_t *  opt )

virtual

Execution

Definition at line 141 of file PndGemDigitize.cxx.

References DigitizeEvent(), DigitizeRealisticEvent(), Double_t, fDataBuffer, fDigis, fHitOutsideArray, fPoints, fRealisticResponse, fSaveOutsideHits, fTNofEvents, fVerbose, and Reset().

                                   {

  if ( fSaveOutsideHits ) {
    if ( ! fHitOutsideArray ) Fatal("Exec", "No fHitOutsideArray");
    fHitOutsideArray->Delete();
  }

  Reset();

  if ( fVerbose )
    cout << "EVENT " << fTNofEvents << endl;

  fTNofEvents++;

  if ( fRealisticResponse ) 
    DigitizeRealisticEvent();
  else
    DigitizeEvent();

  PndGemDigi* adigi;
  Double_t EventTime = FairRootManager::Instance()->GetEventTime();
  // cout << "---------------> Event " << fTNofEvents << endl;
  // cout << "--- at  " << EventTime << " ns" << endl;
  // cout << "--- has " << fDigis->GetEntriesFast() << " digis made from " << fPoints->GetEntriesFast() << " points" << endl;
  // cout << "---------------> Put them in output buffer" << endl;
  if ( fVerbose ) {
    cout << "---PndGemDigitize---> Event " << fTNofEvents
         << " at  " << EventTime << " ns"
         << " has " << fDigis->GetEntriesFast() << " digis made from " << fPoints->GetEntriesFast() << " points" << endl;
  }
  // this loop should be done after each point, and not at the end of the event not to miss any late digis...
  for ( Int_t idigi = 0 ; idigi < fDigis->GetEntriesFast() ; idigi++ ) {
    adigi = dynamic_cast<PndGemDigi*>(fDigis->At(idigi));
    adigi->SetTimeStamp(adigi->GetTimeStamp()+EventTime);
    fDataBuffer->FillNewData(adigi,
                             adigi->GetTimeStamp(),
                             adigi->GetTimeStamp()+100.); // 100 ns dead time
    //                               adigi->GetTimeStamp()+EventTime,
    //                               adigi->GetTimeStamp()+EventTime+100.); // 100 ns dead time
  }
  //  cout << "------------------------------------------" << endl;
}

void PndGemDigitize::Finish ( )

privatevirtual

Finish at the end of each event

Definition at line 624 of file PndGemDigitize.cxx.

References Double_t, fDigiPar, fDigis, fTNofDigis, fTNofEvents, fTNofPoints, PndGemDigiPar::GetNChannels(), and PndGemDigiPar::GetNSensors().

                            {
  if ( fDigis ) fDigis->Delete();

  cout << "-------------------- " << fName.Data() << " : Summary ------------------------" << endl;
  cout << " Events:        " << setw(10) << fTNofEvents << endl;
  cout << " MC Points:     " << setw(10) << fTNofPoints << "    ( " << (Double_t)fTNofPoints/((Double_t)fTNofEvents) << " per event )" << endl;
  cout << " Digis:         " << setw(10) << fTNofDigis  << "    ( " << (Double_t)fTNofDigis /((Double_t)fTNofEvents) << " per event )" << endl;
  cout << "                       -->    ( " << (Double_t)fTNofDigis /((Double_t)fTNofEvents)/((Double_t)fDigiPar->GetNSensors()) << " per sensor )" << endl;
  cout << "                       -->    ( " << (Double_t)fTNofDigis /((Double_t)fTNofEvents)/((Double_t)fDigiPar->GetNChannels())*100. << "% occupancy )" << endl;
  cout << "                       -->    ( 2 x " << (Double_t)fTNofDigis /((Double_t)fTNofPoints)/2. << " per point )" << endl;
  cout << "---------------------------------------------------------------------" << endl; 
}

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(), PndRichHitProducer::Init(), PndSttHitProducerRealFast::Init(), PndDiscTaskReconstruction::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(), PndMdtPointsToWaveform::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(), PndEmcHitProducer::Init(), PndRecoMultiKalmanTask2::Init(), PndDrcHitProducerReal::Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndDrcTimeDigiTask::Init(), PndRecoKalmanTask2::Init(), PndEmcExpClusterSplitter::Init(), PndSdsNoiseProducer::Init(), PndFtsHoughTrackerTask::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsHybridHitProducer::Init(), PndSdsIdealRecoTask::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndSdsStripHitProducer::Init(), 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; }

Bool_t PndGemDigitize::GetTimeBased ( ) const

inline

Definition at line 70 of file PndGemDigitize.h.

References fTimeOrderedDigi.

{ return fTimeOrderedDigi; }

InitStatus PndGemDigitize::Init ( )

privatevirtual

Intialisation

Definition at line 566 of file PndGemDigitize.cxx.

References fDataBuffer, fDigiPar, fDigis, fHitOutsideArray, fPoints, fSaveOutsideHits, fTimeOrderedDigi, PndGemDigiPar::GetNStations(), and PndPersistencyTask::GetPersistency().

                                {

  // Get input array 
  FairRootManager* ioman = FairRootManager::Instance();

  if ( ! ioman ) Fatal("Init", "No FairRootManager");
  fPoints = (TClonesArray*) ioman->GetObject("GEMPoint");

  if ( fSaveOutsideHits ) {
    fHitOutsideArray = new TClonesArray("PndGemHit");
    ioman->Register("GEMOutsideHit", "PndGem Hits in inactive region",
                    fHitOutsideArray, GetPersistency());
  }

  // Register output array StsDigi
  fDigis = new TClonesArray("PndGemDigi",10000);
  ioman->Register("GEMDigiNormal", "Digital response in GEM", fDigis, kFALSE);

  // Register output buffer
  fDataBuffer = new PndGemDigiWriteoutBuffer("GEMDigi", "GEM", GetPersistency());
  fDataBuffer = (PndGemDigiWriteoutBuffer*)ioman->RegisterWriteoutBuffer("GEMDigi", fDataBuffer);
  fDataBuffer->ActivateBuffering(fTimeOrderedDigi);

  cout << "-I- " << fName.Data() << "::Init(). There are " << fDigiPar->GetNStations() << " GEM stations." << endl;
  cout << "-I- " << fName.Data() << "::Init(). Initialization succesfull." << endl;
  return kSUCCESS;

}

void PndGemDigitize::PrintDigis ( )

private

Print Digis

Definition at line 533 of file PndGemDigitize.cxx.

References Double_t, fDigis, PndGemDigi::GetChannelNr(), PndGemDigi::GetCharge(), PndGemDigi::GetDetectorId(), PndGemDigi::GetSensorNr(), PndGemDigi::GetSide(), and PndGemDigi::GetStationNr().

                                {
  Int_t nofDigis = fDigis->GetEntriesFast();
  for ( Int_t idigi = 0 ; idigi < nofDigis ; idigi++ ) {
    PndGemDigi* digiToDraw   = (PndGemDigi*)fDigis->At(idigi);
    Int_t digiStationNr      = digiToDraw->GetStationNr();
    Int_t digiSensorNr       = digiToDraw->GetSensorNr();
    Int_t digiSide           = digiToDraw->GetSide();
    Double_t channelNumber   = digiToDraw->GetChannelNr();

    cout << "HAVE digi " << idigi << " at " << digiStationNr << "." << digiSensorNr << "." << digiSide << "." << channelNumber << " ( " << digiToDraw->GetDetectorId() << " ) with signal = " << digiToDraw->GetCharge() << endl;

  }
}

InitStatus PndGemDigitize::ReInit ( )

privatevirtual

Reinitialisation

Definition at line 599 of file PndGemDigitize.cxx.

                                  {

  // Clear digitisation scheme
  //  fDigiScheme->Clear();

  // Build new digitisation scheme
  //  if ( fDigiScheme->Init(fGeoPar, fDigiPar) ) return kSUCCESS;

  return kSUCCESS;

}

void PndGemDigitize::Reset ( )

private

Reset eventwise counters

Definition at line 615 of file PndGemDigitize.cxx.

References fChannelMap, fDigiMatches, fDigis, fNDigis, fNFailed, fNMulti, fNOutside, and fNPoints.

Referenced by Exec(), PndGemDigitize(), and ~PndGemDigitize().

                           {
  fNPoints = fNFailed = fNOutside = fNMulti = fNDigis = 0;
  fChannelMap.clear();
  if ( fDigis ) fDigis->Delete();
  if ( fDigiMatches ) fDigiMatches->Clear();
}

void PndGemDigitize::RunTimeBased ( Bool_t  bt = kTRUE )

inline

Definition at line 68 of file PndGemDigitize.h.

References fTimeOrderedDigi.

Referenced by TB_gem_digi().

{fTimeOrderedDigi   = bt;}

void PndGemDigitize::SaveOutsideHits ( Bool_t  bt = kTRUE )

inline

Switch saving outside hits

Definition at line 63 of file PndGemDigitize.h.

References fSaveOutsideHits.

{fSaveOutsideHits = bt;}

void PndGemDigitize::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 549 of file PndGemDigitize.cxx.

References fDigiPar, and 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 GEM digitisation parameter container
  fDigiPar = (PndGemDigiPar*)(db->getContainer("PndGemDetectors"));
}

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(), PndGemFindHits::PndGemFindHits(), PndIdealTrackFinder::PndIdealTrackFinder(), PndLmdPixelClusterTask::PndLmdPixelClusterTask(), PndLmdPixelHitProducerFast::PndLmdPixelHitProducerFast(), PndMdtClusterTask::PndMdtClusterTask(), PndMdtDigitization::PndMdtDigitization(), PndMdtHitProducerIdeal::PndMdtHitProducerIdeal(), PndMdtPointsToWaveform::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(), PndMdtPointsToWaveform::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 PndGemDigitize::SetRealisticResponse ( Bool_t  bt = kTRUE )

inline

Definition at line 66 of file PndGemDigitize.h.

References fRealisticResponse.

Referenced by TB_gem_digi().

{fRealisticResponse = bt;}

void PndGemDigitize::SetShowerSigma ( Double_t  dd )

inline

Definition at line 72 of file PndGemDigitize.h.

References fSSigma.

{fSSigma = dd;} //sigam for realistic response

void PndGemDigitize::SimulateGaussianResponse ( PndGemSensor *  sensor, Int_t  side, PndGemMCPoint *  gemPoint, Double_t  showerSigma, Double_t  showerStrength, Int_t  iPoint  )

private

Simulated gaussian response

Definition at line 390 of file PndGemDigitize.cxx.

References ActivateChannel(), CAMath::Cos(), Double_t, fDigiPar, fRand, PndGemSensor::GetChannel(), PndGemSensor::GetChannel2(), PndGemSensor::GetDetectorId(), PndGemSensor::GetDistance(), PndGemDigiPar::GetNodeName(), PndGemMCPoint::GetSensorId(), PndGemSensor::GetStripOrientation(), PndGemMCPoint::GetXOut(), PndGemMCPoint::GetYOut(), PndGemMCPoint::GetZOut(), gGeoManager, Pi, CAMath::Sin(), TString, and CAMath::TwoPi().

Referenced by DigitizeRealisticEvent().

                                                                                                                                                                    {
  Double_t posIn [3] = {gemPoint->GetX(),
                        gemPoint->GetY(),
                        gemPoint->GetZ()};
  Double_t posOut[3] = {gemPoint->GetXOut(),
                        gemPoint->GetYOut(),
                        gemPoint->GetZOut()};
  
  Int_t sensorDetId = sensor->GetDetectorId();

  Int_t sensorId = gemPoint->GetSensorId();
  TString nodeName = fDigiPar->GetNodeName(sensorId);
  gGeoManager->cd(nodeName.Data());
  TGeoNode* curNode = gGeoManager->GetCurrentNode();
  
  Double_t locPosIn [4];
  Double_t locPosOut[4];
  
  curNode->MasterToLocal(posIn ,locPosIn );
  curNode->MasterToLocal(posOut,locPosOut);

  Double_t scanOrient = sensor->GetStripOrientation(locPosIn[0],locPosIn[1],side) + TMath::Pi()/2.;

  Double_t feeDist;
  Int_t minChan = sensor->GetChannel2(locPosIn[0],
                                      locPosIn[1],
                                      side,feeDist);
  Int_t maxChan = sensor->GetChannel2(locPosOut[0],
                                      locPosOut[1],
                                      side,feeDist);

  Int_t nofBBB = 1000+10*sensor->GetDistance(side,minChan,maxChan);

  Double_t unitSig = showerStrength/((Double_t)nofBBB);

  Double_t deltaPos[4];
  for ( Int_t ico = 0 ; ico < 4 ; ico++ ) deltaPos[ico] = (locPosOut[ico] - locPosIn[ico])/((Double_t)(nofBBB));

  Int_t chanNr = sensor->GetChannel(locPosIn[0],locPosIn[1],side);

  Double_t time = gemPoint->GetTime();

  if ( scanOrient > TMath::TwoPi() ) 
    scanOrient -= TMath::TwoPi();

  for ( Int_t ibbb = 0 ; ibbb < nofBBB ; ibbb++ ) {
    Double_t dist = fRand->Gaus(0.,showerSigma);

    Double_t xloc = locPosIn[0]+ibbb*deltaPos[0]-dist*TMath::Sin(scanOrient);
    Double_t yloc = locPosIn[1]+ibbb*deltaPos[1]+dist*TMath::Cos(scanOrient);
    chanNr = sensor->GetChannel2(xloc,yloc,side,feeDist);
    time = gemPoint->GetTime()+fRand->Gaus(5,0.5)+fRand->Gaus(feeDist/30.,0.5);
    //time = fTNofEvents;
    if ( chanNr == -1 ) continue;
    ActivateChannel(sensorDetId,side,chanNr,unitSig,time,iPoint);
  }

}

void PndGemDigitize::SimulateRectangularResponse ( Int_t  sensorDetId, Int_t  side, Double_t  channelInd, Double_t  stripWidth, Double_t  showerSigma, Double_t  showerStrength, Int_t  iPoint  )

private

Simulated rectangular response

Definition at line 451 of file PndGemDigitize.cxx.

References ActivateChannel(), and Double_t.

                                                                                                                                                                                     {
  Int_t channelNumber = (Int_t)channelInd;

  Double_t thisASig = 0.;
  
  Double_t leftSide = showerSigma;
  Double_t leftTBin = (channelInd-(Double_t)channelNumber)*stripWidth;
  Int_t    leftChan = channelNumber;
  
  for ( Int_t ibin = 0 ; ibin < 50 ; ibin++ ) {
    Double_t leftASig = leftTBin*showerStrength/showerSigma;
    if ( leftTBin > leftSide ) {
      leftASig = leftSide*showerStrength/showerSigma;
      ibin = 100; // to finish the loop
    }
    
    //cout << "Lwill activate channel " << leftChan << " with " << leftASig << endl;
    if ( ibin == 0 ) 
      thisASig += leftASig;
    else 
      ActivateChannel(sensorDetId,side,leftChan,leftASig,0.,iPoint);
    
    leftSide  = leftSide - leftTBin;
    leftTBin  = stripWidth;
    leftChan -= 1;
  }
  
  Double_t rightSide = showerSigma;
  Double_t rightTBin = (1.+(Double_t)channelNumber-channelInd)*stripWidth;
  Int_t    rightChan = channelNumber;
  
  for ( Int_t ibin = 0 ; ibin < 50 ; ibin++ ) {
    Double_t rightASig = rightTBin*showerStrength/showerSigma;
    if ( rightTBin > rightSide ) {
      rightASig = rightSide*showerStrength/showerSigma;
      ibin = 100; // to finish the loop
    }
    
    //cout << "Rwill activate channel " << rightChan << " with " << rightASig << endl;
    if ( ibin == 0 ) 
      thisASig += rightASig;
    else 
      ActivateChannel(sensorDetId,side,rightChan,rightASig,0.,iPoint);
    
    rightSide  = rightSide - rightTBin;
    rightTBin  = stripWidth;
    rightChan += 1;
  }
  ActivateChannel(sensorDetId,side,channelNumber,thisASig,0.,iPoint);
}

Member Data Documentation

std::map<std::pair<Int_t, Int_t>, Int_t> PndGemDigitize::fChannelMap

private

Map of active channels (pair detectorId, channel number) to index of PndGemDigi

Definition at line 101 of file PndGemDigitize.h.

Referenced by ActivateChannel(), and Reset().

PndGemDigiWriteoutBuffer* PndGemDigitize::fDataBuffer

private

Definition at line 88 of file PndGemDigitize.h.

Referenced by Exec(), and Init().

TClonesArray* PndGemDigitize::fDigiMatches

private

Output array of PndGemDigi

Definition at line 79 of file PndGemDigitize.h.

Referenced by Reset(), and ~PndGemDigitize().

PndGemDigiPar* PndGemDigitize::fDigiPar

private

Definition at line 76 of file PndGemDigitize.h.

Referenced by DigitizeEvent(), DigitizeRealisticEvent(), Finish(), Init(), SetParContainers(), SimulateGaussianResponse(), and ~PndGemDigitize().

TClonesArray* PndGemDigitize::fDigis

private

Input array of PndGemPoint

Definition at line 78 of file PndGemDigitize.h.

Referenced by ActivateChannel(), Exec(), Finish(), Init(), PrintDigis(), Reset(), and ~PndGemDigitize().

TClonesArray* PndGemDigitize::fHitOutsideArray

private

whether to save the hits

Definition at line 82 of file PndGemDigitize.h.

Referenced by Exec(), and Init().

Int_t PndGemDigitize::fNDigis

private

Definition at line 95 of file PndGemDigitize.h.

Referenced by ActivateChannel(), and Reset().

Int_t PndGemDigitize::fNFailed

private

Definition at line 92 of file PndGemDigitize.h.

Referenced by Reset().

Int_t PndGemDigitize::fNMulti

private

Definition at line 94 of file PndGemDigitize.h.

Referenced by Reset().

Int_t PndGemDigitize::fNOutside

private

Definition at line 93 of file PndGemDigitize.h.

Referenced by Reset().

Int_t PndGemDigitize::fNPoints

private

Definition at line 91 of file PndGemDigitize.h.

Referenced by Reset().

TClonesArray* PndGemDigitize::fPoints

private

Definition at line 77 of file PndGemDigitize.h.

Referenced by DigitizeEvent(), DigitizeRealisticEvent(), Exec(), and Init().

TRandom2* PndGemDigitize::fRand

private

Definition at line 87 of file PndGemDigitize.h.

Referenced by SimulateGaussianResponse().

Bool_t PndGemDigitize::fRealisticResponse

private

Definition at line 83 of file PndGemDigitize.h.

Referenced by Exec(), and SetRealisticResponse().

Bool_t PndGemDigitize::fSaveOutsideHits

private

Definition at line 81 of file PndGemDigitize.h.

Referenced by DigitizeEvent(), Exec(), Init(), and SaveOutsideHits().

Double_t PndGemDigitize::fSSigma

private

Output array of PndGemDigiMatches

Definition at line 80 of file PndGemDigitize.h.

Referenced by DigitizeRealisticEvent(), and SetShowerSigma().

Bool_t PndGemDigitize::fTimeOrderedDigi

private

Definition at line 89 of file PndGemDigitize.h.

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

TStopwatch PndGemDigitize::fTimer

private

Definition at line 97 of file PndGemDigitize.h.

Int_t PndGemDigitize::fTNofDigis

private

Definition at line 86 of file PndGemDigitize.h.

Referenced by ActivateChannel(), and Finish().

Int_t PndGemDigitize::fTNofEvents

private

Definition at line 84 of file PndGemDigitize.h.

Referenced by Exec(), and Finish().

Int_t PndGemDigitize::fTNofPoints

private

Definition at line 85 of file PndGemDigitize.h.

Referenced by DigitizeEvent(), DigitizeRealisticEvent(), and Finish().

---

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

• PndGemDigitize.h
• PndGemDigitize.cxx