PndAnalysisTaskExample Class Reference

#include <PndAnalysisTaskExample.h>

Inheritance diagram for PndAnalysisTaskExample:

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

Public Member Functions
	PndAnalysisTaskExample ()
	~PndAnalysisTaskExample ()
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 (PndAnalysisTaskExample, 1)

Private Attributes
TClonesArray *	fChargedArray
TClonesArray *	fNeutralArray

Detailed Description

Definition at line 21 of file PndAnalysisTaskExample.h.

Member Typedef Documentation

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

Definition at line 25 of file PndAnalysisTaskExample.h.

Constructor & Destructor Documentation

PndAnalysisTaskExample::PndAnalysisTaskExample

(

)

Default constructor

Definition at line 36 of file PndAnalysisTaskExample.cxx.

                                               :
  FairTask("Panda Analysis Task") { 
}

PndAnalysisTaskExample::~PndAnalysisTaskExample

(

)

Destructor

Definition at line 42 of file PndAnalysisTaskExample.cxx.

{ }

Member Function Documentation

PndAnalysisTaskExample::ClassDef ( PndAnalysisTaskExample  , 1    )

private

void PndAnalysisTaskExample::Exec ( Option_t *  opt )

virtual

Virtual method Exec

Definition at line 130 of file PndAnalysisTaskExample.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 PndAnalysisTaskExample::Finish ( )

virtual

Definition at line 219 of file PndAnalysisTaskExample.cxx.

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

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

}

InitStatus PndAnalysisTaskExample::Init ( )

virtual

Virtual method Init

Definition at line 48 of file PndAnalysisTaskExample.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- PndAnalysisTaskExample::Init: "
         << "No PndChargedCandidates && PndNeutralCandidates array!" << endl;
    return kERROR;
  }

  // Create and register output array
  cout << "-I- PndAnalysisTaskExample: 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 PndAnalysisTaskExample::SetParContainers ( )

privatevirtual

Geo file to use Get parameter containers

Definition at line 115 of file PndAnalysisTaskExample.cxx.

References run.

                                              {

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

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

 
}

Member Data Documentation

TH1F* PndAnalysisTaskExample::ds0mass

protected

Definition at line 70 of file PndAnalysisTaskExample.h.

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

TH1F* PndAnalysisTaskExample::dsmass

protected

Definition at line 69 of file PndAnalysisTaskExample.h.

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

RhoMassParticleSelector* PndAnalysisTaskExample::dsMSel

protected

Definition at line 62 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

int PndAnalysisTaskExample::evcount

protected

Definition at line 47 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

TClonesArray* PndAnalysisTaskExample::fChargedArray

private

Input array of LheTrack

Definition at line 78 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

TClonesArray* PndAnalysisTaskExample::fNeutralArray

private

Definition at line 79 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

RhoSimpleKaonSelector* PndAnalysisTaskExample::kSel

protected

Definition at line 63 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

RhoMinusParticleSelector* PndAnalysisTaskExample::minusSel

protected

Definition at line 56 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

RhoNeutralParticleSelector* PndAnalysisTaskExample::neutralSel

protected

Definition at line 54 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

TH1F* PndAnalysisTaskExample::nmult

protected

Definition at line 73 of file PndAnalysisTaskExample.h.

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

TH1F* PndAnalysisTaskExample::phimass

protected

book all the histograms

Definition at line 67 of file PndAnalysisTaskExample.h.

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

RhoMassParticleSelector* PndAnalysisTaskExample::phiMSel

protected

Definition at line 60 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

TH1F* PndAnalysisTaskExample::pi0mass

protected

Definition at line 68 of file PndAnalysisTaskExample.h.

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

RhoMassParticleSelector* PndAnalysisTaskExample::pi0MSel

protected

Definition at line 61 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

RhoSimplePionSelector* PndAnalysisTaskExample::piSel

protected

Definition at line 64 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

RhoPlusParticleSelector* PndAnalysisTaskExample::plusSel

protected

Definition at line 55 of file PndAnalysisTaskExample.h.

Referenced by Exec(), and Init().

TH1F* PndAnalysisTaskExample::ppmass

protected

Definition at line 71 of file PndAnalysisTaskExample.h.

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

---

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

• PndAnalysisTaskExample.h
• PndAnalysisTaskExample.cxx