X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=MUON%2FMUONmassPlot_ESD.C;h=7237e48b3e306f26052ab2a1c8e58518b509f365;hb=e1f6be442d709908d02410c6c8a4a56443b6d618;hp=f2e93bda73dea4ef7871a7c108e37b5d4301853b;hpb=8b4a9b898803945c38c09ed9f603ad74359ef44f;p=u%2Fmrichter%2FAliRoot.git diff --git a/MUON/MUONmassPlot_ESD.C b/MUON/MUONmassPlot_ESD.C index f2e93bda73d..7237e48b3e3 100644 --- a/MUON/MUONmassPlot_ESD.C +++ b/MUON/MUONmassPlot_ESD.C @@ -14,60 +14,55 @@ #include // STEER includes -#include "AliRun.h" #include "AliLog.h" -#include "AliRunLoader.h" -#include "AliHeader.h" -#include "AliLoader.h" -#include "AliStack.h" #include "AliMagFMaps.h" -#include "AliESD.h" +#include "AliESDEvent.h" +#include "AliESDVertex.h" #include "AliTracker.h" +#include "AliESDMuonTrack.h" // MUON includes #include "AliMUONTrackParam.h" #include "AliMUONTrackExtrap.h" -#include "AliESDMuonTrack.h" +#include "AliMUONESDInterface.h" #endif -// -// Macro MUONmassPlot.C for ESD -// Ch. Finck, Subatech, April. 2004 -// - -// macro to make invariant mass plots -// for combinations of 2 muons with opposite charges, -// from root file "MUON.tracks.root" containing the result of track reconstruction. -// Histograms are stored on the "MUONmassPlot.root" file. -// introducing TLorentzVector for parameter calculations (Pt, P,rap,etc...) -// using Invariant Mass for rapidity. - -// Arguments: -// ExtrapToVertex (default -1) -// <0: no extrapolation; -// =0: extrapolation to (0,0,0); -// >0: extrapolation to ESDVertex if available, else to (0,0,0) -// FirstEvent (default 0) -// LastEvent (default 0) -// ResType (default 553) -// 553 for Upsilon, anything else for J/Psi -// Chi2Cut (default 100) -// to keep only tracks with chi2 per d.o.f. < Chi2Cut -// PtCutMin (default 1) -// to keep only tracks with transverse momentum > PtCutMin -// PtCutMax (default 10000) -// to keep only tracks with transverse momentum < PtCutMax -// massMin (default 9.17 for Upsilon) -// & massMax (default 9.77 for Upsilon) -// to calculate the reconstruction efficiency for resonances with invariant mass -// massMin < mass < massMax. - -// Add parameters and histograms for analysis - -Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = -1, char* geoFilename = "geometry.root", + +/// \ingroup macros +/// \file MUONmassPlot_ESD.C +/// \brief Macro MUONefficiency.C for ESD +/// +/// \author Ch. Finck, Subatech, April. 2004 +/// +/// +/// Macro to make invariant mass plots +/// for combinations of 2 muons with opposite charges, +/// from root file "MUON.tracks.root" containing the result of track reconstruction. +/// Histograms are stored on the "MUONmassPlot.root" file. +/// introducing TLorentzVector for parameter calculations (Pt, P,rap,etc...) +/// using Invariant Mass for rapidity. +/// +/// Add parameters and histograms for analysis + +Bool_t MUONmassPlot(Int_t ExtrapToVertex = -1, char* geoFilename = "geometry.root", Int_t FirstEvent = 0, Int_t LastEvent = 10000, char* esdFileName = "AliESDs.root", Int_t ResType = 553, Float_t Chi2Cut = 100., Float_t PtCutMin = 1., Float_t PtCutMax = 10000., Float_t massMin = 9.17,Float_t massMax = 9.77) { +/// \param ExtrapToVertex (default -1) +/// - <0: no extrapolation; +/// - =0: extrapolation to (0,0,0); +/// - >0: extrapolation to ESDVertex if available, else to (0,0,0) +/// \param FirstEvent (default 0) +/// \param LastEvent (default 0) +/// \param ResType 553 for Upsilon, anything else for J/Psi (default 553) +/// \param Chi2Cut to keep only tracks with chi2 per d.o.f. < Chi2Cut (default 100) +/// \param PtCutMin to keep only tracks with transverse momentum > PtCutMin (default 1) +/// \param PtCutMax to keep only tracks with transverse momentum < PtCutMax (default 10000) +/// \param massMin (default 9.17 for Upsilon) +/// \param massMax (default 9.77 for Upsilon); +/// to calculate the reconstruction efficiency for resonances with invariant mass +/// massMin < mass < massMax. + cout << "MUONmassPlot " << endl; cout << "FirstEvent " << FirstEvent << endl; cout << "LastEvent " << LastEvent << endl; @@ -128,6 +123,8 @@ Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = - Double_t fZVertex=0; Double_t fYVertex=0; Double_t fXVertex=0; + Double_t errXVtx=0; + Double_t errYVtx=0; TLorentzVector fV1, fV2, fVtot; @@ -135,7 +132,7 @@ Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = - if (!gGeoManager) { TGeoManager::Import(geoFilename); if (!gGeoManager) { - Error("MUONmass_ESD", "getting geometry from file %s failed", filename); + Error("MUONmass_ESD", "getting geometry from file %s failed", geoFilename); return kFALSE; } } @@ -146,20 +143,6 @@ Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = - AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG); AliTracker::SetFieldMap(field, kFALSE); - // open run loader and load gAlice, kinematics and header - AliRunLoader* runLoader = AliRunLoader::Open(filename); - if (!runLoader) { - Error("MUONmass_ESD", "getting run loader from file %s failed", filename); - return kFALSE; - } -/* - runLoader->LoadgAlice(); - if (!gAlice) { - Error("MUONmass_ESD", "no galice object found"); - return kFALSE; - } -*/ - // open the ESD file TFile* esdFile = TFile::Open(esdFileName); if (!esdFile || !esdFile->IsOpen()) { @@ -167,7 +150,7 @@ Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = - return kFALSE; } - AliESD* esd = new AliESD(); + AliESDEvent* esd = new AliESDEvent(); TTree* tree = (TTree*) esdFile->Get("esdTree"); if (!tree) { Error("CheckESD", "no ESD tree found"); @@ -176,18 +159,13 @@ Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = - // tree->SetBranchAddress("ESD", &esd); esd->ReadFromTree(tree); - - runLoader->LoadHeader(); - nevents = runLoader->GetNumberOfEvents(); + nevents = (Int_t)tree->GetEntries(); AliMUONTrackParam trackParam; // Loop over events for (Int_t iEvent = FirstEvent; iEvent <= TMath::Min(LastEvent, nevents - 1); iEvent++) { - // get current event - runLoader->GetEvent(iEvent); - // get the event summary data tree->GetEvent(iEvent); if (!esd) { @@ -201,6 +179,8 @@ Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = - fZVertex = Vertex->GetZv(); fYVertex = Vertex->GetYv(); fXVertex = Vertex->GetXv(); + errXVtx = Vertex->GetXRes(); + errYVtx = Vertex->GetYRes(); } hPrimaryVertex->Fill(fZVertex); @@ -214,17 +194,20 @@ Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = - // loop over all reconstructed tracks (also first track of combination) for (Int_t iTrack = 0; iTrack < nTracks; iTrack++) { + // skip ghosts + if (!esd->GetMuonTrack(iTrack)->ContainTrackerData()) continue; + AliESDMuonTrack* muonTrack = new AliESDMuonTrack(*(esd->GetMuonTrack(iTrack))); // extrapolate to vertex if required and available if (ExtrapToVertex > 0 && Vertex->GetNContributors()) { - trackParam.GetParamFromUncorrected(*muonTrack); - AliMUONTrackExtrap::ExtrapToVertex(&trackParam, fXVertex, fYVertex, fZVertex); - trackParam.SetParamFor(*muonTrack); // put the new parameters in this copy of AliESDMuonTrack + AliMUONESDInterface::GetParamAtFirstCluster(*muonTrack, trackParam); + AliMUONTrackExtrap::ExtrapToVertex(&trackParam, fXVertex, fYVertex, fZVertex, errXVtx, errYVtx); + AliMUONESDInterface::SetParamAtVertex(trackParam, *muonTrack); // put the new parameters in this copy of AliESDMuonTrack } else if ((ExtrapToVertex > 0 && !Vertex->GetNContributors()) || ExtrapToVertex == 0){ - trackParam.GetParamFromUncorrected(*muonTrack); - AliMUONTrackExtrap::ExtrapToVertex(&trackParam, 0., 0., 0.); - trackParam.SetParamFor(*muonTrack); // put the new parameters in this copy of AliESDMuonTrack + AliMUONESDInterface::GetParamAtFirstCluster(*muonTrack, trackParam); + AliMUONTrackExtrap::ExtrapToVertex(&trackParam, 0., 0., 0., 0., 0.); + AliMUONESDInterface::SetParamAtVertex(trackParam, *muonTrack); // put the new parameters in this copy of AliESDMuonTrack } fCharge1 = Int_t(TMath::Sign(1.,muonTrack->GetInverseBendingMomentum())); @@ -267,17 +250,20 @@ Bool_t MUONmassPlot(char* filename = "galice_sim.root", Int_t ExtrapToVertex = - // loop over second track of combination for (Int_t iTrack2 = iTrack + 1; iTrack2 < nTracks; iTrack2++) { + // skip ghosts + if (!esd->GetMuonTrack(iTrack2)->ContainTrackerData()) continue; + AliESDMuonTrack* muonTrack2 = new AliESDMuonTrack(*(esd->GetMuonTrack(iTrack2))); // extrapolate to vertex if required and available if (ExtrapToVertex > 0 && Vertex->GetNContributors()) { - trackParam.GetParamFromUncorrected(*muonTrack2); - AliMUONTrackExtrap::ExtrapToVertex(&trackParam, fXVertex, fYVertex, fZVertex); - trackParam.SetParamFor(*muonTrack2); // put the new parameters in this copy of AliESDMuonTrack + AliMUONESDInterface::GetParamAtFirstCluster(*muonTrack2, trackParam); + AliMUONTrackExtrap::ExtrapToVertex(&trackParam, fXVertex, fYVertex, fZVertex, errXVtx, errYVtx); + AliMUONESDInterface::SetParamAtVertex(trackParam, *muonTrack2); // put the new parameters in this copy of AliESDMuonTrack } else if ((ExtrapToVertex > 0 && !Vertex->GetNContributors()) || ExtrapToVertex == 0){ - trackParam.GetParamFromUncorrected(*muonTrack2); - AliMUONTrackExtrap::ExtrapToVertex(&trackParam, 0., 0., 0.); - trackParam.SetParamFor(*muonTrack2); // put the new parameters in this copy of AliESDMuonTrack + AliMUONESDInterface::GetParamAtFirstCluster(*muonTrack2, trackParam); + AliMUONTrackExtrap::ExtrapToVertex(&trackParam, 0., 0., 0., 0., 0.); + AliMUONESDInterface::SetParamAtVertex(trackParam, *muonTrack2); // put the new parameters in this copy of AliESDMuonTrack } fCharge2 = Int_t(TMath::Sign(1.,muonTrack2->GetInverseBendingMomentum()));