AliPhysics  fc9925d (fc9925d)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Mid-rapidity tracklet code for dN/deta
Collaboration diagram for Mid-rapidity tracklet code for dN/deta:

Files

file  AliAODSimpleHeader.C
 A simplified AOD header.
 
file  AliAODTracklet.C
 Tracklet AOD object.
 
file  AliSimpleHeaderTask.C
 A task to make a simplified AOD header.
 
file  AliTrackletAODdNdeta.C
 AOD tasks to do final dN/deta in midrapidity.
 
file  AliTrackletAODTask.C
 Tasks to make tracklet AOD output.
 
file  AliTrackletAODUtils.C
 Utilities for midrapidity analysis.
 
file  AliTrackletdNdeta.C
 Post processing.
 
file  AliTrackletdNdeta2.C
 To post processing 2nd version.
 
file  AliTrackletWeights.C
 Encode simulation weights for 2nd pass.
 
file  BackOfTheEnvelope.C
 Do back of the envelope calculation of effect of reweighing.
 
file  CompareResults.C
 Compare different results.
 
file  DetailsCompare.C
 
file  BackOfTheEnvelope.C
 Do back of the envelope calculation of effect of reweighing.
 
file  ExtractGSE.C
 Extract GraphSysErr from results.
 
file  ExtractGSE2.C
 Extract GraphSysErr from results.
 
file  ExtractMCWeights.C
 
file  MakeDeltas.C
 
file  MakeDeltaWeights.C
 Make weight based on delta, eta, ipz.
 
file  MakeEposWeight.C
 
file  MakeK0S.C
 
file  MakeRealistic.C
 
file  MakeRuben.C
 
file  MakeStrange.C
 
file  MakeWeight.C
 
file  PureMCWeights.C
 
file  ReadWeights.C
 
file  TrackletAODdNdeta.C
 A tracklet dNdeta train.
 
file  TrackletAODTrain.C
 Tracklet AOD train.
 
file  TryBuild.C
 

Classes

class  AliAODSimpleHeader
 
class  AliAODTracklet
 
class  AliAODMCTracklet
 
class  AliSimpleHeaderTask
 
class  AliTrackletAODdNdeta
 
class  AliTrackletAODMCdNdeta
 
class  AliTrackletAODWeightedMCdNdeta
 
class  AliTrackletAODTask
 
class  AliTrackletAODMCTask
 
class  AliTrackletAODUtils
 
struct  AliTrackletdNdeta
 
struct  AliTrackletdNdeta2
 
class  AliTrackletPtPidStrWeights
 
struct  Calculation
 
struct  Calculation::Particle
 
struct  DetailsComparer
 
struct  FixPaths
 
struct  MakeDeltaWeights
 
struct  PureMCWeights
 
struct  TrackletAODdNdeta
 
struct  TrackletAODTrain
 

Functions

void BackOfTheEnvelope (const char *file, Double_t c1=0, Double_t c2=0, Bool_t mid=true)
 
void DetailsCompare ()
 
void Expectations (Bool_t mid=true, Double_t c1=0, Double_t c2=5)
 
void ExtractMCWeights (const char *nFile="dpmjet.root", const char *dFile="hijing.root", const char *oFile=0)
 
void MakeDeltas (const char *realFileName, const char *simFileName)
 
void TryBuild ()
 

Variables

const char * fmt = "%-10s | %5.3f | %5.2f | %5.2f | %5.2f %5.2f |"
 
const char * parFmt = "%-10s | %5.3f |"
 

Detailed Description

Overview

The code in this module constitutes tools for analysing SPD tracklet data for the charged-particle pseudorapidity density. It is based on Ruben's original code (see pwglf_fwd_spd_tracklet_1), but differs in some important aspects.

To use

AOD production

To produce the AODs with the tracklet information in, one needs to run a train with a task of the class AliTrackletAODTask (or AliTrackletAODMCTask for simulated data) and a task of the class AliSimpleHeaderTask in it. This is most easily done using the TrainSetup (Using the TrainSetup facility) derived class TrackletAODTrain.

For example for real data from run 245064 of LHC15o using the first physics pass

runTrain --name=LHC15o_245064_fp_AOD \
         --class=TrackletAODTrain.C \
         --url="alien:///alice/data/2015/LHC15o?run=245064&pattern=pass_lowint_firstphys/*/AliESDs.root&aliphysics=last,regular#esdTree"

or for simulated data from the LHC15k1a1 production anchored to run 245064

runTrain --name=LHC15k1a1_245064_fp_AOD \
         --class=TrackletAODTrain.C \
         --url=alien:///alice/sim/2015/LHC15k1a1?run=245064&pattern=*/AliESDs.root&aliphysics=last,regular&mc#esdTree

(note the addition of the option "&mc" to the URL argument)

In both cases a sub-directory - named of the name argument - of the current directory is created. In that sub-directory there are scripts for merging the output, downloading results, and downloading the generated AODs.

It is highly recommended to download the generated AODs to your local work station to allow fast second step analysis. To download the AODs, go to the generated sub-directory an run the DownloadAOD.C script. For example, for the real data analysis of run 245064 of LHC15o, one would do

(cd LHC15o_245064_fp_AOD && root -l -b -q DownloadAOD.C) 

and similar for the analysis of the simulated data.

Histogram production

To produce the histograms for the final charged-particle pseudorapidity density , one needs to run a train with a task of the class AliTrackletAODdNdeta (or AliTrackletAODMCdNdeta for simulated data) in it. This is most easily done using the TrainSetup derived class TrackletAODdNdeta.

For example for real data from run 245064 of LHC15o where we store the AODs generated above on the grid

runTrain --name=LHC15o_245064_fp_dNdeta \
         --class=TrackletAODdNdeta.C \
         --url="alien:///alice/cern.ch/user/a/auser/LHC15o_245064_fp_dNdeta/output?run=245064&pattern=* /AliAOD.root&aliphysics=last,regular#aodTree"

or for simulated data from the LHC15k1a1 production anchored to run 245064

runTrain --name=LHC15k1a1_245064_fp_dNdeta \
         --class=TrackletAODdNdeta.C \
         --url=alien:///alice/cern.ch/user/a/auser/LHC15k1a1_245064_fp_AOD/output?run=245064&pattern=* /AliAOD.root&aliphysics=last,regular&mc#esdTree

(note the addition of the option "&mc" to the URL argument)

If we had downloaded the AODs, we can use ProofLite to do this step

runTrain --name=LHC15o_245064_fp_dNdeta \
         --class=TrackletAODdNdeta.C \
         --url="lite:///${PWD}/LHC15o_245064_fp_dNdeta?pattern=AliAOD_*.root#aodTree"

and similar for simulated data

runTrain --name=LHC15k1a1_245064_fp_dNdeta \
         --class=TrackletAODdNdeta.C \
         --url="lite:///${PWD}/LHC15k1a1_245064_fp_dNdeta?pattern=AliAOD_*.root&mc#aodTree"

(note the addition of the option "&mc" to the URL argument)

Final result

The final result is obtained by runnin the class AliTrackletdNdeta2 over the histograms from both real data and simulations. As an example, suppose we ran or histogram production on the Grid and have downloaded the merged results into LHC15o_245064_fp_dNdeta/root_archive_000245064/AnalysisResult.root (using LHC15o_245064_fp_dNdeta/Download.C) and LHC15k1a1_245064_fp_dNdeta/root_archive_245064/AnalysisResult.root (using LHC15k1a1_245064_fp_dNdeta/Download.C). Then we should do

gROOT->LoadMacro("AliTrackletAODUtils.C+g");
gROOT->LoadMacro("AliTrackletdNdeta2.C+g");
p->Run(0x8,0x3FF,9,
"LHC15o_245064_fp_dNdeta/root_archive_000245064/AnalysisResult.root",
"LHC15k1a1_245064_fp_dNdeta/root_archive_245064/AnalysisResult.root",
"result.root");

(see AliTrackletdNdeta2::Run for more information on arguments)

By default, each plot will be made and the process paused. To advance, simple press the space-bar.

If we had made the histograms using ProofLite, we should do

gROOT->LoadMacro("AliTrackletAODUtils.C+g");
gROOT->LoadMacro("AliTrackletdNdeta2.C+g");
p->Run(0x8,0x3FF,9,
"LHC15o_245064_fp_dNdeta/AnalysisResult.root",
"LHC15k1a1_245064_fp_dNdeta/AnalysisResult.root",
"result.root");

Function Documentation

void BackOfTheEnvelope ( const char *  file,
Double_t  c1 = 0,
Double_t  c2 = 0,
Bool_t  mid = true 
)

Do back of the envoloope calculation

Parameters
file
c1
c2
mid

Definition at line 356 of file BackOfTheEnvelope.C.

void DetailsCompare ( )

Definition at line 604 of file DetailsCompare.C.

void Expectations ( Bool_t  mid = true,
Double_t  c1 = 0,
Double_t  c2 = 5 
)

Do back of the envoloope calculation

Parameters
file
c1
c2
mid

Definition at line 505 of file Expectations.C.

void ExtractMCWeights ( const char *  nFile = "dpmjet.root",
const char *  dFile = "hijing.root",
const char *  oFile = 0 
)
Parameters
nFile
dFile
oFile

Definition at line 22 of file ExtractMCWeights.C.

void MakeDeltas ( const char *  realFileName,
const char *  simFileName 
)
Parameters
realFileName
simFileName

Definition at line 21 of file MakeDeltas.C.

void TryBuild ( )

Try to build the code

Definition at line 17 of file TryBuild.C.

Variable Documentation

const char* fmt = "%-10s | %5.3f | %5.2f | %5.2f | %5.2f %5.2f |"

Print format

Definition at line 41 of file BackOfTheEnvelope.C.

Referenced by GridRailway::AuxSave(), Calculation::Particle::Print(), and Calculation::Print().

const char* parFmt = "%-10s | %5.3f |"

Print format

Definition at line 44 of file Expectations.C.

Referenced by Calculation::Particle::Print(), and Calculation::Print().