for(ievent = 0; ievent < nevents; ievent++){
+ gAlice->GetEvent(ievent) ;
+ gAlice->SetEvent(ievent) ;
+
if(!ReadDigits(ievent)) //reads digits for event fEvent
- return;
+ continue;
MakeClusters() ;
fNumberOfEmcClusters = 0 ;
fNumberOfCpvClusters = 0 ;
- // Get Digits Tree header from file
- gAlice->GetEvent(event) ;
- gAlice->SetEvent(event) ;
-
if ( gAlice->TreeD()==0) {
cerr << "ERROR: AliPHOSClusterizerv1::ReadDigits There is no Digit Tree" << endl;
return kFALSE;
TClonesArray * digits = gime->Digits() ;
Int_t index ;
- //Evaluate poisition, dispersion and other RecPoint properties...
+ //Evaluate position, dispersion and other RecPoint properties...
for(index = 0; index < emcRecPoints->GetEntries(); index++)
((AliPHOSEmcRecPoint *)emcRecPoints->At(index))->EvalAll(fW0,digits) ;
TString branchName(GetName()) ;
branchName.ReplaceAll(Version(), "") ;
- //First EMC
Int_t bufferSize = 32000 ;
Int_t splitlevel = 0 ;
+
+ //First EMC
TBranch * emcBranch = gAlice->TreeR()->Branch("PHOSEmcRP","TObjArray",&emcRecPoints,bufferSize,splitlevel);
emcBranch->SetTitle(branchName);
if (filename) {
while ((sb=(TBranch*)next())) {
sb->SetFile(filename);
}
-
+
cwd->cd();
}
}
cwd->cd();
}
-
- emcBranch->Fill() ;
- cpvBranch->Fill() ;
+ emcBranch ->Fill() ;
+ cpvBranch ->Fill() ;
clusterizerBranch->Fill() ;
gAlice->TreeR()->Write(0,kOverwrite) ;
Float_t * emcEnergies = emcRP->GetEnergiesList() ;
- AliPHOSGeometry * geom = AliPHOSGeometry::GetInstance() ;
-
+ const AliPHOSGeometry * geom = AliPHOSGetter::GetInstance()->PHOSGeometry() ;
fret = 0. ;
Int_t iparam ;
TObjArray * emcRecPoints = AliPHOSGetter::GetInstance()->EmcRecPoints() ;
TObjArray * cpvRecPoints = AliPHOSGetter::GetInstance()->CpvRecPoints() ;
- cout << "AliPHOSClusterizerv1: " << endl ;
+ cout << "AliPHOSClusterizerv1: : event "<<gAlice->GetEvNumber() << endl ;
cout << " Found "<< emcRecPoints->GetEntriesFast() << " EMC Rec Points and "
<< cpvRecPoints->GetEntriesFast() << " CPV RecPoints" << endl ;
for(ievent = 0; ievent < nevents; ievent++){
+ gAlice->GetEvent(ievent);
+ gAlice->SetEvent(ievent) ;
if(!ReadSDigits(ievent)) //read sdigits event(s) evaluated by Combinator() from file(s)
- return ;
+ continue ;
Digitize(ievent) ; //Add prepared SDigits to digits and add the noise
WriteDigits(ievent) ;
- }
- if(strstr(option,"deb"))
- PrintDigits(option);
-
+ if(strstr(option,"deb"))
+ PrintDigits(option);
+ }
if(strstr(option,"tim")){
gBenchmark->Stop("PHOSDigitizer");
AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ;
TClonesArray * digits = gime->Digits() ;
- cout << "AliPHOSDigitiser:"<< endl ;
+ cout << "AliPHOSDigitiser: event " << gAlice->GetEvNumber() << endl ;
cout << " Number of entries in Digits list " << digits->GetEntriesFast() << endl ;
cout << endl ;
if(strstr(option,"all")){
TString treeName("TreeS") ;
treeName += event ;
TTree * treeS = (TTree*)file->Get(treeName.Data());
-
+
if(treeS==0){
cerr << "ERROR: AliPHOSDigitizer::ReadSDigits There is no SDigit Tree" << endl;
return kFALSE;
TClonesArray * sdigits = gime->SDigits() ;
sdigitsBranch->SetAddress(&sdigits) ;
+ sdigits->Clear();
AliPHOSSDigitizer * sdigitizer = gime->SDigitizer() ;
sdigitizerBranch->SetAddress(&sdigitizer) ;
fPedestal = sdigitizer->GetPedestalParameter() ;
fSlope = sdigitizer->GetCalibrationParameter() ;
-
- return kTRUE ;
+
+ return kTRUE ;
}
AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ;
TClonesArray * digits = gime->Digits() ;
-
- gAlice->GetEvent(event) ;
if(gAlice->TreeD()==0)
gAlice->MakeTree("D") ;
gAlice->TreeD()->Write(0,kOverwrite) ;
}
-
-
//Fist read Track Segment Branch and extract RecPointsBranch from fTSMaker
- gAlice->GetEvent(event) ;
- gAlice->SetEvent(event) ;
-
// Get TreeR header from file
if(gAlice->TreeR()==0){
cerr << "ERROR: AliPHOSPIDv1::ReadTrackSegments -> There is no Reconstruction Tree" << endl;
Int_t ievent ;
for(ievent = 0; ievent < nevents; ievent++){
+ gAlice->GetEvent(ievent) ;
+ gAlice->SetEvent(ievent) ;
+
if(!ReadTrackSegments(ievent)) //reads TrackSegments for event ievent
- return ;
+ continue ;
MakeRecParticles() ;
- cout << "MAKE" << endl ;
+
WriteRecParticles(ievent);
if(strstr(option,"deb"))
TObjArray * cpvRecPoints = gime->CpvRecPoints() ;
TClonesArray * trackSegments = gime->TrackSegments() ;
TClonesArray * recParticles = gime->RecParticles() ;
+ recParticles->Clear();
TIter next(trackSegments) ;
AliPHOSTrackSegment * ts ;
AliPHOSGetter *gime = AliPHOSGetter::GetInstance() ;
TClonesArray * recParticles = gime->RecParticles() ;
- gAlice->GetEvent(event) ;
-
//Make branch in TreeR for RecParticles
char * filename = 0;
if(gSystem->Getenv("CONFIG_SPLIT_FILE")!=0){ //generating file name
TClonesArray * recParticles = gime->RecParticles() ;
- cout << "AliPHOSPIDv1: " << endl ;
+ cout << "AliPHOSPIDv1: event "<<gAlice->GetEvNumber() << endl ;
cout << " found " << recParticles->GetEntriesFast() << " RecParticles " << endl ;
if(strstr(option,"all")) { // printing found TS
AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ;
TClonesArray * sdigits = gime->SDigits(GetName()) ;
- cout << "AliPHOSSDigitizer: " << endl ;
+ cout << "AliPHOSSDigitizer: event "<<gAlice->GetEvNumber() << endl ;
cout << " Number of entries in SDigits list " << sdigits->GetEntriesFast() << endl ;
cout << endl ;
(((AliPHOSRecPoint *) emcRecPoints->At(fEmcLast))->GetPHOSMod() == fModule );
fEmcLast ++) ;
-
//Now CPV clusters
Int_t totalCpv = cpvRecPoints->GetEntriesFast() ;
TObjArray * cpvRecPoints = gime->CpvRecPoints() ;
TClonesArray * trackSegments = gime->TrackSegments() ;
- //Make arrays to mark clusters already chousen
+ //Make arrays to mark clusters already chosen
Int_t * emcExist = 0;
if(fEmcLast > fEmcFirst)
emcExist = new Int_t[fEmcLast-fEmcFirst] ;
Int_t ievent ;
for(ievent = 0; ievent < nevents; ievent++){
+ gAlice->GetEvent(ievent) ;
+ gAlice->SetEvent(ievent) ;
+
if(!ReadRecPoints(ievent)) //reads RecPoints for event ievent
- return;
+ continue;
for(fModule = 1; fModule <= geom->GetNModules() ; fModule++ ){
fPpsdFirst= 0 ;
fPpsdLast = 0 ;
-
- gAlice->GetEvent(event) ;
- gAlice->SetEvent(event) ;
-
// Get TreeR header from file
if(gAlice->TreeR()==0){
cerr << "ERROR: AliPHOSTrackSegmentMakerv1::ReadRecPoints -> There is no Reconstruction Tree" << endl;
TObjArray * cpvRecPoints = gime->CpvRecPoints() ;
cpvRecPoints->Clear() ;
cpvbranch->SetAddress(&cpvRecPoints) ;
-
+
+ TClonesArray * trackSegments = gime->TrackSegments() ;
+ trackSegments->Clear() ;
AliPHOSClusterizer * clusterizer = 0 ;
clusterizerbranch->SetAddress(&clusterizer) ;
emcbranch->GetEntry(0) ;
cpvbranch->GetEntry(0) ;
+
clusterizerbranch->GetEntry(0) ;
+
+ printf("ReadRecPoint: EMC=%d, CPV=%d\n",emcRecPoints->GetEntries(),cpvRecPoints->GetEntries());
return kTRUE ;
AliPHOSGetter *gime = AliPHOSGetter::GetInstance() ;
TClonesArray * trackSegments = gime->TrackSegments() ;
- gAlice->GetEvent(event) ;
-
//Make branch in TreeR for TrackSegments
char * filename = 0;
if(gSystem->Getenv("CONFIG_SPLIT_FILE")!=0){ //generating file name
TClonesArray * trackSegments = AliPHOSGetter::GetInstance()->TrackSegments() ;
- cout << "AliPHOSTrackSegmentMakerv1: " << endl ;
+ cout << "AliPHOSTrackSegmentMakerv1: event "<<gAlice->GetEvNumber() << endl ;
cout << " Found " << trackSegments->GetEntriesFast() << " trackSegments " << endl ;
if(strstr(option,"all")) { // printing found TS