WriteRecPoints(ievent) ;
- if(strstr(option,"deb")) PrintRecPoints(option) ;
- }
+ if(strstr(option,"deb"))
+ PrintRecPoints(option) ;
+
+ //increment the total number of digits per run
+ fRecPointsInRun += gime->EmcRecPoints()->GetEntriesFast() ;
+ fRecPointsInRun += gime->CpvRecPoints()->GetEntriesFast() ;
+ }
if(strstr(option,"tim")){
gBenchmark->Stop("PHOSClusterizer");
if(strstr(option,"deb"))
PrintDigits(option);
+
+ //increment the total number of Digits per run
+ fDigitsInRun += gime->Digits()->GetEntriesFast() ;
}
if(strstr(option,"tim")){
if(strstr(option,"deb"))
PrintRecParticles(option) ;
+
+ //increment the total number of rec particles per run
+ fRecParticlesInRun += gime->RecParticles()->GetEntriesFast() ;
+
}
if(strstr(option,"tim")){
cout << "AliPHOSPIDv1: event "<<gAlice->GetEvNumber() << endl ;
cout << " found " << recParticles->GetEntriesFast() << " RecParticles " << endl ;
-
- fRecParticlesInRun += recParticles->GetEntriesFast() ;
if(strstr(option,"all")) { // printing found TS
if(strstr(option,"deb"))
PrintTrackSegments(option) ;
+
+ //increment the total number of track segments per run
+ fTrackSegmentsInRun += gime->TrackSegments()->GetEntriesFast() ;
}
cout << "AliPHOSTrackSegmentMakerv1: event "<<gAlice->GetEvNumber() << endl ;
cout << " Found " << trackSegments->GetEntriesFast() << " trackSegments " << endl ;
- fTrackSegmentsInRun += trackSegments->GetEntriesFast() ;
-
if(strstr(option,"all")) { // printing found TS
cout << "TrackSegment # " << " EMC RP# " << " CPV RP# " << endl ;
gMC->CurrentVolOffID(10, moduleNumber) ; // get the PHOS module number ;
- // fill the relevant QA Checkables
- fQATotEner->Update( xyze[4] ) ; // total energy in PHOS
- (static_cast<AliPHOSQAFloatCheckable*>((*fQATotEnerB)[moduleNumber-1]))->Update( xyze[4] ) ; // energy in this block
-
Int_t strip ;
gMC->CurrentVolOffID(3, strip);
Int_t cell ;
// add current hit to the hit list
AddHit(fIshunt, primary,tracknumber, absid, xyze);
-
+ // fill the relevant QA Checkables
+ fQATotEner->Update( xyze[4] ) ; // total energy in PHOS
+ (static_cast<AliPHOSQAFloatCheckable*>((*fQATotEnerB)[moduleNumber-1]))->Update( xyze[4] ) ; // energy in this block
+
} // there is deposited energy
} // we are inside a PHOS Xtal
-############################### PHOS Makefile #################################
+ ############################### PHOS Makefile #################################
# Include machine specific definitions
@echo " "
@echo " Starting the test of the simulation/reconstruction software. Please don't take the warning messages into account. "
@echo " ____________________________________________________________ "
- @aliroot -b -q "testsim.C(100)" > out
- @aliroot -b -q testsimglobal.C > out
+ @aliroot -b -q "test.C(100)" > out
@rm out
@rm testPHOS.root
@echo " "
@echo " Starting the test of the simulation/reconstruction software. Please don't take the warning messages into account. "
@echo " ____________________________________________________________ "
- @aliroot -b -q "testsim.C(10)" > out
- @aliroot -b -q testsimglobal.C > out
+ @aliroot -b -q "test.C(10)" > out
@rm out
@rm testPHOS.root
--- /dev/null
+void test (Int_t nevent=1)
+{
+ gSystem->SetIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include -I.");
+ gROOT->ProcessLine(".x testsim.C(nevent)") ;
+ gROOT->ProcessLine(".x testsimglobal.C++") ;
+}
void testsimglobal()
{
+ gSystem->SetIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include -I.");
Bool_t error = kFALSE ;
TString mess("") ;
// Definition of the alarm bounds for the examination
const Float_t maxAlaHitsM = 12.79 ; // total multiplicity
- const Float_t maxAlaTotEn = 19.34 ; // total energy
- const Float_t maxAlaTotEnB = 18.35 ; // per block multiplicity
+ const Float_t maxAlaTotEn = 19.34 ; // total energy
+ const Float_t maxAlaTotEnB = 18.35 ; // per block multiplicity
const Float_t maxAlaHitsMB = 11.1 ; // per block energy
TString mess("") ;
//Digits process
AliPHOSGetter * gime = AliPHOSGetter::GetInstance("testPHOS.root") ;
TString reconame = "test suite" ;
- const Float_t maxDigits = 3489.41 ;
+ const Float_t maxDigits = 2860. ;
const Float_t widDigits = TMath::Sqrt(maxDigits) ;
TString mess("") ;
AliPHOSGetter * gime = AliPHOSGetter::GetInstance("testPHOS.root") ;
TString reconame = "test suite" ;
- const Float_t maxRecPoints = 222.83 ;
+ const Float_t maxRecPoints = 1.0 ;
const Float_t widRecPoints = TMath::Sqrt(maxRecPoints) ;
TString mess("") ;