Expand all   Collapse all Weekly Summaries

Daily Entries

• projects : https://twiki.cern.ch/twiki/bin/view/Main/CSGammaIsolation

todo

• put Firas presentaiton on grenoble twiki
• review cms gamma paper for 12 avril meeting group
• site internet update
• QA trigger

11 juillet

• slides pwgga

Sinon, a priori, pour une demande de papier, pas besoin de presenter la semaien prochaine mais on peut aller direct au forum fin juillet (à confirmer dans la journee par Yuri), ce qui laisserait + de temps pour finir les checks et les propagations au xe + updater la note.

Voici ce que j'ai noté de la presentation d'auj :

slide 4 - renforcer un peu le message physique - attention c'est la fct de fragmentation du parton et pas du photon, mesuré par les directs photons

slide 9 - montrer que la limite inf du fit à 0.27 n'influe pas (changer à qq chose de + loin de la region du signal) : syst negligeable je crois

slide 13 - syst lié à l'efficacité de reconstruction ? si il y en a, varier dans l'analyse certains parametre du MC (j'ai pas trop compris lesquels mais peut etre qu'on peut smearer le pt selon la resolution ou le shifter ?)

questions : 1) est ce que la difference entre les analyses de Mauro et de Marie pour le spectre des photons isolés pourrait avoir un impact sur cette analyse ? la difference est deja au niveau du spectre raw de clusters isolés... => ici ca doit moins affecter de tt maniere parce qu'on utilise des correlations... mais il faudrait le demontrer ? ou alors est ce que ca suffit de montrer qu'on a le meme xE avec la method threshold et sum pt ? (peut etre comparer alors aussi le spectre raw des clusters isolé avec les 2 methodes...) 2) est ce que la future reconstruction des données de 7 TeV avec une meilleurs description des calibrations des detecteurs aurait un impact sur les resultats ? A: ben on a des grandes barres d'erreurs stat là... pas sur qu'on soit sensible... Gustavo ?

Dans les trucs à discuter - pq on retrouve pas la pureté MC du bin pt<0.5 quand on regarde les bin [0.2,0.3][0.3,0.4]et [0.4,0.5] - pourquoi on a une eff differente quand on utilise gj et jj (+ pq vs pt on a gen<reco) - un tour des systematiques

Dans les trucs à faire ou faire tourner (ou pas) - finaliser les plots - finaliser les systematiques - propager au xE - l'analyse avec sum pt <1GeV - l'analyse MC pour pi0 triggué et le bin [0.5,0.6] pour comprendre pq ya un saut ? - d'autres choses pour les systematiques ? - updater la note

9 juillet

analyse complete

pi0

• vieille, 1 run, 159582, LHC11d /Users/catherinesilvestretello/Documents/Alice/work/svn/macros/copyExtractScaleMerge/simu_pp7_JJ_Pi0_AOD/simu_pp7_JJ_Pi0_AOD_Merged.root
• 153371, 153733, 154211, 154383, LHC11c /Users/catherinesilvestretello/Documents/Alice/work/svn/macros/copyExtractScaleMerge/simu_pp7_JJ_Pi0_AOD_New/simu_pp7_JJ_Pi0_AOD_New_MergedAll.root

• trigger 5 GeV (lhc11d) inputData/Simu/simu_pp7_GJ_AOD_new_PtThreshold_MultipleCuts_159582_Merged_RC50_TM1.root
• trigger 5 GeV (lhc11d) /Users/catherinesilvestretello/Documents/Alice/work/svn/macros/copyExtractScaleMerge/simu_pp7_GJ_EMCAL_AOD/simu_pp7_GJ_EMCAL_AOD_Merged.root * trigger 8 GeV (lhc11c) /Users/catherinesilvestretello/Documents/Alice/work/svn/macros/copyExtractScaleMerge/simu_pp7_GJ_AOD_new/simu_pp7_GJ_AOD_new_MergedAll.root

analyse pureté

trigger 5 GeV (lhc11d)

• more bins [0.2,0.3], [0.3,0.4], [0.4,0.5], 2,5] ,[5,10], [10,100]
  • Users/catherinesilvestretello/Documents/Alice/work/svn/macros/copyExtractScaleMerge/simu_pp_7TeV_Pi005GeV_AOD_MultipleCuts_Merged.root
• initial bins >0.5 [0.5,0.7], [0.7,0.9], [0.9,1.3], [1.3,2]
  • inputData/Simu/simu_pp7_JJ_Pi0_AOD_PtThreshold_MultipleCuts_sum_Merged_RC50_TM1.root

trigger 8 GeV (lhc11c?)

• more bins [0.2,0.3], [0.3,0.4], [0.4,0.5], 2,5] ,[5,10], [10,100]
  • /Users/catherinesilvestretello/Documents/Alice/work/svn/macros/copyExtractScaleMerge/simu_pp_7TeV_Pi008GeV_AOD_MultipleCuts_MergedAll.root
• initial bins >0.5 [0.5,0.7], [0.7,0.9], [0.9,1.3], [1.3,2]
  • 8 TeV inputData/Simu/simu_pp7_JJ_Pi0_AOD_new_PtThreshold_MultipleCuts_159582_Merged_RC50_TM1.root
  • 8 TeV inputData/Simu/simu_pp7_JJ_Pi0_AOD_new_PtThreshold_MultipleCuts_sum_Merged_RC50_TM1.root

27 juin

git

• make a clone : http://aliceinfo.cern.ch/Offline/node/2972#patch
• update locally : git pull --rebase
• git status
• ecraser un fichier changé : git checkout -- AliAnaParticleIsolation.cxx

merging Gustavo

• /Users/catherinesilvestretello/Documents/Alice/work/svn/macros/copyExtractScaleMerge

voici un fichier tar avec les suivants fichiers:

* copier les outputs de la grille cpBinsJJPi05GeV.sh : simulations JJ et JJ+Pi0> 5 (LHC11d), pour l’executer regards dedans (executer 2 fois une pour JJ et une autre pour JJ+Pi0) cpBinsJJPi08GeV.sh : simulations JJ+Pi0 > 8 GeV (LHC11d et LHC11c) (executer 5 fois, une fois avec input un run number)

* execute la macro ScaleMerge.C pour chaque un des bins, applique normalisation, merge des bins extractScaleJJPi05GeV.sh (executer 2 fois une pour JJ et une autre pour JJ+Pi0) extractScaleJJPi08GeV.sh (executer 5 fois, une fois avec input un run number)

* Macro merge des bins ScaleMerged.C ScaleMergedNew.C cpBinsJJPi05GeV.sh cpBinsJJPi08GeV.sh

Pour les données, les fichiers merges se trouveront dans: /alice/cern.ch/user/g/gconesab/PartIsoCorr/outputMerged/data_LHC11c_pp7_Pass1_AOD106_PtThreshold_MultipleCuts/mergedall /alice/cern.ch/user/g/gconesab/PartIsoCorr/outputMerged/data_LHC11d_pp7_Pass1_AOD106_PtThreshold_MultipleCuts/mergedall

Seulement JJ et JJ+Pi0 > 5 GeV sont finis, les autres le seront dans l’âpres midi/soiree (soiree pour les données).

Gustavo

20 juin

plots/Purity3_background_MC.pdf

Purity2Fits_Final_IsoPt_NewMC_HightPt_1Ebins.pdf

ChangingBackground_MCMultipleCuts159582_SignalIso05_BackgroundNoIso_Purity3_vsE.pd Purity_multi_ChangingBackground_MCMultipleCuts159582_SignalIso05_BackgroundNoIso_6Ebins.pdf

• check all plots from 20th

• make presentation
  • different purity comparison : ok to go for purity 3
  • hypothesis on the background and purity : extrapolation to 0.25 in data and in MC
  • comparison with MC : alpha value
  • conclusion : value of the purity for our analysis bins
  • conclusion : xe + systematics

17 avril

• compile git : cmake "$ALICE_ROOT" ; make clean;make

20 Fev

Analyse

• pureté 3 slides 7 fev
• analysis status slides

• finir check MC de la pureté
• valider method de la pureté
• checker xE distrib
• systematics

• note

16 May 2014

pureté

• Mc finir :
  • étude sur MC des zones
  • si on se trompait de x% de signal sur le bruit, quel impact sur la pureté ?
  • refaire MC en bin de non iso

• 3 : pureté difference 5%
  • tracer la pureté en fonction de pt range
  • extrapoler à [0.1,0.5]GeV : depenence lineaire
  • suppose que c'est vraiment du bruit (ce qui semble vrai puisque on n'a pas un écart qui grandi à grand pt)
  • verifier que dans [0.5, 0.7] pas de signal en MC
• 4 : à bas pt, photons de décroissance
  • faire ratio lambda02 pour voir à haut lambda02 ce qui se passe
  • zommer à frand lambda02
• 5 ; valeur de scaling NLM2 ? NLM1 ? appliqué 2 à 1 ?
• 6 : haut pt on isole moins
  • bas pt + faux qu'avec les chargés car pureté pas moins bonne

• message qu'on est dominé par la stat sur cette analyse
• la simu ne donne pas d'évolution brutale : données donne
• syst : difference entre vert et extrapolation; valeur central extrapo

21 janv

submit ana.jdl  158285_bin6 simu_pp7_JJ_Pi0_ESD

submit MergeHistos.jdl  simu_pp7_JJ_Pi0_ESD/158285_bin5 /alice/cern.ch/user/c/csilvest/KinematicGen/output/simu_pp7_JJ_Pi0_ESD/158285_bin5/003/AnalysisResults.root

17 Dec

Utiliser la grille (et installer alien sur ton ordi pour lancer les jobs de la maison). Tu as tous mes scripts sur /alice/cern.ch/user/g/gconesab/PartIsoCorr/jdl/

il y a beaucoup de scripts … les principaux sont

• submitRun.sh
• submitMergeRuns.sh
• submitMergeRunsAll.sh (ne pas utiliser dans des simulations avec bins en pTHard)

et le input de ces fichiers est

• run submitRuns.sh firstRunIndex lastRunIndex productionName jdlName
• run submitMergeRuns.sh firstRunIndex lastRunIndex productionName
• run submitMergeRunsAll.sh productionName

Tu as des jdls aussi dans ce répertoire, le plus recent ana.jdl, il faudra changer la version de aliroot/root pour la plus récente qui contienne toutes les motidifications et éviter d'utiliser fichiers par, mais ça aussi marche bien.

tu as aussi des xml sur /alice/cern.ch/user/g/gconesab/PartIsoCorr/collections et les dernières productions s'appellent

• simu_pp7_JJ_Pi0_AOD_new
• simu_pp7_GJ_AOD_new

Tu peux copier les scripts (pas besoin des collections) et les modifier de façon qu'ils marchent chez toi. Si tu ne réussis pas a les faire marcher on verra ça ensemble.

9 dec

http://cern.ch/go/NcS7

sudo curl https://raw.github.com/gerrywastaken/git-new-workdir/master/git-new-workdir /usr/bin/git-new-workdir

sudo chmod +x /usr/bin/git-new-workdir

• installer root avec ./configure  --with-f77=${FORTRAN_COMPILER}  --with-pythia6-uscore=SINGLE --enable-minuit2 --enable-roofit --enable-soversion --with-ssl-shared=yes --with-alien-incdir=${ALIEN_ROOT}/api/include --with-alien-libdir=${ALIEN_ROOT}/api/lib --with-monalisa-libdir=${ALIEN_ROOT}/api/lib  --with-monalisa-incdir=${ALIEN_ROOT}/api/include --with-xrootd-incdir=${ALIEN_ROOT}/api/include/xrootd  --with-xrootd-libdir=${ALIEN_ROOT}/api/lib --enable-gsl-shared --disable-globus --build=debug --enable-builtin-freetype

28 Nov

• lancer l'analyse sur la grille des pi0 nouvelle simu et rajouter la tache AliAnaGeneratoKine sans la leading

19 Nov

• Grid certificate :
  • export from key chain : https://ca.cern.ch/ca/Help/?kbid=060111
  • put in .globus : https://ca.cern.ch/ca/Help/?kbid=024010 and do chmod 400 $HOME/.globus/userkey.pem

9 aout

• je pense qu'il faut pas faire de bin en lambda02 max different
• refaire les macros xE seul lundi et computer les systematiques

29 juiillet

• run PbPb avec trigger || comme il faut
• https://twiki.cern.ch/twiki/bin/viewauth/ALICE/ALICEIsoPhotonHadronPaper

20 juillet

• /cebaf/arbor/ALICE/correlation/DataAnalysis/Gustavo/correlations/Mix_NoL0Cut
• /cebaf/arbor/ALICE/correlation/ThesisPlot/pp/Isolation
• finir PurityRootFile_FitSS.C pour plotter les fit SS et les chi2 pour chaque bin
  • attention binning [ ] ou [ [ !
  • attention la pureté vs. pt a changé...
• finir V2 pour renvoie à Mauro

• check contamination PbPb : merger job + regarder
  • presentation lundi d'apres au PWGGA ?

• check job astrid + lancer simu PbPb : grille !! :S

• rapport soumya

4 Juillet

• create svn PhotonHadronCorrelationsPP macro area in : /cebaf/usr_cebaf/alice/packages/aliroot/trunk/PWGGA/CaloTrackCorrelations/macro and document https://twiki.cern.ch/twiki/bin/viewauth/ALICE/ALICEIsoPhotonHadronPaper

2 juilllet

• The rules for requesting a production are published here

14 juin

• CorrelationsGammaHadron twiki
• la note en l'etat : CorrelationsGammaHadrons_pp7TeV.pdf
• Je l'ai mis la note sous svn : aliPhotonCorr/GHCorrelations7TeVpp
• Je pensais la soumettre dans le systeme une fois que je rajoute les 3 plots qui me manquait, la biblio.
• J'ai mis un plot pour comparer à QM12 (celui de ta derniere presentaiton Nicolas). Une fois que t'aura fini la redaction de ta these Nicolas, on pourra completer en traduisant les plots en anglais (mettant les macros sous svn), rajoutant les plots de Gustavo, ...

10 Juin

• Nicolas :
• Il me manque les figures suivantes Efficiency_Lambda0.pdf , UE_regions.pdf, et faudrait aussi qu'on mette la figure fig:L0_EffPurityv2} je crois.
• Et est ce que tu pourrais m'envoyer ton fichier bib ?

• Est ce qu'il ne faudrait pas rajouter "prompt" dans le titre "Prompt photon-hadron and neutral pion-hadron correlations with the EMCAL in pp [...]"

• section 2 : Est ce qu'on ne doit pas mettre ici les courbes d'efficacité obtenues pour les deux périodes (celles présentées sur la figure \ref{fig:L0Eff_pp} de la section \ref{sec:TriggerPerf} de ta these) ?

• j'ai mis la description sur le MC avant les paragraph sur les particules chargées et neutre parce qu'on utilise le MC dans le texte sur les particules neutres.

• section 2.3 : est ce qu'on montre la resolution des traces ? (see Sec~\ref{subsec:trajectographie} de ta these)

• j'ai introduit le fait qu'on inclura les differences entre les sub samples dans les systematics

• section 3 : du coup en pp on prend R=0.4 et pt thresh <0.5GeV/c ou sum pt <0.1GeV/c

• section 3.2 fin : ne devrait on pas discuter de systematique ici sur l'efficacité d'identification et d'isolation des photons prompts (comme on en a une ensuite dans la partie purity) ? Par exemple des difference entre le lambda02 MC et data ne devrait pas apparaitre ici plutot que dans la partie purity (ou alors c'est redondant avec ce quie st fait pour la fenetre de fit du signal dans l'estimation de la pureté) ? ou en comparant l'efficacité d'isolation obtenu à ce qu'il y a dans la génération (MC thruth) ? Ou juste prendre la difference entre 2 methodes d'isolation ou cones comme une erreur syst additionnelle pour la suite ?

• section 3.3 paragraphe Systematic uncertainties

• pour les systematiques sur la pureté : Ne doit on pas estimer la purité avec d'autre methode d'isolation et cone du coup ? et refaire les variations systematiques decrites ici ?

• Du coup je montrerai bien le fit pour plusieurs bins (quite à les mettres en annexe)

• mettre le plot de beam test lambda02 (see section \ref{sec:SSG4} fig:Lambda0_DataMC)

• table 2 : je crois qu'il faudrait aussi faire ca en bin de Xe puisque c'est ce qui sera utilisé ensuite

• section 4 : pi0 selection

• section 5 : correlations

• partie resultat 5.2 : la coupure en lambda02 c'est 0.1< ou 0.1<=

• changer la soustraction de la contamination par tout le bruit de lambda02>0.5 et pas juste les pi0 (comme c'est deja fait pour l'estimation de la pureté)

• on dit qu'on peut pas trop comparer les correlations obtenues prompt photon-h avec pi0-h à cause des barres d'erreurs... Et en rebinnant c'est pas un peu plus jolie (plot fig 24 droite) ?

• section 6.1 : combinatoric background subtraction in the fragmentaiton function

• section 6. 3 contamination

• table 3: rajouter la systematique de la contamination pour le cas des photon-h

• rajouter un tableau avec les valeurs de la fig 30 ?

3 Juin

• tourner plus de stat MB V1 et V2 : en cours
• tourner analyse purity SS > 0.5 sur MB data pour R=0.3 sumpt>2 sumpt>3, R=0.2 sumpt >2 sumpt >1

30 mai

hadd AnalysisResults_cent010_mergedScaled.root AnalysisResults_Bin*_cent010_mergedScaled.root
hadd AnalysisResults_cent6090_mergedScaled.root AnalysisResults_Bin*_cent6090_mergedScaled.root

hadd AnalysisResults_Bin1_merged.root 000GJ_EMCAL_Bin1*/AnalysisResults_ok.root
hadd AnalysisResults_Bin2_merged.root 000GJ_EMCAL_Bin2*/AnalysisResults_ok.root
hadd AnalysisResults_Bin3_merged.root 000GJ_EMCAL_Bin3*/AnalysisResults_ok.root
hadd AnalysisResults_Bin4_merged.root 000GJ_EMCAL_Bin4*/AnalysisResults_ok.root
hadd AnalysisResults_Bin5_merged.root 000GJ_EMCAL_Bin5*/AnalysisResults_ok.root
hadd AnalysisResults_Bin6_merged.root 000GJ_EMCAL_Bin6*/AnalysisResults_ok.root
hadd AnalysisResults_Bin7_merged.root 000GJ_EMCAL_Bin7*/AnalysisResults_ok.root

28 mai

• -> regarder biais traces, cluster separer
• -> faire pour g-jet embeddé, et lambda02>0.5 bruit donnees relles : R=0.1, 0.2, 0.3, sumpt>1, sumpt>2, sumpt>4
  • j'ai deja g-jet embeddé smpt >1 pour R=0.1, R=0.2, R=0.3
  • g-jet embeddé smpt >3 pour R=0.2
• -> regarder V1 vs V2
• slide 10 : + de bin pour les gaussiennes vs lambda0>0.5
• exprimer le biais en % de l'E soustraite
• varier le pt [10-25 GeV] faire eff vs pt
• tourner l'analyse (avec lambda) sur lambda-jet
• slide 3 : separer cluster et trces
• slide 4 : rouge pas gaussienne fit, plutot poisson a la limite
• seuil sum pt 1, 2, 3, 4 ? regarder la pureté vs. la stat qu'il reste

21 mai

• enlever flag no fQuality ! et +++++GetEventCentrality() -1.0000
• faire azur coloque
• pp : merger 0.2, tourner 0.1 et 0.3 et 0.4

17 mai

• faire le plot mergé scalé bin par bin (couleur) et rejeter les pt qui deccroche (trouver une coupure par bin) puis merger

• • faudra peut etre appliquer une coupure par bin lorsque je tourne... ou plutot filtrer les donnees une fois pour toute ?

• faire multi isolation et tourner sur pp aussi + R=0.4, 0.3, 0.2 pour gg et jj
• faire azur colloque

14 mai

• check twiki julien pour calculatino HV
• pourquoi R=1 pour photon iso alors que pas pour pi0 ?? + check si histo rempli avant ou apres iso -> retourner !
• debugguer code fragmentation
• ecrire AN

26 avril

• working dir : >cd /cebaf/cebaf/EMCAL/cath/PbPb
• macros
  • >source execute_aod : changer le nom du dir set passNb = "PbPb_backgroundSub_embeddingJJ_cone04_fiducialEtaPhi"; et uncommenter ou commenter si on veux tourner sur JJ ou PbPb
  • appel runDistribPC_Embedding_aod_JJ.pl ou runDistribPC_Embedding_aod.pl (GJ) ou runDistribPC_aod.pl (PbPb)
  • tourne /cebaf/csilvestre/work/CaloTrackCorrelations/PbPb/lpsc_running/macros/aodpar/AddTaskCaloTrackCorr_embedding.C où on peux changer
    • la fiducial cut ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.2, 1.8*TMath::RadToDeg(), 2.9*TMath::RadToDeg()) ;
    • le rayon du cone : rechercher ConeSize dans le fichier pour PbPb (2 endroits à changer)
  • appelé par /cebaf/csilvestre/work/CaloTrackCorrelations/PbPb/lpsc_running/macros/aodpar/ana_embedding.C

• puis merge par bin : cf /Users/catherinesilvestretello/Documents/Alice/work/isolation/macros/merging.txt

cd PbPb_backgroundSub_embeddingJJ_cone04_fiducialEtaPhi_try2/
hadd AnalysisResults_Bin2_merged.root 000JJ_Bin2*/AnalysisResults_ok.root
hadd AnalysisResults_Bin3_merged.root 000JJ_Bin3*/AnalysisResults_ok.root
hadd AnalysisResults_Bin4_merged.root 000JJ_Bin4*/AnalysisResults_ok.root
hadd AnalysisResults_Bin5_merged.root 000JJ_Bin5*/AnalysisResults_ok.root
hadd AnalysisResults_Bin6_merged.root 000JJ_Bin6*/AnalysisResults_ok.root
hadd AnalysisResults_Bin7_merged.root 000JJ_Bin7*/AnalysisResults_ok.root
hadd AnalysisResults_Bin8_merged.root 000JJ_Bin8*/AnalysisResults_ok.root
hadd AnalysisResults_Bin9_merged.root 000JJ_Bin9*/AnalysisResults_ok.root

• puis hscale
  • >cd cd PbPb_backgroundSub_embeddingJJ_cone02_fiducialEtaPhi
  • >root -l /data2/alice/csilvestre/CaloTrackCorrelations/PbPb/ScaleMerged.C // check if the right scaling factors are used of JJ or GJ uncommenting the right line
• et puis remerge : /data2/alice/csilvestre/CaloTrackCorrelations/PbPb/PbPb_backgroundSub_embeddingJJ_cone02_fiducialEtaPhi>hadd AnalysisResults_mergedScaled.root AnalysisResults_Bin*_mergedScaled.root

• peut etre uncomment those lines in the AddTask quand on tourne jet jet pour avoir un pt spectra scaled plus plat :

 // Event rejection cuts for jet-jet simulations
  /* reader->SetPtHardAndJetPtComparison(kTRUE);
   reader->SetPtHardAndJetPtFactor(4);
 */  
  // reader->SetPtHardAndClusterPtComparison(kTRUE);
  // reader->SetPtHardAndClusterPtFactor(1.5);
   

PbPb ongoing or todo

• RAJOUTER tach pi0!!
• je tourne sur JJ et GJ embeddés
• coupure fiduciaire eta[-0.2,0.2] phi [1.6, 2.9]
  • JJ R=0.2 /data2/alice/csilvestre/CaloTrackCorrelations/PbPb/PbPb_backgroundSub_embeddingJJ_cone02_fiducialEtaPhi/AnalysisResults_mergedScaled.root
  • JJ R=0.4 /data2/alice/csilvestre/CaloTrackCorrelations/PbPb/PbPb_backgroundSub_embeddingJJ_cone04_fiducialEtaPhi/AnalysisResults_mergedScaled.root
  • GJ R=0.2 /data2/alice/csilvestre/CaloTrackCorrelations/PbPb/PbPb_backgroundSub_embeddingGJ_cone02_fiducialEtaPhi/AnalysisResults_mergedScaled.root
  • GJ R=0.4 /data2/alice/csilvestre/CaloTrackCorrelations/PbPb/PbPb_backgroundSub_embeddingGJ_cone04_fiducialEtaPhi/AnalysisResults_mergedScaled.root

• par contre j'ai les histo vs centrality vides...

• je compte tourner l'analyse sur gamma-jet et jet-jet embeddé
  • pour eta [-0.6,-0.2] [0.2,0.6] autres bandes
    • [on peut pas donner un coupure non symetrique pour eta ? ana->GetFiducialCut()->SetSimpleEMCALFiducialCut(0.2, 1.8*TMath::RadToDeg(), 2.9*TMath::RadToDeg()) ;]
  • et 0.3
  • et pour changer le Eclust à 500MeV et 1GeV il faut reclusteriser non ? Je pense pas que j'aurai le temps... we'll see

23 avril

• rajouter coupure fiduciaire en phi symetrique !
• regarder (sum pt - sum bruit)/sum bruit
• en phi : 2/60GeV -> combien ca fait en surface
• quelle est la somme pt avant correction ?
• comparer plot eta:phi avec les histo AliAnaPhoton, et ceux de la tach QA pour voir si ya aussi le biais pour les clusters
• faire eta:phi pour des energies de clusters differents : est ce que c'est bien uniforme à haut pt
• faire la projection en eta et en phi
• faire 3 bandes en eta: [-0.6, -0.2][-0.2,0.2][0.2,0.6] pour voir
• trouver une methode pour corriger de l'inhomogeneité de detection : regarder dans la tach QA la hit map moyenne des cellules et des clusters, et voir comment ca chance vs E_cluster =300MeV, E_cluster=500MeV, E_cluster=1GeV (hit/accumulated R) -> ya un histo comme ca
• regarder comment il font pour les jets en PbPb (cath va emailer)
• embedding: distinguer les photons et les pi0 : peut on utiliserune coupure sur la somme des pt ?

22 Avril

1. Gamma-jet /data2/alice/arbor/Embedding/Embedding_GJ/simu_pp2.76_GJ_EMCAL_AOD/MB0/000GJ_EMCAL_Bin1/ ... /data2/alice/arbor/Embedding/Embedding_GJ/simu_pp2.76_GJ_EMCAL_AOD/MB0/000GJ_EMCAL_Bin6/

Pour comparaison avec / sans embedding il faut utiliser la simu gamma-jet : /data2/alice/emcal/aliAOD/simu_pp2.76_GJ_EMCAL_AOD/000GJ_EMCAL_Bin1/ ... /data2/alice/emcal/aliAOD/simu_pp2.76_GJ_EMCAL_AOD/000GJ_EMCAL_Bin6/

*1. jet-jet */data2/alice/arbor/Embedding/Embedding_JJ/simu_pp2.76_JJ_AOD/MB0/000JJ_Bin2/ ... /data2/alice/arbor/Embedding/Embedding_JJ/simu_pp2.76_JJ_AOD/MB0/000JJ_Bin9/

Pour comparaison avec / sans embedding il faut utiliser la simu jet-jet : /data2/alice/real/aliAOD/simu_pp2.76_JJ_AOD/000JJ_Bin2/ ... /data2/alice/real/aliAOD/simu_pp2.76_JJ_AOD/000JJ_Bin9/

12 avril

• rajouterhisto ligne 527 pour les traces dans et hors cone (eta,phi), retourner.

Calib

• tourner sur le MC
• tourner pass 0 avec T corrections
• cf plots pass7-pass8 vs runnumber
  • plotter S/N vs runnumber
  • correlation avec fill ? bunchs, magntet ?
  • condition background ?
  • coeffs a, b, c du background vs runnumber

11 Mars

• StageAlban

5 Mars

• StageStefano

4 Mars

• tourner code pour Stefano
• master classes : presentation, doc, ...
• calib : pass6_bis

22 Fev

• trigger QA run : 2.76 TeV pp check QA
• calib : check pass6_bis, refarie plot vs. runnumber, faire nouveau path, relancer encore, check vs. runnumber, soumettre ocdb.

15 Fev

• Pi0Calib :
  • i dont find the right ocdb root file I should sue to creaste the new one with the recalibrator factors

    //Take OCDB used to calibrate existing data, that needs to be updated in
      // alien:/alice/data/year/EMCAL/Calib/Data, take the latest version
      TFile * forg = new TFile("org/Run113461_999999999_v4_s0.root","READ");
     

11 Feb

• jet workshop link
  • yenjie
    • background subtraction
    • slide 20 : c'est quoi la jet shape ? smeared pp ref...
    • slide 22 : check how they do simu
    • slide 23 : syst purity determincation
    • question slide 17 : differentially in R ?
    • slide 29 what do we learn if theory is wrong ?
    • background subtraction sensitive to structures
    • bjets
    • side 15 : which energy ? dependance ?
  • atlas
    • compare jet v2 and single inclusive hadron as they are not expected to be similar

todo

• master classes: faire exo Z, temps, cf twiki
• Calib : faire database file, ordre directories
• QA : check latest runs
• AN read gustavos' note
• preparer stage
• trigger : read

7 Feb 2013

• Presentation end of pi0 calib pdf
• comments :
  • slide 2 : check periodic structure
  • add numbers in slide 3
  • check number of masked towers : if pass1 or production should be more that 150 not 40 !
  • slide 5: sigma = width of mean right ?
  • 6) ask trd people about the asymetry : 2D maps send
  • 7) T corrections applied so ok to produce ocdb
  • wait for gustavo

25 Janv 2013

Jet Framework slides

Port release commit

Salut,

je vais le faire, mais la procédure est la suivant pour la prochain fois :

1) Tu demandes via savannah que le port soi fait

https://savannah.cern.ch/bugs/?group=aliroot

Regards si tu peux accéder a savannah et soumettre un bug, sinon il faut demander a Peter ou Latchezar de t'inclure dans le savannah du projet aliRoot (meme si savannah va être remplace par Jirah a un moment donne)

Tu cherches "Bugs" et tu click sur submit. Tu dois assigner le ticket a Peter Hristov Je t'ai mis en copie dans la requête, comme ça tu as un exemple de comment faire la demande.

2) La requête est prise en compte dans le offline meeting de lundi et si pas de commentaire, elle sera incluse dans le prochaine tag d'aliroot, qui se fait normalement le mercredi.

Parfois, comme les modifs de Michael, le commit et requête est fait juste en temps de faire le tag, et Peter les inclus, ou bien il y a un tag d'urgence, et alors on peut avoir la chance.

3) Une note, dans la requête tu dois mettre tous les commits précédents qui ne sont pas inclus dans le release, et en ordre, dans notre cas, non seulement 60491 mais aussi 60452

8 Nov

• faire plot STU vs V0 pour run p Pb
• retrouver QA results : cf dans aliensh /alice/data/2012/LHC12g/000188362/ESDs/pass2/

31 Oct

• Master Class
• write up offline calibration paragraph
• acces LAF + check jobs
• calibration : check runs and issues, summary of pi0 pass3 status
• site web update
• L1 trigger catch up
• print poster
• retour de mission

oct 24

• daw indico overview training
• training part 1, training part 2, HLT training, CTP trigger training, its page, alice offline

oct 19

• nedit change remplace
  • Dans nedit, rechercher " 23 [0-9]\.[0-9].*$", remplacer par " 23 1.0" (j'espere avec mis les bons backslashs !).

oct 16

• do pass 3 and launch reco
• add integral peak - hors peak tower 0 and 23, echel en z
• oadb bas channel : check fees issu qa warm channels ?
• lhc12d : check in logbook stat triggered

oct 15

• MasterClass2013

oct 12 :

• Fete de la sciences 2012 doc
• ouvrir compte LAF : doc
• pense bete sabine ici
• Calib EMCalPi0Calib2012

oct 8

Correlation

• summarize summer studies
• launch prod of isolation on the new MC : γ (EMCal) – jet : LHC12a15g, jet – jet : LHC12a15f (not biased), Analysis of this production shown before summer, jet – jet : triggered on a particle with pT > 5 GeV/c in EMCal, π0 : LHC12f2a, Hadron (not π0 nor η) : LHC12f2b
• cf HowToRunCaloTrackCorr : cf gustavo's code

sept 20

• todo :
  • ppt avec plot de l'ete sur les corelations
  • liste calib prob pour emcal + illustration
  • calibration
  • site web

sept 10

• rsa : http://www.idris.fr/data/faq/ssh_keys.html

Sept 7

• Calib2012 : twiki
  • tourner sur donnees locales : jsp check si tout a ete transferé; regarder macro pour tourner code /cebaf/cebaf//EMCAL/calibPi0/run_local/runDistribPCNew.pl; tourner code d'analyse de julien ensuite pour voir si 118M event de c ca suffit ou si il faut reconstruire d aussi.

• site internet : actualisé stage et these !

• email Jiri sur simu trigger EMCAl => voir rachid sur ce qu'il manque et coupures

Sept 6

• renew certif
  • sign VO here
  • renew certificate here in firefox, puis export and put in .globus/ with right permissions

$ openssl pkcs12 -nocerts -in mycert.p12 -out userkey.pem
$ openssl pkcs12 -clcerts -nokeys -in mycert.p12 -out usercert.pem
$ chmod 0400 userkey.pem
$ chmod 0600 usercert.pem

Aout 16

• update twiki correlation
• actualise cuts
• list tasks

• Isolation (todo in data and in MC)
  • why is there some differences between isolation fraction between pion identified using the shower shape+split technic and with the inv. mass ?
  • Pile up :
    • Check how the distrib varies in the cone with and without pile up.
    • How to select pile up properly ? Also see the pile up thread on the gamma e-list.
  • for data and MC Check the isolation of pi0 and inclusive neutral clusters as a function of
    • lambda0 cuts
    • the number of local maxima (NLM) (if higher, there is a jet -> + isolated : check if we see this behavior)
  • purity
    • compare isolation with SS purity estimate
• simulation
  • background gamma and triggered => purity
  • stat, bin in pt
  • compare 2.76 and 7 TeV MC
    • UE
    • XE in trigger
    • Isolation
    • distrib in cone, in 0.5 GeV
• systematics : shower shape
• kt charge slope
• eff pp tracks, purity => 7 TeV, trigger bias to correct ?
• pp corrections tracks secondaries, tracks purity, cross-check with Oliviers'
• eff(pi0) vs. pt drops :
  • how does it affect XE ?
  • study generator level : # trigger
• cross-check with Mauro purity compatible ? List of runs ?
• Diphox isolation ?

PbPb

• isolation
  • gamma complicated...
  • pi0
  • subtract background
  • function of Rxnp : * in the cone as a function of the reaction plane (distrib part sum(Econe) vs rxnp) * add also the centrality in the available parameters (distrib part sum(Econe) vs. centrality)
• simulation
  • embedding in PbPb : which centrality ?
• IAA:
  • gamma
  • which purity
  • study sample
  • averaged activity ala Mauro

Aout 2

Voici quelques elements d'information sur la comparaison des FF les plus récentes AKK08, DSS et HKNS. Je propose d'inclure 1 ou 2 lignes d'explications dans la note d'analyse.

Dans l'article (arXiv 0810.1193), F. Arléo présente une comparaison precise des 3 fonctions de fragmentation AKK08, DSS et HKNS.

Les différences essentielles sont resumees dans le tableau 1 p.7. - Les analyses HKK08 et DSS reposent sur les données e+e-/pp et pp_bar ainsi que H_s^0 et Lambda+Lambda_bar. De son cote HKNS, repose uniquement sur l'analyse des données e+e-. - DSS et HKNS possède un traitement des erreurs qui permet d'évaluer la précision des FF ainsi extraites pour un domaine en x donne. La figure 3 (p9) montre les incertitudes sur la FF HKNS. On peut clairement voir que la FF de gluon n'est pas connue avec une grande précision. - Seul AKK08 possed un traitement des corrections des masses de hadrons et de effets de resommation.

En p11, il compare directement les FF : - Les FF des quarks en pi+- sont assez semblables pour les 3 analyses et en particulier pour z>0.3 et en particulier pour es plus grandes valeurs de z. - Pour contre les FF de gluons en pi+- données par AKK08 et DSS sont plus dures que celle donnée par HKNS, ce qui résulte des contraintes introduites par les données de RHIC et du Tevatron. --> Il en conclue que AKK08 et DSS devraient a priori donner des prédictions plus fiable au LHC. - Cependant certaines différences existent a grands z entre AKK08 et DSS, ce qu'il interprète par des données moins précises en pp/ppbar qu'en e+e-. --> Ainsi meme si les mesures a RHIC et au Tevatron aide pour contraindre les FF de gluon, elles introduisent des erreurs significatives.

Remarque supplementaire : - Toutes les analyses les plus récentes montrent que l'accord entre les differents jeux de FF est meilleur a l'ordre NLO qu'a l'ordre LO.

Juillet 24

• comprendre les differences d'eff sur l'isolation des pi0 selon les suggestions de gustavo
  • pile up, tester pileup() dans les donnees 2011 : Ana.C, event selection reader : reject event with pileup : check N event, isolation...
  • MC dans chaque pt hard bin,
  • NLM et shower shape dans données et simulation, ...
• calculate purity (Eff_performance.C) in bins of pt as Nicol slide 4 ici. Ya tj rien sur S/B avec l'isolation du coup. Faudrait que j'essaye qd mm de faire une figure comme la 34 (purity vs. pt) pour voir comment ca se compare en bin de pt plus fin. (d'ailleurs sur cette fig. les box c'est les erreurs syst ? Ya pas d'erreurs stat ? Je ne vois pas où cette figure est commenté dans la note...

Je vais continuer de chercher à Pour tourner avec pileup sur les donnees reelles, vous me conseiller de faire ca sur la grille ou sur les donnees au lpsc ou ?

Juillet 10

• Checker eff pi0 MC (merger, cp ordi, plotter...)
• Checker calculation of errors in S/B, purity, ... Comparer avec SS method, ...
• Corriger note. Comment on fait pour les commentaires ? Tu vas tout repondre et corrigé seul Nicolas ?
  • https://aliceinfo.cern.ch/Notes/node/86
• Je vais aussi refaire les plots eff vs purity fig 23 qui ne sont pas dans les bon bin ni avec le bon MC.
• Et puis updater les syst pour 2.76TeV vs. 7TeV.

Packages = {
    "VO_ALICE@AliRoot::v5-03-38-AN",
    "VO_ALICE@GEANT3::v1-14-2",
    "VO_ALICE@ROOT::v5-33-02b"
};

• dans l'add task ajouter :

  reader->SetPtHardAndJetPtComparison(kTRUE);
  reader->SetPtHardAndJetPtFactor(4);
  
  reader->SetPtHardAndClusterPtComparison(kTRUE);
  reader->SetPtHardAndClusterPtFactor(0.9);

June 30

EMCal

• https://twiki.cern.ch/twiki/bin/view/ALICE/EMCalOnCall

June 15

• comparer isolation Jet-Jet, Gamma-Jet, donnee relles, pour les photons et les pions
• rajouter la decomposition N photon des simu Jet-Jet
• checker pq bin 3 est pas là
• checker pourquoi run 17* n'est pas complet

June 4

• tourner multi isolation sur Jet-Jet MC prod
• tourner Corr sur Real Data pp7 AOD 106
• tourner QA sur runs
• regarder output PbPb
• faire plot ratio denum {0.4, 0.5} sur reste pour quantifier.

• ttaché tu trouveras le par file avec les modifs de Rachid. Dans le ana.C tu dois changer la ligne
  • gSystem->Load("libPWGGAEMCALTasks.so");
  • par celui-ci SetupPar("PWGGAEMCALTasks");
• La façon de créer un par file est
  • make nomLibrairie.par dans le repertoir build

May 31

• jet jet sur grille MC + correlation charged only sur data
  • check :
    • /alice/cern.ch/user/g/gconesab/PartIsoCorr/macros/GJ
    • /alice/cern.ch/user/g/gconesab/PartIsoCorr/jdl/ana.jdl (il faut le modifier pour prender les macros/GJ)
    • mais rappelle toi de activer les analyses pi0 et ses correlations. Cette example est pour Gamma-jet ou faire correlation avec pi0 n'a pas de sense

May 29

• harmoniser filling histo iso et iso cut cxx
• rajouter pt : ptleading et sum pt : pt leading
• re run sur pp
• run grille jet - jet
• check pbpb
• cone 0.2?

Isolation

• Pourquoi relacher la coupure enleve t il des particles isolees ? prob code ou ?

May 20

• QA : running
• slides Firas
• prepa reunion demain (isolation, deadlines): last summary prez, twiki.
• PbPb: QA, plot vs. centrality (cf macro Amaury)
• corriger rapport Firas

Correlation

• /cebaf/balbastre/CorrAnalysis/21_05_2012
• 2 fichiers comprimês avec les histogrammes de l'analyse pp et PbPb
  • ppCorr.tgz :
    • 3 repertoirs per 3 periodes LHC11a, LHC11c, LHC11d, la correction en T est apliquê (mais les corrections manquent pour un dizaine de runs, et même après correction LHC11d a une mass du pi0 trop basse)
    • Dans chaque periode vous trouverez 4 repertoirs indicant diverses coupures de R et pt threshold pour l'isolation R = 0.4, 0.5, pT = 0.5, 1 GeV. Dans R0.4_Pt0.5 vous trouverez les histogrammes de correlation avec et sans isolation, dans les autres que avec isolation
  • PbPbCorr.tgz :
    • Vous trouverez 8 repertoires, chaque un correspondant aux subdivisions annoncêes sur le twiki. Dans chaque repertoire vous trouverez diverses fichers correspondant a divers trigger ou centratlitê
      • TrigCen_ClV1_CenCentral.root : Trigger Central, centralitê 0-20 (no selection)
      • TrigEMC_ClV1_CenCentral.root : Trigger EGA, centralitê 0-20
      • TrigEMC_ClV1_CenPeripheral.root : Trigger EGA, centralitê 20-40
      • TrigEMC_ClV1_CenSemiCentral.root : Trigger EGA, centralitê 60-80
      • TrigINT_ClV1_CenCentral.root : Trigger INT, centralitê 0-20
      • TrigINT_ClV1_CenPeripheral.root : Trigger INT, centralitê 20-40
      • TrigINT_ClV1_CenSemiCentral.root : Trigger INT, centralitê 60-80
    • tout est fait avec clusterizer V1. Dans le Trigger INT ou Central, l'analyse des fake triggers est fait (pour estimer l'efficacitê de la methode). ici la correction de temperature n'est pas faite.

Isolation

• /cebaf/csilvestre/work/CaloTrackCorrelations/macros/ChargedNeutral/Multi/AnalysisResults_ChargedNeutral_multi_PbPb_merged.root
• make QA selection cf twiki gustavo QA
• prob dans la compilation si j'utilise un par file... Surement mes lib etaient trop vielle. Penser á updater regulierement.
• avec ca ca devrait marcher sans par file :

Packages = {
    "VO_ALICE@AliRoot::v5-03-24-AN",
    "VO_ALICE@GEANT3::v1-14-2",
    "VO_ALICE@ROOT::v5-33-02b"
};

May 15

• par files
  • attachê tu trouveras le par file avec les modifs de Rachid. Dans le ana.C tu dois changer la ligne gSystem->Load("libPWGGAEMCALTasks"); par celui-ci SetupPar("PWGGAEMCALTasks");. La façon de crêer un par file est make nomLibrairie.par dans le repertoir build
• QA : je tourne avec une lib d'avril pour avoir si ca va mieu les index, et je teste le par file pour voir si les modifs sont ok.
• EMCAL_OnCallGuide.pdf
• Analisis notes location
• physic forum

May 11

• check jobs (8 qui tournent), puis merger, puis email Romain

• QA
  • lancer la qa pour les runs suivant:

    -177 858
    -177 860
    -177 861
    -177 864
    -177 869
    -177 799
    -177 804 
177 849

• HowToRunQA

May 10 2012

• After meeting
  • generalitês JPL jean-philipe
  • caracterisitques cynthia

April 28 2012

• HowToRunCaloTrackCorr

April 25 2012

• lettre reco firas : ok
• Corr task ca plante ??
• merger QA : ok

April 24 2012

• Verifier # events -> j'ai relancê du coup sur tous les run vec cone 0.5 et pt<0.5
• run correlation for pt threshold 0.5, 1 GeV, cone 0.4, cone 0.5, sum fraction 10% cone 0.4
• donner coupure á Gustavo : ok

• cp outputNoPhysSel/PHYSICS1_17/178025/001/AnalysisResults.root file:AnalysisResults_178025_001.root
• submit MergeHistos.jdl PHYSICS1_17/178030 /alice/cern.ch/user/c/csilvest/QA_trigger/outputNoPhysSel/PHYSICS1_17/178030/001/AnalysisResults.root

Correlation

• V1 SS -> split clusters
• V2 Mass Inv. haut pt
• h-h traces L electrons non coupês
  • correlation change avec contamination

Isolation

• contamination pi0 vs. pt
• comtamination vs. method
• vs. R
• vs. seuil
• changeg legende y : Nphoton
• V2+Mass Inv. et V1+SS
• PbPb : avec ou sens soustraction bruit de fond
• comment estimer la contamination pi0 ?
• combien de gama survivent ?
• eff -> MC

QA

• check reco train if QA task ran properly or not
• objectif : avec run de ref et autres run conclure avec macro qui compare systematiquement

• Isolation cf plots : /cebaf/abedelnabi/Desktop/4Firas/Multi_Cut/Multi_cut_merged/PPR
• QA task : ca tourne le QA_gustavo.jdl... va savoir pourquoi...
  • tourne sur 2 runs :
  • merger : /alice/cern.ch/user/c/csilvest/QA_trigger/jdl/ >run submitMergeRuns.sh 1 1 PHYSICS1_17

• /alice/cern.ch/user/g/gconesab/QA_Trigger/
• dans le repertoir macros tu as le AddTaskQA et ana.C, sur le repertoir NoPhysSel avec les macros qui ont fourni des histos remplis.
• dans le repertoir jdl tu as les scripts pour lancer
  • TriggerQA.jdl : lance le job pour un run pour l'executer : submit TriggerQA.jdl numeroRun production
  • submitRuns.sh : script qui lance les jobs run par run, execute le precedent
  • pour l'executer : run submitRuns.sh first last production
• production : un repertoir que j'ai dans /alice/cern.ch/user/g/gconesab/collection, avec les collections de donnêes par run
  • (first, last sont le premier et dernier run que tu veux analyser en ordre dans le fichier list.txt que tu trouves dans ce repertoir avec les collections)
• Il faudra que tu modifies le TriggerQA.jdl et le submitRuns.sh afin qu'ils pointent vers tes repertoirs.

April 23 2012

• tutorial aliensh grille
• grid certificate : > openssl x509 -noout -in $HOME/.globus/usercert.pem -dates
  • Do I now basically have to go through all steps at "How To get a certificate at CERN" on the alien webpage?
  • get a cerificate from your local CA
  • Do not forget to register with your new certificate in VOMRS https://lcg-voms.cern.ch:8443/vo/alice/vomrs

pp

• preparer des fichiers root avec different parametres d'isolation pour les correlations : attention á bien choisir les pt hadron, cuts, ...
• (Amaury) etudier l'impact de la contamination sur le xE, systematiques
• (Firas) conclure sur l'isolation en pp

• tourner jobs avec coupure bruit de fond de 2GeV
• tourner avec soustraciton de bruit de fond : mais induit plus de perte uncertitude que de gain (up to 12GeV !); tester cone á meme etat
• etudier l'isolation

• tourner tache QA pour avoir histos de ref et aboutir á macros de systematisation des resultats QA
• simu

April 20 2012

• todo : checker job sur grille et 3 interactif. mettre la macro de nico en mode grille pour l'executer + vite ou alors faire un script qui lance tout seul.
  • interactif : trop long ! ca a plantê a cause du token
  • j'ai relancê sur 100 segments...
• check correlations
  • amaury essaye de reproduire les plots de gustavo. qq remarques et questions :
  • photons : V1 vs. V2 change nombre d'entree mais pas la forme des correlations
    • peak away side existe tj ?
  • pi0 : V1 moins de stat que V2, normal
  • normalisation isolê : checker -> essayer avec l'integral partielle
  • SS mieux haut pt ? chaque clusterisers, aussi pour V1 a bas pt
  • pq hh away side + grand qye photon-hadron ?
• hh-> V1 = V2 ?
  • Comment les hadron sont ils reconstruits
  • existe trigger pt info ? oui delta phi vs. pt leading -> faire graph 8-15, 16-40
  • photon-hadron = pi0-hadron ?
  • c'est quoi spllot ?
  • c'est quoi xE ?

• á la page de monalisa á gauche tu as une espèce de vêlocimètre, tu vois que il est presque au maximum ...

• sinon tu peux contrôler tes jobs directement avec

https://alimonitor.cern.ch/users/jobs.jsp

• http://pcalimonitor.cern.ch/show?page=jobStatus.html

April 19 2012

• aliensh
  • submit TriggerQA.jdl 176753 PHYSICS1_17
  • ps
  • top -user csilvest
  • ps -trace 173129386 all
  • ls /alice/data/2012/LHC12a/000176753/raw/ | tail -2 | awk '{print "alien:///alice/data/2012/LHC12a/000176753/raw/"$1 }'
  • find -x collection /alice/data/2012/LHC12b/000178025/ AliESDs.root >../../collection/PHYSICS1_17/178025.xml /alice/cern.ch/user/c/csilvest/QA_trigger/jdl/ >submit TriggerQA.jdl 178025 PHYSICS1_17

• • [lpsc2035x.in2p3.fr]/cebaf/csilvestre/work/QA/macros>aliroot -l -b -q rootlogon.C STU_AMORE_Grid.C\(\"List_177750.txt\",\"RawResults_177750.root\"\) > & List_177750.log &
  • lpsc2030x:~/work/QA/macros> aliroot -l -b -q rootlogon.C STU_AMORE_Grid.C\(\"List_177798.txt\",\"RawResults_177798.root\"\) > & List_177798.log &
  • lpsc2036x:~/work/QA/macros> aliroot -l -b -q rootlogon.C STU_AMORE_Grid.C\(\"List_177792.txt\",\"RawResults_177792.root\"\) > & List_177792.log &
• QA
  • aliroot rootlogon.C STU_AMORE_Grid.C

76753 2h Pas bon (dêcorrêlation Lo//L1)
177750* moins de stat mais L0 vraiment nickel
177792* 2h OK pour reference ?
177798* 1h OK pour reference ?
177942 31 min  Pas bon (1 hotspot en L0 plus marquê et peu de stats)
178056 51 min Peu de stats
178058 2h  Peu de stats
178062 37 Peu de stats
178025 19 min LHC12b  plot L0 bizarre
178030 2h LHC12b 1 TRU bizarre (L1 Gamma Trigger) 

April 17 2012

• Isolation plots grouped here pdf
• Correlation : je sais pas si t'as jetê un oeil á la presentation de Gustavo ici mais peut etre ce que tu pourrais essayer de reproduire dêja comme plots ce sont ceux des slides 12->16 ? Qu'en penses tu ? Tu devrais avoir moins de stat mais je pense que t'as tout ce dont tu as besoin dans le fichier root.
• logbook alice shifter
• alimonitor monalisa

April 12 2012

• presentation firas ici
  • changer Y name histo
  • tester cone 0.5
  • V1 vs. V2 : rejette pi0 avec v2 jusqu'a 20GeV
  • V1 + shower shape: existe sample + pur pour denum
  • comparer V1 et V2 spectres, ou integral apres 10GeV
  • V1 + unfolding ?
  • tester 2D avec seuil á 0.1 GeV
  • faire le 2D avec les barres d'erreur

• reunion de group
• decay empty : coupure mass ?
• dqm seuil :
  • amplitude QA moyen das fastir oiyur ce ryb au niveau 0 et niveau 1
  • regarder dans le train si les donnees sont reconstruites
  • demander reco de ca et ca
  • regarder les donnees 2011 : reco a nouveau dans les runs
  • bon run de l'annee derniere vs mauvaus run

April 10 2012

• Thanks to gustavo : AddTaskIC.C works (j'avais zappê de cp dans le .pl) : removes les photons isolês en trop
• rajouter R=0.5 ?
• Ca tourne pour V1, V2, 4 methods : merger, faire plots 1D, 2D d'eff...
• Firas : presentation en cours. Reflechir aux ccl. Cf PPR

• QA : regarder les fichiers de gustavo : faire une macro pour comparer les runs... vs. run de ref ou moyenne ? ou ?

• resumê cms paper ?

Apr 5

• debugging multi : sum tj empty (cf threshold) -> run sur all runs !
• QA : email Marie OK, doc DQM OK (umm revoir qd mm non ?)
• presentation QA pour le 12 avril

April 2 2012

Comments

• review comment https://aliceinfo.cern.ch/ArtSubmission/node/115 for friday

Isolation

• analyse
  • prob ptSum
  • refaire multi avec kDensity
  • plot eff isolation eta:phi pour chaque methode

• pour l'instant
  • utilisant les infos chargés, neutres, ou chargés+neutre
  • clusturiser V1 vs. V2
  • cone size {0.3,0.4}
  • method {ptThresh, ptFrac}
  • 4 seuils pour chaque methode (ceux qu'il y avait par default... jsais pas si cest les + pertinent)
• Je rajouterai la methode ptSum (des que de trouve pq les histo sont vide), et la nouvelle methode density (la prochaine fois que je tourne).

• regarder la fraction d'isolation et le nombre de photon restant.
• pour chaque clusteriser, cone size, chaque methode, les 4 seuils (cad sur un mm plot 1 clusterizer, 1 cone size, 1 methode avec 4 seuils et les neutres/chargés/neutre+chargé).

• prob plotting colors et canvas when more than 1
• prob for pt sum
• prob for pt frac
  • check isolationCut.cxx

Mars 30

Isolation

• Les fichiers avec les historammes comportant tous les cas sont ici :

/cebaf/csilvestre/work/CaloTrackCorrelations/macros/AnalysisResults_ChargedNeutral_multi.root
/cebaf/csilvestre/work/CaloTrackCorrelations/macros/AnalysisResults_OnlyNeutral_multi.root
/cebaf/csilvestre/work/CaloTrackCorrelations/macros/AnalysisResults_OnlyCharged_multi.root

• faire des plots lisibles : scp FracIsoPartCone_graph.C csilvestre@lpsc2021xNOSPAMPLEASE.in2p3.fr:/cebaf/csilvestre/work/CaloTrackCorrelations/macros/.
• retourner la nuit avec kDensity et checker le ptSum...
• faire aussi les plot 2D

Mars 28

todo

• launch multi analysis : ok

Mars 23 2012

Todo isolation

• checker les barres d'erreurs
• tourner sur toute la stat !
• faire varier tous les paramettres :
  • learn about Multi
  • rajouter les histos qui manques
• comparer les histo 2D eta:phi pour les spectres isoles : neutre vs. chargês vs. both
  • voit t on des zone mortes ?
  • plotter les erreurs stat associês
  • plotter eta:phi 4x4 cellules (ou SM par SM)
  • comparer les points chauds avec la meme cartographie avec les clusters
  • ou varier le threshold (0.5GeV) pour voir l'impact
• V1 : les clusters sont + isolês
• V2 : les clusters de gamma sont moins isolees car mergês plus tard
• pq V2 est plat ? --> voir MC...
• limite en pt des decay (et pion) vient qu'au dela on ne peux pas appliquer la methode de la masse invariante
• les decays ont 1 seuls photons (mass inv)
• les pi0 sont 2 photons (mass inv)
  • avec la shower shape ca donne quoi ?
  • pour les chargês et neutres, + E est grand, + il y a des particules dans le cone

Isolation

• Summary presentation

Mars 20 2012

Group meeting

• Correlation : reactiver l'AN ? (run selectoin , choix des cuts, tuning MC)
• check sur les runs que j'utilise n'ont pas de prob avec la TPC
• check Shower Shape
• plotter le nombre de photons avant et apres isolation pour chaque methode
• Decays : mass invariante pour id que c'est un photon decay de pi0

Isolation

• refaire avec kPtThresIC et sizeCone=0.4
• tourner SumPtThresIC thresh=1.0
• touner multi IC :
  • varier methode cone cut, thresholds, conesize
  • check histograms 2D :
    • maybe need to add histo of photon decay
    • check aliIsolationCut.h -> SumPtThreshold

• sumPtFrac = 10%, et ConeSize = 0.4

Mars 19 2012

Isolation

• refaire avec une autre coupure sumPtFrac, et ConeSize = 0.4 -> ungoing
• V2 marche mieu pour les decay et pion (?), si on utilise un coupure moins strcit -> a checker
• pour l'instant, ya pas bcp de stat dans les decay et pion isoles lorsqu'on utilise l'info neutre : a cause de la coupure trop strict
• conclusion pour l'instant c'est mieu d'utiliser charge+neutral parce qu'on a moins de contamination. V1 remonte a partir de 20GeV parce que ya un truc en plus (cf gustavo)
• faudra confirmer avec la simu, et peut etre regarder pbpb mais alors faut implementer la coupure du bruit de fond en plus.

	kPtThresIC=0.5 coneSize=0.3	kSumPtFracIC coneSize=0.4
photon
	kPtThresIC coneSize=0.3	kSumPtFracIC coneSize=0.4
Decay
	kPtThresIC coneSize=0.3	kSumPtFracIC coneSize=0.4
Pi0

Work

• use liste_pc : 2017x, 2030x, 2035x, 2036x

Mars 16 2012

Isolation

• Les histos qu'on regardaient sont les suivant dans FindObject():

dans le directory
    TList *list = (TList *) f->Get("EMCAL_TrigEMC7_ClV1_Ecell50_Eseed100_DT0_WT0");


    hPhotonNoIso[file]      = (TH1F *) list->FindObject("AnaIsolPhoton_TM0_hPtNoIso");
      hPhotonIso[file]      = (TH1F *) list->FindObject("AnaIsolPhoton_TM0_hPt");
     hPhotonNoIso[file]      = (TH2F *) list->FindObject("AnaIsolPhoton_TM0_hEtaPhiNoIso");
      hPhotonIso[file]      = (TH2F *) list->FindObject("AnaIsolPhoton_TM0_hEtaPhiIso");

      hDecayPhotonNoIso[file]      = (TH1F *) list->FindObject("AnaIsolPhoton_TM0_hPtDecayNoIso");
      hDecayPhotonIso[file]      = (TH1F *) list->FindObject("AnaIsolPhoton_TM0_hPtDecayIso");
      hDecayPhotonNoIso[file]      = (TH2F *) list->FindObject("AnaIsolPhoton_TM0_hEtaPhiDecayNoIso");
      hDecayPhotonIso[file]      = (TH2F *) list->FindObject("AnaIsolPhoton_TM0_hEtaPhiDecayIso");

      hPi0NoIso[file]      = (TH1F *) list->FindObject("AnaIsolPi0_TM0_hPtNoIso");
      hPi0Iso[file]      = (TH1F *) list->FindObject("AnaIsolPi0_TM0_hPt");
      hPi0NoIso[file]      = (TH2F *) list->FindObject("AnaIsolPi0_TM0_hEtaPhiNoIso");
      hPi0Iso[file]      = (TH2F *) list->FindObject("AnaIsolPi0_TM0_hEtaPhiIso");


Je vais rajouter un directory pour V2 qui sera comme ca :
    TList *list = (TList *) f->Get("EMCAL_TrigEMC7_ClV2_Ecell50_Eseed100_DT0_WT0");

• On pourra comparer les fractions d'isolation pour les photons, les pi0 et les decay, pour v1 et v2, en fonction de pt et d'eta:phi pour en parler lundi aprem avec Nicolas et Gustavo et discuter pour savoir ce qu'ils en pensent et ce que tu pourrais regar

• tourner V2 pour avoir plus de decay et pi0 : 3 jobs lances avant de partir dans les directory respectif. Reste a merge et a plotter.
• si neutre only, mettre a 0 la coupure sur delta cone : gustavo a update le code
• dans quoi s'insert l'analyse ? aider isolation... to check with nico et gustavo
• check : what is Neutral, Charged, how is the matching done, how are pi0, decays, iso, photons selected. => regarder dans le code
• faire chtite presentation pour mardi

Mars 15 2012

todo

• reparter lib ordi pour tourner dessus !

Isolation

• je tourne charge+neutral dans 4Firas, onlyNeutral dans /OnlyNeutral/ sur lpsc2006x
• et charged only dans macros/
• reste á
  • merger histograms
  • tourner neutral only
  • bidouiller AddTask pour ne pas avoir á lancer 3 jobs ! arg... add name in the histo,et in the ().
• macro d'analyse : faire avec TGraphError pour voir difference et pouvoir binner

Mars 14 2012

Isolation

• utiliser neutre ou pas dans le cone ou neutre+charged
• cut fidutial
• fraction d'isolation vs. cone neutre, chargê, ou both, en eta:phi et pt pour les photons, les decays, et pi0

• dans 4Firas/ ca tourne ./anaAllAOD.sh LHC11d_pp7_pass1 3 300 avec charged+neutral et photon, decay, pi0, iso and not iso
• todo
  • merger Analysis output
  • tourner pour charged only, et pour neutral only, puis merger

Mars 13 2012

Starting to analyse

• tourner sur + d'evenemnt ana.C qui ne tourne qu'une tache d'analyse et fait moins d'histo
• merger
• regarder les histos

• pour Firas :
  • J'ai modifie le directory : /cebaf/csilvestre/work/CaloTrackCorrelations/macros/4Firas pour y copie les seuls fichie dont tu as besoin pour tourner le script anaAllAOD.sh
  • J'y ai mis le fichier anaAllAOD.sh ou il y a seulement la creation des liens symboliques (ln -s) vwers les fichiers necessaire lorsque on tourne sur chaque input (cad dans le directory LHC11d_pp7_pass1/###/)
  • Tu peux donc copier a nouveau mon dossier 4Firas=/ et essayer d'executer : => ./anaAllAOD.sh LHC11d_pp7_pass1 0 0
  • Tu verras la creation de LHC11d_pp7_pass1/000156797/
  • et le fichie de sortie root : AnalysisResults.root (entre autres)
  • Pour faire les plots, tu peux soit ouvrir dans root le TBrowser comme on a fait hier (mais apres c'est un peu laborieux pour sauvegarder les fichiers je trouve et faire des plots plus compliquer comme des efficacitees ou superposer des courbes). J'ai mis un exemple de macro root que tu peux t'amuser a regarder / modifier a souhait pour y mettre d'autres plots : /cebaf/csilvestre/work/CaloTrackCorrelations/macros/4Firas/Analyse.C
  • Une fois que t'as execute > ./anaAllAOD.sh LHC11d_pp7_pass1 0 0
  • tu peux tourner la macro d'analyse pour voir : root -l Analyse.C
  • Ca plot le pt des photons, et photons isoles, et decay photons que l'on a analyses, et en dernier une efficacite d'isolation.
  • Ca sauve ca ici plots/Photon_Decay_Eff_pt.pdf

Mars 12 2012

Tourner ana.C

• /cebaf/arbor/ALICE/correlation/NewTrunk/ana_Local.C

Mars 9 2012

Getting started with analysis Isolation

• check SetupAlice from arbor
• rot Ana_Local.C

1. Ana.C : lance taches d'analyse

• load CaloTrackCorr : pi0, gamma selection, choix correlation, tracks
• EMCALClusterize : V1 default, ou V2

1. tourner local : ppLHCD11d, ou donnees PbPb
2. recup from Gustavo : /data2/alice/balastre/Analysis/anaAllAOD.sh et anaAOD.sh (et les autres AddTaskCaloTrackCorr par exemple)
3. AddTaskCaloTrackCorr : choix des histos

• regarder eta, phi, energy
  • pour les photon et pi0, isolês ou non isolê
  • correlations delta phi et xe, isolê vs. no isolês
• check chez arbor NewTrunk/AnalyseOutput.C

1. Commenter pour ne regarder qu'une tache d'analyse, ex evenment triggees EMCal dans Ana.C
2. addTask, en haut, check maker->AddAna ... et ne garder que AnaPhoton et AnaIsoPhoton par exemple
3. changer cuts dans configure et voir ce que ca change
   • quelle methode garde le plus de photon direct et rejette le plus de pi0 ?
     • en pbpb faut tenir compte du brruit de fond
   • comment change l'eff en changeant les cuts : comment varie id pi0 ? est elle biaisêe ? check vs eta, pt, ...
   • eff pi0 isolation en pp?
     • phton decroissance 5%
     • haut pt : shower shape eff merdouille ?
     • pi0 mass inv, existe 1 seul cluster
     • eff 20% pi0
     • pq eff isolation photon est meilleur quand le pt>12GeV ?
     • eff photon direct 95% from simu
     • eff photon decroissance : 20% de contamination, 5-10% avec les traces
       • il y a + de pi0 a grand pt qu'on n'arrive pas a remover parce que isolês?
       • proportion des photon de fragmentation ?

Todo

• listen to EMCal meeting
• continue with tutorial and checking code

Mars 8 2012

Todo

• discuss https://twiki.cern.ch/twiki/bin/viewauth/ALICE/PhotonHadronCorrelationTasks
• start looking at code
  • esd aod : code info
  • ex code $ALICE_ROOT/EMCAL/macros/TestESD.C or TestAOD.C
  • script setup env. lscp ?
  • lancer jobs ?
  • se logguer sur ssh -X csilvestre@lpsc2021xNOSPAMPLEASE.in2p3.fr
  • alien-token-init csilvest
  • tutorial
    • alien
    • aliroot manual
• check intro theses : gustavo
• rajouter Firas dans la liste du groupe
• computer Firas ?

Mars 6 2012

todo

• lister taches d'analyse

Correlations

• tasks list twiki

Mars 5 2012

todo

Stagiaire

• tour administratif
• intro par Nico

Analysis

• meeting de vendredi sur evo
• etude du twiki de gustavo
• où est le code ? voir Nico/Gustavo

Correlation and photons

• twiki https://twiki.cern.ch/twiki/bin/viewauth/ALICE/PhotonHadronCorrelations
• issues https://twiki.cern.ch/twiki/bin/viewauth/ALICE/%20PhotonHadronCorrelationsStatusProgressIssues
• from the PWG-GA working meeting https://indico.cern.ch/conferenceDisplay.py?confId=180529

• MC runs for reference
• Need of simulation
• study possible trigger particle
  • Inclusive (isolated) photons : shower shape + track veto (+isolation)
  • Inclusive(isolated) pi0/eta : shower shape + track veto (+isolation)
    • Selection via pairs inside invariant mass window : Need to check the optimal cuts : mass range, time cut, opening angle
    • Selection via pairs outside mass window : Need studies
    • Selection of clusters in Shower shape pi0 band : To be studied, selection band refined
    • Selection of large clusters and splitting (cleaner than previous) : To be studied, refined
• several isolation methods studied
• Need to study the necessity of fiducial cuts on trigger particle, depending on cone size
• Need to work on mixing with AODs (ESDs) containing different triggers
• Need to understand how to estimate systematic errors (check what PWGCF do)
• use energy density instead of just the energy to remove bad channels and acceptance effects
• why isn't the acceptance issue not visible in phi ?
• Using photon isolation method. More details Ref. PRD. 82,072001,2010

• isolation
  • comprendre pourquoi pas d'exces de photons isoles en phi quand on ne met pas de coupure fiducial (alors qu'il y en a un en eta)
  • essayer avec des cells plutot que des cluster
  • essayer une coupure sur la densite d'energie dans le cone (Energie tot / # cells OK) pour tenir compte des bad channels. En mettant une coupure sur la fraction de surface du cone qui doit être "bonne" (par ex : plus de 75% du cone doit contenir des cells good) on peut aussi jouer sur la coupure fiduciale.
• mixed event (xiangrong)
  • si on fait du mixed avec des events triggês par EMCal, il semble qu'une correlation apparaisse a cause de la geometrie limitêe de ces evenements (le mixed n'est pas plat).
  • comparaison du mixed entre bons et mauvais runs.
• Taches attribuêes :
  • dana/xiangrong doivent essayer de faire un QA des runs (tracks) pour proposer une liste
  • Alexis doit regarder si la calibration des runs peut jouer sur l'isolation
  • Je dois regarder le biais d'acceptance en testant le fait d'isoler avec des cells et la methode de la densitê.
• Par contre pour ce qui est d'appliquer l'isolation sur les cells (ce qui me semble être une bonne idêe car on n'est plus biaisêe par la clusterisation) j'aurais besoin de ton aide pour faire ca propre (par exemple ajouter une option dans l'isolation Cells/Clusters). On pourra regarder ca apres le CNRS ?

• checking https://twiki.cern.ch/twiki/bin/viewauth/ALICE/PhotonHadronCorrelations
  • http://aliweb.cern.ch/Offline/Activities/Analysis/PWGDocumentation/PWG4/PartCorr.html

• egroups https://groups.cern.ch/Pages/default.aspx

August 31 2011

Lire les AOD

• voir /cebaf/arbor/ALICE/stagiairesL2/clusterization/analyse_emcalReclusterize.C
  • Il n'y a pas grand chose mais cela montre comment lire les AOD. Apres tu n'as plus qu'a regarder dans la classe AliAODCaloCluster pour voir toutes les infos que tu peux recuperer dans les AOD.
• Un autre exemple interessant est le code de Diane : /cebaf/faivre/labo/stageM2Diane20102011/dianephy/analyseEmcalReclusterize10.C
  • La dedans tu dois pouvoir trouver des infos sur la lecture des clusters, mais aussi celle des tours (la classe est AliAODCaloCells).

August 25 2011

• Trying to run emcalSimReclusterize.C after adding AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler()->SetFillAOD(kTRUE); in AliAnalysisTaskEmbed.cxx (from /cebaf/csilvestre/work)
  • error from the output container :

E-AliAnalysisTaskEmbed::CheckOwnership: ####### IMPORTANT! ####### 
                Task task (AliAnalysisTaskEmbed) posts a container that is not owner at output #1. This may apply for other embedded containers. 
                ####### FIX YOUR CODE, THIS WILL PRODUCE A FATAL ERROR IN FUTURE! ##########
MGR: Processing event #0
 *** Break *** segmentation violation

August 17 2011

• Jet
  • AddTaskJet : macros originale utilisee pour la reconstruction des jets.
  • anaJet.C : macros pour lancer AddTaskJet (que j'ai directement integre a la macro).
  • Il ne reste plus qu'a modifier le directory au debut de la macro pour indiquer un repertoire ou tu auras copie un echantillon d'AliAOD.root.
  • Tu peux aussi aller voir dans : aliroot/PWG4/macros pour trouver les macros officielles.

August 16 2011

• catching up on mails and todo

• jet trigger analysis
• ALiroot, install play

notes from some meeting

• fonctionement algo en pp et pbpb
• discuter avec Magalie et Matseutz : etude Emcal JetAn pour integration Emcal dans reco jet
  • emailer Magalie : code Nantes
• JetAn
  • Existe doc
    • Alexandre, Marco, Adam
    • macro reco jet chargê
  • courbe eff jet
• reco run avec Jet : info tours calo pour simuler patch : pseudo courbe eff
  • mettre l'info trigger dans les AOD
  • modifier code local avec copie des ESD dans AOD
    • cf Rachid ou Gustavo
• rentree meeting photon

1. Simu trigger
2. Simu embedding

July 22 2011

• 14h rdv lpsc
• rapport annuel
• installation aliroot
• check https://indico.cern.ch:443/conferenceDisplay.py?confId=147060
• check mails todo

July 21 2011

• EPS-HEP

July 20 2011

todo

• corriger talk eps y pratcise
• A general rehearsal will be organized on Wednesday 20 July starting early (10:00) in the morning at Alpes-Congres (access map on our site). MAKE SURE TO BE THERE AS THIS WILL BE THE ONLY OCCASION FOR EACH TEAM TO REHEARSE IN REAL CONDITIONS BEFORE THE CONFERENCE.
  • http://eps-hep2011.eu

July 19 2011

• twiki alice grenoble https://twiki.cern.ch/twiki/bin/view/Main/AliceGrenoble

July 18 2011

• eps quarkonia CMS in HI (talk)

Group analysis meeting

• Comment checker les bons runs ?
• c'est quoi le roll back ?
• eff = 2x2 clusters triggê / chaque cluster, par TRU

• study Jet funder
  • info pas dans les aod
  • lancer code sur ESD
• PW4Jet
• ESD trigger
  • reco energy dans patch en pp
  • recalculer index
• mettre EMCal dans reco
• regarder acceptance du jet chargê

• regarder Jet finder
  • code a faire
  • tourner sur un run
  • train de Christian : cf sur la liste jet
  • puis recuperer info L1
    • eff -> existe patch ?
  • changer ref EMCal pour ref Alice
• Magalie : format AOD ? de quelle facon patch info besoin ?
• Demontrer que en pp trigger march + facteur de rejection
• Embedding : lancer job sur la grille

-- CatherineSilvestre - 19-Jul-2011