EmulsionDet Class Reference

#include <EmulsionDet.h>

Inheritance diagram for EmulsionDet:

Collaboration diagram for EmulsionDet:

Public Member Functions
	EmulsionDet (const char *name, Bool_t Active, const char *Title="EmulsionDet")
	EmulsionDet ()
virtual	~EmulsionDet ()
void	ConstructGeometry ()
virtual void	Initialize ()
virtual Bool_t	ProcessHits (FairVolume *v=0)
virtual void	Register ()
virtual TClonesArray *	GetCollection (Int_t iColl) const
virtual void	Reset ()
EmulsionDetPoint *	AddHit (Int_t trackID, Int_t detID, TVector3 pos, TVector3 mom, Double_t time, Double_t length, Double_t eLoss, Int_t pdgcode, TVector3 Lpos, TVector3 Lmom)
virtual void	CopyClones (TClonesArray *cl1, TClonesArray *cl2, Int_t offset)
virtual void	SetSpecialPhysicsCuts ()
virtual void	EndOfEvent ()
virtual void	FinishPrimary ()
virtual void	FinishRun ()
virtual void	BeginPrimary ()
virtual void	PostTrack ()
virtual void	PreTrack ()
virtual void	BeginEvent ()
void	DecodeBrickID (Int_t detID, Int_t &NWall, Int_t &NRow, Int_t &NColumn, Int_t &NPlate)
TString	PathBrickID (Int_t detID)
void	SetConfPar (TString name, Float_t value)
void	SetConfPar (TString name, Int_t value)
void	SetConfPar (TString name, TString value)
Float_t	GetConfParF (TString name)
Int_t	GetConfParI (TString name)
TString	GetConfParS (TString name)
	EmulsionDet (const EmulsionDet &)
EmulsionDet &	operator= (const EmulsionDet &)

Public Attributes
ClassDef(EmulsionDet, 5) private Int_t	fVolumeID
	track index
TLorentzVector	fPos
	volume id
TLorentzVector	fMom
	position at entrance
Double32_t	fTime
	momentum at entrance
Double32_t	fLength
	time
Double32_t	fELoss
	length
TClonesArray *	fEmulsionDetPointCollection
	energy loss
std::map< TString, Float_t >	conf_floats
std::map< TString, Int_t >	conf_ints
std::map< TString, TString >	conf_strings

Protected Member Functions
Int_t	InitMedium (const char *name)

Protected Attributes
Double_t	XDimension
Double_t	YDimension
Double_t	ZDimension
Int_t	fNWall
Int_t	fNRow
Int_t	fNCol
Int_t	fNTarget
Int_t	number_of_plates
Double_t	WallXDim
Double_t	WallYDim
Double_t	WallZDim
Double_t	TotalWallZDim
Double_t	WallZBorder_offset
Double_t	EmulsionThickness
Double_t	EmulsionX
Double_t	EmulsionY
Double_t	PlasticBaseThickness
Double_t	PassiveThickness
Double_t	EmPlateWidth
Double_t	AllPlateWidth
Double_t	BrickPackageX
Double_t	BrickPackageY
Double_t	BrickPackageZ
Double_t	BrickZ
Double_t	BrickY
Double_t	BrickX
Int_t	fPassiveOption
Double_t	TTrackerZ
Double_t	ShiftY
Double_t	ShiftX

Detailed Description

Definition at line 23 of file EmulsionDet.h.

Constructor & Destructor Documentation

EmulsionDet() [1/3]

EmulsionDet::EmulsionDet	(	const char *	name,
		Bool_t	Active,
		const char *	Title = "EmulsionDet"
	)

Definition at line 74 of file EmulsionDet.cxx.

: FairDetector(name, true, kEmulsionDet),
fTrackID(-1),
fVolumeID(-1),
fPos(),
fMom(),
fTime(-1.),
fLength(-1.),
fELoss(-1),
fEmulsionDetPointCollection(new TClonesArray("EmulsionDetPoint"))
{
}

EmulsionDet() [2/3]

EmulsionDet::EmulsionDet

(

)

Definition at line 61 of file EmulsionDet.cxx.

: FairDetector("EmulsionDet", "",kTRUE),
fTrackID(-1),
fVolumeID(-1),
fPos(),
fMom(),
fTime(-1.),
fLength(-1.),
fELoss(-1),
fEmulsionDetPointCollection(new TClonesArray("EmulsionDetPoint"))
{
}

~EmulsionDet()

EmulsionDet::~EmulsionDet ( )

virtual

Definition at line 87 of file EmulsionDet.cxx.

{
    if (fEmulsionDetPointCollection) {
        fEmulsionDetPointCollection->Delete();
        delete fEmulsionDetPointCollection;
    }
}

EmulsionDet() [3/3]

EmulsionDet::EmulsionDet

(

const EmulsionDet &

)

Member Function Documentation

AddHit()

EmulsionDetPoint * EmulsionDet::AddHit	(	Int_t	trackID,
		Int_t	detID,
		TVector3	pos,
		TVector3	mom,
		Double_t	time,
		Double_t	length,
		Double_t	eLoss,
		Int_t	pdgcode,
		TVector3	Lpos,
		TVector3	Lmom
	)

How to add your own point of type EmulsionDetPoint to the clones array

Definition at line 516 of file EmulsionDet.cxx.

{
    TClonesArray& clref = *fEmulsionDetPointCollection;
    Int_t size = clref.GetEntriesFast();
    return new(clref[size]) EmulsionDetPoint(trackID,detID, pos, mom,
                    time, length, eLoss, pdgCode,Lpos,Lmom);
}

BeginEvent()

virtual void EmulsionDet::BeginEvent ( )

inlinevirtual

Definition at line 61 of file EmulsionDet.h.

{;}

BeginPrimary()

virtual void EmulsionDet::BeginPrimary ( )

inlinevirtual

Definition at line 58 of file EmulsionDet.h.

{;}

ConstructGeometry()

void EmulsionDet::ConstructGeometry

(

)

Create the detector geometry

Definition at line 162 of file EmulsionDet.cxx.

{
    // configuration parameters
    XDimension = conf_floats["EmulsionDet/xdim"];
    YDimension = conf_floats["EmulsionDet/ydim"];
    ZDimension = conf_floats["EmulsionDet/zdim"];
    fNCol        = conf_ints["EmulsionDet/col"];
    fNRow     = conf_ints["EmulsionDet/row"];
    fNWall      = conf_ints["EmulsionDet/wall"];
    fNTarget = conf_ints["EmulsionDet/target"];
    WallXDim = conf_floats["EmulsionDet/WallXDim"];
    WallYDim = conf_floats["EmulsionDet/WallYDim"];
    WallZDim = conf_floats["EmulsionDet/WallZDim"];
    TotalWallZDim = conf_floats["EmulsionDet/TotalWallZDim"];
    WallZBorder_offset = conf_floats["EmulsionDet/WallZBorder_offset"];
    EmulsionThickness = conf_floats["EmulsionDet/EmTh"];
    EmulsionX = conf_floats["EmulsionDet/EmX"];
    EmulsionY = conf_floats["EmulsionDet/EmY"];
    PlasticBaseThickness = conf_floats["EmulsionDet/PBTh"];
    EmPlateWidth = conf_floats["EmulsionDet/EPlW"];
    PassiveThickness = conf_floats["EmulsionDet/PassiveTh"];
    AllPlateWidth = conf_floats["EmulsionDet/AllPW"];
    fPassiveOption = conf_ints["EmulsionDet/PassiveOption"];
    BrickPackageX = conf_floats["EmulsionDet/BrPackX"];
    BrickPackageY = conf_floats["EmulsionDet/BrPackY"];
    BrickPackageZ = conf_floats["EmulsionDet/BrPackZ"];
    BrickX = conf_floats["EmulsionDet/BrX"];
    BrickY = conf_floats["EmulsionDet/BrY"];
    BrickZ = conf_floats["EmulsionDet/BrZ"];
    number_of_plates = conf_ints["EmulsionDet/n_plates"];
    TTrackerZ = conf_floats["EmulsionDet/TTz"];
    ShiftX = conf_floats["EmulsionDet/ShiftX"];
    ShiftY = conf_floats["EmulsionDet/ShiftY"];
    
    int NTungstenPlatesTB24 = conf_ints["EmulsionDet/n_tungsten_plates_tb24"];
    bool testbeam_2024_setup = false;
    if (NTungstenPlatesTB24 > 0) {
      testbeam_2024_setup = true;
    }
 
    TGeoVolume *top=gGeoManager->FindVolumeFast("Detector");
    if(!top)  LOG(ERROR) << "no Detector volume found " ;
    gGeoManager->SetVisLevel(10);
 
        LOG(INFO) << "    X: "<< XDimension<< "  "<< YDimension<<"   Z: "<<ZDimension;
    LOG(INFO) <<"  BrickX: "<< BrickX<<"  Y: "<< BrickY<< "  Z: "<< BrickZ;
    //Materials 
    InitMedium("vacuum");
    TGeoMedium *vacuum =gGeoManager->GetMedium("vacuum");
    
    InitMedium("air");
    TGeoMedium *air =gGeoManager->GetMedium("air");
 
    InitMedium("Aluminum");
    TGeoMedium *Al  = gGeoManager->GetMedium("Aluminum");
 
    InitMedium("PlasticBase");
    TGeoMedium *PBase =gGeoManager->GetMedium("PlasticBase");
 
    InitMedium("NuclearEmulsion");
    TGeoMedium *NEmu =gGeoManager->GetMedium("NuclearEmulsion");
    
    TGeoMaterial *NEmuMat = NEmu->GetMaterial(); //I need the materials to build the mixture
    TGeoMaterial *PBaseMat = PBase->GetMaterial();
 
    TGeoMixture * emufilmmixture = new TGeoMixture("EmulsionFilmMixture", 2.00); // two nuclear emulsions separated by the plastic base
    Double_t frac_emu = NEmuMat->GetDensity() * 2 * EmulsionThickness /(NEmuMat->GetDensity() * 2 * EmulsionThickness + PBaseMat->GetDensity() * PlasticBaseThickness);
 
    emufilmmixture->AddElement(NEmuMat,frac_emu);
    emufilmmixture->AddElement(PBaseMat,1. - frac_emu);
 
    TGeoMedium *Emufilm = new TGeoMedium("EmulsionFilm",100,emufilmmixture);
 
    InitMedium("tungstensifon");
    TGeoMedium *tungsten = gGeoManager->GetMedium("tungstensifon");
 
 
    Int_t NPlates = number_of_plates; //Number of doublets emulsion + Pb    
 
 
    //TGeoVolume *top = gGeoManager->MakeBox("Top",vacuum,10.,10.,10.);
    //gGeoManager->SetTopVolume(top);
 
    //Definition of the target box containing emulsion bricks + target trackers (TT) 
         TGeoVolumeAssembly *volTarget  = new TGeoVolumeAssembly("volTarget");
 
    //
    //  //Volumes definition
    //    //
 
        /*For testbeam 2024 there are plastic plates filling the upstream part of the target.
      The last 28 layers have W+plastic plates.
      The plastic layers in Wall1, Brick 1 are the same thickness as the plastic base in TI18 conf.
      The plastic layers in Wall2, Brick 1 are 2 types of different thickness: 1 or 4/3 times the thickness of the TI18 plastic base.
      These layers are alternating starting with the 4/3x-thick one, then 1x-thick one.
      Whenever testbeam 2024 items are added in this function, there will be a note.
    */
    
    TGeoBBox *Walltot = new TGeoBBox("walltot",XDimension/2, YDimension/2, TotalWallZDim/2);
        TGeoBBox *Wallint = new TGeoBBox("wallint",WallXDim/2, WallYDim/2, WallZDim/2);
 
        TGeoTranslation * Wallborderpos = new TGeoTranslation("Walborderpos",0,0,WallZBorder_offset);
        Wallborderpos->RegisterYourself();
 
        TGeoCompositeShape *Wallborder = new TGeoCompositeShape("wallborder","walltot - (wallint:Walborderpos)");
        TGeoVolume *volWallborder = new TGeoVolume("volWallborder", Wallborder, Al);
        volWallborder->SetLineColor(kGray);
 
        TGeoVolume *volWall = new TGeoVolume("Wall",Wallint,air);
        TGeoVolume *volWall_2 = new TGeoVolume("Wall_2",Wallint,air);
    
    //Rows
    TGeoBBox *Row = new TGeoBBox("row",WallXDim/2, BrickY/2, BrickZ/2);
    TGeoVolume *volRow = new TGeoVolume("Row",Row,air);
    TGeoBBox *Row_W2 = new TGeoBBox("row_W2",WallXDim/2, BrickY/2, BrickZ/2);// testbeam 2024
    TGeoVolume *volRow_W2 = new TGeoVolume("Row_W2",Row_W2,air);
 
    //Bricks
    TGeoBBox *Brick = new TGeoBBox("brick", BrickX/2, BrickY/2, BrickZ/2);
    TGeoVolume *volBrick = new TGeoVolume("Brick",Brick,air);
    volBrick->SetLineColor(kCyan);
    volBrick->SetTransparency(1);
    
    TGeoBBox *Brick_W2 = new TGeoBBox("brick_W2", BrickX/2, BrickY/2, BrickZ/2);// testbeam 2024
    TGeoVolume *volBrick_W2 = new TGeoVolume("Brick_W2",Brick_W2,air);
    volBrick_W2->SetLineColor(kOrange);
    volBrick_W2->SetTransparency(1);
 
    TGeoBBox *Passive = new TGeoBBox("Passive", EmulsionX/2, EmulsionY/2, PassiveThickness/2);
    TGeoVolume *volPassive = new TGeoVolume("volPassive",Passive,tungsten);
    volPassive->SetTransparency(1);
    volPassive->SetLineColor(kGray);
    
    // For the testbeam 2024, the upstream part of the brick is filled with plastic
    float UpstreamPlasticThickness = (NPlates-NTungstenPlatesTB24)*AllPlateWidth;
    // it is a little thinner for Wall2, where the pлastic plates have 2 thicknesses
    float correction_W2 = NTungstenPlatesTB24/2*PlasticBaseThickness/3.;
    float UpstreamPlasticThickness_W2 = UpstreamPlasticThickness - correction_W2;
    TGeoBBox *Passive_plastic = new TGeoBBox("Passive_plastic", (EmulsionX-2.)/2, (EmulsionY-2.)/2, UpstreamPlasticThickness/2);
    TGeoVolume *volPassive_plastic = new TGeoVolume("volPassive_plastic",Passive_plastic,PBase);
    volPassive_plastic->SetLineColor(kMagenta);
    volPassive_plastic->SetVisibility(kTRUE);
    
    TGeoBBox *Passive_plastic_W2 = new TGeoBBox("Passive_plastic_W2", (EmulsionX-2.)/2, (EmulsionY-2.)/2, UpstreamPlasticThickness_W2/2);
    TGeoVolume *volPassive_plastic_W2 = new TGeoVolume("volPassive_plastic_W2",Passive_plastic_W2,PBase);
    volPassive_plastic_W2->SetLineColor(kMagenta+2);
    volPassive_plastic_W2->SetVisibility(kTRUE);
    
    float accumulatve_width = UpstreamPlasticThickness_W2; // testbeam 2024
 
    for(Int_t n=0; n<NPlates; n++)
    {
        if (n==0 && testbeam_2024_setup){// testbeam 2024
          volBrick->AddNode(volPassive_plastic, n, 
                       new TGeoTranslation(0,0,-BrickZ/2 + UpstreamPlasticThickness/2));// upstream plastic
          volBrick_W2->AddNode(volPassive_plastic_W2, n, new TGeoTranslation(0,0,-BrickZ/2 + UpstreamPlasticThickness_W2/2));// upstream plastic
        }
        if (n>=(NPlates-NTungstenPlatesTB24) || !testbeam_2024_setup) { // instrumented part for both TI18 and testbeam 2024
          volBrick->AddNode(volPassive, n, 
                       new TGeoTranslation(0,0,-BrickZ/2 + EmPlateWidth + PassiveThickness/2 + n*AllPlateWidth));// default(1x) plastic base thickness
          if (testbeam_2024_setup){ // testbeam 2024
            // interchange plastic plates of 1x or 4/3x default thickness
            volBrick_W2->AddNode(volPassive, n, 
                            new TGeoTranslation(0,0,-BrickZ/2 + (1+((n+1)%2)/3.)*EmPlateWidth + PassiveThickness/2 + accumulatve_width));
            accumulatve_width+=(1+((n+1)%2)/3.)*EmPlateWidth + PassiveThickness;
          }
        }
    }
 
    TGeoBBox *EmulsionFilm = new TGeoBBox("EmulsionFilm", EmulsionX/2, EmulsionY/2, EmPlateWidth/2);
    TGeoVolume *volEmulsionFilm = new TGeoVolume("Emulsion",EmulsionFilm,Emufilm); //TOP
    volEmulsionFilm->SetLineColor(kBlue);
    LOG(INFO) << "EmulsionDet : Passive option (0: all active, 1: all passive) set to " << fPassiveOption ;
    if (fPassiveOption == 0) AddSensitiveVolume(volEmulsionFilm);
    if (!testbeam_2024_setup){ // TI18
      for(Int_t n=0; n<NPlates+1; n++)
      {
        volBrick->AddNode(volEmulsionFilm, n, new TGeoTranslation(0,0,-BrickZ/2+ EmPlateWidth/2 + n*AllPlateWidth));
      }
    }
 
    volBrick->SetVisibility(kTRUE);
    if (testbeam_2024_setup) { // testbeam 2024
      volBrick_W2->SetVisibility(kTRUE);
    }
 
//alignment
    double dx_survey[5] = {conf_floats["EmulsionDet/Xpos0"],conf_floats["EmulsionDet/Xpos1"],conf_floats["EmulsionDet/Xpos2"],conf_floats["EmulsionDet/Xpos3"],conf_floats["EmulsionDet/Xpos4"]};
    double dy_survey[5] = {conf_floats["EmulsionDet/Ypos0"],conf_floats["EmulsionDet/Ypos1"],conf_floats["EmulsionDet/Ypos2"],conf_floats["EmulsionDet/Ypos3"],conf_floats["EmulsionDet/Ypos4"]};
    double dz_survey[5] = {conf_floats["EmulsionDet/Zpos0"],conf_floats["EmulsionDet/Zpos1"],conf_floats["EmulsionDet/Zpos2"],conf_floats["EmulsionDet/Zpos3"],conf_floats["EmulsionDet/Zpos4"]};
 
    top->AddNode(volTarget,1,new TGeoTranslation(0,0,0));
 
    //adding walls
    
    Double_t d_cl_z = - ZDimension/2;
 
    if (!testbeam_2024_setup){
      for(int l = 0; l < fNWall; l++) {
        volTarget->AddNode(volWallborder,l,new TGeoTranslation(-dx_survey[l]-XDimension/2., dz_survey[l]+YDimension/2., dy_survey[l]+TotalWallZDim/2.)); //the survey points refer to the down-left corner
        volTarget->AddNode(volWall,l,new TGeoTranslation(-dx_survey[l]-XDimension/2., dz_survey[l]+YDimension/2., dy_survey[l]+TotalWallZDim/2.+WallZBorder_offset)); //the survey points refer to the down-left corner
        d_cl_z += BrickZ + TTrackerZ;
      }
    }
    else { // testbeam 2024 
      volTarget->AddNode(volWallborder,0,new TGeoTranslation(-dx_survey[0]-XDimension/2., dz_survey[0]+YDimension/2., dy_survey[0]+TotalWallZDim/2.));
      volTarget->AddNode(volWall,0,new TGeoTranslation(-dx_survey[0]-XDimension/2., dz_survey[0]+YDimension/2., dy_survey[0]+TotalWallZDim/2.+WallZBorder_offset));
      d_cl_z += BrickZ + TTrackerZ;
      volTarget->AddNode(volWallborder,1,new TGeoTranslation(-dx_survey[1]-XDimension/2., dz_survey[1]+YDimension/2., dy_survey[1]+TotalWallZDim/2.));
      volTarget->AddNode(volWall_2,1,new TGeoTranslation(-dx_survey[1]-XDimension/2., dz_survey[1]+YDimension/2., dy_survey[1]+TotalWallZDim/2.+WallZBorder_offset));
      d_cl_z += BrickZ + TTrackerZ;
    }
 
    //adding rows
    Double_t d_cl_y = -WallYDim/2;
    
    for(int k= 0; k< fNRow; k++)
      {
      volWall->AddNode(volRow,k,new TGeoTranslation(0, d_cl_y + BrickY/2,0));
      if (testbeam_2024_setup){ // testbeam 2024 
        volWall_2->AddNode(volRow_W2,k,new TGeoTranslation(0, d_cl_y + BrickY/2,0));
      }
 
      // 2mm is the distance for the structure that holds the brick
      d_cl_y += BrickY;
      }
 
    //adding columns of bricks
    Double_t d_cl_x = -WallXDim/2;
    for(int j= 0; j < fNCol; j++)
      {
        volRow->AddNode(volBrick,j,new TGeoTranslation(d_cl_x+BrickX/2, 0, 0));
        if (testbeam_2024_setup){ // testbeam 2024 
          volRow_W2->AddNode(volBrick_W2,j,new TGeoTranslation(d_cl_x+BrickX/2, 0, 0));
        }
        d_cl_x += BrickX;
      }
}

CopyClones()

virtual void EmulsionDet::CopyClones ( TClonesArray *  cl1, TClonesArray *  cl2, Int_t  offset  )

inlinevirtual

Definition at line 53 of file EmulsionDet.h.

{;}

DecodeBrickID()

void EmulsionDet::DecodeBrickID	(	Int_t	detID,
		Int_t &	NWall,
		Int_t &	NRow,
		Int_t &	NColumn,
		Int_t &	NPlate
	)

Obtain info about brick position from detectorID

Definition at line 121 of file EmulsionDet.cxx.

                                                                                                    {
 
  NWall = detID/1E4;
  Int_t NTransverse = (detID - NWall*1E4)/1E3;
  switch (NTransverse){
    case (1):  
      NColumn = 1;
      NRow = 2;
      break;
      
    case (2):
      NColumn = 1;
      NRow = 1;
      break;
    
    case (3):
      NColumn = 2;
      NRow = 2;
      break;
 
    case (4):
      NColumn = 2;
      NRow = 1;
      break;
    }
 
  NPlate = detID - NWall*1E4 - NTransverse*1E3;
  
 
}

EndOfEvent()

void EmulsionDet::EndOfEvent ( )

virtual

Definition at line 485 of file EmulsionDet.cxx.

{
    fEmulsionDetPointCollection->Clear();
}

FinishPrimary()

virtual void EmulsionDet::FinishPrimary ( )

inlinevirtual

Definition at line 56 of file EmulsionDet.h.

{;}

FinishRun()

virtual void EmulsionDet::FinishRun ( )

inlinevirtual

Definition at line 57 of file EmulsionDet.h.

{;}

GetCollection()

TClonesArray * EmulsionDet::GetCollection ( Int_t  iColl ) const

virtual

Gets the produced collections

Definition at line 504 of file EmulsionDet.cxx.

{
    if (iColl == 0) { return fEmulsionDetPointCollection; }
    else { return NULL; }
}

GetConfParF()

Float_t EmulsionDet::GetConfParF ( TString  name )

inline

Definition at line 70 of file EmulsionDet.h.

{return conf_floats[name];} 

GetConfParI()

Int_t EmulsionDet::GetConfParI ( TString  name )

inline

Definition at line 71 of file EmulsionDet.h.

{return conf_ints[name];}

GetConfParS()

TString EmulsionDet::GetConfParS ( TString  name )

inline

Definition at line 72 of file EmulsionDet.h.

{return conf_strings[name];}

Initialize()

void EmulsionDet::Initialize ( )

virtual

Initialization of the detector is done here

Definition at line 95 of file EmulsionDet.cxx.

{
    FairDetector::Initialize();
}

InitMedium()

Int_t EmulsionDet::InitMedium ( const char *  name )

protected

Definition at line 101 of file EmulsionDet.cxx.

{
    static FairGeoLoader *geoLoad=FairGeoLoader::Instance();
    static FairGeoInterface *geoFace=geoLoad->getGeoInterface();
    static FairGeoMedia *media=geoFace->getMedia();
    static FairGeoBuilder *geoBuild=geoLoad->getGeoBuilder();
    
    FairGeoMedium *ShipMedium=media->getMedium(name);
    
    if (!ShipMedium)
    {
        Fatal("InitMedium","Material %s not defined in media file.", name);
        return -1111;
    }
    TGeoMedium* medium=gGeoManager->GetMedium(name);
    if (medium!=NULL)
        return ShipMedium->getMediumIndex();
    return geoBuild->createMedium(ShipMedium);
}

operator=()

EmulsionDet & EmulsionDet::operator=

(

const EmulsionDet &

)

PathBrickID()

TString EmulsionDet::PathBrickID

(

Int_t

detID

)

Definition at line 152 of file EmulsionDet.cxx.

                                           {
  //provide path to the node, from the detectorID
  Int_t NWall, NRow, NColumn, NPlate;
  DecodeBrickID(detID, NWall, NRow, NColumn, NPlate);
  
  TString emulsionpath = TString::Format("/cave_1/Detector_0/volTarget_1/Wall_%i/Row_%i/Brick_%i/Emulsion_%i",NWall-1,NRow-1,NColumn-1,NPlate -1);
 
  return emulsionpath;
}

PostTrack()

virtual void EmulsionDet::PostTrack ( )

inlinevirtual

Definition at line 59 of file EmulsionDet.h.

{;}

PreTrack()

virtual void EmulsionDet::PreTrack ( )

inlinevirtual

Definition at line 60 of file EmulsionDet.h.

{;}

ProcessHits()

Bool_t EmulsionDet::ProcessHits ( FairVolume *  v = 0 )

virtual

Method called for each step during simulation (see FairMCApplication::Stepping())

This method is called from the MC stepping

Definition at line 405 of file EmulsionDet.cxx.

{
    //Set parameters at entrance of volume. Reset ELoss.
    if ( gMC->IsTrackEntering() ) {
        fELoss  = 0.;
        fTime   = gMC->TrackTime() * 1.0e09;
        fLength = gMC->TrackLength();
        gMC->TrackPosition(fPos);
        gMC->TrackMomentum(fMom);
    }
    // Sum energy loss for all steps in the active volume
    fELoss += gMC->Edep();
    
    // Create EmulsionDetPoint at exit of active volume
    if ( gMC->IsTrackExiting()    ||
        gMC->IsTrackStop()       ||
        gMC->IsTrackDisappeared()   ) {
        fTrackID  = gMC->GetStack()->GetCurrentTrackNumber();
    gMC->CurrentVolID(fVolumeID);
 
        //finding number of wall, row and column
        Int_t detID = fVolumeID;
        Int_t MaxLevel = gGeoManager->GetLevel();
        const Int_t MaxL = MaxLevel;
        //cout << "MaxLevel = " << MaxL << endl;
        //cout << gMC->CurrentVolPath()<< endl;
    
 
        Int_t motherV[MaxL];
        Int_t NPlate =0;
        const char *name;
    
        name = gMC->CurrentVolName();
        //cout << name << endl;
        NPlate = detID;
    
        Int_t  NWall = -2, NColumn =-2, NRow =-2;
 
        for(Int_t i = 0; i <= MaxL;i++)
      {  
        gMC->CurrentVolOffID(i, motherV[i]);         
        const char *mumname = gMC->CurrentVolOffName(i);                 
        
        if(strcmp(mumname, "Brick") == 0) NColumn = (1-motherV[i]); //0 or 1 (0 higher x, 1 lower x, x are negative, so higher x are closer to the beam)
        if(strcmp(mumname, "Row") == 0) NRow = motherV[i]; // 0 or 1 (0 lower y, 1 higher y)
        if(strcmp(mumname, "Wall") == 0) NWall = motherV[i]; //0,1,2,3 (increasing along the beam verse)
        
      }
 
        detID = (NWall+1)*1E4+(NRow*2+NColumn+1)*1E3+(NPlate+1);
        fVolumeID = detID;
    //found number of row, column and wall
        //if (NColumn > 2 || NRow > 2 || NWall > 5) cout<<"Debug test for detID "<<detID<<" is "<<NColumn<<" "<<NRow<<" "<<NWall<<" "<<NPlate<<endl;
        TParticle* p=gMC->GetStack()->GetCurrentTrack();
    Int_t pdgCode = p->GetPdgCode();
    if (pdgCode == 22) { return kFALSE; } //discard only photons
    
        TLorentzVector Pos; 
        gMC->TrackPosition(Pos); 
        Double_t xmean = (fPos.X()+Pos.X())/2. ;      
        Double_t ymean = (fPos.Y()+Pos.Y())/2. ;      
        Double_t zmean = (fPos.Z()+Pos.Z())/2. ;     
        
        TLorentzVector Mom;
        gMC->TrackMomentum(Mom);
        AddHit(fTrackID, fVolumeID, TVector3(xmean, ymean,  zmean),
           TVector3(fMom.Px(), fMom.Py(), fMom.Pz()), fTime, fLength,
           fELoss,pdgCode,TVector3(Pos.X(), Pos.Y(), Pos.Z()),TVector3(Mom.Px(), Mom.Py(), Mom.Pz()) );
    
        // Increment number of muon det points in TParticle
        ShipStack* stack = (ShipStack*) gMC->GetStack();
        stack->AddPoint(kEmulsionDet);
    }
    
    return kTRUE;
}

Register()

void EmulsionDet::Register ( )

virtual

Registers the produced collections in FAIRRootManager.

This will create a branch in the output tree called TargetPoint, setting the last parameter to kFALSE means: this collection will not be written to the file, it will exist only during the simulation.

Definition at line 491 of file EmulsionDet.cxx.

{
    
    FairRootManager::Instance()->Register("EmulsionDetPoint", "EmulsionDet",
                                          fEmulsionDetPointCollection, kTRUE);
}

Reset()

void EmulsionDet::Reset ( )

virtual

has to be called after each event to reset the containers

Definition at line 510 of file EmulsionDet.cxx.

{
    fEmulsionDetPointCollection->Clear();
}

SetConfPar() [1/3]

void EmulsionDet::SetConfPar ( TString  name, Float_t  value  )

inline

Definition at line 67 of file EmulsionDet.h.

{conf_floats[name]=value;}

SetConfPar() [2/3]

void EmulsionDet::SetConfPar ( TString  name, Int_t  value  )

inline

Definition at line 68 of file EmulsionDet.h.

{conf_ints[name]=value;}

SetConfPar() [3/3]

void EmulsionDet::SetConfPar ( TString  name, TString  value  )

inline

Definition at line 69 of file EmulsionDet.h.

{conf_strings[name]=value;}

SetSpecialPhysicsCuts()

virtual void EmulsionDet::SetSpecialPhysicsCuts ( )

inlinevirtual

Definition at line 54 of file EmulsionDet.h.

{;}

Member Data Documentation

AllPlateWidth

Double_t EmulsionDet::AllPlateWidth

protected

Definition at line 121 of file EmulsionDet.h.

BrickPackageX

Double_t EmulsionDet::BrickPackageX

protected

Definition at line 124 of file EmulsionDet.h.

BrickPackageY

Double_t EmulsionDet::BrickPackageY

protected

Definition at line 125 of file EmulsionDet.h.

BrickPackageZ

Double_t EmulsionDet::BrickPackageZ

protected

Definition at line 126 of file EmulsionDet.h.

BrickX

Double_t EmulsionDet::BrickX

protected

Definition at line 130 of file EmulsionDet.h.

BrickY

Double_t EmulsionDet::BrickY

protected

Definition at line 129 of file EmulsionDet.h.

BrickZ

Double_t EmulsionDet::BrickZ

protected

Definition at line 128 of file EmulsionDet.h.

conf_floats

std::map<TString,Float_t> EmulsionDet::conf_floats

configuration parameters

Definition at line 93 of file EmulsionDet.h.

conf_ints

std::map<TString,Int_t> EmulsionDet::conf_ints

Definition at line 94 of file EmulsionDet.h.

conf_strings

std::map<TString,TString> EmulsionDet::conf_strings

Definition at line 95 of file EmulsionDet.h.

EmPlateWidth

Double_t EmulsionDet::EmPlateWidth

protected

Definition at line 120 of file EmulsionDet.h.

EmulsionThickness

Double_t EmulsionDet::EmulsionThickness

protected

Definition at line 114 of file EmulsionDet.h.

EmulsionX

Double_t EmulsionDet::EmulsionX

protected

Definition at line 115 of file EmulsionDet.h.

EmulsionY

Double_t EmulsionDet::EmulsionY

protected

Definition at line 116 of file EmulsionDet.h.

fELoss

Double32_t EmulsionDet::fELoss

length

Definition at line 88 of file EmulsionDet.h.

fEmulsionDetPointCollection

TClonesArray* EmulsionDet::fEmulsionDetPointCollection

energy loss

container for data points

Definition at line 91 of file EmulsionDet.h.

fLength

Double32_t EmulsionDet::fLength

time

Definition at line 87 of file EmulsionDet.h.

fMom

TLorentzVector EmulsionDet::fMom

position at entrance

Definition at line 85 of file EmulsionDet.h.

fNCol

Int_t EmulsionDet::fNCol

protected

Definition at line 103 of file EmulsionDet.h.

fNRow

Int_t EmulsionDet::fNRow

protected

Definition at line 103 of file EmulsionDet.h.

fNTarget

Int_t EmulsionDet::fNTarget

protected

Definition at line 104 of file EmulsionDet.h.

fNWall

Int_t EmulsionDet::fNWall

protected

Definition at line 103 of file EmulsionDet.h.

fPassiveOption

Int_t EmulsionDet::fPassiveOption

protected

Definition at line 132 of file EmulsionDet.h.

fPos

TLorentzVector EmulsionDet::fPos

volume id

Definition at line 84 of file EmulsionDet.h.

fTime

Double32_t EmulsionDet::fTime

momentum at entrance

Definition at line 86 of file EmulsionDet.h.

fVolumeID

ClassDef (EmulsionDet,5) private Int_t EmulsionDet::fVolumeID

track index

Track information to be stored until the track leaves the active volume.

Definition at line 83 of file EmulsionDet.h.

number_of_plates

Int_t EmulsionDet::number_of_plates

protected

Definition at line 105 of file EmulsionDet.h.

PassiveThickness

Double_t EmulsionDet::PassiveThickness

protected

Definition at line 119 of file EmulsionDet.h.

PlasticBaseThickness

Double_t EmulsionDet::PlasticBaseThickness

protected

Definition at line 118 of file EmulsionDet.h.

ShiftX

Double_t EmulsionDet::ShiftX

protected

Definition at line 141 of file EmulsionDet.h.

ShiftY

Double_t EmulsionDet::ShiftY

protected

Definition at line 139 of file EmulsionDet.h.

TotalWallZDim

Double_t EmulsionDet::TotalWallZDim

protected

Definition at line 111 of file EmulsionDet.h.

TTrackerZ

Double_t EmulsionDet::TTrackerZ

protected

Definition at line 136 of file EmulsionDet.h.

WallXDim

Double_t EmulsionDet::WallXDim

protected

Definition at line 107 of file EmulsionDet.h.

WallYDim

Double_t EmulsionDet::WallYDim

protected

Definition at line 108 of file EmulsionDet.h.

WallZBorder_offset

Double_t EmulsionDet::WallZBorder_offset

protected

Definition at line 112 of file EmulsionDet.h.

WallZDim

Double_t EmulsionDet::WallZDim

protected

Definition at line 109 of file EmulsionDet.h.

XDimension

Double_t EmulsionDet::XDimension

protected

Definition at line 99 of file EmulsionDet.h.

YDimension

Double_t EmulsionDet::YDimension

protected

Definition at line 100 of file EmulsionDet.h.

ZDimension

Double_t EmulsionDet::ZDimension

protected

Definition at line 101 of file EmulsionDet.h.

---

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

• shipLHC/EmulsionDet.h
• shipLHC/EmulsionDet.cxx