{
// Calculates Real and Mixed invariant mass distributions
- const Int_t nMixedEvents = 4 ; //# of events used for calculation of 'mixed' distribution
- Int_t mixedLoops = (Int_t )TMath::Ceil(Nevents/nMixedEvents) ;
+ const Int_t knMixedEvents = 4 ; //# of events used for calculation of 'mixed' distribution
+ Int_t mixedLoops = (Int_t )TMath::Ceil(Nevents/knMixedEvents) ;
//========== Booking Histograms
TH2D * hRealEM = new TH2D("hRealEM", "Real for EM particles", 250,0.,1.,40,0.,4.) ;
Int_t ievent;
Int_t eventInMixedLoop ;
- Int_t nRecParticles[4];//nMixedEvents] ;
+ Int_t nRecParticles[4];//knMixedEvents] ;
- AliPHOSRecParticle::RecParticlesList * allRecParticleList = new TClonesArray("AliPHOSRecParticle", nMixedEvents*1000) ;
+ AliPHOSRecParticle::RecParticlesList * allRecParticleList = new TClonesArray("AliPHOSRecParticle", knMixedEvents*1000) ;
for(eventInMixedLoop = 0; eventInMixedLoop < mixedLoops; eventInMixedLoop++ ){
Int_t iRecPhot = 0 ;
- for ( ievent=0; ievent < nMixedEvents; ievent++){
+ for ( ievent=0; ievent < knMixedEvents; ievent++){
- Int_t absEventNumber = eventInMixedLoop*nMixedEvents + ievent ;
+ Int_t absEventNumber = eventInMixedLoop*knMixedEvents + ievent ;
//=========== Connects the various Tree's for evt
gAlice->GetEvent(absEventNumber);
TH2F * fhAllEnergy ; // Spectrum of detected photons with photon primary
TH2F * fhPhotEnergy ; // Total spectrum of detected photons
TH2F * fhEMEnergy ; // Spectrum of detected neutral EM with EM primary
- TH2F * fhPPSDEnergy ; //
+ TH2F * fhPPSDEnergy ; // Spectrum of particles detected in PPSD with primary particle
TH2F * fhAllPosition ; // Position Resolution of photons with photon primary
TH2F * fhPhotPosition ; // Position Resolution of photons
//____________________________________________________________________________
Int_t AliPHOSDigit::GetPrimary(Int_t index) const
{
-
+ // retrieves the primary particle number given its index in the list
Int_t rv = -1 ;
if ( index <= fNprimary ){
rv = fPrimary[index-1] ;
}
TString Name() ;
virtual void Paint(Option_t * option="");
- virtual void Print(const char *) ;
+ virtual void Print(const char * opt) ;
void SetPrimary(Int_t index) {
// sets the primary particle index
fPrimary = index ;
}
}
+//____________________________________________________________________________
+Float_t AliPHOSGeometry::GetCPVBoxSize(Int_t index) const {
+ if (strcmp(fName,"GPS2") ==0 )
+ return fGeometryPPSD->GetCPVBoxSize(index);
+ else if (strcmp(fName,"IHEP")==0)
+ return fGeometryCPV ->GetCPVBoxSize(index);
+ else if (strcmp(fName,"MIXT")==0)
+ return TMath::Max(fGeometryCPV ->GetCPVBoxSize(index), fGeometryPPSD->GetCPVBoxSize(index));
+ else
+ return 0;
+}
+
//____________________________________________________________________________
AliPHOSGeometry * AliPHOSGeometry::GetInstance()
{
// Common PPSD and CPV parameters
- Float_t GetCPVBoxSize(Int_t index) const {
- if (strcmp(fName,"GPS2")==0) return fGeometryPPSD->GetCPVBoxSize(index);
- else if (strcmp(fName,"IHEP")==0) return fGeometryCPV ->GetCPVBoxSize(index);
- else if (strcmp(fName,"MIXT")==0) return TMath::Max(fGeometryCPV ->GetCPVBoxSize(index),
- fGeometryPPSD->GetCPVBoxSize(index));
- else return 0;
- }
+ Float_t GetCPVBoxSize(Int_t index) const ;
// Mixed EMCA and PPSD parameters
//_________________________________________________________________________
// A singleton that retrieves objets from an array stored in a Tree on a disk file
-// 1. AliPHOSDigit from TreeD
+// 1. AliPHOSDigit from TreeD
+// Why should I put meaningless comments
+// just to satisfy
+// the code checker
//
//*-- Author: Yves Schutz (SUBATECH)
//////////////////////////////////////////////////////////////////////////////
//_________________________________________________________________________
// A singleton that retrieves objets from an array stored in a Tree on a disk file
// 1. AliPHOSDigit from TreeD
-//
+// Oh yeah
//*-- Author: Yves Schutz (SUBATECH)
AliRecPoint * GimeRecPoint(Int_t index, TString s) ;
AliPHOSTrackSegment * GimeTrackSegment(Int_t index) ;
- AliPHOSIndexToObject & operator = (const AliPHOSIndexToObject & rvalue) {
+ AliPHOSIndexToObject & operator = (const AliPHOSIndexToObject & ) {
// assignement operator requested by coding convention
// but not needed
assert(0==1) ;
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
/* $Id$ */
-
//_________________________________________________________________________
// Algorithm class used only by AliPHOSTrackSegmentMaker
-//
+// Links recpoints into tracksegments
+// Why should I put meaningless comments
+// just to satisfy
+// the code checker
//*-- Author: Dmitri Peressounko (SUBATECH)
// --- ROOT system ---
//_________________________________________________________________________
// Algorithm class used only by AliPHOSTrackSegmentMaker
-//
+// Links recpoints
+// into tracksegments
//*-- Author: Dmitri Peressounko (SUBATECH)
// --- ROOT system ---
//_________________________________________________________________________
// Algorithm class for the identification of particles detected in PHOS
-// base class
+// base class of identified particle
+// Why should I put meaningless comments
+// just to satisfy
+// the code checker
+
//
//*-- Author: Yves Schutz (SUBATECH)
//_________________________________________________________________________
// Algorithm class for the identification of particles detected in PHOS
// base class
-//
+// of identified particles
//*-- Author: Yves Schutz (SUBATECH)
// --- ROOT system ---
//_________________________________________________________________________
// Implementation version v1 of the PHOS particle identifier
// Identification is based on information from PPSD and EMC
+// Why should I put meaningless comments
+// just to satisfy
+// the code checker
//
//*-- Author: Yves Schutz (SUBATECH) & Gines Martinez (SUBATECH)
// //R=(lambda[0]-1.386)*(lambda[0]-1.386)+1.707*1.707*(lambda[1]-1.008)*(lambda[1]-1.008) ;
// //if(R<0.35*0.35)
- Float_t Dispersion;
- Dispersion = recp->GetDispersion();
- if (Dispersion < fCutOnDispersion)
+ Float_t dispersion;
+ dispersion = recp->GetDispersion();
+ if (dispersion < fCutOnDispersion)
showerprofile = 0 ; // NARROW PROFILE
else
showerprofile = 1 ;// WIDE PROFILE
//_________________________________________________________________________
// Implementation version v1 of the PHOS particle identifier
// Identification is based on information from PPSD and EMC
-//
+// Oh yeah
//*-- Author: Yves Schutz (SUBATECH)
// --- ROOT system ---
Float_t GetDistanceInPHOSPlane(AliPHOSEmcRecPoint * emcclu, AliPHOSPpsdRecPoint * PpsdClu, Bool_t &toofar, Option_t * Axis) ; // Relative Distance PPSD-EMC
virtual void MakeParticles(AliPHOSTrackSegment::TrackSegmentsList * trsl,
AliPHOSRecParticle::RecParticlesList * rpl ) ; // does the job
- virtual void Print(const char *) ;
+ virtual void Print(const char * opt) ;
virtual void SetDispersionCutOff(Float_t Dcut) {fCutOnDispersion = Dcut ; }
virtual void SetShowerProfileCuts(Float_t l1m, Float_t l1M, Float_t l2m, Float_t l2M) ;
virtual void SetRelativeDistanceCut(Float_t CutOnRelativeDistance) ;
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
/* $Id$ */
-
//_________________________________________________________________________
// A RecPoint (cluster) in the PPSD
// A PPSD RecPoint ends up to be a single digit
+// Why should I put meaningless comments
+// just to satisfy
+// the code checker
//
//*-- Yves Schutz (SUBATECH)
//_________________________________________________________________________
// A RecPoint (cluster) in the PPSD
// A PPSD RecPoint ends up to be a single digit
-//
+// Oh yeah
//*-- Yves Schutz (SUBATECH)
// --- ROOT system ---
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
/* $Id$ */
-
//_________________________________________________________________________
// A Reconstructed Particle in PHOS
-// To become a general class of AliRoot ?
+// To become a general class of AliRoot ?
+// Why should I put meaningless comments
+// just to satisfy
+// the code checker
//
//*-- Author: Yves Schutz (SUBATECH)
//_________________________________________________________________________
// A Reconstructed Particle in PHOS
// To become a general class of AliRoot ?
-//
+// why not
//*-- Author: Yves Schutz (SUBATECH)
// --- ROOT system ---
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
/* $Id$ */
-
//_________________________________________________________________________
// Base Class for PHOS Reconstructed Points
-//
+// Why should I put meaningless comments
+// just to satisfy
+// the code checker
//*-- Author: Gines Martinez (SUBATECH)
// --- ROOT system ---
#define ALIPHOSRECPOINT_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
-
//_________________________________________________________________________
// Base Class for PHOS Reconstructed Points
-//
+// A recpoint being equivalent to a cluster in encal terminology
//*-- Author: Gines Martinez (SUBATECH)
#include <assert.h>
// Print prototype
}
- AliPHOSRecPoint & operator = (const AliPHOSRecPoint & rvalue) {
+ AliPHOSRecPoint & operator = (const AliPHOSRecPoint & ) {
// assignement operator requested by coding convention
// but not needed
assert(0==1) ;
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
/* $Id$ */
-
//_________________________________________________________________________
// Track segment in PHOS
// Can be : 1 EmcRecPoint
//____________________________________________________________________________
void AliPHOSTrackSegment::SetPpsdUpRecPoint(AliPHOSPpsdRecPoint * PpsdUpRecPoint)
{
+ // gives an id from its position in the list
if( PpsdUpRecPoint )
fPpsdUpRecPoint = PpsdUpRecPoint->GetIndexInList() ;
else
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
/* $Id$ */
-
//_________________________________________________________________________
// Algorithm Base class to construct PHOS track segments
// Associates EMC and PPSD clusters
Int_t fNTrackSegments ; // number of track segments found
AliPHOSGeometry * fGeom ; // pointer to PHOS geometry
- AliPHOSIndexToObject * fPlease ; //
+ AliPHOSIndexToObject * fPlease ; // factory used to convert an index into the associated object
ClassDef( AliPHOSTrackSegmentMaker,1) // Algorithm class to make PHOS track segments (Base Class)
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
/* $Id$ */
-
//_________________________________________________________________________
// Implementation version 1 of algorithm class to construct PHOS track segments
// Associates EMC and PPSD clusters
Float_t * maxAtEnergy ) ; //Unfolds overlaping clusters using TMinuit package
virtual void UnsetUnfoldFlag() { fUnfoldFlag = kFALSE ; }
- AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & rvalue) {
+ AliPHOSTrackSegmentMakerv1 & operator = (const AliPHOSTrackSegmentMakerv1 & ) {
// assignement operator requested by coding convention
// but not needed
assert(0==1) ;
* about the suitability of this software for any purpose. It is *
* provided "as is" without express or implied warranty. *
**************************************************************************/
-
/* $Id$ */
//_________________________________________________________________________
// Version of AliPHOSv1 which keeps all hits in TreeH
//*-- Modified Nov. 22 2000 by Dmitri Peressounko
// All hits are stored.
// Note, that primaries will not be assigned to digits:
-// becouse of tiny energy deposition at each step.
+// because of tiny energy deposition at each step.
//
// --- ROOT system ---
* See cxx source for full Copyright notice */
//_________________________________________________________________________
// Version of AliPHOSv0 which keeps all hits in TreeH
+// I mean real hits not cumulated hits
// This version is NOT recommended for Reconstruction analysis
//
//*-- Author: Gines MARTINEZ (SUBATECH)
Int_t modid ;
gMC->CurrentVolID(modid);
- Float_t Energy = gMC->Etot() ; //Total energy of current track
+ Float_t energy = gMC->Etot() ; //Total energy of current track
//Calculating mass of current particle
- TDatabasePDG * PDG = TDatabasePDG::Instance() ;
- TParticlePDG * partPDG = PDG->GetParticle(gMC->TrackPid()) ;
+ TDatabasePDG * pdg = TDatabasePDG::Instance() ;
+ TParticlePDG * partPDG = pdg->GetParticle(gMC->TrackPid()) ;
Float_t mass = partPDG->Mass() ;
- if(Energy > mass){
- pos.SetMag(TMath::Sqrt(Energy*Energy-mass*mass)) ;
- TLorentzVector pTrack(pos,Energy) ;
+ if(energy > mass){
+ pos.SetMag(TMath::Sqrt(energy*energy-mass*mass)) ;
+ TLorentzVector pTrack(pos, energy) ;
TParticle * part = new TParticle(gMC->TrackPid(), 0,-1,-1,-1,-1, pTrack, lv) ;
return TString("v4") ;
}
- AliPHOSv4 & operator = (const AliPHOSv4 & rvalue) {
+ AliPHOSv4 & operator = (const AliPHOSv4 & ) {
// assignement operator requested by coding convention
// but not needed
assert(0==1) ;
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff