+++ /dev/null
-#!/bin/bash
-
-# ana.sh
-#
-#
-# Created by schutz on 09.04.08.
-# Copyright 2008 ALICE. All rights reserved.
-
-#Analysis mode, 0 local, 1 localCAF, 2 PROOF, 3 GRID
-export MODE=$1
-
-#Local mode parameters
-#Data root directory
-export INDIR=$2
-#Pattern of data directory name
-export PATTERN=Run
-#Number of files to analyze
-export NEVENT=$3
-
-#GRID mode parameters
-#name of xml file
-export XML=$2
-
-#Do the analysis
-alienroot -b <<EOF
-.L ana.C
-ana($MODE)
-EOF
+++ /dev/null
-//------------------------------------\r
-// Configuration macro example:\r
-//\r
-// Configure EMCal electron analysis\r
-//\r
-// Modified by: K. Read\r
-//\r
-//------------------------------------\r
-\r
-AliAnaPartCorrMaker* ConfigAnalysis()\r
-{\r
- //\r
- // Configuration goes here\r
- // \r
- printf("======================== \n");\r
- printf("ConfigAnalysisElectron() \n");\r
- printf("======================== \n");\r
- Bool_t kInputIsESD = kTRUE; //uncomment for input ESD\r
-//Bool_t kInputIsESD = kFALSE; //uncomment for input AODs\r
- Bool_t kFollowsFilter = kTRUE; //uncomment if follows ESD filter task\r
-//Bool_t kFollowsFilter = kFALSE; //uncomment if no ESD filter task\r
- Bool_t kMC = kTRUE; //set to kFALSE for data\r
-\r
- //enum for the different electron cut sets\r
- //defined for dR and p/E\r
- //kTight2 is the default standard cuts\r
- enum kCutTypes {kTight2, kLooseTight, kTightLoose, kLoose2};\r
- Int_t kCutSet = kTight2;\r
- Double_t pOverEmin = 0.8; //tight\r
- Double_t pOverEmax = 1.2; //tight\r
- Double_t dRmax = 0.02; //tight\r
- if (gSystem->Getenv("ELECUTSET")){\r
- kCutSet = atoi(gSystem->Getenv("ELECUTSET"));\r
- }\r
- if(kCutSet == kLooseTight) {\r
- pOverEmin = 0.6; //loose\r
- pOverEmax = 1.4; //loose\r
- dRmax = 0.02; //tight\r
- }\r
- if(kCutSet == kTightLoose) {\r
- pOverEmin = 0.8; //tight\r
- pOverEmax = 1.2; //tight\r
- dRmax = 0.05; //loose\r
- }\r
- if(kCutSet == kLoose2) {\r
- pOverEmin = 0.6; //loose\r
- pOverEmax = 1.4; //loose\r
- dRmax = 0.05; //loose\r
- } \r
-\r
- //Alternatively, select input via anaInputData environment variable.\r
- if (gSystem->Getenv("anaInputData")){\r
- TString kInputData = gSystem->Getenv("anaInputData");\r
- if( kInputData == "AOD" ){\r
- kInputIsESD = kFALSE;\r
- kFollowsFilter = kFALSE;\r
- }\r
- }\r
-\r
- //Alternatively, adjust for real data based on kMC value.\r
- if (gSystem->Getenv("anakMC")){\r
- kMC = atoi(gSystem->Getenv("anakMC"));\r
- }\r
-\r
- //Detector Fiducial Cuts\r
- AliFiducialCut * fidCut = new AliFiducialCut();\r
- fidCut->DoCTSFiducialCut(kFALSE) ;\r
- fidCut->DoEMCALFiducialCut(kFALSE) ;\r
- fidCut->DoPHOSFiducialCut(kFALSE) ;\r
-\r
- //fidCut->SetSimpleCTSFiducialCut(0.9,0.,360.);\r
- //fidCut->SetSimpleEMCALFiducialCut(0.7,80.,190.);\r
- //fidCut->SetSimplePHOSFiducialCut(0.13,220.,320.);\r
-\r
- fidCut->Print("");\r
-\r
- //----------------------------------------------------------- \r
- // Reader\r
- //-----------------------------------------------------------\r
- if(kInputIsESD && !kFollowsFilter)AliCaloTrackESDReader *reader = new AliCaloTrackESDReader();\r
- else AliCaloTrackAODReader *reader = new AliCaloTrackAODReader();\r
- reader->SetDebug(-1);//10 for lots of messages\r
-\r
- //Switch on or off the detectors information that you want\r
- reader->SwitchOnEMCAL();\r
- reader->SwitchOnCTS();\r
- //reader->SwitchOffEMCALCells(); \r
- reader->SwitchOffPHOS();\r
- //reader->SwitchOffPHOSCells(); \r
-\r
- //Kine\r
- if(kMC && !kInputIsESD){\r
- reader->SwitchOffStack(); // On by default, remember to SwitchOnMCData() in analysis classes\r
- reader->SwitchOnAODMCParticles(); // Off by default, remember to SwitchOnMCData() in analysis classes\r
- }\r
-\r
- //Select Trigger Class for real data\r
- if(!kMC) reader->SetFiredTriggerClassName("CINT1B-ABCE-NOPF-ALL");\r
-\r
- //Min particle pT\r
- reader->SetCTSPtMin(0.0); //new\r
- reader->SetEMCALPtMin(0.0); //new\r
- if(kFollowsFilter)reader->SetTrackStatus(0); //to prevent automatic TPC and ITS refit\r
-\r
- //In case of generating jet events (with PYTHIA), if pt hard bin is small\r
- //reject events with large difference between ptHard and triggered jet\r
- //reader->SetPtHardAndJetPtComparison(kTRUE);\r
-\r
- reader->SetFiducialCut(fidCut);\r
-\r
- if(!kInputIsESD){\r
- // Analysis with tracks, select only tracks with\r
- // following bits\r
-\r
- // //We want tracks fitted in the detectors:\r
- // ULong_t status=AliAODTrack::kTPCrefit;\r
- // status|=AliAODTrack::kITSrefit;\r
- \r
- // We want tracks whose PID bit is set:\r
- // ULong_t status =AliAODTrack::kITSpid;\r
- // status|=AliAODTrack::kTPCpid; \r
-\r
- // reader->SetTrackStatus(status);\r
- }\r
-\r
- reader->Print("");\r
-\r
-\r
- //Detector Fiducial Cuts\r
- AliFiducialCut * fidCut2 = new AliFiducialCut();\r
- fidCut2->DoEMCALFiducialCut(kTRUE) ;\r
- fidCut2->SetSimpleEMCALFiducialCut(0.7,80.,190.);\r
-\r
- fidCut2->DoCTSFiducialCut(kTRUE) ;\r
- fidCut2->SetSimpleCTSFiducialCut(0.9,0.,360.); \r
-\r
- fidCut2->Print("");\r
-\r
- //---------------------------------------------------------------------\r
- // Analysis algorithm\r
- //---------------------------------------------------------------------\r
-\r
- AliAnaElectron *anaelectron = new AliAnaElectron();\r
- anaelectron->SetDebug(-1); //10 for lots of messages\r
- anaelectron->SetCalorimeter("EMCAL");\r
- if(kMC){\r
- anaelectron->SwitchOnDataMC();\r
- anaelectron->SetMinPt(1.);\r
- }\r
- anaelectron->SetOutputAODName("Electrons");\r
- anaelectron->SetOutputAODClassName("AliAODPWG4Particle");\r
- anaelectron->SetWriteNtuple(kFALSE);\r
- //Determine which cuts to use based on enum\r
- anaelectron->SetpOverEmin(pOverEmin);\r
- anaelectron->SetpOverEmax(pOverEmax);\r
- anaelectron->SetResidualCut(dRmax);\r
- //Set Histrograms bins and ranges\r
- anaelectron->SetHistoPtRangeAndNBins(0, 100, 100) ;\r
- anaelectron->SetHistoPhiRangeAndNBins(0, TMath::TwoPi(), 100) ;\r
- anaelectron->SetHistoEtaRangeAndNBins(-0.7, 0.7, 100) ;\r
- anaelectron->Print("");\r
-\r
- //---------------------------------------------------------------------\r
- // Set analysis algorithm and reader\r
- //---------------------------------------------------------------------\r
- maker = new AliAnaPartCorrMaker();\r
- maker->SetReader(reader);//pointer to reader\r
- maker->AddAnalysis(anaelectron,0);\r
- maker->SetAnaDebug(-1) ;\r
- maker->SwitchOnHistogramsMaker() ;\r
- maker->SwitchOnAODsMaker() ;\r
-\r
- maker->Print("");\r
- //\r
- printf("============================ \n");\r
- printf("END ConfigAnalysisElectron() \n");\r
- printf("============================ \n");\r
- return maker ;\r
-}\r
+++ /dev/null
-//------------------------------------\r
-// Configuration macro:\r
-//\r
-// Configure JETAN FastJet analysis.\r
-//\r
-// Modified by: K. Read\r
-//\r
-//------------------------------------\r
-\r
-AliJetFinder* ConfigJetAnalysis()\r
-{\r
- //\r
- // Configuration goes here\r
- // \r
- printf("========================== \n");\r
- printf("ConfigJetAnalysisFastJet() \n");\r
- printf("========================== \n");\r
-\r
- Bool_t kInputIsESD = kTRUE; //uncomment for input ESD\r
- //Bool_t kInputIsESD = kFALSE; //uncomment for input AODs\r
- Bool_t kFollowsFilter = kTRUE; //uncomment if follows ESD filter task\r
- //Bool_t kFollowsFilter = kFALSE; //uncomment if no ESD filter task\r
-\r
- //Alternatively, select input via anaInputData environment variable.\r
- if (gSystem->Getenv("anaInputData")){\r
- TString kInputData = gSystem->Getenv("anaInputData");\r
- if( kInputData == "AOD" ){\r
- kInputIsESD = kFALSE;\r
- kFollowsFilter = kFALSE;\r
- }\r
- }\r
-\r
- // Define the grids\r
- AliJetGrid *grid = new AliJetGrid(419,119,0.,2*TMath::Pi(),-0.9,0.9); \r
- grid->SetGridType(1);\r
- grid->InitParams(80.*TMath::Pi()/180,190.*TMath::Pi()/180,-0.7,0.7); \r
- grid->SetMatrixIndexes();\r
- grid->SetIndexIJ();\r
- AliJetGrid *grid2 = new AliJetGrid(131,95,80.*TMath::Pi()/180.,190.*TMath::Pi()/180.,-0.7,0.7); \r
- grid2->SetGridType(0);\r
- grid2->SetMatrixIndexes();\r
- grid2->SetIndexIJ();\r
-\r
- // Define reader header\r
- if(kInputIsESD && !kFollowsFilter) AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader();\r
- else AliJetAODReaderHeader *jrh = new AliJetAODReaderHeader();\r
- jrh->SetComment("Testing");\r
- if(kInputIsESD && !kFollowsFilter) jrh->SetReadSignalOnly(kFALSE);\r
-\r
- // Detector options: 0 = Charged particles only (MomentumArray)\r
- // 1 = Charged particles only (UnitArray)\r
- // 2 = Neutral cells only (UnitArray)\r
- // 3 = Charged particles + neutral cells (UnitArray)\r
- jrh->SetDetector(3);\r
- //jrh->SetDebug(-1);\r
- //jrh->SetFiducialEta(-0.7,0.7);\r
- //jrh->SetFiducialPhi(80.*TMath::Pi()/180,190.*TMath::Pi()/180);\r
- jrh->SetPtCut(0.1);\r
- jrh->SetFiducialEta(-0.9,0.9); //fiducial range used by AliJetFillUnitArrayTracks\r
- jrh->SetFiducialPhi(0,2*TMath::Pi()); //fiducial range used by AliJetFillUnitArrayTracks\r
- \r
- // Define reader and set its header\r
- if(kInputIsESD && !kFollowsFilter) AliJetESDReader *er = new AliJetESDReader();\r
- else AliJetAODReader *er = new AliJetAODReader();\r
- er->SetReaderHeader(jrh);\r
- er->SetTPCGrid(grid);\r
- er->SetEMCalGrid(grid2);\r
- er->SetApplyMIPCorrection(kFALSE);\r
- //hadronic correction\r
- er->SetApplyFractionHadronicCorrection(kTRUE);\r
- er->SetFractionHadronicCorrection(1.0);\r
-\r
- // Define jet header\r
- AliFastJetHeaderV1 *jh=new AliFastJetHeaderV1();\r
- Double_t R=0.4;\r
- Double_t Rbkg=0.2;\r
-\r
- // AliFastJetHeaderV1 *jh=new AliFastJetHeaderV1();\r
- jh->SetComment("Fast jet code with default parameters");\r
- //jh->SetDebug(-1);\r
- //jh->SetBGMode(1); //Do BG Subtraction\r
- jh->SetBGMode(0); //No BG Subtraction. Store AOD track refs.\r
- jh->SetRparam(R); // setup parameters\r
- jh->SetRparamBkg(Rbkg); // setup parameters\r
- jh->SetPtMin(0.2);\r
- //jh->SetGhostEtaMax(0.9);\r
- jh->SetGhostArea(0.01);\r
- jh->SetGhostEtaMax(0.7);//used to set the rap_min and rap_max, that are then used by FJ\r
- jh->SetPhiRange(80.*TMath::Pi()/180+R,190.*TMath::Pi()/180-R);//used in AliFastJetFinder for the range\r
-\r
- // Define jet finder. Set its header and reader\r
- jetFinder = new AliFastJetFinder();\r
- jetFinder->SetJetHeader(jh);\r
- jetFinder->SetJetReader(er);\r
- //jetFinder->SetPlotMode(kTRUE);\r
-\r
- printf("============================== \n");\r
- printf("END ConfigJetAnalysisFastJet() \n");\r
- printf("============================== \n");\r
- \r
- return jetFinder;\r
-}\r
+++ /dev/null
-{
- gEnv->SetValue("XSec.GSI.DelegProxy","2");
- logs=TProof::Mgr("kread@localhost")->GetSessionLogs(0);
- logs->Display();
- //logs->Grep("segmentation violation");
- logs->Save("*","logs.txt");
-}
+++ /dev/null
-\r
-//\r
-// Author: G. Balbastre\r
-// Modified by K. Read\r
-//\r
-// Usage: \r
-// After previously using mymerger.C to merge each PYHIA hard bin,\r
-// do the following to download one file per pt bin:\r
-// aliensh; cd <production name>/output/merged; cp histoss* file: \r
-// Then, adjust binlist below to reflect which histograms are non-empty.\r
-// Then, locally:\r
-// root -b -x MergeFileInBins.C\r
-//\r
-// Note: This macro overwrites the downloaded histograms with the same names.\r
-// It may be appropriate to make a backup of the histograms before running\r
-// this macro.\r
-//\r
-\r
-void MergeFileInBins()\r
-{\r
- char name[128] ;\r
- TFile * mfile;\r
- TFile * sfile;\r
- TList* list;\r
- TList* newlist;\r
- TString binlist = "00 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16";\r
-\r
- if (binlist.Length()) {\r
- arr = binlist.Tokenize(" ");\r
- TObjString *os;\r
- TIter next(arr);\r
- while ((os=(TObjString*)next())) {\r
- printf("Rescaling histosscaled-merged%s.root\n",os->GetString().Data());\r
- \r
- //Rescale histograms\r
- sprintf(name,"histosscaled-merged%s.root",os->GetString().Data());\r
- mfile = new TFile(name,"read");\r
- list = (TList*) mfile->Get("histosscaled");\r
- mfile->Close();\r
- \r
- TObject * h ; \r
- Int_t split = ((TH1F*) list->FindObject("hCount"))->GetEntries();\r
- cout<<"scale with factor "<<split<< " histograms "<<list->GetEntries()<<endl;\r
- newlist = new TList();\r
- sprintf(name,"%s2",list->GetName());\r
- newlist->SetName(name);\r
-\r
- for(Int_t iter = 0; iter < list->GetEntries(); iter++){\r
- h = list->At(iter);\r
- if(h && (h->GetName()!="hCount")){\r
- if ( !strncmp(h->ClassName(),"TH",2) ) {\r
- char name[128] ; \r
- sprintf(name, "%s", h->GetName()) ; \r
- //cout<<iter<<" histo scaled : "<<name<<endl;\r
- TH1 * hout = dynamic_cast<TH1*> (h->Clone(name)) ; \r
- //if(fSumw2) hout->Sumw2();\r
- hout->Scale(1./split) ; \r
- newlist->Add(hout) ; \r
- }\r
- } \r
- }\r
-\r
- sprintf(name,"histosscaled-merged%s.root",os->GetString().Data());\r
- sfile = new TFile(name,"recreate");\r
- newlist->Write();\r
- sfile->Close(); \r
- cout<<name<<" has been recreated "<<endl;\r
- }//while\r
- }\r
-\r
-\r
- //Merge all histos\r
-\r
- cout<<"Merge all bins "<<endl;\r
-\r
- TFileMerger m;\r
- sprintf(name,"TOTALhistosscaled.root");\r
-\r
- if (binlist.Length()) {\r
- m.OutputFile(name);\r
- arr = binlist.Tokenize(" ");\r
- TObjString *os;\r
- TIter next(arr);\r
- while ((os=(TObjString*)next())) {\r
- printf("Adding bin histosscaled-merged%s.root\n",os->GetString().Data());\r
- sprintf(name,"histosscaled-merged%s.root",os->GetString().Data());\r
- m.AddFile(name);\r
- }\r
- m.Merge();\r
- }\r
-\r
-}\r
+++ /dev/null
-Documentation of these macros and goodies is located at:
-
-http://aliceinfo.cern.ch/Offline/Activities/Analysis/PWGDocumentation/PWG4/ElectronMacros.html
-
-https://twiki.cern.ch/twiki/bin/view/ALICE/PWG4JetElectron
+++ /dev/null
-void ReadAODJete(){\r
- gSystem->Load("libTree.so");\r
- gSystem->Load("libPhysics.so");\r
- gSystem->Load("libGeom.so");\r
- gSystem->Load("libVMC.so");\r
- gSystem->Load("libCGAL");\r
- gSystem->Load("libfastjet");\r
- gSystem->Load("libSISConePlugin");\r
- gSystem->Load("libSTEERBase");\r
- gSystem->Load("libESD");\r
- gSystem->Load("libAOD");\r
- gSystem->Load("libANALYSIS");\r
- gSystem->Load("libANALYSISalice");\r
- gSystem->Load("libJETAN");\r
- gSystem->Load("libFASTJETAN");\r
- gSystem->Load("libPWG4PartCorrBase");\r
- gSystem->Load("libPWG4PartCorrDep");\r
-\r
-\r
- TFile* fin=new TFile("aodoutput.root");\r
-\r
- TTree *aodTree = (TTree*)fin->Get("aodTree");\r
- AliAODEvent *ev = new AliAODEvent();\r
- ev->ReadFromTree(aodTree);\r
- Int_t nEvents = aodTree->GetEntries();\r
- cout<<nEvents<<endl;\r
-\r
- //header information\r
- AliFastJetHeaderV1* hd=( AliFastJetHeaderV1*)aodTree->GetUserInfo()->FindObject("AliJetHeader_jets");\r
- cout<<"Rparam = "<<hd->GetRparam()<<" BGmode = "<<hd->GetBGMode()<<endl;\r
- TH1F* h1=new TH1F("h1","",100,0,1000);\r
- TH1F* h2=new TH1F("h2","",100,0,1000);\r
- TH1F* h3=new TH1F("h3","",100,0,1000);\r
- TH1F* h4=new TH1F("h4","",100,0,1000);\r
-\r
- for (Int_t nEv = 0; nEv < nEvents; nEv++) {\r
-\r
- aodTree->GetEvent(nEv);\r
-\r
- //getting jets\r
- Int_t nJets = ev->GetNJets();\r
- for(Int_t iJet=0; iJet<nJets; iJet++){ \r
- jet = (AliAODJet*)ev->GetJet(iJet); \r
- cout<<jet->Pt()<<endl;\r
- cout<<"*****************************START JET PRINTOUT**************************"<<endl;\r
- jet->Print();\r
- cout<<"*****************************END JET PRINTOUT**************************"<<endl;\r
- }\r
-\r
- //getting bkg info\r
- AliAODJetEventBackground* evBkg=(AliAODJetEventBackground*)ev->FindListObject("jeteventbackground");\r
-\r
- cout<<"*******************************evBkg = "<<evBkg<<endl;\r
-\r
- if(evBkg){\r
- cout<<"*****************************START BACKGROUND PRINTOUT********************"<<endl;\r
- evBkg->Print();\r
-\r
- Float_t bkg1=evBkg->GetBackground(0);\r
- Float_t bkg2=evBkg->GetBackground(1);\r
- Float_t bkg3=evBkg->GetBackground(2);\r
- Float_t bkg4=evBkg->GetBackground(3);\r
-\r
- cout<<"ev = "<<nEv<<" bkg FastJet ="<<bkg1<<" bkg FastJet Charg="<<bkg2<<" bkg OutOfCone ="<<bkg3<<" bkg OutOfJet ="<<bkg4<<endl;\r
- cout<<"*****************************END BACKGROUND PRINTOUT********************"<<endl;\r
- h1->Fill(bkg1);\r
- h2->Fill(bkg2);\r
- h3->Fill(bkg3);\r
- h4->Fill(bkg4);\r
- }\r
-\r
- }\r
- h1->GetXaxis()->SetTitle("#rho");\r
- h2->GetXaxis()->SetTitle("#rho");\r
- h3->GetXaxis()->SetTitle("#rho");\r
- h4->GetXaxis()->SetTitle("#rho");\r
- h1->GetYaxis()->SetTitle("Entries");\r
- h2->GetYaxis()->SetTitle("Entries");\r
- h3->GetYaxis()->SetTitle("Entries");\r
- h4->GetYaxis()->SetTitle("Entries");\r
-\r
- h1->SetTitle("FastJet R=0.2 All");\r
- h2->SetTitle("FastJet R=0.4 Charged");\r
- h3->SetTitle("Out-of-cone");\r
- h4->SetTitle("Out-of-jet");\r
-\r
-\r
- TCanvas* c1=new TCanvas();\r
- c1->Divide(2,2);\r
- c1->cd(1);\r
- h1->Draw();\r
- c1->cd(2);\r
- h2->Draw();\r
- c1->cd(3);\r
- h3->Draw();\r
- c1->cd(4);\r
- h4->Draw();\r
-\r
-}\r
+++ /dev/null
-/* $Id: $ */\r
-//--------------------------------------------------\r
-// Example macro to do multi-platform JETAN/FASTJET analysis\r
-// Can be executed with Root and AliRoot\r
-//\r
-// Configured by options and definitions set in the lines below and\r
-// additional external configuration files and environement variables.\r
-//\r
-// Author: K. Read\r
-//\r
-//-------------------------------------------------\r
-enum anaModes {mLocal, mLocalCAF, mPROOF, mGRID, mPLUGIN};\r
-//mLocal = 0: Analyze locally files in your computer\r
-//mLocalCAF = 1: Analyze locally CAF files\r
-//mPROOF = 2: Analyze CAF files with PROOF\r
-//mGRID = 3: Analyze files on GRID\r
-//mPLUGIN = 4: Analyze files on GRID with AliEn plugin\r
-\r
-//---------------------------------------------------------------------------\r
-//Settings to read locally several files, only for "mLocal" mode\r
-//The different values are default, they can be set with environmental \r
-//variables: INDIR, PATTERN, NFILES, respectively\r
-//char * kInDir = "/afs/cern.ch/user/k/kread/public/data/"; \r
-char * kInDir = "/user/data/files";\r
-char * kPattern = ""; // Data are in files kInDir/kPattern+i\r
-Int_t kFile = 1; // Number of files\r
-//---------------------------------------------------------------------------\r
-//Collection file for grid analysis\r
-char * kXML = "collection.xml";\r
-//---------------------------------------------------------------------------\r
-//Data directory for PROOF analysis\r
-char * kmydataset = "/COMMON/COMMON/LHC09a4_run8101X";\r
-//char * kmydataset = "/PWG4/mcosenti/LHC08d10_ppElectronB_Jets#esdTree";\r
-//---------------------------------------------------------------------------\r
-//Scale histograms from file. Change to kTRUE when xsection file exists\r
-//Put name of file containing xsection \r
-//Put number of events per ESD file\r
-//This is an specific case for normalization of Pythia files.\r
-Bool_t kGetXSectionFromFileAndScale = kTRUE;\r
-const char * kXSFileName = "pyxsec.root";\r
-const Int_t kNumberOfEventsPerFile = 200; \r
-//---------------------------------------------------------------------------\r
-\r
-Bool_t kMC = kTRUE; //With real data kMC = kFALSE\r
-TString kInputData = "ESD";//ESD, AOD, MC\r
-TString kTreeName = "esdTree";\r
-//const Bool_t kMergeAODs = kTRUE; //uncomment for AOD merging\r
-const Bool_t kMergeAODs = kFALSE; //uncomment for no AOD merging\r
-const Bool_t kUsePAR = kTRUE; //set to kFALSE for libraries\r
-const Bool_t kDoJetTask = kTRUE; //set to kFALSE to skip JETAN task\r
-Int_t sevent = 0;\r
-\r
-Int_t mode = mLocal;\r
-char sconfig1[1024] = "ConfigPWG4AODtoAOD"; //"ConfigAnalysis";\r
-char sconfig2[1024] = "ConfigJetAnalysisFastJet.C";//"ConfigAnalysis";\r
-char sconfig3[1024] = "ConfigAnalysisElectron"; //"ConfigAnalysis";\r
-\r
-//Initialize the cross section and ntrials values. Do not modify.\r
-Double_t xsection = 0;\r
-Float_t ntrials = 0;\r
-\r
-void anaJete()\r
-{\r
- // Main\r
-\r
- //Process environmental variables from command line:\r
- ProcessEnvironment(); \r
- printf("Final Variables: kInputData %s, kMC %d, mode %d, config2 %s, config3 %s, sevent %d\n",kInputData.Data(),kMC,mode,sconfig2,sconfig3,sevent);\r
-\r
- if(!kMC) {\r
- kGetXSectionFromFileAndScale = kFALSE;\r
- }\r
-\r
- //--------------------------------------------------------------------\r
- // Load analysis libraries\r
- // Look at the method below, \r
- // change whatever you need for your analysis case\r
- // ------------------------------------------------------------------\r
- LoadLibraries(mode) ;\r
- \r
- //-------------------------------------------------------------------------------------------------\r
- //Create chain from ESD and from cross sections files, look below for options.\r
- //------------------------------------------------------------------------------------------------- \r
- if(kInputData == "ESD") kTreeName = "esdTree" ;\r
- else if(kInputData == "AOD") kTreeName = "aodTree" ;\r
- else if (kInputData == "MC") kTreeName = "TE" ;\r
- else {\r
- cout<<"Wrong data type "<<kInputData<<endl;\r
- break;\r
- }\r
-\r
- TChain * chain = new TChain(kTreeName) ;\r
- TChain * chainxs = new TChain("Xsection") ;\r
-\r
- if (mode==mLocal || mode==mLocalCAF || mode == mGRID) {\r
- CreateChain(mode, chain, chainxs); \r
- cout<<"Chain created"<<endl;\r
- }\r
-\r
- if( chain || mode==mPROOF || mode==mPLUGIN ){\r
- AliLog::SetGlobalLogLevel(AliLog::kError);//Minimum prints on screen\r
- \r
- //--------------------------------------\r
- // Make the analysis manager\r
- //-------------------------------------\r
- AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager");\r
-\r
- if( mode == mPLUGIN ){\r
- // Create and configure the alien handler plugin\r
- if (!AliAnalysisGrid::CreateToken()) return NULL;\r
- AliAnalysisAlien *plugin = new AliAnalysisAlien();\r
- plugin->SetRunMode("submit");\r
- //Uncomment the following 3 lines to permit auto xml creation\r
- //plugin->SetGridDataDir("/alice/sim/PDC_08b/LHC08d10/"); //dummy\r
- //plugin->SetDataPattern("AliESDs.root"); //dummy\r
- //plugin->AddRunNumber(30010); //dummy\r
- plugin->AddDataFile("mycollect.xml");\r
- plugin->SetGridWorkingDir("work3");\r
- plugin->SetAdditionalLibs("anaJet.C ConfigJetAnalysisFastJet.C ConfigAnalysisElectron.C ANALYSIS.par ANALYSISalice.par AOD.par EMCALUtils.par ESD.par PHOSUtils.par STEERBase.par JETAN.par FASTJETAN.par");\r
- plugin->SetJDLName("anaJet.jdl");\r
- plugin->SetExecutable("anaJet.sh");\r
- plugin->SetOutputFiles("histos.root");\r
- AliAnalysisGrid *alienHandler = plugin;\r
- if (!alienHandler) return;\r
-\r
- // Connect plug-in to the analysis manager\r
- mgr->SetGridHandler(alienHandler);\r
- }\r
-\r
- // MC handler\r
- if( (kMC && (kInputData == "ESD")) || kInputData == "MC"){\r
- AliMCEventHandler* mcHandler = new AliMCEventHandler();\r
- mcHandler->SetReadTR(kFALSE);//Do not search TrackRef file\r
- mgr->SetMCtruthEventHandler(mcHandler);\r
- if( kInputData == "MC") mgr->SetInputEventHandler(NULL);\r
- }\r
-\r
- // AOD output handler\r
- AliAODHandler* aodoutHandler = new AliAODHandler();\r
- aodoutHandler->SetOutputFileName("aodoutput.root");\r
- if(kMergeAODs)aodoutHandler->SetCreateNonStandardAOD();\r
- mgr->SetOutputEventHandler(aodoutHandler);\r
-\r
- //input\r
- if(kInputData == "ESD"){\r
- // ESD handler\r
- AliESDInputHandler *esdHandler = new AliESDInputHandler();\r
- mgr->SetInputEventHandler(esdHandler);\r
- }\r
- if(kInputData == "AOD"){\r
- // AOD handler\r
- AliAODInputHandler *aodHandler = new AliAODInputHandler();\r
- mgr->SetInputEventHandler(aodHandler);\r
- if(kMergeAODs){\r
- char path[1024];\r
- sprintf(path,"AliAOD.root");\r
- if(gSystem->Getenv("SIMPATH"))\r
- sprintf(path,"%s/AliAOD.root",gSystem->Getenv("SIMPATH"));\r
- cout<<"Config: Second input file: "<<path<<endl;\r
- aodHandler->SetMergeEvents(kTRUE);\r
- aodHandler->AddFriend(path);\r
- cout<<"Config: Starting event for second input file: "<<sevent<<endl;\r
- aodHandler->SetMergeOffset(sevent);\r
- }\r
- }\r
-\r
- mgr->SetDebugLevel(3); // For debugging, do not uncomment if you want no messages.\r
-\r
-\r
- const Bool_t kDoESDFilter = kTRUE; //need this for JETAN with input ESDs\r
- //const Bool_t kDoESDFilter = kFALSE;\r
- if(kInputData == "ESD" && kDoESDFilter){\r
- printf("Applying ESD filter cuts appropriate for jet analysis\n");\r
- //\r
- // Set of cuts\r
- // \r
- // standard\r
- AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Loose");\r
- //esdTrackCutsL->SetMinNClustersTPC(50);\r
- //esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);\r
- //esdTrackCutsL->SetRequireTPCRefit(kTRUE);\r
- //esdTrackCutsL->SetMaxDCAToVertexXY(2.4);\r
- //esdTrackCutsL->SetMaxDCAToVertexZ(3.2);\r
- //esdTrackCutsL->SetDCAToVertex2D(kTRUE);\r
- //esdTrackCutsL->SetRequireSigmaToVertex(kFALSE);\r
- //esdTrackCutsL->SetAcceptKinkDaughters(kFALSE);\r
- //\r
- // hard\r
- AliESDtrackCuts* esdTrackCutsH = new AliESDtrackCuts("AliESDtrackCuts", "Hard");\r
- esdTrackCutsH->SetMinNClustersTPC(100);\r
- esdTrackCutsH->SetMaxChi2PerClusterTPC(2.0);\r
- esdTrackCutsH->SetRequireTPCRefit(kTRUE);\r
- esdTrackCutsH->SetMaxDCAToVertexXY(2.4);\r
- esdTrackCutsH->SetMaxDCAToVertexZ(3.2);\r
- esdTrackCutsH->SetDCAToVertex2D(kTRUE);\r
- esdTrackCutsH->SetRequireSigmaToVertex(kFALSE);\r
- esdTrackCutsH->SetAcceptKinkDaughters(kFALSE);\r
- //\r
- AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");\r
- trackFilter->AddCuts(esdTrackCutsL);\r
- // trackFilter->AddCuts(esdTrackCutsH);\r
- //\r
- AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");\r
- esdfilter->SetTrackFilter(trackFilter);\r
- esdfilter->SetDebugLevel(10);\r
- mgr->AddTask(esdfilter);\r
- }\r
-\r
-\r
- //-------------------------------------------------------------------------\r
- //Define task, put here any other task that you want to use.\r
- //-------------------------------------------------------------------------\r
-\r
- //\r
- // Jet analysis\r
- //\r
- if( kDoJetTask ){\r
- AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis",chain);\r
- jetana->SetDebugLevel(2);\r
- jetana->SetConfigFile(sconfig2); //Default ConfigJetAnalysisFastJet\r
- //Uncommenting the following line produces too many AddAtAndExpand warnings for now.\r
- if(kMergeAODs)jetana->ReadAODFromOutput(); //Uncomment when AOD merging\r
- mgr->AddTask(jetana);\r
- }\r
-\r
- //\r
- // electron analysis\r
- //\r
- AliAnalysisTaskParticleCorrelation * taskpwg4 = new AliAnalysisTaskParticleCorrelation ("Particle");\r
- taskpwg4->SetConfigFileName("ConfigAnalysisElectron"); //Default name is ConfigAnalysisElectron\r
- mgr->AddTask(taskpwg4);\r
- \r
- // Create containers for input/output\r
- AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();\r
- AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();\r
- AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),\r
- AliAnalysisManager::kOutputContainer, "histos.root");\r
-\r
-\r
- if(kInputData == "ESD" && kDoESDFilter){\r
- mgr->ConnectInput (esdfilter, 0, cinput1 );\r
- mgr->ConnectOutput (esdfilter, 0, coutput1 );\r
- }\r
-\r
- if( kDoJetTask ){\r
- mgr->ConnectInput (jetana, 0, cinput1 );\r
- mgr->ConnectOutput (jetana, 0, coutput1 );\r
- mgr->ConnectOutput (jetana, 1, coutput2 );\r
- }\r
-\r
- mgr->ConnectInput (taskpwg4, 0, cinput1 );\r
- mgr->ConnectOutput (taskpwg4, 0, coutput1 );\r
- mgr->ConnectOutput (taskpwg4, 1, coutput2 );\r
- \r
-\r
- //------------------------ \r
- //Scaling task\r
- //-----------------------\r
- cout<<">>> Scaling Task"<<endl;\r
- Int_t nfiles = chain->GetListOfFiles()->GetEntriesFast();//chainxs->GetEntries();\r
- Int_t nevents = chain->GetEntries();\r
- cout<<"Get? "<<kGetXSectionFromFileAndScale<<" nfiles "<<nfiles<<" nevents "<<nevents<<endl;\r
- if(kGetXSectionFromFileAndScale && nfiles > 0){\r
- cout<<"Init AnaScale"<<endl;\r
- AliAnaScale * scale = new AliAnaScale("scale") ;\r
- cout<<"Summed xsection "<<xsection<<" Summed ntrials "<<ntrials<<" total events "<<nevents<<endl;\r
- //Calculate the average\r
- xsection/=nfiles;\r
- ntrials/=nfiles;\r
- cout<<"Average xsection "<<xsection<<" Average ntrials "<<ntrials<<" total events "<<nevents<<endl;\r
- Double_t scaleFactor = xsection/ntrials/nevents ;\r
- cout<<"Scale factor "<<scaleFactor<<endl;\r
- scale->Set(scaleFactor) ;\r
- scale->MakeSumw2(kTRUE);//If you want histograms with error bars set to kTRUE\r
- scale->SetDebugLevel(2);\r
- mgr->AddTask(scale);\r
- \r
- AliAnalysisDataContainer *coutput3 = mgr->CreateContainer("histosscaled", \r
- TList::Class(), AliAnalysisManager::kOutputContainer, "histosscaled.root");\r
- mgr->ConnectInput (scale, 0, coutput2);\r
- mgr->ConnectOutput (scale, 0, coutput3 );\r
- }\r
- \r
- //-----------------------\r
- // Run the analysis\r
- //----------------------- \r
- TString smode = "";\r
- if (mode==mLocal || mode == mLocalCAF) \r
- smode = "local";\r
- else if (mode==mPROOF) \r
- smode = "proof";\r
- else if (mode==mGRID) \r
- smode = "local";\r
- else if (mode==mPLUGIN) \r
- smode = "grid";\r
- \r
- //mgr->ResetAnalysis();\r
- mgr->InitAnalysis();\r
- mgr->PrintStatus();\r
- if (mode==mPROOF)\r
- mgr->StartAnalysis(smode.Data(),kmydataset,1500,0);\r
- else if (mode==mPLUGIN)\r
- mgr->StartAnalysis(smode.Data());\r
- else\r
- mgr->StartAnalysis(smode.Data(),chain);\r
-\r
- cout <<" Analysis ended sucessfully "<< endl ;\r
-\r
- }\r
- else cout << "Chain was not produced ! "<<endl;\r
- \r
-}\r
-\r
-void LoadLibraries(const anaModes mode) {\r
- \r
- \r
- //----------------------------------------------------------\r
- // >>>>>>>>>>> Local mode <<<<<<<<<<<<<< \r
- //----------------------------------------------------------\r
- if (mode==mLocal || mode == mLocalCAF || mode == mGRID || mode == mPLUGIN) {\r
-\r
- //--------------------------------------\r
- // Load the needed libraries most of them already loaded by aliroot\r
- //--------------------------------------\r
- gSystem->Load("libTree.so");\r
- gSystem->Load("libGeom.so");\r
- gSystem->Load("libVMC.so");\r
- gSystem->Load("libXMLIO.so");\r
- if( kDoJetTask ){\r
- gSystem->Load("libCGAL");\r
- gSystem->Load("libfastjet");\r
- gSystem->Load("libSISConePlugin");\r
- }\r
-\r
- if(kUsePAR){\r
- //--------------------------------------------------------\r
- //If you want to use root and par files from aliroot\r
- //-------------------------------------------------------- \r
- SetupPar("STEERBase");\r
- SetupPar("ESD");\r
- SetupPar("AOD");\r
- SetupPar("ANALYSIS");\r
- SetupPar("ANALYSISalice");\r
- SetupPar("PHOSUtils");\r
- SetupPar("EMCALUtils");\r
- if( kDoJetTask ){\r
- cerr<<"Now Loading JETAN"<<endl;\r
- SetupPar("JETAN");\r
- cerr<<"Done Loading JETAN"<<endl;\r
- cerr<<"Now Loading FASTJETAN"<<endl;\r
- SetupPar("FASTJETAN");\r
- cerr<<"Done Loading FASTJETAN"<<endl;\r
- }\r
- SetupPar("PWG4PartCorrBase");\r
- SetupPar("PWG4PartCorrDep");\r
- }\r
- else{\r
- //--------------------------------------------------------\r
- // If you want to use already compiled libraries \r
- // in the aliroot distribution\r
- //--------------------------------------------------------\r
- gSystem->Load("libSTEERBase");\r
- gSystem->Load("libESD");\r
- gSystem->Load("libAOD");\r
- gSystem->Load("libANALYSIS");\r
- gSystem->Load("libANALYSISalice");\r
- gSystem->Load("libPHOSUtils");\r
- gSystem->Load("libEMCALUtils");\r
- if( kDoJetTask ){\r
- gSystem->Load("libJETAN");\r
- gSystem->Load("libFASTJETAN");\r
- }\r
- gSystem->Load("libPWG4PartCorrBase");\r
- gSystem->Load("libPWG4PartCorrDep");\r
- }\r
-\r
- \r
- }\r
-\r
- //---------------------------------------------------------\r
- // <<<<<<<<<< PROOF mode >>>>>>>>>>>>\r
- //---------------------------------------------------------\r
- else if (mode==mPROOF) {\r
- //\r
- // Connect to proof\r
- // Put appropriate username here\r
- // char* myproofname = "alicecaf";\r
- char* myproofname = "kread@localhost";\r
-\r
- //TProof::Reset("proof://kread@lxb6046.cern.ch");\r
- //TProof::Reset("proof://myproofname);\r
- //TProof::Reset("myproofname",kTRUE);\r
- gEnv->SetValue("XSec.GSI.DelegProxy","2"); \r
- //TProof::Mgr(myproofname)->ShowROOTVersions();\r
- //TProof::Mgr(myproofname)->SetROOTVersion("v5-24-00");\r
- TProof::Open(myproofname);\r
-\r
- // gProof->ClearPackages();\r
- // gProof->SetLogLevel(5);\r
- // gProof->ClearPackage("STEERBase");\r
- // gProof->ClearPackage("ESD");\r
- // gProof->ClearPackage("AOD");\r
- // gProof->ClearPackage("ANALYSIS");\r
- // gProof->ClearPackage("ANALYSISalice");\r
- // gProof->ClearPackage("PHOSUtils");\r
- // gProof->ClearPackage("EMCALUtils");\r
- // if( kDoJetTask ){\r
- // gProof->ClearPackage("JETAN");\r
- // gProof->ClearPackage("FASTJETAN");\r
- // }\r
- // gProof->ClearPackage("PWG4PartCorrBase");\r
- // gProof->ClearPackage("PWG4PartCorrDep");\r
- // gProof->ShowEnabledPackages();\r
-\r
- // Enable the STEERBase Package\r
- gProof->UploadPackage("STEERBase.par");\r
- gProof->EnablePackage("STEERBase");\r
- // Enable the ESD Package\r
- gProof->UploadPackage("ESD.par");\r
- gProof->EnablePackage("ESD");\r
- // Enable the AOD Package\r
- gProof->UploadPackage("AOD.par");\r
- gProof->EnablePackage("AOD");\r
- // Enable the Analysis Package\r
- gProof->UploadPackage("ANALYSIS.par");\r
- gProof->EnablePackage("ANALYSIS");\r
- // Enable the Analysis Package\r
- gProof->UploadPackage("ANALYSISalice.par");\r
- gProof->EnablePackage("ANALYSISalice");\r
- // Enable the PHOSUtils Package\r
- gProof->UploadPackage("PHOSUtils.par");\r
- gProof->EnablePackage("PHOSUtils");\r
- // Enable the EMCALUtils Package\r
- gProof->UploadPackage("EMCALUtils.par");\r
- gProof->EnablePackage("EMCALUtils");\r
- if( kDoJetTask ){\r
- // Enable JETAN analysis\r
- gProof->UploadPackage("JETAN.par");\r
- gProof->EnablePackage("JETAN");\r
- // Enable FASTJETAN analysis\r
- gProof->UploadPackage("FASTJETAN.par");\r
- gProof->EnablePackage("FASTJETAN");\r
- }\r
- // Enable the PWG4PartCorrBase Package\r
- gProof->UploadPackage("PWG4PartCorrBase.par");\r
- gProof->EnablePackage("PWG4PartCorrBase");\r
- // Enable the PWG4PartCorrDep Package\r
- gProof->UploadPackage("PWG4PartCorrDep.par");\r
- gProof->EnablePackage("PWG4PartCorrDep");\r
-\r
- gProof->ShowEnabledPackages();\r
- } \r
- \r
-}\r
-\r
-void SetupPar(char* pararchivename)\r
-{\r
- //Load par files, create analysis libraries\r
- //For testing, if par file already decompressed and modified\r
- //classes then do not decompress.\r
- \r
- TString cdir(Form("%s", gSystem->WorkingDirectory() )) ; \r
- TString parpar(Form("%s.par", pararchivename)) ; \r
- if ( gSystem->AccessPathName(pararchivename) ) { \r
- TString processline = Form(".! tar xvzf %s",parpar.Data()) ;\r
- gROOT->ProcessLine(processline.Data());\r
- }\r
- \r
- TString ocwd = gSystem->WorkingDirectory();\r
- gSystem->ChangeDirectory(pararchivename);\r
- \r
- // check for BUILD.sh and execute\r
- if (!gSystem->AccessPathName("PROOF-INF/BUILD.sh")) {\r
- printf("*******************************\n");\r
- printf("*** Building PAR archive ***\n");\r
- cout<<pararchivename<<endl;\r
- printf("*******************************\n");\r
- \r
- if (gSystem->Exec("PROOF-INF/BUILD.sh")) {\r
- Error("runProcess","Cannot Build the PAR Archive! - Abort!");\r
- return -1;\r
- }\r
- }\r
- // check for SETUP.C and execute\r
- if (!gSystem->AccessPathName("PROOF-INF/SETUP.C")) {\r
- printf("*******************************\n");\r
- printf("*** Setup PAR archive ***\n");\r
- cout<<pararchivename<<endl;\r
- printf("*******************************\n");\r
- gROOT->Macro("PROOF-INF/SETUP.C");\r
- }\r
- \r
- gSystem->ChangeDirectory(ocwd.Data());\r
- printf("Current dir: %s\n", ocwd.Data());\r
-}\r
-\r
-\r
-\r
-void CreateChain(const anaModes mode, TChain * chain, TChain * chainxs){\r
- //Fills chain with data\r
-\r
- TString datafile="";\r
- if(kInputData == "ESD") datafile = "AliESDs.root" ;\r
- else if(kInputData == "AOD") datafile = "AliAOD.root" ;\r
- else if(kInputData == "MC") datafile = "galice.root" ;\r
-\r
- if(kInputData == "AOD") kXSFileName = "pyxsec_hists.root";\r
-\r
- char * kGener = gSystem->Getenv("GENER");\r
- if(kGener) {\r
- cout<<"GENER "<<kGener<<endl;\r
- if(!strcmp(kGener,"PYTHIA")) kXSFileName = "pyxsec.root";\r
- else if(!strcmp(kGener,"HERWIG")) kXSFileName = "hexsec.root";\r
- else cout<<" UNKNOWN GENER, use default: "<<kXSFileName<<endl;\r
- }\r
- else cout<<" GENER not set, use default xs file name: "<<kXSFileName<<endl;\r
-\r
-\r
- TString ocwd = gSystem->WorkingDirectory();\r
- \r
- //-----------------------------------------------------------\r
- //Analysis of CAF data locally\r
- //-----------------------------------------------------------\r
- if(mode == mLocalCAF){\r
- // Read the input list of files and add them to the chain\r
- TString line;\r
- ifstream in;\r
- in.open("ESDlist.txt");\r
- while (in.good()) {\r
- in >> line;\r
- if (line.Length() == 0) continue;\r
- // cout << " line = " << line << endl;\r
- chain->Add(line);\r
- }\r
- }\r
- \r
- //---------------------------------------\r
- //Local files analysis\r
- //---------------------------------------\r
- else if(mode == mLocal){\r
- //If you want to add several ESD files sitting in a common directory INDIR\r
- //Specify as environmental variables the directory (INDIR), the number of files \r
- //to analyze (NFILES) and the pattern name of the directories with files (PATTERN)\r
-\r
- if(gSystem->Getenv("INDIR")) \r
- kInDir = gSystem->Getenv("INDIR") ; \r
- else cout<<"INDIR not set, use default: "<<kInDir<<endl; \r
- \r
- if(gSystem->Getenv("PATTERN")) \r
- kPattern = gSystem->Getenv("PATTERN") ; \r
- else cout<<"PATTERN not set, use default: "<<kPattern<<endl;\r
- \r
- if(gSystem->Getenv("NFILES"))\r
- kFile = atoi(gSystem->Getenv("NFILES")) ;\r
- else cout<<"NFILES not set, use default: "<<kFile<<endl;\r
- \r
- //Check if env variables are set and are correct\r
- if ( kInDir && kFile) {\r
- printf("Get %d files from directory %s\n",kFile,kInDir);\r
- if ( ! gSystem->cd(kInDir) ) {//check if ESDs directory exist\r
- printf("%s does not exist\n", kInDir) ;\r
- return ;\r
- }\r
-\r
- //if(gSystem->Getenv("XSFILE")) \r
- //kXSFileName = gSystem->Getenv("XSFILE") ; \r
- //else cout<<" XS file name not set, use default: "<<kXSFileName<<endl; \r
-\r
- cout<<"INDIR : "<<kInDir<<endl;\r
- cout<<"NFILES : "<<kFile<<endl;\r
- cout<<"PATTERN: " <<kPattern<<endl;\r
- cout<<"XSFILE : "<<kXSFileName<<endl;\r
- \r
-\r
- //Loop on ESD files, add them to chain\r
- Int_t event =0;\r
- Int_t skipped=0 ; \r
- char file[120] ;\r
- char filexs[120] ;\r
- \r
- for (event = 0 ; event < kFile ; event++) {\r
- sprintf(file, "%s/%s%d/%s", kInDir,kPattern,event,datafile.Data()) ; \r
- sprintf(filexs, "%s/%s%d/%s", kInDir,kPattern,event,kXSFileName) ; \r
- TFile * dataFile = 0 ; \r
- //Check if file exists and add it, if not skip it\r
- if ( dataFile = TFile::Open(file)) {\r
- if ( dataFile->Get(kTreeName) ) { \r
- Int_t nEventsPerFile = ((TTree*) dataFile->Get(kTreeName)) ->GetEntries();\r
- printf("++++ Adding %s, with %d events \n", file, nEventsPerFile) ;\r
- chain->AddFile(file);\r
- if(kGetXSectionFromFileAndScale) GetXsection(nEventsPerFile, filexs); \r
- }\r
- }\r
- else { \r
- printf("---- Skipping %s\n", file) ;\r
- skipped++ ;\r
- }\r
- }\r
- printf("number of entries # %lld, skipped %d\n", chain->GetEntries(), skipped*100) ; \r
- }\r
- else {\r
- TString input = "AliESDs.root" ;\r
- cout<<">>>>>> No list added, take a single file <<<<<<<<< "<<input<<endl;\r
- chain->AddFile(input);\r
- }\r
- \r
- }// local files analysis\r
- \r
- //------------------------------\r
- //GRID xml files\r
- //-----------------------------\r
- else if(mode == mGRID){\r
- //Get colection file. It is specified by the environmental\r
- //variable XML\r
-\r
- if(gSystem->Getenv("XML") )\r
- kXML = gSystem->Getenv("XML");\r
- else\r
- sprintf(kXML, "collection.xml") ; \r
- \r
- if (!TFile::Open(kXML)) {\r
- printf("No collection file with name -- %s -- was found\n",kXML);\r
- return ;\r
- }\r
- else cout<<"XML file "<<kXML<<endl;\r
-\r
- //Load necessary libraries and connect to the GRID\r
- gSystem->Load("libNetx.so") ; \r
- gSystem->Load("libRAliEn.so"); \r
- TGrid::Connect("alien://") ;\r
-\r
- //Feed Grid with collection file\r
- //TGridCollection * collection = (TGridCollection*)gROOT->ProcessLine(Form("TAlienCollection::Open(\"%s\", 0)", kXML));\r
- TGridCollection * collection = (TGridCollection*) TAlienCollection::Open(kXML);\r
- if (! collection) {\r
- AliError(Form("%s not found", kXML)) ; \r
- return kFALSE ; \r
- }\r
- TGridResult* result = collection->GetGridResult("",0 ,0);\r
- \r
- // Makes the ESD chain \r
- printf("*** Getting the Chain ***\n");\r
- Int_t nEventsPerFile = 0;\r
- for (Int_t index = 0; index < result->GetEntries(); index++) {\r
- TString alienURL = result->GetKey(index, "turl") ; \r
- cout << "================== " << alienURL << endl ; \r
- chain->Add(alienURL) ; \r
-\r
- if(kGetXSectionFromFileAndScale){\r
- //Get the number of events per file.\r
- //Do it only once, no need to open all the files.\r
- if(index == 0) {\r
- TFile * df = TFile::Open(alienURL);\r
- nEventsPerFile = ((TTree*) df->Get(kTreeName)) ->GetEntries();\r
- df->Close();\r
- } \r
- alienURL.ReplaceAll(datafile,kXSFileName);\r
- GetXsection(nEventsPerFile, alienURL);//chainxs->Add(alienURL) ; \r
- }\r
- }\r
- }// xml analysis\r
-\r
- gSystem->ChangeDirectory(ocwd.Data());\r
-}\r
-\r
-//________________________________________________\r
-void GetXsection(Int_t nEventsPerFile, TString filexs)\r
-{\r
- // Get the cross section from the corresponding file in the directory\r
- // where the data sits.\r
- // The xsection and ntrials global variables are updated per each file.\r
- // The average of these cuantities should be calculated after.\r
-\r
- TFile *fxs = TFile::Open(filexs);\r
- if(kInputData =="AOD") { //needs improvement, in case of train with ESDs, reading output AODs for example this is wrong.\r
- TList *l = (TList*) fxs->Get("cFilterList");\r
- TH1F * hxs = l->FindObject("h1Xsec") ;\r
- TH1F * htrial = l->FindObject("h1Trials") ;\r
- if(htrial->GetEntries()!=hxs->GetEntries() || htrial->GetEntries()==0){\r
- cout<<"Careful!!! Entries in histo for cross section "<<hxs->GetEntries()<< ", for trials "<<htrial->GetEntries()<<endl;\r
- continue;\r
- } \r
- xsection += hxs->GetBinContent(1);\r
- ntrials += htrial->GetBinContent(1)/nEventsPerFile;\r
- cout << "Chain: xsection " <<hxs->GetBinContent(1)<<" ntrials "<<htrial->GetBinContent(1)<<endl; \r
- cout << "nEventsPerFile = " << nEventsPerFile <<endl; \r
- cout << "Accumulating ntrials/event = " << ntrials << endl;\r
- }\r
- else {\r
- Double_t xs = 0;\r
- UInt_t ntr = 0;\r
- TTree * xstree = (TTree*)fxs->Get("Xsection");\r
- xstree->SetBranchAddress("xsection",&xs);\r
- xstree->SetBranchAddress("ntrials",&ntr);\r
- xstree->GetEntry(0);\r
- cout << "Chain: xsection " <<xs<<" ntrials "<<ntr<<endl;\r
- xsection += xs ;\r
- ntrials += (float) ntr/nEventsPerFile;\r
- cout << "Latest values read are ntr = " << ntr << " and nEventsPerFile = " << nEventsPerFile <<endl; \r
- cout << "Accumulating ntrials/event = " << ntrials <<endl;\r
- } \r
-}\r
-\r
-void ProcessEnvironment(){\r
-\r
- if (gSystem->Getenv("anaInputData"))\r
- kInputData = gSystem->Getenv("anaInputData");\r
-\r
- if (gSystem->Getenv("anakMC"))\r
- kMC = atoi(gSystem->Getenv("anakMC"));\r
-\r
- if (gSystem->Getenv("MODE"))\r
- mode = atoi(gSystem->Getenv("MODE"));\r
-\r
- if (gSystem->Getenv("CONFIG1"))\r
- sprintf(sconfig1,gSystem->Getenv("CONFIG1"));\r
-\r
- if (gSystem->Getenv("CONFIG2"))\r
- sprintf(sconfig2,gSystem->Getenv("CONFIG2"));\r
-\r
- if (gSystem->Getenv("CONFIG3"))\r
- sprintf(sconfig3,gSystem->Getenv("CONFIG3"));\r
-\r
- if (gSystem->Getenv("SEVENT"))\r
- sevent = atoi (gSystem->Getenv("SEVENT"));\r
- \r
- printf("Process: Variables: kInputData %s, kMC %d, mode %d, config2 %s, config3 %s, sevent %d\n",kInputData.Data(),kMC,mode,sconfig2,sconfig3,sevent);\r
-\r
-}\r
+++ /dev/null
-Executable="root";
-Jobtag ={"Analysis: Processing collection for analysis "};
-Packages ={"VO_ALICE@ROOT::v5-24-00","VO_ALICE@APISCONFIG::V2.4","VO_ALICE@fastjet::v2.4.0"};
-
-TTL=15000;
-
-InputFile={"LF:/alice/cern.ch/user/k/kread/work12/anaJete.C",
- "LF:/alice/cern.ch/user/k/kread/work12/run.C",
- "LF:/alice/cern.ch/user/k/kread/work12/ConfigJetAnalysisFastJet.C",
- "LF:/alice/cern.ch/user/k/kread/work12/ConfigAnalysisElectron.C",
- "LF:/alice/cern.ch/user/k/kread/work12/STEERBase.par",
- "LF:/alice/cern.ch/user/k/kread/work12/ESD.par",
- "LF:/alice/cern.ch/user/k/kread/work12/AOD.par",
- "LF:/alice/cern.ch/user/k/kread/work12/ANALYSIS.par",
- "LF:/alice/cern.ch/user/k/kread/work12/ANALYSISalice.par",
- "LF:/alice/cern.ch/user/k/kread/work12/PHOSUtils.par",
- "LF:/alice/cern.ch/user/k/kread/work12/EMCALUtils.par",
- "LF:/alice/cern.ch/user/k/kread/work12/JETAN.par",
- "LF:/alice/cern.ch/user/k/kread/work12/FASTJETAN.par",
- "LF:/alice/cern.ch/user/k/kread/work12/PWG4PartCorrBase.par",
- "LF:/alice/cern.ch/user/k/kread/work12/PWG4PartCorrDep.par"};
-OutputArchive={"log_archive:stdout,stderr,*.log@ALICE::CERN::ALICEDISK",
- "root_archive:*.root@ALICE::CERN::ALICEDISK"};
-
-OutputDir = "/alice/cern.ch/user/k/kread/work12/output/#alien_counter_03i#";
-
-InputDataCollection = "LF:/alice/cern.ch/user/k/kread/work12/mycollect.xml,nodownload";
-
-InputDataListFormat = "xml-single";
-InputDataList = "collection.xml";
-
-Split = "se";
-SplitMaxInputFileNumber = "15";
-
-splitarguments="run.C --anaInputData ESD --mode 3 --config2 ConfigJetAnalysisFastJet.C --config3 ConfigAnalysisElectron --sevent 0";
-
-
-# Validation script to be run for each subjob
-Validationcommand = "/alice/cern.ch/user/k/kread/work12/validate.sh";
-
-OutputFile="histos.root";
-MergeOutputDir = "/alice/cern.ch/user/k/kread/work12/output/merged";
-Merge = {"histos.root:/alice/cern.ch/user/k/kread/work12/mergeout.jdl:histos-merged.root"};
+++ /dev/null
-#!/bin/bash
-
-#
-# Author: K. Read
-#
-
- export ELECUTSET=0 #kTight2 - default
-#export ELECUTSET=1 #kLooseTight - loose p/E, tight dR
-#export ELECUTSET=2 #kTightLoose - tight p/E, loose dR
-#export ELECUTSET=3 #kLoose2 - both loose
-
-# Set anaKMC=0 for data or 1 for MC
- export anakMC=1
-
- export anaInputData=ESD
- export INDIR=/work2/data/test
- export PATTERN=
- export NFILES=1
-
-#export anaInputData=AOD
-#export INDIR=/work2/data/LHC08d6/AOD/000
-#export PATTERN=01
-#export NFILES=1
-
-#export anaInputData=AOD
-#export INDIR=/work2/data/bjetfilter/AOD/117005
-#export PATTERN=00
-#export NFILES=5
-
- export MODE=0
-
- export CONFIG2=ConfigJetAnalysisFastJet.C
- export CONFIG3=ConfigAnalysisElectron
-
- export SIMPATH=/work2/data/bjetfilter/AOD/117005/004
-#export SIMPATH=/work2/data/LHC09a5/000042
-#export SIMPATH=/work2/data/LHC08d6/AOD/000/010
-
- export SEVENT=0
-
-root -q -b -l anaJete.C
+++ /dev/null
-Executable="root";
-Jobtag ={"Analysis: Processing $1 collection number $2 for analysis "};
-Packages ={"VO_ALICE@ROOT::v5-24-00","VO_ALICE@APISCONFIG::V2.4","VO_ALICE@fastjet::v2.4.0"};
-
-TTL=15000;
-
-InputFile={"LF:/alice/cern.ch/user/k/kread/$1/anaJete.C",
- "LF:/alice/cern.ch/user/k/kread/$1/run.C",
- "LF:/alice/cern.ch/user/k/kread/$1/ConfigJetAnalysisFastJet.C",
- "LF:/alice/cern.ch/user/k/kread/$1/ConfigAnalysisElectron.C",
- "LF:/alice/cern.ch/user/k/kread/$1/STEERBase.par",
- "LF:/alice/cern.ch/user/k/kread/$1/ESD.par",
- "LF:/alice/cern.ch/user/k/kread/$1/AOD.par",
- "LF:/alice/cern.ch/user/k/kread/$1/ANALYSIS.par",
- "LF:/alice/cern.ch/user/k/kread/$1/ANALYSISalice.par",
- "LF:/alice/cern.ch/user/k/kread/$1/PHOSUtils.par",
- "LF:/alice/cern.ch/user/k/kread/$1/EMCALUtils.par",
- "LF:/alice/cern.ch/user/k/kread/$1/JETAN.par",
- "LF:/alice/cern.ch/user/k/kread/$1/FASTJETAN.par",
- "LF:/alice/cern.ch/user/k/kread/$1/PWG4PartCorrBase.par",
- "LF:/alice/cern.ch/user/k/kread/$1/PWG4PartCorrDep.par"};
-OutputArchive={"log_archive:stdout,stderr,*.log@ALICE::CERN::ALICEDISK",
- "root_archive:*.root@ALICE::CERN::ALICEDISK"};
-
-OutputDir = "/alice/cern.ch/user/k/kread/$1/output/$2/#alien_counter_03i#";
-
-InputDataCollection = "LF:/alice/cern.ch/user/k/kread/$1/mycollect$2.xml,nodownload";
-
-InputDataListFormat = "xml-single";
-InputDataList = "collection.xml";
-
-Split = "se";
-SplitMaxInputFileNumber = "1";
-
-splitarguments="run.C --anaInputData AOD --mode 3 --config2 ConfigJetAnalysisFastJet.C --config3 ConfigAnalysisElectron --sevent 0";
-
-# Validation script to be run for each subjob
-Validationcommand = "/alice/cern.ch/user/k/kread/$1/validate.sh";
-
-OutputFile="histos.root";
+++ /dev/null
-/////////////////////////////////////////
-// A set of plotting utilities for the
-// electron chapter of the PPR
-//
-// J.L.Klay (Cal Poly)
-// 28-Oct-2009
-////////////////////////////////////////
-
-TH1F *alltte, *alltrk, *allemc;
-TH1F *sumtte, *sumtrk, *sumemc; //all but misid'ed hadrons
-TH1F *btte, *btrk, *bemc;
-TH1F *ctte, *ctrk, *cemc;
-TH1F *cbtte, *cbtrk, *cbemc;
-TH1F *convtte, *convtrk, *convemc;
-TH1F *daltte, *daltrk, *dalemc;
-TH1F *wztte, *wztrk, *wzemc;
-TH1F *htte, *htrk, *hemc;
-
-TH1F *allMC;
-TH1F *bMC, *cMC, *cbMC;
-TH1F *convMC, *dalMC, *wzMC;
-TH1F *mchad;
-TH1F *allheratio, *behratio;
-
-TF1* fpow;
-
-void makeData(char* jjfname = "data/scaled25Oct09/TOTALhistosscaled-LHC09b2-0.root",
- char* bfname = "data/scaled25Oct09/histosscaledLHC09b4AODc.root",
- char* wfname = "data/scaled25Oct09/histosWboson.root") {
-
- //TO AVOID DOUBLE-COUNTING ELECTRONS:
- //NOTE: Jet-Jet events are "minimum bias" which means that there
- //are B-Jets included and for the higher pThard bins the B-jets will
- //play a more important role. So to AVOID DOUBLE-COUNTING
- //B-electrons, we will suppress the electrons in the Jet-Jet events
- //that come from B->e or B->C->e when we combine the productions
- //into the total histograms
-
- //TO AVOID DOUBLE-COUNTING HADRONS:
- //For the hadron yields, we need to use only the yields from the
- //Jet-Jet events + W events and NOT USE the hadrons from the B-Jet events
-
- //For bjet and jet-jet events
- double pyscale = (1.E6)*0.5*208*208*100/360; //seconds*lumi*Pb*Pb*acceptance
- double bscale = pyscale; //Do we need to scale by Branching ratio
- //for forced semi-leptonic decays?
- double wscale = pyscale;
-
- TFile* jjfile = new TFile(jjfname);
- if(!jjfile) { printf("NO JET-JET FILE\n"); return; }
- TH2F* jjtte = (TH2F*)jjfile->Get("AnaElectron_hPtNPEleTTEScaled");
- TH2F* jjemc = (TH2F*)jjfile->Get("AnaElectron_hPtNPEleEMCALScaled");
- TH2F* jjtrk = (TH2F*)jjfile->Get("AnaElectron_hPtNPEleTPCTRDScaled");
- TH1F* jjmult = (TH1F*)jjfile->Get("AnaElectron_hRefMultScaled");
- Int_t nEvtJJ = jjmult->GetEntries();
- jjtte->Scale(pyscale);
- jjemc->Scale(pyscale);
- jjtrk->Scale(pyscale);
- TH2F* jjmcele = (TH2F*)jjfile->Get("AnaElectron_hPtMCElectronScaled");
- TH1F* jjmchad = (TH1F*)jjfile->Get("AnaElectron_hPtMCHadronScaled");
- jjmcele->Scale(pyscale);
- jjmchad->Scale(pyscale);
-
- TFile* bfile = new TFile(bfname);
- if(!bfile) { printf("NO B-JET FILE\n"); return; }
- TH2F* bjtte = (TH2F*)bfile->Get("AnaElectron_hPtNPEleTTEScaled");
- TH2F* bjemc = (TH2F*)bfile->Get("AnaElectron_hPtNPEleEMCALScaled");
- TH2F* bjtrk = (TH2F*)bfile->Get("AnaElectron_hPtNPEleTPCTRDScaled");
- TH1F* bjmult = (TH1F*)bfile->Get("AnaElectron_hRefMultScaled");
- Int_t nEvtB = bjmult->GetEntries();
- bjtte->Scale(bscale);
- bjemc->Scale(bscale);
- bjtrk->Scale(bscale);
- TH2F* bjmcele = (TH2F*)bfile->Get("AnaElectron_hPtMCElectronScaled");
- TH1F* bjmchad = (TH1F*)bfile->Get("AnaElectron_hPtMCHadronScaled");
- bjmcele->Scale(bscale);
- bjmchad->Scale(bscale);
-
- TFile* wfile = new TFile(wfname);
- if(!wfile) { printf("NO W-BOSON FILE\n"); return; }
- TH2F* wjtte = (TH2F*)wfile->Get("AnaElectron_hPtNPEleTTE");
- TH2F* wjemc = (TH2F*)wfile->Get("AnaElectron_hPtNPEleEMCAL");
- TH2F* wjtrk = (TH2F*)wfile->Get("AnaElectron_hPtNPEleTPCTRD");
- TH1F* wjmult = (TH1F*)wfile->Get("AnaElectron_hRefMult");
- Int_t nEvtW = wjmult->GetEntries();
- wjtte->Scale(wscale); //already scaled by evts
- wjemc->Scale(wscale); //already scaled by evts
- wjtrk->Scale(wscale); //already scaled by evts
- TH2F* wjmcele = (TH2F*)wfile->Get("AnaElectron_hPtMCElectron");
- TH1F* wjmchad = (TH1F*)wfile->Get("AnaElectron_hPtMCHadron");
- wjmcele->Scale(wscale);
- wjmchad->Scale(wscale);
-
- printf("Event statistics: %d (JET-JET) %d (B-JET) %d (W-Boson)\n",nEvtJJ,nEvtB,nEvtW);
-
- makeMCElectrons(jjmcele,bjmcele,wjmcele,jjmchad,wjmchad);
- makeTTEElectrons(jjtte,bjtte,wjtte);
- makeEMCElectrons(jjemc,bjemc,wjemc);
- makeTRKElectrons(jjtrk,bjtrk,wjtrk);
-
- TFile* fout = new TFile("CombinedCocktailHistograms.root","RECREATE");
- fout->cd();
- alltte->Write();
- alltrk->Write();
- allemc->Write();
-
- sumtte->Write();
- sumtrk->Write();
- sumemc->Write();
-
- btte->Write();
- btrk->Write();
- bemc->Write();
-
- ctte->Write();
- ctrk->Write();
- cemc->Write();
-
- cbtte->Write();
- cbtrk->Write();
- cbemc->Write();
-
- convtte->Write();
- convtrk->Write();
- convemc->Write();
-
- daltte->Write();
- daltrk->Write();
- dalemc->Write();
-
- wztte->Write();
- wztrk->Write();
- wzemc->Write();
-
- htte->Write();
- htrk->Write();
- hemc->Write();
-
- allMC->Write();
- bMC->Write();
- cMC->Write();
- cbMC->Write();
- convMC->Write();
- dalMC->Write();
- wzMC->Write();
- mchad->Write();
-
- allheratio->Write();
- behratio->Write();
-
- fout->Close();
-
-}
-
-
-void ScaleAndConfigure(TH1F* hist,Double_t scale, Int_t color,Bool_t keepErr)
-{
- hist->Scale(scale);
- hist->SetLineColor(color);
- hist->SetMarkerColor(color);
- hist->SetLineWidth(2);
- if(keepErr == kFALSE) {
- //remove the error bars - useful for MC rates
- for(Int_t i = 1; i <= hist->GetNbinsX(); i++) {
- if(hist->GetBinContent(i) > 0.) {
- if(hist->GetBinError(i)/hist->GetBinContent(i) > 0.5) {
- Double_t avg = 0.;
- if(i > 1 && i < hist->GetNbinsX())
- avg = (hist->GetBinContent(i-1) + hist->GetBinContent(i+1))/2.;
- hist->SetBinContent(i,avg);
- }
- }
- hist->SetBinError(i,0.);
- }
- } else {
- //Set the error bars to statistics of the bin content
- for(Int_t i = 1; i <= hist->GetNbinsX(); i++) {
- if(hist->GetBinContent(i) > 0.) {
- if(hist->GetBinError(i)/hist->GetBinContent(i) > 0.5) {
- Double_t avg = 0;
- if(i > 1 && i < hist->GetNbinsX())
- avg = (hist->GetBinContent(i-1) + hist->GetBinContent(i+1))/2.;
- hist->SetBinContent(i,avg);
- }
- }
- hist->SetBinError(i,TMath::Sqrt(hist->GetBinContent(i)));
- }
- }
-}
-
-TH1F* GetPtCutHisto(TH1F* input)
-{
- //Given a rate histogram vs pt, return the histogram with yield
- //above a given pTcut
-
- TH1F* result = (TH1F*)input->Clone();
- char name[100];
- sprintf(name,"%s_ptCut",result->GetName());
- // printf("####### %s #########\n",name);
- result->SetNameTitle(name,name);
- for(Int_t i = 1; i <= result->GetNbinsX(); i++) {
- Double_t val = input->Integral(i,result->GetNbinsX());
- //printf("<Bin %d> value %.2f integral above %.2f\n",i,input->GetBinContent(i),val);
- result->SetBinContent(i,val);
- result->SetBinError(i,0.);
- }
- //printf("################\n");
- return result;
-
-}
-
-void makeMCElectrons(TH2F* jjmcele, TH2F* bjmcele, TH2F* wjmcele, TH1F* jjmchad, TH1F* wjmchad) {
-
- //Jet-Jet MC electrons
- TH1F* jjallmc = (TH1F*)jjmcele->ProjectionX("jjallmc",1,1);
- TH1F* jjbmc = (TH1F*)jjmcele->ProjectionX("jjbmc",2,2);
- TH1F* jjcmc = (TH1F*)jjmcele->ProjectionX("jjcmc",3,3);
- TH1F* jjcandbmc = (TH1F*)jjmcele->ProjectionX("jjcandbmc",4,4);
- TH1F* jjconvmc = (TH1F*)jjmcele->ProjectionX("jjconvmc",5,5);
- TH1F* jjdalmc = (TH1F*)jjmcele->ProjectionX("jjdalmc",6,6);
-
- //Bottom-Jet MC electrons
- TH1F* bjallmc = (TH1F*)bjmcele->ProjectionX("bjallmc",1,1);
- TH1F* bjbmc = (TH1F*)bjmcele->ProjectionX("bjbmc",2,2);
- TH1F* bjcmc = (TH1F*)bjmcele->ProjectionX("bjcmc",3,3);
- TH1F* bjcandbmc = (TH1F*)bjmcele->ProjectionX("bjcandbmc",4,4);
- TH1F* bjconvmc = (TH1F*)bjmcele->ProjectionX("bjconvmc",5,5);
- TH1F* bjdalmc = (TH1F*)bjmcele->ProjectionX("bjdalmc",6,6);
-
- //W-Jet MC electrons
- TH1F* wjallmc = (TH1F*)wjmcele->ProjectionX("wjallmc",1,1);
- TH1F* wjbmc = (TH1F*)wjmcele->ProjectionX("wjbmc",2,2);
- TH1F* wjcmc = (TH1F*)wjmcele->ProjectionX("wjcmc",3,3);
- TH1F* wjcandbmc = (TH1F*)wjmcele->ProjectionX("wjcandbmc",4,4);
- TH1F* wjconvmc = (TH1F*)wjmcele->ProjectionX("wjconvmc",5,5);
- TH1F* wjdalmc = (TH1F*)wjmcele->ProjectionX("wjdalmc",6,6);
- TH1F* wjwzmc = (TH1F*)wjmcele->ProjectionX("wjwzmc",7,7);
-
- //MC Hadrons (from jj events only)
- TCanvas *ctemp = new TCanvas("ctemp");
- ctemp->Divide(2,3);
- ctemp->cd(1); gPad->SetLogy(); jjmchad->Draw();
- mchad = (TH1F*)jjmchad->Clone(); mchad->SetName("mchad");
- TH1F* temphad = (TH1F*)mchad->Clone(); temphad->SetName("temphad");
- smoothWithFit(temphad,1e5,-3,10,60,100);
- for(Int_t i = 10; i<= mchad->GetNbinsX(); i++) {
- Double_t pt = mchad->GetBinCenter(i);
- mchad->SetBinContent(i,temphad->GetFunction("fpow")->Eval(pt));
- }
- jjmchad->Draw();
- temphad->Draw("same");
- mchad->Draw("same");
- mchad->SetTitle("MC hadrons in Pb+Pb in EMCAL acceptance");
- mchad->SetName("mchad");
- mchad->SetXTitle("p_T (GeV/c)");
-
- bMC = (TH1F*)bjbmc->Clone(); bMC->SetName("bMC"); //B-Jet + W-jet
- bMC->Add(wjbmc);
- ctemp->cd(2); gPad->SetLogy(); bMC->Draw();
- TH1F* foob = (TH1F*)bMC->Clone(); foob->SetName("foob");
- smoothWithFit(foob,1e6,-3,15,50,100);
- for(Int_t i = 10; i<= bMC->GetNbinsX(); i++) bMC->SetBinContent(i,foob->GetBinContent(i));
- bMC->Draw();
- foob->Draw("same");
-
- cMC = (TH1F*)jjcmc->Clone(); cMC->SetName("cMC"); //Jet-Jet + W-jet
- cMC->Add(wjcmc);
- ctemp->cd(3); gPad->SetLogy(); cMC->Draw();
- TH1F* fooc = (TH1F*)cMC->Clone(); fooc->SetName("fooc");
- smoothWithFit(fooc,1e6,-3,5,30,100);
- for(Int_t i = 5; i<= cMC->GetNbinsX(); i++) cMC->SetBinContent(i,fooc->GetBinContent(i));
- cMC->Draw();
- fooc->Draw("same");
-
- cbMC = (TH1F*)bjcandbmc->Clone(); cbMC->SetName("cbMC"); //B-Jet + W-jet
- cbMC->Add(wjcandbmc);
- ctemp->cd(4); gPad->SetLogy(); cbMC->Draw();
- TH1F* foocb = (TH1F*)cbMC->Clone(); foocb->SetName("foocb");
- smoothWithFit(foocb,1e6,-3,8,35,100);
- for(Int_t i = 8; i<= cbMC->GetNbinsX(); i++) cbMC->SetBinContent(i,foocb->GetBinContent(i));
- cbMC->Draw();
- foocb->Draw("same");
-
- convMC = (TH1F*)jjconvmc->Clone(); convMC->SetName("convMC"); //Jet-Jet + W-jet
- convMC->Add(wjconvmc);
- ctemp->cd(5); gPad->SetLogy(); convMC->Draw();
- TH1F* fooconv = (TH1F*)convMC->Clone(); fooconv->SetName("fooconv");
- smoothWithFit(fooconv,1e6,-3,6,40,100);
- // for(Int_t i = 6; i<= convMC->GetNbinsX(); i++) convMC->SetBinContent(i,fooconv->GetBinContent(i));
- convMC->Draw();
- fooconv->Draw("same");
-
- dalMC = (TH1F*)jjdalmc->Clone(); dalMC->SetName("dalMC"); //Jet-Jet + W-jet
- dalMC->Add(wjdalmc);
- ctemp->cd(6); gPad->SetLogy(); dalMC->Draw();
- TH1F* foodal = (TH1F*)dalMC->Clone(); foodal->SetName("foodal");
- smoothWithFit(foodal,1e6,-3,3,20,100);
- for(Int_t i = 3; i<= dalMC->GetNbinsX(); i++) dalMC->SetBinContent(i,foodal->GetBinContent(i));
- dalMC->Draw();
- foodal->Draw("same");
-
- wzMC = (TH1F*)wjwzmc->Clone(); wzMC->SetName("wzMC"); //W-jet only
- TCanvas* cw= new TCanvas("cw");
- cw->cd(); gPad->SetLogy(); wzMC->GetYaxis()->SetRangeUser(1,1e3); wzMC->Draw();
- TH1F* foowz = (TH1F*)wzMC->Clone(); foowz->SetName("foowz");
- TF1* fws = new TF1("fws","[0]*(1+exp((x-[1])/[2]))^-1",39,100);
- fws->SetParameters(100,30,5);
- foowz->Fit(fws,"R");
- TF1* fwzexp = new TF1("fwzexp","[0]+[1]*log(x/[2])^2",4,40);
- fwzexp->SetParameters(100,10,3);
- foowz->Fit(fwzexp,"R");
- for(Int_t i = 8; i<= wzMC->GetNbinsX(); i++) {
- Double_t pt = wzMC->GetBinCenter(i);
- if(pt < 40) wzMC->SetBinContent(i,fwzexp->Eval(pt));
- if(pt > 40) wzMC->SetBinContent(i,fws->Eval(pt));
- }
- wzMC->GetYaxis()->SetRangeUser(1,1e3);
- wzMC->Draw();
- fws->Draw("same");
- foowz->Draw("same");
-
- //All mc electrons is the sum of
- //Jet-Jet: conversions + direct charm + dalitz
- //Bottom-Jet: direct bottom + indirect bottom
- //W-Jet: all (because these events are exclusive of the others)
- allMC = (TH1F*)wzMC->Clone(); allMC->SetNameTitle("allMC","All MC Electrons");
- allMC->Add(convMC); allMC->Add(cMC); allMC->Add(dalMC);
- allMC->Add(bMC); allMC->Add(cbMC);
-
- //Hadron/electron ratios
- allheratio = (TH1F*)allMC->Clone(); allheratio->SetName("allheratio");
- behratio = (TH1F*)bMC->Clone(); behratio->SetName("behratio");
- allheratio->SetTitle("MC hadrons and electrons in Pb+Pb, 5.5 TeV");
- allheratio->SetXTitle("p_{T} (GeV/c)");
- allheratio->SetYTitle("Hadrons/Electrons");
- for(Int_t i = 1; i <= allheratio->GetNbinsX(); i++) {
- Double_t vale = allMC->GetBinContent(i);
- Double_t valb = bMC->GetBinContent(i);
- Double_t valh = mchad->GetBinContent(i);
- //printf("pT %.2f, Hadron %.1f, Electron %.1f, B-electron
- //%.1f\n",all->GetBinCenter(i),valh,vale,valb);
- if(vale>0) allheratio->SetBinContent(i,valh/vale);
- else allheratio->SetBinContent(i,0.);
-
- if(valb>0) behratio->SetBinContent(i,valh/valb);
- else behratio->SetBinContent(i,0.);
-
- allheratio->SetBinError(i,0.);
- behratio->SetBinError(i,0.);
- }
-
- double myscale = 1.; //we already scaled them
- ScaleAndConfigure(allMC,myscale,kBlack,kFALSE);
- ScaleAndConfigure(bMC,myscale,kRed,kFALSE);
- ScaleAndConfigure(cMC,myscale,kBlue,kFALSE);
- ScaleAndConfigure(cbMC,myscale,kViolet,kFALSE);
- ScaleAndConfigure(convMC,myscale,kOrange-3,kFALSE);
- ScaleAndConfigure(dalMC,myscale,kGreen-3,kFALSE);
- ScaleAndConfigure(wzMC,myscale,kOrange-7,kFALSE);
- ScaleAndConfigure(mchad,myscale,kGreen+2,kFALSE);
-
- return;
-}
-
-void makeTTEElectrons(TH2F* jjtte, TH2F* bjtte, TH2F* wjtte) {
-
- //Jet-Jet TTE Electrons
- TH1F* jjalltte = (TH1F*)jjtte->ProjectionX("jjalltte",1,1);
- TH1F* jjbtte = (TH1F*)jjtte->ProjectionX("jjbtte",2,2);
- TH1F* jjctte = (TH1F*)jjtte->ProjectionX("jjctte",3,3);
- TH1F* jjcbtte = (TH1F*)jjtte->ProjectionX("jjcbtte",4,4);
- TH1F* jjconvtte = (TH1F*)jjtte->ProjectionX("jjconvtte",5,5);
- TH1F* jjdaltte = (TH1F*)jjtte->ProjectionX("jjdaltte",6,6);
- TH1F* jjwztte = (TH1F*)jjtte->ProjectionX("jjwztte",7,7);
- TH1F* jjhtte = (TH1F*)jjtte->ProjectionX("jjhtte",9,9);
-
- //B-Jet TTE Electrons
- TH1F* bjalltte = (TH1F*)bjtte->ProjectionX("bjalltte",1,1);
- TH1F* bjbtte = (TH1F*)bjtte->ProjectionX("bjbtte",2,2);
- TH1F* bjctte = (TH1F*)bjtte->ProjectionX("bjctte",3,3);
- TH1F* bjcbtte = (TH1F*)bjtte->ProjectionX("bjcbtte",4,4);
- TH1F* bjconvtte = (TH1F*)bjtte->ProjectionX("bjconvtte",5,5);
- TH1F* bjdaltte = (TH1F*)bjtte->ProjectionX("bjdaltte",6,6);
- TH1F* bjwztte = (TH1F*)bjtte->ProjectionX("bjwztte",7,7);
- TH1F* bjhtte = (TH1F*)bjtte->ProjectionX("bjhtte",9,9);
-
- //W-Jet TTE Electrons
- TH1F* wjalltte = (TH1F*)wjtte->ProjectionX("wjalltte",1,1);
- TH1F* wjbtte = (TH1F*)wjtte->ProjectionX("wjbtte",2,2);
- TH1F* wjctte = (TH1F*)wjtte->ProjectionX("wjctte",3,3);
- TH1F* wjcbtte = (TH1F*)wjtte->ProjectionX("wjcbtte",4,4);
- TH1F* wjconvtte = (TH1F*)wjtte->ProjectionX("wjconvtte",5,5);
- TH1F* wjdaltte = (TH1F*)wjtte->ProjectionX("wjdaltte",6,6);
- TH1F* wjwztte = (TH1F*)wjtte->ProjectionX("wjwztte",7,7);
- TH1F* wjhtte = (TH1F*)wjtte->ProjectionX("wjhtte",9,9);
-
- btte = (TH1F*)bjbtte->Clone(); btte->SetName("btte"); //B-Jet + W-jet
- btte->Add(wjbtte)
- TCanvas * ctemptte = new TCanvas("ctemptte");
- ctemptte->Divide(2,3);
- ctemptte->cd(1); gPad->SetLogy(); btte->Draw();
- TH1F* foobtte = (TH1F*)btte->Clone(); foobtte->SetName("foobtte");
- smoothWithFit(foobtte,1e5,-3,8,40,100);
- btte->Draw(); foobtte->Draw("same");
-
- ctte = (TH1F*)jjctte->Clone(); ctte->SetName("ctte"); //Jet-Jet + W-jet
- ctte->Add(wjctte);
- ctemptte->cd(2); gPad->SetLogy(); ctte->Draw();
- TH1F* fooctte = (TH1F*)ctte->Clone(); fooctte->SetName("fooctte");
- smoothWithFit(fooctte,1e5,-3,3,12,100);
- ctte->Draw(); fooctte->Draw("same");
-
- cbtte = (TH1F*)bjcbtte->Clone(); cbtte->SetName("cbtte"); //B-Jet + W-jet
- cbtte->Add(wjcbtte);
- ctemptte->cd(3); gPad->SetLogy(); cbtte->Draw();
- TH1F* foocbtte = (TH1F*)cbtte->Clone(); foocbtte->SetName("foocbtte");
- smoothWithFit(foocbtte,1e5,-3,5,20,100);
- cbtte->Draw(); foocbtte->Draw("same");
-
- convtte = (TH1F*)jjconvtte->Clone(); convtte->SetName("convtte"); //Jet-Jet + W-jet
- convtte->Add(wjconvtte);
- ctemptte->cd(4); gPad->SetLogy(); convtte->Draw();
- TH1F* fooconvtte = (TH1F*)convtte->Clone(); fooconvtte->SetName("fooconvtte");
- smoothWithFit(fooconvtte,1e5,-3,6,13,100);
- convtte->Draw(); fooconvtte->Draw("same");
-
- daltte = (TH1F*)jjdaltte->Clone(); daltte->SetName("daltte"); //Jet-Jet + W-jet
- daltte->Add(wjdaltte);
- ctemptte->cd(5); gPad->SetLogy(); daltte->Draw();
- TH1F* foodaltte = (TH1F*)daltte->Clone(); foodaltte->SetName("foodaltte");
- smoothWithFit(foodaltte,1e2,-5,2,10,100);
- daltte->Draw(); foodaltte->Draw("same");
-
- htte = (TH1F*)jjhtte->Clone(); htte->SetName("htte");
- htte->Add(wjhtte);
- ctemptte->cd(6); gPad->SetLogy(); htte->Draw();
- TH1F* foohtte = (TH1F*)htte->Clone(); foohtte->SetName("foohtte");
- smoothWithFit(foohtte,1e5,-3,10,40,100);
- htte->Draw(); foohtte->Draw("same");
-
- wztte = (TH1F*)wjwztte->Clone(); wztte->SetName("wztte"); //W-jet only
- TCanvas* ctempwztte = new TCanvas("ctempwztte");
- ctempwztte->cd(); gPad->SetLogy(); wztte->Draw();
- TH1F* foowztte = (TH1F*)wztte->Clone(); foowztte->SetName("foowztte");
- smoothWithFit(foowztte,1e5,-3,10,40,100);
- wztte->Draw(); foowztte->Draw("same");
-
- //All TTE electrons is the sum of
- //Jet-Jet: conversions + direct charm + dalitz + misid
- //Bottom-Jet: direct bottom + indirect bottom
- //W-Jet: all (because these events are exclusive of the others)
- alltte = (TH1F*)wjalltte->Clone(); alltte->SetName("alltte");
- alltte->Add(jjconvtte); alltte->Add(jjctte); alltte->Add(jjdaltte); alltte->Add(jjhtte);
- alltte->Add(bjbtte); alltte->Add(bjcbtte);
- sumtte = (TH1F*)wjalltte->Clone(); sumtte->SetName("sumtte");
- sumtte->Add(jjconvtte); sumtte->Add(jjctte); sumtte->Add(jjdaltte);
- sumtte->Add(bjbtte); sumtte->Add(bjcbtte);
-
- double myscale = 1.; //we already scaled them
- ScaleAndConfigure(alltte,myscale,kBlack,kFALSE);
- ScaleAndConfigure(sumtte,myscale,kBlack,kFALSE);
- ScaleAndConfigure(btte,myscale,kRed,kFALSE);
- ScaleAndConfigure(ctte,myscale,kBlue,kFALSE);
- ScaleAndConfigure(cbtte,myscale,kViolet,kFALSE);
- ScaleAndConfigure(convtte,myscale,kOrange-3,kFALSE);
- ScaleAndConfigure(daltte,myscale,kGreen-3,kFALSE);
- ScaleAndConfigure(htte,myscale,kGreen+2,kFALSE);
- ScaleAndConfigure(wztte,myscale,kOrange-7,kFALSE);
-
- return;
-}
-
-void makeEMCElectrons(TH2F* jjemc, TH2F* bjemc, TH2F* wjemc) {
-
- //Jet-Jet EMC Electrons
- TH1F* jjallemc = (TH1F*)jjemc->ProjectionX("jjallemc",1,1);
- TH1F* jjbemc = (TH1F*)jjemc->ProjectionX("jjbemc",2,2);
- TH1F* jjcemc = (TH1F*)jjemc->ProjectionX("jjcemc",3,3);
- TH1F* jjcbemc = (TH1F*)jjemc->ProjectionX("jjcbemc",4,4);
- TH1F* jjconvemc = (TH1F*)jjemc->ProjectionX("jjconvemc",5,5);
- TH1F* jjdalemc = (TH1F*)jjemc->ProjectionX("jjdalemc",6,6);
- TH1F* jjwzemc = (TH1F*)jjemc->ProjectionX("jjwzemc",7,7);
- TH1F* jjhemc = (TH1F*)jjemc->ProjectionX("jjhemc",9,9);
-
- //B-Jet EMC Electrons
- TH1F* bjallemc = (TH1F*)bjemc->ProjectionX("bjallemc",1,1);
- TH1F* bjbemc = (TH1F*)bjemc->ProjectionX("bjbemc",2,2);
- TH1F* bjcemc = (TH1F*)bjemc->ProjectionX("bjcemc",3,3);
- TH1F* bjcbemc = (TH1F*)bjemc->ProjectionX("bjcbemc",4,4);
- TH1F* bjconvemc = (TH1F*)bjemc->ProjectionX("bjconvemc",5,5);
- TH1F* bjdalemc = (TH1F*)bjemc->ProjectionX("bjdalemc",6,6);
- TH1F* bjwzemc = (TH1F*)bjemc->ProjectionX("bjwzemc",7,7);
- TH1F* bjhemc = (TH1F*)bjemc->ProjectionX("bjhemc",9,9);
-
- //W-Jet EMC Electrons
- TH1F* wjallemc = (TH1F*)wjemc->ProjectionX("wjallemc",1,1);
- TH1F* wjbemc = (TH1F*)wjemc->ProjectionX("wjbemc",2,2);
- TH1F* wjcemc = (TH1F*)wjemc->ProjectionX("wjcemc",3,3);
- TH1F* wjcbemc = (TH1F*)wjemc->ProjectionX("wjcbemc",4,4);
- TH1F* wjconvemc = (TH1F*)wjemc->ProjectionX("wjconvemc",5,5);
- TH1F* wjdalemc = (TH1F*)wjemc->ProjectionX("wjdalemc",6,6);
- TH1F* wjwzemc = (TH1F*)wjemc->ProjectionX("wjwzemc",7,7);
- TH1F* wjhemc = (TH1F*)wjemc->ProjectionX("wjhemc",9,9);
-
- bemc = (TH1F*)bjbemc->Clone(); bemc->SetName("bemc"); //B-Jet + W-jet
- bemc->Add(wjbemc);
- TCanvas * ctempemc = new TCanvas("ctempemc","",0,0,800,800);
- ctempemc->Divide(2,3);
- ctempemc->cd(1); gPad->SetLogy(); bemc->Draw();
- TH1F* foobemc = (TH1F*)bemc->Clone(); foobemc->SetName("foobemc");
- smoothWithFit(foobemc,1e5,-3,8,40,100);
- TRandom rand;
- for(Int_t i = 8; i<= bemc->GetNbinsX(); i++) {
- Double_t dither = rand.Gaus(0.,foobemc->GetBinContent(i)/2.);
- if(dither + foobemc->GetBinContent(i) < foobemc->GetBinContent(i)/100.) dither = 0.;
- bemc->SetBinContent(i,foobemc->GetBinContent(i)+dither);
- bemc->SetBinError(i,sqrt(foobemc->GetBinContent(i)));
- }
- bemc->Draw(); foobemc->Draw("same");
-
- cemc = (TH1F*)jjcemc->Clone(); cemc->SetName("cemc"); //Jet-Jet + W-jet
- cemc->Add(wjcemc);
- ctempemc->cd(2); gPad->SetLogy(); cemc->Draw();
- TH1F* foocemc = (TH1F*)cemc->Clone(); foocemc->SetName("foocemc");
- smoothWithFit(foocemc,1e5,-3,5,14,100);
- for(Int_t i = 5; i<= cemc->GetNbinsX(); i++) {
- Double_t dither = rand.Gaus(0.,foocemc->GetBinContent(i)/2.);
- if(dither + foocemc->GetBinContent(i) < foocemc->GetBinContent(i)/100.) dither = 0.;
- cemc->SetBinContent(i,foocemc->GetBinContent(i)+dither);
- cemc->SetBinError(i,sqrt(foocemc->GetBinContent(i)));
- }
- cemc->Draw(); foocemc->Draw("same");
-
- cbemc = (TH1F*)bjcbemc->Clone(); cbemc->SetName("cbemc"); //B-Jet + W-jet
- cbemc->Add(wjcbemc);
- ctempemc->cd(3); gPad->SetLogy(); cbemc->Draw();
- TH1F* foocbemc = (TH1F*)cbemc->Clone(); foocbemc->SetName("foocbemc");
- smoothWithFit(foocbemc,1e5,-3,8,20,100);
- for(Int_t i = 8; i<= cbemc->GetNbinsX(); i++) {
- Double_t dither = rand.Gaus(0.,foocbemc->GetBinContent(i)/2.);
- if(dither + foocbemc->GetBinContent(i) < foocbemc->GetBinContent(i)/100.) dither = 0.;
- cbemc->SetBinContent(i,foocbemc->GetBinContent(i)+dither);
- cbemc->SetBinError(i,sqrt(foocbemc->GetBinContent(i)));
- }
- cbemc->Draw(); foocbemc->Draw("same");
-
- convemc = (TH1F*)jjconvemc->Clone(); convemc->SetName("convemc"); //Jet-Jet + W-jet
- convemc->Add(wjconvemc);
- ctempemc->cd(4); gPad->SetLogy(); convemc->Draw();
- TH1F* fooconvemc = (TH1F*)convemc->Clone(); fooconvemc->SetName("fooconvemc");
- smoothWithFit(fooconvemc,1e6,-3,5,15,100);
- for(Int_t i = 5; i<= convemc->GetNbinsX(); i++) {
- Double_t dither = rand.Gaus(0.,fooconvemc->GetBinContent(i)/2.);
- if(dither + fooconvemc->GetBinContent(i) < fooconvemc->GetBinContent(i)/100.) dither = 0.;
- convemc->SetBinContent(i,fooconvemc->GetBinContent(i)+dither);
- convemc->SetBinError(i,sqrt(fooconvemc->GetBinContent(i)));
- }
- convemc->Draw(); fooconvemc->Draw("same");
-
- dalemc = (TH1F*)jjdalemc->Clone(); dalemc->SetName("dalemc"); //Jet-Jet + W-jet
- dalemc->Add(wjdalemc);
- ctempemc->cd(5); gPad->SetLogy(); dalemc->Draw();
- TH1F* foodalemc = (TH1F*)dalemc->Clone(); foodalemc->SetName("foodalemc");
- for(Int_t i = 18; i <= dalemc->GetNbinsX(); i++) {
- dalemc->SetBinContent(i,0);
- dalemc->SetBinError(i,0);
- foodalemc->SetBinContent(i,0);
- foodalemc->SetBinError(i,0);
- }
- // smoothWithFit(foodalemc,1e6,-3,8.,12.,100);
- dalemc->Draw(); //foodalemc->Draw("same");
-
- hemc = (TH1F*)jjhemc->Clone(); hemc->SetName("hemc");
- hemc->Add(wjhemc);
- ctempemc->cd(6); gPad->SetLogy(); hemc->Draw();
- TH1F* foohemc = (TH1F*)hemc->Clone(); foohemc->SetName("foohemc");
- smoothWithFit(foohemc,1e5,-3,10,40,100);
- for(Int_t i = 10; i<= hemc->GetNbinsX(); i++) {
- Double_t dither = rand.Gaus(0.,foohemc->GetBinContent(i)/2.);
- if(dither + foohemc->GetBinContent(i) < foohemc->GetBinContent(i)/100.) dither = 0.;
- hemc->SetBinContent(i,foohemc->GetBinContent(i)+dither);
- hemc->SetBinError(i,sqrt(foohemc->GetBinContent(i)));
- }
- hemc->Draw(); foohemc->Draw("same");
-
- wzemc = (TH1F*)wjwzemc->Clone(); wzemc->SetName("wzemc"); //W-jet only
- TCanvas* ctempwzemc = new TCanvas("ctempwzemc");
- ctempwzemc->cd(); gPad->SetLogy(); wzemc->Draw();
- TH1F* foowzemc = (TH1F*)wzemc->Clone(); foowzemc->SetName("foowzemc");
- TF1* fwsemc = new TF1("fwsemc","[0]*(1+exp((x-[1])/[2]))^-1",30,50);
- fwsemc->SetParameters(10,30,3);
- foowzemc->Fit(fwsemc,"R");
- TF1* fwzexpemc = new TF1("fwzexpemc","[0]+[1]*log(x/[2])^2",5,20);
- fwzexpemc->SetParameters(10,10,3);
- // foowzemc->Fit(fwzexpemc,"R");
- /*
- for(Int_t i = 8; i<= wzMC->GetNbinsX(); i++) {
- Double_t pt = wzMC->GetBinCenter(i);
- if(pt < 40) wzMC->SetBinContent(i,fwzexp->Eval(pt));
- if(pt > 40) wzMC->SetBinContent(i,fws->Eval(pt));
- }
- */
- wzemc->Draw(); foowzemc->Draw("same");
- fwzexpemc->Draw("same");
-
- //All EMC electrons is the sum of
- //Jet-Jet: conversions + direct charm + dalitz + misid
- //Bottom-Jet: direct bottom + indirect bottom
- //W-Jet: all (because these events are exclusive of the others)
- allemc = (TH1F*)wjallemc->Clone(); allemc->SetName("allemc");
- allemc->Add(jjconvemc); allemc->Add(jjcemc); allemc->Add(jjdalemc); allemc->Add(jjhemc);
- allemc->Add(bjbemc); allemc->Add(bjcbemc);
- sumemc = (TH1F*)wjallemc->Clone(); sumemc->SetName("sumemc");
- sumemc->Add(jjconvemc); sumemc->Add(jjcemc); sumemc->Add(jjdalemc);
- sumemc->Add(bjbemc); sumemc->Add(bjcbemc);
-
- double myscale = 1.; //we already scaled them
- ScaleAndConfigure(allemc,myscale,kBlack,kFALSE);
- ScaleAndConfigure(sumemc,myscale,kBlack,kFALSE);
- ScaleAndConfigure(bemc,myscale,kRed,kFALSE);
- ScaleAndConfigure(cemc,myscale,kBlue,kFALSE);
- ScaleAndConfigure(cbemc,myscale,kViolet,kFALSE);
- ScaleAndConfigure(convemc,myscale,kOrange-3,kFALSE);
- ScaleAndConfigure(dalemc,myscale,kGreen-3,kFALSE);
- ScaleAndConfigure(hemc,myscale,kGreen+2,kFALSE);
- ScaleAndConfigure(wzemc,myscale,kOrange-7,kFALSE);
-
- return;
-}
-
-void makeTRKElectrons(TH2F* jjtrk, TH2F* bjtrk, TH2F* wjtrk) {
-
- //Jet-Jet TRK Electrons
- TH1F* jjalltrk = (TH1F*)jjtrk->ProjectionX("jjalltrk",1,1);
- TH1F* jjbtrk = (TH1F*)jjtrk->ProjectionX("jjbtrk",2,2);
- TH1F* jjctrk = (TH1F*)jjtrk->ProjectionX("jjctrk",3,3);
- TH1F* jjcbtrk = (TH1F*)jjtrk->ProjectionX("jjcbtrk",4,4);
- TH1F* jjconvtrk = (TH1F*)jjtrk->ProjectionX("jjconvtrk",5,5);
- TH1F* jjdaltrk = (TH1F*)jjtrk->ProjectionX("jjdaltrk",6,6);
- TH1F* jjwztrk = (TH1F*)jjtrk->ProjectionX("jjwztrk",7,7);
- TH1F* jjhtrk = (TH1F*)jjtrk->ProjectionX("jjhtrk",9,9);
-
- //B-Jet TRK Electrons
- TH1F* bjalltrk = (TH1F*)bjtrk->ProjectionX("bjalltrk",1,1);
- TH1F* bjbtrk = (TH1F*)bjtrk->ProjectionX("bjbtrk",2,2);
- TH1F* bjctrk = (TH1F*)bjtrk->ProjectionX("bjctrk",3,3);
- TH1F* bjcbtrk = (TH1F*)bjtrk->ProjectionX("bjcbtrk",4,4);
- TH1F* bjconvtrk = (TH1F*)bjtrk->ProjectionX("bjconvtrk",5,5);
- TH1F* bjdaltrk = (TH1F*)bjtrk->ProjectionX("bjdaltrk",6,6);
- TH1F* bjwztrk = (TH1F*)bjtrk->ProjectionX("bjwztrk",7,7);
- TH1F* bjhtrk = (TH1F*)bjtrk->ProjectionX("bjhtrk",9,9);
-
- //W-Jet TRK Electrons
- TH1F* wjalltrk = (TH1F*)wjtrk->ProjectionX("wjalltrk",1,1);
- TH1F* wjbtrk = (TH1F*)wjtrk->ProjectionX("wjbtrk",2,2);
- TH1F* wjctrk = (TH1F*)wjtrk->ProjectionX("wjctrk",3,3);
- TH1F* wjcbtrk = (TH1F*)wjtrk->ProjectionX("wjcbtrk",4,4);
- TH1F* wjconvtrk = (TH1F*)wjtrk->ProjectionX("wjconvtrk",5,5);
- TH1F* wjdaltrk = (TH1F*)wjtrk->ProjectionX("wjdaltrk",6,6);
- TH1F* wjwztrk = (TH1F*)wjtrk->ProjectionX("wjwztrk",7,7);
- TH1F* wjhtrk = (TH1F*)wjtrk->ProjectionX("wjhtrk",9,9);
-
- btrk = (TH1F*)bjbtrk->Clone(); btrk->SetName("btrk"); //B-Jet + W-jet
- btrk->Add(wjbtrk);
- TCanvas * ctemptrk = new TCanvas("ctemptrk","",0,0,800,800);
- ctemptrk->Divide(2,3);
- ctemptrk->cd(1); gPad->SetLogy(); btrk->Draw();
- TH1F* foobtrk = (TH1F*)btrk->Clone(); foobtrk->SetName("foobtrk");
- smoothWithFit(foobtrk,1e5,-3,8,40,100);
- TRandom rand2;
- for(Int_t i = 8; i<= btrk->GetNbinsX(); i++) {
- Double_t dither = rand2.Gaus(0.,foobtrk->GetBinContent(i)/2.);
- if(dither + foobtrk->GetBinContent(i) < foobtrk->GetBinContent(i)/100.) dither = 0.;
- btrk->SetBinContent(i,foobtrk->GetBinContent(i)+dither);
- btrk->SetBinError(i,sqrt(foobtrk->GetBinContent(i)));
- }
- btrk->Draw(); foobtrk->Draw("same");
-
- ctrk = (TH1F*)jjctrk->Clone(); ctrk->SetName("ctrk"); //Jet-Jet + W-jet
- ctrk->Add(wjctrk);
- ctemptrk->cd(2); gPad->SetLogy(); ctrk->Draw();
- TH1F* fooctrk = (TH1F*)ctrk->Clone(); fooctrk->SetName("fooctrk");
- smoothWithFit(fooctrk,1e5,-3,5,20,100);
- for(Int_t i = 5; i<= ctrk->GetNbinsX(); i++) {
- Double_t dither = rand2.Gaus(0.,fooctrk->GetBinContent(i)/2.);
- if(dither + fooctrk->GetBinContent(i) < fooctrk->GetBinContent(i)/100.) dither = 0.;
- ctrk->SetBinContent(i,fooctrk->GetBinContent(i)+dither);
- ctrk->SetBinError(i,sqrt(fooctrk->GetBinContent(i)));
- }
- ctrk->Draw(); fooctrk->Draw("same");
-
- cbtrk = (TH1F*)bjcbtrk->Clone(); cbtrk->SetName("cbtrk"); //B-Jet + W-jet
- cbtrk->Add(wjcbtrk);
- ctemptrk->cd(3); gPad->SetLogy(); cbtrk->Draw();
- TH1F* foocbtrk = (TH1F*)cbtrk->Clone(); foocbtrk->SetName("foocbtrk");
- smoothWithFit(foocbtrk,1e5,-3,5,20,100);
- for(Int_t i = 5; i<= cbtrk->GetNbinsX(); i++) {
- Double_t dither = rand2.Gaus(0.,foocbtrk->GetBinContent(i)/2.);
- if(dither + foocbtrk->GetBinContent(i) < foocbtrk->GetBinContent(i)/100.) dither = 0.;
- cbtrk->SetBinContent(i,foocbtrk->GetBinContent(i)+dither);
- cbtrk->SetBinError(i,sqrt(foocbtrk->GetBinContent(i)));
- }
- cbtrk->Draw(); foocbtrk->Draw("same");
-
- convtrk = (TH1F*)jjconvtrk->Clone(); convtrk->SetName("convtrk"); //Jet-Jet + W-jet
- convtrk->Add(wjconvtrk);
- ctemptrk->cd(4); gPad->SetLogy(); convtrk->Draw();
- TH1F* fooconvtrk = (TH1F*)convtrk->Clone(); fooconvtrk->SetName("fooconvtrk");
- smoothWithFit(fooconvtrk,1e6,-3,5,20,100);
- for(Int_t i = 5; i<= convtrk->GetNbinsX(); i++) {
- Double_t dither = rand2.Gaus(0.,fooconvtrk->GetBinContent(i)/2.);
- if(dither + fooconvtrk->GetBinContent(i) < fooconvtrk->GetBinContent(i)/100.) dither = 0.;
- convtrk->SetBinContent(i,fooconvtrk->GetBinContent(i)+dither);
- convtrk->SetBinError(i,sqrt(fooconvtrk->GetBinContent(i)));
- }
- convtrk->Draw(); fooconvtrk->Draw("same");
-
- daltrk = (TH1F*)jjdaltrk->Clone(); daltrk->SetName("daltrk"); //Jet-Jet + W-jet
- daltrk->Add(wjdaltrk);
- ctemptrk->cd(5); gPad->SetLogy(); daltrk->Draw();
- TH1F* foodaltrk = (TH1F*)daltrk->Clone(); foodaltrk->SetName("foodaltrk");
- smoothWithFit(foodaltrk,1e6,-3,5.,40.,100);
- for(Int_t i = 5; i<= daltrk->GetNbinsX(); i++) {
- Double_t dither = rand2.Gaus(0.,foodaltrk->GetBinContent(i)/2.);
- if(dither + foodaltrk->GetBinContent(i) < foodaltrk->GetBinContent(i)/100.) dither = 0.;
- daltrk->SetBinContent(i,foodaltrk->GetBinContent(i)+dither);
- daltrk->SetBinError(i,sqrt(foodaltrk->GetBinContent(i)));
- }
- daltrk->Draw(); foodaltrk->Draw("same");
-
- htrk = (TH1F*)jjhtrk->Clone(); htrk->SetName("htrk");
- htrk->Add(wjhtrk);
- ctemptrk->cd(6); gPad->SetLogy(); htrk->Draw();
- TH1F* foohtrk = (TH1F*)htrk->Clone(); foohtrk->SetName("foohtrk");
- smoothWithFit(foohtrk,1e5,-3,5,80,100);
- for(Int_t i = 5; i<= htrk->GetNbinsX(); i++) {
- Double_t dither = rand2.Gaus(0.,foohtrk->GetBinContent(i)/2.);
- if(dither + foohtrk->GetBinContent(i) < foohtrk->GetBinContent(i)/100.) dither = 0.;
- htrk->SetBinContent(i,foohtrk->GetBinContent(i)+dither);
- htrk->SetBinError(i,sqrt(foohtrk->GetBinContent(i)));
- }
- htrk->Draw(); foohtrk->Draw("same");
-
- wztrk = (TH1F*)wjwztrk->Clone(); wztrk->SetName("wztrk"); //W-jet only
- TCanvas* ctempwztrk = new TCanvas("ctempwztrk");
- ctempwztrk->cd(); gPad->SetLogy(); wztrk->Draw();
- TH1F* foowztrk = (TH1F*)wztrk->Clone(); foowztrk->SetName("foowztrk");
- TF1* fwstrk = new TF1("fwstrk","[0]*(1+exp((x-[1])/[2]))^-1",30,50);
- fwstrk->SetParameters(10,30,3);
- foowztrk->Fit(fwstrk,"R");
- TF1* fwzexptrk = new TF1("fwzexptrk","[0]+[1]*log(x/[2])^2",8,40);
- fwzexptrk->SetParameters(10,10,3);
- // foowztrk->Fit(fwzexptrk,"R");
- /*
- for(Int_t i = 8; i<= wzMC->GetNbinsX(); i++) {
- Double_t pt = wzMC->GetBinCenter(i);
- if(pt < 40) wzMC->SetBinContent(i,fwzexp->Eval(pt));
- if(pt > 40) wzMC->SetBinContent(i,fws->Eval(pt));
- }
- */
- wztrk->Draw(); foowztrk->Draw("same");
-
- //All TRK electrons is the sum of
- //Jet-Jet: conversions + direct charm + dalitz + misid
- //Bottom-Jet: direct bottom + indirect bottom
- //W-Jet: all (because these events are exclusive of the others)
- alltrk = (TH1F*)wjalltrk->Clone(); alltrk->SetName("alltrk");
- alltrk->Add(jjconvtrk); alltrk->Add(jjctrk); alltrk->Add(jjdaltrk); alltrk->Add(jjhtrk);
- alltrk->Add(bjbtrk); alltrk->Add(bjcbtrk);
- sumtrk = (TH1F*)wjalltrk->Clone(); sumtrk->SetName("sumtrk");
- sumtrk->Add(jjconvtrk); sumtrk->Add(jjctrk); sumtrk->Add(jjdaltrk);
- sumtrk->Add(bjbtrk); sumtrk->Add(bjcbtrk);
-
- double myscale = 1.; //we already scaled them
- ScaleAndConfigure(alltrk,myscale,kBlack,kFALSE);
- ScaleAndConfigure(sumtrk,myscale,kBlack,kFALSE);
- ScaleAndConfigure(btrk,myscale,kRed,kFALSE);
- ScaleAndConfigure(ctrk,myscale,kBlue,kFALSE);
- ScaleAndConfigure(cbtrk,myscale,kViolet,kFALSE);
- ScaleAndConfigure(convtrk,myscale,kOrange-3,kFALSE);
- ScaleAndConfigure(daltrk,myscale,kGreen-3,kFALSE);
- ScaleAndConfigure(htrk,myscale,kGreen+2,kFALSE);
- ScaleAndConfigure(wztrk,myscale,kOrange-7,kFALSE);
-
- return;
-}
-
-void smoothWithFit(TH1F* hist, Double_t p0, Double_t p1, Double_t min, Double_t max,Double_t remax) {
-
- fpow = new TF1("fpow","[0]*pow(x,[1])",min,max);
- fpow->SetParameters(p0,p1);
- hist->Fit(fpow,"R");
- for(Int_t i = (Int_t)min; i <= (Int_t)remax; i++) {
- Double_t pt = hist->GetBinCenter(i);
- Double_t val = fpow->Eval(pt);
- hist->SetBinContent(i,val);
- hist->SetBinError(i,sqrt(val));
- }
-
- return;
-}
-
-
+++ /dev/null
-
-
-
-void makeSysErr() {
-
- TFile* f0 = new TFile("data/scaled25Oct09/TOTALhistosscaled-LHC09b2-0.root");
- TFile* f1 = new TFile("data/scaled25Oct09/TOTALhistosscaled-LHC09b2-1.root");
- TFile* f2 = new TFile("data/scaled25Oct09/TOTALhistosscaled-LHC09b2-2.root");
- TFile* f3 = new TFile("data/scaled25Oct09/TOTALhistosscaled-LHC09b2-3.root");
-
- double pyscale = (1.E6)*0.5*208*208*100/360; //seconds*lumi*Pb*Pb*acceptance
-
- TH2F* tte0 = (TH2F*)f0->Get("AnaElectron_hPtNPEleTTEScaled");
- TH2F* emc0 = (TH2F*)f0->Get("AnaElectron_hPtNPEleEMCALScaled");
-
- TH2F* tte1 = (TH2F*)f1->Get("AnaElectron_hPtNPEleTTEScaled");
- TH2F* emc1 = (TH2F*)f1->Get("AnaElectron_hPtNPEleEMCALScaled");
-
- TH2F* tte2 = (TH2F*)f2->Get("AnaElectron_hPtNPEleTTEScaled");
- TH2F* emc2 = (TH2F*)f2->Get("AnaElectron_hPtNPEleEMCALScaled");
-
- TH2F* tte3 = (TH2F*)f3->Get("AnaElectron_hPtNPEleTTEScaled");
- TH2F* emc3 = (TH2F*)f3->Get("AnaElectron_hPtNPEleEMCALScaled");
-
- tte0->Scale(pyscale);
- tte1->Scale(pyscale);
- tte2->Scale(pyscale);
- tte3->Scale(pyscale);
- emc0->Scale(pyscale);
- emc1->Scale(pyscale);
- emc2->Scale(pyscale);
- emc3->Scale(pyscale);
-
- alltte0 = (TH1F*)tte0->ProjectionX("alltte0",1,1);
- allemc0 = (TH1F*)emc0->ProjectionX("allemc0",1,1);
- alltte0->Rebin(5);
- allemc0->Rebin(5);
-
- alltte1 = (TH1F*)tte1->ProjectionX("alltte1",1,1);
- allemc1 = (TH1F*)emc1->ProjectionX("allemc1",1,1);
- alltte1->Rebin(5);
- allemc1->Rebin(5);
-
- alltte2 = (TH1F*)tte2->ProjectionX("alltte2",1,1);
- allemc2 = (TH1F*)emc2->ProjectionX("allemc2",1,1);
- alltte2->Rebin(5);
- allemc2->Rebin(5);
-
- alltte3 = (TH1F*)tte3->ProjectionX("alltte3",1,1);
- allemc3 = (TH1F*)emc3->ProjectionX("allemc3",1,1);
- alltte3->Rebin(5);
- allemc3->Rebin(5);
-
- Double_t tsum = 0.;
- Double_t esum = 0.;
- for(Int_t i = 1; i <= alltte0->GetNbinsX(); i++) {
- Double_t t0 = alltte0->GetBinContent(i);
- Double_t t1 = alltte1->GetBinContent(i);
- Double_t t2 = alltte2->GetBinContent(i);
- Double_t t3 = alltte3->GetBinContent(i);
- Double_t e0 = allemc0->GetBinContent(i);
- Double_t e1 = allemc1->GetBinContent(i);
- Double_t e2 = allemc2->GetBinContent(i);
- Double_t e3 = allemc3->GetBinContent(i);
-
- Double_t td01 = (t0 - t1)/(t0+t1+t2+t3);
- Double_t td02 = (t0 - t2)/(t0+t1+t2+t3);
- Double_t td03 = (t0 - t3)/(t0+t1+t2+t3);
-
- Double_t tave = (TMath::Abs(td01) + TMath::Abs(td02) + TMath::Abs(td03))/2.;
- tsum += tave;
-
- Double_t ed01 = (e0 - e1)/(e0+e1+e2+e3);
- Double_t ed02 = (e0 - e2)/(e0+e1+e2+e3);
- Double_t ed03 = (e0 - e3)/(e0+e1+e2+e3);
-
- Double_t eave = (TMath::Abs(ed01) + TMath::Abs(ed02) + TMath::Abs(ed03))/2.;
- printf("%d Average unc = %2.2f\n",i,eave);
- esum += eave;
- }
- Double_t tfinal = tsum/alltte0->GetNbinsX();
- Double_t efinal = esum/allemc0->GetNbinsX();
- printf("tfinal %f, efinal %f\n",tfinal,efinal);
-
- TCanvas *c1 = new TCanvas();
- gPad->SetLogy();
- alltte0->SetLineWidth(2);
- alltte0->Draw();
-
- TCanvas *c2 = new TCanvas();
- gPad->SetLogy();
- allemc0->SetLineWidth(2);
- allemc0->Draw();
-
- TGraphErrors* terr = new TGraphErrors();
- terr->SetName("tteErr");
- TGraphErrors* eerr = new TGraphErrors();
- eerr->SetName("emcErr");
- for(Int_t i = 1; i <= alltte0->GetNbinsX(); i++) {
- terr->SetPoint(i-1,alltte0->GetBinCenter(i),alltte0->GetBinContent(i));
- terr->SetPointError(i-1,0.,(tfinal)*alltte0->GetBinContent(i));
-
- eerr->SetPoint(i-1,allemc0->GetBinCenter(i),allemc0->GetBinContent(i));
- eerr->SetPointError(i-1,0.,(efinal)*allemc0->GetBinContent(i));
- }
- c1->cd();
- terr->SetFillColor(kRed-8);
- terr->Draw("3same");
- alltte0->Draw("same");
-
- c2->cd();
- eerr->SetFillColor(kRed-8);
- eerr->Draw("3same");
- allemc0->Draw("same");
-}
+++ /dev/null
-#! /bin/bash
-
-echo "Making parfiles..."
-
-source /home/aron/root/setup.sh
-
-/home/aron/root/alice/AliRoot/PWG4/macros/electrons/./parmaker all remote
-
-tar xvf EMCALUtils.par
-
-echo "...copying"
-#cp /home/aron/root/alice/AliRoot/EMCAL/AliCaloPeakFinderConstants.h EMCALUtils/.
-cp -v /home/aron/root/alice/AliRoot/EMCAL/AliCaloConstants.h EMCALUtils/.
-
-echo "...removing"
-rm -v EMCALUtils.par
-
-echo "...remaking EMCALUtils.par"
-/home/aron/root/alice/AliRoot/PWG4/macros/electrons/./parmaker EMCALUtils
-
-echo "...removing EMCALUtils dir"
-rm -rfv EMCALUtils
-
-echo "...making directory GRID and copying important stuff there"
-mkdir GRID
-
-cp -v *.par parmaker anaJete.C run.C ConfigAnalysisElectron.C ConfigJetAnalysisFastJet.C anaJete.jdl validate.sh mergeout.jdl GRID/.
-
+++ /dev/null
-Packages ={"VO_ALICE@ROOT::v5-24-00","VO_ALICE@APISCONFIG::V2.4"};
-
-Arguments = "27498605 histos.root histos-merged.root kread /alice/cern.ch/user/k/kread/work12/output";
-Executable = "/alice/bin/mergerootfile";
-Email = "";
-OutputFile = "histos-merged.root@ALICE::CERN::ALICEDISK";
-InputFile = "LF:/alice/macros/mergerootfile.C";
-OutputDir = "/alice/cern.ch/user/k/kread/work12/output/merged";
+++ /dev/null
-Packages ={"VO_ALICE@ROOT::v5-24-00","VO_ALICE@APISCONFIG::V2.4"};
-
-Arguments = "27498605 histosscaled.root histosscaled-merged.root kread /alice/cern.ch/user/k/kread/work12/output";
-Executable = "/alice/bin/mergerootfile";
-Email = "";
-OutputFile = "histosscaled-merged.root@ALICE::CERN::ALICEDISK";
-InputFile = "LF:/alice/macros/mergerootfile.C";
-OutputDir = "/alice/cern.ch/user/k/kread/work12/output/merged";
+++ /dev/null
-Packages ={"VO_ALICE@ROOT::v5-24-00","VO_ALICE@APISCONFIG::V2.4"};
-
-Arguments = "27498605 histosscaled.root histosscaled-merged$2.root kread /alice/cern.ch/user/k/kread/$1/output/$2";
-Executable = "/alice/bin/mergerootfile";
-Email = "";
-OutputFile = "histosscaled-merged$2.root@ALICE::CERN::ALICEDISK";
-InputFile = "LF:/alice/macros/mergerootfile.C";
-OutputDir = "/alice/cern.ch/user/k/kread/$1/output/merged";
+++ /dev/null
-<?xml version="1.0"?>
-<alien>
- <collection name="mycollect">
- <event name="1">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40671374" entryId="40671390" expiretime="" gowner="aliprod" guid="B597B85E-4FF3-11DE-8001-91B0F682BEEF" guidtime="11DE4FF3" lfn="/alice/sim/PDC_08b/LHC08d10/50036/001/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13333579" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/001/AliESDs.root" type="f" />
- </event>
- <event name="10">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40747801" entryId="40747817" expiretime="" gowner="aliprod" guid="A68AC094-4FF4-11DE-8001-810A010ABEEF" guidtime="11DE4FF4" lfn="/alice/sim/PDC_08b/LHC08d10/50036/010/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13486614" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/010/AliESDs.root" type="f" />
- </event>
- <event name="11">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40775061" entryId="40775072" expiretime="" gowner="aliprod" guid="0307DB22-4FF5-11DE-9717-6C5530C1BEEF" guidtime="11DE4FF5" lfn="/alice/sim/PDC_08b/LHC08d10/50036/011/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12813014" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/011/AliESDs.root" type="f" />
- </event>
- <event name="12">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40691174" entryId="40691190" expiretime="" gowner="aliprod" guid="7F8DE4E0-500C-11DE-9717-F1B28E80BEEF" guidtime="11DE500C" lfn="/alice/sim/PDC_08b/LHC08d10/50036/012/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12956949" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/012/AliESDs.root" type="f" />
- </event>
- <event name="13">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40774647" entryId="40774662" expiretime="" gowner="aliprod" guid="05198E42-4FFF-11DE-9717-7A6AB58CBEEF" guidtime="11DE4FFF" lfn="/alice/sim/PDC_08b/LHC08d10/50036/013/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12852838" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/013/AliESDs.root" type="f" />
- </event>
- <event name="14">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40608228" entryId="40608239" expiretime="" gowner="aliprod" guid="2E83C03C-4FFE-11DE-8001-D2C59A83BEEF" guidtime="11DE4FFE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/014/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13293740" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/014/AliESDs.root" type="f" />
- </event>
- <event name="15">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40702799" entryId="40702815" expiretime="" gowner="aliprod" guid="186B04CA-500A-11DE-9717-EC05A8C0BEEF" guidtime="11DE500A" lfn="/alice/sim/PDC_08b/LHC08d10/50036/015/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12745825" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/015/AliESDs.root" type="f" />
- </event>
- <event name="16">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40866536" entryId="40866542" expiretime="" gowner="aliprod" guid="9AA293AE-503C-11DE-8001-F5FFFF0ABEEF" guidtime="11DE503C" lfn="/alice/sim/PDC_08b/LHC08d10/50036/016/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13211762" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/016/AliESDs.root" type="f" />
- </event>
- <event name="17">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40752699" entryId="40752715" expiretime="" gowner="aliprod" guid="61D68EC4-4FEF-11DE-8001-9818010ABEEF" guidtime="11DE4FEF" lfn="/alice/sim/PDC_08b/LHC08d10/50036/017/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13084293" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/017/AliESDs.root" type="f" />
- </event>
- <event name="18">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40724580" entryId="40724596" expiretime="" gowner="aliprod" guid="FBDD6F60-5003-11DE-9717-D305A8C0BEEF" guidtime="11DE5003" lfn="/alice/sim/PDC_08b/LHC08d10/50036/018/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13135179" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/018/AliESDs.root" type="f" />
- </event>
- <event name="19">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40726035" entryId="40726051" expiretime="" gowner="aliprod" guid="663469C6-500A-11DE-8001-AF50A8C0BEEF" guidtime="11DE500A" lfn="/alice/sim/PDC_08b/LHC08d10/50036/019/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13014898" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/019/AliESDs.root" type="f" />
- </event>
- <event name="2">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40661719" entryId="40661735" expiretime="" gowner="aliprod" guid="FC12FA00-4FEE-11DE-8001-9718010ABEEF" guidtime="11DE4FEE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/002/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13072511" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/002/AliESDs.root" type="f" />
- </event>
- <event name="20">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40745174" entryId="40745190" expiretime="" gowner="aliprod" guid="4B6A5970-5000-11DE-9717-15E05D9FBEEF" guidtime="11DE5000" lfn="/alice/sim/PDC_08b/LHC08d10/50036/020/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12822516" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/020/AliESDs.root" type="f" />
- </event>
- <event name="21">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40866968" entryId="40866983" expiretime="" gowner="aliprod" guid="651CB574-503D-11DE-9717-C2E05D9FBEEF" guidtime="11DE503D" lfn="/alice/sim/PDC_08b/LHC08d10/50036/021/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13367820" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/021/AliESDs.root" type="f" />
- </event>
- <event name="22">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40797195" entryId="40797211" expiretime="" gowner="aliprod" guid="082A8E4E-4FF8-11DE-9717-1E000A0ABEEF" guidtime="11DE4FF8" lfn="/alice/sim/PDC_08b/LHC08d10/50036/022/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13337182" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/022/AliESDs.root" type="f" />
- </event>
- <event name="23">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40859027" entryId="40859039" expiretime="" gowner="aliprod" guid="4586BB44-5031-11DE-9717-BF74A8C0BEEF" guidtime="11DE5031" lfn="/alice/sim/PDC_08b/LHC08d10/50036/023/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12723670" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/023/AliESDs.root" type="f" />
- </event>
- <event name="24">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40786046" entryId="40786062" expiretime="" gowner="aliprod" guid="8074E9DA-5002-11DE-8001-70C59A83BEEF" guidtime="11DE5002" lfn="/alice/sim/PDC_08b/LHC08d10/50036/024/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12721991" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/024/AliESDs.root" type="f" />
- </event>
- <event name="25">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40697471" entryId="40697474" expiretime="" gowner="aliprod" guid="E0B754E0-5002-11DE-8001-C1C59A83BEEF" guidtime="11DE5002" lfn="/alice/sim/PDC_08b/LHC08d10/50036/025/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12740213" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/025/AliESDs.root" type="f" />
- </event>
- <event name="26">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40603957" entryId="40603970" expiretime="" gowner="aliprod" guid="F7C92D88-4FF9-11DE-8001-0A62ABC2BEEF" guidtime="11DE4FF9" lfn="/alice/sim/PDC_08b/LHC08d10/50036/026/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13015648" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/026/AliESDs.root" type="f" />
- </event>
- <event name="27">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40603975" entryId="40603988" expiretime="" gowner="aliprod" guid="1066BFB4-5003-11DE-8001-D7C59A83BEEF" guidtime="11DE5003" lfn="/alice/sim/PDC_08b/LHC08d10/50036/027/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12940124" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/027/AliESDs.root" type="f" />
- </event>
- <event name="28">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40603993" entryId="40604006" expiretime="" gowner="aliprod" guid="E5799554-4FFF-11DE-8001-66C59A83BEEF" guidtime="11DE4FFF" lfn="/alice/sim/PDC_08b/LHC08d10/50036/028/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13078050" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/028/AliESDs.root" type="f" />
- </event>
- <event name="29">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40604030" entryId="40604043" expiretime="" gowner="aliprod" guid="43ACD9C2-4FEE-11DE-8001-3DAD9E86BEEF" guidtime="11DE4FEE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/029/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12792236" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/029/AliESDs.root" type="f" />
- </event>
- <event name="3">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40681607" entryId="40681623" expiretime="" gowner="aliprod" guid="4F6076C8-5004-11DE-8001-DBB0F682BEEF" guidtime="11DE5004" lfn="/alice/sim/PDC_08b/LHC08d10/50036/003/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13394655" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/003/AliESDs.root" type="f" />
- </event>
- <event name="30">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40604048" entryId="40604061" expiretime="" gowner="aliprod" guid="F3DF7E66-4FEA-11DE-9717-99C68E80BEEF" guidtime="11DE4FEA" lfn="/alice/sim/PDC_08b/LHC08d10/50036/030/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13056035" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/030/AliESDs.root" type="f" />
- </event>
- <event name="31">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40724526" entryId="40724542" expiretime="" gowner="aliprod" guid="7B23D0AC-4FFA-11DE-9717-D2E05D9FBEEF" guidtime="11DE4FFA" lfn="/alice/sim/PDC_08b/LHC08d10/50036/031/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12956389" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/031/AliESDs.root" type="f" />
- </event>
- <event name="32">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40698947" entryId="40698963" expiretime="" gowner="aliprod" guid="0CA7AA34-5010-11DE-8001-7402010ABEEF" guidtime="11DE5010" lfn="/alice/sim/PDC_08b/LHC08d10/50036/032/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12629822" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/032/AliESDs.root" type="f" />
- </event>
- <event name="33">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40694105" entryId="40694121" expiretime="" gowner="aliprod" guid="70E5F638-4FF0-11DE-8001-9919010ABEEF" guidtime="11DE4FF0" lfn="/alice/sim/PDC_08b/LHC08d10/50036/033/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13268281" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/033/AliESDs.root" type="f" />
- </event>
- <event name="34">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40849966" entryId="40849982" expiretime="" gowner="aliprod" guid="2C03A336-501C-11DE-8001-19539E86BEEF" guidtime="11DE501C" lfn="/alice/sim/PDC_08b/LHC08d10/50036/034/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12763894" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/034/AliESDs.root" type="f" />
- </event>
- <event name="35">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40637244" entryId="40637258" expiretime="" gowner="aliprod" guid="8393CFE4-4FFF-11DE-8001-7309010ABEEF" guidtime="11DE4FFF" lfn="/alice/sim/PDC_08b/LHC08d10/50036/035/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13280927" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/035/AliESDs.root" type="f" />
- </event>
- <event name="36">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40673424" entryId="40673440" expiretime="" gowner="aliprod" guid="BCB51728-4FF1-11DE-8001-15AA9E86BEEF" guidtime="11DE4FF1" lfn="/alice/sim/PDC_08b/LHC08d10/50036/036/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13352897" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/036/AliESDs.root" type="f" />
- </event>
- <event name="37">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40859894" entryId="40859905" expiretime="" gowner="aliprod" guid="95084CD2-5027-11DE-8001-3019E793BEEF" guidtime="11DE5027" lfn="/alice/sim/PDC_08b/LHC08d10/50036/037/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12930117" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/037/AliESDs.root" type="f" />
- </event>
- <event name="38">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40622958" entryId="40622974" expiretime="" gowner="aliprod" guid="2819F70A-500F-11DE-9717-D704A8C0BEEF" guidtime="11DE500F" lfn="/alice/sim/PDC_08b/LHC08d10/50036/038/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12798199" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/038/AliESDs.root" type="f" />
- </event>
- <event name="39">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40788565" entryId="40788581" expiretime="" gowner="aliprod" guid="6878251A-4FF6-11DE-8001-6F6E010ABEEF" guidtime="11DE4FF6" lfn="/alice/sim/PDC_08b/LHC08d10/50036/039/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13709293" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/039/AliESDs.root" type="f" />
- </event>
- <event name="4">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40646904" entryId="40646920" expiretime="" gowner="aliprod" guid="BFDC19F4-4FFD-11DE-8001-6A08010ABEEF" guidtime="11DE4FFD" lfn="/alice/sim/PDC_08b/LHC08d10/50036/004/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13334006" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/004/AliESDs.root" type="f" />
- </event>
- <event name="40">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40769379" entryId="40769391" expiretime="" gowner="aliprod" guid="AB163FD0-4FF4-11DE-8001-7B72010ABEEF" guidtime="11DE4FF4" lfn="/alice/sim/PDC_08b/LHC08d10/50036/040/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13014676" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/040/AliESDs.root" type="f" />
- </event>
- <event name="41">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40804301" entryId="40804312" expiretime="" gowner="aliprod" guid="88D2B850-4FF8-11DE-8001-92B0F682BEEF" guidtime="11DE4FF8" lfn="/alice/sim/PDC_08b/LHC08d10/50036/041/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13501676" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/041/AliESDs.root" type="f" />
- </event>
- <event name="42">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40804337" entryId="40804353" expiretime="" gowner="aliprod" guid="CBC86BD4-4FF1-11DE-8001-65AA9E86BEEF" guidtime="11DE4FF1" lfn="/alice/sim/PDC_08b/LHC08d10/50036/042/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13447298" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/042/AliESDs.root" type="f" />
- </event>
- <event name="43">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40804355" entryId="40804371" expiretime="" gowner="aliprod" guid="86951C26-4FEF-11DE-8001-6719010ABEEF" guidtime="11DE4FEF" lfn="/alice/sim/PDC_08b/LHC08d10/50036/043/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12829045" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/043/AliESDs.root" type="f" />
- </event>
- <event name="44">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40804373" entryId="40804389" expiretime="" gowner="aliprod" guid="91188864-5007-11DE-8001-C5FAB796BEEF" guidtime="11DE5007" lfn="/alice/sim/PDC_08b/LHC08d10/50036/044/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13669059" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/044/AliESDs.root" type="f" />
- </event>
- <event name="45">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40693889" entryId="40693905" expiretime="" gowner="aliprod" guid="8D7C3516-4FFD-11DE-8001-7402010ABEEF" guidtime="11DE4FFD" lfn="/alice/sim/PDC_08b/LHC08d10/50036/045/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13323293" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/045/AliESDs.root" type="f" />
- </event>
- <event name="46">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40829575" entryId="40829591" expiretime="" gowner="aliprod" guid="90AABE5A-4FFE-11DE-8001-CFC59A83BEEF" guidtime="11DE4FFE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/046/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13318485" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/046/AliESDs.root" type="f" />
- </event>
- <event name="47">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40874200" entryId="40874209" expiretime="" gowner="aliprod" guid="7DB37CEE-5083-11DE-9717-06050B0ABEEF" guidtime="11DE5083" lfn="/alice/sim/PDC_08b/LHC08d10/50036/047/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12567924" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/047/AliESDs.root" type="f" />
- </event>
- <event name="48">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40779429" entryId="40779445" expiretime="" gowner="aliprod" guid="059E1244-4FFE-11DE-9717-0B020A0ABEEF" guidtime="11DE4FFE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/048/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12747565" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/048/AliESDs.root" type="f" />
- </event>
- <event name="49">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40786136" entryId="40786152" expiretime="" gowner="aliprod" guid="AE0A9EFE-4FF7-11DE-8001-371A9E86BEEF" guidtime="11DE4FF7" lfn="/alice/sim/PDC_08b/LHC08d10/50036/049/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12799950" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/049/AliESDs.root" type="f" />
- </event>
- <event name="5">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40646922" entryId="40646938" expiretime="" gowner="aliprod" guid="3DCBDE64-5001-11DE-8001-8008010ABEEF" guidtime="11DE5001" lfn="/alice/sim/PDC_08b/LHC08d10/50036/005/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13196292" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/005/AliESDs.root" type="f" />
- </event>
- <event name="50">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40866662" entryId="40866663" expiretime="" gowner="aliprod" guid="F129D6AC-503B-11DE-8001-537B4B86BEEF" guidtime="11DE503B" lfn="/alice/sim/PDC_08b/LHC08d10/50036/050/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13054347" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/050/AliESDs.root" type="f" />
- </event>
- <event name="51">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40747711" entryId="40747727" expiretime="" gowner="aliprod" guid="ADF3D1F8-4FF5-11DE-8001-DCCDBDD4BEEF" guidtime="11DE4FF5" lfn="/alice/sim/PDC_08b/LHC08d10/50036/051/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12995014" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/051/AliESDs.root" type="f" />
- </event>
- <event name="52">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40669144" entryId="40669160" expiretime="" gowner="aliprod" guid="CF8A1582-5001-11DE-9717-2AE05D9FBEEF" guidtime="11DE5001" lfn="/alice/sim/PDC_08b/LHC08d10/50036/052/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13779601" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/052/AliESDs.root" type="f" />
- </event>
- <event name="53">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40617740" entryId="40617753" expiretime="" gowner="aliprod" guid="88A07D68-4FEE-11DE-8001-B11A010ABEEF" guidtime="11DE4FEE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/053/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13283218" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/053/AliESDs.root" type="f" />
- </event>
- <event name="54">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40731755" entryId="40731771" expiretime="" gowner="aliprod" guid="882D3D2E-4FF0-11DE-8001-35AA9E86BEEF" guidtime="11DE4FF0" lfn="/alice/sim/PDC_08b/LHC08d10/50036/054/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13025824" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/054/AliESDs.root" type="f" />
- </event>
- <event name="55">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40794555" entryId="40794564" expiretime="" gowner="aliprod" guid="1832DA18-4FFC-11DE-8001-A7C59A83BEEF" guidtime="11DE4FFC" lfn="/alice/sim/PDC_08b/LHC08d10/50036/055/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13307629" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/055/AliESDs.root" type="f" />
- </event>
- <event name="56">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40657406" entryId="40657422" expiretime="" gowner="aliprod" guid="5A76798E-4FF7-11DE-9717-EF05A8C0BEEF" guidtime="11DE4FF7" lfn="/alice/sim/PDC_08b/LHC08d10/50036/056/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13422545" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/056/AliESDs.root" type="f" />
- </event>
- <event name="57">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40829593" entryId="40829609" expiretime="" gowner="aliprod" guid="BC47C880-4FF6-11DE-8001-62539E86BEEF" guidtime="11DE4FF6" lfn="/alice/sim/PDC_08b/LHC08d10/50036/057/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12685116" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/057/AliESDs.root" type="f" />
- </event>
- <event name="58">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40754139" entryId="40754155" expiretime="" gowner="aliprod" guid="ABC11A54-500D-11DE-8001-6F04010ABEEF" guidtime="11DE500D" lfn="/alice/sim/PDC_08b/LHC08d10/50036/059/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13054486" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/059/AliESDs.root" type="f" />
- </event>
- <event name="59">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40662581" entryId="40662597" expiretime="" gowner="aliprod" guid="445B1A24-5001-11DE-8001-CFC59A83BEEF" guidtime="11DE5001" lfn="/alice/sim/PDC_08b/LHC08d10/50036/060/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13252131" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/060/AliESDs.root" type="f" />
- </event>
- <event name="6">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40633557" entryId="40633572" expiretime="" gowner="aliprod" guid="CE8B0E22-4FF5-11DE-8001-DFA99E86BEEF" guidtime="11DE4FF5" lfn="/alice/sim/PDC_08b/LHC08d10/50036/006/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="11194332" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/006/AliESDs.root" type="f" />
- </event>
- <event name="60">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40735536" entryId="40735552" expiretime="" gowner="aliprod" guid="DA3E3C84-5000-11DE-8001-51C59A83BEEF" guidtime="11DE5000" lfn="/alice/sim/PDC_08b/LHC08d10/50036/061/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13462787" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/061/AliESDs.root" type="f" />
- </event>
- <event name="61">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40864562" entryId="40864566" expiretime="" gowner="aliprod" guid="673CA0A2-502F-11DE-9717-E3FB14ACBEEF" guidtime="11DE502F" lfn="/alice/sim/PDC_08b/LHC08d10/50036/062/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13657967" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/062/AliESDs.root" type="f" />
- </event>
- <event name="62">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40795902" entryId="40795918" expiretime="" gowner="aliprod" guid="778E6FD6-5006-11DE-9717-4C5530C1BEEF" guidtime="11DE5006" lfn="/alice/sim/PDC_08b/LHC08d10/50036/063/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="11882407" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/063/AliESDs.root" type="f" />
- </event>
- <event name="63">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40671302" entryId="40671318" expiretime="" gowner="aliprod" guid="2F8CF5F4-5006-11DE-8001-DDC59A83BEEF" guidtime="11DE5006" lfn="/alice/sim/PDC_08b/LHC08d10/50036/064/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13966029" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/064/AliESDs.root" type="f" />
- </event>
- <event name="64">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40693133" entryId="40693149" expiretime="" gowner="aliprod" guid="EBCE4D84-4FEE-11DE-8001-73AD9E86BEEF" guidtime="11DE4FEE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/065/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13342053" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/065/AliESDs.root" type="f" />
- </event>
- <event name="65">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40693151" entryId="40693167" expiretime="" gowner="aliprod" guid="79043A10-4FF4-11DE-8001-830A010ABEEF" guidtime="11DE4FF4" lfn="/alice/sim/PDC_08b/LHC08d10/50036/066/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12755188" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/066/AliESDs.root" type="f" />
- </event>
- <event name="66">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40785092" entryId="40785108" expiretime="" gowner="aliprod" guid="85579CC0-4FFF-11DE-8001-7307010ABEEF" guidtime="11DE4FFF" lfn="/alice/sim/PDC_08b/LHC08d10/50036/067/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12923247" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/067/AliESDs.root" type="f" />
- </event>
- <event name="67">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40654545" entryId="40654562" expiretime="" gowner="aliprod" guid="FA6BE908-4FF0-11DE-8001-AF17010ABEEF" guidtime="11DE4FF0" lfn="/alice/sim/PDC_08b/LHC08d10/50036/068/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13458172" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/068/AliESDs.root" type="f" />
- </event>
- <event name="68">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40688567" entryId="40688583" expiretime="" gowner="aliprod" guid="36AAD23A-4FF1-11DE-8001-9D18010ABEEF" guidtime="11DE4FF1" lfn="/alice/sim/PDC_08b/LHC08d10/50036/069/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13134088" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/069/AliESDs.root" type="f" />
- </event>
- <event name="69">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40624847" entryId="40624863" expiretime="" gowner="aliprod" guid="8F4BFA60-4FEF-11DE-8001-7518010ABEEF" guidtime="11DE4FEF" lfn="/alice/sim/PDC_08b/LHC08d10/50036/070/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13452011" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/070/AliESDs.root" type="f" />
- </event>
- <event name="7">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40633574" entryId="40633590" expiretime="" gowner="aliprod" guid="D7D5DE84-500A-11DE-8001-AF50A8C0BEEF" guidtime="11DE500A" lfn="/alice/sim/PDC_08b/LHC08d10/50036/007/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13281651" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/007/AliESDs.root" type="f" />
- </event>
- <event name="70">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40697921" entryId="40697937" expiretime="" gowner="aliprod" guid="ECCC7EEE-4FF4-11DE-9717-A074A8C0BEEF" guidtime="11DE4FF4" lfn="/alice/sim/PDC_08b/LHC08d10/50036/072/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12979534" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/072/AliESDs.root" type="f" />
- </event>
- <event name="71">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40787396" entryId="40787411" expiretime="" gowner="aliprod" guid="4059A672-4FF0-11DE-8001-AA18010ABEEF" guidtime="11DE4FF0" lfn="/alice/sim/PDC_08b/LHC08d10/50036/073/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13790285" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/073/AliESDs.root" type="f" />
- </event>
- <event name="72">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40766806" entryId="40766820" expiretime="" gowner="aliprod" guid="4F2839CA-4FF0-11DE-8001-6C18010ABEEF" guidtime="11DE4FF0" lfn="/alice/sim/PDC_08b/LHC08d10/50036/074/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13145552" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/074/AliESDs.root" type="f" />
- </event>
- <event name="73">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40766824" entryId="40766830" expiretime="" gowner="aliprod" guid="75EB3A92-5001-11DE-8001-7B02010ABEEF" guidtime="11DE5001" lfn="/alice/sim/PDC_08b/LHC08d10/50036/075/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13636333" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/075/AliESDs.root" type="f" />
- </event>
- <event name="74">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40835800" entryId="40835816" expiretime="" gowner="aliprod" guid="19A41D2A-4FFD-11DE-9717-3C000A0ABEEF" guidtime="11DE4FFD" lfn="/alice/sim/PDC_08b/LHC08d10/50036/076/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13180712" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/076/AliESDs.root" type="f" />
- </event>
- <event name="75">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40862705" entryId="40862722" expiretime="" gowner="aliprod" guid="763F8E82-502C-11DE-8001-3519E793BEEF" guidtime="11DE502C" lfn="/alice/sim/PDC_08b/LHC08d10/50036/077/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12953949" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/077/AliESDs.root" type="f" />
- </event>
- <event name="76">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40676013" entryId="40676029" expiretime="" gowner="aliprod" guid="57A72B1C-5007-11DE-8001-C5FAB796BEEF" guidtime="11DE5007" lfn="/alice/sim/PDC_08b/LHC08d10/50036/078/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13123362" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/078/AliESDs.root" type="f" />
- </event>
- <event name="77">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40748880" entryId="40748896" expiretime="" gowner="aliprod" guid="712FF322-4FF1-11DE-9717-B874A8C0BEEF" guidtime="11DE4FF1" lfn="/alice/sim/PDC_08b/LHC08d10/50036/079/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13109639" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/079/AliESDs.root" type="f" />
- </event>
- <event name="78">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40704182" entryId="40704198" expiretime="" gowner="aliprod" guid="381DAD86-4FF1-11DE-9717-B574A8C0BEEF" guidtime="11DE4FF1" lfn="/alice/sim/PDC_08b/LHC08d10/50036/080/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13380836" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/080/AliESDs.root" type="f" />
- </event>
- <event name="79">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40769811" entryId="40769827" expiretime="" gowner="aliprod" guid="81719CB0-4FFE-11DE-8001-6E01010ABEEF" guidtime="11DE4FFE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/081/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13338504" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/081/AliESDs.root" type="f" />
- </event>
- <event name="8">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40829557" entryId="40829573" expiretime="" gowner="aliprod" guid="9BD4B5AA-4FF5-11DE-8001-96B0F682BEEF" guidtime="11DE4FF5" lfn="/alice/sim/PDC_08b/LHC08d10/50036/008/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13188874" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/008/AliESDs.root" type="f" />
- </event>
- <event name="80">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40769829" entryId="40769842" expiretime="" gowner="aliprod" guid="1D8CE508-5001-11DE-8001-8903010ABEEF" guidtime="11DE5001" lfn="/alice/sim/PDC_08b/LHC08d10/50036/082/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12516746" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/082/AliESDs.root" type="f" />
- </event>
- <event name="81">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40751097" entryId="40751112" expiretime="" gowner="aliprod" guid="4C277C82-5003-11DE-8001-D7C59A83BEEF" guidtime="11DE5003" lfn="/alice/sim/PDC_08b/LHC08d10/50036/083/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13244517" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/083/AliESDs.root" type="f" />
- </event>
- <event name="82">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40726233" entryId="40726249" expiretime="" gowner="aliprod" guid="1BEB08FA-4FF3-11DE-8001-20A99E86BEEF" guidtime="11DE4FF3" lfn="/alice/sim/PDC_08b/LHC08d10/50036/084/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12421638" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/084/AliESDs.root" type="f" />
- </event>
- <event name="83">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40684072" entryId="40684088" expiretime="" gowner="aliprod" guid="BDE665A8-5004-11DE-8001-7003010ABEEF" guidtime="11DE5004" lfn="/alice/sim/PDC_08b/LHC08d10/50036/085/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13456866" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/085/AliESDs.root" type="f" />
- </event>
- <event name="84">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40681805" entryId="40681821" expiretime="" gowner="aliprod" guid="D7BC60E6-4FFE-11DE-8001-7901010ABEEF" guidtime="11DE4FFE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/086/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12985667" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/086/AliESDs.root" type="f" />
- </event>
- <event name="85">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40689377" entryId="40689393" expiretime="" gowner="aliprod" guid="0F099A5A-500E-11DE-8001-29539E86BEEF" guidtime="11DE500E" lfn="/alice/sim/PDC_08b/LHC08d10/50036/087/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13427938" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/087/AliESDs.root" type="f" />
- </event>
- <event name="86">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40676409" entryId="40676425" expiretime="" gowner="aliprod" guid="3047E370-4FF0-11DE-8001-7E18010ABEEF" guidtime="11DE4FF0" lfn="/alice/sim/PDC_08b/LHC08d10/50036/088/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13012063" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/088/AliESDs.root" type="f" />
- </event>
- <event name="87">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40659401" entryId="40659417" expiretime="" gowner="aliprod" guid="1A51B6B4-4FFE-11DE-8001-4B7B4B86BEEF" guidtime="11DE4FFE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/089/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13417110" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/089/AliESDs.root" type="f" />
- </event>
- <event name="88">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40868309" entryId="40868322" expiretime="" gowner="aliprod" guid="DDED5DE4-5049-11DE-8001-B2B4F682BEEF" guidtime="11DE5049" lfn="/alice/sim/PDC_08b/LHC08d10/50036/090/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13131600" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/090/AliESDs.root" type="f" />
- </event>
- <event name="89">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40829611" entryId="40829627" expiretime="" gowner="aliprod" guid="76665B70-4FF0-11DE-8001-9C18010ABEEF" guidtime="11DE4FF0" lfn="/alice/sim/PDC_08b/LHC08d10/50036/091/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12743482" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/091/AliESDs.root" type="f" />
- </event>
- <event name="9">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40631092" entryId="40631108" expiretime="" gowner="aliprod" guid="8E1127F6-5003-11DE-9717-BCE05D9FBEEF" guidtime="11DE5003" lfn="/alice/sim/PDC_08b/LHC08d10/50036/009/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13359744" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/009/AliESDs.root" type="f" />
- </event>
- <event name="90">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40781693" entryId="40781709" expiretime="" gowner="aliprod" guid="065E0BD2-4FF7-11DE-8001-9BB0F682BEEF" guidtime="11DE4FF7" lfn="/alice/sim/PDC_08b/LHC08d10/50036/092/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12857299" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/092/AliESDs.root" type="f" />
- </event>
- <event name="91">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40653340" entryId="40653355" expiretime="" gowner="aliprod" guid="945F2148-4FF0-11DE-8001-7CAA9E86BEEF" guidtime="11DE4FF0" lfn="/alice/sim/PDC_08b/LHC08d10/50036/093/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12787987" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/093/AliESDs.root" type="f" />
- </event>
- <event name="92">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40727601" entryId="40727617" expiretime="" gowner="aliprod" guid="0ED8A1FE-5007-11DE-8001-7003010ABEEF" guidtime="11DE5007" lfn="/alice/sim/PDC_08b/LHC08d10/50036/094/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13134432" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/094/AliESDs.root" type="f" />
- </event>
- <event name="93">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40676499" entryId="40676515" expiretime="" gowner="aliprod" guid="1CB73D38-5009-11DE-8001-A1B0F682BEEF" guidtime="11DE5009" lfn="/alice/sim/PDC_08b/LHC08d10/50036/095/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13209577" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/095/AliESDs.root" type="f" />
- </event>
- <event name="94">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40710424" entryId="40710440" expiretime="" gowner="aliprod" guid="27B0624C-4FF4-11DE-8001-78A99E86BEEF" guidtime="11DE4FF4" lfn="/alice/sim/PDC_08b/LHC08d10/50036/096/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13257355" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/096/AliESDs.root" type="f" />
- </event>
- <event name="95">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40866554" entryId="40866557" expiretime="" gowner="aliprod" guid="267D547E-5036-11DE-9717-FB05A8C0BEEF" guidtime="11DE5036" lfn="/alice/sim/PDC_08b/LHC08d10/50036/097/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13591633" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/097/AliESDs.root" type="f" />
- </event>
- <event name="96">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40700227" entryId="40700243" expiretime="" gowner="aliprod" guid="F3B62D2C-4FEF-11DE-8001-B21A010ABEEF" guidtime="11DE4FEF" lfn="/alice/sim/PDC_08b/LHC08d10/50036/098/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="13268494" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/098/AliESDs.root" type="f" />
- </event>
- <event name="97">
- <file name="AliESDs.root" aclId="" broken="0" ctime="2009-08-03 23:13:24" dir="40614699" entryId="40614714" expiretime="" gowner="aliprod" guid="DD5DC93C-4FEE-11DE-8001-03AE9E86BEEF" guidtime="11DE4FEE" lfn="/alice/sim/PDC_08b/LHC08d10/50036/099/AliESDs.root" md5="" owner="aliprod" perm="755" replicated="0" size="12847272" turl="alien:///alice/sim/PDC_08b/LHC08d10/50036/099/AliESDs.root" type="f" />
- </event>
- <info command="[/alice/cern.ch/user/k/kread/]: find -z -x mycollect /alice/sim/PDC_08b/LHC08d10/50036/0* AliESDs.root" creator="admin" date="Sun Sep 27 07:21:33 CEST 2009" timestamp="1254028893" />
- </collection>
-</alien>
+++ /dev/null
-
-//
-// AliEn Multi MasterJob Coordination
-// Author: K. Read
-//
-
-void mylauncher()
-{
- TString worksubdir = "work12"; //already existing work directory
-//TString worksubdir = "LHC09b4"; //a better name for MultiMode
- // Name of JDL file to upload. Leave blank for no upload and no submit.
- TString jdlfilename = "anaJete.jdl";//ignored for MultiMode
- // Name of executable to upload. Leave blank for no upload.
- TString collectfilename = "mycollect.xml";//ignored for MultiMode
- // List any other files to upload in string filenames separated by blanks.
- TString execfilename = "";
- // List any other files to upload in string filenames separated by blanks.
- TString filelist = "run.C validate.sh anaJete.C ConfigJetAnalysisFastJet.C ConfigAnalysisElectron.C mergeout.jdl mergeoutscaled.jdl ANALYSIS.par ANALYSISalice.par AOD.par EMCALUtils.par ESD.par PHOSUtils.par STEERBase.par JETAN.par FASTJETAN.par PWG4PartCorrBase.par PWG4PartCorrDep.par";
- TString filename;
-
- Bool_t kMultiMode = kFALSE;//uncomment for regular mode
-//Bool_t kMultiMode = kTRUE; //uncomment for coordinated multiple masterjobs
-//MultiMode Variables. All ignored for regular mode.
- TString multijdlfilename = "anaJetemulti.jdl";
-//TString datatype = "AliESDs.root";//uncomment for ESDs
- TString datatype = "AliAOD.root";//uncomment for AODs
- TString findlist = "
-/alice/sim/LHC09b4/AOD/000
-/alice/sim/LHC09b4/AOD/001
-/alice/sim/LHC09b4/AOD/002"; //line-feed separated list of desired collections
-
-
- gSystem->Load("libNetx.so") ;
- gSystem->Load("libRAliEn.so");
-
- TGrid::Connect("alien://") ;
- if (gGrid && gGrid->IsConnected()) {
- TString homedir = gGrid->GetHomeDirectory(); // has a trailing slash
- TString workdir = homedir + worksubdir;
- if (gGrid->Cd(workdir)) {
-
- // Upload files listed
- if (filelist.Length()) {
- arr = filelist.Tokenize(" ");
- TObjString *os;
- TIter next(arr);
- while ((os=(TObjString*)next())) {
- Info("Grid Upload", "Copying %s to your AliEn work directory", os->GetString().Data());
- if (FileExists(os->GetString())) gGrid->Rm(os->GetString());
- TFile::Cp(Form("file:%s",os->GetString().Data()), Form("alien://%s/%s", workdir.Data(), os->GetString().Data()));
- }
- delete arr;
- }
-
- // Prepare collection(s)
- if(kMultiMode){
- // Find collections if desired
- if (findlist.Length()) {
- arr = findlist.Tokenize("\n");
- TObjString *os;
- TIter next(arr);
- Int_t count=0;
- while ((os=(TObjString*)next())) {
- TGridResult *res;
- Info("Find", "Making collection for %s", os->GetString().Data());
- res = gGrid->Command(Form("find -x collection %s %s > mycollect%02d.xml", os->GetString().Data(),datatype.Data(),count));
- if(!res) Info("Find", "Failed making collection for %s \n",os->GetString().Data());
- delete res;
- count++;
- }
- }
- }
- else{
- if (collectfilename.Length()) {
- Info("Grid Upload", "Copying %s to your AliEn work directory", collectfilename.Data());
- if (FileExists(collectfilename)) gGrid->Rm(collectfilename);
- TFile::Cp(Form("file:%s",collectfilename.Data()), Form("alien://%s/%s", workdir.Data(), collectfilename.Data()));
- }
- }
-
- // Upload executable if listed
- if (execfilename.Length()) {
- filename = Form("%sbin/%s", homedir.Data(), execfilename.Data());
- if (FileExists(filename)) gGrid->Rm(filename);
- Info("Grid Upload", "Copying executable file %s to your AliEn bin directory", execfilename.Data());
- TFile::Cp(Form("file:%s",execfilename.Data()), Form("alien://%s", filename.Data()));
- }
-
- // Upload and submit JDL if listed
- if (kMultiMode) jdlfilename = multijdlfilename;
- if (jdlfilename.Length()) {
- filename = Form("%s/%s", workdir.Data(), jdlfilename.Data());
- if (FileExists(filename)) gGrid->Rm(filename);
- Info("Grid Upload", "Copying JDL file %s to your AliEn work directory", jdlfilename.Data());
- TFile::Cp(Form("file:%s",jdlfilename.Data()), Form("alien://%s", filename.Data()));
-
- if (!kMultiMode) Int_t count=1;
- for (Int_t index = 0; index < count; index++) {
- TGridResult *res;
- TString jobID = "";
- if(kMultiMode){
- res = gGrid->Command(Form("submit %s %s %02d", jdlfilename.Data(),worksubdir.Data(),index));
- Info("Launcher:", "Submitting %s %s %02d", jdlfilename.Data(),worksubdir.Data(),index);
- }
- else{
- res = gGrid->Command(Form("submit %s", jdlfilename.Data()));
- Info("Launcher:", "Submitting %s", jdlfilename.Data());
- }
- if (res) {
- const char *cjobId = res->GetKey(0,"jobId");
- if (!cjobId) {
- Error("Launcher:", "Your JDL %s could not be submitted", jdlfilename.Data());
- return;
- }
- else {
- Info("Launcher:", "Your JDL %s was successfully submitted.\n\n\t\t\t THE JOB ID IS: %s\n",
- jdlfilename.Data(), cjobId);
- }
- delete res;
- }
- }
- }
-
- // Launch alien shell
- gSystem->Exec("aliensh");
-
- }
- }
-}
-
-Bool_t FileExists(const char *lfn) const
-{
-// Returns true if file exists.
- if (!gGrid) {
- Error("FileExists", "No connection to grid");
- return kFALSE;
- }
- TGridResult *res = gGrid->Ls(lfn);
- if (!res) return kFALSE;
- TMap *map = dynamic_cast<TMap*>(res->At(0));
- if (!map) {
- delete res;
- return kFALSE;
- }
- TObjString *objs = dynamic_cast<TObjString*>(map->GetValue("name"));
- if (!objs || !objs->GetString().Length()) {
- delete res;
- return kFALSE;
- }
- delete res;
- return kTRUE;
-}
+++ /dev/null
-
-//
-// AliEn Initial PYTHIA pt hard bin-by-bin merging to prepare for rescaling
-// Author: K. Read
-//
-
-void mymerger()
-{
- TString filename;
- TString jdlfilename = "mergeoutscaledi.jdl";
- //TString worksubdir = "LHC09b4b";
- //Int_t maxbin = 15;
- TString worksubdir = "LHC09b2ESDb";
- Int_t maxbin = 16;
-
- //Later do: aliensh; cd worksubdir/output/merged; cp histoss* file:
- //And then process files locally with MergeFileInBins.C.
-
- gSystem->Load("libNetx.so") ;
- gSystem->Load("libRAliEn.so");
-
- TGrid::Connect("alien://") ;
- if (gGrid && gGrid->IsConnected()) {
- TString homedir = gGrid->GetHomeDirectory(); // has a trailing slash
- TString workdir = homedir + worksubdir;
- if (gGrid->Cd(workdir)) {
-
- // Upload and submit JDL if listed
- if (jdlfilename.Length()) {
- filename = Form("%s/%s", workdir.Data(), jdlfilename.Data());
- if (FileExists(filename)) gGrid->Rm(filename);
- Info("Grid Upload", "Copying JDL file %s to your AliEn work directory", jdlfilename.Data());
- TFile::Cp(Form("file:%s",jdlfilename.Data()), Form("alien://%s", filename.Data()));
-
- for (Int_t index = 0; index <= maxbin; index++) {
- TGridResult *res;
- TString jobID = "";
- res = gGrid->Command(Form("submit %s %s %02d", jdlfilename.Data(),worksubdir.Data(),index));
- Info("Launcher:", "Submitting %s %s %02d", jdlfilename.Data(),worksubdir.Data(),index);
- if (res) {
- const char *cjobId = res->GetKey(0,"jobId");
- if (!cjobId) {
- Error("Launcher:", "Your JDL %s could not be submitted", jdlfilename.Data());
- return;
- }
- else {
- Info("Launcher:", "Your JDL %s was successfully submitted.\n\n\t\t\t THE JOB ID IS: %s\n",
- jdlfilename.Data(), cjobId);
- }
- delete res;
- }
- }
- }
-
- // Launch alien shell
- gSystem->Exec("aliensh");
-
- }
- }
-}
-
-Bool_t FileExists(const char *lfn) const
-{
-// Returns true if file exists.
- if (!gGrid) {
- Error("FileExists", "No connection to grid");
- return kFALSE;
- }
- TGridResult *res = gGrid->Ls(lfn);
- if (!res) return kFALSE;
- TMap *map = dynamic_cast<TMap*>(res->At(0));
- if (!map) {
- delete res;
- return kFALSE;
- }
- TObjString *objs = dynamic_cast<TObjString*>(map->GetValue("name"));
- if (!objs || !objs->GetString().Length()) {
- delete res;
- return kFALSE;
- }
- delete res;
- return kTRUE;
-}
+++ /dev/null
-{
- TFile *fin = TFile::Open("histos.root");
- TCanvas *myc1 = new TCanvas("myc1","myc1",1);
-// myc1->SetLogx(1);
-// myc1->SetLogy(1);
- histos.Print();
- histos->FindObject("AnaElectron_hPtElectron")->Draw("pe");
-//AnaElectron_hPtElectron->GetXaxis()->SetRangeUser(1,200);
-//AnaElectron_hPtElectron->GetYaxis()->SetRangeUser(1,4000);
-//AnaElectron_hPtElectron->SetTitle("pT of Electron");
-//AnaElectron_hPtElectron->SetStats(kFALSE);
-//AnaElectron_hPtElectron->Draw("pe");
-
-// myc1->Print("electron.pdf");
-}
+++ /dev/null
-{
- TFile *fin = TFile::Open("histosscaled.root");
- TCanvas *myc1 = new TCanvas("myc1","myc1",1);
-// myc1->SetLogx(1);
-// myc1->SetLogy(1);
- histosscaled.Print();
- histosscaled->FindObject("AnaElectron_hPtElectronScaled")->Draw("pe");
-//AnaElectron_hPtElectronScaled->GetXaxis()->SetRangeUser(1,200);
-// AnaElectron_hPtElectronScaled->GetYaxis()->SetRangeUser(1,4000);
-// AnaElectron_hPtElectronScaled->SetTitle("pT of Electron");
-// AnaElectron_hPtElectronScaled->SetStats(kFALSE);
-// AnaElectron_hPtElectronScaled->Draw("pe");
-
-// myc1->Print("electron.pdf");
-}
+++ /dev/null
-#!/bin/bash
-
-#
-# Author: K. Read
-#
-# This script makes PAR files without requiring write access to
-# $ALICE_ROOT. It is not necessary to first make AliRoot after doing
-# an svn checkout or update. Execute this script in your local area
-# with write access.
-#
-# Usage:
-# parmaker JETAN
-# parmaker JETAN remote
-# parmaker JETAN remote patch
-# parmaker all remote
-#
-# So far only available for these par files:
-# ANALYSIS ANALYSISalice AOD EMCALUtils ESD FASTJETAN JETAN
-# PHOSUtils PWG4PartCorrBase PWG4PartCorrDep STEERBase
-#
-
-case $2 in
-
- "remote")
-
- case $1 in
- "all")
- ./parmaker ANALYSIS remote
- ./parmaker ANALYSISalice remote
- ./parmaker AOD remote
- ./parmaker EMCALUtils remote
- ./parmaker ESD remote
- ./parmaker FASTJETAN remote
- ./parmaker JETAN remote
- ./parmaker PHOSUtils remote
- ./parmaker PWG4PartCorrBase remote
- ./parmaker PWG4PartCorrDep remote
- ./parmaker STEERBase remote
- echo " "
- echo "finished making complete set of par files"
- exit
- ;;
- "ANALYSIS")
- parmaker_input_basedir="ANALYSIS"
- parmaker_input_dir_subdir="ANALYSIS"
- parmaker_output_dir_subdir="ANALYSIS"
- ;;
- "ANALYSISalice")
- parmaker_input_basedir="ANALYSIS"
- parmaker_input_dir_subdir="ANALYSIS"
- parmaker_output_dir_subdir="ANALYSISalice"
- ;;
- "AOD")
- parmaker_input_basedir="STEER"
- parmaker_input_dir_subdir="STEER"
- parmaker_output_dir_subdir="AOD"
- ;;
- "EMCALUtils")
- parmaker_input_basedir="EMCAL"
- parmaker_input_dir_subdir="EMCAL"
- parmaker_output_dir_subdir="EMCALUtils"
- ;;
- "ESD")
- parmaker_input_basedir="STEER"
- parmaker_input_dir_subdir="STEER"
- parmaker_output_dir_subdir="ESD"
- ;;
- "PHOSUtils")
- parmaker_input_basedir="PHOS"
- parmaker_input_dir_subdir="PHOS"
- parmaker_output_dir_subdir="PHOSUtils"
- ;;
- "PWG4PartCorrBase")
- parmaker_input_basedir="PWG4"
- parmaker_input_dir_subdir="PWG4/PartCorrBase"
- parmaker_output_dir_subdir="PWG4PartCorrBase/PartCorrBase"
- ;;
- "PWG4PartCorrDep")
- parmaker_input_basedir="PWG4"
- parmaker_input_dir_subdir="PWG4/PartCorrDep"
- parmaker_output_dir_subdir="PWG4PartCorrDep/PartCorrDep"
- ;;
- "STEERBase")
- parmaker_input_basedir="STEER"
- parmaker_input_dir_subdir="STEER"
- parmaker_output_dir_subdir="STEERBase"
- ;;
- "JETAN")
- parmaker_input_basedir="JETAN"
- parmaker_input_dir_subdir="JETAN"
- parmaker_output_dir_subdir="JETAN"
- ;;
- "FASTJETAN")
- parmaker_input_basedir="JETAN"
- parmaker_input_dir_subdir="JETAN"
- parmaker_output_dir_subdir="FASTJETAN"
- ;;
- *)
- echo "parmaker: I'm sorry Dave, I'm afraid I can't do that."
- exit
- esac
-
-
- echo " "
- echo "parmaker to use source $ALICE_ROOT/$parmaker_input_dir_subdir"
- if [ -e "$ALICE_ROOT/$parmaker_input_dir_subdir" ]
- then
- echo "parmaker creating $1.par"
- mkdir $1
- if [ $parmaker_input_basedir != $parmaker_input_dir_subdir ]
- then
- mkdir $parmaker_output_dir_subdir
- fi
-
- list=`grep Ali ${ALICE_ROOT}/${parmaker_input_basedir}/lib${1}.pkg | sed -e '/EXPORT/d' -e '/EINCLUDE/d' -e 's:.cxx::g' -e 's:SRCS::' -e 's:=::' -e 's:+::' -e 's:\\\::'`
- for i in $list; do
-# echo "adding:" $i
- cp $ALICE_ROOT/$parmaker_input_basedir/$i.cxx $parmaker_output_dir_subdir
- cp $ALICE_ROOT/$parmaker_input_basedir/$i.h $parmaker_output_dir_subdir
- done
-
- mkdir $1/PROOF-INF
- cp -r $ALICE_ROOT/$parmaker_input_basedir/PROOF-INF.$1/* $1/PROOF-INF
- cp $ALICE_ROOT/$parmaker_input_basedir/${1}LinkDef.h $1
- cp $ALICE_ROOT/$parmaker_input_basedir/lib${1}.pkg $1
- cp $ROOTSYS/test/Makefile.arch $1
- cp $ALICE_ROOT/$parmaker_input_basedir/Makefile $1/Makefiletemp
- sed -e 's:include \$(ROOTSYS)\/test\/Makefile.arch:include Makefile.arch:' -e "s:PACKAGE = .*:PACKAGE = ${1}:" $1/Makefiletemp > $1/Makefile
-
- /bin/rm $1/Makefiletemp
-
-
-# Apply temporary patches.
- if [ -z $3 ]
- then
- echo "no non-standard patches applied"
- else
- if [ $3 = "patch" ]
- then
- if [ -e "patches/$1" ]
- then
- echo "applying temporary $1 patch"
- cp patches/$1/* $parmaker_output_dir_subdir
- echo "completed temporary $1 patch"
- else
- echo "no patches found for $1"
- fi
- else
- echo "argument $3 not processed"
- fi
- fi
-
-# Apply further FASTJETAN tweaks.
- if [ $1 = "FASTJETAN" ]
- then
- echo "Making temporary adjustments to FASTJETAN.par."
- cp -r $ALICE_ROOT/JETAN/fastjet FASTJETAN
-
- sed -e "s:CXXFLAGS +:\nALICEINC += -I../JETAN -I./fastjet -I./fastjet/fastjet\n\n\nCXXFLAGS +:" FASTJETAN/Makefile > FASTJETAN/Makefile2
- mv FASTJETAN/Makefile2 FASTJETAN/Makefile
- fi
-
- tar cfzh $1.par $1
- /bin/rm -rf $1
- fi
- ;;
-
- *)
- echo "parmaker to use local source $1"
- if [ -e "$1" ]
- then
- echo "parmaker creating $1.par"
- tar cfzh $1.par $1
- else
- echo "local subdirectory $1 not found"
- fi
-
-esac
+++ /dev/null
-
-
-{
- //=========Macro generated from canvas: c7/7 efficiency per pt
- //========= (Sat Oct 24 13:36:23 2009) by ROOT version5.14/00
- TCanvas *c7 = new TCanvas("c7", "7 efficiency per pt",158,233,600,400);
- c7->Range(-5,-0.08875,55,0.89875);
- c7->SetFillColor(0);
- c7->SetBorderMode(0);
- c7->SetBorderSize(2);
- c7->SetFrameBorderMode(0);
- c7->SetFrameBorderMode(0);
-
- TGraphErrors *gre = new TGraphErrors(5);
- gre->SetName("Graph");
- gre->SetTitle("PID efficiency");
- gre->SetFillColor(1);
- gre->SetMarkerStyle(20);
- gre->SetPoint(0,5,0.140749);
- gre->SetPointError(0,0,0);
- gre->SetPoint(1,15,0.559721);
- gre->SetPointError(1,0,0);
- gre->SetPoint(2,25,0.465536);
- gre->SetPointError(2,0,0);
- gre->SetPoint(3,35,0.418012);
- gre->SetPointError(3,0,0);
- gre->SetPoint(4,45,0.341103);
- gre->SetPointError(4,0,0);
-
- TH1 *Graph1 = new TH1F("Graph1","PID efficiency",100,1,49);
- Graph1->SetMinimum(0.01);
- Graph1->SetMaximum(0.8);
- Graph1->SetDirectory(0);
- Graph1->SetStats(0);
- Graph1->GetXaxis()->SetTitle("pT (GeV/c)");
- Graph1->GetYaxis()->SetTitle("efficiency");
- gre->SetHistogram(Graph1);
-
- gre->Draw("apl");
-
- gre = new TGraphErrors(5);
- gre->SetName("Graph2");
- gre->SetTitle("Graph2");
- gre->SetFillColor(1);
- gre->SetMarkerColor(2);
- gre->SetMarkerStyle(20);
- gre->SetPoint(0,5,0.136597);
- gre->SetPointError(0,0,0);
- gre->SetPoint(1,15,0.525284);
- gre->SetPointError(1,0,0);
- gre->SetPoint(2,25,0.425577);
- gre->SetPointError(2,0,0);
- gre->SetPoint(3,35,0.397009);
- gre->SetPointError(3,0,0);
- gre->SetPoint(4,45,0.316323);
- gre->SetPointError(4,0,0);
- gre->Draw("pl");
-
- TLegend *leg = new TLegend(0.5,0.6,0.85,0.85,NULL,"brNDC");
- leg->SetBorderSize(1);
- leg->SetLineColor(1);
- leg->SetLineStyle(1);
- leg->SetLineWidth(1);
- leg->SetFillColor(10);
- leg->SetFillStyle(1001);
- TLegendEntry *entry=leg->AddEntry("Graph","Identified Electrons (EMCAL)","p");
- entry->SetLineColor(1);
- entry->SetLineStyle(1);
- entry->SetLineWidth(1);
- entry->SetMarkerColor(1);
- entry->SetMarkerStyle(20);
- entry->SetMarkerSize(1);
- entry->SetTextAlign(12);
- entry->SetTextColor(1);
- entry=leg->AddEntry("Graph2","Identified Electrons (EMCAL+TRD+TPC)","p");
- entry->SetLineColor(1);
- entry->SetLineStyle(1);
- entry->SetLineWidth(1);
- entry->SetMarkerColor(2);
- entry->SetMarkerStyle(20);
- entry->SetMarkerSize(1);
- entry->SetTextAlign(12);
- entry->SetTextColor(1);
- leg->Draw();
- TLatex * tex = new TLatex(0.4,0.91,"LHC09b4, 5.5Tev, b-jet");
- tex->SetNDC();
- tex->SetTextSize(0.04);
- tex->SetLineWidth(2);
- tex->Draw();
- tex = new TLatex(0.4,0.85,"oct20");
- tex->SetNDC();
- tex->SetTextSize(0.04);
- tex->SetLineWidth(2);
- tex->Draw();
- tex = new TLatex(0.15,0.85,"version = 99");
- tex->SetNDC();
- tex->SetTextSize(0.04);
- tex->SetLineWidth(2);
- tex->Draw();
- tex = new TLatex(0.15,0.75,"nEvents=406580");
- tex->SetNDC();
- tex->SetTextSize(0.04);
- tex->SetLineWidth(2);
- tex->Draw();
-
- TPaveText *pt = new TPaveText(0.01,0.940161,0.217919,0.995,"blNDC");
- pt->SetName("title");
- pt->SetBorderSize(1);
- text = pt->AddText("PID efficiency");
- pt->Draw();
- c7->Modified();
- c7->cd();
- c7->SetSelected(c7);
-}
+++ /dev/null
-{
- //=========Macro generated from canvas: c3/cr
- //========= (Sat Oct 24 14:06:41 2009) by ROOT version5.23/01
- TCanvas *c3 = new TCanvas("c3", "cr",0,0,1000,600);
- gStyle->SetOptStat(0);
- c3->Range(0,0,1,1);
- c3->SetFillColor(0);
- c3->SetBorderMode(0);
- c3->SetBorderSize(2);
- c3->SetFrameBorderMode(0);
-
- TGraphBentErrors *grbe = new TGraphBentErrors(4);
- grbe->SetName("Graph");
- grbe->SetTitle("Graph");
- grbe->SetFillColor(1);
- grbe->SetMarkerColor(2);
- grbe->SetMarkerStyle(21);
- grbe->SetPoint(0,10,1163);
- grbe->SetPointError(0,0,0,319.8685,242.6435,0,0,0,0);
- grbe->SetPoint(1,20,1473);
- grbe->SetPointError(1,0,0,355.2512,271.8543,0,0,0,0);
- grbe->SetPoint(2,40,2172);
- grbe->SetPointError(2,0,0,490.1327,376.0158,0,0,0,0);
- grbe->SetPoint(3,80,3329);
- grbe->SetPointError(3,0,0,800.9482,426.2069,0,0,0,0);
-
- TH1 *Graph1 = new TH1F("Graph1","",100,3,87);
- Graph1->SetMinimum(10);
- Graph1->SetMaximum(5000);
- Graph1->SetDirectory(0);
- Graph1->SetStats(0);
- Graph1->GetXaxis()->SetTitle("p_{T}^{track} (GeV/c)");
- // Graph1->GetXaxis()->SetTitleOffset(1.5);
- Graph1->GetYaxis()->SetTitle("Rejection Power");
- //Graph1->GetYaxis()->SetTitleOffset(1.5);
- grbe->SetHistogram(Graph1);
-
- grbe->Draw("apl");
-
- grbe = new TGraphBentErrors(4);
- grbe->SetName("Grapha");
- grbe->SetTitle("Grapha");
- grbe->SetFillColor(1);
- grbe->SetMarkerColor(4);
- grbe->SetMarkerStyle(21);
- grbe->SetPoint(0,10,214);
- grbe->SetPointError(0,0,0,23.04028,20.36199,0,0,0,0);
- grbe->SetPoint(1,20,600);
- grbe->SetPointError(1,0,0,86.54205,72.02112,0,0,0,0);
- grbe->SetPoint(2,40,651);
- grbe->SetPointError(2,0,0,77.49317,70.38893,0,0,0,0);
- grbe->SetPoint(3,80,1705);
- grbe->SetPointError(3,0,0,249.8343,221.2924,0,0,0,0);
-
- TH1 *Graph2 = new TH1F("Graph2","Rejection Power",100,3,87);
- Graph2->SetMinimum(17.42646);
- Graph2->SetMaximum(2099.826);
- Graph2->SetDirectory(0);
- Graph2->SetStats(0);
- Graph2->GetXaxis()->SetTitle("p_{T}^{track} (GeV/c)");
- Graph2->GetYaxis()->SetTitle("Rejection Power");
- grbe->SetHistogram(Graph2);
-
- grbe->Draw("pl");
-
- grbe = new TGraphBentErrors(2);
- grbe->SetName("Graphb");
- grbe->SetTitle("Graphb");
- grbe->SetFillColor(1);
- grbe->SetMarkerColor(6);
- grbe->SetMarkerStyle(29);
- grbe->SetMarkerSize(2);
- grbe->SetPoint(0,40,686);
- grbe->SetPointError(0,0,0,0,0,0,0,0,0);
- grbe->SetPoint(1,80,1210);
- grbe->SetPointError(1,0,0,0,0,0,0,0,0);
-
- TH1 *Graph3 = new TH1F("Graph3","Rejection Power",100,36,84);
- Graph3->SetMinimum(633.6);
- Graph3->SetMaximum(1262.4);
- Graph3->SetDirectory(0);
- Graph3->SetStats(0);
- Graph3->GetXaxis()->SetTitle("p_{T}^{track} (GeV/c)");
- Graph3->GetYaxis()->SetTitle("Rejection Power");
- grbe->SetHistogram(Graph3);
-
- grbe->Draw("pl");
-
- TLegend *leg = new TLegend(0.15,0.7,0.7,0.9,NULL,"brNDC");
- leg->SetBorderSize(1);
- leg->SetTextFont(62);
- leg->SetTextSize(0.03);
- leg->SetLineColor(1);
- leg->SetLineStyle(1);
- leg->SetLineWidth(1);
- leg->SetFillColor(0);
- leg->SetFillStyle(1001);
- TLegendEntry *entry=leg->AddEntry("Graph","80% e^{-} Efficiency - Simulation","p");
- entry->SetLineColor(1);
- entry->SetLineStyle(1);
- entry->SetLineWidth(1);
- entry->SetMarkerColor(1);
- entry->SetMarkerStyle(21);
- entry->SetMarkerSize(1);
- entry=leg->AddEntry("Grapha","90% e^{-} Efficiency - Simulation","p");
- entry->SetLineColor(1);
- entry->SetLineStyle(1);
- entry->SetLineWidth(1);
- entry->SetMarkerColor(1);
- entry->SetMarkerStyle(21);
- entry->SetMarkerSize(1);
- entry=leg->AddEntry("Graphb","90% e^{-} Efficiency - Test Beam","p");
- entry->SetLineColor(1);
- entry->SetLineStyle(1);
- entry->SetLineWidth(1);
- entry->SetMarkerColor(1);
- entry->SetMarkerStyle(21);
- entry->SetMarkerSize(1);
- leg->Draw();
- // c3->Modified();
- //c3->cd();
- //c3->SetSelected(c3);
-}
+++ /dev/null
-/////////////////////////////////////////////////\r
-//\r
-// Macro for plotting MC rates of electrons\r
-// for the EMCAL PPR\r
-//\r
-// J.L. Klay (Cal Poly)\r
-//\r
-/////////////////////////////////////////////////\r
-\r
-TLegend* leg;\r
-\r
-void plotMCRates() {\r
-\r
- gROOT->LoadMacro("makeCombinedData.C");\r
- makeData("data/scaled25Oct09/TOTALhistosscaled-LHC09b2-0.root",\r
- "data/scaled25Oct09/histosscaledLHC09b4AODc.root",\r
- "data/scaled25Oct09/histosWboson.root");\r
-\r
- gStyle->SetOptStat(0);\r
- drawAnnualYields();\r
- drawPtCutRates();\r
- drawHadEleRatios();\r
-\r
-}\r
-\r
-void drawAnnualYields() {\r
-\r
- TCanvas* crates = new TCanvas();\r
- crates->SetFillColor(0);\r
- crates->SetBorderMode(0);\r
- crates->SetBorderSize(2);\r
- crates->SetFrameBorderMode(0);\r
- crates->SetFrameBorderMode(0);\r
-\r
- crates->cd();\r
- gPad->SetLogy();\r
- allMC->SetXTitle("p_{T} (GeV/c)");\r
- allMC->SetTitle("MC electrons in Pb+Pb, 5.5 TeV");\r
- allMC->SetYTitle("Annual yield in EMCAL dN/dp_{T} (GeV/c)^{-1}");\r
- allMC->GetYaxis()->SetRangeUser(1,2.E6);\r
- allMC->GetXaxis()->SetRangeUser(10.,50.);\r
- allMC->Draw();\r
- bMC->Draw("same");\r
- cMC->Draw("same"); \r
- cbMC->Draw("same"); \r
- convMC->Draw("same"); \r
- dalMC->Draw("same"); \r
- wzMC->Draw("same");\r
-\r
- leg = new TLegend(0.6,0.6,0.9,0.9);\r
- leg->SetFillColor(0);\r
- leg->SetTextSize(leg->GetTextSize()*1.2);\r
- leg->AddEntry(allMC,"All MC electrons","l");\r
- leg->AddEntry(bMC,"Bottom e","l");\r
- leg->AddEntry(cMC,"Charm e","l");\r
- leg->AddEntry(cbMC,"B-->C e","l");\r
- leg->AddEntry(dalMC,"Dalitz e","l");\r
- leg->AddEntry(convMC,"Conversion e","l");\r
- leg->AddEntry(wzMC,"W Boson e","l");\r
- leg->Draw();\r
- crates->Print("MCRates_all.pdf");\r
-\r
- TCanvas* crates2 = new TCanvas();\r
- crates2->Divide(2,4);\r
- crates2->cd(1); gPad->SetLogy(); allMC->Draw();\r
- crates2->cd(2); gPad->SetLogy(); bMC->Draw();\r
- crates2->cd(3); gPad->SetLogy(); cMC->Draw();\r
- crates2->cd(4); gPad->SetLogy(); cbMC->Draw();\r
- crates2->cd(5); gPad->SetLogy(); convMC->Draw();\r
- crates2->cd(6); gPad->SetLogy(); dalMC->Draw();\r
- crates2->cd(7); gPad->SetLogy(); wzMC->Draw();\r
- crates2->cd(8); gPad->SetLogy(); mchad->Draw();\r
-\r
-}\r
-\r
-void drawPtCutRates() {\r
-\r
- TCanvas* cptcut = new TCanvas();\r
- cptcut->SetFillColor(0);\r
- cptcut->SetBorderMode(0);\r
- cptcut->SetBorderSize(2);\r
- cptcut->SetFrameBorderMode(0);\r
- cptcut->SetFrameBorderMode(0);\r
-\r
- cptcut->cd();\r
- gPad->SetLogy();\r
- TH1F* alleptcut = GetPtCutHisto(allMC);\r
- TH1F* beleptcut = GetPtCutHisto(bMC);\r
- TH1F* celeptcut = GetPtCutHisto(cMC);\r
- TH1F* cbeleptcut = GetPtCutHisto(cbMC);\r
- TH1F* dalitzptcut = GetPtCutHisto(dalMC);\r
- TH1F* convptcut = GetPtCutHisto(convMC);\r
- TH1F* wzptcut = GetPtCutHisto(wzMC);\r
- alleptcut->GetXaxis()->SetRangeUser(10,50);\r
- alleptcut->GetYaxis()->SetRangeUser(10,2.e6);\r
- alleptcut->SetXTitle("p_{T}^{cut} (GeV/c)");\r
- alleptcut->SetYTitle("Annual Yield in EMCAL for p_{T}>p_{T}^{cut}");\r
- alleptcut->SetTitle("MC electrons in Pb+Pb, 5.5 TeV");\r
- alleptcut->Draw();\r
- beleptcut->Draw("same");\r
- celeptcut->Draw("same");\r
- cbeleptcut->Draw("same");\r
- dalitzptcut->Draw("same");\r
- convptcut->Draw("same");\r
- wzptcut->Draw("same");\r
- leg->Draw();\r
- cptcut->Print("MCRates_ptcut_all.pdf");\r
-\r
-}\r
-\r
-void drawHadEleRatios() {\r
-\r
- TCanvas* ceh = new TCanvas();\r
- ceh->SetFillColor(0);\r
- ceh->SetBorderMode(0);\r
- ceh->SetBorderSize(2);\r
- ceh->SetFrameBorderMode(0);\r
- ceh->SetFrameBorderMode(0);\r
-\r
- ceh->cd();\r
- gPad->SetLogy();\r
- gStyle->SetOptStat(0);\r
- allheratio->Rebin(2); allheratio->Scale(1./2.);\r
- behratio->Rebin(2); behratio->Scale(1./2.);\r
- allheratio->SetLineWidth(2);\r
- allheratio->GetYaxis()->SetRangeUser(10,2e3);\r
- allheratio->GetXaxis()->SetRangeUser(10.,49.);\r
- behratio->GetXaxis()->SetRangeUser(10.,49.);\r
- allheratio->SetMarkerStyle(20);\r
- behratio->SetMarkerStyle(24);\r
- behratio->SetMarkerColor(1);\r
- allheratio->Draw();\r
- behratio->Draw("psame");\r
-\r
- TLegend *heleg = new TLegend(0.15,0.15,0.5,0.35);\r
- heleg->SetFillColor(0);\r
- heleg->SetTextSize(heleg->GetTextSize()*1.5);\r
- heleg->AddEntry(allheratio,"All electrons","l");\r
- heleg->AddEntry(behratio,"Bottom electrons","p");\r
- heleg->Draw();\r
- ceh->Print("MCRates_heratio.pdf");\r
-}\r
-\r
+++ /dev/null
-/////////////////////////////////////////////////\r
-//\r
-// Macro for plotting rates of identified Non-photonic electrons\r
-// for the EMCAL PPR\r
-//\r
-// J.L. Klay (Cal Poly)\r
-//\r
-/////////////////////////////////////////////////\r
-\r
-TLegend* leg;\r
-\r
-void plotNPERates(const char* which = "EMC") {\r
-\r
- gROOT->LoadMacro("makeCombinedData.C");\r
- makeData("data/scaled25Oct09/TOTALhistosscaled-LHC09b2-0.root",\r
- "data/scaled25Oct09/histosscaledLHC09b4AODc.root",\r
- "data/scaled25Oct09/histosWboson.root");\r
-\r
- //define common legend\r
- leg = new TLegend(0.5,0.6,0.9,0.9);\r
- leg->SetFillColor(0);\r
- leg->SetTextSize(leg->GetTextSize()*1.2);\r
- //leg->AddEntry(alltte,"All N-P e candidates","l");\r
- leg->AddEntry(sumtte,"All N-P electrons","l");\r
- leg->AddEntry(btte,"Bottom e","l");\r
- leg->AddEntry(ctte,"Charm e","l");\r
- leg->AddEntry(cbtte,"B-->C e","l");\r
- leg->AddEntry(daltte,"Dalitz e","l");\r
- leg->AddEntry(convtte,"Conversion e","l");\r
- leg->AddEntry(wztte,"W Boson e","l");\r
- // leg->AddEntry(htte,"Misidentified hadrons","l");\r
-\r
- gStyle->SetOptStat(0);\r
- // drawAnnualYields(which);\r
- // drawPtCutRates(which);\r
- //drawCompareTruth();\r
-\r
-}\r
-\r
-void drawAnnualYields(char* which = "EMC") {\r
-\r
- TCanvas* crates = new TCanvas();\r
- crates->SetFillColor(0);\r
- crates->SetBorderMode(0);\r
- crates->SetBorderSize(2);\r
- crates->SetFrameBorderMode(0);\r
- crates->SetFrameBorderMode(0);\r
-\r
- crates->cd();\r
- gPad->SetLogy();\r
-\r
- if(strcmp(which,"EMC")==0) {\r
- allemc->SetXTitle("p_{T} (GeV/c)");\r
- allemc->SetYTitle("Annual yield");\r
- allemc->SetTitle("Annual yield of non-phot. electron candidates (EMCAL pid)");\r
- allemc->Rebin(5); allemc->Scale(0.2);\r
-\r
- sumemc->SetXTitle("p_{T} (GeV/c)");\r
- sumemc->SetYTitle("Annual yield");\r
- sumemc->SetTitle("Annual yield of non-phot. electrons (EMCAL pid)");\r
- sumemc->Rebin(5); sumemc->Scale(0.2);\r
-\r
- bemc->Rebin(5); bemc->Scale(0.2);\r
- cemc->Rebin(5); cemc->Scale(0.2);\r
- cbemc->Rebin(5); cbemc->Scale(0.2);\r
- convemc->Rebin(5); convemc->Scale(0.2);\r
- dalemc->Rebin(5); dalemc->Scale(0.2);\r
- wzemc->Rebin(5); wzemc->Scale(0.2);\r
- hemc->Rebin(5); hemc->Scale(0.2);\r
-\r
- allemc->GetYaxis()->SetRangeUser(1.,2.e6);\r
- allemc->GetXaxis()->SetRangeUser(10.,49.);\r
- //allemc->Draw();\r
- sumemc->GetYaxis()->SetRangeUser(1.,2.e6);\r
- sumemc->GetXaxis()->SetRangeUser(10.,49.);\r
- sumemc->Draw();\r
- bemc->Draw("same");\r
- cemc->Draw("same");\r
- cbemc->Draw("same");\r
- convemc->Draw("same");\r
- dalemc->Draw("same");\r
- wzemc->Draw("same");\r
- // hemc->Draw("same");\r
- leg->Draw();\r
- crates->Print("NPERates_EMC_all.pdf");\r
- }\r
- if(strcmp(which,"TRK")==0) {\r
- alltrk->SetXTitle("p_{T} (GeV/c)");\r
- alltrk->SetYTitle("Annual yield");\r
- alltrk->SetTitle("Annual yield of non-phot. electron candidates (TPC+TRD pid)");\r
- alltrk->Rebin(5); alltrk->Scale(0.2);\r
- sumtrk->SetXTitle("p_{T} (GeV/c)");\r
- sumtrk->SetYTitle("Annual yield");\r
- sumtrk->SetTitle("Annual yield of non-phot. electrons (TPC+TRD pid)");\r
- sumtrk->Rebin(5); sumtrk->Scale(0.2);\r
-\r
- btrk->Rebin(5); btrk->Scale(0.2);\r
- ctrk->Rebin(5); ctrk->Scale(0.2);\r
- cbtrk->Rebin(5); cbtrk->Scale(0.2);\r
- convtrk->Rebin(5); convtrk->Scale(0.2);\r
- daltrk->Rebin(5); daltrk->Scale(0.2);\r
- wztrk->Rebin(5); wztrk->Scale(0.2);\r
- htrk->Rebin(5); htrk->Scale(0.2);\r
-\r
- alltrk->GetYaxis()->SetRangeUser(1.,6.e6);\r
- alltrk->GetXaxis()->SetRangeUser(10.,49.);\r
- // alltrk->Draw();\r
- sumtrk->GetYaxis()->SetRangeUser(1.,6.e6);\r
- sumtrk->GetXaxis()->SetRangeUser(10.,49.);\r
- sumtrk->Draw();\r
- btrk->Draw("same");\r
- ctrk->Draw("same");\r
- cbtrk->Draw("same");\r
- convtrk->Draw("same");\r
- daltrk->Draw("same");\r
- wztrk->Draw("same");\r
- // htrk->Draw("same");\r
- leg->Draw();\r
- crates->Print("NPERates_TRK_all.pdf");\r
- }\r
- if(strcmp(which,"TTE")==0) {\r
- alltte->SetXTitle("p_{T} (GeV/c)");\r
- alltte->SetYTitle("Annual yield");\r
- alltte->SetTitle("Annual yield of non-phot. electron candidates (Tracking+EMCAL pid)");\r
- alltte->Rebin(5); alltte->Scale(0.2);\r
- sumtte->SetXTitle("p_{T} (GeV/c)");\r
- sumtte->SetYTitle("Annual yield");\r
- sumtte->SetTitle("Annual yield of non-phot. electrons (Tracking+EMCAL pid)");\r
- sumtte->Rebin(5); sumtte->Scale(0.2);\r
-\r
- btte->Rebin(5); btte->Scale(0.2);\r
- ctte->Rebin(5); ctte->Scale(0.2);\r
- cbtte->Rebin(5); cbtte->Scale(0.2);\r
- convtte->Rebin(5); convtte->Scale(0.2);\r
- daltte->Rebin(5); daltte->Scale(0.2);\r
- wztte->Rebin(5); wztte->Scale(0.2);\r
- htte->Rebin(5); htte->Scale(0.2);\r
-\r
- alltte->GetYaxis()->SetRangeUser(1.,2.e6);\r
- alltte->GetXaxis()->SetRangeUser(10.,49.);\r
- // alltte->Draw();\r
- sumtte->GetYaxis()->SetRangeUser(1.,2.e6);\r
- sumtte->GetXaxis()->SetRangeUser(10.,49.);\r
- sumtte->Draw();\r
- btte->Draw("same");\r
- ctte->Draw("same");\r
- cbtte->Draw("same");\r
- convtte->Draw("same");\r
- daltte->Draw("same");\r
- wztte->Draw("same");\r
- htte->Draw("same");\r
- leg->Draw();\r
- crates->Print("NPERates_TTE_all.pdf");\r
- }\r
-\r
-}\r
-\r
-void drawPtCutRates(char* which = "EMC") {\r
-\r
- TCanvas* cptcut = new TCanvas();\r
- cptcut->SetFillColor(0);\r
- cptcut->SetBorderMode(0);\r
- cptcut->SetBorderSize(2);\r
- cptcut->SetFrameBorderMode(0);\r
- cptcut->SetFrameBorderMode(0);\r
- cptcut->cd();\r
- gPad->SetLogy();\r
- if(strcmp(which,"EMC")==0) {\r
- // TH1F* alleptcut = GetPtCutHisto(allemc);\r
- TH1F* alleptcut = GetPtCutHisto(sumemc);\r
- TH1F* beleptcut = GetPtCutHisto(bemc);\r
- TH1F* celeptcut = GetPtCutHisto(cemc);\r
- TH1F* cbeleptcut = GetPtCutHisto(cbemc);\r
- TH1F* dalitzptcut = GetPtCutHisto(dalemc);\r
- TH1F* convptcut = GetPtCutHisto(convemc);\r
- TH1F* wzptcut = GetPtCutHisto(wzemc);\r
- TH1F* misidptcut = GetPtCutHisto(hemc);\r
- alleptcut->SetXTitle("p_{T}^{cut} (GeV/c)");\r
- alleptcut->SetYTitle("Annual Yield in EMCAL for p_{T}>p_{T}^{cut}");\r
- alleptcut->SetTitle("Pb+Pb, 5.5 TeV reconstructed N-P electrons (EMCAL pid)");\r
- alleptcut->GetXaxis()->SetRangeUser(10.,49.);\r
- alleptcut->GetYaxis()->SetRangeUser(1,4.e5);\r
- alleptcut->Draw();\r
- beleptcut->Draw("same");\r
- celeptcut->Draw("same");\r
- cbeleptcut->Draw("same");\r
- dalitzptcut->Draw("same");\r
- convptcut->Draw("same");\r
- wzptcut->Draw("same");\r
- //misidptcut->Draw("same");\r
- leg->Draw();\r
- cptcut->Print("NPERates_EMC_ptcut_all.pdf");\r
- }\r
- if(strcmp(which,"TRK")==0) {\r
- // TH1F* alleptcut = GetPtCutHisto(alltrk);\r
- TH1F* alleptcut = GetPtCutHisto(sumtrk);\r
- TH1F* beleptcut = GetPtCutHisto(btrk);\r
- TH1F* celeptcut = GetPtCutHisto(ctrk);\r
- TH1F* cbeleptcut = GetPtCutHisto(cbtrk);\r
- TH1F* dalitzptcut = GetPtCutHisto(daltrk);\r
- TH1F* convptcut = GetPtCutHisto(convtrk);\r
- TH1F* wzptcut = GetPtCutHisto(wztrk);\r
- TH1F* misidptcut = GetPtCutHisto(htrk);\r
- alleptcut->SetXTitle("p_{T}^{cut} (GeV/c)");\r
- alleptcut->SetYTitle("Annual Yield in EMCAL for p_{T}>p_{T}^{cut}");\r
- alleptcut->SetTitle("Pb+Pb, 5.5 TeV reconstructed N-P electrons (TPC+TRD pid)");\r
- alleptcut->GetXaxis()->SetRangeUser(10.,49.);\r
- alleptcut->GetYaxis()->SetRangeUser(1,6.e6);\r
- alleptcut->Draw();\r
- beleptcut->Draw("same");\r
- celeptcut->Draw("same");\r
- cbeleptcut->Draw("same");\r
- dalitzptcut->Draw("same");\r
- convptcut->Draw("same");\r
- wzptcut->Draw("same");\r
- //misidptcut->Draw("same");\r
- leg->Draw();\r
- cptcut->Print("NPERates_TRK_ptcut_all.pdf");\r
- }\r
- if(strcmp(which,"TTE")==0) {\r
- // TH1F* alleptcut = GetPtCutHisto(alltte);\r
- TH1F* alleptcut = GetPtCutHisto(sumtte);\r
- TH1F* beleptcut = GetPtCutHisto(btte);\r
- TH1F* celeptcut = GetPtCutHisto(ctte);\r
- TH1F* cbeleptcut = GetPtCutHisto(cbtte);\r
- TH1F* dalitzptcut = GetPtCutHisto(daltte);\r
- TH1F* convptcut = GetPtCutHisto(convtte);\r
- TH1F* wzptcut = GetPtCutHisto(wztte);\r
- TH1F* misidptcut = GetPtCutHisto(htte);\r
- alleptcut->SetXTitle("p_{T}^{cut} (GeV/c)");\r
- alleptcut->SetYTitle("Annual Yield in EMCAL for p_{T}>p_{T}^{cut}");\r
- alleptcut->SetTitle("Pb+Pb, 5.5 TeV reconstructed N-P electrons (Tracking+EMCAL pid)");\r
- alleptcut->GetXaxis()->SetRangeUser(10.,49.);\r
- alleptcut->GetYaxis()->SetRangeUser(1,4.e5);\r
- alleptcut->Draw();\r
- beleptcut->Draw("same");\r
- celeptcut->Draw("same");\r
- cbeleptcut->Draw("same");\r
- dalitzptcut->Draw("same");\r
- convptcut->Draw("same");\r
- wzptcut->Draw("same");\r
- misidptcut->Draw("same");\r
-\r
- leg->Draw();\r
- cptcut->Print("NPERates_TTE_ptcut_all.pdf");\r
- }\r
-\r
-}\r
-\r
-TH1F* drawCompareTruth() {\r
-\r
- TH1F* mctruth = (TH1F*)belemc->Clone();\r
- mctruth->SetName("mctruth");\r
- mctruth->Add(celemc);\r
- mctruth->Add(candbmc);\r
- mctruth->Add(wzmc);\r
- mctruth->Rebin(2); mctruth->Scale(0.5);\r
-\r
- TFile* effic = new TFile("elec_eff.root");\r
- TH1F* heff = (TH1F*)effic->Get("h111");\r
-\r
- TH1F* hcorr = (TH1F*)sumHFemc->Clone();\r
- hcorr->SetName("hcorr");\r
- for(Int_t i = 1; i < heff->GetNbinsX(); i++) {\r
- Double_t pt = heff->GetBinCenter(i);\r
- Double_t eff = heff->GetBinContent(i); \r
- Double_t corr = 0.;\r
- if(eff > 0.) corr = hcorr->GetBinContent(i)/eff;\r
- hcorr->SetBinContent(i,corr);\r
- }\r
- hcorr->Rebin(2); hcorr->Scale(0.5);\r
- sumHFemc->Rebin(2); sumHFemc->Scale(0.5);\r
-\r
- Double_t efinal = 0.258;\r
- TGraphErrors* eerr = new TGraphErrors();\r
- eerr->SetName("emcErr");\r
- for(Int_t i = 1; i <= hcorr->GetNbinsX(); i++) {\r
- eerr->SetPoint(i-1,hcorr->GetBinCenter(i),hcorr->GetBinContent(i));\r
- eerr->SetPointError(i-1,0.,efinal*hcorr->GetBinContent(i));\r
- }\r
- eerr->SetFillColor(kRed-8);\r
-\r
- TCanvas* ctruth = new TCanvas();\r
- ctruth->SetFillColor(0);\r
- ctruth->SetBorderMode(0);\r
- ctruth->SetBorderSize(2);\r
- ctruth->SetFrameBorderMode(0);\r
- ctruth->SetFrameBorderMode(0);\r
-\r
- ctruth->cd();\r
- gPad->SetLogy();\r
- mctruth->SetTitle("Comparison of MC and reco HF+W electrons");\r
- mctruth->SetMarkerColor(kBlack); mctruth->SetLineColor(kBlack);\r
- mctruth->SetYTitle("Annual yield in EMCAL dN/dp_{T} (GeV/c)^{-1}");\r
- mctruth->SetXTitle("p_{T} (GeV/c)");\r
- mctruth->GetXaxis()->SetRangeUser(10.,49.);\r
- mctruth->Draw();\r
- eerr->Draw("3same");\r
- hcorr->SetMarkerColor(kRed);\r
- hcorr->SetLineColor(kRed);\r
- hcorr->Draw("same");\r
- sumHFemc->SetMarkerColor(kBlue);\r
- sumHFemc->SetLineColor(kBlue);\r
- sumHFemc->Draw("same");\r
- TLegend *legy = new TLegend(0.3,0.7,0.9,0.9);\r
- legy->SetFillColor(0);\r
- legy->AddEntry(mctruth,"MC HF+W electrons","l");\r
- legy->AddEntry(sumHFemc,"Rec (EMCAL) HF+W electrons","l");\r
- legy->AddEntry(hcorr,"Eff. corrected reco HF+W electrons","l");\r
- legy->AddEntry(eerr,"Systematic uncertainty","f");\r
- legy->Draw();\r
-\r
- ctruth->Print("NPERates_TruthComparison.pdf");\r
-\r
-}\r
+++ /dev/null
-/////////////////////////////////////////////////\r
-//\r
-// Macro for plotting rates of identified Non-photonic electrons\r
-// for the EMCAL PPR\r
-//\r
-// J.L. Klay (Cal Poly)\r
-//\r
-/////////////////////////////////////////////////\r
-\r
-void plotPIDCompare(char* which = "EMC") {\r
-\r
- gROOT->LoadMacro("makeCombinedData.C");\r
- makeData("data/scaled25Oct09/TOTALhistosscaled-LHC09b2-0.root",\r
- "data/scaled25Oct09/histosscaledLHC09b4AODc.root",\r
- "data/scaled25Oct09/histosWboson.root");\r
-\r
- //define common legend\r
- TLegend* leg = new TLegend(0.5,0.6,0.9,0.9);\r
- leg->SetFillColor(0);\r
- leg->SetTextSize(leg->GetTextSize()*1.2);\r
- leg->AddEntry(alltte,"All N-P e candidates","l");\r
- //leg->AddEntry(sumtte,"All N-P electrons","l");\r
- leg->AddEntry(btte,"Bottom e","l");\r
- leg->AddEntry(ctte,"Charm e","l");\r
- leg->AddEntry(cbtte,"B-->C e","l");\r
- leg->AddEntry(daltte,"Dalitz e","l");\r
- leg->AddEntry(convtte,"Conversion e","l");\r
- leg->AddEntry(wztte,"W Boson e","l");\r
- //leg->AddEntry(htte,"Mis-identified hadrons","l");\r
-\r
- gStyle->SetOptStat(0);\r
- \r
- TCanvas* crates = new TCanvas("crates","",0,0,1600,600);\r
- crates->SetFillColor(0);\r
- crates->SetBorderMode(0);\r
- crates->SetBorderSize(2);\r
- crates->SetFrameBorderMode(0);\r
- crates->SetFrameBorderMode(0);\r
- crates->Divide(2,1);\r
- crates->cd(1);\r
- gPad->SetLogy();\r
- alltrk->SetXTitle("p_{T} (GeV/c)");\r
- alltrk->SetYTitle("Annual yield in EMCAL dN/dp_{T} (GeV/c)^{-1}");\r
- alltrk->SetTitle("PID comparison: Tracking only vs. EMCAL only");\r
- alltrk->Rebin(5); alltrk->Scale(0.2);\r
- alltrk->GetYaxis()->SetRangeUser(10.,alltrk->GetMaximum()*2.);\r
- alltrk->GetXaxis()->SetRangeUser(10.,49.);\r
- alltrk->Draw();\r
- htrk->Rebin(5); htrk->Scale(0.2);\r
- htrk->Draw("same");\r
- TH1F* tempallemc = (TH1F*)allemc->Clone();\r
- tempallemc->SetNameTitle("tempallemc","tempallemc");\r
- tempallemc->SetLineColor(kBlue);\r
- tempallemc->Rebin(5); tempallemc->Scale(0.2);\r
- tempallemc->Draw("same");\r
- \r
- TH1F* temphemc = (TH1F*)hemc->Clone();\r
- temphemc->SetNameTitle("temphemc","temphemc");\r
- temphemc->SetLineColor(kOrange-3);\r
- temphemc->Rebin(5); temphemc->Scale(0.2);\r
- temphemc->Draw("same");\r
-\r
- TLegend* leg2 = new TLegend(0.35,0.6,0.9,0.9);\r
- leg2->SetFillColor(0);\r
- leg2->SetTextSize(leg->GetTextSize()*1.2);\r
- leg2->AddEntry(alltrk,"Electron Candidates (Tracking PID only)","l");\r
- leg2->AddEntry(htrk,"Hadron Contamination (Tracking PID only)","l");\r
- leg2->AddEntry(tempallemc,"Electron Candidates (EMCAL PID)","l");\r
- leg2->AddEntry(temphemc,"Hadron Contamination (EMCAL PID)","l");\r
- leg2->Draw();\r
-\r
- crates->cd(2);\r
- gPad->SetLogy();\r
- TH1F* subtrk = (TH1F*)alltrk->Clone();\r
- for(Int_t i = 1; i <= alltrk->GetNbinsX(); i++) {\r
- Double_t diff = alltrk->GetBinContent(i) - htrk->GetBinContent(i);\r
- Double_t unc = 0.;\r
- if(diff < 0) diff = 0.;\r
- if(diff > 0) \r
- unc = diff/TMath::Sqrt(diff+2.*htrk->GetBinContent(i));\r
- printf("<%d> Cand %d, Contam %d, diff %d, unc %d\n",i,alltrk->GetBinContent(i),htrk->GetBinContent(i), diff, unc);\r
- subtrk->SetBinContent(i,diff);\r
- subtrk->SetBinError(i,unc);\r
- }\r
- subtrk->SetYTitle("Annual yield in EMCAL dN/dp_{T} (GeV/c)^{-1}");\r
- subtrk->SetLineColor(kRed);\r
- subtrk->SetMarkerStyle(20); subtrk->SetMarkerColor(kRed);\r
- subtrk->Draw();\r
-\r
- TH1F* subemc = (TH1F*)tempallemc->Clone();\r
- for(Int_t i = 1; i <= tempallemc->GetNbinsX(); i++) {\r
- Double_t diff = tempallemc->GetBinContent(i) - temphemc->GetBinContent(i);\r
- Double_t unc = 0.;\r
- if(diff < 0.) diff = subemc->GetBinContent(i-1);\r
- if(diff > 0)\r
- unc = diff/TMath::Sqrt(diff+2.*temphemc->GetBinContent(i));\r
- // printf("Cand %d, Contam %d, diff %d, unc %d\n",tempallemc->GetBinContent(i),temphemc->GetBinContent(i), diff, unc);\r
- subemc->SetBinContent(i,diff);\r
- subemc->SetBinError(i,unc);\r
- }\r
-\r
- TFile* effic = new TFile("elec_eff.root");\r
- TH1F* heff = (TH1F*)effic->Get("h111");\r
- heff->Rebin(5); heff->Scale(0.2);\r
- TH1F* hcorr = (TH1F*)subemc->Clone();\r
- hcorr->SetName("hcorr");\r
- for(Int_t i = 1; i < heff->GetNbinsX(); i++) {\r
- Double_t pt = heff->GetBinCenter(i);\r
- Double_t eff = heff->GetBinContent(i);\r
- Double_t corr = 0.;\r
- if(eff > 0.) corr = hcorr->GetBinContent(i)/eff;\r
- hcorr->SetBinContent(i,corr);\r
- }\r
-\r
- Double_t efinal = 0.4;\r
- TGraphErrors* eerr = new TGraphErrors();\r
- eerr->SetName("emcErr");\r
- int count=0;\r
- for(Int_t i = 1; i <= hcorr->GetNbinsX(); i++) {\r
- if (hcorr->GetBinCenter(i) <10.0) continue;\r
- if (hcorr->GetBinCenter(i) >50.0) break;\r
- cout <<"bin:"<< i << ", bin-center:"<< hcorr->GetBinCenter(i)<< endl;\r
- eerr->SetPoint(count,hcorr->GetBinCenter(i),hcorr->GetBinContent(i));\r
- eerr->SetPointError(count,0.,efinal*hcorr->GetBinContent(i));\r
- // if(hcorr->GetBinCenter(i) <20.0) eerr->SetPointError(count,0.,1.5*efinal*hcorr->GetBinContent(i));\r
- count++;\r
- }\r
- eerr->SetFillColor(kRed-8);\r
- eerr->Draw("3same");\r
- subtrk->Draw("same");\r
- hcorr->SetMarkerStyle(20); hcorr->SetMarkerColor(kBlue);\r
- hcorr->SetLineColor(kBlue);\r
- hcorr->Draw("same");\r
- \r
- allMC->Draw("same");\r
- mchad->Draw("same");\r
-\r
- TLegend *legx = new TLegend(0.3,0.7,0.9,0.9);\r
- legx->SetFillColor(0);\r
- legx->AddEntry(subtrk,"Signal (candidates - contam) (Tracking PID only)","pl");\r
- legx->AddEntry(hcorr,"Eff. corrected signal (EMCAL PID)","pl");\r
- legx->AddEntry(eerr,"Systematic uncertainty","f");\r
- legx->AddEntry(allMC,"MC Electrons","l");\r
- legx->AddEntry(mchad,"MC Hadrons","l");\r
- legx->Draw();\r
-\r
- \r
- TLatex* latex = new TLatex(0.5,0.6,"Unc = #frac{S}{#sqrt{S+2B}}");\r
- latex->SetNDC();\r
- // latex->Draw();\r
-\r
- crates->Print("NPERates_PIDCompare_all.pdf");\r
-\r
- TCanvas* ccomp = new TCanvas("ccomp","",0,0,600,800);\r
- ccomp->SetFillColor(0);\r
- ccomp->SetBorderMode(0);\r
- ccomp->SetBorderSize(2);\r
- ccomp->SetFrameBorderMode(0);\r
- ccomp->SetFrameBorderMode(0);\r
- TPad* upperPad = new TPad("upperPad", "upperPad", \r
- .005, .2525, .995, .995);\r
- upperPad->SetFillColor(0);\r
- upperPad->SetBorderMode(0);\r
- upperPad->SetBorderSize(2);\r
- upperPad->SetFrameBorderMode(0);\r
- upperPad->SetFrameBorderMode(0);\r
- TPad* lowerPad = new TPad("lowerPad", "lowerPad", \r
- .005, .005, .995, .2525);\r
- lowerPad->SetFillColor(0);\r
- lowerPad->SetFillColor(0);\r
- lowerPad->SetBorderMode(0);\r
- lowerPad->SetBorderSize(2);\r
- lowerPad->SetFrameBorderMode(0);\r
- lowerPad->SetFrameBorderMode(0);\r
- upperPad->Draw(); \r
- lowerPad->Draw(); \r
- upperPad->cd();\r
- gPad->SetLogy();\r
- hcorr->SetTitle("Efficiency corrected non-photonic electrons");\r
- hcorr->SetYTitle("Annual yield in EMCAL dN/dp_{T} (GeV/c)^{-1}");\r
- hcorr->SetXTitle("p_{T} (GeV/c)");\r
- hcorr->GetXaxis()->SetRangeUser(10.,49.);\r
- hcorr->GetYaxis()->SetTitleOffset(1.2);\r
- hcorr->Draw();\r
- eerr->Draw("3same");\r
- hcorr->Draw("same");\r
- allMC->Draw("same");\r
- TLegend *myleg = new TLegend(0.3,0.7,0.9,0.9);\r
- myleg->SetFillColor(0);\r
- myleg->AddEntry(allMC,"MC Electrons","l");\r
- myleg->AddEntry(hcorr,"Eff. corrected EMCAL N-P electrons","pl");\r
- myleg->AddEntry(eerr,"Systematic uncertainty","f");\r
- myleg->Draw();\r
- lowerPad->cd();\r
- TH1F* hmc = (TH1F*)allMC->Clone(); hmc->SetName("hmc");\r
- hmc->Rebin(5); hmc->Scale(0.2);\r
- TH1F* hmcratio = (TH1F*)hcorr->Clone(); hmcratio->SetName("hmcratio");\r
- hmcratio->Divide(hmc);\r
- hmcratio->SetTitle();\r
- hmcratio->SetYTitle("MC/Data");\r
- hmcratio->GetYaxis()->SetRangeUser(0.,2.);\r
- hmcratio->GetYaxis()->SetNdivisions(505);\r
- hmcratio->GetYaxis()->SetLabelSize(2.*hmcratio->GetYaxis()->GetLabelSize());\r
- hmcratio->GetYaxis()->SetTitleSize(3.*hmcratio->GetYaxis()->GetTitleSize());\r
- hmcratio->GetYaxis()->SetTitleOffset(0.3);\r
- hmcratio->SetXTitle("");\r
- hmcratio->GetXaxis()->SetLabelSize(2.*hmcratio->GetXaxis()->GetLabelSize());\r
- hmcratio->Draw();\r
- TGraphErrors* rerr = new TGraphErrors();\r
- rerr->SetName("ratioErr");\r
- count=0;\r
- for(Int_t i = 1; i <= hmcratio->GetNbinsX(); i++) {\r
- if (hmcratio->GetBinCenter(i) <10.0) continue;\r
- if (hmcratio->GetBinCenter(i) >50.0) break;\r
- rerr->SetPoint(count,hmcratio->GetBinCenter(i),hmcratio->GetBinContent(i));\r
- rerr->SetPointError(count,0.,efinal*hmcratio->GetBinContent(i));\r
- //if(hmcratio->GetBinCenter(i) <20.0) rerr->SetPointError(count,0.,1.5*efinal*hmcratio->GetBinContent(i));\r
- count++;\r
- }\r
- rerr->SetFillColor(kRed-8);\r
- rerr->Draw("3same");\r
- hmcratio->Draw("same");\r
- TLine* line = new TLine(10.,1.,50.,1.);\r
- line->SetLineStyle(2);\r
- line->Draw();\r
- ccomp->Print("EMCAL_NPE.pdf");\r
- crates->cd();\r
-}\r
-\r
+++ /dev/null
-//#define VERBOSEARGS
-{
- //Process environmental variables from command line:
- char anainputdata[1024];
- char smode[1024];
- char sconfig1[1024];
- char sconfig2[1024];
- char sconfig3[1024];
- char sevent[1024];
-
- sprintf(anainputdata,"");
- sprintf(smode,"");
- sprintf(sconfig1,"");
- sprintf(sconfig2,"");
- sprintf(sconfig3,"");
- sprintf(sevent,"");
-
-
- for (int i=0; i< gApplication->Argc();i++){
-#ifdef VERBOSEARGS
- printf("Arg %d: %s\n",i,gApplication->Argv(i));
-#endif
-
- if (!(strcmp(gApplication->Argv(i),"--anaInputData")))
- sprintf(anainputdata,"%s",gApplication->Argv(i+1));
- if (!(strcmp(gApplication->Argv(i),"--mode")))
- sprintf(smode,"%s",gApplication->Argv(i+1));
- if (!(strcmp(gApplication->Argv(i),"--config1")))
- sprintf(sconfig1,"%s",gApplication->Argv(i+1));
- if (!(strcmp(gApplication->Argv(i),"--config2")))
- sprintf(sconfig2,"%s",gApplication->Argv(i+1));
- if (!(strcmp(gApplication->Argv(i),"--config3")))
- sprintf(sconfig3,"%s",gApplication->Argv(i+1));
- if (!(strcmp(gApplication->Argv(i),"--sevent")))
- sprintf(sevent,"%s",gApplication->Argv(i+1));
- }
-
- gSystem->Setenv("anaInputData", anainputdata);
- gSystem->Setenv("MODE" , smode);
- gSystem->Setenv("CONFIG1", sconfig1);
- gSystem->Setenv("CONFIG2", sconfig2);
- gSystem->Setenv("CONFIG3", sconfig3);
- gSystem->Setenv("SEVENT" , sevent);
-
-
- printf("Run.C: Variables: anaInputData %s, mode %s, config1 %s, config2 %s, config3 %s, first event %s\n", anainputdata,smode,sconfig1,sconfig2,sconfig3,sevent);
-
- gSystem->Exec("root -b -q -l anaJete.C > anaJete.log 2>&1");
-
-}
+++ /dev/null
-#!/bin/bash
-##################################################
-validateout=`dirname $0`
-validatetime=`date`
-validated="0";
-error=0
-if [ -z $validateout ]
-then
- validateout="."
-fi
-
-cd $validateout;
-validateworkdir=`pwd`;
-
-echo "*******************************************************" >> stdout
-echo "* Automatically generated validation script *" >> stdout
-
-echo "* Time: $validatetime " >> stdout
-echo "* Dir: $validateout" >> stdout
-echo "* Workdir: $validateworkdir" >> stdout
-echo "* ----------------------------------------------------*" >> stdout
-ls -la ./ >> stdout
-echo "* ----------------------------------------------------*" >> stdout
-
-##################################################
-if ! [ -f histos.root ] ; then
- error=1
- echo "Output file(s) not found. Job FAILED !" >> stdout
- echo "Output file(s) not found. Job FAILED !" >> stderr
-fi
-if [ $error = 0 ] ; then
- echo "* ---------------- Job Validated ------------------*" >> stdout
-fi
-echo "* ----------------------------------------------------*" >> stdout
-echo "*******************************************************" >> stdout
-cd -
-exit $error