PndGemDigiAna Class Reference

analyze GEM digis More...

#include <PndGemDigiAna.h>

Inheritance diagram for PndGemDigiAna:

Public Member Functions
	PndGemDigiAna ()
	PndGemDigiAna (Int_t iVerbose)
	PndGemDigiAna (TString taskName, Int_t iVerbose)
virtual	~PndGemDigiAna ()
virtual void	Exec (Option_t *opt)
void	SetVerbose (const Int_t &verbose)

Private Member Functions
void	CreateHistos ()
void	AnaHistos ()
virtual void	SetParContainers ()
virtual void	Finish ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
	ClassDef (PndGemDigiAna, 1)

Private Attributes
PndGemDigiPar *	fDigiPar
TClonesArray *	fGemDigiArray
Int_t	fNofEvents
	event counter More...
TList *	fHistoList
TClonesArray *	fHistoArray

Detailed Description

analyze GEM digis

Author

R. Karabowicz r.kar.nosp@m.abow.nosp@m.icz@g.nosp@m.si.d.nosp@m.e

Date

19.03.2009 Compare the found digis. First version compares only the time between digis. For each strip the time of last digi on that strip is stored: DLTS By default the DLTS is set to -1. If strip is fired and the DLTS is different from -1, the histogram with time between digis is filled. Moreover, at any given time the DLTS will show for how long was the strip inactive.

Definition at line 33 of file PndGemDigiAna.h.

Constructor & Destructor Documentation

PndGemDigiAna::PndGemDigiAna

(

)

Default constructor

Definition at line 33 of file PndGemDigiAna.cxx.

  : FairTask("GEM Digis Ana", 1),
    fDigiPar(NULL),
    fGemDigiArray(NULL),
    fNofEvents(0),
    fHistoList(NULL),
    fHistoArray(NULL)
{
}

PndGemDigiAna::PndGemDigiAna

(

Int_t

iVerbose

)

Definition at line 45 of file PndGemDigiAna.cxx.

  : FairTask("GEM Digis Ana", iVerbose) ,
    fDigiPar(NULL),
    fGemDigiArray(NULL),
    fNofEvents(0),
    fHistoList(NULL),
    fHistoArray(NULL)
{
}

PndGemDigiAna::PndGemDigiAna	(	TString	taskName,
		Int_t	iVerbose
	)

Definition at line 57 of file PndGemDigiAna.cxx.

  : FairTask(taskName, iVerbose) ,
    fDigiPar(NULL),
    fGemDigiArray(NULL),
    fNofEvents(0),
    fHistoList(NULL),
    fHistoArray(NULL)
{
}

PndGemDigiAna::~PndGemDigiAna ( )

virtual

Destructor

Definition at line 69 of file PndGemDigiAna.cxx.

{ }

Member Function Documentation

void PndGemDigiAna::AnaHistos ( )

private

Definition at line 245 of file PndGemDigiAna.cxx.

Referenced by Finish().

                              {

  //Int_t nofStations = fDigiPar->GetNStations(); //[R.K. 01/2017] unused variable?
}

PndGemDigiAna::ClassDef ( PndGemDigiAna  , 1    )

private

void PndGemDigiAna::CreateHistos ( )

private

Definition at line 192 of file PndGemDigiAna.cxx.

References fDigiPar, fHistoArray, fHistoList, fVerbose, PndGemStation::GetNSensors(), PndGemDigiPar::GetNStations(), PndGemStation::GetSensor(), PndGemSensor::GetSideChannels(), PndGemDigiPar::GetStation(), and sensor.

Referenced by Init().

                                 {
  fHistoList = new TList();

  if ( fVerbose )
    cout << "ought to create histogram for " << fDigiPar->GetNStations() << " stations." << endl;

  Int_t nofStations = fDigiPar->GetNStations();
  Int_t nofHists = 0;

  fHistoArray = new TClonesArray("TH1F", 5000);

  for ( Int_t istat = 0 ; istat < nofStations ; istat++ ) {
    PndGemStation* station = (PndGemStation*)fDigiPar->GetStation(istat);

    if ( station->GetNSensors() != 2 ) {
      std::cout << "will have problems with different numbers of sensors per station (" << station->GetNSensors() << ")" << std::endl;
    }

    for ( Int_t isens = 0 ; isens < station->GetNSensors() ; isens++ ) {
      PndGemSensor* sensor = (PndGemSensor*)station->GetSensor(isens);

      for ( Int_t iside = 0 ; iside < 2 ; iside++ ) {

        new ((*fHistoArray)[nofHists++]) TH1F(Form("fhLDTS_s%d_s%d_s%d",istat,isens,iside),
                                              Form("Last digi time stamp, station %d, sensor %d, %s side",istat,isens,(iside?"back":"front")),
                                              sensor->GetSideChannels(iside)+1,-0.5,sensor->GetSideChannels(iside)+0.5);
      }
    }
  }

  for ( Int_t istat = 0 ; istat < nofStations ; istat++ ) {
    PndGemStation* station = (PndGemStation*)fDigiPar->GetStation(istat);

    for ( Int_t isens = 0 ; isens < station->GetNSensors() ; isens++ ) {
      //PndGemSensor* sensor = (PndGemSensor*)station->GetSensor(isens); //[R.K. 01/2017] unused variable?

      for ( Int_t iside = 0 ; iside < 2 ; iside++ ) {

        new ((*fHistoArray)[nofHists++]) TH1F(Form("fhTimeBetweenDigis_s%d_s%d_s%d",istat,isens,iside),
                                              Form("Time between digis, station %d, sensor %d, %s side",istat,isens,(iside?"back":"front")),
                                              100000,0,100000);
      }
    }
  }

  if ( fVerbose ) {
    cout << "HM, seems that fHistoArray is filled with histograms, and there are " << nofHists << " of them" << endl;
  }

}

void PndGemDigiAna::Exec ( Option_t *  opt )

virtual

Execution

Definition at line 142 of file PndGemDigiAna.cxx.

References Double_t, fDigiPar, fGemDigiArray, fHistoArray, fNofEvents, fVerbose, PndGemDigi::GetChannelNr(), PndGemDigiPar::GetNStations(), PndGemDigi::GetSensorNr(), PndGemDigi::GetSide(), and PndGemDigi::GetStationNr().

                                  {
  
  if ( fVerbose ) {
    cout << "IN EVENT GOT " << fGemDigiArray->GetEntries() << " hits." << endl;
    cout << " - - - - - > " << fHistoArray->GetEntries() << " histograms" << endl;
  }

  Int_t nofGemDigis = fGemDigiArray->GetEntries();

  PndGemDigi* gemDigi = NULL;

  for ( Int_t idigi = 0 ; idigi < nofGemDigis ; idigi++ ) {
    gemDigi = (PndGemDigi*)fGemDigiArray->At(idigi);

    Int_t digiCode = 
      4*(gemDigi->GetStationNr()-1)+
      2*(gemDigi->GetSensorNr()-1)+
      gemDigi->GetSide();

    if ( digiCode < 0 || digiCode > 11 ) 
      cout << "oops, wrong digiCode " << digiCode
           << " from station " << gemDigi->GetStationNr()
           << " sensor " << gemDigi->GetSensorNr()
           << " side " << gemDigi->GetSide()
           << endl;

    Double_t lastTime = ((TH1F*)fHistoArray->At(digiCode))->GetBinContent(gemDigi->GetChannelNr());

    if ( lastTime > 0. ) {
      // cout << " DIGI on station " << gemDigi->GetStationNr()
      //           << " sensor " << gemDigi->GetSensorNr()
      //           << " side " << gemDigi->GetSide()
      //           << " channel " << gemDigi->GetChannelNr()
      //           << " came again at " << gemDigi->GetTimeStamp()
      //           << " after " << lastTime 
      //           << endl;
      ((TH1F*)fHistoArray->At(fDigiPar->GetNStations()*4+digiCode))->Fill(gemDigi->GetTimeStamp()-lastTime);

    }

    ((TH1F*)fHistoArray->At(digiCode))->SetBinContent(gemDigi->GetChannelNr(),gemDigi->GetTimeStamp());
  }

  fNofEvents++;

}

void PndGemDigiAna::Finish ( )

privatevirtual

Finish

Definition at line 252 of file PndGemDigiAna.cxx.

References AnaHistos(), fHistoArray, fHistoList, fNofEvents, and next.

                           {
  cout << "-------------------- PndGemDigiAna : Summary ------------------" << endl;
  cout << " Events:        " << setw(10) << fNofEvents << endl;
  cout << "---------------------------------------------------------------------" << endl; 

  AnaHistos();

  TFile* temp = gFile;
  FairRootManager* ioman = FairRootManager::Instance();
  gFile = ioman->GetOutFile();
  gDirectory = (TDirectory*)gFile;

  gDirectory->mkdir("GemDigiAna");
  gDirectory->cd("GemDigiAna");
  TIter next(fHistoList);
  for ( Int_t ihist = 0 ; ihist < fHistoArray->GetEntries() ; ihist++ ) {
    fHistoArray->At(ihist)->Write();
  }
  gDirectory->cd("..");

  gFile = temp;
}

InitStatus PndGemDigiAna::Init ( )

privatevirtual

Intialisation

Definition at line 72 of file PndGemDigiAna.cxx.

References CreateHistos(), fDigiPar, fGemDigiArray, PndGemDigiPar::GetNStations(), and rtdb.

                               {
  
  // Get and check FairRootManager
  FairRootManager* ioman = FairRootManager::Instance();
  if( !ioman ) {
    cout << "-E- "<< GetName() <<"::Init: "
         << "RootManager not instantised!" << endl;
    return kERROR;
  }
  
  // Get the pointer to the singleton FairRunAna object
  FairRunAna* ana = FairRunAna::Instance();
  if(NULL == ana) {
    cout << "-E- "<< GetName() <<"::Init :"
         <<" no FairRunAna object!" << endl;
    return kERROR;
  }
  // Get the pointer to run-time data base
  FairRuntimeDb* rtdb = ana->GetRuntimeDb();
  if(NULL == rtdb) {
    cout << "-E- "<< GetName() <<"::Init :"
         <<" no runtime database!" << endl;
    return kERROR;
  }
  
  // Get Gem digi Array
  fGemDigiArray  = (TClonesArray*) ioman->GetObject("GEMDigi");
  if ( !fGemDigiArray ) {
    cout << "-E- " << GetName() << "::Init: No PndGemDigi array!" << endl;
    return kERROR;
  }

  // Get GEM digitisation parameter container
  fDigiPar = (PndGemDigiPar*)(rtdb->getContainer("PndGemDetectors"));

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

  CreateHistos();

  return kSUCCESS;
}

InitStatus PndGemDigiAna::ReInit ( )

privatevirtual

Reinitialisation

Definition at line 117 of file PndGemDigiAna.cxx.

                                 {

  return kSUCCESS;
}

void PndGemDigiAna::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 124 of file PndGemDigiAna.cxx.

References fDigiPar, fNofEvents, 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"));

  fNofEvents = 0;
}

void PndGemDigiAna::SetVerbose ( const Int_t &  verbose )

inline

Public modifiers

Definition at line 51 of file PndGemDigiAna.h.

References fVerbose, and verbose.

{ fVerbose = verbose; };

Member Data Documentation

PndGemDigiPar* PndGemDigiAna::fDigiPar

private

Definition at line 51 of file PndGemDigiAna.h.

Referenced by CreateHistos(), Exec(), Init(), and SetParContainers().

TClonesArray* PndGemDigiAna::fGemDigiArray

private

Definition at line 57 of file PndGemDigiAna.h.

Referenced by Exec(), and Init().

TClonesArray* PndGemDigiAna::fHistoArray

private

Definition at line 65 of file PndGemDigiAna.h.

Referenced by CreateHistos(), Exec(), and Finish().

TList* PndGemDigiAna::fHistoList

private

Definition at line 62 of file PndGemDigiAna.h.

Referenced by CreateHistos(), and Finish().

Int_t PndGemDigiAna::fNofEvents

private

event counter

Event counter

Definition at line 60 of file PndGemDigiAna.h.

Referenced by Exec(), Finish(), and SetParContainers().

---

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

• PndGemDigiAna.h
• PndGemDigiAna.cxx