PndEmcHitProducer Class Reference

creates PndEmcHits from PndEmcPoints More...

#include <PndEmcHitProducer.h>

Inheritance diagram for PndEmcHitProducer:

Public Types
typedef std::map< Int_t, Float_t >	mapper

Public Member Functions
	PndEmcHitProducer ()
	PndEmcHitProducer (Bool_t val)
	~PndEmcHitProducer ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
PndEmcHit *	AddHit (Int_t trackID, Int_t detID, Float_t energy, Float_t time, std::vector< Int_t > &mctruth, FairMultiLinkedData entering, FairMultiLinkedData exiting)
void	SetStorageOfData (Bool_t val)
void	SetNonuniformityFile (const char *filename)
void	SetDayOne (bool d=true)
void	FinishTask ()
void	SetPersistency (Bool_t val=kTRUE)
Bool_t	GetPersistency ()

Private Member Functions
virtual void	SetParContainers ()
void	cleansortmclist (std::vector< Int_t > &newlist, TClonesArray *mcTrackArray)
	PndEmcHitProducer (const PndEmcHitProducer &L)
PndEmcHitProducer &	operator= (const PndEmcHitProducer &)
	ClassDef (PndEmcHitProducer, 1)
bool	AcceptDayOne (PndEmcPoint *p)

Private Attributes
Int_t	fUse_nonuniformity
TString	fNonuniformityFile
TClonesArray *	fPointArray
TClonesArray *	fMCTrackArray
TClonesArray *	fHitArray
TObjArray *	fVolumeArray
Int_t	fMapVersion
Float_t	fEnergyThreshold
mapper	emcX
mapper	emcY
mapper	emcZ
PndEmcStructure *	fEmcStr
PndEmcMapper *	fMapper
PndEmcDigiPar *	fDigiPar
PndEmcGeoPar *	fGeoPar
PndEmcDigiNonuniformityPar *	fNonuniformityPar
Bool_t	fStoreHits
map< Int_t, Float_t >	fTrackEnergy
map< Int_t, Float_t >	fTrackTime
map< Int_t, std::vector< Int_t > >	fTrackMcTruth
map< Int_t, std::vector< Int_t > >	fPointMatch
map< Int_t, FairMultiLinkedData >	fTrackEntering
map< Int_t, FairMultiLinkedData >	fTrackExiting
bool	fDayOne

Detailed Description

creates PndEmcHits from PndEmcPoints

Definition at line 37 of file PndEmcHitProducer.h.

Member Typedef Documentation

typedef std::map<Int_t, Float_t> PndEmcHitProducer::mapper

Definition at line 41 of file PndEmcHitProducer.h.

Constructor & Destructor Documentation

PndEmcHitProducer::PndEmcHitProducer

(

)

Default constructor

Definition at line 44 of file PndEmcHitProducer.cxx.

References fNonuniformityFile, and PndPersistencyTask::SetPersistency().

                                     :
        PndPersistencyTask("Ideal EMC hit Producer"),
        fUse_nonuniformity(0), fNonuniformityFile(""), fPointArray(), fMCTrackArray(), fHitArray(), fVolumeArray(), fMapVersion(0), fEnergyThreshold(0), emcX(), emcY(), emcZ(), fEmcStr(), fMapper(), fDigiPar(), fGeoPar(), fNonuniformityPar(), fDayOne(false)
{
        fNonuniformityFile=gSystem->Getenv("VMCWORKDIR");
        fNonuniformityFile+="/input/EmcDigiNoniformityPars.root";
        SetPersistency(kTRUE);
}

PndEmcHitProducer::PndEmcHitProducer

(

Bool_t

val

)

Definition at line 54 of file PndEmcHitProducer.cxx.

References fNonuniformityFile, and PndPersistencyTask::SetPersistency().

                                               :
        PndPersistencyTask("Ideal EMC hit Producer"),
        fUse_nonuniformity(0), fNonuniformityFile(""), fPointArray(), fMCTrackArray(), fHitArray(), fVolumeArray(), fMapVersion(0), fEnergyThreshold(0), emcX(), emcY(), emcZ(), fEmcStr(), fMapper(), fDigiPar(), fGeoPar(), fNonuniformityPar(), fDayOne(false)
{
        fNonuniformityFile=gSystem->Getenv("VMCWORKDIR");
        fNonuniformityFile+="/input/EmcDigiNoniformityPars.root";
        SetPersistency(val);
}

PndEmcHitProducer::~PndEmcHitProducer

(

)

Destructor

Definition at line 64 of file PndEmcHitProducer.cxx.

References fEmcStr.

{ delete fEmcStr;}

PndEmcHitProducer::PndEmcHitProducer ( const PndEmcHitProducer &  L )

private

Member Function Documentation

bool PndEmcHitProducer::AcceptDayOne ( PndEmcPoint *  p )

private

Definition at line 424 of file PndEmcHitProducer.cxx.

References fDayOne, PndEmcPoint::GetPhi(), and phi.

Referenced by Exec().

{
  if(!fDayOne) return true;
  if(p->GetDetectorID()>250e6) return true; // afaik this is convention for the barrel
  float phi=p->GetPhi();
  if(abs(phi-90)<22.5) return false;
  if(abs(phi-270)<22.5) return false;
  if(abs(phi+90)<22.5) return false;
  return true;
}

PndEmcHit * PndEmcHitProducer::AddHit	(	Int_t	trackID,
		Int_t	detID,
		Float_t	energy,
		Float_t	time,
		std::vector< Int_t > &	mctruth,
		FairMultiLinkedData	entering,
		FairMultiLinkedData	exiting
	)

Definition at line 389 of file PndEmcHitProducer.cxx.

References emcX, emcY, emcZ, fHitArray, and hit.

Referenced by Exec().

{
        // It fills the PndEmcHit category

        //cout << "PndEmcHitProducer: track " << trackID << " evt " << eventID
        //<< " sec " << sec << " plane " << pla << " strip " << strip << "box
        //" << box << " tube " << tub << endl;
        TClonesArray& clref = *fHitArray;
        Int_t size = clref.GetEntriesFast();
        PndEmcHit* hit = new(clref[size]) PndEmcHit(trackID, detID, energy, time, emcX[detID],
                        emcY[detID], emcZ[detID], mctruth, entering, exiting);
        //hit->Print();
        return  hit;
}

PndEmcHitProducer::ClassDef ( PndEmcHitProducer  , 1    )

private

void PndEmcHitProducer::cleansortmclist ( std::vector< Int_t > &  newlist, TClonesArray *  mcTrackArray  )

private

Definition at line 169 of file PndEmcHitProducer.cxx.

References PndMCTrack::GetMotherID(), PndMCTrack::GetStartVertex(), gGeoManager, pt(), and TString.

Referenced by Exec().

{
        std::vector <Int_t> tmplist;
        std::vector <Int_t> tmplist2;
        // Sort list...
        std::sort( newlist.begin(), newlist.end());
        // and copy every id only once (even so it might be in the list several times)
        std::unique_copy( newlist.begin(), newlist.end(), std::back_inserter( tmplist ) );

        // Now check if mother or (grand)^x-mother are already in the list
        // (which means i am a secondary)... if so, remove myself
        for(Int_t j=tmplist.size()-1; j>=0; j--){
                bool flag = false;
                PndMCTrack *pt;
                Int_t id = tmplist[j];
                // if -1 index put it in the list and continue
                if(id < 0) {
                        tmplist2.push_back(id);
                        continue;
                }
                while(!flag){
                        pt=((PndMCTrack*)mcTrackArray->At(id));
                        // If the particle is primary store it and stop
                        if(pt->GetMotherID()<0) {
                                tmplist2.push_back(id);
                                break;
                        }
                        // Stop when it finds the first MCTrack not produced in emc
                        TString node =  gGeoManager->FindNode(pt->GetStartVertex().X(),pt->GetStartVertex().Y(),pt->GetStartVertex().Z())->GetName();
                        if ( !(node.BeginsWith("emc") || node.BeginsWith("CrystalVol") || node.BeginsWith("Fsc") ) ) {
                                tmplist2.push_back(id);
                                break;
                        }
                        // not exactly clear this part of the code, it needs to be checked (SS, 18/03/2015)
                        id = pt->GetMotherID();
                        for(Int_t k=j-1; k>=0; k--){
                                if(tmplist[k]==id){
                                        tmplist.erase(tmplist.begin()+j);
                                        flag=true;
                                        break;
                                }
                        }
                }
        }
        newlist.clear();
        std::unique_copy( tmplist2.begin(), tmplist2.end(), std::back_inserter( newlist) );
}

void PndEmcHitProducer::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 218 of file PndEmcHitProducer.cxx.

References AcceptDayOne(), AddHit(), c, cleansortmclist(), Double_t, fEmcStr, fEnergyThreshold, fHitArray, fMapper, fMCTrackArray, fNonuniformityPar, fPointArray, fPointMatch, PndEmcXtal::frontCentre(), fTrackEnergy, fTrackEntering, fTrackExiting, fTrackMcTruth, fTrackTime, fUse_nonuniformity, fVerbose, PndEmcPoint::GetEntering(), PndEmcPoint::GetExiting(), PndEmcPoint::GetModule(), PndEmcDigiNonuniformityPar::GetNonuniformityParameters(), PndEmcMapper::GetTCI(), PndEmcStructure::GetTciXtalMap(), HowManyHitsAboveThreshold, HowManyHitsAll, HowManyPoints, idx, max(), PndEmcXtal::normalToFrontFace(), p, point, printf(), push_back(), and zpos.

{
        if (fVerbose>1) cout << " -I- PndEmcHitProducer POINT EXECUTION *********************" << endl;
        // Reset output array
        if (! fHitArray ) Fatal("Exec", "No DigiArray");

        fHitArray->Delete();

        // Declare some variables
        //PndEmcPoint* point  = NULL;
        Int_t DetId;

        fTrackEnergy.clear();
        fTrackTime.clear();
        fTrackMcTruth.clear();
        fPointMatch.clear();
        fTrackEntering.clear();
        fTrackExiting.clear();

        map<Int_t, Float_t>::const_iterator p;

        std::vector<PndEmcPoint*> fPointList;// to pass to EmcHit
        const PndEmcTciXtalMap &XtalMap = fEmcStr->GetTciXtalMap();
        TVector3 frontvec;
        TVector3 normvec;
        TVector3 pointvec;
        TVector3 distvec;
        Double_t zpos;
        Double_t energyscalefactor=1.0;
        Double_t c[3];
        PndEmcXtal *tmpXtal;
        PndEmcTwoCoordIndex *tmpTCI;
        // Loop over EmcPoints
        Int_t nPoints = fPointArray->GetEntriesFast();

        Double_t point_time = 0.00;
        //------- init containers ---

        for (Int_t iPoint = 0; iPoint < nPoints; iPoint++){
                PndEmcPoint* point  = (PndEmcPoint*) fPointArray->At(iPoint);
    if ( ! AcceptDayOne(point) ) continue;
                fTrackEnergy[point->GetDetectorID()] = 0.00;
                fTrackTime  [point->GetDetectorID()] = std::numeric_limits<float>::max();
        }

        //----------------------------

        //Int_t Counter[3] = {0};
        HowManyPoints += nPoints;

        for (Int_t iPoint=0; iPoint<nPoints; iPoint++)
        {
                PndEmcPoint* point  = (PndEmcPoint*) fPointArray->At(iPoint);
    if ( ! AcceptDayOne(point) ) continue;
                DetId = point->GetDetectorID();

                if (point->GetEntering()){
                        fTrackEntering[DetId].AddLinks(point->GetLinksWithType(FairRootManager::Instance()->GetBranchId("MCTrack")));
                }
                if (point->GetExiting()){
                        fTrackExiting[DetId].AddLinks(point->GetLinksWithType(FairRootManager::Instance()->GetBranchId("MCTrack")));
                }
                if(point->GetEnergyLoss() == 0 ) continue;
                if(point->GetModule() == 10 )
                  {
                    cout << " -I-  PndEmcHitProducer::Exec" << "\t" << "Skipping Module 10 (FscFiber)" << endl;
                    continue;
                  }
                //if(point->GetTrackID() < 0)
                //      std::cout<<"negative track id #"<<point->GetTrackID()<<std::endl;

                //std::cout<<"point belongs to track #"<<point->GetTrackID()<<std::endl;
                //if(point->GetTrackID() == 0) ++Counter[0];
                //if(point->GetTrackID() == 1) ++Counter[1];
                //if(point->GetTrackID() == 2) ++Counter[2];

                if(fUse_nonuniformity !=0 ){
                        //light output is z-dependent, so calculate z

                        tmpTCI = fMapper->GetTCI(DetId);
                        if(tmpTCI == NULL){
                                printf("no TCI found for DetectorID %d\n",DetId);
                                continue;
                        }
                        tmpXtal =XtalMap.find(tmpTCI)->second;
                        point->Position(pointvec);
                        frontvec = tmpXtal->frontCentre();
                        normvec = tmpXtal->normalToFrontFace();
                        distvec = pointvec-frontvec;
                        zpos = distvec.Dot(normvec);
                        fNonuniformityPar->GetNonuniformityParameters(DetId,c);
                        energyscalefactor=c[0]+zpos*(c[1]+zpos*c[2]);
                        fTrackEnergy[DetId] += point->GetEnergyLoss() * energyscalefactor;
                        fPointMatch[DetId].push_back(iPoint);
                        //        printf("point with detID %d has z Position %f and energyloss %f scaled with %f\n",DetId,zpos, point->GetEnergyLoss(),energyscalefactor);
                        //        printf("front is at x: %f y: %f z: %f\n", frontvec.X(),frontvec.Y(),frontvec.Z());
                } else {
                        fTrackEnergy[DetId] += point->GetEnergyLoss();
                        fPointMatch[DetId].push_back(iPoint);
                        //        printf("point with detID %d has z Position %f and energyloss %f not scaled\n",DetId,zpos, point->GetEnergyLoss());
                }
                point_time=point->GetTime();

                if (point_time < fTrackTime[point->GetDetectorID()]){
                        fTrackTime[point->GetDetectorID()] = point_time;
                }

                // Check and save MC truth information
                // Eloss==0 tracks are only stored in point, if track is entering detector from outside
                // and thats what we are interested in...
                //std::cout<<"track id #"<<point->GetTrackID()<<", Energyloss #"<<point->GetEnergyLoss()<<endl;
                if(point->GetEnergyLoss() >0)
                        (fTrackMcTruth[point->GetDetectorID()]).push_back(point->GetTrackID());
                //if(point->GetEnergyLoss() == 0 ){
                //      cout << "ELoss== 0 : " <<point->GetEnergyLoss()<<", ID "<<
                //              point->GetTrackID()<<","<<point->GetDetectorID()<<","<<point->GetXPad()<<","<<point->GetYPad()<<endl;
                //}else{
                //      cout << "ELoss>0 : " <<point->GetEnergyLoss()<<", ID "<<
                //              point->GetTrackID()<<","<<point->GetDetectorID()<<","<<point->GetXPad()<<","<<point->GetYPad()<<endl;
                //}
        }

        // Loop over EmcPoint

        // Loop to register EmcHit
        Int_t idx = 0;
        for( p = fTrackEnergy.begin(); p != fTrackEnergy.end(); ++p){
                ++HowManyHitsAll;
                ++idx;
                if ((*p).second > fEnergyThreshold){
                        ++ HowManyHitsAboveThreshold;
                        // Check and save MC truth information B.S.
                        // remove MC Truth particles which are not needed (eg grand^x-daugherts)
                        if( fMCTrackArray){
                                std::vector<Int_t>& plist = fTrackMcTruth[(*p).first];
                                cleansortmclist(fTrackMcTruth[(*p).first],fMCTrackArray);
                                //std::cout<<"The "<<(idx)<<" hit produced by track #";
                                //for(Int_t ip=0;ip<plist.size();++ip){
                                //      cout<<plist[ip]<<", ";
                                //}
                                //cout<<endl;
                        }
                        //std::vector<Int_t>& track = fTrackMcTruth[(*p).first];
                        //std::cout<<track.size()<<" tracks contributes energy to this hit, ";
                        //for(size_t i=0;i<track.size();++i){
                        //      if(i < track.size() - 1){
                        //              std::cout<<track[i]<<", ";
                        //      }else{
                        //              std::cout<<track[i];
                        //      }
                        //}
                        //std::cout<<std::endl;
                        AddHit(1, (*p).first, (*p).second, fTrackTime[(*p).first], fTrackMcTruth[(*p).first], fTrackEntering[(*p).first], fTrackExiting[(*p).first]);
                        //myHit->AddLinks(FairMultiLinkedData("EmcPoint", fPointMatch[p->first]));
                        //myHit->AddLinks(FairMultiLinkedData("MCTrack", fTrackMcTruth[(*p).first));
                }
        }
        //check
        //Int_t nTrack = fMCTrackArray->GetEntriesFast();
        //for(Int_t itrack = 0; itrack < nTrack; ++itrack){
        //      PndMCTrack* pt1 =((PndMCTrack*)fMCTrackArray->At(itrack));
        //      if(pt1->IsGeneratorCreated()){
        //              pt1->Print(itrack);
        //      }
        //}


}

void PndEmcHitProducer::FinishTask

(

)

Definition at line 412 of file PndEmcHitProducer.cxx.

References fDayOne, fEnergyThreshold, HowManyHitsAboveThreshold, HowManyHitsAll, and HowManyPoints.

{
        std::cout<<"========================================================="<<std::endl;
        std::cout<<"PndEmcHitProducer::FinishTask"<<std::endl;
        std::cout<<"*********************************************************"<<std::endl;
  if(fDayOne)   std::cout<<" DAY 1 Setup active, only 12/16 Slices available "<<std::endl;
        std::cout<<"Read points # "<<HowManyPoints<<std::endl;
        std::cout<<"Produc hits# "<<HowManyHitsAll<<", threshold# "<<fEnergyThreshold<<std::endl;
        std::cout<<"Hits above threshhod#"<<HowManyHitsAboveThreshold<<std::endl;
        std::cout<<"*********************************************************"<<std::endl;
}

Bool_t PndPersistencyTask::GetPersistency ( )

inlineinherited

Definition at line 32 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency.

Referenced by PndLmdPixelHitProducerFast::GetPersistance(), PndMdtDigitization::Init(), PndMdtHitProducerIdeal::Init(), PndMdtClusterTask::Init(), PndFtsHitProducerRealFast::Init(), PndDiscTaskReconstruction::Init(), PndRichHitProducer::Init(), PndSttHitProducerRealFast::Init(), PndSttHelixHitProducer::Init(), PndDiscTaskPID::Init(), PndIdealTrackFinder::Init(), PndSttMvdGemTracking::Init(), PndMdtTrkProducer::Init(), PndFtsHitProducerRealFull::Init(), PndLmdPixelClusterTask::Init(), PndSttHitProducerRealFull::Init(), PndLmdStripClusterTask::Init(), PndEmcApdHitProducer::Init(), PndMissingPzCleanerTask::Init(), PndEmcMakeRecoHit::Init(), PndEmcMakeClusterOnline::Init(), PndTrackSmearTask::Init(), PndEmcFWEndcapTimebasedWaveforms::Init(), PndSttHitProducerIdeal::Init(), PndEmcFWEndcapDigi::Init(), PndFtsHitProducerIdeal::Init(), PndEmcMakeCluster::Init(), PndMdtPointsToWaveform::Init(), PndDiscTaskDigitization::Init(), PndEmcMakeDigi::Init(), PndSdsTimeWalkCorrTask::Init(), PndLmdPixelHitProducerFast::Init(), PndDrcHitFinder::Init(), PndRichHitFinder::Init(), PndEmcMakeCorr::Init(), PndFtofHitProducerIdeal::Init(), PndEmcHitsToWaveform::Init(), PndSciTDigiTask::Init(), PndDrcHitProducerIdeal::Init(), PndSdsHitProducerIdeal::Init(), PndSciTHitProducerIdeal::Init(), Init(), PndRecoMultiKalmanTask2::Init(), PndDrcHitProducerReal::Init(), PndDskFLGHitProducerIdeal::Init(), PndEmcTmpWaveformToDigi::Init(), PndDrcDigiTask::Init(), PndEmcWaveformToDigi::Init(), PndSttMatchTracks::Init(), PndEmcWaveformToCalibratedDigi::Init(), PndTrkTracking2::Init(), PndSttFindTracks::Init(), PndEmcMultiWaveformToCalibratedDigi::Init(), PndDrcTimeDigiTask::Init(), PndRecoKalmanTask2::Init(), PndEmcExpClusterSplitter::Init(), PndSdsNoiseProducer::Init(), PndFtsHoughTrackerTask::Init(), PndEmcPhiBumpSplitter::Init(), PndSdsHybridHitProducer::Init(), PndSdsIdealRecoTask::Init(), PndRecoMultiKalmanTask::Init(), PndSdsIdealClusterTask::Init(), PndRecoKalmanTask::Init(), PndSdsStripHitProducerDif::Init(), PndSdsStripHitProducer::Init(), PndGemDigitize::Init(), PndGemFindHits::Init(), PndSdsPixelClusterTask::Init(), PndSdsStripClusterTask::Init(), PndMvdGemTrackFinderOnHits::Init(), PndBarrelTrackFinder::Init(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcMakeBump::PndEmcMakeBump(), PndUnassignedHitsTask::RegisterBranches(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndEmcMakeBump::SetStorageOfData(), and PndEmcFullDigiTask::StoreDigi().

{ return fPersistency; }

InitStatus PndEmcHitProducer::Init ( )

virtual

Virtual method Init

Definition at line 69 of file PndEmcHitProducer.cxx.

References emcX, emcY, emcZ, fDigiPar, fEmcStr, fEnergyThreshold, fGeoPar, fHitArray, fMapper, fMCTrackArray, fNonuniformityFile, fNonuniformityPar, fPointArray, fUse_nonuniformity, PndEmcStructure::GetEmcX(), PndEmcStructure::GetEmcY(), PndEmcStructure::GetEmcZ(), PndEmcDigiPar::GetEnergyHitThreshold(), PndPersistencyTask::GetPersistency(), PndEmcDigiPar::GetUse_nonuniformity(), PndEmcGeoPar::InitEmcMapper(), PndEmcMapper::Instance(), PndEmcStructure::Instance(), printf(), and PndEmcDigiNonuniformityPar::SetNonuniParObject().

                                  {

        cout << " -I- PndEmcHitProducer INITIALIZATION *********************" << endl;

        //FairDetector::Initialize();
        //FairRun* sim = FairRun::Instance();
        //FairRuntimeDb* rtdb=sim->GetRuntimeDb();

        // Get RootManager
        FairRootManager* ioman = FairRootManager::Instance();
        if (! ioman ){
                cout << "-E- PndEmcHitProducer::Init: "
                        << "RootManager not instantiated!" << endl;
                return kFATAL;
        }

        // Get input array
        fPointArray = (TClonesArray*) ioman->GetObject("EmcPoint");
        if (! fPointArray ){
                cout << "-W- PndEmcHitProducer::Init: "
                        << "No EmcPoint array!" << endl;
                return kERROR;
        }

        // Get input array
        fMCTrackArray = (TClonesArray*) ioman->GetObject("MCTrack");
        if (! fMCTrackArray ){
                cout << "-W- PndEmcMakeCluster::Init: "
                        << "No MCTrack array! Needed for MC Truth" << endl;
                //return kERROR;
        }

        // Create and register output array
        fHitArray = new TClonesArray("PndEmcHit");

        ioman->Register("EmcHit","Emc",fHitArray, GetPersistency());

        fGeoPar->InitEmcMapper();
        fMapper=PndEmcMapper::Instance();
        fEmcStr=PndEmcStructure::Instance();

        emcX=fEmcStr->GetEmcX();
        emcY=fEmcStr->GetEmcY();
        emcZ=fEmcStr->GetEmcZ();;

        fEnergyThreshold =fDigiPar->GetEnergyHitThreshold();
        fUse_nonuniformity = fDigiPar->GetUse_nonuniformity();

        if(fUse_nonuniformity){
                cout << "-I- PndEmcHitProducer: Using nonuniform lightoutput" << endl;
        }
        if(fUse_nonuniformity && fNonuniformityFile.Length()>0){
                TFile *nonuniformityfile = new TFile(fNonuniformityFile);
                if(nonuniformityfile==NULL){
                        cout << "-E- PndEmcHitProducer: Could not open file " << fNonuniformityFile.Data() << " for Nonuniformity Information" << endl;
                } else {
                        PndEmcDigiNonuniParObject *parObject;
                        nonuniformityfile->GetObject("PndEmcDigiNonuniParObject",parObject);
                        if(parObject == NULL){
                                cout << "-E- PndEmcHitProducer: Could not get Nonuniformity information from file " << fNonuniformityFile.Data() << endl;
                        } else {
                                fNonuniformityPar->SetNonuniParObject(parObject);
                        }
                }
        }


        printf("HitProducer has EnergyHitThreshold of %f GeV and Use_nonuniformity %i\n", fEnergyThreshold, fUse_nonuniformity);
        cout << "-I- PndEmcHitProducer: Intialization successfull" << endl;

        return kSUCCESS;
}

PndEmcHitProducer& PndEmcHitProducer::operator= ( const PndEmcHitProducer &  )

inlineprivate

Definition at line 106 of file PndEmcHitProducer.h.

{return *this;}

void PndEmcHitProducer::SetDayOne ( bool  d = true )

inline

Definition at line 67 of file PndEmcHitProducer.h.

References d, and fDayOne.

Referenced by PndMasterSimTask::PndMasterSimTask().

{fDayOne=d;};

void PndEmcHitProducer::SetNonuniformityFile ( const char *  filename )

inline

Definition at line 66 of file PndEmcHitProducer.h.

References filename, and fNonuniformityFile.

{fNonuniformityFile = filename;};

void PndEmcHitProducer::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 142 of file PndEmcHitProducer.cxx.

References fDigiPar, fGeoPar, fNonuniformityPar, and run.

                                        {

        // Get run and runtime database
        FairRun* run = FairRun::Instance();
        if (! run ) Fatal("SetParContainers", "No analysis run");

        FairRuntimeDb* db = run->GetRuntimeDb();
        if (! db ) Fatal("SetParContainers", "No runtime database");

        // Get Emc geometry parameter container
        fGeoPar = (PndEmcGeoPar*) db->getContainer("PndEmcGeoPar");

        // Get Emc digitisation parameter container
        fDigiPar = (PndEmcDigiPar*) db->getContainer("PndEmcDigiPar");

        fNonuniformityPar = (PndEmcDigiNonuniformityPar*) db->getContainer("PndEmcDigiNonuniformityPar");

        fDigiPar->setChanged();
        fDigiPar->setInputVersion(run->GetRunId(),1);

        fNonuniformityPar->setChanged();
        fNonuniformityPar->setInputVersion(run->GetRunId(),1);
}

void PndPersistencyTask::SetPersistency ( Bool_t  val = kTRUE )

inlineinherited

Definition at line 31 of file PndPersistencyTask.h.

References PndPersistencyTask::fPersistency, and val.

Referenced by barrelTrackFinder(), digi_complete(), digi_complete_newSTT(), digiOnly_complete(), PndBarrelTrackFinder::PndBarrelTrackFinder(), PndCATracking::PndCATracking(), PndDrcHitFinder::PndDrcHitFinder(), PndEmc2DLocMaxFinder::PndEmc2DLocMaxFinder(), PndEmcExpClusterSplitter::PndEmcExpClusterSplitter(), PndEmcFullDigiTask::PndEmcFullDigiTask(), PndEmcFWEndcapDigi::PndEmcFWEndcapDigi(), PndEmcFWEndcapTimebasedWaveforms::PndEmcFWEndcapTimebasedWaveforms(), PndEmcHitProducer(), PndEmcHitsToWaveform::PndEmcHitsToWaveform(), PndEmcMakeBump::PndEmcMakeBump(), PndEmcMakeCluster::PndEmcMakeCluster(), PndEmcMakeClusterOnline::PndEmcMakeClusterOnline(), PndEmcMakeDigi::PndEmcMakeDigi(), PndEmcMakeRecoHit::PndEmcMakeRecoHit(), PndEmcMultiWaveformToCalibratedDigi::PndEmcMultiWaveformToCalibratedDigi(), PndEmcPhiBumpSplitter::PndEmcPhiBumpSplitter(), PndEmcTmpWaveformToDigi::PndEmcTmpWaveformToDigi(), PndEmcWaveformToCalibratedDigi::PndEmcWaveformToCalibratedDigi(), PndEmcWaveformToDigi::PndEmcWaveformToDigi(), PndFtofHitProducerIdeal::PndFtofHitProducerIdeal(), PndFtsCATracking::PndFtsCATracking(), PndFtsHitProducerIdeal::PndFtsHitProducerIdeal(), PndFtsHitProducerRealFast::PndFtsHitProducerRealFast(), PndFtsHitProducerRealFull::PndFtsHitProducerRealFull(), PndFtsHoughTrackerTask::PndFtsHoughTrackerTask(), PndGemDigitize::PndGemDigitize(), PndGemFindHits::PndGemFindHits(), PndIdealTrackFinder::PndIdealTrackFinder(), PndLmdPixelClusterTask::PndLmdPixelClusterTask(), PndLmdPixelHitProducerFast::PndLmdPixelHitProducerFast(), PndMdtClusterTask::PndMdtClusterTask(), PndMdtDigitization::PndMdtDigitization(), PndMdtHitProducerIdeal::PndMdtHitProducerIdeal(), PndMdtPointsToWaveform::PndMdtPointsToWaveform(), PndMdtTrkProducer::PndMdtTrkProducer(), PndMissingPzCleanerTask::PndMissingPzCleanerTask(), PndMvdGemTrackFinderOnHits::PndMvdGemTrackFinderOnHits(), PndMvdHitProducerIdeal::PndMvdHitProducerIdeal(), PndMvdPixelClusterTask::PndMvdPixelClusterTask(), PndMvdTimeWalkCorrTask::PndMvdTimeWalkCorrTask(), PndMvdToPix4ClusterTask::PndMvdToPix4ClusterTask(), PndRecoKalmanTask::PndRecoKalmanTask(), PndRecoKalmanTask2::PndRecoKalmanTask2(), PndRecoMultiKalmanTask::PndRecoMultiKalmanTask(), PndRecoMultiKalmanTask2::PndRecoMultiKalmanTask2(), PndRichHitFinder::PndRichHitFinder(), PndRichHitProducer::PndRichHitProducer(), PndSciTDigiTask::PndSciTDigiTask(), PndSciTHitProducerIdeal::PndSciTHitProducerIdeal(), PndSdsHitProducerIdeal::PndSdsHitProducerIdeal(), PndSdsHybridHitProducer::PndSdsHybridHitProducer(), PndSdsIdealClusterTask::PndSdsIdealClusterTask(), PndSdsIdealRecoTask::PndSdsIdealRecoTask(), PndSdsNoiseProducer::PndSdsNoiseProducer(), PndSdsPixelClusterTask::PndSdsPixelClusterTask(), PndSdsStripClusterTask::PndSdsStripClusterTask(), PndSdsStripHitProducer::PndSdsStripHitProducer(), PndSdsTimeWalkCorrTask::PndSdsTimeWalkCorrTask(), PndSttFindTracks::PndSttFindTracks(), PndSttHelixHitProducer::PndSttHelixHitProducer(), PndSttHitProducerIdeal::PndSttHitProducerIdeal(), PndSttHitProducerRealFast::PndSttHitProducerRealFast(), PndSttHitProducerRealFull::PndSttHitProducerRealFull(), PndSttMatchTracks::PndSttMatchTracks(), PndSttMvdGemTracking::PndSttMvdGemTracking(), PndTrackSmearTask::PndTrackSmearTask(), PndTrkTracking2::PndTrkTracking2(), reco(), reco_complete(), reco_complete_gf2(), reco_complete_newSTT(), reco_complete_sec(), recoideal_complete(), PndMvdClusterTask::SetPersistance(), PndMvdDigiTask::SetPersistance(), PndLmdPixelHitProducerFast::SetPersistance(), PndSdsHitProducerIdeal::SetPersistance(), PndSttMvdGemTracking::SetPersistenc(), PndMdtClusterTask::SetPersistence(), PndSttHelixHitProducer::SetPersistence(), PndMissingPzCleanerTask::SetPersistence(), PndFtsHitProducerRealFast::SetPersistence(), PndFtsHitProducerRealFull::SetPersistence(), PndSttHitProducerIdeal::SetPersistence(), PndSttHitProducerRealFull::SetPersistence(), PndSttHitProducerRealFast::SetPersistence(), PndFtsHitProducerIdeal::SetPersistence(), PndTrackSmearTask::SetPersistence(), PndSciTHitProducerIdeal::SetPersistence(), PndIdealTrackFinder::SetPersistence(), PndSttMatchTracks::SetPersistence(), PndSttFindTracks::SetPersistence(), PndFtsHoughTrackerTask::SetPersistence(), PndTrkTracking2::SetPersistence(), PndEmcMakeRecoHit::SetStorageOfData(), PndEmcMakeClusterOnline::SetStorageOfData(), PndEmcFWEndcapDigi::SetStorageOfData(), PndEmcFWEndcapTimebasedWaveforms::SetStorageOfData(), PndEmcMakeDigi::SetStorageOfData(), PndMdtPointsToWaveform::SetStorageOfData(), PndEmc2DLocMaxFinder::SetStorageOfData(), PndEmcMakeCluster::SetStorageOfData(), PndEmcHitsToWaveform::SetStorageOfData(), PndEmcMakeBump::SetStorageOfData(), PndEmcTmpWaveformToDigi::SetStorageOfData(), PndEmcWaveformToDigi::SetStorageOfData(), PndEmcWaveformToCalibratedDigi::SetStorageOfData(), PndEmcMultiWaveformToCalibratedDigi::SetStorageOfData(), PndEmcExpClusterSplitter::SetStorageOfData(), PndEmcPhiBumpSplitter::SetStorageOfData(), standard_tracking(), and PndEmcFullDigiTask::StoreDigi().

{ fPersistency = val; }

void PndEmcHitProducer::SetStorageOfData

(

Bool_t

val

)

Definition at line 406 of file PndEmcHitProducer.cxx.

References fStoreHits, and val.

Referenced by digi_sttcombi(), emc(), emc_complete(), emc_complete_corr(), emc_correction_data_production(), emc_correction_QA_data_production(), and error_matrix_data_production().

{
        fStoreHits=val;
        return;
}

Member Data Documentation

mapper PndEmcHitProducer::emcX

private

Definition at line 87 of file PndEmcHitProducer.h.

Referenced by AddHit(), and Init().

mapper PndEmcHitProducer::emcY

private

Definition at line 88 of file PndEmcHitProducer.h.

Referenced by AddHit(), and Init().

mapper PndEmcHitProducer::emcZ

private

Definition at line 89 of file PndEmcHitProducer.h.

Referenced by AddHit(), and Init().

bool PndEmcHitProducer::fDayOne

private

Definition at line 117 of file PndEmcHitProducer.h.

Referenced by AcceptDayOne(), FinishTask(), and SetDayOne().

PndEmcDigiPar* PndEmcHitProducer::fDigiPar

private

Definition at line 94 of file PndEmcHitProducer.h.

Referenced by Init(), and SetParContainers().

PndEmcStructure* PndEmcHitProducer::fEmcStr

private

Definition at line 91 of file PndEmcHitProducer.h.

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

Float_t PndEmcHitProducer::fEnergyThreshold

private

Definition at line 85 of file PndEmcHitProducer.h.

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

PndEmcGeoPar* PndEmcHitProducer::fGeoPar

private

Digitisation parameter container

Definition at line 95 of file PndEmcHitProducer.h.

Referenced by Init(), and SetParContainers().

TClonesArray* PndEmcHitProducer::fHitArray

private

Output array of PndEmcHit

Definition at line 78 of file PndEmcHitProducer.h.

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

PndEmcMapper* PndEmcHitProducer::fMapper

private

Definition at line 92 of file PndEmcHitProducer.h.

Referenced by Exec(), and Init().

Int_t PndEmcHitProducer::fMapVersion

private

Geo file to use

Definition at line 84 of file PndEmcHitProducer.h.

TClonesArray* PndEmcHitProducer::fMCTrackArray

private

Definition at line 75 of file PndEmcHitProducer.h.

Referenced by Exec(), and Init().

TString PndEmcHitProducer::fNonuniformityFile

private

Definition at line 72 of file PndEmcHitProducer.h.

Referenced by Init(), PndEmcHitProducer(), and SetNonuniformityFile().

PndEmcDigiNonuniformityPar* PndEmcHitProducer::fNonuniformityPar

private

Geometry parameter container

Definition at line 96 of file PndEmcHitProducer.h.

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

TClonesArray* PndEmcHitProducer::fPointArray

private

Input array of PndEmcPoints

Definition at line 74 of file PndEmcHitProducer.h.

Referenced by Exec(), and Init().

map<Int_t, std::vector <Int_t> > PndEmcHitProducer::fPointMatch

private

Definition at line 113 of file PndEmcHitProducer.h.

Referenced by Exec().

Bool_t PndEmcHitProducer::fStoreHits

private

Definition at line 101 of file PndEmcHitProducer.h.

Referenced by SetStorageOfData().

map<Int_t, Float_t> PndEmcHitProducer::fTrackEnergy

private

Definition at line 110 of file PndEmcHitProducer.h.

Referenced by Exec().

map<Int_t, FairMultiLinkedData> PndEmcHitProducer::fTrackEntering

private

Definition at line 114 of file PndEmcHitProducer.h.

Referenced by Exec().

map<Int_t, FairMultiLinkedData> PndEmcHitProducer::fTrackExiting

private

Definition at line 115 of file PndEmcHitProducer.h.

Referenced by Exec().

map<Int_t, std::vector <Int_t> > PndEmcHitProducer::fTrackMcTruth

private

Definition at line 112 of file PndEmcHitProducer.h.

Referenced by Exec().

map<Int_t, Float_t> PndEmcHitProducer::fTrackTime

private

Definition at line 111 of file PndEmcHitProducer.h.

Referenced by Exec().

Int_t PndEmcHitProducer::fUse_nonuniformity

private

Definition at line 71 of file PndEmcHitProducer.h.

Referenced by Exec(), and Init().

TObjArray* PndEmcHitProducer::fVolumeArray

private

Definition at line 81 of file PndEmcHitProducer.h.

---

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

• PndEmcHitProducer.h
• PndEmcHitProducer.cxx