PndJpsi2PiAna Class Reference

#include <PndJpsi2PiAna.h>

Inheritance diagram for PndJpsi2PiAna:

Public Member Functions
	PndJpsi2PiAna ()
	~PndJpsi2PiAna ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *opt)
virtual void	Finish ()

Protected Member Functions
void	InitArrays ()
int	GetEvent ()
bool	FillList (RhoCandList &l, std::string listkey="All")

Protected Attributes
RhoMassParticleSelector *	jpsiMSel
TH1F *	jpsimass
TH1F *	jpsimassf
TH1F *	ppmass
TH1F *	ppmassf
TH1F *	mcmass

Private Member Functions
	ClassDef (PndJpsi2PiAna, 1)

Private Attributes
UInt_t	fEvtCount

Detailed Description

Definition at line 13 of file PndJpsi2PiAna.h.

Constructor & Destructor Documentation

PndJpsi2PiAna::PndJpsi2PiAna

(

)

Default constructor

Definition at line 34 of file PndJpsi2PiAna.cxx.

                             : PndAnaTask()
{

}

PndJpsi2PiAna::~PndJpsi2PiAna

(

)

Destructor

Definition at line 41 of file PndJpsi2PiAna.cxx.

{
}

Member Function Documentation

PndJpsi2PiAna::ClassDef ( PndJpsi2PiAna  , 1    )

private

void PndJpsi2PiAna::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Reimplemented from PndAnaTask.

Definition at line 71 of file PndJpsi2PiAna.cxx.

References RhoCandList::Combine(), fEvtCount, PndAnaTask::FillList(), Rho4CFitter::FitConserveMasses(), PndAnaTask::GetEvent(), RhoCandList::GetLength(), i, jpsimass, jpsimassf, RhoCandidate::M(), mcmass, RhoCandidate::P4(), ppmass, and ppmassf.

{
        if (0==(++fEvtCount)%100) cout <<"evt "<<fEvtCount<<endl;

    int i,j;

        GetEvent();

    RhoCandList ep, em, pip, pim, jpsi, pp, mc;

        FillList(ep,"ElectronVeryLoosePlus");
        FillList(em,"ElectronVeryLooseMinus");
        FillList(pip,"PionVeryLoosePlus");
        FillList(pim,"PionVeryLooseMinus");

        FillList(mc,"McTruth");

        TLorentzVector ini(0, 0, 6.231711, 7.240222);

/*      cout<<" e+:" <<ep.GetLength();
        cout<<" e-:" <<em.GetLength();
        cout<<" pi+:" <<pip.GetLength();
        cout<<" pi-:" <<pim.GetLength();
        cout <<endl;*/

        jpsi.Combine(ep,em);

/*      cout <<jpsi.GetLength();*/

        for (j=0;j<jpsi.GetLength();++j) jpsimass->Fill(jpsi[j]->M());
        //jpsi.Select(jpsiMSel);

        for (j=0;j<mc.GetLength();++j) mcmass->Fill(mc[j]->M());

        pp.Combine(jpsi,pip,pim);

        for (j=0;j<pp.GetLength();++j)
        {
                        ppmass->Fill(pp[j]->M());

                        //do the 4C fit on the pbar p System
                        Rho4CFitter fitter(pp[j],ini);

                        fitter.FitConserveMasses();

            RhoCandidate *ppfit=pp[j]->GetFit();
                        ppmassf->Fill(ppfit->M());

            RhoCandidate *epfit=(pp[j]->Daughter(0)->Daughter(0))->GetFit();
            RhoCandidate *emfit=(pp[j]->Daughter(0)->Daughter(1))->GetFit();

                        TLorentzVector sum=epfit->P4()+emfit->P4();

                        jpsimassf->Fill(sum.Mag());

        }


}

bool PndAnaTask::FillList ( RhoCandList &  l, std::string  listkey = "All"  )

protectedinherited

Definition at line 93 of file PndAnaTask.cxx.

References RhoCandList::Add(), RhoCandList::Cleanup(), PndAnaTask::fChargedCands, PndAnaTask::fChargedProbability, PndPidListMaker::FillList(), PndAnaTask::fLAllCands, PndAnaTask::fLChargedCands, PndAnaTask::fLNeutralCands, PndAnaTask::fMcCands, PndAnaTask::fNeutralCands, PndAnaTask::fNeutralProbability, PndAnaTask::fPidListMaker, PndPidProbability::GetElectronPidProb(), GetEntriesFast(), PndPidProbability::GetKaonPidProb(), RhoCandList::GetLength(), PndPidProbability::GetMuonPidProb(), PndPidProbability::GetPionPidProb(), PndPidProbability::GetProtonPidProb(), PndPidListMaker::SetBaseList(), RhoCandidate::SetPidInfo(), and uid().

Referenced by Exec().

{
  // Reads the specified List for the current event

  l.Cleanup();
  UInt_t uid=1;

  // when the first list is requested read in the event
  /*
  if (!fEventRead)
  {
    fRootManager->ReadEvent(fEvtCount-1);
    fEventRead=true;
  }
  */

  if (listkey=="McTruth") {
    if (fMcCands) {
      //cout <<"PndAnaTask::FillList: "<<fMcCands->GetEntriesFast()<<" mc cands"<<endl;
      for (int i1=0; i1<fMcCands->GetEntriesFast(); i1++) {
        RhoCandidate* tc = (RhoCandidate*)fMcCands->At(i1);
        l.Add(tc);
      }
      return true;
    } else { return false; }
  }

  // fill all, neutral and charged from the PndPidCandidate Array
//  if (fMicroCands && allCands.GetLength()==0)
//  {
//    for (int i1=0; i1<fMicroCands->GetEntriesFast(); i1++)
//    {
//      VAbsMicroCandidate *mic = (VAbsMicroCandidate *)fMicroCands->At(i1);
//      RhoCandidate tc(*mic,i1+1);
//
//      allCands.Add(tc);
//
//      if (fabs(tc.Charge())>0.01)
//        chargedCands.Add(tc);
//      else
//        neutralCands.Add(tc);
//    }
//  }
  //else
  if (fLAllCands.GetLength() == 0) { // do only when we didn't read something yet.
    // removed now compatibility to RhoCandidate readin ... instead read PndPidCandidates
    if (fNeutralCands && fLNeutralCands.GetLength()==0)
      for (int i1=0; i1<fNeutralCands->GetEntriesFast(); i1++) {
        PndPidCandidate* mic = (PndPidCandidate*)fNeutralCands->At(i1);
        RhoCandidate tc(*mic,uid++);

        // TODO: Do we want to set something here? It is neutrals anyway.
        if(0!=fNeutralProbability && i1<fNeutralProbability->GetEntriesFast()) {
          PndPidProbability* neuProb = (PndPidProbability*)fNeutralProbability->At(i1);
          // numbering see PndPidListMaker
          tc.SetPidInfo(0,neuProb->GetElectronPidProb());
          tc.SetPidInfo(1,neuProb->GetMuonPidProb());
          tc.SetPidInfo(2,neuProb->GetPionPidProb());
          tc.SetPidInfo(3,neuProb->GetKaonPidProb());
          tc.SetPidInfo(4,neuProb->GetProtonPidProb());
        }
        fLNeutralCands.Add(&tc);
        fLAllCands.Add(&tc);
      }

    if (fChargedCands && fLChargedCands.GetLength()==0) {
      for (int i1=0; i1<fChargedCands->GetEntriesFast(); i1++) {
        PndPidCandidate* mic = (PndPidCandidate*)fChargedCands->At(i1);
        RhoCandidate tc(*mic,uid++);

        if(0!=fChargedProbability && i1<fChargedProbability->GetEntriesFast()) {
          PndPidProbability* chProb = (PndPidProbability*)fChargedProbability->At(i1);
          // numbering see PndPidListMaker
          tc.SetPidInfo(0,chProb->GetElectronPidProb());
          tc.SetPidInfo(1,chProb->GetMuonPidProb());
          tc.SetPidInfo(2,chProb->GetPionPidProb());
          tc.SetPidInfo(3,chProb->GetKaonPidProb());
          tc.SetPidInfo(4,chProb->GetProtonPidProb());
        }
        fLChargedCands.Add(&tc);
        fLAllCands.Add(&tc);
      }
    }
  }


  // set the base list for the PID list maker
  fPidListMaker->SetBaseList(fLChargedCands);

  if (listkey=="All" ) {
    l=fLAllCands;
    return true;
  }

  if (listkey=="Neutral") {
    l=fLNeutralCands;
    return true;
  }

  if (listkey=="Charged") {
    l=fLChargedCands;
    return true;
  }

  return fPidListMaker->FillList(l,listkey);
}

void PndJpsi2PiAna::Finish ( )

virtual

Reimplemented from PndAnaTask.

Definition at line 133 of file PndJpsi2PiAna.cxx.

References jpsimass, jpsimassf, mcmass, ppmass, and ppmassf.

{
  jpsimass->Write();
  jpsimassf->Write();
  ppmass->Write();
  ppmassf->Write();
  mcmass->Write();
}

int PndAnaTask::GetEvent ( )

protectedinherited

Definition at line 81 of file PndAnaTask.cxx.

References RhoCandList::Cleanup(), PndAnaTask::fLAllCands, PndAnaTask::fLChargedCands, PndAnaTask::fLMcCands, PndAnaTask::fLNeutralCands, RhoFactory::Instance(), and RhoFactory::Reset().

Referenced by Exec().

{
  RhoFactory::Instance()->Reset();

  fLAllCands.Cleanup();
  fLChargedCands.Cleanup();
  fLNeutralCands.Cleanup();
  fLMcCands.Cleanup();

  return 0;
}

InitStatus PndJpsi2PiAna::Init ( )

virtual

Virtual method Init

Reimplemented from PndAnaTask.

Definition at line 47 of file PndJpsi2PiAna.cxx.

References fEvtCount, PndAnaTask::InitArrays(), jpsimass, jpsimassf, jpsiMSel, mcmass, ppmass, and ppmassf.

{
        std::cout << "-I- PndJpsi2PiAna: Initialization" << std::endl;

        InitArrays();


        jpsimass=new TH1F("mjpsi","jpsi",200,0,3.5);
        jpsimassf=new TH1F("mjpsif","jpsi fit",200,0,3.5);

        ppmass=new TH1F("mpp","mpp",200,0,4.5);
        ppmassf=new TH1F("mppf","mpp",200,0,4.5);

        mcmass=new TH1F("mmc","mmc",200,0,4.0);

        jpsiMSel=new RhoMassParticleSelector("jpsi",3.1,0.2);

    fEvtCount=0;

    return kSUCCESS;
}

void PndAnaTask::InitArrays ( )

protectedinherited

Definition at line 59 of file PndAnaTask.cxx.

References PndAnaTask::fChargedCands, PndAnaTask::fChargedProbability, PndAnaTask::fMcCands, PndAnaTask::fNeutralCands, PndAnaTask::fNeutralProbability, and PndAnaTask::fPidListMaker.

Referenced by Init().

{
  std::cout << "-I- PndAnaTask: Initialization" << std::endl;

  FairRootManager* fRootManager = FairRootManager::Instance();

  if ( 0 == fRootManager ) {
    std::cout << "-E- PndAnaTask: No FairRootManager found!" << std::endl;
    return;
  }

  //read
  fChargedCands = (TClonesArray*) fRootManager->GetObject("PidChargedCand");
  fNeutralCands = (TClonesArray*) fRootManager->GetObject("PidNeutralCand");
  fChargedProbability = (TClonesArray*) fRootManager->GetObject("PidChargedProbability");
  fNeutralProbability = (TClonesArray*) fRootManager->GetObject("PidNeutralProbability");
  fMcCands = (TClonesArray*) fRootManager->GetObject("PndMcTracks");

  fPidListMaker = new PndPidListMaker();

}

Member Data Documentation

UInt_t PndJpsi2PiAna::fEvtCount

private

Definition at line 46 of file PndJpsi2PiAna.h.

Referenced by Exec(), and Init().

TH1F* PndJpsi2PiAna::jpsimass

protected

book all the histograms

Definition at line 39 of file PndJpsi2PiAna.h.

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

TH1F* PndJpsi2PiAna::jpsimassf

protected

Definition at line 40 of file PndJpsi2PiAna.h.

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

RhoMassParticleSelector* PndJpsi2PiAna::jpsiMSel

protected

Definition at line 36 of file PndJpsi2PiAna.h.

Referenced by Init().

TH1F* PndJpsi2PiAna::mcmass

protected

Definition at line 43 of file PndJpsi2PiAna.h.

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

TH1F* PndJpsi2PiAna::ppmass

protected

Definition at line 41 of file PndJpsi2PiAna.h.

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

TH1F* PndJpsi2PiAna::ppmassf

protected

Definition at line 42 of file PndJpsi2PiAna.h.

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

---

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

• PndJpsi2PiAna.h
• PndJpsi2PiAna.cxx