#include <PndGemFindClusters.h>

Inheritance diagram for PndGemFindClusters:

Public Member Functions
	PndGemFindClusters ()

	PndGemFindClusters (Int_t iVerbose)

	PndGemFindClusters (const char *name, Int_t iVerbose=0)

virtual	~PndGemFindClusters ()

virtual void	Exec (Option_t *opt)

Private Member Functions
virtual void	SetParContainers ()

virtual InitStatus	Init ()

virtual InitStatus	ReInit ()

Int_t	SortDigis ()

Int_t	CreateClusters ()

Int_t	WriteClusters ()

Bool_t	CompareDigiToClustersDigis (Int_t digiNumber)

void	SortClusters ()

void	SortCluster (Int_t clus)

void	PrintClusters ()

void	PrintCluster (Int_t clus)

void	JoinTwoClusters (Int_t clus1, Int_t clus2)

void	AddDigiToCluster (Int_t digiNr, Int_t clusNr)

virtual void	Finish ()

	ClassDef (PndGemFindClusters, 1)

Private Attributes
PndGemDigiPar *	fDigiPar

TClonesArray *	fDigis

TClonesArray *	fClusters

std::map< Double_t, Int_t >	fTimeOrderedDigis

std::vector< DigiCluster >	fDigiClusters

Int_t	fTNofEvents

Int_t	fTNofDigis

Int_t	fTNofClusters

TStopwatch	fTimer

Detailed Description

Definition at line 51 of file PndGemFindClusters.h.

Constructor & Destructor Documentation

PndGemFindClusters::PndGemFindClusters ( )

Default constructor

Definition at line 50 of file PndGemFindClusters.cxx.

References fClusters, fDigiPar, fDigis, fTNofClusters, fTNofDigis, and fTNofEvents.

   : FairTask("GEMFindClusters", 1) {
   fDigiPar  = NULL;
   fDigis    = NULL;
   fClusters = NULL;
 
   fTNofEvents   = 0;
   fTNofDigis    = 0;
   fTNofClusters = 0;
 }

PndGemFindClusters::PndGemFindClusters ( Int_t iVerbose )

Standard constructor

Definition at line 65 of file PndGemFindClusters.cxx.

References fClusters, fDigiPar, fDigis, fTNofClusters, fTNofDigis, and fTNofEvents.

   : FairTask("GEMFindClusters", iVerbose) {
   fDigiPar  = NULL;
   fDigis    = NULL;
   fClusters = NULL;
 
   fTNofEvents   = 0;
   fTNofDigis    = 0;
   fTNofClusters = 0;
 }

PndGemFindClusters::PndGemFindClusters	(	const char *	name,
		Int_t	iVerbose = `0`
	)

Constructor with task name

Definition at line 80 of file PndGemFindClusters.cxx.

References fClusters, fDigiPar, fDigis, fTNofClusters, fTNofDigis, and fTNofEvents.

   : FairTask(name, iVerbose) {
   fDigiPar  = NULL;
   fDigis    = NULL;
   fClusters = NULL;
 
   fTNofEvents   = 0;
   fTNofDigis    = 0;
   fTNofClusters = 0;
 }

PndGemFindClusters::~PndGemFindClusters ( )

virtual

Destructor

Definition at line 95 of file PndGemFindClusters.cxx.

References fClusters.

                                         {
   if ( fClusters ) {
     fClusters->Delete();
     delete fClusters;
   }
 }

Member Function Documentation

void PndGemFindClusters::AddDigiToCluster	(	Int_t	digiNr,
		Int_t	clusNr
	)

private

Definition at line 948 of file PndGemFindClusters.cxx.

References digi, Double_t, fDigiClusters, fDigiPar, fDigis, fVerbose, PndGemDigi::GetChannelNr(), PndGemDigi::GetCharge(), PndGemDigi::GetDetectorId(), PndGemSensor::GetMeanChannel(), PndGemDigiPar::GetSensor(), PndGemDigi::GetSensorNr(), PndGemDigi::GetSide(), PndGemDigi::GetStationNr(), and sensor.

Referenced by CompareDigiToClustersDigis(), and JoinTwoClusters().

                                                                     {
   if ( fVerbose > 0 ) cout << "PndGemFindClusters::AddDigiToCluster()" << endl;
   
   PndGemDigi* digi = (PndGemDigi*) fDigis->At(digiNr);
   Int_t stationNr = digi->GetStationNr();
   Int_t sensorNr  = digi->GetSensorNr();
   Int_t iSide     = digi->GetSide();
   
   PndGemSensor* sensor = fDigiPar->GetSensor(stationNr, sensorNr);
 
   if ( fDigiClusters[clusNr].digiNr.size() == 0 ) {
     //    cout << "NEW CLUSTER " << clusNr << " CREATED WITH AddDigiToCluster(" << digiNr << "," << clusNr << ") IN DETECTOR " << digi->GetDetectorId() << endl;
     fDigiClusters[clusNr].detId  = digi->GetDetectorId();
     fDigiClusters[clusNr].cluPos = digi->GetChannelNr();
     fDigiClusters[clusNr].cluADC = digi->GetCharge();
     fDigiClusters[clusNr].cluTDC = digi->GetTimeStamp();
     fDigiClusters[clusNr].cluPMn = digi->GetChannelNr();
     fDigiClusters[clusNr].cluPMx = digi->GetChannelNr();
     fDigiClusters[clusNr].cluMVl = digi->GetCharge();
     fDigiClusters[clusNr].cluMPs = digi->GetChannelNr();
   }
   else {
     fDigiClusters[clusNr].cluTDC = (fDigiClusters[clusNr].cluTDC*((Double_t)(fDigiClusters[clusNr].sigADC.size()-1.))+digi->GetTimeStamp())/((Double_t)(fDigiClusters[clusNr].sigADC.size()));
     fDigiClusters[clusNr].cluPos = sensor->GetMeanChannel(iSide,fDigiClusters[clusNr].cluPos,fDigiClusters[clusNr].cluADC,digi->GetChannelNr(),digi->GetCharge());
     if ( fDigiClusters[clusNr].cluPMn > digi->GetChannelNr() ) fDigiClusters[clusNr].cluPMn = digi->GetChannelNr();
     if ( fDigiClusters[clusNr].cluPMx < digi->GetChannelNr() ) fDigiClusters[clusNr].cluPMx = digi->GetChannelNr();
     fDigiClusters[clusNr].cluADC =  fDigiClusters[clusNr].cluADC+digi->GetCharge();
     if ( digi->GetCharge() > fDigiClusters[clusNr].cluMVl ) {
       fDigiClusters[clusNr].cluMVl = digi->GetCharge();
       fDigiClusters[clusNr].cluMPs = digi->GetChannelNr();
     }
   }
 
   fDigiClusters[clusNr].digiNr.push_back(digiNr);
   fDigiClusters[clusNr].chanNr.push_back((Int_t)digi->GetChannelNr());
   fDigiClusters[clusNr].sigADC.push_back(digi->GetCharge());
 }

PndGemFindClusters::ClassDef	(	PndGemFindClusters	,
		1
	)

private

Bool_t PndGemFindClusters::CompareDigiToClustersDigis ( Int_t digiNumber )

private

Compare digi to clusters Compare digi to digis

Definition at line 877 of file PndGemFindClusters.cxx.

References AddDigiToCluster(), digi, Double_t, fDigiClusters, fDigiPar, fDigis, fVerbose, PndGemDigi::GetChannelNr(), PndGemDigi::GetDetectorId(), PndGemSensor::GetDistance2(), PndGemDigiPar::GetSensor(), PndGemDigi::GetSensorNr(), PndGemDigi::GetSide(), PndGemDigi::GetStationNr(), JoinTwoClusters(), and sensor.

Referenced by CreateClusters().

                                                                       {
   if ( fVerbose > 0 ) cout << "PndGemFindClusters::CompareDigiToClustersDigis" << endl;
   
   PndGemSensor* sensor;
   PndGemDigi* digi;
   digi = (PndGemDigi*) fDigis->At(digiNumber);
 
   Int_t iDigiUsed = -1;
   for ( size_t idc = 0 ; idc < fDigiClusters.size() ; idc++ ) {
     if ( fDigiClusters[idc].cluPos < -0.5 ) continue;
     if ( digi->GetDetectorId() != fDigiClusters[idc].detId ) continue;
 
     Int_t stationNr = digi->GetStationNr();
     Int_t sensorNr  = digi->GetSensorNr();
     Int_t iSide     = digi->GetSide();
 
     sensor = fDigiPar->GetSensor(stationNr, sensorNr);
 
     Double_t digiDigiDist=999999999, testDist;
     for ( size_t id = 0 ; id < fDigiClusters[idc].digiNr.size() ; id++ ) {
       testDist = sensor->GetDistance2(iSide,digi->GetChannelNr(),fDigiClusters[idc].chanNr[id]);
       if ( testDist < digiDigiDist  ) digiDigiDist=testDist;
     }
     
     if ( digiDigiDist < 0. ) continue; 
     if ( digiDigiDist > 1 ) continue;//1 for next strip to accept
 
     if ( fVerbose > 1 ) 
       cout << "digi " << digiNumber << " ( " << digi->GetDetectorId() << " / " << digi->GetChannelNr() << " ) could attach to cluster " 
            << idc << " ( " << fDigiClusters[idc].cluPos << " ) with distance " << digiDigiDist << " . " << fDigiClusters[idc].digiNr.size() << " digis in cluster" << flush;
     
     AddDigiToCluster(digiNumber,idc);
 
     if ( fVerbose > 1 ) 
       cout << " >>> mean channel = " << fDigiClusters[idc].cluPos << endl;
 
     if ( iDigiUsed >= 0 ) {
       //cout << "DIGI " << digiNumber << " MATCHES TO CLUSTER " << iDigiUsed << " and " << idc << endl;
       JoinTwoClusters(idc,iDigiUsed);
     }
     
     iDigiUsed = idc;
   }
   return (iDigiUsed==-1?kFALSE:kTRUE);
 
 }

Int_t PndGemFindClusters::CreateClusters ( )

private

Find clusters Look for clusters

Definition at line 322 of file PndGemFindClusters.cxx.

References DigiCluster::chanNr, DigiCluster::cluADC, DigiCluster::cluMPs, DigiCluster::cluMVl, DigiCluster::cluPMn, DigiCluster::cluPMx, DigiCluster::cluPos, DigiCluster::cluTDC, CompareDigiToClustersDigis(), DigiCluster::detId, digi, DigiCluster::digiNr, fDigiClusters, fDigis, fVerbose, PndGemDigi::GetChannelNr(), PndGemDigi::GetCharge(), PndGemDigi::GetDetectorId(), and DigiCluster::sigADC.

Referenced by Exec().

                                          {
   if ( fVerbose > 0 ) {
     cout << "-I-  PndGemFindClusters::CreateClusters()" <<endl; 
   }
   
   fDigiClusters.clear();
 
   PndGemDigi* digi = NULL;
   for ( Int_t idigi = 0 ; idigi < fDigis->GetEntries() ; idigi++ ) {
     digi = (PndGemDigi*) fDigis->At(idigi);
     
     if ( fVerbose ) {
       //    cout << "digi " << idigi << " in " << digi->GetDetectorId() << ", there are already " << fDigiClusters.size() << " clusters" << endl;
       //    cout << idigi << " ---> " << digi->GetDetectorId() << " @ " << digi->GetChannelNr() << " @ " << digi->GetTimeStamp() << endl;
     }
     if ( CompareDigiToClustersDigis(idigi) == kTRUE ) continue;
 
     DigiCluster tempDC;
     tempDC.detId  = digi->GetDetectorId();
     tempDC.digiNr.push_back(idigi);
     tempDC.chanNr.push_back((Int_t)digi->GetChannelNr());
     tempDC.sigADC.push_back(digi->GetCharge());
     tempDC.cluTDC = digi->GetTimeStamp();
     tempDC.cluADC = digi->GetCharge();
     tempDC.cluPos = digi->GetChannelNr();
     tempDC.cluPMn = digi->GetChannelNr();
     tempDC.cluPMx = digi->GetChannelNr();
     tempDC.cluMPs = digi->GetChannelNr();
     tempDC.cluMVl = digi->GetCharge();
     fDigiClusters.push_back(tempDC);
 
   }
   return fDigiClusters.size();
 }

void PndGemFindClusters::Exec ( Option_t * opt )

virtual

Execution

Definition at line 106 of file PndGemFindClusters.cxx.

References CreateClusters(), fClusters, fDigiClusters, fTimer, fTNofClusters, fTNofEvents, fVerbose, PrintClusters(), SortClusters(), SortDigis(), and WriteClusters().

                                        {
   if ( fVerbose > 0 ) cout <<"-I- PndGemFindClusters::Exec()"<<endl;
   
   fTNofEvents++;
 
   fTimer.Start();
   //Bool_t warn = kFALSE; //[R.K. 01/2017] unused variable?
 
   // Clear output array
   fClusters->Clear();
 
   //Int_t nDigis    = 0; //[R.K.03/2017] unused variable
   Int_t nClusters = 0;
 
   // Sort digis with respect to time
   SortDigis();//nDigis =  //[R.K.03/2017] unused variable
 
   //FindClusters();
   
   CreateClusters();
 
   if ( fVerbose > 1 ) {
     Int_t nofC2W = 0;
     for ( size_t idc = 0 ; idc < fDigiClusters.size() ; idc++ )
       if ( fDigiClusters[idc].cluPos > -0.5 )
         nofC2W++;
     cout << ">>>> created  " << fDigiClusters.size() << " clusters, " << nofC2W << " to write" << endl;
   }
   
   //Skip Analyze... the join is already done in the FindClusters()..
   //AnalyzeClusters();
   
 //   if ( fVerbose ) {
 //     Int_t nofC2W = 0;
 //     for ( Int_t idc = 0 ; idc < fDigiClusters.size() ; idc++ )
 //       if ( fDigiClusters[idc].cluPos > -0.5 )
 //      nofC2W++;
 //     cout << ">>>> analyzed " << fDigiClusters.size() << " clusters, " << nofC2W << " to write" << endl;
 //   }
 
   //ClearClusters();
   //ClearClusters2();//sort and print only.. 
   SortClusters(); 
 
   if ( fVerbose > 1 ) {
     PrintClusters();
   }
   
 //   if ( fVerbose ) {
 //     Int_t nofC2W = 0;
 //     for ( Int_t idc = 0 ; idc < fDigiClusters.size() ; idc++ )
 //       if ( fDigiClusters[idc].cluPos > -0.5 )
 //      nofC2W++;
 //     cout << ">>>> cleared  " << fDigiClusters.size() << " clusters, " << nofC2W << " to write" << endl;
 //   }
 
   nClusters = WriteClusters();
 
   fTimer.Stop();  
 
   //  cout << "CREATED " << nClusters << " CLUSTERS OUT OF " << nDigis << " DIGIS IN " << fTimer.RealTime() << "s" << endl;
 
   fTNofClusters += nClusters;
 }

void PndGemFindClusters::Finish ( )

privatevirtual

Finish at the end of each event

Definition at line 1063 of file PndGemFindClusters.cxx.

References Double_t, fClusters, fTNofClusters, fTNofDigis, and fTNofEvents.

                                 {
   if ( fClusters ) fClusters->Clear();
 
   cout << "-------------------- " << fName.Data() << " : Summary -----------------------" << endl;
   cout << " Events:        " << setw(10) << fTNofEvents << endl;
   cout << " Digis:         " << setw(10) << fTNofDigis      << "    ( " << (Double_t)fTNofDigis   /((Double_t)fTNofEvents) << " per event )" << endl;
   cout << " Clusters:      " << setw(10) << fTNofClusters   << "    ( " << (Double_t)fTNofClusters/((Double_t)fTNofEvents) << " per event )" << endl;
   cout << "                       -->    ( " << (Double_t)fTNofDigis  /((Double_t)fTNofClusters) << " digis per cluster )" << endl;
   cout << "---------------------------------------------------------------------" << endl; 
 
 }

InitStatus PndGemFindClusters::Init ( )

privatevirtual

Intialisation

Definition at line 193 of file PndGemFindClusters.cxx.

References fClusters, fDigiPar, fDigis, PndGemDigiPar::GetNStations(), PndGemSensor::GetOuterRadius(), PndGemStation::GetSensor(), PndGemDigiPar::GetStation(), and sensor.

                                     {
 
   // Get input array
   FairRootManager* ioman = FairRootManager::Instance();
   if ( ! ioman ) Fatal("Init", "No FairRootManager");
   fDigis = (TClonesArray*) ioman->GetObject("GEMDigi");
 
   // Register output array
   //fClusters = new TClonesArray("PndGemDigi", 1000);
   fClusters = new TClonesArray("PndGemCluster", 1000);
   ioman->Register("GEMCluster", "Cluster in GEM", fClusters, kTRUE);
 
   cout << "-I- " << fName.Data() << "::Init(). There are " << fDigiPar->GetNStations() << " GEM stations." << endl;
   cout << "-I- " << fName.Data() << "::Init(). Initialization succesfull." << endl;  
   
   PndGemStation* station = (PndGemStation*)fDigiPar->GetStation(0);
   PndGemSensor* sensor = (PndGemSensor*)station->GetSensor(0);
   cout << "sensor out rad is " << sensor->GetOuterRadius() << endl;
 
   return kSUCCESS;
 }

void PndGemFindClusters::JoinTwoClusters	(	Int_t	clus1,
		Int_t	clus2
	)

private

Definition at line 926 of file PndGemFindClusters.cxx.

References AddDigiToCluster(), fDigiClusters, and fVerbose.

Referenced by CompareDigiToClustersDigis().

                                                                  {
   if ( fVerbose > 0 ) cout << "PndGemFindClusters::JointTwoClusters()" << endl;
   
   //PndGemDigi* digi; //[R.K. 01/2017] unused variable?
 
   for ( size_t id2 = 0 ; id2 < fDigiClusters[clus2].digiNr.size() ; id2++ ) { 
     Int_t digi2 = fDigiClusters[clus2].digiNr[id2];
     Int_t chan2 = fDigiClusters[clus2].chanNr[id2];
 
     Int_t addDigi = kTRUE;
     for ( size_t id1 = 0 ; id1 < fDigiClusters[clus1].digiNr.size() ; id1++ ) {
       if ( chan2 == fDigiClusters[clus1].chanNr[id1] ) addDigi = kFALSE;
     }
     
     if ( addDigi ) AddDigiToCluster(digi2,clus1);
   }
   fDigiClusters[clus2].cluPos = -1.;
 }

void PndGemFindClusters::PrintCluster ( Int_t clus )

private

Definition at line 1033 of file PndGemFindClusters.cxx.

References Double_t, fDigiClusters, and fVerbose.

Referenced by PrintClusters().

                                                 {
   if ( fVerbose > 0 ) cout << "PndGemFindClusters::PrintCluster()" << endl;
   
   cout << "CLUSTER " << clus << " IN " << fDigiClusters[clus].detId << " AT " << fDigiClusters[clus].cluPos 
        << " CREATED FROM " << fDigiClusters[clus].digiNr.size() << " DIGIS, FROM " 
        << fDigiClusters[clus].cluPMn << " TO " << fDigiClusters[clus].cluPMx << " WITH MAXIMUM OF "
        << fDigiClusters[clus].cluMVl << " AT " << fDigiClusters[clus].cluMPs << endl;
   
   for ( Int_t icol = 0 ; icol < fDigiClusters[clus].cluPMx-fDigiClusters[clus].cluPMn+3 ; icol++ ) cout << "-" << flush; cout<<endl;
   for ( Double_t fh = 0.99 ; fh >= 0. ; fh-=0.247475 ) {
     cout << "|" << flush;
     Int_t    thisDigi = 0;
     Double_t thisADC = 0.;
     for ( Int_t ichan = fDigiClusters[clus].cluPMn ; ichan <= fDigiClusters[clus].cluPMx ; ichan++ ) {
       thisADC = 0.;
       if ( fDigiClusters[clus].chanNr[thisDigi] == ichan ) {
         thisADC = fDigiClusters[clus].sigADC[thisDigi];
         thisDigi++;
       }
       if ( thisADC > fh*fDigiClusters[clus].cluMVl ) cout << "*" << flush;
       else                                          cout << " " << flush;
     }
     cout << "|" << endl;
   }
   for ( Int_t icol = 0 ; icol < fDigiClusters[clus].cluPMx-fDigiClusters[clus].cluPMn+3 ; icol++ ) cout << "-" << flush; cout<<endl;
 }

void PndGemFindClusters::PrintClusters ( )

private

Definition at line 1021 of file PndGemFindClusters.cxx.

References fDigiClusters, fVerbose, and PrintCluster().

Referenced by Exec().

                                        {
   if ( fVerbose > 0 ) cout << "PndGemFindClusters::PrintClusters()" << endl;
 
   for ( size_t idc = 0 ; idc < fDigiClusters.size() ; idc++ ) {
     //    if ( fDigiClusters[idc].detId != 568329024 ) continue;
     if ( fDigiClusters[idc].cluPos < -0.5 ) continue;
     PrintCluster(idc);
   }
 }

InitStatus PndGemFindClusters::ReInit ( )

privatevirtual

Reinitialisation

Definition at line 220 of file PndGemFindClusters.cxx.

                                       {
 
   // Create sectorwise digi sets
   //  MakeSets();
 
   return kSUCCESS;
 }

void PndGemFindClusters::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 175 of file PndGemFindClusters.cxx.

References fDigiPar, and run.

                                           {
 
   // Get run and runtime database
   FairRunAna* run = FairRunAna::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 PndGemFindClusters::SortCluster ( Int_t clus )

private

Definition at line 997 of file PndGemFindClusters.cxx.

References Double_t, fDigiClusters, and fVerbose.

Referenced by SortClusters().

                                                {
   if ( fVerbose > 0 ) cout << "PndGemFindClusters::SortCluster()" << endl;
   Int_t channelPosition = 0;
 
   for ( Int_t ich = fDigiClusters[clus].cluPMn ; ich <= fDigiClusters[clus].cluPMx ; ich++ ) {
     for ( size_t idigi = channelPosition ; idigi < fDigiClusters[clus].digiNr.size() ; idigi++ ) {
       if ( fDigiClusters[clus].chanNr[idigi] == ich ) {
         Int_t    dN = fDigiClusters[clus].digiNr[idigi];
         Int_t    cN = fDigiClusters[clus].chanNr[idigi];
         Double_t sA = fDigiClusters[clus].sigADC[idigi];
         fDigiClusters[clus].digiNr[idigi] = fDigiClusters[clus].digiNr[channelPosition];
         fDigiClusters[clus].chanNr[idigi] = fDigiClusters[clus].chanNr[channelPosition];
         fDigiClusters[clus].sigADC[idigi] = fDigiClusters[clus].sigADC[channelPosition];
         fDigiClusters[clus].digiNr[channelPosition] = dN;
         fDigiClusters[clus].chanNr[channelPosition] = cN;
         fDigiClusters[clus].sigADC[channelPosition] = sA;
         channelPosition++;
       }
     }
   }
 }

void PndGemFindClusters::SortClusters ( )

private

Analyze clusters Clear clusters

Definition at line 988 of file PndGemFindClusters.cxx.

References fDigiClusters, fVerbose, and SortCluster().

Referenced by Exec().

                                       {
   if ( fVerbose > 0 ) cout << "PndGemFindClusters::SortClusters()" << endl;
   for ( size_t idc = 0 ; idc < fDigiClusters.size() ; idc++ ) {
     SortCluster(idc);
   }
 }

Int_t PndGemFindClusters::SortDigis ( )

private

Sort digis timewise

Definition at line 231 of file PndGemFindClusters.cxx.

References Double_t, fDigis, fTimeOrderedDigis, fTNofDigis, and fVerbose.

Referenced by Exec().

                                     {
    if ( fVerbose ) cout << "-I-  PndGemFindClusters::SortDigis()" <<endl;
 
   // Check input array
   if ( ! fDigis ) {
     cout << "-E- " << fName << "::SortDigis: No input array!" << endl;
     return -1;
   }
   TRandom2* randNumber = new TRandom2();
 
   //PndGemDigi* digi; //[R.K. 01/2017] unused variable?
   Int_t nDigis = fDigis->GetEntriesFast();
 
   fTimeOrderedDigis.clear();
 
   fTNofDigis += nDigis;
 
   for (Int_t iDigi=0; iDigi<nDigis; iDigi++) {
     //digi = (PndGemDigi*) fDigis->At(iDigi); //[R.K. 01/2017] unused variable?
     //    Double_t time = digi->GetTimeStamp();
     Double_t time = randNumber->Gaus(0.,5.);
 
     fTimeOrderedDigis[time] = iDigi;
   }
 
   //  cout << "have " << nDigis << " digis ( " << fTimeOrderedDigis.size() << " ) " << endl;
 
   map<Double_t, Int_t>::iterator mapIt;
   Int_t nd = 0;
   for (mapIt=fTimeOrderedDigis.begin(); mapIt!=fTimeOrderedDigis.end(); mapIt++) {
     //    cout << mapIt->first << " --> " << mapIt->second << endl;
     nd++;
   }
   // cout << " ---> " << nd << endl;
 
   return nDigis;
 }

Int_t PndGemFindClusters::WriteClusters ( )

private

Write clusters

Definition at line 790 of file PndGemFindClusters.cxx.

References fDigiClusters, and fVerbose.

Referenced by Exec().

                                         {
   if ( fVerbose > 0 ) cout << "-I- PndGemFindClusters::WriteClusters()" <<endl;
   
   Int_t nClusters = 0;
   //PndGemDigi* digi; //[R.K. 01/2017] unused variable?
 
   for ( size_t idc = 0 ; idc < fDigiClusters.size() ; idc++ ) {
     
     //    if ( fDigiClusters[idc].cluADC < 1. ) continue;
     if ( fDigiClusters[idc].cluPos < -0.5 ) continue;
     
     if ( fVerbose > 1 ) {      
       cout << "cluster at " << fDigiClusters[idc].cluPos << ", h = " << fDigiClusters[idc].cluADC << ", for detId = " << fDigiClusters[idc].detId << " from " << fDigiClusters[idc].digiNr.size() << " digis (from " << fDigiClusters[idc].cluPMn << " to " << fDigiClusters[idc].cluPMx << "): " << endl;
       
       for ( size_t idigi = 0 ; idigi < fDigiClusters[idc].digiNr.size() ; idigi++ ) {
         cout <<" digiNr="  << fDigiClusters[idc].digiNr[idigi] 
              <<" chanNr=" << fDigiClusters[idc].chanNr[idigi] 
              <<" sigADC=" << fDigiClusters[idc].sigADC[idigi] << endl;
       }
     }
     
     new ((*fClusters)[nClusters]) PndGemCluster(fDigiClusters[idc].detId, 
                                                           fDigiClusters[idc].cluPos,
                                                           fDigiClusters[idc].cluPMn,
                                                           fDigiClusters[idc].cluPMx,
                                                           fDigiClusters[idc].cluADC,
                                                           fDigiClusters[idc].cluTDC,
                                                           fDigiClusters[idc].digiNr);
     
     if ( fVerbose > 1 ) 
       cout << "creating cluster at " << fDigiClusters[idc].cluPos << " with height of " << fDigiClusters[idc].cluADC << " /// compare maximum at " << fDigiClusters[idc].cluMPs << endl;
     nClusters++;
   }
   return nClusters;
 }