ecalAnalysisMaterial Class Reference

#include <ecalAnalysisMaterial.h>

Inheritance diagram for ecalAnalysisMaterial:

Collaboration diagram for ecalAnalysisMaterial:

Public Member Functions
	ecalAnalysisMaterial (const char *name, const Int_t iVerbose=0)
	ecalAnalysisMaterial ()
virtual InitStatus	Init ()
virtual void	Exec (Option_t *option)
virtual void	Finish ()
virtual	~ecalAnalysisMaterial ()

Private Member Functions
void	InitTree ()
	ecalAnalysisMaterial (const ecalAnalysisMaterial &)
ecalAnalysisMaterial &	operator= (const ecalAnalysisMaterial &)

Private Attributes
TTree *	fTree
Int_t	fEvent
Double_t	fX
Double_t	fY
Double_t	fZ
Double_t	fX0
Double_t	fX0z [400]
Double_t	fIntL
TClonesArray *	fMC
TClonesArray *	fRadLen

Detailed Description

Definition at line 14 of file ecalAnalysisMaterial.h.

Constructor & Destructor Documentation

ecalAnalysisMaterial() [1/3]

ecalAnalysisMaterial::ecalAnalysisMaterial	(	const char *	name,
		const Int_t	iVerbose = 0
	)

Definition at line 96 of file ecalAnalysisMaterial.cxx.

  : FairTask(name, iVerbose),
    fTree(NULL),
    fEvent(0),
    fX(0.),
    fY(0.),
    fZ(0.),
    fX0(0.),
    fIntL(0.),
    fMC(NULL),
    fRadLen(NULL)
{
}

ecalAnalysisMaterial() [2/3]

ecalAnalysisMaterial::ecalAnalysisMaterial

(

)

Default constructor

Definition at line 82 of file ecalAnalysisMaterial.cxx.

  : FairTask(),
    fTree(NULL),
    fEvent(0),
    fX(0.),
    fY(0.),
    fZ(0.),
    fX0(0.),
    fIntL(0.),
    fMC(NULL),
    fRadLen(NULL)
{
}

~ecalAnalysisMaterial()

virtual ecalAnalysisMaterial::~ecalAnalysisMaterial ( )

inlinevirtual

Destructor

Definition at line 32 of file ecalAnalysisMaterial.h.

{};

ecalAnalysisMaterial() [3/3]

ecalAnalysisMaterial::ecalAnalysisMaterial ( const ecalAnalysisMaterial &  )

private

Member Function Documentation

Exec()

void ecalAnalysisMaterial::Exec ( Option_t *  option )

virtual

Loop procedure

Definition at line 20 of file ecalAnalysisMaterial.cxx.

{
  fEvent++;
  if (fVerbose>0)
    Info("Exec", "Event %d.", fEvent);
  InitTree();
  ShipMCTrack* tr;
  FairRadLenPoint* p;
  Int_t i;
  Int_t k;
  Int_t n;
  Double_t z=3535;
  TVector3 newp;
  TVector3 oldp(0.0, 0.0, 0.0);
  TVector3 d;
  if (fMC->GetEntriesFast()!=1)
  {
    cout << "Size of MC points arrays differs from one" << endl;
  }
  tr=(ShipMCTrack*)fMC->At(0);
  fX=tr->GetPx();
  fY=tr->GetPy();
  fZ=tr->GetPz();
  fX*=z/fZ;
  fY*=z/fZ;
  fZ*=z/fZ;
//  if (TMath::Abs(fX)>594||TMath::Abs(fY)>474) return;
//  if (TMath::Abs(fX)<42&&TMath::Abs(fY)<42) return;
  n=fRadLen->GetEntriesFast();
  fX0=0; fIntL=0;
  for(i=0;i<400;i++) fX0z[i]=0;
  for(i=0;i<n;i++)
  {
    p=(FairRadLenPoint*)fRadLen->At(i);
    newp=p->GetPositionOut();
//    cout << d.Mag() << "  " << p->GetRadLength() << " " << newp.Z() << "  " << fX0 << endl;
    if (newp.Z()>3535) break;
    k=newp.Z()/10;
    d=newp-oldp;
//    cout << p->GetLength() << "   " << p->GetRadLength() << " " << fEvent << "    " <<p->GetXOut() << "   " << p->GetYOut() << "  " << p->GetZOut() << endl;
    fX0+=d.Mag()/p->GetRadLength();
    fX0z[k]+=d.Mag()/p->GetRadLength();
    oldp=newp;
  }
  fTree->Fill();
}

Finish()

void ecalAnalysisMaterial::Finish ( )

virtual

Finishing routine

Definition at line 137 of file ecalAnalysisMaterial.cxx.

{
  if (fTree)
    fTree->Write();
}

Init()

InitStatus ecalAnalysisMaterial::Init ( )

virtual

Initing routine

Definition at line 111 of file ecalAnalysisMaterial.cxx.

{
  FairRootManager* fManager=FairRootManager::Instance();
  if (!fManager)
  {
    Fatal("Init", "Can't find a Root Manager.");
    return kFATAL;
  }
  fMC=(TClonesArray*)fManager->GetObject("MCTrack");
  if (!fMC)
  {
    Fatal("Init", "Can't find an array of MC tracks.");
    return kFATAL;
  }
  fRadLen=(TClonesArray*)fManager->GetObject("RadLen");
  if (!fRadLen)
  {
    Fatal("Init", "Can't find an array of radiation length points.");
    return kFATAL;
  }
  fTree=NULL;
 
  return kSUCCESS;
}

InitTree()

void ecalAnalysisMaterial::InitTree ( )

private

Initialize a tree

Initializes tree

Definition at line 69 of file ecalAnalysisMaterial.cxx.

{
  if (fTree) return;
  fTree=new TTree("rad","Radiation lenght tree");
  fTree->Branch("ev", &fEvent, "ev/I");
  fTree->Branch("x", &fX, "x/D");
  fTree->Branch("y", &fY, "Y/D");
  fTree->Branch("z", &fZ, "z/D");
  fTree->Branch("x0", &fX0, "x0/D");
  fTree->Branch("x0z", &fX0z, "x0z[400]/D");
  fTree->Branch("intl", &fIntL, "intl/D");
}

operator=()

ecalAnalysisMaterial & ecalAnalysisMaterial::operator= ( const ecalAnalysisMaterial &  )

private

Member Data Documentation

fEvent

Int_t ecalAnalysisMaterial::fEvent

private

Event number

Definition at line 40 of file ecalAnalysisMaterial.h.

fIntL

Double_t ecalAnalysisMaterial::fIntL

private

Definition at line 47 of file ecalAnalysisMaterial.h.

fMC

TClonesArray* ecalAnalysisMaterial::fMC

private

MC tracks array

Definition at line 51 of file ecalAnalysisMaterial.h.

fRadLen

TClonesArray* ecalAnalysisMaterial::fRadLen

private

Array of radlen points

Definition at line 53 of file ecalAnalysisMaterial.h.

fTree

TTree* ecalAnalysisMaterial::fTree

private

Out tree

Definition at line 37 of file ecalAnalysisMaterial.h.

fX

Double_t ecalAnalysisMaterial::fX

private

Information particle

Definition at line 42 of file ecalAnalysisMaterial.h.

fX0

Double_t ecalAnalysisMaterial::fX0

private

Definition at line 45 of file ecalAnalysisMaterial.h.

fX0z

Double_t ecalAnalysisMaterial::fX0z[400]

private

Definition at line 46 of file ecalAnalysisMaterial.h.

fY

Double_t ecalAnalysisMaterial::fY

private

Definition at line 43 of file ecalAnalysisMaterial.h.

fZ

Double_t ecalAnalysisMaterial::fZ

private

Definition at line 44 of file ecalAnalysisMaterial.h.

---

The documentation for this class was generated from the following files:

• ecal/ecalAnalysisMaterial.h
• ecal/ecalAnalysisMaterial.cxx