PndAnalysisTask Class Reference

#include <PndAnalysisTask.h>

Inheritance diagram for PndAnalysisTask:

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

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

Protected Attributes
int	evcount
RhoNeutralParticleSelector *	neutralSel
RhoPlusParticleSelector *	plusSel
RhoMinusParticleSelector *	minusSel
RhoMassParticleSelector *	phiMSel
RhoMassParticleSelector *	pi0MSel
RhoMassParticleSelector *	dsMSel
RhoSimpleKaonSelector *	kSel
RhoSimplePionSelector *	piSel
TH1F *	phimass
TH1F *	pi0mass
TH1F *	dsmass
TH1F *	ds0mass
TH1F *	ppmass
TH1F *	nmult

Private Member Functions
virtual void	SetParContainers ()
	ClassDef (PndAnalysisTask, 1)

Private Attributes
TClonesArray *	fChargedArray
TClonesArray *	fNeutralArray

Detailed Description

Definition at line 23 of file PndAnalysisTask.h.

Member Typedef Documentation

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

Definition at line 27 of file PndAnalysisTask.h.

Constructor & Destructor Documentation

PndAnalysisTask::PndAnalysisTask

(

)

Default constructor

Definition at line 35 of file PndAnalysisTask.cxx.

                                 :
  FairTask("Panda Analysis Task")
{
}

PndAnalysisTask::~PndAnalysisTask

(

)

Destructor

Definition at line 42 of file PndAnalysisTask.cxx.

{ }

Member Function Documentation

PndAnalysisTask::ClassDef ( PndAnalysisTask  , 1    )

private

void PndAnalysisTask::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 132 of file PndAnalysisTask.cxx.

References RhoCandList::Add(), RhoCandList::Cleanup(), RhoCandList::Combine(), ds0mass, dsmass, dsMSel, evcount, fChargedArray, fNeutralArray, RhoCandList::GetLength(), RhoFactory::Instance(), kSel, RhoCandidate::M(), minusSel, neutralSel, RhoCandListIterator::Next(), nmult, phimass, phiMSel, pi0mass, pi0MSel, piSel, plusSel, ppmass, RhoCandList::RemoveClones(), RhoFactory::Reset(), and RhoCandList::Select().

{

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

  if (!(++evcount%100)) { cout <<"evt "<<evcount<<endl; }

  // **** create all the particle lists we'll need for rebuilding the decay tree
  //
  RhoCandList neutralCands,chargedCands, plusCands,minusCands;

  RhoCandList kpCands,kmCands,piCands;

  RhoCandList phiCands,pi0Cands,dsCands,ds0Cands,ppCands;


  RhoCandidate* tc;

  // **** loop over all Candidates and add them to the list allCands
  //
  chargedCands.Cleanup();
  neutralCands.Cleanup();

  for (Int_t i1=0; i1<fChargedArray->GetEntriesFast(); i1++) {
    tc = (RhoCandidate*)fChargedArray->At(i1);
    chargedCands.Add(tc);
  }

  for (Int_t i1=0; i1<fNeutralArray->GetEntriesFast(); i1++) {
    tc = (RhoCandidate*)fNeutralArray->At(i1);
    neutralCands.Add(tc);
  }

  //cout <<"c:"<<chargedCands.GetLength()<<" n:"<<neutralCands.GetLength()<<endl;


  // **** select all the basic lists
  //
  nmult->Fill(neutralCands.GetLength());

  plusCands.Select(chargedCands  ,plusSel);
  minusCands.Select(chargedCands ,minusSel);

  // **** pid selection
  //
  kpCands.Select(plusCands    ,kSel);
  kmCands.Select(minusCands   ,kSel);
  piCands.Select(chargedCands ,piSel);

  // **** now start combining all composits; inbetween plot masses
  //      before using the mass selectors
  //
  phiCands.Combine(kpCands,kmCands);

  RhoCandListIterator iterPhi(phiCands);
  while ((tc=iterPhi.Next())) { phimass->Fill(tc->M()); }
  phiCands.Select(phiMSel);

  dsCands.Combine(phiCands,piCands);

  RhoCandListIterator iterDs(dsCands);
  while ((tc=iterDs.Next())) { dsmass->Fill(tc->M()); }
  dsCands.Select(dsMSel);

  pi0Cands.Combine(neutralCands,neutralCands);

  RhoCandListIterator iterPi0(pi0Cands);
  while ((tc=iterPi0.Next())) { pi0mass->Fill(tc->M()); }
  pi0Cands.Select(pi0MSel);


  ds0Cands.Combine(dsCands,pi0Cands);

  RhoCandListIterator iterDs0(ds0Cands);
  while ((tc=iterDs0.Next())) { ds0mass->Fill(tc->M()); }

  ppCands.Combine(ds0Cands,dsCands);
  ppCands.Select(neutralSel);

  // **** since we didn't care about the D_s charge, combinations might appear in two
  //      different ways; RemoveClones removes double candidates based on the same final states
  ppCands.RemoveClones();

  RhoCandListIterator iterPp(ppCands);
  while ((tc=iterPp.Next())) { ppmass->Fill(tc->M()); }


}

void PndAnalysisTask::Finish ( )

virtual

Definition at line 222 of file PndAnalysisTask.cxx.

References ds0mass, dsmass, nmult, phimass, pi0mass, and ppmass.

{
  phimass->Write();
  pi0mass->Write();
  dsmass->Write();
  ds0mass->Write();
  ppmass->Write();
  nmult->Write();

}

InitStatus PndAnalysisTask::Init ( )

virtual

Virtual method Init

Definition at line 48 of file PndAnalysisTask.cxx.

References ds0mass, dsmass, dsMSel, evcount, fChargedArray, fNeutralArray, kSel, minusSel, neutralSel, nmult, phimass, phiMSel, pi0mass, pi0MSel, piSel, plusSel, ppmass, and RhoParticleSelectorBase::SetCriterion().

{

  //cout << " Inside the Init function****" << endl;

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

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

// Get input array
  fChargedArray = (TClonesArray*) ioman->GetObject("PndChargedCandidates");
  fNeutralArray = (TClonesArray*) ioman->GetObject("PndNeutralCandidates");

  if ( !fChargedArray && !fNeutralArray) {
    cout << "-W- PndAnalysisTask::Init: "
         << "No PndChargedCandidates && PndNeutralCandidates array!" << endl;
    return kERROR;
  }

  // Create and register output array
  cout << "-I- PndAnalysisTask: Intialization successfull" << endl;

  phimass = new TH1F("phimass","phi cands",100,0.95,1.1);
  pi0mass = new TH1F("pi0mass","pi0 cands",100,0.135-0.03,0.135+0.03);
  dsmass  = new TH1F("dsmass","Ds cands",100,1.968-0.03,1.968+0.03);
  ds0mass = new TH1F("ds0mass","Ds0 cands",100,2.317-0.05,2.317+0.05);
  ppmass = new TH1F("ppmass","pbarp cands",100,4.306-0.1,4.306+0.1);

  nmult=new TH1F("nmult","# neutrals",15,0,15);

  phimass->SetMinimum(0);
  pi0mass->SetMinimum(0);
  dsmass->SetMinimum(0);
  ds0mass->SetMinimum(0);
  ppmass->SetMinimum(0);

  // **** create and configure the selectors/filters we'd like to use later
  //
  //chargedSel = new RhoChargedParticleSelector;
  neutralSel = new RhoNeutralParticleSelector;
  plusSel    = new RhoPlusParticleSelector;
  minusSel   = new RhoMinusParticleSelector;

  // **** mass selectors for the resonances/composites
  //
  phiMSel       = new RhoMassParticleSelector("phiSelector" , 1.0195 , 0.01);
  pi0MSel       = new RhoMassParticleSelector("pi0Selector" , 0.135  , 0.005);
  dsMSel        = new RhoMassParticleSelector("dsSelector"  , 1.9685 , 0.01);

  kSel    = new RhoSimpleKaonSelector();
  kSel->SetCriterion("loose");
  piSel   = new RhoSimplePionSelector();
  piSel->SetCriterion("veryLoose");

  evcount=0;

  return kSUCCESS;

}

void PndAnalysisTask::SetParContainers ( )

privatevirtual

Geo file to use Get parameter containers

Definition at line 116 of file PndAnalysisTask.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");


}

Member Data Documentation

TH1F* PndAnalysisTask::ds0mass

protected

Definition at line 72 of file PndAnalysisTask.h.

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

TH1F* PndAnalysisTask::dsmass

protected

Definition at line 71 of file PndAnalysisTask.h.

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

RhoMassParticleSelector* PndAnalysisTask::dsMSel

protected

Definition at line 64 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

int PndAnalysisTask::evcount

protected

Definition at line 49 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndAnalysisTask::fChargedArray

private

Input array of TpcLheTrack

Definition at line 80 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndAnalysisTask::fNeutralArray

private

Definition at line 81 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

RhoSimpleKaonSelector* PndAnalysisTask::kSel

protected

Definition at line 65 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

RhoMinusParticleSelector* PndAnalysisTask::minusSel

protected

Definition at line 58 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

RhoNeutralParticleSelector* PndAnalysisTask::neutralSel

protected

Definition at line 56 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

TH1F* PndAnalysisTask::nmult

protected

Definition at line 75 of file PndAnalysisTask.h.

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

TH1F* PndAnalysisTask::phimass

protected

book all the histograms

Definition at line 69 of file PndAnalysisTask.h.

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

RhoMassParticleSelector* PndAnalysisTask::phiMSel

protected

Definition at line 62 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

TH1F* PndAnalysisTask::pi0mass

protected

Definition at line 70 of file PndAnalysisTask.h.

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

RhoMassParticleSelector* PndAnalysisTask::pi0MSel

protected

Definition at line 63 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

RhoSimplePionSelector* PndAnalysisTask::piSel

protected

Definition at line 66 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

RhoPlusParticleSelector* PndAnalysisTask::plusSel

protected

Definition at line 57 of file PndAnalysisTask.h.

Referenced by Exec(), and Init().

TH1F* PndAnalysisTask::ppmass

protected

Definition at line 73 of file PndAnalysisTask.h.

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

---

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

• PndAnalysisTask.h
• PndAnalysisTask.cxx