#include <PndGemEventBuilderOnTracks.h>

Inheritance diagram for PndGemEventBuilderOnTracks:

Public Member Functions
	PndGemEventBuilderOnTracks ()

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

virtual	~PndGemEventBuilderOnTracks ()

virtual void	StoreEventData (FairRecoEventHeader *recoEvent)

std::vector< std::pair< double, FairRecoEventHeader * > >	FindEvents ()

std::vector< std::pair< double, PndTrack * > >	Modify (std::pair< double, PndTrack * > oldData, std::pair< double, PndTrack * > newData)

void	AddNewDataToTClonesArray (FairTimeStamp *data)

double	FindTimeForData (FairTimeStamp *data)

void	FillDataMap (FairTimeStamp *data, double activeTime)

void	EraseDataFromDataMap (FairTimeStamp *data)

Private Member Functions
virtual void	SetParContainers ()

virtual Bool_t	Init ()

virtual Bool_t	ReInit ()

virtual void	Print ()

virtual void	Finish ()

Int_t	CompareTrackToPreviousEvents (Int_t trackId, PndTrack *tempTrack)

	ClassDef (PndGemEventBuilderOnTracks, 1)

Private Attributes
TClonesArray *	fGemTracks

TClonesArray *	fGemOutTracks

Int_t	fTNofEvents

Int_t	fTNofTracks

Int_t	fTNofRecoEvents

Double_t	fGemTrackDelay

TStopwatch	fTimer

Double_t	fExecTime

std::vector< RecoEvent >	fRecoEvents

std::map< FairRecoEventHeader, double >	fEvent_map

std::map< PndTrack, double >	fData_map

Detailed Description

Definition at line 42 of file PndGemEventBuilderOnTracks.h.

Constructor & Destructor Documentation

PndGemEventBuilderOnTracks::PndGemEventBuilderOnTracks ( )

Default constructor

Definition at line 51 of file PndGemEventBuilderOnTracks.cxx.

   : FairEventBuilder("GEMEvTrack","PndTrack","GEM",kTRUE),
     fGemTracks     (NULL),
     fTNofEvents    (0),
     fTNofTracks    (0),
     fTNofRecoEvents(0),
     fExecTime      (0.),
     fGemTrackDelay (0.)
 {
 }

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

Constructor with task name

Definition at line 66 of file PndGemEventBuilderOnTracks.cxx.

References fVerbose, and iVerbose.

   : FairEventBuilder("GEMEvTrack","PndTrack","GEM",kTRUE),
     fGemTracks     (NULL),
     fTNofEvents    (0),
     fTNofTracks    (0),
     fTNofRecoEvents(0),
     fExecTime      (0.),
     fGemTrackDelay (0.)
 {
   SetBuilderName(name);
   fVerbose = iVerbose;
 }

PndGemEventBuilderOnTracks::~PndGemEventBuilderOnTracks ( )

virtual

Destructor

Definition at line 83 of file PndGemEventBuilderOnTracks.cxx.

83 {

84 }

Member Function Documentation

void PndGemEventBuilderOnTracks::AddNewDataToTClonesArray ( FairTimeStamp * data )

Definition at line 212 of file PndGemEventBuilderOnTracks.cxx.

References fVerbose.

 {
   FairRootManager* ioman = FairRootManager::Instance();
   TClonesArray* myArray = ioman->GetTClonesArray(fBranchName);
   if (fVerbose > 1) std::cout << "Data Inserted: "  <<  *(PndTrack*)(data) << std::endl;
   new ((*myArray)[myArray->GetEntries()]) PndTrack(*(PndTrack*)(data));
 }

PndGemEventBuilderOnTracks::ClassDef	(	PndGemEventBuilderOnTracks	,
		1
	)

private

Int_t PndGemEventBuilderOnTracks::CompareTrackToPreviousEvents	(	Int_t	trackId,
		PndTrack *	tempTrack
	)

private

Definition at line 244 of file PndGemEventBuilderOnTracks.cxx.

References CAMath::Abs(), Double_t, fRecoEvents, fVerbose, RecoEvent::meanTime, RecoEvent::nofTracks, and RecoEvent::trackIndex.

Referenced by FindEvents().

                                                                                                  {
   Double_t trackTime = tempTrack->GetTimeStamp();
   Int_t    recoEvent = -1;
   for ( Int_t irev = 0 ; irev < fRecoEvents.size() ; irev++ ) { 
     RecoEvent iterREV = fRecoEvents[irev];
 
     // SHOULD BE SOMETHING MORE SERIOUS THAN "2."
     // 26.08 - just checked the reconstructed track time difference,
     // got a distribution with sigma of 9.49795e-01
     // Setting to 3.0
     if ( TMath::Abs(trackTime-iterREV.meanTime) < 3. ) { 
       if ( fVerbose && recoEvent != -1 ) {
         cout << "track @ " << trackTime << " matches to two reco events: " << irev << " and " << recoEvent << endl;
       }
       else {
         recoEvent = irev;
       }
     }
   }
 
   // matching event found
   if ( recoEvent != -1 ) { 
     fRecoEvents[recoEvent].meanTime = (fRecoEvents[recoEvent].meanTime*fRecoEvents[recoEvent].nofTracks+trackTime)/((Double_t)(fRecoEvents[recoEvent].nofTracks+1));
     fRecoEvents[recoEvent].nofTracks = fRecoEvents[recoEvent].nofTracks+1;
     fRecoEvents[recoEvent].trackIndex.push_back(trackId);
   }
   else {
     RecoEvent tempREV;
     tempREV.meanTime   = trackTime;
     tempREV.nofTracks  = 1;
     tempREV.trackIndex.push_back(trackId);
 
     fRecoEvents.push_back(tempREV);
   }
 
   return recoEvent;
 }

void PndGemEventBuilderOnTracks::EraseDataFromDataMap ( FairTimeStamp * data )

Definition at line 236 of file PndGemEventBuilderOnTracks.cxx.

References fData_map.

 {
   PndTrack myData = *(PndTrack*)data;
   if (fData_map.find(myData) != fData_map.end())
     fData_map.erase(fData_map.find(myData));
 }

void PndGemEventBuilderOnTracks::FillDataMap	(	FairTimeStamp *	data,
		double	activeTime
	)

Definition at line 230 of file PndGemEventBuilderOnTracks.cxx.

References fData_map.

 {
   PndTrack myData = *(PndTrack*)data;
   fData_map[myData] = activeTime;
 }

std::vector< std::pair< double, FairRecoEventHeader * > > PndGemEventBuilderOnTracks::FindEvents ( )

Definition at line 90 of file PndGemEventBuilderOnTracks.cxx.

References CompareRecoEvents(), CompareTrackToPreviousEvents(), fEvent_map, fExecTime, fGemTrackDelay, fGemTracks, fRecoEvents, fTimer, fTNofEvents, fTNofRecoEvents, fTNofTracks, and fVerbose.

                                                                                          {
   if ( fVerbose )
     cout << endl << "======== PndGemEventBuilderOnTracks::FindEvents(Event = " << fTNofEvents << " ) ====================" << endl;
   
   fTimer.Start();
   
   fRecoEvents.clear();
   
   fTNofTracks += fGemTracks->GetEntries();
 
   PndTrack* tempTrack = NULL;
 
   if ( fVerbose )
     cout << "There are " << fGemTracks->GetEntries() << " tracks." << endl;
 
   for ( Int_t itrack = 0 ; itrack < fGemTracks->GetEntries() ; itrack++ ) {
     tempTrack = (PndTrack*)fGemTracks->At(itrack);
     
     if ( fVerbose > 1 ) {
       cout << "    " << fTNofEvents << " . " << itrack << " @ " << tempTrack->GetTimeStamp() << "ns" << endl
            << "         trying to match to one of " << fRecoEvents.size() << " reco events" << endl;
     }
 
     Int_t revNo = CompareTrackToPreviousEvents(itrack,tempTrack);
 
     FillNewData(tempTrack,
                 tempTrack->GetTimeStamp(),
                 tempTrack->GetTimeStamp());
 
   }
 
   fTNofEvents += 1;
 
   fTNofRecoEvents += fRecoEvents.size();
 
   std::sort(fRecoEvents.begin(),fRecoEvents.end(),CompareRecoEvents);
 
   std::vector<std::pair<double, FairRecoEventHeader*> > result;
   std::pair<double,FairRecoEventHeader*> singleResult;
 
   if ( fVerbose )
     cout << "There are " << fRecoEvents.size() << " events." << endl;
 
   for ( Int_t irev = 0 ; irev < fRecoEvents.size() ; irev++ ) {
     FairRecoEventHeader recoEvent;
     recoEvent.SetEventTime(fRecoEvents[irev].meanTime-fGemTrackDelay,5.);
     recoEvent.SetIdentifier(GetIdentifier());
     SetMaxAllowedTime(recoEvent.GetEventTime());
 
     fEvent_map[recoEvent] = recoEvent.GetEventTime();
 
     singleResult.first  = recoEvent.GetEventTime();
     singleResult.second = (FairRecoEventHeader*)&fEvent_map.rbegin()->first;
 
     if ( fVerbose )
       cout << "FOUND EVENT @ " << singleResult.first << endl;
     
     result.push_back(singleResult);
   }
 
   fExecTime += fTimer.RealTime();
   fTimer.Stop();  
 
   return result;
 }

double PndGemEventBuilderOnTracks::FindTimeForData ( FairTimeStamp * data )

Definition at line 220 of file PndGemEventBuilderOnTracks.cxx.

References fData_map.

 {
   std::map<PndTrack, double>::iterator it;
   PndTrack myData = *(PndTrack*)data;
   it = fData_map.find(myData);
   if (it == fData_map.end())
     return -1;
   else
     return it->second;
 }

void PndGemEventBuilderOnTracks::Finish ( )

privatevirtual

Finish at the end of each event

Definition at line 334 of file PndGemEventBuilderOnTracks.cxx.

References Double_t, fExecTime, fTNofEvents, fTNofRecoEvents, and fTNofTracks.

                                         {
 
   cout << "-------------------- " << GetBuilderName() << " : Summary -----------------------" << endl;
   cout << " Events:        " << setw(10) << fTNofEvents     << endl;
   cout << " Tracks:        " << setw(10) << fTNofTracks     << "    ( " << (Double_t)fTNofTracks   /((Double_t)fTNofEvents) << " per event )" << endl;
   cout << " Reco Events:   " << setw(10) << fTNofRecoEvents << endl;
   cout << " . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . " << endl;
   cout << " >>> EB >>> all   time  = " << fExecTime << "s" << endl;
   cout << "---------------------------------------------------------------------" << endl; 
 }

Bool_t PndGemEventBuilderOnTracks::Init ( )

privatevirtual

Intialisation

Definition at line 299 of file PndGemEventBuilderOnTracks.cxx.

References fGemTrackDelay, and fGemTracks.

                                         {
 
   // Get input array
   FairRootManager* ioman = FairRootManager::Instance();
   if ( ! ioman ) Fatal("Init", "No FairRootManager");
 
   fGemTracks = (TClonesArray*) ioman->GetObject("GEMTrack");
 
   // fGemOutTracks = new TClonesArray("PndTrack",100);
   // ioman->Register("GEMTrackQQQ", "Gem Tracks", fGemOutTracks, kTRUE);
   // ioman->RegisterWriteoutBuffer("GEMDigi", fDataBuffer);
   // fDataBuffer->ActivateBuffering(fTimeOrderedDigi);
 
   fGemTrackDelay = 4.96229;
 
   return kTRUE;
 }

std::vector< std::pair< double, PndTrack * > > PndGemEventBuilderOnTracks::Modify	(	std::pair< double, PndTrack * >	oldData,
		std::pair< double, PndTrack * >	newData
	)

Definition at line 205 of file PndGemEventBuilderOnTracks.cxx.

 {
   std::vector<std::pair<double, PndTrack*> > result;
   result.push_back(newData);
   return result;
 }

void PndGemEventBuilderOnTracks::Print ( )

privatevirtual

Definition at line 158 of file PndGemEventBuilderOnTracks.cxx.

References fEvent_map.

                                        {
   cout << " print " << endl;
   cout << " fEventMap has " << fEvent_map.size() << " entries" << endl;
   std::map<FairRecoEventHeader,double>::iterator iter;
   for ( iter = fEvent_map.begin() ; iter != fEvent_map.end() ; iter++ ) {
     cout << " // for time " << iter->second << " got " << iter->first.GetEventTime() << " at an address " << &iter->first << endl;
   }
   cout << "-------" << endl;
 }

Bool_t PndGemEventBuilderOnTracks::ReInit ( )

privatevirtual

Reinitialisation

Definition at line 322 of file PndGemEventBuilderOnTracks.cxx.

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

void PndGemEventBuilderOnTracks::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 284 of file PndGemEventBuilderOnTracks.cxx.

References 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");
 }

void PndGemEventBuilderOnTracks::StoreEventData ( FairRecoEventHeader * recoEvent )

virtual

Definition at line 170 of file PndGemEventBuilderOnTracks.cxx.

References CAMath::Abs(), fData_map, fGemTrackDelay, and fVerbose.

                                                                               {
   typedef std::multimap<double, std::pair<double, FairTimeStamp*> >::iterator startTimeMapIter;
 
   if ( fVerbose ) {
     cout << "* " << fData_map.size() << " in fData_map " 
          << "* " << fStartTime_map.size() << " in fStartTime_map "
          << " TO PUT IN EVENT AT " << recoEvent->GetEventTime() << endl;
     cout << "CHECK TRACKS WITH TimeStamps:" << endl;
     for (startTimeMapIter jter = fStartTime_map.begin() ; jter != fStartTime_map.end(); ++jter) {
       std::pair<double, FairTimeStamp*> datb = jter->second;
       cout << "@@@ " << datb.second->GetTimeStamp() << endl;
     }
   }
   
   startTimeMapIter iter = fStartTime_map.begin();
   while ( iter != fStartTime_map.end() ) {
     std::pair<double, FairTimeStamp*> data = iter->second;
     if ( fVerbose ) {
       cout << " ---> check data at " << data.second->GetTimeStamp() 
            << " /// " << data.second->GetTimeStamp()-recoEvent->GetEventTime()-fGemTrackDelay << endl;
     }
     if ( TMath::Abs(data.second->GetTimeStamp()-recoEvent->GetEventTime()-fGemTrackDelay) < 3. ) {
       if ( fVerbose ) cout << "WILL ADD THIS ONE" << endl;
       FillDataToDeadTimeMap(data.second, data.first);
 
       startTimeMapIter save = iter;
       save++;
       fStartTime_map.erase(iter);
       iter = save;
     }
     else 
       iter++;
   }
 }