PndLmdDigiProducer Class Reference

#include <PndLmdDigiProducer.h>

Inheritance diagram for PndLmdDigiProducer:

Public Member Functions
	PndLmdDigiProducer ()
	PndLmdDigiProducer (Int_t k, SensorSegment segment=kLEFT)
	PndLmdDigiProducer (Double_t k, SensorSide side=kTOP)
	~PndLmdDigiProducer ()
virtual void	SetParContainers ()
virtual InitStatus	Init ()
virtual InitStatus	ReInit ()
virtual void	Exec (Option_t *opt)
void	AddDigi (Int_t &iStrip, Int_t iPoint, Int_t detID, TString detName, Int_t fe, Int_t channel, Double_t charge)
void	SetParamSet ()

Private Member Functions
void	Register ()
void	ReSet ()
void	ProduceHits ()
Bool_t	SelectSensorParams (TString detname)
	ClassDef (PndLmdDigiProducer, 5)

Private Attributes
PndLmdDigiPara *	fDigiParTrap
PndLmdDigiPara *	fCurrentDigiPar
PndLmdDigiStrips *	temp
PndLmdCalStrip *	fStripCalTrapLeft
PndLmdCalStrip *	fStripCalTrapRight
PndLmdCalStrip *	fCurrentStripCalLeft
PndLmdCalStrip *	fCurrentStripCalRight
PndLmdCalStrip *	fStripCalTrapCircle
PndLmdCalStrip *	fCurrentStripCalCircle
SensorSegment	fSegment
SensorSide	fSide
TString	fBranchName
TClonesArray *	fPointArray
TClonesArray *	fStripArray
PndLmdGeoHandling *	fGeoH

Detailed Description

Definition at line 32 of file PndLmdDigiProducer.h.

Constructor & Destructor Documentation

PndLmdDigiProducer::PndLmdDigiProducer

(

)

Definition at line 35 of file PndLmdDigiProducer.cxx.

References fBranchName.

                                       : FairTask("Pnd LMD Digi Producer")
{
  fBranchName = "PndLmdMCPoint";
}

PndLmdDigiProducer::PndLmdDigiProducer	(	Int_t	k,
		SensorSegment	segment = kLEFT
	)

Definition at line 40 of file PndLmdDigiProducer.cxx.

References fBranchName, and fSegment.

                    :FairTask("Pnd LMD Digi Producer")
{
        fBranchName = "PndLmdMCPoint";
//      fDigiParTrap = digipar;
        fSegment = segment;
        }

PndLmdDigiProducer::PndLmdDigiProducer	(	Double_t	k,
		SensorSide	side = kTOP
	)

Definition at line 48 of file PndLmdDigiProducer.cxx.

References fBranchName, and fSide.

                   :FairTask("Pnd LMD Digi Producer")
 {
        fBranchName = "PndLmdMCPoint";
  //  fDigiParTrap = digipar;
    fSide = side;
  }

PndLmdDigiProducer::~PndLmdDigiProducer

(

)

Definition at line 57 of file PndLmdDigiProducer.cxx.

References fGeoH.

{
        delete fGeoH;
        }

Member Function Documentation

void PndLmdDigiProducer::AddDigi	(	Int_t &	iStrip,
		Int_t	iPoint,
		Int_t	detID,
		TString	detName,
		Int_t	fe,
		Int_t	channel,
		Double_t	charge
	)

Definition at line 241 of file PndLmdDigiProducer.cxx.

References Bool_t, fe, fStripArray, PndLmdDigiStrips::GetChannel(), PndLmdDigiStrips::GetDetID(), PndLmdDigiStrips::GetDetName(), and PndLmdDigiStrips::GetFe().

Referenced by Exec().

{
    Bool_t found = kFALSE;
    PndLmdDigiStrips* aDigi= 0;
    for(Int_t skr=0;skr<iStrip && found==kFALSE;skr++){
        aDigi = (PndLmdDigiStrips*)fStripArray->At(skr);
        if(aDigi->GetDetID()==detID,
           aDigi->GetDetName()==detName,
           aDigi->GetFe()==fe,
           aDigi->GetChannel()==channel)
        {
        //      aDigi->AddCharge(charge);
        //      aDigi->AddIndex(iPoint);
                found=kTRUE;
        }
    }

    if(found == kFALSE){
        new((*fStripArray)[iStrip]) PndLmdDigiStrips(iPoint,detID,detName,fe,channel,charge);
        iStrip++;
    }
}

PndLmdDigiProducer::ClassDef ( PndLmdDigiProducer  , 5    )

private

void PndLmdDigiProducer::Exec ( Option_t *  opt )

virtual

Definition at line 127 of file PndLmdDigiProducer.cxx.

References AddDigi(), PndLmdCalStrip::CalChannelFromStrip(), PndLmdCalStrip::CalFeFromStrip(), PndLmdCalStrip::DigiStripFromPoint(), Double_t, fCurrentDigiPar, fCurrentStripCalLeft, fCurrentStripCalRight, fDigiParTrap, fGeoH, fPointArray, fStripArray, fStripCalTrapLeft, fStripCalTrapRight, fVerbose, PndLmdDigiPara::GetNrFeChannels(), PndLmdDigiPara::GetNrLeftFe(), point, and pos.

{
        //reset the output array
        fStripArray->Clear();

        //define variables
        PndLmdMCPoint* point = NULL;

        Int_t nPoint = fPointArray->GetEntriesFast();
        if(fVerbose>0){
                std::cout<<"the Nr of points is"<<nPoints<<std::endl;
        }

        //loop all points in an event
        Int_t iStrip = 0;
        for(Int_t iPoint = 0;iPoint<nPoints;iPoint++)
        {
                point = (PndLmdMCPoint*)fPointArray->At(iPoint);
                if(! point){
                        std::cout<<"No point !"<<std::endl;
                        continue;
                }

                fCurrentStripCalLeft = fStripCalTrapLeft;
                fCurrentStripCalRight = fStripCalTrapRight;
                fCurrentDigiPar = fDigiParTrap;

                //transform the point from global coordinate to local coordinate
                TVector3 posInL = fGeoH->MasterToLocalId(point->GetEntryPoint(),point->GetDetName());
                TVector3 posOutL = fGeoH->MasterToLocalId(point->GetExitPoint(),point->GetDetName());
                //std::cout<<"the point in global coordinate is "<<point->GetEntryPoint().X()<<","<<point->GetEntryPoint().Y()
                //<<","<<point->GetEntryPoint().Z()<<std::endl;
                std::cout<<"the Entry point in local coordinate is "<<posInL.X()<<","<<posInL.Y()<<","<<posInL.Z()<<std::endl;
                std::cout<<"the Exit point in local coordinate is "<<posOutL.X()<<","<<posOutL.Y()<<","<<posOutL.Z()<<std::endl;

                Double_t eLoss = point->GetEnergyLoss();
                std::cout<<"the eLoss of this point is: "<<eLoss<<endl;

                TVector3 pos,pos1 = posInL + posOutL;
                pos.SetXYZ(pos1.X()/2.,pos1.Y()/2.,pos1.Z()/2.);
                //for Left segment
                if(pos.X()>0)
                {
                //      std::cout<<"Hits on Left segment"<<std::endl;

                //calculate the cluster strips fired
                std::vector<PndLmdStrips> leftStrips =
                        fCurrentStripCalLeft->DigiStripFromPoint(posInL.X(),posInL.Y(),posInL.Z(),
                                                                        posOutL.X(),posOutL.Y(),posOutL.Z(),
                                                                        eLoss);
                      if(leftStrips.size()!=0){
                          for(std::vector<PndLmdStrips>::const_iterator kit=leftStrips.begin();kit!=leftStrips.end();++kit)
                          {//register the strips fired
                                  AddDigi(iStrip,iPoint,kLMDHitsStrip,point->GetDetName(),
                                          fCurrentStripCalLeft->CalFeFromStrip(kit->GetIndex()),
                                          fCurrentStripCalLeft->CalChannelFromStrip(kit->GetIndex()),
                                          kit->GetCharge());
                                  Int_t chan1 = (fCurrentStripCalLeft->CalFeFromStrip(kit->GetIndex()))*(fCurrentDigiPar->GetNrFeChannels())
                                    +( fCurrentStripCalLeft->CalChannelFromStrip(kit->GetIndex()));
                                  
                                  cout<<"the iPoint is: "<<iPoint<<endl;
                                  cout<<"the iStrip is :"<<iStrip-1<<endl;
                                  cout<<"the channel is :"<<chan1<<endl;
                                  cout<<"the information of digis: "<<point->GetDetName()<<","<<kit->GetCharge()<<endl;
                                  cout<<endl;
                                  
                          }
                      }else std::cout<<"the left segment is empty"<<std::endl;

                }else{// std::cout<<"Hits on Right segment"<<std::endl;
              std::vector<PndLmdStrips> rightStrips =
                  fCurrentStripCalRight->DigiStripFromPoint(posInL.X(),posInL.Y(),posInL.Z(),
                                                                        posOutL.X(),posOutL.Y(),posOutL.Z(),
                                                                        eLoss);
        //      cout<<"Call the calStrip class to get digi successfully"<<endl;

              if(rightStrips.size()!=0){
                  for(std::vector<PndLmdStrips>::const_iterator kit=rightStrips.begin();
                       kit!=rightStrips.end();++kit){
                            //register the strips fired
                            AddDigi(iStrip,iPoint,kLMDHitsStrip,point->GetDetName(),
                                        fCurrentStripCalRight->CalFeFromStrip(kit->GetIndex())+fCurrentDigiPar->GetNrLeftFe(),
                                        fCurrentStripCalRight->CalChannelFromStrip(kit->GetIndex()),
                                        kit->GetCharge());
                         Int_t chan2 = (fCurrentStripCalRight->CalFeFromStrip(kit->GetIndex())+fCurrentDigiPar->GetNrLeftFe())*(fCurrentDigiPar->GetNrFeChannels())
                                       +(fCurrentStripCalLeft->CalChannelFromStrip(kit->GetIndex()));

/*
                            new((*fStripArray)[iStrip]) PndLmdDigiStrips(iPoint,kLMDHitsStrip,point->GetDetName(),
                                        fCurrentStripCalRight->CalFeFromStrip(kit->GetIndex())+fCurrentDigiPar->GetNrLeftFe(),
                                        fCurrentStripCalRight->CalChannelFromStrip(kit->GetIndex()),
                                        kit->GetCharge());
                                        iStrip++;
*/
                        //    cout<<endl;
                            cout<<"the iPoint is: "<<iPoint<<endl;
                            cout<<"the iStrip is :"<<iStrip-1<<endl;
                            cout<<"the channel is :"<<chan2<<endl;
                            cout<<"the information of digis: "<<point->GetDetName()<<","<<kit->GetCharge()<<endl;
                            cout<<endl;

                      }
              }else std::cout<<"the right segment is empty"<<std::endl;

                     }//judge the hits on left segment or right segment, then calculate strips

        }//loop all points for one event

 //   if(fVerbose>0)
        std::cout<<"-I- PndLmdDigiProducer: "<<nPoints<<" PndLmdMCPoint, "
                    <<iStrip<<" PndLmdDigis"<<std::endl<<std::endl;

}//public method for Exec

InitStatus PndLmdDigiProducer::Init ( )

virtual

Definition at line 79 of file PndLmdDigiProducer.cxx.

References fBranchName, fDigiParTrap, fGeoH, fPointArray, fStripArray, fStripCalTrapLeft, fStripCalTrapRight, gGeoManager, PndLmdDigiPara::print(), and PndLmdCalStrip::print().

{
    FairRun* ana = FairRun::Instance();
    FairRootManager* ioman = FairRootManager::Instance();

    fGeoH = new PndLmdGeoHandling(gGeoManager);

    if(! ioman)
    {
      std::cout<<"-E- PndLmdDigiProducer::Init :"
      <<" The FariRootManager is not initialized!"<<std::endl;
      return kFATAL;
    }

    //input array initialization
    fPointArray = (TClonesArray*)ioman->GetObject(fBranchName);
    if(! fPointArray)
    {
        std::cout<<"-E- PndLmdDigiProducer::Init :"
        <<"No PndLmdMCPoint Array"<<std::endl;
     return kERROR;
    }

    //creat and register output array
    fStripArray = new TClonesArray("PndLmdDigiStrips");
    ioman->Register("LMDStripDigis","LMD",fStripArray,kTRUE);

    if(! fDigiParTrap)
    {
       std::cout<<"-E- PndLmdDigiProducer::DigiParTrap Container doesn't exist!"<<std::endl;
       return kERROR;
    }
//print the digipar parameters for check
    fDigiParTrap->print();

    fStripCalTrapLeft = new PndLmdCalStrip(fDigiParTrap,kLEFT);
    cout<<"the digi paras for left segment!"<<endl;
    fStripCalTrapLeft->print();

    fStripCalTrapRight = new PndLmdCalStrip(fDigiParTrap,kRIGHT);
    cout<<"the digi paras for right segment!"<<endl;
    fStripCalTrapRight->print();

    std::cout<<"the PndLmdDigiProcuder initialized successfully!"<<std::endl;

    return kSUCCESS;
}

void PndLmdDigiProducer::ProduceHits ( )

private

void PndLmdDigiProducer::Register ( )

private

InitStatus PndLmdDigiProducer::ReInit ( )

virtual

Definition at line 72 of file PndLmdDigiProducer.cxx.

References SetParContainers().

{
   SetParContainers();
   return kSUCCESS;
}

void PndLmdDigiProducer::ReSet ( )

private

Bool_t PndLmdDigiProducer::SelectSensorParams ( TString  detname )

private

void PndLmdDigiProducer::SetParamSet

(

)

void PndLmdDigiProducer::SetParContainers ( )

virtual

Definition at line 62 of file PndLmdDigiProducer.cxx.

References fDigiParTrap, and rtdb.

Referenced by ReInit().

{
        FairRun* ana = FairRun::Instance();
        FairRuntimeDb* rtdb = ana->GetRuntimeDb();
        fDigiParTrap = (PndLmdDigiPara*)(rtdb->getContainer("LmdStripDigiParaTrap"));
        //if(!fDigiParTrap) cout<<"The digi parameters has been set to container!"<<endl;
        //fDigiParTrap->print();
        }

Member Data Documentation

TString PndLmdDigiProducer::fBranchName

private

Definition at line 76 of file PndLmdDigiProducer.h.

Referenced by Init(), and PndLmdDigiProducer().

PndLmdDigiPara* PndLmdDigiProducer::fCurrentDigiPar

private

Definition at line 60 of file PndLmdDigiProducer.h.

Referenced by Exec().

PndLmdCalStrip* PndLmdDigiProducer::fCurrentStripCalCircle

private

Definition at line 71 of file PndLmdDigiProducer.h.

PndLmdCalStrip* PndLmdDigiProducer::fCurrentStripCalLeft

private

Definition at line 67 of file PndLmdDigiProducer.h.

Referenced by Exec().

PndLmdCalStrip* PndLmdDigiProducer::fCurrentStripCalRight

private

Definition at line 68 of file PndLmdDigiProducer.h.

Referenced by Exec().

PndLmdDigiPara* PndLmdDigiProducer::fDigiParTrap

private

Definition at line 59 of file PndLmdDigiProducer.h.

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

PndLmdGeoHandling* PndLmdDigiProducer::fGeoH

private

Definition at line 89 of file PndLmdDigiProducer.h.

Referenced by Exec(), Init(), and ~PndLmdDigiProducer().

TClonesArray* PndLmdDigiProducer::fPointArray

private

Definition at line 79 of file PndLmdDigiProducer.h.

Referenced by Exec(), and Init().

SensorSegment PndLmdDigiProducer::fSegment

private

Definition at line 73 of file PndLmdDigiProducer.h.

Referenced by PndLmdDigiProducer().

SensorSide PndLmdDigiProducer::fSide

private

Definition at line 74 of file PndLmdDigiProducer.h.

Referenced by PndLmdDigiProducer().

TClonesArray* PndLmdDigiProducer::fStripArray

private

Definition at line 81 of file PndLmdDigiProducer.h.

Referenced by AddDigi(), Exec(), and Init().

PndLmdCalStrip* PndLmdDigiProducer::fStripCalTrapCircle

private

Definition at line 70 of file PndLmdDigiProducer.h.

PndLmdCalStrip* PndLmdDigiProducer::fStripCalTrapLeft

private

Definition at line 65 of file PndLmdDigiProducer.h.

Referenced by Exec(), and Init().

PndLmdCalStrip* PndLmdDigiProducer::fStripCalTrapRight

private

Definition at line 66 of file PndLmdDigiProducer.h.

Referenced by Exec(), and Init().

PndLmdDigiStrips* PndLmdDigiProducer::temp

private

Definition at line 62 of file PndLmdDigiProducer.h.

---

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

• PndLmdDigiProducer.h
• PndLmdDigiProducer.cxx