PndEmcAnalysis Class Reference

#include <PndEmcAnalysis.h>

Inheritance diagram for PndEmcAnalysis:

Public Member Functions
	PndEmcAnalysis (Int_t verbose=0, Bool_t storeRootFile=kTRUE)
virtual	~PndEmcAnalysis ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
void	SetOutput (TString name)
void	RunTimeBased ()
virtual void	FinishTask ()

Public Attributes
Bool_t	fSaveWave
Bool_t	fSaveDigi
Bool_t	fSaveBump
Bool_t	fSaveTask
Bool_t	fSaveHits
Bool_t	fSaveMcTruth
Bool_t	fSaveReco
Bool_t	fSavePidCharged
Bool_t	fSavePidNeutral
Bool_t	fSaveJpsiTo3g

Private Member Functions
void	SetParContainers ()
	PndEmcAnalysis (const PndEmcAnalysis &L)
PndEmcAnalysis &	operator= (const PndEmcAnalysis &)
	ClassDef (PndEmcAnalysis, 1)

Private Attributes
TClonesArray *	fEmcHitsArray
TClonesArray *	fWaveformArray
TClonesArray *	fDigiArray
TClonesArray *	fSharedDigiArray
TClonesArray *	fClusterArray
TClonesArray *	fBumpArray
TClonesArray *	fMcTrackArray
TClonesArray *	fEvtHeaderArray
TClonesArray *	fChargedCand
TClonesArray *	fNeutralCand
Int_t	fVerbose
Bool_t	fStoreRooTFile
BinaryFunctor *	fFunctor
Bool_t	fTimeOrderedDigi
	set to kTRUE to use the time ordering of the output data. More...
std::map< Int_t, std::pair < Int_t, Int_t > >	evtMap
std::set< Int_t >	pevtset
std::set< Int_t >	evtset
TString	fileName
TFile *	fRootFile
TTree *	tHit
TTree *	tDigi
TTree *	tWave
TTree *	tBump
TTree *	tTask
TTree *	tEvtPileup
TTree *	tTaskPileup
TTree *	tMcTruth
TTree *	tRec
Double_t	fHitTime
Double_t	fHitEnergy
Double_t	fTheta
Double_t	fPhi
Int_t	fNPoint
Double_t	fEvtTime
Int_t	fModule
Int_t	fHitIndex
Double_t	fWaveTime
Double_t	fWavePeak
Int_t	fWaveLen
Int_t	fWaveIdx
Int_t	fEventNo
Double_t	fBaseline
Int_t	fXcor
Int_t	fYcor
Int_t	fPileupCount
Double_t	fDigiTime
Double_t	fDigiEnergy
Double_t	fDistance
TVector3	fPosition
TVector3	fSeedPosition
Double_t	fEnergy1C
Double_t	fEnergy2C
Double_t	fBumpTime1
Double_t	fBumpTime2
Double_t	fBumpTime3
Double_t	fSeedTime
Double_t	fBumpEnergy
Int_t	fNumDigi
Int_t	fNumBump
Int_t	fClusterIdx
Int_t	fNumOfSeed
Int_t	fDigiMapSize
Int_t	fNumOfMatch
Int_t	fNumOfDismatch
Double_t	totDigiEnergy
Double_t	totSharedDigiEnergy
Int_t	fNumOfHits
Int_t	fNumWaveform
Int_t	fNumSharedDigi
Int_t	fNumCluster
Int_t	fNumMcTrack
Int_t	fTotWave
Int_t	fPileupWave
Int_t	totNumOfWave
Int_t	totNumOfSubWave
Int_t	totNumOfDigi
Int_t	totNumOfBump
Int_t	totNumOfSharedDigi
Int_t	totNumOfEvents
Int_t	totNumOfPileupEvents
Double_t	fp4 [4 *NElement]
Double_t	fMcp4 [4 *NElement]
Double_t	fPositionV [3 *NElement]
Double_t	fEnergy [NElement]
Double_t	fEnergyC1 [NElement]
Double_t	fEnergyC2 [NElement]
Double_t	fZ201 [NElement]
Double_t	fZ531 [NElement]
Double_t	fLat1 [NElement]
Int_t	fModule1 [NElement]
Int_t	fPDGCode [NElement]
Double_t	fTheta1 [NElement]
Double_t	fPhi1 [NElement]
Double_t	fRestEnergy1
Double_t	fRestEnergy2
Double_t	fRestEnergy3
Int_t	fNumGT
Int_t	fDet1 [NElement]
Int_t	fHitEmc1 [NElement]
Int_t	fHitEmc2 [NElement]
Int_t	fNPointV [NElement]
Int_t	fchar [NElement]
Double_t	fpos [NElement]
Double_t	fmom [NElement]
Int_t	fmcidx [NElement]
Double_t	frawe [NElement]
Double_t	fcale [NElement]
Int_t	fndigi [NElement]
Int_t	fmod [NElement]
Int_t	fcluIdx [NElement]
Double_t	fe25 [NElement]
Double_t	fe9 [NElement]
Int_t	fnbump [NElement]
TTree *	tCharged
TTree *	tNeutral
TTree *	tAna3g
Int_t	fNumberOfGoodPhoton
Double_t	fExternalEnergy
Double_t	f4p1 [4]
Double_t	fMc4p1 [4]
Double_t	fpos1 [3]
Double_t	fmom1
Double_t	fenergy1
Double_t	femcraw1
Int_t	fndigi1
Int_t	fmod1
Int_t	fMcTrack1
Double_t	fe251
Double_t	flat1
Double_t	fz201
Double_t	fz531
Double_t	f4p2 [4]
Double_t	fMc4p2 [4]
Double_t	fpos2 [3]
Double_t	fmom2
Double_t	fenergy2
Double_t	femcraw2
Int_t	fndigi2
Int_t	fmod2
Int_t	fMcTrack2
Double_t	fe252
Double_t	flat2
Double_t	fz202
Double_t	fz532
Double_t	f4p3 [4]
Double_t	fMc4p3 [4]
Double_t	fpos3 [3]
Double_t	fmom3
Double_t	fenergy3
Double_t	femcraw3
Int_t	fndigi3
Int_t	fmod3
Int_t	fMcTrack3
Double_t	fe253
Double_t	flat3
Double_t	fz203
Double_t	fz533
PndEmcDigiCalibrator	digiCalibrator

Static Private Attributes
static const Int_t	NElement = 30

Detailed Description

Definition at line 49 of file PndEmcAnalysis.h.

Constructor & Destructor Documentation

PndEmcAnalysis::PndEmcAnalysis	(	Int_t	verbose = 0,
		Bool_t	storeRootFile = kTRUE
	)

Definition at line 81 of file PndEmcAnalysis.cxx.

References fileName, fSaveBump, fSaveDigi, fSaveHits, fSaveJpsiTo3g, fSaveMcTruth, fSavePidCharged, fSavePidNeutral, fSaveReco, fSaveTask, fSaveWave, and fTimeOrderedDigi.

                                                             :
        fEmcHitsArray(0)
        ,fWaveformArray(0)
        ,fDigiArray(0)
        ,fSharedDigiArray(0)
        ,fClusterArray(0)
        ,fBumpArray(0)
        ,fMcTrackArray(0)
        ,fEvtHeaderArray(0)
        ,fVerbose(verbose)
        ,fStoreRooTFile(storeFile)
{
        fileName="PndEmcAnalysis.root";
        fTimeOrderedDigi = kFALSE;
        fSaveWave = kFALSE;
        fSaveDigi = kFALSE;
        fSaveBump = kFALSE;
        fSaveTask = kFALSE;
        fSaveHits = kFALSE;
        fSaveMcTruth = kFALSE;
        fSaveReco = kFALSE;
        fSavePidCharged = kFALSE;
        fSavePidNeutral = kFALSE;
        fSaveJpsiTo3g = kFALSE;
}

PndEmcAnalysis::~PndEmcAnalysis ( )

virtual

Definition at line 111 of file PndEmcAnalysis.cxx.

{
}

PndEmcAnalysis::PndEmcAnalysis ( const PndEmcAnalysis &  L )

private

Member Function Documentation

PndEmcAnalysis::ClassDef ( PndEmcAnalysis  , 1    )

private

void PndEmcAnalysis::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 428 of file PndEmcAnalysis.cxx.

References Bool_t, CAMath::Cos(), ctime, PndEmcCluster::DigiList(), Double_t, PndEmcAbsClusterCalibrator::Energy(), PndEmcCluster::energy(), evtMap, evtset, f4p1, f4p2, f4p3, fBaseline, fBumpArray, fBumpEnergy, fBumpTime1, fcale, fchar, fChargedCand, fcluIdx, fClusterArray, fClusterIdx, fDet1, fDigiArray, fDigiEnergy, fDigiMapSize, fDigiTime, fDistance, fe25, fe251, fe252, fe253, fe9, fEmcHitsArray, femcraw1, femcraw2, femcraw3, fEnergy, fenergy1, fEnergy1C, fenergy2, fenergy3, fEnergyC1, fEnergyC2, fEventNo, PndEmcDigi::fEvtNo, fEvtTime, fExternalEnergy, fFunctor, fHitEmc1, fHitEmc2, fHitEnergy, fHitIndex, fHitTime, fLat1, flat1, flat2, flat3, fMc4p1, fMc4p2, fMc4p3, fmcidx, fMcp4, fMcTrack1, fMcTrack2, fMcTrack3, fMcTrackArray, fmod, fmod1, fmod2, fmod3, fModule, fModule1, fmom, fnbump, fndigi, fndigi1, fndigi2, fndigi3, fNeutralCand, fNPoint, fNPointV, fNumberOfGoodPhoton, fNumBump, fNumCluster, fNumDigi, fNumGT, fNumMcTrack, fNumOfDismatch, fNumOfHits, fNumOfMatch, fNumOfSeed, fNumSharedDigi, fNumWaveform, fp4, fPDGCode, fPhi, fPhi1, fPileupCount, fpos, fpos1, fpos2, fpos3, fPosition, fPositionV, frawe, fRestEnergy1, fRestEnergy2, fRestEnergy3, fSaveBump, fSaveDigi, fSaveHits, fSaveJpsiTo3g, fSaveMcTruth, fSavePidCharged, fSavePidNeutral, fSaveReco, fSaveTask, fSaveWave, fSeedPosition, fSeedTime, fSharedDigiArray, fTheta, fTheta1, fTimeOrderedDigi, fVerbose, fWaveformArray, fWaveIdx, fWaveLen, fWavePeak, fWaveTime, fXcor, fYcor, fZ201, fz201, fz202, fz203, fZ531, fz531, fz532, fz533, PndMCTrack::Get4Momentum(), PndEmcWaveform::GetBaseline(), PndPidCandidate::GetCharge(), PndEmcBump::GetClusterIndex(), PndPidCandidate::GetEmcCalEnergy(), PndPidCandidate::GetEmcClusterE25(), PndPidCandidate::GetEmcClusterE9(), PndPidCandidate::GetEmcClusterLat(), PndPidCandidate::GetEmcClusterZ20(), PndPidCandidate::GetEmcClusterZ53(), PndPidCandidate::GetEmcIndex(), PndPidCandidate::GetEmcModule(), PndPidCandidate::GetEmcNumberOfBumps(), PndPidCandidate::GetEmcNumberOfCrystals(), PndPidCandidate::GetEmcRawEnergy(), PndEmcHit::GetEnergy(), PndPidCandidate::GetEnergy(), PndEmcDigi::GetEnergy(), PndEmcWaveform::GetEvtList(), PndEmcWaveform::GetHitIndex(), PndEmcDigi::GetHitIndex(), PndPidCandidate::GetLorentzVector(), PndPidCandidate::GetMcIndex(), PndEmcHit::GetModule(), PndEmcCluster::GetModule(), PndEmcWaveform::GetModule(), PndEmcDigi::GetModule(), PndPidCandidate::GetMomentum(), PndMCTrack::GetMomentum(), PndMCTrack::GetMotherID(), PndMCTrack::GetNPoints(), PndMCTrack::GetPdgCode(), PndEmcHit::GetPhi(), PndEmcWaveform::GetPileupCount(), PndEmcHit::GetPointList(), PndPidCandidate::GetPosition(), PndEmcWaveform::GetTCI(), PndEmcHit::GetTheta(), PndEmcHit::GetTime(), PndEmcWaveform::GetWaveformLength(), PndEmcHit::GetXPad(), PndEmcHit::GetYPad(), i, idx, PndMCTrack::IsGeneratorCreated(), kEMC, PndEmcCluster::LatMom(), PndEmcCluster::LocalMaxMap(), PndEmcClusterCalibrator::MakeEmcClusterCalibrator(), PndEmcWaveform::Max(), PndEmcCluster::MemberDigiMap(), PndEmcCluster::NBumps(), NElement, nHits, PndEmcCluster::NumberOfDigis(), p1, p2, pevtset, pos, PndMCTrack::Print(), rtime, CAMath::Sin(), tAna3g, tBump, tCharged, tDigi, tHit, timer, tMcTruth, tNeutral, totDigiEnergy, totNumOfBump, totNumOfDigi, totNumOfSharedDigi, totNumOfSubWave, totNumOfWave, totSharedDigiEnergy, tRec, tTask, tWave, PndEmcAbsClusterCalibrator::Where(), PndEmcCluster::where(), PndEmcDigi::where(), PndEmcTwoCoordIndex::XCoord(), PndEmcTwoCoordIndex::YCoord(), PndEmcCluster::Z20(), and PndEmcCluster::Z53().

{
        TStopwatch timer;
        if (fVerbose>2){
                timer.Start();
        }

 
        fEvtTime = FairRootManager::Instance()->GetEventTime();
        cout<<"=======================PndEmcAnalysis======================="<<endl;
        if(fSaveHits){
                Int_t nHits = fEmcHitsArray->GetEntriesFast();
                for(Int_t iHit =0; iHit < nHits; ++iHit){
                        PndEmcHit* theHit = (PndEmcHit*)fEmcHitsArray->At(iHit);
                        fHitTime   = theHit->GetTime(); 
                        fModule    = theHit->GetModule();
                        fHitEnergy = theHit->GetEnergy();
                        fTheta     = theHit->GetTheta(); 
                        fPhi       = theHit->GetPhi();
                        fNPoint    = theHit->GetPointList().size();
                        fXcor      = theHit->GetXPad();
                        fYcor      = theHit->GetYPad();
                        tHit->Fill();
                }
                if(fVerbose>0)
                        cout<<"Save "<<nHits<<" PndEmcHit done"<<endl;
        }
        if(fTimeOrderedDigi){
                Double_t time_length = 40.;//99.98%
                if(FairRunAna::Instance()->IsTimeStamp()){
                        fWaveformArray  = FairRootManager::Instance()->GetData("EmcSortedWaveform"
                                        , fFunctor
                                        , FairRootManager::Instance()->GetEventTime() + time_length);
                }
        }

        Int_t nWaveforms = fWaveformArray->GetEntriesFast();
        totNumOfWave += nWaveforms;
        if(fSaveWave){
                for (Int_t iWaveform=0; iWaveform<nWaveforms; iWaveform++) 
                {
                        PndEmcWaveform* theWaveform = (PndEmcWaveform*) fWaveformArray->At(iWaveform);
                        if(fTimeOrderedDigi){
                                const std::vector<Int_t>& evtList = theWaveform->GetEvtList();
                                for(size_t i=0; i< evtList.size(); ++i){
                                        evtset.insert(evtList[i]);
                                        ++ evtMap[evtList[i]].first;//all waveforms counter
                                        if(evtList.size() > 1){
                                                pevtset.insert(evtList[i]);
                                                ++ evtMap[evtList[i]].second;//pileup waveforms counter
                                        }
                                }
                        }

                        fPileupCount= theWaveform->GetPileupCount();
                        totNumOfSubWave += (fPileupCount+1);
                        fHitIndex   = theWaveform->GetHitIndex();
                        fModule                 = theWaveform->GetModule();             
                        fWaveTime               = theWaveform->GetTimeStamp();  
                        fWavePeak         = theWaveform->Max();
                        fWaveLen    = theWaveform->GetWaveformLength();
                        fWaveIdx          = iWaveform;
                        fBaseline   = theWaveform->GetBaseline();
                        fXcor           = theWaveform->GetTCI()->XCoord();
                        fYcor       = theWaveform->GetTCI()->YCoord();              
                        if(!fTimeOrderedDigi)
                                fHitEnergy  = ((PndEmcHit*)fEmcHitsArray->At(fHitIndex))->GetEnergy();
                        tWave->Fill();
                }
                if(fVerbose>0)
                        cout<<"Save "<<nWaveforms<<" PndEmcWaveform done"<<endl;
        }

        totDigiEnergy = 0.;
        totSharedDigiEnergy = 0.;
        Int_t nDigi = fDigiArray->GetEntriesFast();
        Int_t nSharedDigi = fSharedDigiArray->GetEntriesFast();
        if(!fTimeOrderedDigi){
                totNumOfDigi += nDigi;
                totNumOfSharedDigi += nSharedDigi;
        }
        if(fSaveDigi){
                for (Int_t iDigi=0; iDigi<nDigi; iDigi++) 
                {
                        PndEmcDigi* theDigi = (PndEmcDigi*) fDigiArray->At(iDigi);
                        totDigiEnergy += theDigi->GetEnergy();
                        fDigiTime                 = theDigi->GetTimeStamp();
                        fModule                                 = theDigi->GetModule();
                        fWaveIdx                                = theDigi->GetHitIndex();
                        if(!fTimeOrderedDigi){
                                fWaveTime                               = ((PndEmcWaveform*)fWaveformArray->At(fWaveIdx))->GetTimeStamp();
                                fHitEnergy      = ((PndEmcHit*)fEmcHitsArray->At(fWaveIdx))->GetEnergy();
                        }
                        fDigiEnergy         = theDigi->GetEnergy();
                        fDistance                               = theDigi->where().Mag();
                        fPosition               = theDigi->where();
                        tDigi->Fill();
                }
                for (Int_t iDigi=0; iDigi<nSharedDigi; iDigi++) 
                {
                        PndEmcDigi* theDigi = (PndEmcDigi*) fSharedDigiArray->At(iDigi);
                        totSharedDigiEnergy += theDigi->GetEnergy();
                }
                if(fVerbose>0)
                        cout<<"Save "<<nDigi<<" PndEmcDigi done"<<endl;
        }
        Int_t nBump = fBumpArray->GetEntriesFast();
        totNumOfBump += nBump;
        std::vector<key> goodHits;
        if(fSaveBump){
                PndEmcAbsClusterCalibrator * calibrator1= PndEmcClusterCalibrator::MakeEmcClusterCalibrator(1);
                PndEmcClusterCalibrator::MakeEmcClusterCalibrator(2); //PndEmcAbsClusterCalibrator * calibrator2= //[R.K.03/2017] unused variable
                for (Int_t iBump=0; iBump<nBump; iBump++) 
                {
                        PndEmcBump* theBump = (PndEmcBump*) fBumpArray->At(iBump);
                        if((fEnergy1C=calibrator1->Energy(theBump)) > 0.050){//50MeV
                                goodHits.push_back(key(iBump, fEnergy1C));
                        }
                        fModule                         =       theBump->GetModule();
                        fBumpEnergy       = theBump->energy();
                        fPosition       =       theBump->where();
                        fNumDigi      = theBump->NumberOfDigis(); 
                        fNumBump                        =       theBump->NBumps();
                        fClusterIdx       = theBump->GetClusterIndex();
                        fNumOfSeed    = theBump->LocalMaxMap().size();
                        fDigiMapSize  = theBump->MemberDigiMap().size();
                        const std::vector<Int_t>& digis = theBump->DigiList();
                        if(fTimeOrderedDigi){
                                totNumOfDigi += digis.size();
                                totNumOfSharedDigi += digis.size();
                        }
                        PndEmcDigi* seedDigi(0);
                        Double_t maxE(-1.);
                        for(size_t id=0; id < digis.size(); ++id){
                                PndEmcDigi* theDigi = (PndEmcDigi*)fSharedDigiArray->At(digis[id]);
                                if(theDigi->GetEnergy() > maxE){
                                        maxE = theDigi->GetEnergy();
                                        seedDigi = theDigi;
                                }
                                theDigi->fEvtNo == fEventNo ? fNumOfMatch++ : fNumOfDismatch ++;
                        }
                        fSeedTime                       = seedDigi == 0? 0 :/* digiCalibrator.CalibrationEvtTimeByDigi(seedDigi, kFALSE);*/seedDigi->GetTimeStamp();
                        fBumpTime1                      = theBump->GetTimeStamp();
                        //fBumpTime2                    = theBump->GetTimeStamp1();
                        fSeedPosition   = seedDigi->where();
                        //fBumpTime3                    = theBump->GetTimeStamp3();
                        tBump->Fill();
                }
                if(fVerbose>0)
                        cout<<"Save "<<nBump<<" PndEmcBump done"<<endl;
        }
        if(fSaveReco){
                std::vector<key> goodTrack;
                //std::vector<IndexKey> indicator;
                sort(goodHits.begin(), goodHits.end(), std::greater<key>());
                if(goodHits.size() >=1){
                        //reset
                        for(Int_t idx=0; idx < NElement; ++idx){
                                fEnergy[idx] = fEnergyC1[idx] = fEnergyC2[idx] = fTheta1[idx] = fPhi1[idx]
                                        = fZ201[idx] = fZ531[idx] = fLat1[idx] = fModule1[idx] = 0;
                                fPositionV[idx*3] = fPositionV[idx*3+1] = fPositionV[idx*3+2] = 0;
                                fp4[idx*4] = fp4[idx*4+1] = fp4[idx*4+2] = fp4[idx*4+3] = 0;
                        }
                        if(fMcTrackArray){
                                Int_t nTrack = fMcTrackArray->GetEntriesFast() > NElement ? NElement : fMcTrackArray->GetEntriesFast();
                                for(Int_t i=0; i < nTrack; ++i){
                                        PndMCTrack* p1 = (PndMCTrack*)fMcTrackArray->At(i);
                                        if(p1->IsGeneratorCreated()){
                                                goodTrack.push_back(key(i, p1->Get4Momentum().E()));
                                                //TLorentzVector p4(p1->GetMomentum(), p1->GetE());
                                                //memcpy(&fMcp4[i*4], &p4[0], 4*sizeof(Double_t)); 
                                        }else
                                                break;
                                }
                                sort(goodTrack.begin(), goodTrack.end(), std::greater<key>());
                                for(size_t i=0; i< goodTrack.size(); ++i){
                                        PndMCTrack* p1 = (PndMCTrack*)fMcTrackArray->At(goodTrack[i].iHit);
                                        TLorentzVector p4(p1->GetMomentum(), p1->Get4Momentum().E());
                                        memcpy(&fMcp4[i*4], &p4[0], 4*sizeof(Double_t)); 
                                        fPDGCode[i] = p1->GetPdgCode();
                                }
                        }

                        PndEmcAbsClusterCalibrator * calibrator1= PndEmcClusterCalibrator::MakeEmcClusterCalibrator(1);
                        PndEmcAbsClusterCalibrator * calibrator2= PndEmcClusterCalibrator::MakeEmcClusterCalibrator(2);
                        Int_t NEL = goodHits.size() > NElement ? NElement : goodHits.size();
                        fRestEnergy1 = fRestEnergy2 = fRestEnergy3 = 0.;
                        for(size_t i=NEL;i<goodHits.size();++i){
                                PndEmcBump* theBump = (PndEmcBump*) fBumpArray->At(goodHits[i].iHit);
                                fRestEnergy1 += theBump->energy();
                                fRestEnergy2 += calibrator1->Energy(theBump);
                                fRestEnergy3 += calibrator1->Energy(theBump);
                        }
                        fNumBump = goodHits.size();
                        Int_t locIdx =0;

                        for(Int_t idx=0; idx < NEL ; ++idx){
                                PndEmcBump* theBump = (PndEmcBump*)fBumpArray->At(goodHits[idx].iHit);
                                //PndEmcBump* theBump = (PndEmcBump*)fBumpArray->At(index);
                                //if(idx < goodTrack.size())
                                //      locIdx = goodTrack[idx].iHit;
                                //else
                                //      locIdx = idx;
                                locIdx = idx;
                                fEnergy[locIdx]   = theBump->energy();
                                fEnergyC1[locIdx] = calibrator1->Energy(theBump)/1.018;
                                fEnergyC2[locIdx] = calibrator2->Energy(theBump);
                                TVector3 pos(calibrator1->Where(theBump));
                                memcpy(&fPositionV[locIdx*3], &pos[0], 3*sizeof(Double_t));
                                fTheta1[locIdx]   = pos.Theta();
                                fPhi1[locIdx]     = pos.Phi();
                                fZ201[locIdx]     = theBump->Z20();
                                fZ531[locIdx]     = theBump->Z53();
                                fLat1[locIdx]     = theBump->LatMom();
                                fModule1[locIdx]  = theBump->GetModule();
                                TLorentzVector p4(fEnergyC1[locIdx]*TMath::Sin(fTheta1[locIdx])*TMath::Cos(fPhi1[locIdx])
                                                , fEnergyC1[locIdx]*TMath::Sin(fTheta1[locIdx])*TMath::Sin(fPhi1[locIdx])
                                                , fEnergyC1[locIdx]*TMath::Cos(fTheta1[locIdx]), fEnergyC1[locIdx]);
                                memcpy(&fp4[locIdx*4], &p4[0], 4*sizeof(Double_t)); 
                        }
                        tRec->Fill();
                        if(fVerbose>0)
                                cout<<"Save "<<NEL<<" reconstruction data done"<<endl;
                }
        }
        if(fSaveMcTruth){
                Bool_t isFill = kFALSE; 
                Int_t nTrack = fMcTrackArray->GetEntriesFast();// > NElement ? NElement : fMcTrackArray->GetEntriesFast();
                std::vector< std::map<Int_t, Int_t> > trackMap;//top node, and its tree
                fNumGT = 0;
                for(Int_t idx=0; idx < NElement; ++idx){
                        fTheta1[idx] = fPhi1[idx] = fDet1[idx] = fNPointV[idx]
                                = fModule1[idx] = fHitEmc1[idx] = fHitEmc2[idx]= 0;
                        fPositionV[idx*3] = fPositionV[idx*3+1] = fPositionV[idx*3+2] = 0;
                        fMcp4[idx*4] = fMcp4[idx*4+1] = fMcp4[idx*4+2] = fMcp4[idx*4+3] = 0;
                }
                for(Int_t itrack = 0; itrack < nTrack; ++itrack){
                        PndMCTrack* pt1 =((PndMCTrack*)fMcTrackArray->At(itrack));
                        if(!pt1) break;
                        if(pt1->IsGeneratorCreated()){
                                pt1->Print(itrack);
                                std::map<Int_t, Int_t> newMap;
                                newMap.insert(std::pair<Int_t,Int_t>(itrack, pt1->GetMotherID()));
                                trackMap.push_back(newMap);
                                isFill = kTRUE;
                                TLorentzVector p4(pt1->GetMomentum(), pt1->Get4Momentum().E());
                                memcpy(&fMcp4[itrack*4], &p4[0], 4*sizeof(Double_t)); 
                                fTheta1[itrack]    = p4.Theta();
                                fPhi1[itrack]      = p4.Phi();
                                fHitEmc1[itrack]   = fMcTrackArray->GetEntriesFast();
                                fHitEmc2[itrack]   = pt1->GetNPoints(kEMC);
                                fNPointV[itrack]   = pt1->GetNPoints(kEMC);
                                //fDet1[itrack] = -1;
                                //for(Int_t i=0;i<15;++i){
                                //      if(detIds[i] == kEMC){
                                //              if(pt1->GetNPoints(detIds[i])>0)
                                //              {fDet1[itrack] = detIds[i];break;}
                                //      }else{
                                //              if(pt1->GetNPoints(detIds[i])>0)
                                //              {fDet1[itrack] = detIds[i];break;}
                                //      }
                                //}
                        }else{
                                if( itrack < 6){
                                        pt1->Print(itrack);
                                        cout<<"MC points #"<<pt1->GetNPoints(kEMC)<<", #"<<trackMap.size()<<endl;
                                }
                                for(size_t im=0; im < trackMap.size(); ++im){
                                        std::map<Int_t, Int_t>& theMap = trackMap[im];
                                        if(theMap.find(pt1->GetMotherID()) != theMap.end()){
                                                theMap.insert(std::pair<Int_t,Int_t>(itrack, pt1->GetMotherID()));
                                                fNPointV[im]  += pt1->GetNPoints(kEMC);
                                                break;
                                        }
                                }
                                if(pt1->GetNPoints(kEMC) > 0)  ++ fNumGT;
                        }
                }
                if(isFill)
                        tMcTruth->Fill();
                if(fVerbose>0)
                        cout<<"Save MCTruth data done"<<endl;
        }
        if(fSavePidCharged){
                for(Int_t idx=0; idx < NElement; ++idx){
                        fchar[idx] = fmcidx[idx] = 0;
                        fmom[idx] = 0.;
                        fpos[idx*3] = fpos[idx*3+1] = fpos[idx*3+2] = 0;
                        fp4[idx*4] = fp4[idx*4+1] = fp4[idx*4+2] = fp4[idx*4+3] = 0;
                }
                Int_t nCand = fChargedCand->GetEntriesFast() > NElement ? NElement : fChargedCand->GetEntriesFast();
                for(Int_t ic=0;ic < nCand; ++ic){
                        PndPidCandidate* theCand = (PndPidCandidate*)fChargedCand->At(ic);
                        fchar[ic] = theCand->GetCharge();
                        fmom[ic]  = theCand->GetMomentum().Mag();
                        fmcidx[ic] = theCand->GetMcIndex();
                        TVector3 pos = theCand->GetPosition();
                        memcpy(&fpos[ic*3], &pos[0], 3*sizeof(Double_t));
                        TLorentzVector p4 = theCand->GetLorentzVector();
                        memcpy(&fp4[ic*4], &p4[0], 4*sizeof(Double_t)); 
                }
                tCharged->Fill();
                if(fVerbose>0)
                        cout<<"Save "<<nCand<<" charged candidates done"<<endl;
        }
        if(fSavePidNeutral){
                Int_t nCand = fNeutralCand->GetEntriesFast() > NElement ? NElement : fNeutralCand->GetEntriesFast();
                for(Int_t idx=0; idx < NElement; ++idx){
                        fchar[idx] = fmcidx[idx] =  fnbump[idx] = fndigi[idx] = fmod[idx] = fcluIdx[idx] = 0;
                        frawe[idx] = fcale[idx] = fe25[idx] = fe9[idx] = 0;
                        fpos[idx*3] = fpos[idx*3+1] = fpos[idx*3+2] = 0;
                        fp4[idx*4] = fp4[idx*4+1] = fp4[idx*4+2] = fp4[idx*4+3] = 0;
                }
                for(Int_t ic=0;ic < nCand; ++ic){
                        PndPidCandidate* theCand = (PndPidCandidate*)fNeutralCand->At(ic);
                        fchar[ic] = theCand->GetCharge();
                        frawe[ic] = theCand->GetEmcRawEnergy();
                        fcale[ic] = theCand->GetEmcCalEnergy();
                        fndigi[ic]= theCand->GetEmcNumberOfCrystals();
                        fnbump[ic] = theCand->GetEmcNumberOfBumps();
                        fmod[ic]  = theCand->GetEmcModule();
                        fcluIdx[ic] = theCand->GetEmcIndex();
                        fe25[ic]  = theCand->GetEmcClusterE25();
                        fe9[ic]  = theCand->GetEmcClusterE9();
                        TVector3 pos = theCand->GetPosition();
                        memcpy(&fpos[ic*3], &pos[0], 3*sizeof(Double_t));
                        TLorentzVector p4 = theCand->GetLorentzVector();
                        memcpy(&fp4[ic*4], &p4[0], 4*sizeof(Double_t)); 
                }
                tNeutral->Fill();
                if(fVerbose>0)
                        cout<<"Save "<<nCand<<" neutral candidates done"<<endl;
        }
        if(fSaveJpsiTo3g && fChargedCand->GetEntriesFast() == 0){
                Int_t nCand = fNeutralCand->GetEntriesFast();
                if(fVerbose>0) 
                        cout<<"begin write #"<<nCand<<" candidates."<<endl;
                std::vector<key> goodCand;
                for(Int_t ic=0;ic<nCand;++ic){
                        PndPidCandidate* theCand = (PndPidCandidate*)fNeutralCand->At(ic);
                        if(theCand->GetEnergy() > 0.02)
                                goodCand.push_back(key(ic, theCand->GetEnergy()));
                }
                if(goodCand.size() >= 3){
                        sort(goodCand.begin(), goodCand.end(), std::greater<key>());
                        PndPidCandidate* cand1 = (PndPidCandidate*)fNeutralCand->At(goodCand[0].iHit);
                        fenergy1 = cand1->GetEnergy();
                        femcraw1 = cand1->GetEmcRawEnergy();
                        fndigi1 = cand1->GetEmcNumberOfCrystals();
                        fmod1   = cand1->GetEmcModule();
                        fMcTrack1 = cand1->GetMcIndex();
                        fe251   = cand1->GetEmcClusterE25();
                        flat1   = cand1->GetEmcClusterLat();
                        fz201   = cand1->GetEmcClusterZ20();
                        fz531   = cand1->GetEmcClusterZ53();
                        TLorentzVector p4 = cand1->GetLorentzVector();
                        TVector3 pos = cand1->GetPosition();
                        memcpy(&f4p1[0], &p4[0], 4*sizeof(Double_t));
                        memcpy(&fpos1[0], &pos[0], 3*sizeof(Double_t));
                        //f4p1    = cand1->GetLorentzVector();
                        //fpos1  = cand1->GetPosition();
                        PndPidCandidate* cand2 = (PndPidCandidate*)fNeutralCand->At(goodCand[1].iHit);
                        fenergy2 = cand2->GetEnergy();
                        femcraw2 = cand2->GetEmcRawEnergy();
                        fndigi2 = cand2->GetEmcNumberOfCrystals();
                        fmod2   = cand2->GetEmcModule();
                        fe252   = cand2->GetEmcClusterE25();
                        flat2   = cand2->GetEmcClusterLat();
                        fz202   = cand2->GetEmcClusterZ20();
                        fz532   = cand2->GetEmcClusterZ53();
                        fMcTrack2 = cand2->GetMcIndex();
                        p4 = cand2->GetLorentzVector();
                        pos = cand2->GetPosition();
                        memcpy(&f4p2[0], &p4[0], 4*sizeof(Double_t));
                        memcpy(&fpos2[0], &pos[0], 3*sizeof(Double_t));
                        //f4p2    = cand2->GetLorentzVector();
                        //fpos2  = cand2->GetPosition();
                        PndPidCandidate* cand3 = (PndPidCandidate*)fNeutralCand->At(goodCand[2].iHit);
                        fenergy3 = cand3->GetEnergy();
                        femcraw3 = cand3->GetEmcRawEnergy();
                        fndigi3 = cand3->GetEmcNumberOfCrystals();
                        fmod3   = cand3->GetEmcModule();
                        fe253   = cand3->GetEmcClusterE25();
                        flat3   = cand3->GetEmcClusterLat();
                        fz203   = cand3->GetEmcClusterZ20();
                        fz533   = cand3->GetEmcClusterZ53();
                        fMcTrack3 = cand3->GetMcIndex();
                        p4 = cand3->GetLorentzVector();
                        pos = cand3->GetPosition();
                        memcpy(&f4p3[0], &p4[0], 4*sizeof(Double_t));
                        memcpy(&fpos3[0], &pos[0], 3*sizeof(Double_t));
                        //f4p3    = cand3->GetLorentzVector();
                        //fpos3  = cand3->GetPosition();
                        fNumberOfGoodPhoton = goodCand.size();
                        fExternalEnergy = 0.;
                        for(size_t jc=3;jc<goodCand.size();++jc){
                                PndPidCandidate* cand = (PndPidCandidate*)fNeutralCand->At(goodCand[jc].iHit);
                                fExternalEnergy += cand->GetEnergy();
                        }
                        memset(&fMc4p1[0], 0, 4*sizeof(Double_t));
                        memset(&fMc4p2[0], 0, 4*sizeof(Double_t));
                        memset(&fMc4p3[0], 0, 4*sizeof(Double_t));
                        if(fMcTrackArray){
                                if(fMcTrack1 >= 0){
                                        PndMCTrack* p1 = (PndMCTrack*)fMcTrackArray->At(fMcTrack1);
                                        p4 = TLorentzVector(p1->GetMomentum(), p1->Get4Momentum().E());
                                        memcpy(&fMc4p1[0], &p4[0], 4*sizeof(Double_t)); 
                                }
                                if(fMcTrack2 >= 0){
                                        PndMCTrack* p2 = (PndMCTrack*)fMcTrackArray->At(fMcTrack2);
                                        p4 = TLorentzVector(p2->GetMomentum(), p2->Get4Momentum().E());
                                        memcpy(&fMc4p2[0], &p4[0], 4*sizeof(Double_t)); 
                                }
                                if(fMcTrack3 >= 0){
                                        PndMCTrack* p3 = (PndMCTrack*)fMcTrackArray->At(fMcTrack3);
                                        p4 = TLorentzVector(p3->GetMomentum(), p3->Get4Momentum().E());
                                        memcpy(&fMc4p3[0], &p4[0], 4*sizeof(Double_t)); 
                                }
                        }
                        tAna3g->Fill();
                }
        }


        if(fSaveTask){
                fNumOfHits    = fEmcHitsArray->GetEntriesFast();
                fNumWaveform  = fWaveformArray->GetEntriesFast();
                fNumDigi      = fDigiArray->GetEntriesFast();
                fNumSharedDigi= fSharedDigiArray->GetEntriesFast();
                fNumCluster   = fClusterArray->GetEntriesFast();
                fNumBump      = fBumpArray->GetEntriesFast();
                fNumMcTrack   = fMcTrackArray->GetEntriesFast();
                tTask->Fill();
                if(fVerbose>0)
                        cout<<"Save task data done"<<endl;
        }

        if (fVerbose>2){
                timer.Stop();
                Double_t rtime = timer.RealTime();
                Double_t ctime = timer.CpuTime();
                cout << "PndEmcAnalysis, Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
        }
        //fEvtHeaderArray->Delete();
        fEventNo ++;
        cout<<"================================================================="<<endl;
}

void PndEmcAnalysis::FinishTask ( )

virtual

Definition at line 895 of file PndEmcAnalysis.cxx.

References evtMap, evtset, fEventNo, fPileupWave, fRootFile, fSaveBump, fSaveDigi, fSaveHits, fSaveJpsiTo3g, fSaveMcTruth, fSavePidCharged, fSavePidNeutral, fSaveReco, fSaveTask, fSaveWave, fStoreRooTFile, fTimeOrderedDigi, fTotWave, pevtset, tAna3g, tBump, tCharged, tDigi, tEvtPileup, tHit, tMcTruth, tNeutral, totNumOfBump, totNumOfDigi, totNumOfEvents, totNumOfPileupEvents, totNumOfSharedDigi, totNumOfSubWave, totNumOfWave, tRec, tTask, tTaskPileup, and tWave.

{

        std::cout<<"*********************************************************"<<std::endl;
        std::cout<<"PndEmcAnalysis::FinishTask"<<std::endl;
        std::cout<<"*********************************************************"<<std::endl;
        std::cout<<"Total waveforms# "<<totNumOfWave<<std::endl;
        std::cout<<"Total subwaveforms# "<<totNumOfSubWave<<std::endl;
        std::cout<<"Total fpga digis# "<<totNumOfDigi<<std::endl;
        std::cout<<"Total shared digis# "<<totNumOfSharedDigi<<std::endl;
        std::cout<<"Total bumps# "<<totNumOfBump<<std::endl;
        //std::cout<<"Digis to events they belong to#"<<matchCount<<std::endl;
        //std::cout<<"Digis to events they dont belong to#"<<unmatchCount<<std::endl;
        std::cout<<"*********************************************************"<<std::endl;
        if(fTimeOrderedDigi){
                std::map<Int_t, std::pair<Int_t,Int_t> >::iterator it = evtMap.begin();
                for(; it != evtMap.end(); ++it){
                        fEventNo = (*it).first;
                        fTotWave = (*it).second.first;
                        fPileupWave = (*it).second.second;
                        tEvtPileup->Fill();
                }
                totNumOfEvents = evtset.size();
                totNumOfPileupEvents = pevtset.size();
                tTaskPileup->Fill();
        }


        if(fStoreRooTFile){
                TFile* gOldFile = gFile;
                gFile = fRootFile;
                gFile->cd();
                if(fSaveHits){
                        cout<<"write "<<tHit->GetTitle()<<endl;
                        tHit->Write();
                }
                if(fSaveWave){
                        cout<<"write "<<tWave->GetTitle()<<endl;
                        tWave->Write();
                }
                if(fSaveDigi){
                        cout<<"write "<<tDigi->GetTitle()<<endl;
                        tDigi->Write();
                }
                if(fSaveBump){
                        cout<<"write "<<tBump->GetTitle()<<endl;
                        tBump->Write();
                }
                if(fSaveTask){
                        cout<<"write "<<tTask->GetTitle()<<endl;
                        tTask->Write();
                }
                if(fSaveReco){
                        cout<<"write "<<tRec->GetTitle()<<endl;
                        tRec->Write();
                }
                if(fSaveMcTruth){
                        cout<<"write "<<tMcTruth->GetTitle()<<endl;
                        tMcTruth->Write();
                }
                if(fTimeOrderedDigi){
                        cout<<"write "<<tTaskPileup->GetTitle()<<endl;
                        cout<<"write "<<tEvtPileup->GetTitle()<<endl;
                        tTaskPileup->Write();
                        tEvtPileup->Write();
                }
                if(fSavePidCharged){
                        cout<<"write "<<tCharged->GetTitle()<<endl;
                        tCharged->Write();
                }
                if(fSavePidNeutral){
                        cout<<"write "<<tNeutral->GetTitle()<<endl;
                        tNeutral->Write();
                }
                if(fSaveJpsiTo3g){
                        cout<<"write "<<tAna3g->GetTitle()<<endl;
                        tAna3g->Write();
                }

                gFile = gOldFile;
        }

        evtset.clear();
        pevtset.clear();
        evtMap.clear();
}

InitStatus PndEmcAnalysis::Init ( )

virtual

Virtual method Init

Definition at line 116 of file PndEmcAnalysis.cxx.

References evtMap, evtset, f4p1, f4p2, f4p3, fBaseline, fBumpArray, fBumpEnergy, fBumpTime1, fBumpTime2, fBumpTime3, fcale, fchar, fChargedCand, fcluIdx, fClusterArray, fClusterIdx, fDigiArray, fDigiEnergy, fDigiMapSize, fDigiTime, fDistance, fe25, fe251, fe252, fe253, fe9, fEmcHitsArray, femcraw1, femcraw2, femcraw3, fEnergy, fenergy1, fenergy2, fenergy3, fEnergyC1, fEnergyC2, fEventNo, fEvtTime, fExternalEnergy, fFunctor, fHitEmc1, fHitEmc2, fHitEnergy, fHitIndex, fHitTime, fileName, fLat1, flat1, flat2, flat3, fMc4p1, fMc4p2, fMc4p3, fmcidx, fMcp4, fMcTrack1, fMcTrack2, fMcTrack3, fMcTrackArray, fmod, fmod1, fmod2, fmod3, fModule, fModule1, fmom, fnbump, fndigi, fndigi1, fndigi2, fndigi3, fNeutralCand, fNPoint, fNPointV, fNumberOfGoodPhoton, fNumBump, fNumCluster, fNumDigi, fNumGT, fNumMcTrack, fNumOfDismatch, fNumOfHits, fNumOfMatch, fNumOfSeed, fNumSharedDigi, fNumWaveform, fp4, fPDGCode, fPhi, fPhi1, fPileupCount, fPileupWave, fpos, fpos1, fpos2, fpos3, fPosition, fPositionV, frawe, fRestEnergy1, fRestEnergy2, fRestEnergy3, fRootFile, fSaveBump, fSaveDigi, fSaveHits, fSaveJpsiTo3g, fSaveMcTruth, fSavePidCharged, fSavePidNeutral, fSaveReco, fSaveTask, fSaveWave, fSeedPosition, fSeedTime, fSharedDigiArray, fStoreRooTFile, fTheta, fTheta1, fTimeOrderedDigi, fTotWave, fWaveformArray, fWaveIdx, fWaveLen, fWavePeak, fWaveTime, fXcor, fYcor, fZ201, fz201, fz202, fz203, fZ531, fz531, fz532, fz533, pevtset, tAna3g, tBump, tCharged, tDigi, tEvtPileup, tHit, tMcTruth, tNeutral, totDigiEnergy, totNumOfBump, totNumOfDigi, totNumOfEvents, totNumOfPileupEvents, totNumOfSharedDigi, totNumOfSubWave, totNumOfWave, totSharedDigiEnergy, tRec, tTask, tTaskPileup, and tWave.

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

        fEmcHitsArray = (TClonesArray*) ioman->GetObject("EmcHit");
        if ( ! fEmcHitsArray ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No EmcHit array!" << endl;
                fSaveHits = kFALSE;
        }
        // Get input array
        if(fTimeOrderedDigi){
                fWaveformArray = (TClonesArray*) ioman->GetObject("EmcSortedWaveform");
        }else{
                fWaveformArray = (TClonesArray*) ioman->GetObject("EmcWaveform");
        }
        //fEvtHeaderArray = (TClonesArray*) ioman->GetObject("EventHeader.");
        if ( ! fWaveformArray ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No PndEmcWaveform array!" << endl;
                fSaveWave = kFALSE;
        }

        fDigiArray = (TClonesArray*) ioman->GetObject("EmcDigi");
        if ( ! fDigiArray ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No PndEmcDigi array!" << endl;
                fSaveDigi = kFALSE;
        }
        fSharedDigiArray = (TClonesArray*) ioman->GetObject("EmcSharedDigi");
        if ( ! fSharedDigiArray ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No PndEmcSharedDigi array!" << endl;
                fSaveDigi = kFALSE;
        }
        fMcTrackArray = (TClonesArray*) ioman->GetObject("MCTrack");
        if ( ! fMcTrackArray ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No MCTrack array!" << endl;
                fSaveMcTruth = kFALSE;
        }
        fClusterArray = (TClonesArray*) ioman->GetObject("EmcCluster");
        if ( ! fClusterArray ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No PndEmcCluster array!" << endl;
                //fSaveBump = kFALSE;
        }
        fBumpArray = (TClonesArray*) ioman->GetObject("EmcBump");
        if ( ! fBumpArray ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No PndEmcBump array!" << endl;
                fSaveBump = kFALSE;
        }
        fChargedCand = (TClonesArray*)ioman->GetObject("PidChargedCand");
        if ( ! fChargedCand ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No PidChargedCand array!" << endl;
                fSavePidCharged = kFALSE;
        }
        fNeutralCand = (TClonesArray*)ioman->GetObject("PidNeutralCand");
        if ( ! fNeutralCand ) {
                cout << "-W- PndEmcAnalysis::Init: "
                        << "No PidNeutralCand array!" << endl;
                fSavePidNeutral = kFALSE;
        }


        fFunctor = new StopTime();
        fStoreRooTFile = fSaveDigi || fSaveWave || fSaveBump || fSaveTask || fTimeOrderedDigi || fSaveHits || fSaveMcTruth || fSaveReco;

        if(fStoreRooTFile)
                fRootFile = new TFile(fileName,"RECREATE");

        if(fSaveHits){
                tHit = new TTree("hit", "hit data");
                tHit->Branch("ie"    ,&fEventNo                         ,"EventNO/I");
                tHit->Branch("ht"    ,&fHitTime                         ,"HitTime/D");
                tHit->Branch("mod"   ,&fModule                          ,"Module/I");
                tHit->Branch("he"    ,&fHitEnergy     ,"HitEnergy/D");
                tHit->Branch("theta" ,&fTheta         ,"Theta/D");
                tHit->Branch("phi"   ,&fPhi           ,"Phi/D");
                tHit->Branch("np"    ,&fNPoint        ,"NumOfPoint/I");
                tHit->Branch("xcor"  ,&fXcor              ,"XCoor/I");
                tHit->Branch("ycor"  ,&fYcor              ,"XCoor/I");
        }

        if(fSaveWave){
                tWave = new TTree("wave", "wave data");
                tWave->Branch("et"    ,&fEvtTime                        ,"EventTime/D");
                tWave->Branch("mod"   ,&fModule                                 ,"Module/I");
                tWave->Branch("hitid" ,&fHitIndex       ,"HitIndex/I");
                tWave->Branch("wavet" ,&fWaveTime                               ,"WaveTime/D");
                tWave->Branch("peak"  ,&fWavePeak               ,"WavePeak/D");
                tWave->Branch("wl"    ,&fWaveLen        ,"WaveLength/I");
                tWave->Branch("iw"    ,&fWaveIdx                                ,"WaveIndex/I");
                tWave->Branch("ie"    ,&fEventNo                                ,"EventNO/I");
                tWave->Branch("baseE" ,&fBaseline       ,"BaseLine/D");
                tWave->Branch("xcor"  ,&fXcor               ,"XCoor/I");
                tWave->Branch("ycor"  ,&fYcor               ,"YCoor/I");
                tWave->Branch("wpc"   ,&fPileupCount    ,"Pileup/I");
                tWave->Branch("he"    ,&fHitEnergy       ,"HitEnergy/D");
        }
        if(fSaveDigi){
                tDigi = new TTree("digi", "digi data");
                tDigi->Branch("digit" ,&fDigiTime                 ,"DigiTime/D");
                tDigi->Branch("et"    ,&fEvtTime                        ,"EventTime/D");
                tDigi->Branch("mod"   ,&fModule                                 ,"Module/I");
                tDigi->Branch("wt"    ,&fWaveTime                               ,"WaveTime/D");
                tDigi->Branch("hite"  ,&fHitEnergy                      ,"HitEnergy/D");
                tDigi->Branch("e"     ,&fDigiEnergy         ,"Energy/D");
                tDigi->Branch("r"     ,&fDistance                               ,"Path/D");
                tDigi->Branch("iw"    ,&fWaveIdx                                ,"WaveIndex/I");
                tDigi->Branch("ie"    ,&fEventNo                                ,"EventNO/I");
                tDigi->Branch("pos"   ,&fPosition[0]    ,"Position[3]/D");
        }
        if(fSaveBump){
                tBump = new TTree("bump", "bump data");
                tBump->Branch("et"    ,&fEvtTime                        ,"EventTime/D");
                tBump->Branch("ie"    ,&fEventNo                                ,"EventNO/I");
                tBump->Branch("mod"   ,&fModule                                 ,"Module/I");
                tBump->Branch("bt1"   ,&fBumpTime1                      ,"WeightedTimeByTimeError/D");
                tBump->Branch("bt2"   ,&fBumpTime2                      ,"WeightedTimeByTimeErrorAndEnergy/D");
                tBump->Branch("bt3"   ,&fBumpTime3                      ,"WeightedTimeByEnergy/D");
                tBump->Branch("seedt" ,&fSeedTime                               ,"SeedTime/D");
                tBump->Branch("e"     ,&fBumpEnergy             ,"BumpEnergy/D");
                tBump->Branch("pos"   ,&fPosition[0]      ,"Position[3]/D");
                tBump->Branch("seedpos",&fSeedPosition[0]         ,"SeedPosition[3]/D");
                tBump->Branch("ndigi" ,&fNumDigi        ,"NumOfDigi/I");
                tBump->Branch("nbump" ,&fNumBump                                ,"NumofBump/I");
                tBump->Branch("ic"    ,&fClusterIdx         ,"ClusterIndex/I");
                tBump->Branch("nseed" ,&fNumOfSeed      ,"NumOfSeed/I");
                tBump->Branch("nmap"  ,&fDigiMapSize      ,"NumDigiMap/I");
        }
        if(fSaveTask){
                tTask = new TTree("task", "task data");
                tTask->Branch("ie"      ,&fEventNo                       ,"EventNO/I");
                tTask->Branch("nhits"   ,&fNumOfHits     ,"NumOfEmcHit/I");
                tTask->Branch("nwave"   ,&fNumWaveform   ,"NumOfWaveform/I");
                tTask->Branch("ndigi"   ,&fNumDigi       ,"NumOfDigi/I");
                tTask->Branch("sdigi"   ,&fNumSharedDigi ,"NumOfSharedDigi/I");
                tTask->Branch("nclus"   ,&fNumCluster    ,"NumOfCluster/I");
                tTask->Branch("nbump"   ,&fNumBump       ,"NumOfBump/I");
                tTask->Branch("ntrack"  ,&fNumMcTrack    ,"NumOfMcTrack/I");
                tTask->Branch("totde"   ,&totDigiEnergy   ,"TotalDigiEnergy/D");
                tTask->Branch("totse"   ,&totSharedDigiEnergy,"TotalSharedDigiEnergy/D");
        }
        if(fTimeOrderedDigi){
                tEvtPileup= new TTree("evtpu","event pileup data");
                tEvtPileup->Branch("evtNo"   ,&fEventNo    ,"evtNo/I");
                tEvtPileup->Branch("totw"    ,&fTotWave    ,"totw/I");
                tEvtPileup->Branch("pilew"   ,&fPileupWave ,"pilew/I");

                tTaskPileup= new TTree("taskpu","task pileup data");
                tTaskPileup->Branch("nw"      ,&totNumOfWave       ,"TotNumOfWave/I");
                tTaskPileup->Branch("nsw"     ,&totNumOfSubWave    ,"TotNumOfSubWave/I");
                tTaskPileup->Branch("nd"      ,&totNumOfDigi       ,"TotNumOfDigi/I");
                tTaskPileup->Branch("nsd"     ,&totNumOfSharedDigi ,"TotNumOfSharedDigi/I");
                tTaskPileup->Branch("nevt"    ,&totNumOfEvents     ,"TotNumOfEvents/I");
                tTaskPileup->Branch("npevt"   ,&totNumOfPileupEvents,"TotNumOfPileupEvents/I");
                tTaskPileup->Branch("match"   ,&fNumOfMatch        ,"NumOfMatch/I");
                tTaskPileup->Branch("mismatch",&fNumOfDismatch     ,"NumOfDismatch/I");
        }
        if(fSaveReco){
                tRec = new TTree("ResA", "Reconstruction");
                tRec->Branch("E1"    ,&fEnergy[0]               ,"E1[30]/D");
                tRec->Branch("E1C"   ,&fEnergyC1[0]             ,"E1C[30]/D");
                tRec->Branch("E2C"   ,&fEnergyC2[0]                     ,"E2C[30]/D");
                tRec->Branch("p1"    ,&fp4[0]                               ,"p1[120]/D");
                tRec->Branch("Mcp1"  ,&fMcp4[0]                         ,"Mcp1[120]/D");
                tRec->Branch("pdg"  ,&fPDGCode[0]                                       ,"pdg[30]/I");
                tRec->Branch("Z201"  ,&fZ201[0]                                 ,"Z201[30]/D");
                tRec->Branch("Z531"  ,&fZ531[0]                                 ,"Z531[30]/D");
                tRec->Branch("Lat1"  ,&fLat1[0]                                 ,"Lat1[30]/D");
                tRec->Branch("pos"    ,&fPositionV[0]   ,"Position[30]/D");
                tRec->Branch("Mod1"  ,&fModule1[0]                      ,"Mod1[30]/I");
                tRec->Branch("Theta1",&fTheta1[0]                               ,"Theta1[30]/D");
                tRec->Branch("Phi1"  ,&fPhi1[0]                                 ,"Phi1[30]/D");
                tRec->Branch("ExtE1" ,&fRestEnergy1             ,"ExtE1/D");
                tRec->Branch("ExtE2" ,&fRestEnergy2             ,"ExtE2/D");
                tRec->Branch("ExtE3" ,&fRestEnergy3             ,"ExtE3/D");
                tRec->Branch("NBump" ,&fNumBump                 ,"NBump/I");
        }
        if(fSaveMcTruth){
                tMcTruth = new TTree("Res1A", "Mc Truth");
                tMcTruth->Branch("Mcp1"         ,&fMcp4[0]              ,"Mcp1[120]/D");
                tMcTruth->Branch("Theta1" ,&fTheta1[0]          ,"Theta1[30]/D");
                tMcTruth->Branch("Phi1"         ,&fPhi1[0]                      ,"Phi1[30]/D");
                //tMcTruth->Branch("Det1"               ,&fDet1[0]                      ,"Det1[30]/I");
                tMcTruth->Branch("isHit1"       ,&fHitEmc1[0]           ,"SizeOfMcTrackArray[30]/I");
                tMcTruth->Branch("isHit2"       ,&fHitEmc2[0]           ,"NumHitOfPrimaryTrack[30]/I");
                tMcTruth->Branch("isHit3"       ,&fNPointV[0]           ,"TotalNumHit[30]/I");
                tMcTruth->Branch("isHit4"       ,&fNumGT                        ,"NumOfSecondaryHit/I");
                //tMcTruth->Branch("X1"                 ,&fX1                                   ,"X1/D");
                //tMcTruth->Branch("Y1"                 ,&fY1                                   ,"Y1/D");
                //tMcTruth->Branch("Z1"                 ,&fZ1                                   ,"Z1/D");
        }
        if(fSavePidCharged)
        {
                tCharged= new TTree("char", "charged cand");
                tCharged->Branch("char"         ,&fchar[0]      ,"char[30]/I");
                tCharged->Branch("p4"           ,&fp4[0]                        ,"p4[120]/D");
                tCharged->Branch("pos"          ,&fpos[0]                       ,"pos[90]/D");
                tCharged->Branch("mom"          ,&fmom[0]                       ,"mom[30]/D");
                tCharged->Branch("mcIdx"        ,&fmcidx[0]             ,"mctrackid[30]/I");

        }
        if(fSavePidNeutral)
        {
                tNeutral= new TTree("neut", "neutral cand");
                tNeutral->Branch("char"         ,&fchar[0]      ,"char[30]/I");
                tNeutral->Branch("mcIdx"        ,&fmcidx[0]             ,"mctrackid[30]/I");
                tNeutral->Branch("p4"           ,&fp4[0]                        ,"p4[120]/D");
                tNeutral->Branch("pos"          ,&fpos[0]                       ,"pos[90]/D");
                tNeutral->Branch("rawe"         ,&frawe[0]              ,"rawe[30]/D");
                tNeutral->Branch("cale"         ,&fcale[0]              ,"cale[30]/D");
                tNeutral->Branch("ndigi"        ,&fndigi[0]             ,"ndigi[30]/I");
                tNeutral->Branch("mod"          ,&fmod[0]                       ,"EmcModule[30]/I");
                tNeutral->Branch("cluIdx"       ,&fcluIdx[0]    ,"clusterIdx[30]/I");
                tNeutral->Branch("e25"          ,&fe25[0]                       ,"e25[30]/D");
                tNeutral->Branch("e9"                   ,&fe9[0]                        ,"e9[30]/D");
                tNeutral->Branch("nbump"        ,&fnbump[0]             ,"nbump[30]/I");
        }
        if(fSaveJpsiTo3g){
                tAna3g = new TTree("Ana", "Jpsi-->3photon");
                tAna3g->Branch("ng"               ,&fNumberOfGoodPhoton         ,"GoodPhoton/I");
                tAna3g->Branch("eE"         ,&fExternalEnergy               ,"ExtralEnergy/D");
                tAna3g->Branch("4p1"    ,&f4p1[0]               ,"Photon1[4]/D");
                tAna3g->Branch("mc4p1"  ,&fMc4p1[0]             ,"Truth1[4]/D");
                tAna3g->Branch("pos1"   ,&fpos1[0]              ,"pos1[3]/D");
                tAna3g->Branch("e1"     ,&fenergy1              ,"e1/D");
                tAna3g->Branch("rawe1"  ,&femcraw1              ,"rawe1/D");
                tAna3g->Branch("nd1"    ,&fndigi1               ,"NumOfDigi1/I");
                tAna3g->Branch("mcTrk1" ,&fMcTrack1            ,"McTrackID1/I");
                tAna3g->Branch("mod1"   ,&fmod1                ,"mod1/I");
                tAna3g->Branch("e251"   ,&fe251                ,"e251/D");
                tAna3g->Branch("lat1"   ,&flat1                ,"lat1/D");
                tAna3g->Branch("z201"   ,&fz201                ,"z201/D");
                tAna3g->Branch("z531"   ,&fz531                ,"z531/D");

                tAna3g->Branch("4p2"    ,&f4p2[0]               ,"Photon2[4]/D");
                tAna3g->Branch("mc4p2"  ,&fMc4p2[0]             ,"Truth2[4]/D");
                tAna3g->Branch("pos2"   ,&fpos2[0]              ,"pos2[3]/D");
                tAna3g->Branch("e2"     ,&fenergy2              ,"e2/D");
                tAna3g->Branch("rawe2"  ,&femcraw2              ,"rawe2/D");
                tAna3g->Branch("nd2"    ,&fndigi2               ,"NumOfDigi2/I");
                tAna3g->Branch("mcTrk2" ,&fMcTrack2            ,"McTrackID2/I");
                tAna3g->Branch("mod2"   ,&fmod2                ,"mod2/I");
                tAna3g->Branch("e252"   ,&fe252                ,"e252/D");
                tAna3g->Branch("lat2"   ,&flat2                ,"lat2/D");
                tAna3g->Branch("z202"   ,&fz202                ,"z202/D");
                tAna3g->Branch("z532"   ,&fz532                ,"z532/D");

                tAna3g->Branch("4p3"    ,&f4p3[0]               ,"Photon3[4]/D");
                tAna3g->Branch("mc4p3"  ,&fMc4p3[0]             ,"Truth3[4]/D");
                tAna3g->Branch("pos3"   ,&fpos3[0]              ,"pos3[3]/D");
                tAna3g->Branch("e3"     ,&fenergy3              ,"e3/D");
                tAna3g->Branch("rawe3"  ,&femcraw3              ,"rawe3/D");
                tAna3g->Branch("nd3"    ,&fndigi3               ,"NumOfDigi3/I");
                tAna3g->Branch("mcTrk3" ,&fMcTrack3            ,"McTrackID3/I");
                tAna3g->Branch("mod3"   ,&fmod3                ,"mod3/I");
                tAna3g->Branch("e253"   ,&fe253                ,"e253/D");
                tAna3g->Branch("lat3"   ,&flat3                ,"lat3/D");
                tAna3g->Branch("z203"   ,&fz203                ,"z203/D");
                tAna3g->Branch("z533"   ,&fz533                ,"z533/D");

        }


        evtset.clear();
        pevtset.clear();
        evtMap.clear();

        totNumOfWave         = 0;
        totNumOfSubWave      = 0;
        totNumOfDigi         = 0;
        totNumOfSharedDigi   = 0;
        totNumOfBump         = 0;
        totNumOfEvents       = 0;
        totNumOfPileupEvents = 0;
        fNumOfMatch          = 0;
        fNumOfDismatch       = 0;

        //vector resize
        //      const Int_t Num = 1;
        //      fEnergy.resize(Num);
        //      fEnergyC1.resize(Num);
        //      fEnergyC2.resize(Num);
        //      fPositionV.resize(3*Num);
        //      fp4.resize(4*Num);
        //      fMcp4.resize(4*Num);
        //      fZ201.resize(Num);
        //      fZ531.resize(Num);
        //      fLat1.resize(Num);
        //      fModule1.resize(Num);
        //      fTheta1.resize(Num);
        //      fPhi1.resize(Num);
        //      fDet1.resize(Num);
        //      fHitEmc1.resize(Num);
        //      fHitEmc2.resize(Num);
        //      fNPointV.resize(Num);

        cout << "-I- PndEmcAnalysis: Intialization successfull" << endl;
        return kSUCCESS;
}

PndEmcAnalysis& PndEmcAnalysis::operator= ( const PndEmcAnalysis &  )

inlineprivate

Definition at line 105 of file PndEmcAnalysis.h.

{return *this;};

void PndEmcAnalysis::RunTimeBased ( )

inline

Definition at line 70 of file PndEmcAnalysis.h.

{fTimeOrderedDigi = kTRUE;}

void PndEmcAnalysis::SetOutput ( TString  name )

inline

Definition at line 68 of file PndEmcAnalysis.h.

References name.

{ fileName = name; }

void PndEmcAnalysis::SetParContainers ( )

private

Definition at line 876 of file PndEmcAnalysis.cxx.

References 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 digitisation parameter container
        //fDigiPar = (PndEmcDigiPar*) db->getContainer("PndEmcDigiPar");

        // Get Emc reconstruction parameter container
        //fRecoPar = (PndEmcRecoPar*) db->getContainer("PndEmcRecoPar");

        // Get Emc FPGA parameter container
        //fFpgaPar = (PndEmcFpgaPar*) db->getContainer("PndEmcFpgaPar");
}

Member Data Documentation

PndEmcDigiCalibrator PndEmcAnalysis::digiCalibrator

private

Definition at line 296 of file PndEmcAnalysis.h.

std::map<Int_t, std::pair<Int_t, Int_t> > PndEmcAnalysis::evtMap

private

Definition at line 114 of file PndEmcAnalysis.h.

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

std::set<Int_t> PndEmcAnalysis::evtset

private

Definition at line 116 of file PndEmcAnalysis.h.

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

Double_t PndEmcAnalysis::f4p1[4]

private

Definition at line 253 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::f4p2[4]

private

Definition at line 267 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::f4p3[4]

private

Definition at line 281 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fBaseline

private

Definition at line 144 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcAnalysis::fBumpArray

private

Definition at line 93 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fBumpEnergy

private

Definition at line 170 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fBumpTime1

private

Definition at line 166 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fBumpTime2

private

Definition at line 167 of file PndEmcAnalysis.h.

Referenced by Init().

Double_t PndEmcAnalysis::fBumpTime3

private

Definition at line 168 of file PndEmcAnalysis.h.

Referenced by Init().

Double_t PndEmcAnalysis::fcale[NElement]

private

Definition at line 237 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fchar[NElement]

private

Definition at line 231 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcAnalysis::fChargedCand

private

Definition at line 96 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fcluIdx[NElement]

private

Definition at line 240 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcAnalysis::fClusterArray

private

Definition at line 92 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fClusterIdx

private

Definition at line 175 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fDet1[NElement]

private

Definition at line 223 of file PndEmcAnalysis.h.

Referenced by Exec().

TClonesArray* PndEmcAnalysis::fDigiArray

private

Definition at line 90 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fDigiEnergy

private

Definition at line 154 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fDigiMapSize

private

Definition at line 177 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fDigiTime

private

Definition at line 149 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fDistance

private

Definition at line 155 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fe25[NElement]

private

Definition at line 241 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fe251

private

Definition at line 262 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fe252

private

Definition at line 276 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fe253

private

Definition at line 290 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fe9[NElement]

private

Definition at line 242 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcAnalysis::fEmcHitsArray

private

Input arrays

Definition at line 88 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::femcraw1

private

Definition at line 258 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::femcraw2

private

Definition at line 272 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::femcraw3

private

Definition at line 286 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fEnergy[NElement]

private

Definition at line 208 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fenergy1

private

Definition at line 257 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fEnergy1C

private

Definition at line 161 of file PndEmcAnalysis.h.

Referenced by Exec().

Double_t PndEmcAnalysis::fenergy2

private

Definition at line 271 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fEnergy2C

private

Definition at line 162 of file PndEmcAnalysis.h.

Double_t PndEmcAnalysis::fenergy3

private

Definition at line 285 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fEnergyC1[NElement]

private

Definition at line 209 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fEnergyC2[NElement]

private

Definition at line 210 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fEventNo

private

Definition at line 143 of file PndEmcAnalysis.h.

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

TClonesArray* PndEmcAnalysis::fEvtHeaderArray

private

Definition at line 95 of file PndEmcAnalysis.h.

Double_t PndEmcAnalysis::fEvtTime

private

Definition at line 136 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fExternalEnergy

private

Definition at line 251 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

BinaryFunctor* PndEmcAnalysis::fFunctor

private

Definition at line 109 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fHitEmc1[NElement]

private

Definition at line 224 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fHitEmc2[NElement]

private

Definition at line 225 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fHitEnergy

private

Definition at line 131 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fHitIndex

private

Definition at line 138 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fHitTime

private

Definition at line 130 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TString PndEmcAnalysis::fileName

private

Definition at line 118 of file PndEmcAnalysis.h.

Referenced by Init(), and PndEmcAnalysis().

Double_t PndEmcAnalysis::fLat1[NElement]

private

Definition at line 213 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::flat1

private

Definition at line 263 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::flat2

private

Definition at line 277 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::flat3

private

Definition at line 291 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fMc4p1[4]

private

Definition at line 254 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fMc4p2[4]

private

Definition at line 268 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fMc4p3[4]

private

Definition at line 282 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fmcidx[NElement]

private

Definition at line 235 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fMcp4[4 *NElement]

private

Definition at line 206 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fMcTrack1

private

Definition at line 261 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fMcTrack2

private

Definition at line 275 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fMcTrack3

private

Definition at line 289 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcAnalysis::fMcTrackArray

private

Definition at line 94 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fmod[NElement]

private

Definition at line 239 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fmod1

private

Definition at line 260 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fmod2

private

Definition at line 274 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fmod3

private

Definition at line 288 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fModule

private

Definition at line 137 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fModule1[NElement]

private

Definition at line 214 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fmom[NElement]

private

Definition at line 234 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fmom1

private

Definition at line 256 of file PndEmcAnalysis.h.

Double_t PndEmcAnalysis::fmom2

private

Definition at line 270 of file PndEmcAnalysis.h.

Double_t PndEmcAnalysis::fmom3

private

Definition at line 284 of file PndEmcAnalysis.h.

Int_t PndEmcAnalysis::fnbump[NElement]

private

Definition at line 243 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fndigi[NElement]

private

Definition at line 238 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fndigi1

private

Definition at line 259 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fndigi2

private

Definition at line 273 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fndigi3

private

Definition at line 287 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcAnalysis::fNeutralCand

private

Definition at line 97 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNPoint

private

Definition at line 134 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNPointV[NElement]

private

Definition at line 226 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumberOfGoodPhoton

private

Definition at line 250 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumBump

private

Definition at line 174 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumCluster

private

Definition at line 188 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumDigi

private

Definition at line 173 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumGT

private

Definition at line 221 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumMcTrack

private

Definition at line 190 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumOfDismatch

private

Definition at line 179 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumOfHits

private

Definition at line 184 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumOfMatch

private

Definition at line 178 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumOfSeed

private

Definition at line 176 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumSharedDigi

private

Definition at line 187 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fNumWaveform

private

Definition at line 185 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fp4[4 *NElement]

private

Definition at line 205 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fPDGCode[NElement]

private

Definition at line 215 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fPhi

private

Definition at line 133 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fPhi1[NElement]

private

Definition at line 217 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fPileupCount

private

Definition at line 147 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fPileupWave

private

Definition at line 194 of file PndEmcAnalysis.h.

Referenced by FinishTask(), and Init().

Double_t PndEmcAnalysis::fpos[NElement]

private

Definition at line 233 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fpos1[3]

private

Definition at line 255 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fpos2[3]

private

Definition at line 269 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fpos3[3]

private

Definition at line 283 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TVector3 PndEmcAnalysis::fPosition

private

Definition at line 158 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fPositionV[3 *NElement]

private

Definition at line 207 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::frawe[NElement]

private

Definition at line 236 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fRestEnergy1

private

Definition at line 218 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fRestEnergy2

private

Definition at line 219 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fRestEnergy3

private

Definition at line 220 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TFile* PndEmcAnalysis::fRootFile

private

Definition at line 119 of file PndEmcAnalysis.h.

Referenced by FinishTask(), and Init().

Bool_t PndEmcAnalysis::fSaveBump

Definition at line 76 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSaveDigi

Definition at line 75 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSaveHits

Definition at line 78 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSaveJpsiTo3g

Definition at line 83 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSaveMcTruth

Definition at line 79 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSavePidCharged

Definition at line 81 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSavePidNeutral

Definition at line 82 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSaveReco

Definition at line 80 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSaveTask

Definition at line 77 of file PndEmcAnalysis.h.

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

Bool_t PndEmcAnalysis::fSaveWave

Definition at line 74 of file PndEmcAnalysis.h.

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

TVector3 PndEmcAnalysis::fSeedPosition

private

Definition at line 159 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fSeedTime

private

Definition at line 169 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TClonesArray* PndEmcAnalysis::fSharedDigiArray

private

Definition at line 91 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Bool_t PndEmcAnalysis::fStoreRooTFile

private

Definition at line 102 of file PndEmcAnalysis.h.

Referenced by FinishTask(), and Init().

Double_t PndEmcAnalysis::fTheta

private

Definition at line 132 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fTheta1[NElement]

private

Definition at line 216 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Bool_t PndEmcAnalysis::fTimeOrderedDigi

private

set to kTRUE to use the time ordering of the output data.

Definition at line 112 of file PndEmcAnalysis.h.

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

Int_t PndEmcAnalysis::fTotWave

private

Definition at line 193 of file PndEmcAnalysis.h.

Referenced by FinishTask(), and Init().

Int_t PndEmcAnalysis::fVerbose

private

Definition at line 101 of file PndEmcAnalysis.h.

Referenced by Exec().

TClonesArray* PndEmcAnalysis::fWaveformArray

private

Definition at line 89 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fWaveIdx

private

Definition at line 142 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fWaveLen

private

Definition at line 141 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fWavePeak

private

Definition at line 140 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fWaveTime

private

Definition at line 139 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fXcor

private

Definition at line 145 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::fYcor

private

Definition at line 146 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fZ201[NElement]

private

Definition at line 211 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fz201

private

Definition at line 264 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fz202

private

Definition at line 278 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fz203

private

Definition at line 292 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fZ531[NElement]

private

Definition at line 212 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fz531

private

Definition at line 265 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fz532

private

Definition at line 279 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Double_t PndEmcAnalysis::fz533

private

Definition at line 293 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

const Int_t PndEmcAnalysis::NElement = 30

staticprivate

Definition at line 204 of file PndEmcAnalysis.h.

Referenced by Exec().

std::set<Int_t> PndEmcAnalysis::pevtset

private

Definition at line 115 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tAna3g

private

Definition at line 248 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tBump

private

Definition at line 123 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tCharged

private

Definition at line 245 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tDigi

private

Definition at line 121 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tEvtPileup

private

Definition at line 125 of file PndEmcAnalysis.h.

Referenced by FinishTask(), and Init().

TTree* PndEmcAnalysis::tHit

private

Definition at line 120 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tMcTruth

private

Definition at line 127 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tNeutral

private

Definition at line 246 of file PndEmcAnalysis.h.

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

Double_t PndEmcAnalysis::totDigiEnergy

private

Definition at line 181 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

Int_t PndEmcAnalysis::totNumOfBump

private

Definition at line 198 of file PndEmcAnalysis.h.

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

Int_t PndEmcAnalysis::totNumOfDigi

private

Definition at line 197 of file PndEmcAnalysis.h.

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

Int_t PndEmcAnalysis::totNumOfEvents

private

Definition at line 200 of file PndEmcAnalysis.h.

Referenced by FinishTask(), and Init().

Int_t PndEmcAnalysis::totNumOfPileupEvents

private

Definition at line 201 of file PndEmcAnalysis.h.

Referenced by FinishTask(), and Init().

Int_t PndEmcAnalysis::totNumOfSharedDigi

private

Definition at line 199 of file PndEmcAnalysis.h.

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

Int_t PndEmcAnalysis::totNumOfSubWave

private

Definition at line 196 of file PndEmcAnalysis.h.

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

Int_t PndEmcAnalysis::totNumOfWave

private

Definition at line 195 of file PndEmcAnalysis.h.

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

Double_t PndEmcAnalysis::totSharedDigiEnergy

private

Definition at line 182 of file PndEmcAnalysis.h.

Referenced by Exec(), and Init().

TTree* PndEmcAnalysis::tRec

private

Definition at line 128 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tTask

private

Definition at line 124 of file PndEmcAnalysis.h.

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

TTree* PndEmcAnalysis::tTaskPileup

private

Definition at line 126 of file PndEmcAnalysis.h.

Referenced by FinishTask(), and Init().

TTree* PndEmcAnalysis::tWave

private

Definition at line 122 of file PndEmcAnalysis.h.

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

---

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

• PndEmcAnalysis.h
• PndEmcAnalysis.cxx