PndMvdConvertApv Class Reference

Convert meassured data into PndSdsDigiStrip. More...

#include <PndMvdConvertApv.h>

Public Member Functions
	PndMvdConvertApv ()
	PndMvdConvertApv (const TString &CalibFileName, const TString &HitFileName)
virtual	~PndMvdConvertApv ()
	PndMvdConvertApv (const PndMvdConvertApv &)=delete
PndMvdConvertApv &	operator= (const PndMvdConvertApv &)=delete
long int	GetNofEvents ()
std::vector< PndSdsDigiStrip >	ReadAll ()
std::vector< PndSdsDigiStrip >	ReadNext ()
Bool_t	Init ()
void	SetFakePair (Int_t TopModuleID, Int_t BottomModuleID)

Private Member Functions
void	LoadCalibration (TString CalibFileName)
std::vector< PndSdsDigiStrip >	Calc (std::vector< PndMvdApvHit > hitlist)
	ClassDef (PndMvdConvertApv, 1)
	Gives Access to the Path info of a hit. More...

Private Attributes
std::map< Int_t, std::map < Int_t, double > >	fCalibPars
	Calib Parameters: <moduleID<FE<channel> > > More...
long int	fNofEvents
	stored last eventID of the hitfile More...
long int	fEvent
	stored current eventID More...
long int	fLastEvent
	last eventID More...
bool	fNoCalib
	true if calib was succesfull loaded More...
TString	fHitFileName
	name of hitfile More...
std::ifstream	fDataFile
	fstream of hitfile More...
std::vector< PndMvdApvHit >	fhitlist
	stored readed events More...
Int_t	fTopModuleID
	stored fake module top side More...
Int_t	fBottomModuleID
	stored fake module bottom side More...
Bool_t	fFake
	knows if fake is allowed More...
TFile *	f
TTree *	t
PndMvdTsEvent *	tsEv
TClonesArray *	arr
PndGeoHandling *	fGeoH

Detailed Description

Convert meassured data into PndSdsDigiStrip.

PndMvdConvertApv.h

Author

L.Ackermann lars..nosp@m.acke.nosp@m.rmann.nosp@m.@phy.nosp@m.sik.t.nosp@m.u-dr.nosp@m.esden.nosp@m..de

Load calibration parameter and keeps strean to meassured data. Event wise read in and converting to PndSdsDigiStrip

Author

Lars Ackermann

Date

11.03.2009

Definition at line 39 of file PndMvdConvertApv.h.

Constructor & Destructor Documentation

PndMvdConvertApv::PndMvdConvertApv ( )

inline

default constructor

Definition at line 43 of file PndMvdConvertApv.h.

                           :
    fCalibPars(),
          fNofEvents(0),
    fEvent(-1),
          fLastEvent(0),
          fNoCalib(kFALSE),
          fHitFileName(""),
          fDataFile(),
          fhitlist(),
          fTopModuleID(0),
          fBottomModuleID(0),
          fFake(kFALSE),
          f(NULL),
          t(NULL),
          tsEv(NULL),
          arr(NULL),
    fGeoH(NULL)
  {};

PndMvdConvertApv::PndMvdConvertApv	(	const TString &	CalibFileName,
		const TString &	HitFileName
	)

main constructor, call all function to be ready for converting hits from hitfile

Parameters

CalibFileName	name of file where the calibration is stored
HitFileName	name of hitfile

Definition at line 20 of file PndMvdConvertApv.cxx.

References arr, f, fGeoH, fNofEvents, PndGeoHandling::Instance(), LoadCalibration(), t, and tsEv.

                                                                                          :
    fCalibPars(),
          fNofEvents(0),
    fEvent(-1),
          fLastEvent(0),
          fNoCalib(kFALSE),
          fHitFileName(""),
          fDataFile(),
          fhitlist(),
          fTopModuleID(0),
          fBottomModuleID(0),
          fFake(kFALSE),
          f(NULL),
          t(NULL),
          tsEv(NULL),
          arr(NULL),
    fGeoH(NULL)
{

  f = new TFile(HitFileName);
  fGeoH = PndGeoHandling::Instance();
  t = (TTree*) f->Get("T");
  tsEv = new PndMvdTsEvent();
  arr = new TClonesArray("PndMvdSiHit");

  cout << "---------------------------------------" << endl;
  cout << "Number of events: " << t->GetEntries() << endl;
  cout << "---------------------------------------" << endl;
        
  t->SetBranchAddress("events",&tsEv);
  cout << t->GetEntries() << " events in File" << endl;

  //LoadCalibration(CalibFileName, fes);
  LoadCalibration(CalibFileName);

  fNofEvents=t->GetEntries();
 
  cout<<"** end of PndMvdConvertApv::PndMvdConvertApv(const TString& , const TString&) **"<<endl;
}

virtual PndMvdConvertApv::~PndMvdConvertApv ( )

inlinevirtual

Destructor

Definition at line 70 of file PndMvdConvertApv.h.

References fDataFile.

        { fDataFile.close(); };

PndMvdConvertApv::PndMvdConvertApv ( const PndMvdConvertApv &  )

delete

Member Function Documentation

std::vector< PndSdsDigiStrip > PndMvdConvertApv::Calc ( std::vector< PndMvdApvHit >  hitlist )

private

Convert the readed hit and store them in the PndSdsDigiStrip

Parameters

hitlist

vector of PndMvdApvHit hits of one event from hitfile

Returns

list of converted PndSdsDigiStrip

Definition at line 114 of file PndMvdConvertApv.cxx.

References Double_t, fCalibPars, fFake, fNoCalib, and kMVDHitsStrip.

Referenced by ReadNext().

{
  std::vector<PndSdsDigiStrip> result;
  for(UInt_t hitnumber=0;hitnumber<hitlist.size();hitnumber++)
  {
    Double_t q=0.;
    if(fNoCalib)
    {
      q=1.*hitlist[hitnumber].GetADC();                                 // no calib adc -> e !!!
    }else{
      if (fCalibPars[hitlist[hitnumber].GetModuleID()].size())
      {
        if (fCalibPars[hitlist[hitnumber].GetModuleID()][hitlist[hitnumber].GetChannel()])
        {
          q=fCalibPars[hitlist[hitnumber].GetModuleID()][hitlist[hitnumber].GetChannel()]*(hitlist[hitnumber].GetADC())*1000.; // in electrons
        }
      }
    }
    
    //TString detPath="Module";
    //Int_t modId=-1; //[R.K.02/2017] Unused variable?
    //if(fFake) //[R.K.02/2017] Unused variable?
    //{ //[R.K.02/2017] Unused variable?
      //if(fTopModuleID==hitlist[hitnumber].GetModuleID() || fBottomModuleID==hitlist[hitnumber].GetModuleID())
        //modId = 1; //[R.K.02/2017] Unused variable?
      //if(fBottomModuleID==hitlist[hitnumber].GetModuleID()) hitlist[hitnumber].SetFeID(hitlist[hitnumber].GetFeID()+3); //[R.K.02/2017] Unused variable?
    //}else{ //[R.K.02/2017] Unused variable?
      //modId = hitlist[hitnumber].GetModuleID(); //[R.K.02/2017] Unused variable?
    //} //[R.K.02/2017] Unused variable?
   /* detPath+=modId;
    detPath+="Rect";
    //     std::cout<<detPath.Data()<<"   |    "<<modId<<std::endl;
    //TGeoVolume* Vol=gGeoManager->FindVolumeFast(detPath);
    TGeoVolume* Vol=0;
    if(Vol!=0) {
      //                std::cout<<Vol->GetName()<<std::endl;
      //                Vol->GetNode(-1)->cd();
      detPath="/SiliconTestStation_1/DummysensorAss_0/";
      detPath+=Vol->GetName();
      detPath+="_0";
      gGeoManager->cd(detPath.Data());
    }
    else {
      //                std::cout<<" -E- PndMvdConvertApv::Calc(): "<<detPath.Data()<<" does not exist"<<std::endl;
      detPath+="_nonexistent";
    }*/
    //     std::cout<<detPath.Data()<<std::endl;
    //     std::cout<<gGeoManager->GetPath()<<std::endl;
    //     if (0==Vol) std::cout<<"0"; 
    //     else std::cout<<"1";
    //     std::cout<<Vol->GetName();
    
    //     detPath = Vol->GetName();
    //  std::cout << "write Digi with "<< detPath.Data()<<" ( "<<fGeoH->GetID(detPath)<<" )"<<std::endl;
    if(fFake)
    {
      PndSdsDigiStrip DigiHit(hitlist[hitnumber].GetEventID(), // index
                              kMVDHitsStrip,                   // panda detID 
                              -1,                              // No  SensorID (from geopath) 
                              hitlist[hitnumber].GetFeID(),    // fe
                              hitlist[hitnumber].GetChannel(), // chan
                              q/*/1000/1000*/,                 // charge
                              hitlist[hitnumber].GetTimestamp()// timestamp
                              );
          DigiHit.SetTimeStampError(hitlist[hitnumber].GetTriggerID());                             
                              
      //      cout << "Ev. " << hitlist[hitnumber].GetEventID() << ", FE: " << hitlist[hitnumber].GetFeID() << ", ch: " << hitlist[hitnumber].GetChannel() << ", FAKE" << endl;

      result.push_back(DigiHit);
    }else{
      PndSdsDigiStrip DigiHit(hitlist[hitnumber].GetEventID(), 
                              kMVDHitsStrip,
                              hitlist[hitnumber].GetModuleID(),
                              hitlist[hitnumber].GetFeID(),
                              hitlist[hitnumber].GetChannel(),
                              q/*/1000/1000*/,
                              hitlist[hitnumber].GetTimestamp() 
                              );
      DigiHit.SetTimeStampError(hitlist[hitnumber].GetTriggerID());
      
                                                        
      //   cout << "Ev. " << hitlist[hitnumber].GetEventID() << ", FE: " << hitlist[hitnumber].GetChannel()/128 << ", ch: " << hitlist[hitnumber].GetChannel()%128 << ", sens" << hitlist[hitnumber].GetModuleID()  << endl;
      result.push_back(DigiHit);
    }
  }
  return result;
}

PndMvdConvertApv::ClassDef ( PndMvdConvertApv  , 1    )

private

Gives Access to the Path info of a hit.

long int PndMvdConvertApv::GetNofEvents

(

)

Returns

long integer of last eventID from hitfile

Definition at line 204 of file PndMvdConvertApv.cxx.

References fNofEvents.

{
  return fNofEvents;
}

Bool_t PndMvdConvertApv::Init

(

)

Initialize global geometry

Returns

success

Definition at line 60 of file PndMvdConvertApv.cxx.

Referenced by PndMvdConvertApvTask::Init().

{
  return kTRUE;
}

void PndMvdConvertApv::LoadCalibration ( TString  CalibFileName )

private

Read the calibration of the modules

Parameters

CalibFileName	Filename where the calibration of the modules are listed
modules	vector of moduleIDs of modules in the hitfile

Returns

void load calibration from database

Definition at line 69 of file PndMvdConvertApv.cxx.

References c, fCalibPars, and fNoCalib.

Referenced by PndMvdConvertApv().

{
        
  std::ifstream calibfile(CalibFileName);
  if(!calibfile)
  {
    cout<<"Calibration file not found"<<endl;
    fNoCalib=true;
    return;
  }
  while (!calibfile.eof())                                              // read data
  {
    char c;
    calibfile>>c;
    if (calibfile.eof()) break;
    if (!isdigit(c))    
    {
      char str[256];
      calibfile.getline(str,256);
      continue;
    }
    calibfile.putback(c);
    
    int boxID, channel;
    double value;
    
    calibfile >> boxID >> channel >> value;
    
    //for(unsigned int vec=0;vec<fes.size();vec++)
    //{
      //if(feID==fes[vec])
      //{
    fCalibPars[boxID][channel]=value;
      //}
    //}
  }
  calibfile.close();
  fNoCalib=false;
  cout<<"Calibration succesfully read"<<endl;
  return;
         
}

PndMvdConvertApv& PndMvdConvertApv::operator= ( const PndMvdConvertApv &  )

delete

std::vector< PndSdsDigiStrip > PndMvdConvertApv::ReadAll

(

)

read all events from hitfile

Returns

vector of PndSdsDigiStrip of the event

Definition at line 279 of file PndMvdConvertApv.cxx.

References fEvent, fNofEvents, i, and ReadNext().

{
  std::vector<PndSdsDigiStrip> result;
  while(fEvent!=fNofEvents)
  {
    std::vector<PndSdsDigiStrip> dummy=ReadNext();
    for(unsigned int i=0;i<dummy.size();++i) result.push_back(dummy[i]);
  }
  return result;
}

std::vector< PndSdsDigiStrip > PndMvdConvertApv::ReadNext

(

)

read the next event from hitfile

Returns

vector of PndSdsDigiStrip of the event

Definition at line 211 of file PndMvdConvertApv.cxx.

References arr, Calc(), ev, PndMvdSiHit::fAdc, PndMvdSiHit::fBox, PndMvdSiHit::fChannel, fe, fEvent, fhitlist, fLastEvent, fNofEvents, PndMvdSiHit::fNumFrames, PndMvdTsEvent::GetEventId(), PndMvdTsEvent::GetExtClockTimeStamp(), PndMvdTsEvent::GetSiHitList(), hit, t, and tsEv.

Referenced by PndMvdConvertApvTask::Exec(), for(), and ReadAll().

{
        
        
        std::vector<PndSdsDigiStrip> digiList;
        //bool work=true;
         
        //int triggID=0;
        int fe=0;
        //int ts=0;
        int frame=0;
        int ch=0;
        int l=0;
        int moduleID=0;
        double q=0.;
        long int ev=0;
        UInt_t ClockReset = 0.;
        ULong64_t ClockCounts = 0.;

        if (fEvent >= -1 && fEvent <= fNofEvents)
        {
        
                t->GetEvent(fEvent);
        
                ClockCounts = tsEv->GetExtClockTimeStamp(ClockReset); 
                // Storing ClockReset in fTriggerID
            // Storing ClockCounts in fTimestamp
        
                arr = tsEv->GetSiHitList();
                fhitlist.clear();
                
                for (Int_t kk = 0 ; kk < arr->GetEntries() ; kk++)
                {
         
                        PndMvdSiHit *hit = (PndMvdSiHit*) arr->At(kk);
                        ev = tsEv->GetEventId();
                        //fe = (Int_t) (hit->fChannel)/128;
                        //ch = (Int_t) (hit->fChannel)%128;
                        ch = (Int_t) (hit->fChannel);
                        q = hit->fAdc;
                        l = hit->fNumFrames;
                        moduleID = hit->fBox;
                
                        //if(fhitlist.size()>20) fhitlist.clear();
                        //                      digiList = Calc(fhitlist);
                        
                        fLastEvent=ev;
                        //PndMvdApvHit Hit(ev, moduleID, fe, triggID, ts, frame, ch, q, l);
                        PndMvdApvHit Hit(ev, moduleID, fe, ClockReset, ClockCounts, frame, ch, q, l);
                        //cout << "Ev. " << fEvent << ", trigg: " << triggID << ", sens: " << moduleID << ", ch: " << ch << ", fe: " << (Int_t) (hit->fChannel)/128 << ", channel: " << (Int_t) (hit->fChannel)%128 << endl;
                        fhitlist.push_back(Hit);
                }

                //if(fhitlist.size()>20) fhitlist.clear();
                digiList = Calc(fhitlist);
                fEvent++;
                        
        }
        
        if(!(fEvent%10000)) cout<<"[ "<<(fEvent*100)/fNofEvents<<" %] "<<fEvent<<" events converted ..."<<endl;
        if(fEvent==fNofEvents) cout<<"[100 %] "<<fEvent<<" events converted"<<endl;
                
        return digiList;
         
}

void PndMvdConvertApv::SetFakePair	(	Int_t	TopModuleID,
		Int_t	BottomModuleID
	)

Set two moduleIDs to merge them as one double sided sensor

Parameters

TopModuleID	moduleID for top side of fake sensor
BottomModuleID	moduleID for bottom side of fake sensor

Returns

void

Definition at line 290 of file PndMvdConvertApv.cxx.

References fBottomModuleID, fFake, and fTopModuleID.

{
  fFake=true;
  fTopModuleID=TopModuleID;
  fBottomModuleID=BottomModuleID;
  return;
}

Member Data Documentation

TClonesArray* PndMvdConvertApv::arr

private

Definition at line 184 of file PndMvdConvertApv.h.

Referenced by PndMvdConvertApv(), and ReadNext().

TFile* PndMvdConvertApv::f

private

Definition at line 180 of file PndMvdConvertApv.h.

Referenced by PndMvdConvertApv().

Int_t PndMvdConvertApv::fBottomModuleID

private

stored fake module bottom side

Definition at line 175 of file PndMvdConvertApv.h.

Referenced by SetFakePair().

std::map<Int_t, std::map<Int_t, double> > PndMvdConvertApv::fCalibPars

private

Calib Parameters: <moduleID<FE<channel> > >

Definition at line 148 of file PndMvdConvertApv.h.

Referenced by Calc(), and LoadCalibration().

std::ifstream PndMvdConvertApv::fDataFile

private

fstream of hitfile

Definition at line 166 of file PndMvdConvertApv.h.

Referenced by ~PndMvdConvertApv().

long int PndMvdConvertApv::fEvent

private

stored current eventID

Definition at line 154 of file PndMvdConvertApv.h.

Referenced by ReadAll(), and ReadNext().

Bool_t PndMvdConvertApv::fFake

private

knows if fake is allowed

Definition at line 178 of file PndMvdConvertApv.h.

Referenced by Calc(), and SetFakePair().

PndGeoHandling* PndMvdConvertApv::fGeoH

private

Definition at line 186 of file PndMvdConvertApv.h.

Referenced by PndMvdConvertApv().

TString PndMvdConvertApv::fHitFileName

private

name of hitfile

Definition at line 163 of file PndMvdConvertApv.h.

std::vector<PndMvdApvHit> PndMvdConvertApv::fhitlist

private

stored readed events

Definition at line 169 of file PndMvdConvertApv.h.

Referenced by ReadNext().

long int PndMvdConvertApv::fLastEvent

private

last eventID

Definition at line 157 of file PndMvdConvertApv.h.

Referenced by ReadNext().

bool PndMvdConvertApv::fNoCalib

private

true if calib was succesfull loaded

Definition at line 160 of file PndMvdConvertApv.h.

Referenced by Calc(), and LoadCalibration().

long int PndMvdConvertApv::fNofEvents

private

stored last eventID of the hitfile

Definition at line 151 of file PndMvdConvertApv.h.

Referenced by GetNofEvents(), PndMvdConvertApv(), ReadAll(), and ReadNext().

Int_t PndMvdConvertApv::fTopModuleID

private

stored fake module top side

Definition at line 172 of file PndMvdConvertApv.h.

Referenced by SetFakePair().

TTree* PndMvdConvertApv::t

private

Definition at line 181 of file PndMvdConvertApv.h.

Referenced by PndMvdConvertApv(), and ReadNext().

PndMvdTsEvent* PndMvdConvertApv::tsEv

private

Definition at line 182 of file PndMvdConvertApv.h.

Referenced by PndMvdConvertApv(), and ReadNext().

---

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

• PndMvdConvertApv.h
• PndMvdConvertApv.cxx