1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
19 // Macro for checking trigger integrated efficiency for dimuons.
20 // The efficiency is calculated with respect to the 3/4 coincidence.
21 // Author: Fabien Guerin, LPC Clermont-Ferrand, Jan. 2006
24 #if !defined(__CINT__) || defined(__MAKECINT__)
27 #include "TClonesArray.h"
33 #include "TParticle.h"
37 #include "AliRunLoader.h"
38 #include "AliHeader.h"
39 #include "AliLoader.h"
44 #include "AliMUONHit.h"
45 #include "AliMUONConstants.h"
46 #include "AliMUONDigit.h"
47 #include "AliMUONGlobalTrigger.h"
48 #include "AliMUONTrack.h"
50 #include "AliMUONDataInterface.h"
51 #include "AliMUONMCDataInterface.h"
52 #include "AliMUONVHitStore.h"
53 #include "AliMUONVTriggerStore.h"
60 void MUONTriggerEfficiency(const char* filenameSim="galice_sim.root",
61 const char* filenameRec="galice.root",
62 Bool_t readFromRP = 0)
67 AliMUONMCDataInterface diSim(filenameSim);
68 AliMUONDataInterface diRec(filenameRec);
72 cerr << "Cannot access " << filenameSim << endl;
78 cerr << "Cannot access " << filenameRec << endl;
82 FILE* fdat = fopen("MUONTriggerEfficiency.out","w");
85 cerr << "Cannot create output file" << endl;
89 Int_t coincmuon,muonlpt,muonhpt;
90 Int_t CoincMuPlus,CoincMuMinus;
95 Int_t nevents = diSim.NumberOfEvents();
97 for (Int_t ievent=0; ievent<nevents; ++ievent)
102 if (ievent%1000==0) printf("\t Event = %d\n",ievent);
107 for (Int_t j=0; j<2; j++)
109 for (Int_t jj=0; jj<4; jj++)
115 Int_t ntracks = (Int_t) diSim.NumberOfTracks(ievent);
117 for ( Int_t itrack=0; itrack<ntracks; ++itrack )
119 AliMUONVHitStore* hitStore = diSim.HitStore(ievent,itrack);
121 TIter next(hitStore->CreateIterator());
123 while ( ( mHit = static_cast<AliMUONHit*>(next()) ) )
125 Int_t Nch = mHit->Chamber();
126 Int_t hittrack = mHit->Track();
127 Float_t IdPart = mHit->Particle();
129 for (Int_t j=0;j<4;j++)
132 if (hittrack==1 || hittrack==2)
134 if (Nch==kch && IdPart==-13 && NbHits[0][j]==0) NbHits[0][j]++;
135 if (Nch==kch && IdPart==13 && NbHits[1][j]==0) NbHits[1][j]++;
142 Int_t SumNbHits=NbHits[0][0]+NbHits[0][1]+NbHits[0][2]+NbHits[0][3];
143 if (SumNbHits==3 || SumNbHits==4) CoincMuPlus=1;
145 SumNbHits=NbHits[1][0]+NbHits[1][1]+NbHits[1][2]+NbHits[1][3];
146 if (SumNbHits==3 || SumNbHits==4) CoincMuMinus=1;
148 if (CoincMuPlus==1 && CoincMuMinus==1) coincmuon++;
151 if ( readFromRP ) tree = "R";
153 AliMUONVTriggerStore* triggerStore = diRec.TriggerStore(ievent,tree.Data());
155 AliMUONGlobalTrigger* gloTrg = triggerStore->Global();
157 if (gloTrg->PairUnlikeLpt()>=1) muonlpt++;
158 if (gloTrg->PairUnlikeHpt()>=1) muonhpt++;
160 } // end loop on event
162 // calculate efficiency with as a ref. at least 3/4 planes fired
163 Float_t efficiencylpt,efficiencyhpt;
164 Double_t coincmu,lptmu,hptmu;
166 coincmu=Double_t(coincmuon);
167 lptmu=Double_t(muonlpt);
168 hptmu=Double_t(muonhpt);
172 fprintf(fdat," Number of events = %d \n",nevents);
173 fprintf(fdat," Number of events with 3/4 coinc = %d \n",coincmuon);
174 fprintf(fdat," Nomber of dimuons with 3/4 coinc Lpt cut = %d \n",muonlpt);
175 fprintf(fdat," Number of dimuons with 3/4 coinc Hpt cut = %d \n",muonhpt);
178 efficiencylpt=lptmu/coincmu;
179 error=efficiencylpt*TMath::Sqrt((lptmu+coincmu)/(lptmu*coincmu));
180 fprintf(fdat," Efficiency Lpt cut = %4.4f +/- %4.4f\n",efficiencylpt,error);
182 efficiencyhpt=hptmu/coincmu;
183 error=efficiencyhpt*TMath::Sqrt((hptmu+coincmu)/(hptmu*coincmu));
184 fprintf(fdat," Efficiency Hpt cut = %4.4f +/- %4.4f\n",efficiencyhpt,error);