if(fNhits==0){
// First hit -> setting flag for primary or secondary particle
- /*TParticle * p = gAlice->GetMCApp()->Particle(track);
+ TParticle * p = gAlice->GetMCApp()->Particle(track);
Int_t imo = p->GetFirstMother();
//
if(track != imo){
}
else if(track == imo){
newquad->SetSFlag(0); // PRIMARY particle entering the ZDC
- }*/
+ }
//
sFlag = newquad->GetSFlag();
primKinEn = newquad->GetPrimKinEn();
motPDGcode = newquad->GetMotherPDGCode();
trackTime = newquad->GetTrackTOF();
trackEta = newquad->GetTrackEta();
+
}
else{
newquad->SetPrimKinEn(primKinEn);
newquad->SetTrackTOF(trackTime);
newquad->SetTrackEta(trackEta);
}
-
+
Int_t j;
for(j=0; j<fNhits; j++){
// If hits are equal (same track, same volume), sum them.
return;
}
}
+ //printf( "PDG from hits[10] = %f\n", hits[10]);
//Otherwise create a new hit
new(lhits[fNhits]) AliZDCHit(*newquad);
void AliZDCHit::Print(Option_t *) const
{
// Print method
- printf("\t **** AliZDCHit: track %d eta %f PDGcode %d TOF %1.1f ns E_prim = %1.2f GeV \n"
+ printf(" **** AliZDCHit: track %d eta %f PDGcode %d TOF %1.1f ns E_prim = %1.2f GeV \n"
"\t DETECTOR (%d, %d) (X, Y)_impact (%1.2f, %1.2f) cm Secflag %1.0f\n"
"\t PMQLight %1.0f, PMCLight %1.0f, E_dep %1.2f\n ",
fTrack,fTrackEta,fPDGCode,fTrackTOF,fPrimKinEn,
// Setters
virtual void SetVolume(Int_t i, Int_t val) {fVolume[i]=val;}
virtual void SetPDGCode(Int_t code) {fPDGCode=code;}
- virtual void SetMotherPDGCode(Int_t code) {fPDGCode=code;}
+ virtual void SetMotherPDGCode(Int_t code) {fMotherPDGCode=code;}
virtual void SetLightPMQ(Float_t value) {fLightPMQ=value;}
virtual void SetLightPMC(Float_t value) {fLightPMC=value;}
virtual void SetSFlag(Float_t value) {fSFlag=value;}