,fHistGammasFoundOutOfAccAltCent(0)
,fHistGammasGeneratedCent(0)
,fHistGammasFoundRecoEnergyCent(0)
+ ,fHistAllGammasFoundRecoEnergyCent(0)
,fHistGammasFoundOutOfAccRecoEnergyCent(0)
+ ,fHistAllGammasFoundOutOfAccRecoEnergyCent(0)
,fHistChargedTracksCut(0)
,fHistChargedTracksAccepted(0)
,fHistGammasCut(0)
delete fHistGammasFoundAltCent; // enter comment here
delete fHistGammasGeneratedCent; // enter comment here
delete fHistGammasFoundRecoEnergyCent;
+ delete fHistAllGammasFoundRecoEnergyCent;
delete fHistGammasFoundOutOfAccRecoEnergyCent;
+ delete fHistAllGammasFoundOutOfAccRecoEnergyCent;
delete fHistChargedTracksCut;
delete fHistChargedTracksAccepted;
delete fHistGammasCut;
{
//Printf("Particle with eta: %f, pid: %d", part->Eta(), code);
// calculate E_T
- if (
+
+ //fHistGammasFoundRecoEnergyCent->Fill(fReconstructedE,fCentClass);
+ //cout<<"Messed up filling would have filled with "<<fReconstructedE<<" cent "<<fCentClass<<endl;
+ if (
TMath::Abs(code) == fgProtonCode ||
TMath::Abs(code) == fgNeutronCode ||
TMath::Abs(code) == fgLambdaCode ||
if(!written && (fDepositedCode==fgGammaCode || fDepositedCode==fgEMinusCode || fDepositedCode ==fgEPlusCode)){//if the particle hitting the calorimeter is gamma, electron and not from a kaon
// cout<<"Signal ";
// PrintFamilyTreeShort(iPart, stack);
+//efficiencies without cutting track matched gammas
+ if(containsGamma){
+ if(!fSelector->CutGeometricalAcceptance(*primPart)){
+ fHistAllGammasFoundOutOfAccRecoEnergyCent->Fill(fReconstructedE,fCentClass);
+ }
+ else{
+ fHistAllGammasFoundRecoEnergyCent->Fill(fReconstructedE,fCentClass);
+ }
+ //cout<<"contains gamma"<<endl;
+ }
+ //else{cout<<"does not contain gamma"<<endl;}
+
fClusterMultGammas++;
written = kTRUE;
fHistGammasFoundOutOfAccRecoEnergyCent->Fill(fReconstructedE,fCentClass);
}
else{
+ //cout<<"filling reference histogram"<<endl;
if(gammaEnergy>0) fHistGammaEnergyCrossCheckAlt->Fill(gammaEnergy,(fReconstructedEt-gammaEnergy)/fReconstructedEt);
fHistGammasFoundRecoEnergyCent->Fill(fReconstructedE,fCentClass);
etGammaCrossCheckAlt += clEt;
fHistGammasFoundCent = new TH2F("fHistGammasFoundCent", "fHistGammasFoundCent",200, 0, 10,20,-0.5,19.5);
fHistGammasFoundOutOfAccAltCent = new TH2F("fHistGammasFoundOutOfAccAltCent", "fHistGammasFoundOutOfAccCent",200, 0, 10,20,-0.5,19.5);
fHistGammasFoundRecoEnergyCent = new TH2F("fHistGammasFoundRecoEnergyCent", "fHistGammasFoundRecoEnergyCent",200, 0, 10,20,-0.5,19.5);
+ fHistAllGammasFoundRecoEnergyCent = new TH2F("fHistAllGammasFoundRecoEnergyCent", "fHistAllGammasFoundRecoEnergyCent",200, 0, 10,20,-0.5,19.5);
fHistGammasFoundOutOfAccRecoEnergyCent = new TH2F("fHistGammasFoundOutOfAccRecoEnergyCent", "fHistGammasFoundOutOfAccRecoEnergyCent",200, 0, 10,20,-0.5,19.5);
+ fHistAllGammasFoundOutOfAccRecoEnergyCent = new TH2F("fHistAllGammasFoundOutOfAccRecoEnergyCent", "fHistAllGammasFoundOutOfAccRecoEnergyCent",200, 0, 10,20,-0.5,19.5);
fHistGammasFoundAltCent = new TH2F("fHistGammasFoundAltCent", "fHistGammasFoundAltCent",200, 0, 10,20,-0.5,19.5);
fHistGammasGeneratedCent = new TH2F("fHistGammasGeneratedCent", "fHistGammasGeneratedCent",200, 0, 10,20,-0.5,19.5);
fHistChargedTracksCut = new TH1F("fHistChargedTracksCut", "fHistChargedTracksCut",100, 0, 5);
list->Add(fHistGammasFoundAltCent);
list->Add(fHistGammasGeneratedCent);
list->Add(fHistGammasFoundRecoEnergyCent);
+ list->Add(fHistAllGammasFoundRecoEnergyCent);
list->Add(fHistGammasFoundOutOfAccRecoEnergyCent);
+ list->Add(fHistAllGammasFoundOutOfAccRecoEnergyCent);
list->Add(fHistChargedTracksCut);
list->Add(fHistChargedTracksAccepted);
list->Add(fHistGammasCut);