Particle identification and Particle identification and Hadron/Jet - - PowerPoint PPT Presentation

▶
particle identification and particle identification and
SMART_READER_LITE
LIVE PREVIEW

Particle identification and Particle identification and Hadron/Jet - - PowerPoint PPT Presentation

Dec 06, 2022 6 likes •171 views

Particle identification and Particle identification and Hadron/Jet correlations analysis f framework k Description, Status, HOW TO use Description, Status, HOW TO use Description, Status, HOW TO use Description, Status, HOW TO use and TO

slide-1
SLIDE 1

Particle identification and Particle identification and Hadron/Jet correlations analysis f k framework

Description, Status, HOW TO use Description, Status, HOW TO use Description, Status, HOW TO use Description, Status, HOW TO use and TO DO list and TO DO list

03/06/08 Gustavo Conesa Balbastre @ PWG4

1/14

slide-2
SLIDE 2

Recent changes in ESDs and AODs Recent changes in ESDs and AODs

  • AliESDCaloCluster and AliAODCaloCluster:
  • Removed fM11 data member, added fTOF data member. Also added/removed related Setters and

Getters Getters.

  • Modifications needed in AliEMCALReconstructor, AliPHOSReconstructor and HLT
  • Backward compatible (as far as I know)
  • AliAODCluster:
  • Particle type list changed now is the same as in AliPID plus kCharged and kNeutral
  • Particle type list changed, now is the same as in AliPID plus kCharged and kNeutral
  • Added Getter of the labels array
  • AliAODEvent:
  • Added methods to get the list of PHOS or EMCAL AODs like for the ESDs.
  • Actual version of AliAnalysisTaskESDfilter is complete in the part of filling
  • Actual version of AliAnalysisTaskESDfilter is complete in the part of filling

AliAODCaloClusters

  • Until 2 months ago the matched tracks, the labels, PID and shower shape information were not filled
  • Last week 2 bugs corrected.
  • REMINDER!!!!
  • EMCAL Pseudoclusters do not exist anymore.
  • Any method in ESDCaloCluster with the name digits in it is doing nothing, it is just there for backward

compatibility reasons.

  • In AliESD(AOD)CaloClusters we keep for each cluster the list of AbsId of the digits that created the

cluster and the list of the fraction of the digits energy used to create the cluster (used only by PHOS right now).

  • In AliESD(AOD)CaloCells we keep the full list of digits amplitude, absId and time (time is not kept in

AODs)

03/06/08 Gustavo Conesa Balbastre @ PWG4

2/14

slide-3
SLIDE 3

What has changed in the PWG4 frame What has changed in the PWG4 frame

Frame completely redesigned since last Alice

  • ffline meeting, now:

g,

AODs are produced and read. More modular: initially developed for the case of

y p direct gamma and correlations physics. Now more kind of analysis can be added easily.

Library name: PWG4Gamma to PWG4PartCorr

03/06/08 Gustavo Conesa Balbastre @ PWG4

3/14

slide-4
SLIDE 4

Particle identification and correlations analysis frame description: Schema frame description: Schema

AliAnalysisTaskSE

AliAnalisisTaskParticleCorrelation

Data

y

Output:

AliAnaMaker AliAODCaloTrackReader ESD, AOD, MC

Data aodEMCal,

Output: AliAODParticleCorrelation and/or AliAODCaloClusters and/or Hi t

aodPHOS, aodCTS

Histograms

aodEMCal, aodPHOS, aodCTS aodParticleCorr AliAnaNNN

derives from AliAnaBaseClass

AliAnaYYY

derives from AliAnaBaseClass

AliAnaXXX

derives from AliAnaBaseClass

…

03/06/08 Gustavo Conesa Balbastre @ PWG4

4/14

slide-5
SLIDE 5

Analysis manager y g

Steerin class t c nnect analysis

Bool_t AliAnaMaker::ProcessEvent(Int_t iEntry){

  • Steering class to connect analysis

and data readers : AliAnaMaker AliAnaMaker

  • Connects the data (ESD/AOD/MC)

h h l l h

… //Each event needs an empty branch fAODBranch->Clear(); //Tell the reader to fill the data in the detector lists ReadInputEvent();

with the analysis algorithm

  • Executes all analysis that user has

provided in the appropriate order.

//Loop on analysis algorithms Int_t nana = fAnalysisContainer->GetEntries() ; for(Int_t iana = 0; iana < nana; iana++){ AliAnaBaseClass * ana = ((AliAnaBaseClass *) fAnalysisContainer->At(iana)) ;

  • First analysis generating AODs (if

needed)

  • Second analysis generating

histograms from AODs (if needed)

(( ) y ( )) ; //Set reader and aod branch for each analysis ana->SetReader(fReader); ana->SetAODBranch(fAODBranch); //Make analysis, create aods in aod branch or AODCaloClusters if(fMakeAOD) ana >MakeAnalysisFillAOD() ;

histograms from AODs (if needed)

  • Several analysis can be executed one

after the other, taking as input the AOD output of the precedent

if(fMakeAOD) ana->MakeAnalysisFillAOD() ; //Make further analysis with aod branch and fill histograms if(fMakeHisto) ana->MakeAnalysisFillHistograms() ; } fReader->ResetLists();

analysis.

  • Initializes all parameters and

histograms connects them with

03/06/08 Gustavo Conesa Balbastre @ PWG4

5/14

  • utput of general analysis frame.
slide-6
SLIDE 6

Particle identification and correlation analysis frame description: Data Readers description: Data Readers

  • Classes to read the ESD, AOD or MC derive from: AliCaloTrackReader

AliCaloTrackReader

  • Output 3 TClonesArray: EMCAL AODCaloClusters, PHOS AODCaloCluster,

CTS AODTrack.

  • Arrays filled if user wants (by default all filled): SwitchOnCTS();

SwitchOnEMCAL(); SwitchOnPHOS()

  • Recover the lists in analysis: GetAODCTS(); GetAODEMCAL(); GetAODPHOS();
  • Also provides the corresponding CaloCells (not in MC case). If reader is ESD

p p g it returns AliESDCaloCells, if AOD it returns AliAODCaloCells:

  • SwitchOnEMCALCells(); SwitchOnEMCALCells()
  • GetEMCALCells(); GetPHOSCells()

(); ()

  • You can get Generator information:

virtual AliStack* GetStack() const ;

i l AliH d * G H d ()

virtual AliHeader* GetHeader() const ; virtual AliGenEventHeader* GetGenEventHeader() const ; More if needed

03/06/08 Gustavo Conesa Balbastre @ PWG4

6/14

slide-7
SLIDE 7

Particle identification and correlation analysis frame description: Data Readers

  • AliCaloTrackESDReader: Done and tested locally and on grid.

description: Data Readers

  • AliCaloTrackAODReader: Done, test in progress.
  • AliCaloTrackMCReader: Done and tested locally.

Output 3 TClonesArray: EMCAL PHOS and CTS TParticles or standard Output 3 TClonesArray: EMCAL, PHOS and CTS TParticles or standard

AODs whatever the user wants.

  • SetClonesArrayType(Int_t type)
  • 1 (default) is AOD, 0 is TParticle

S l i l i h i d

Select particles with a given status code:

  • SwitchOnStatusSelection(); KeepParticleWithStatus(Int_t status)
  • Default keeps status 1

In CTS array keep only charged particles in CTS acceptance In CTS array keep only charged particles in CTS acceptance In Calo arrays keep only neutral particles in Calo acceptances

  • Skip some neutral particles (default skip neutrinos):

SkipNeutralParticles(Int_t pdg)

  • Add some charged particles (electrons): KeepChargedParticles(Int_t pdg)

In PYTHIA by default pi0 has status 1 and is not decayed. Option to

decay it by “hand”:

SwitchOnPi0Decay()

03/06/08 Gustavo Conesa Balbastre @ PWG4

7/14

  • SwitchOnPi0Decay()
slide-8
SLIDE 8

Analysis algorithm base class

  • Analysis base class:

AliAnaBaseClass AliAnaBaseClass

//To fill AODs virtual void AddAODCaloCluster(AliAODCaloCluster calo) ; virtual void AddAODParticleCorrelation(AliAODParticleCorrelation pc) ;

  • All analysis classes must

derive from this one.

  • External modules that can be

called by the analysis classes

virtual void ConnectAODCaloClusters(); //Init Histograms virtual TList * GetCreateOutputObjects(){ return (new TList) ;} //MAIN ANALYSIS METHODS virtual void MakeAnalysisFillAOD() {;}//Creates AODPhotons or AODCaloClusters

called by the analysis classes

  • AliCaloPID
  • Given the PID weights

returns probable pdg id

  • Use as base class for PID

virtual void MakeAnalysisFillAOD() {;}//Creates AODPhotons or AODCaloClusters that will be kept in file virtual void MakeAnalysisFillHistograms() {;} //Final analysis with AODPhotons and production of histograms //Reader virtual AliCaloTrackReader * GetReader() const {return fReader ; } i t l id S tR d (AliG R d * d ) { fR d d }

  • Use as base class for PID

recalculation?

  • AliFidutialCut:
  • Acceptance selection
  • Allows selection in n regions

virtual void SetReader(AliGammaReader * reader) { fReader = reader ; } //AOD lists, outputs //AOD lists, outputs virtual TClonesArray* GetAODBranch() const {return fAODBranch ;} virtual TClonesArray* GetAODCaloClusters() const {return fAODCaloClusters ;} //AOD lists, data for the analysis virtual TClonesArray* GetAODCTS() const ;

Allows selection in n regions in the detectors

  • AliNeutralMesonSelection
  • Given a pair of photons

select event by event as π0 b i i t

y () virtual TClonesArray* GetAODEMCAL() const ; virtual TClonesArray* GetAODPHOS() const ; virtual AliStack * GetMCStack() const ; virtual void SetAODBranch(TClonesArray * tca) { fAODBranch = tca ; }

  • r η by invariant mass or
  • pening angle
  • AliIsolationCut
  • Selects particle if isolated.

virtual AliCaloPID * GetCaloPID() const {return fCaloPID ;} virtual void SetCaloPID(AliCaloPID * pid) { fCaloPID = pid ;} virtual AliIsolationCut * GetIsolationCut() const {return fIC ;} virtual void SetIsolationCut(AliIsolationCut * fc) { fIC = fc ;} virtual AliFidutialCut * GetFidutialCut() const {return fFidCut ;} 03/06/08 Gustavo Conesa Balbastre @ PWG4

8/14

  • more if needed.

virtual void SetFidutialCut(AliFidutialCut * fc) { fFidCut = fc ;}

slide-9
SLIDE 9

Output AODs: AliAODParticleCorrelation Output AODs: AliAODParticleCorrelation

  • Output can be « standard » AODs, AODCaloClusters, or « non

standard » AODs AliAODParticleCorrelation standard » AODs, AliAODParticleCorrelation.

  • AliAODParticleCorrelation is basically a list of TLorentzVectors

h i h h f ll i i f i each one with the following extra information:

  • Several TRefArrays:

Array of cluster/tracks in the isolation cone Array of cluster/track used for hadron correlation or used for jet

reconstruction (if JETAN not used)

Array of cluster/track used for background estimation

PID of the TLVector (photon pi0 hadron )

  • PID of the TLVector (photon, pi0, hadron …).
  • Monte Carlo label.
  • Another label to be set by user

i ith th h l f MC I l b l th h t f d

  • in my case, with the help of MC, I label the photons from decay,

fragmentation and conversion, for example.

  • TLorentzVectors with found Leading particle, Jet and Background

kinematics

03/06/08 Gustavo Conesa Balbastre @ PWG4

9/14

kinematics.

  • Reference to jet found with JETAN
slide-10
SLIDE 10

Analysis class example

//_________________________________________________ TList * AliAnaExample::GetCreateOutputObjects() { // Create histograms to be saved in output file and // store them in fOutputContainer

y p

AliDebug(1,"Init parton histograms"); TList * outputContainer = new TList() ;

  • utputContainer->SetName("GammaExampleHistos") ;

//__________________________________________________________________ void AliAnaExample::MakeAnalysisFillAOD() { fhPt = new TH1F ("hPt","p_T distribution", 100,0,100); fhPt->SetXTitle("p_{T} (GeV/c)");

  • utputContainer->Add(fhPt);

fhPhi = new TH1F ("hPhi","#phi distribution", 100,0,TMath::TwoPi()); { //Do analysis and fill aods ConnectAODCaloClusters(); //Do Only when filling AODCaloClusters //Get vertex for photon momentum calculation Double_t v[3] ; //vertex ; 100,0,TMath::TwoPi()); fhPhi->SetXTitle("#phi (rad)");

  • utputContainer->Add(fhPhi);

… return outputContainer; } GetReader()->GetVertex(v); //Fill AODCaloClusters and AODParticleCorrelation with aods PHOS for(Int_t i = 0; i < GetAODPHOS()->GetEntries(); i++){ AliAODCaloCluster * calo = dynamic_cast<AliAODCaloCluster*> (GetAODPHOS()->At(i)); (GetAODPHOS()->At(i)); … do whatever selection with the clusters //Fill AODCaloClusters AddAODCaloCluster(*calo); //Fill AODPhoton after some selection //________________________________________________________ void AliAnaExample::MakeAnalysisFillHistograms() { //Do analysis and fill histograms TLorentzVector mom ; calo->GetMomentum(mom,v); Double_t pid[10]; calo->GetPID(pid); Int t pdg = GetCaloPID()->GetPdg("PHOS",pid,mom.E()); //Loop on stored AOD photons Int_t naod = GetAODBranch()->GetEntriesFast(); if(GetDebug() > 0) printf("histo aod branch entries %d\n", naod); for(Int_t iaod = 0; iaod < naod ; iaod++){ AliAODPhoton* ph = _ p g () g( ,p , ()); Bool_t in = GetFidutialCut()->IsInFidutialCut(mom,"PHOS") ; if(GetDebug() > 0) printf("pT %3.2f, in %d, pdg %d\n",mom.Pt(), in,pdg); if(mom.Pt() > fPtMin && in && pdg ==AliCaloPID::kPhoton) { AliAODParticleCorrelation ph = AliAODParticleCorrelation(mom); ph SetLabel(calo >GetLabel(0)); dynamic_cast<AliAODPhoton*> (GetAODBranch()- >At(iaod)); fhPt->Fill(ph->Pt()); fhPhi->Fill(ph->Phi()); 03/06/08 Gustavo Conesa Balbastre @ PWG4

10/14

ph.SetLabel(calo->GetLabel(0)); ph.SetPdg(pdg); AddAODParticleCorrelation(ph); } … } }

slide-11
SLIDE 11

Configuration file example Configuration file example

AliAnaMaker* ConfigAnalysis(){ //Detector Fidutial Cuts AliFidutialCut * fidCut = new AliFidutialCut(); fidCut >DoCTSFidutialCut(kFALSE) ; fidCut->DoCTSFidutialCut(kFALSE) ; fidCut->DoEMCALFidutialCut(kTRUE) ; fidCut->DoPHOSFidutialCut(kTRUE) ; fidCut->SetSimpleCTSFidutialCut(0.9,0.,360.); fidCut->SetSimpleEMCALFidutialCut(0.7,80.,190.); fidCut->SetSimplePHOSFidutialCut(0.13,220.,320.); // Reader AliCaloTrackReader *reader = new AliCaloTrackESDReader(); //Min particle pT reader->SetEMCALPtMin(2.); reader->SetPHOSPtMin(2.); e de Se OS ( .); reader->SetCTSPtMin(.2); reader->SetFidutialCut(fidCut); reader->SwitchOnEMCAL(); reader->SwitchOnPHOS(); reader->SwitchOffCTS(); // Analysis algorithm AliCaloPID * pid = new AliCaloPID(); AliCaloPID * pid = new AliCaloPID(); // use selection with simple weights pid->SetPHOSPhotonWeight(0.7); pid->SetPHOSPi0Weight(0.7); pid->SetEMCALPhotonWeight(0.7); pid->SetEMCALPi0Weight(0.7); AliAnaExample *ana = new AliAnaExample(); ana->SetCaloPID(pid); ana->SetFidutialCut(fidCut2); //AliAnaExample *ana2 = new AliAnaExample(); // Set analysis algorithm and reader maker = new AliAnaMaker(); maker->SetReader(reader);//pointer to reader maker->AddAnalysis(ana,0); 03/06/08 Gustavo Conesa Balbastre @ PWG4

11/14

maker >AddAnalysis(ana,0); //maker->AddAnalysis(ana2,1); maker->SetAODBranchName(« MyAOD"); maker->SwitchOnHistogramsMaker() ; maker->SwitchOnAODsMaker() ; return maker ;}

slide-12
SLIDE 12

Direct Gamma analysis code and hadron/jet correlation correlation

  • 2 different analysis modules to be executed one after the other
  • r even independently if the AODParticleCorrelations with all the
  • r even independently if the AODParticleCorrelations with all the

necessary information are already created

  • AliAnaGammaDirect
  • Looks for high pT trigger particle: photons identified

Looks for high pT trigger particle photons identified

  • Check if they are direct: Isolation Cut
  • Fill AODParticleCorrelation and/or histograms
  • Test in progress

p g

  • Correlations: With the trigger particle (direct photon or

anything) correlate

  • AliAnaParticlePartonCorrelation: Delta pT, eta and phi plots. Done,

t t t d not tested.

  • AliAnaParticleHadronCorrelation: Azimuthal correlation plots and

momentum imbalance correlation function. Done, not tested.

  • AliAnaParticleLeadJetCone: Non standard jet finder with the
  • AliAnaParticleLeadJetCone: Non standard jet finder with the

trigger particle as seed for the back-jet. Done, not tested.

  • AliAnaParticleJetFinder: Connects trigger particle and back-jet

found with JETAN (to be done)

03/06/08 Gustavo Conesa Balbastre @ PWG4

12/14

  • Fill AODParticleCorrelations and/or histograms
slide-13
SLIDE 13

HOW TO W

  • Recipe, ingredients needed:

“ f l l f E

  • “par” file containing analysis frame: ESD.par, AOD.par,

STEERBase.par, ANALYSIS.par, ANALYSISalice.par, PWG4PartCorr.par

$

Do in $ALICE_ROOT, “make PWG4PartCorr.par”, or let the analysis

macro do it for you.

  • Configuration file: ConfigXXX.C ()

Examples ConfigAnalysisGammaDirect.C and ConfigAnalysisExample.C

inside $ALICE_ROOT/PWG4/macros

  • Analysis execution file: ana.C

Example inside $ALICE_ROOT/PWG4/macros

  • Exec “root ana.C”
  • Read carefully ana.C, performs analysis locally, on grid or CAF (last
  • Read carefully ana.C, performs analysis locally, on grid or CAF (last

not tested) and different options can be set.

03/06/08 Gustavo Conesa Balbastre @ PWG4

13/14

slide-14
SLIDE 14

Summary and TO DO Summary and TO DO

  • New particle and correlations analysis frame
  • New particle and correlations analysis frame,

1.

Originally designed for photon tagging hadrons/jets

2.

More modularity, new analysis algorithms can be added easily P d AOD d/ hi t

3.

Produces AODs and/or histograms

4.

Can be used by everyone, you are invited to help me improve it.

  • To Do

1.

Needs full test locally and with GRID and as a wagon of the analysis frame.

2.

Add new analysis cases, specially analysis relevant for the first h i physics.

3.

Implement event mixing.

4.

Write tutorial

5.

…

03/06/08 Gustavo Conesa Balbastre @ PWG4

14/14

slide-15
SLIDE 15

Back-up Back up

03/06/08 Gustavo Conesa Balbastre @ PWG4

15/14

slide-16
SLIDE 16

General Analysis Frame short description General Analysis Frame, short description

  • In order to do an analysis you need Root or AliRoot and the libraries:
  • STEERbase

ESD

  • ESD
  • AOD
  • ANALYSIS
  • ANALYSISalice, contains:

Th m s f i ndl AliAn l i T kSE

  • The more user friendly AliAnalysisTaskSE.
  • It is recommended that all the user analysis code derives from this class.
  • AliAnalysisTaskESDfilter: Main wagon of the analysis train, in charge of creating the standard

AODs: AODtracks and AODCaloClusters.

  • Your analysis library

Th h l

  • The way to execute the analysis:

1)

Load the libraries: Create a par file with them (cd $ALICE_ROOT; make ESD.par;) and compile the code.

2)

Init the manager (AliAnalysisManager) and data handlers (MCHandler, AODHandler, ESDHandler) your analysis task and output containers ESDHandler), your analysis task and output containers

  • Output can be AODs or histograms

3)

Create the chain of data, ESDs or AODs

4)

Feed the manager with the chain of data and execute the analysis.

  • Example of analysis task deriving from AliAnalysisTaskSE in
  • Example of analysis task deriving from AliAnalysisTaskSE in

$ALICE_ROOT/PWG4/AliAnalysisTaskPHOSExample.cxx, .h

  • Example of macro to execute analysis in

$ALICE ROOT/PWG4/macros/anaExampleTask C

03/06/08 Gustavo Conesa Balbastre @ PWG4

16/14

$ALICE_ROOT/PWG4/macros/anaExampleTask.C