//____________________________________________________________________________
AliTOFPID::AliTOFPID(char* headerFile, char *cutsFile, const Option_t* opt):TTask("AliTOFPID","")
{
+ felectron = 0;
+ fpion = 0;
+ fkaon = 0;
+ fproton = 0;
fhfile = TFile::Open(headerFile); // connect file with ntuple
fcut = TFile::Open(cutsFile); // connect file for cuts
foutfileName=headerFile;
+ fNtuple = 0;
+ fgen = 0;
Init(opt);
// add Task to //root/Tasks folder
protected:
- ClassDef(AliTOFPID,1) // Task class for TOF pid
+ ClassDef(AliTOFPID,0) // Task class for TOF pid
};
}//event loop
+ // free used memory for ftail
+ if (ftail)
+ {
+ delete ftail;
+ ftail = 0;
+ }
// writing ntuple on output file
foutputfile->cd();
}
}
+ // speed-up the code
+ treehits->SetBranchStatus("*",0); // switch off all branches
+ treehits->SetBranchStatus("TOF*",1); // switch on only TOF
for (Int_t track=0; track<ntracks;track++) { // starting loop on primary tracks for the current event
itrack=0; // itrack: total number of selected TPC tracks
+ // speed-up the code
+ treehits->SetBranchStatus("*",0); // switch off all branches
+ treehits->SetBranchStatus("TPC*",1); // switch on only TPC
+
for (Int_t track=0; track<ntracks;track++) {
gAlice->ResetHits();
nbytes += treehits->GetEvent(track);
Bool_t operator == (const AliTOFReconstructioner & tofrec) const ;
private:
- TGeant3 *fg3; // pointer to GEANT geometry
- TFile *foutputfile; // pointer to output file
- TNtuple *foutputntuple; // pointer to output ntuple
+ TGeant3 *fg3; //! pointer to GEANT geometry
+ TFile *foutputfile; //! pointer to output file
+ TNtuple *foutputntuple; //! pointer to output ntuple
TF1 *fZnoise; // pointer to formula giving the noise along z direction
TF1 *ftail; // pointer to formula for time with tail
Int_t fdbg; // Flag for debug, 0 no debug, 1 debug
Int_t lasttrack=-1;
Int_t nset=0;
+
+ TH->SetBranchStatus("*",0); // switch off all branches
+ TH->SetBranchStatus("TOF*",1); // switch on only TOF
+
// Start loop on primary tracks in the hits containers
Int_t ntracks = static_cast<Int_t>(TH->GetEntries());