MadCluAnalysis Class Reference

#include <MadCluAnalysis.h>

Inheritance diagram for MadCluAnalysis:

List of all members.

Public Member Functions
	MadCluAnalysis (TChain *chainSR=0, TChain *chainMC=0, TChain *chainTH=0, TChain *chainEM=0)
	MadCluAnalysis (JobC *, string, int)
	~MadCluAnalysis ()
void	Plot (char *fileName, Int_t startEnt=0)
void	Plotter (char *fileName)
Double_t	FOM ()
void	DataDistributions (char *filename)
void	DataMCComp (char *, char *)
Bool_t	PassDataCleaningCuts (Int_t shower)
Bool_t	IsTrackInSlice (Int_t slice)
Bool_t	IsNueFidVtxEvt (Int_t event)
Double_t *	PHWProbEM (Int_t event)
Double_t *	PHWCluID (Int_t event)
Double_t *	PHWAvgDev (Int_t event)
Double_t *	ChargeFracRMS (Int_t event, Int_t method=0)
Double_t *	AveNumSSPlane (Int_t event)
Bool_t	FillPHWProbEM (Int_t event)
Bool_t	FillPHWCluID (Int_t event)
Bool_t	FillPHWAvgDev (Int_t event)
Bool_t	FillChargeFracRMS (Int_t event, Int_t method=0)
Bool_t	FillAveNumSSPlane (Int_t event)
Bool_t	FillNCluster (Int_t event)
Bool_t	FillBestProbEM (Int_t event)
Bool_t	ReadPIDFile (char *filename)
Bool_t	ReadNNFile (char *filename)
Public Attributes
Double_t	PAN_PHWCluIDU
Double_t	PAN_PHWCluIDV
Double_t	PAN_PHWCluIDBest
Double_t	PAN_PHWProbEMU
Double_t	PAN_PHWProbEMV
Double_t	PAN_PHWProbEMBest
Double_t	PAN_PHWAvgDevU
Double_t	PAN_PHWAvgDevV
Double_t	PAN_PHWAvgDevBest
Double_t	PAN_ChargeFracRMSU
Double_t	PAN_ChargeFracRMSV
Double_t	PAN_ChargeFracRMSBest
Double_t	PAN_AveNumSSPlaneU
Double_t	PAN_AveNumSSPlaneV
Double_t	PAN_AveNumSSPlaneBest
Int_t	PAN_NShowerStripU
Int_t	PAN_NShowerStripV
Int_t	PAN_NShowerStripBest
Int_t	PAN_NPhysShowerStripU
Int_t	PAN_NPhysShowerStripV
Double_t	PAN_BestProbEM
Double_t	PAN_RatioProbEM
Double_t	PAN_Weight
Float_t	fPID3
Float_t	fPID4
Float_t	fPID5
Int_t	fis_nue_pid
Int_t	PAN_NClusterU
Int_t	PAN_NClusterV
Int_t	PAN_NPhysClusterU
Int_t	PAN_NPhysClusterV
TH1F *	flikelihoodDists [12]
TMultiLayerPerceptron *	fneuralNet
Protected Member Functions
Bool_t	PassTruthSignal (Int_t mcevent=0)
Bool_t	PassAnalysisCuts (Int_t event=0)
Bool_t	PassBasicCuts (Int_t event=0)
Float_t	PID (Int_t event=0, Int_t method=0)
Float_t	RecoAnalysisEnergy (Int_t event=0)
Float_t	GetWeight (Int_t event=0)
Bool_t	AddUserBranches (TTree *)
void	CallUserFunctions (Int_t)

---

Constructor & Destructor Documentation

MadCluAnalysis::MadCluAnalysis (  TChain *  chainSR = 0, TChain *  chainMC = 0, TChain *  chainTH = 0, TChain *  chainEM = 0 )

Definition at line 36 of file MadCluAnalysis.cxx. References MadBase::Clear(), fis_nue_pid, fneuralNet, fPID3, fPID4, fPID5, MadBase::InitChain(), PAN_AveNumSSPlaneBest, PAN_AveNumSSPlaneU, PAN_AveNumSSPlaneV, PAN_BestProbEM, PAN_ChargeFracRMSBest, PAN_ChargeFracRMSU, PAN_ChargeFracRMSV, PAN_NClusterU, PAN_NClusterV, PAN_NPhysClusterU, PAN_NPhysClusterV, PAN_NPhysShowerStripU, PAN_NPhysShowerStripV, PAN_NShowerStripBest, PAN_NShowerStripU, PAN_NShowerStripV, PAN_PHWAvgDevBest, PAN_PHWAvgDevU, PAN_PHWAvgDevV, PAN_PHWCluIDBest, PAN_PHWCluIDU, PAN_PHWCluIDV, PAN_PHWProbEMBest, PAN_PHWProbEMU, PAN_PHWProbEMV, PAN_RatioProbEM, and PAN_Weight. { if(!chainSR) { record = 0; mcrecord = 0; threcord = 0; emrecord = 0; Clear(); whichSource = -1; isMC = true; isTH = true; isEM = true; Nentries = -1; return; } InitChain(chainSR,chainMC,chainTH,chainEM); whichSource = 0; PAN_PHWCluIDU = 0; PAN_PHWCluIDV = 0; PAN_PHWCluIDBest = 0; PAN_PHWProbEMU = 0; PAN_PHWProbEMV = 0; PAN_PHWProbEMBest = 0; PAN_PHWAvgDevU = 0; PAN_PHWAvgDevV = 0; PAN_PHWAvgDevBest = 0; PAN_ChargeFracRMSU = 0; PAN_ChargeFracRMSV = 0; PAN_ChargeFracRMSBest = 0; PAN_AveNumSSPlaneU = 0; PAN_AveNumSSPlaneV = 0; PAN_AveNumSSPlaneBest = 0; PAN_NShowerStripU = 0; PAN_NShowerStripV = 0; PAN_NShowerStripBest = 0; PAN_NPhysShowerStripU = 0; PAN_NPhysShowerStripV = 0; PAN_BestProbEM = 0; PAN_RatioProbEM = 0; PAN_Weight = 0; fPID3 = 0; fPID4 = 0; fPID5 = 0; fis_nue_pid = 0; PAN_NClusterU = 0; PAN_NClusterV = 0; PAN_NPhysClusterU = 0; PAN_NPhysClusterV = 0; fneuralNet = NULL; }

MadCluAnalysis::MadCluAnalysis (  JobC *  , string  , int  )

Definition at line 100 of file MadCluAnalysis.cxx. References MadBase::Clear(), fis_nue_pid, fneuralNet, fPID3, fPID4, fPID5, PAN_AveNumSSPlaneBest, PAN_AveNumSSPlaneU, PAN_AveNumSSPlaneV, PAN_BestProbEM, PAN_ChargeFracRMSBest, PAN_ChargeFracRMSU, PAN_ChargeFracRMSV, PAN_NClusterU, PAN_NClusterV, PAN_NPhysClusterU, PAN_NPhysClusterV, PAN_NPhysShowerStripU, PAN_NPhysShowerStripV, PAN_NShowerStripBest, PAN_NShowerStripU, PAN_NShowerStripV, PAN_PHWAvgDevBest, PAN_PHWAvgDevU, PAN_PHWAvgDevV, PAN_PHWCluIDBest, PAN_PHWCluIDU, PAN_PHWCluIDV, PAN_PHWProbEMBest, PAN_PHWProbEMU, PAN_PHWProbEMV, PAN_RatioProbEM, and PAN_Weight. { record = 0; mcrecord = 0; threcord = 0; emrecord = 0; Clear(); isMC = true; isTH = true; isEM = true; Nentries = entries; whichSource = 1; jcPath = path; fJC = j; fChain = NULL; PAN_PHWCluIDU = 0; PAN_PHWCluIDV = 0; PAN_PHWCluIDBest = 0; PAN_PHWProbEMU = 0; PAN_PHWProbEMV = 0; PAN_PHWProbEMBest = 0; PAN_PHWAvgDevU = 0; PAN_PHWAvgDevV = 0; PAN_PHWAvgDevBest = 0; PAN_ChargeFracRMSU = 0; PAN_ChargeFracRMSV = 0; PAN_ChargeFracRMSBest = 0; PAN_AveNumSSPlaneU = 0; PAN_AveNumSSPlaneV = 0; PAN_AveNumSSPlaneBest = 0; PAN_NShowerStripU = 0; PAN_NShowerStripV = 0; PAN_NShowerStripBest = 0; PAN_NPhysShowerStripU = 0; PAN_NPhysShowerStripV = 0; PAN_BestProbEM = 0; PAN_RatioProbEM = 0; PAN_Weight = 0; fPID3 = 0; fPID4 = 0; fPID5 = 0; fis_nue_pid = 0; PAN_NClusterU = 0; PAN_NClusterV = 0; PAN_NPhysClusterU = 0; PAN_NPhysClusterV = 0; fneuralNet = NULL; }

MadCluAnalysis::~MadCluAnalysis (   )

Definition at line 162 of file MadCluAnalysis.cxx. { }

---

Member Function Documentation

Bool_t MadCluAnalysis::AddUserBranches (  TTree *   )  [protected, virtual]

Reimplemented from MadAnalysis. Definition at line 2712 of file MadCluAnalysis.cxx. References fis_nue_pid, fPID3, fPID4, fPID5, PAN_AveNumSSPlaneBest, PAN_AveNumSSPlaneU, PAN_AveNumSSPlaneV, PAN_BestProbEM, PAN_ChargeFracRMSBest, PAN_ChargeFracRMSU, PAN_ChargeFracRMSV, PAN_NClusterU, PAN_NClusterV, PAN_NPhysClusterU, PAN_NPhysClusterV, PAN_NPhysShowerStripU, PAN_NPhysShowerStripV, PAN_NShowerStripBest, PAN_NShowerStripU, PAN_NShowerStripV, PAN_PHWAvgDevBest, PAN_PHWAvgDevU, PAN_PHWAvgDevV, PAN_PHWCluIDBest, PAN_PHWCluIDU, PAN_PHWCluIDV, PAN_PHWProbEMBest, PAN_PHWProbEMU, PAN_PHWProbEMV, PAN_RatioProbEM, and PAN_Weight. { if(!tree) return false; tree->Branch("ss_PhwIDU",&PAN_PHWCluIDU,"ss_PhwIDU/D",32000); tree->Branch("ss_PhwIDV",&PAN_PHWCluIDV,"ss_PhwIDV/D",32000); tree->Branch("ss_PhwIDBest",&PAN_PHWCluIDBest,"ss_PhwIDBest/D",32000); tree->Branch("ss_PhwProbEMU",&PAN_PHWProbEMU,"ss_PhwProbEMU/D",32000); tree->Branch("ss_PhwProbEMV",&PAN_PHWProbEMV,"ss_PhwProbEMV/D",32000); tree->Branch("ss_PhwProbEMBest",&PAN_PHWProbEMBest,"ss_PhwProbEMBest/D",32000); tree->Branch("ss_PHWAvgDevU",&PAN_PHWAvgDevU,"ss_PHWAvgDevU/D",32000); tree->Branch("ss_PHWAvgDevV",&PAN_PHWAvgDevV,"ss_PHWAvgDevV/D",32000); tree->Branch("ss_PHWAvgDevBest",&PAN_PHWAvgDevBest,"ss_PHWAvgDevBest/D",32000); tree->Branch("ss_ChargeFracRMSU",&PAN_ChargeFracRMSU,"ss_ChargeFracRMSU/D",32000); tree->Branch("ss_ChargeFracRMSV",&PAN_ChargeFracRMSV,"ss_ChargeFracRMSV/D",32000); tree->Branch("ss_ChargeFracRMSBest",&PAN_ChargeFracRMSBest,"ss_ChargeFracRMSBest/D",32000); tree->Branch("ss_AveNumSSPlaneU",&PAN_AveNumSSPlaneU,"ss_AveNumSSPlaneU/D",32000); tree->Branch("ss_AveNumSSPlaneV",&PAN_AveNumSSPlaneV,"ss_AveNumSSPlaneV/D",32000); tree->Branch("ss_AveNumSSPlaneBest",&PAN_AveNumSSPlaneBest,"ss_AveNumSSPlaneBest/D",32000); tree->Branch("ss_NShowerStripU",&PAN_NShowerStripU,"ss_NShowerStripU/I",32000); tree->Branch("ss_NShowerStripV",&PAN_NShowerStripV,"ss_NShowerStripV/I",32000); tree->Branch("ss_NShowerStripBest",&PAN_NShowerStripBest,"ss_NShowerStripBest/I",32000); tree->Branch("ss_NPhysShowerStripU",&PAN_NPhysShowerStripU,"ss_NPhysShowerStripU/I",32000); tree->Branch("ss_NPhysShowerStripV",&PAN_NPhysShowerStripV,"ss_NPhysShowerStripV/I",32000); tree->Branch("ss_BestProbEM",&PAN_BestProbEM,"ss_BestProbEM/D",32000); tree->Branch("ss_RatioProbEM",&PAN_RatioProbEM,"ss_RatioProbEM/D",32000); tree->Branch("ss_Weight",&PAN_Weight,"ss_Weight/D",32000); tree->Branch("pid3",&fPID3,"pid3/F",32000); tree->Branch("pid4",&fPID4,"pid4/F",32000); tree->Branch("pid5",&fPID5,"pid5/F",32000); tree->Branch("is_nue_fid",&fis_nue_pid,"is_nue_fid/I",32000); tree->Branch("NClusterU",&PAN_NClusterU,"NClusterU/I",32000); tree->Branch("NClusterV",&PAN_NClusterV,"NClusterV/I",32000); tree->Branch("NPhysClusterU",&PAN_NPhysClusterU,"NPhysClusterU/I",32000); tree->Branch("NPhysClusterV",&PAN_NPhysClusterV,"NPhysClusterV/I",32000); return true; }

Double_t * MadCluAnalysis::AveNumSSPlane (  Int_t  event  )

Definition at line 2664 of file MadCluAnalysis.cxx. References NtpSRPlane::beg, NtpSRShower::clu, NtpSRPlane::end, NtpSRCluster::id, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), MadBase::LoadStrip(), NtpSRShower::ncluster, NtpSRCluster::nstrip, NtpSRPlane::nu, NtpSRPlane::nv, NtpSRStrip::plane, NtpSRShower::plane, NtpSRCluster::planeview, and NtpSRCluster::stp. Referenced by FillAveNumSSPlane(). { Double_t *aveNumSSPlane = new Double_t[2]; aveNumSSPlane[0] = 0; aveNumSSPlane[1] = 0; if(!LoadEvent(event)) { aveNumSSPlane[0] = -10.; aveNumSSPlane[1] = -10.; return aveNumSSPlane; } Int_t shower = -1; if(!LoadLargestShowerFromEvent(event,shower)) { aveNumSSPlane[0] = -10.; aveNumSSPlane[1] = -10.; return aveNumSSPlane; } Int_t *clusters = ntpShower->clu; for(int k=ntpShower->plane.beg; k<=ntpShower->plane.end; k++){ Double_t planeAveU = 0; Double_t planeAveV = 0; for(int l=0;l<ntpShower->ncluster;l++){ Int_t index = clusters[l]; if(!LoadCluster(index)) continue; if(!(ntpCluster->id==0||ntpCluster->id==1||ntpCluster->id==3)) continue; Int_t *strip = ntpCluster->stp; for(int m=0;m<ntpCluster->nstrip;m++){ if(!LoadStrip(strip[m])) continue; if(ntpStrip->plane==k){ if(ntpCluster->planeview==2) planeAveU += 1; else if(ntpCluster->planeview==3) planeAveV += 1; break; } } } aveNumSSPlane[0] += planeAveU; aveNumSSPlane[1] += planeAveV; } aveNumSSPlane[0]/=ntpShower->plane.nu; aveNumSSPlane[1]/=ntpShower->plane.nv; return aveNumSSPlane; }

void MadCluAnalysis::CallUserFunctions (  Int_t   )  [protected, virtual]

Reimplemented from MadAnalysis. Definition at line 2757 of file MadCluAnalysis.cxx. References abs(), FillAveNumSSPlane(), FillBestProbEM(), FillChargeFracRMS(), FillNCluster(), FillPHWAvgDev(), FillPHWCluID(), FillPHWProbEM(), fis_nue_pid, fPID3, fPID4, fPID5, MadQuantities::GetNShowerStrip(), NtpMCTruth::iaction, NtpSREvent::index, NtpMCTruth::inu, NtpMCTruth::inunoosc, IsNueFidVtxEvt(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), MadBase::LoadTruthForRecoTH(), max, min, PAN_AveNumSSPlaneBest, PAN_AveNumSSPlaneU, PAN_AveNumSSPlaneV, PAN_BestProbEM, PAN_ChargeFracRMSBest, PAN_ChargeFracRMSU, PAN_ChargeFracRMSV, PAN_NClusterU, PAN_NClusterV, PAN_NPhysClusterU, PAN_NPhysClusterV, PAN_NPhysShowerStripU, PAN_NPhysShowerStripV, PAN_NShowerStripBest, PAN_NShowerStripU, PAN_NShowerStripV, PAN_PHWAvgDevBest, PAN_PHWAvgDevU, PAN_PHWAvgDevV, PAN_PHWCluIDBest, PAN_PHWCluIDU, PAN_PHWCluIDV, PAN_PHWProbEMBest, PAN_PHWProbEMU, PAN_PHWProbEMV, PAN_RatioProbEM, PAN_Weight, and PID(). { PAN_PHWCluIDU = 0; PAN_PHWCluIDV = 0; PAN_PHWCluIDBest = 0; PAN_PHWProbEMU = 0; PAN_PHWProbEMV = 0; PAN_PHWProbEMBest = 0; PAN_PHWAvgDevU = 0; PAN_PHWAvgDevV = 0; PAN_PHWAvgDevBest = 0; PAN_ChargeFracRMSU = 0; PAN_ChargeFracRMSV = 0; PAN_ChargeFracRMSBest = 0; PAN_AveNumSSPlaneU = 0; PAN_AveNumSSPlaneV = 0; PAN_AveNumSSPlaneBest = 0; PAN_NShowerStripU = 0; PAN_NShowerStripV = 0; PAN_NShowerStripBest = 0; PAN_NPhysShowerStripU = 0; PAN_NPhysShowerStripV = 0; PAN_BestProbEM = 0; PAN_RatioProbEM = 0; PAN_Weight = 0; fPID3 = 0; fPID4 = 0; fPID5 = 0; fis_nue_pid = 0; PAN_NClusterU = 0; PAN_NClusterV = 0; PAN_NPhysClusterU = 0; PAN_NPhysClusterV = 0; if(!LoadEvent(event)) return; fis_nue_pid = this->IsNueFidVtxEvt(event); this->FillNCluster(event); this->FillBestProbEM(event); this->FillPHWCluID(event); this->FillPHWProbEM(event); this->FillPHWAvgDev(event); this->FillChargeFracRMS(event,1); this->FillAveNumSSPlane(event); //set ratio if(PAN_PHWProbEMU>0 && PAN_PHWProbEMV>0) { if(PAN_PHWProbEMU>PAN_PHWProbEMV) PAN_RatioProbEM = PAN_PHWProbEMV/PAN_PHWProbEMU; else PAN_RatioProbEM = PAN_PHWProbEMU/PAN_PHWProbEMV; } Int_t shower = -1; LoadLargestShowerFromEvent(event,shower); Int_t *nshowerstrip = this->GetNShowerStrip(shower); PAN_NShowerStripU = nshowerstrip[0]; PAN_NShowerStripV = nshowerstrip[1]; delete [] nshowerstrip; Int_t *nphysshowerstrip = this->GetNShowerStrip(shower,1); PAN_NPhysShowerStripU = nphysshowerstrip[0]; PAN_NPhysShowerStripV = nphysshowerstrip[1]; delete [] nphysshowerstrip; double prob = 1.; Int_t mcevent = -1; if(LoadTruthForRecoTH(ntpEvent->index,mcevent)) { Double_t nNuMuFiles = 40; Double_t nNuEFiles = 39; Double_t nNuTauFiles = 39; Double_t POT = 7.4; //prob = Oscillate(ntpTruth->inu,ntpTruth->inunoosc, // ntpTruth->p4neu[3],735.,0.002175, // TMath::ASin(TMath::Sqrt(0.925))/2., // 0.01); //prob = Oscillate(ntpTruth->inu,ntpTruth->inunoosc, // ntpTruth->p4neu[3],735.,0.0025, // TMath::ASin(TMath::Sqrt(1.0))/2., // 0.01); if(ntpTruth->iaction==1){ if(abs(ntpTruth->inunoosc)==12 && abs(ntpTruth->inu)==12) PAN_Weight = prob*(POT/((nNuEFiles+nNuMuFiles)*6.5)); else if(abs(ntpTruth->inu)==12) PAN_Weight = prob*(POT/(nNuEFiles*6.5)); else if(abs(ntpTruth->inu)==14) PAN_Weight = prob*(POT/(nNuMuFiles*6.5)); else if(abs(ntpTruth->inu)==16) PAN_Weight = prob*(POT/(nNuTauFiles*6.5)); } else if(ntpTruth->iaction==0) { if(abs(ntpTruth->inu)==12) PAN_Weight = prob*(POT/((nNuEFiles)*6.5)); else if(abs(ntpTruth->inu)==14) PAN_Weight = prob*(POT/((nNuMuFiles)*6.5)); else if(abs(ntpTruth->inu)==16) PAN_Weight = prob*(POT/((nNuTauFiles)*6.5)); } } else PAN_Weight = 1.; fPID3 = PID(event,3); fPID4 = PID(event,4); fPID5 = PID(event,5); PAN_ChargeFracRMSBest = max(PAN_ChargeFracRMSU,PAN_ChargeFracRMSV); PAN_PHWProbEMBest = max(PAN_PHWProbEMU,PAN_PHWProbEMV); PAN_PHWCluIDBest = min(PAN_PHWCluIDU,PAN_PHWCluIDV); PAN_NShowerStripBest = min(PAN_NShowerStripU,PAN_NShowerStripV); PAN_AveNumSSPlaneBest = min(PAN_AveNumSSPlaneU,PAN_AveNumSSPlaneV); PAN_PHWAvgDevBest = min(fabs(PAN_PHWAvgDevU-0.02), fabs(PAN_PHWAvgDevV-0.02)); }

Double_t * MadCluAnalysis::ChargeFracRMS (  Int_t  event, Int_t  method = 0 )

Definition at line 2610 of file MadCluAnalysis.cxx. References NtpSRShower::clu, NtpSRStripPulseHeight::gev, NtpSRCluster::id, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), NtpSRShower::ncluster, NtpSRShower::nUcluster, NtpSRShower::nVcluster, NtpSRCluster::ph, and NtpSRCluster::planeview. Referenced by FillChargeFracRMS(). { Double_t *chargeFracRMS = new Double_t[2]; chargeFracRMS[0] = 0; chargeFracRMS[1] = 0; if(!LoadEvent(event)) { chargeFracRMS[0] = -10.; chargeFracRMS[1] = -10.; return chargeFracRMS; } Int_t shower = -1; if(!LoadLargestShowerFromEvent(event,shower)) { chargeFracRMS[0] = -10.; chargeFracRMS[1] = -10.; return chargeFracRMS; } Double_t totE_U = 0; Double_t totE_V = 0; Int_t *clusters = ntpShower->clu; for(int k=0; k<ntpShower->ncluster; k++){ Int_t index = clusters[k]; if(!LoadCluster(index)) continue; if(method==1 && !(ntpCluster->id==0 || ntpCluster->id==1 || ntpCluster->id==3)) continue; if(ntpCluster->planeview==2){ chargeFracRMS[0] += ntpCluster->ph.gev*ntpCluster->ph.gev; totE_U += ntpCluster->ph.gev; } else if(ntpCluster->planeview==3){ chargeFracRMS[1] += ntpCluster->ph.gev*ntpCluster->ph.gev; totE_V += ntpCluster->ph.gev; } } if(totE_U!=0) { chargeFracRMS[0] /= (totE_U*totE_U*ntpShower->nUcluster); chargeFracRMS[0] -= TMath::Power(1./ntpShower->nUcluster,2); chargeFracRMS[0] = TMath::Sqrt(chargeFracRMS[0]); } if(totE_V!=0) { chargeFracRMS[1] /= (totE_V*totE_V*ntpShower->nVcluster); chargeFracRMS[1] -= TMath::Power(1./ntpShower->nVcluster,2); chargeFracRMS[1] = TMath::Sqrt(chargeFracRMS[1]); } return chargeFracRMS; }

void MadCluAnalysis::DataDistributions (  char *  filename  )

Definition at line 866 of file MadCluAnalysis.cxx. References NtpSRShower::clu, NtpSREventSummary::date, VldContext::GetDetector(), MadBase::GetEntry(), SpillInfoBlock::GetHorn(), SpillInfoBlock::GetHpos(), SpillInfoBlock::GetMagnet(), VldTimeStamp::GetNanoSec(), SpillInfoBlock::GetPot(), VldTimeStamp::GetSec(), VldContext::GetSimFlag(), SpillInfo::GetSpillInfo(), SpillInfoBlock::GetTgtpos(), VldContext::GetTimeStamp(), RecHeader::GetVldContext(), SpillInfoBlock::GetVpos(), NtpSRStripPulseHeight::gev, NtpMCTruth::iaction, NtpSRCluster::id, NtpSREvent::index, NtpMCTruth::inu, NtpMCTruth::inunoosc, MadQuantities::IsFidVtxEvt(), MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), MadBase::LoadLargestTrackFromEvent(), MadBase::LoadTruthForRecoTH(), NtpSRDate::month, month, NtpSRShower::ncluster, NtpSREventSummary::nevent, NtpSREvent::nshower, PassDataCleaningCuts(), MadQuantities::PassTrackCuts(), NtpSRShower::ph, NtpSREvent::ph, NtpSRCluster::ph, NtpSRCluster::planeview, NtpSRCluster::probem, MadQuantities::RecoEMFrac(), MadQuantities::RecoMuEnergy(), and NtpSRDate::year. { SpillInfo pppinfo; Double_t datapot = 0; char savename[256]; sprintf(savename,"DataClu_%s.root",filename); TFile *save = new TFile(savename,"RECREATE"); save->cd(); string SplitName[6] = {"beamnue","nue","numu","nutau","nc","unknown"}; TH1F *IDHistTrack = new TH1F("IDHistTrack","Cluster ID for Track Events",6,0,5); IDHistTrack->SetXTitle("Cluster ID"); TH1F *IDHistShower = new TH1F("IDHistShower","Cluster ID for Shower Events",6,0,5); IDHistShower->SetXTitle("Cluster ID"); TH1F *IDHistTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",IDHistTrack->GetName(),SplitName[i].c_str()); IDHistTrackSplit[i] = new TH1F(name,"Cluster ID for Track Events",6,0,5); } TH1F *IDHistShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",IDHistShower->GetName(),SplitName[i].c_str()); IDHistShowerSplit[i] = new TH1F(name,"Cluster ID for Shower Events",6,0,5); } TH1F *RecoyTrack = new TH1F("RecoyTrack","Reco y Distribution for Track Events", 101,-0.005,1.005); RecoyTrack->SetXTitle("Reco y"); TH1F *RecoyShower = new TH1F("RecoyShower","Reco y Distribution for Shower Events", 101,-0.005,1.005); RecoyShower->SetXTitle("Reco y"); TH1F *RecoyTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoyTrack->GetName(),SplitName[i].c_str()); RecoyTrackSplit[i] = new TH1F(name,"Reco y Distribution for Track Events", 101,-0.005,1.005); } TH1F *RecoyShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoyShower->GetName(),SplitName[i].c_str()); RecoyShowerSplit[i] = new TH1F(name,"Reco y Distribution for Shower Events", 101,-0.005,1.005); } TH1F *RecoEMFracTrack = new TH1F("RecoEMFracTrack","Reco EMFrac for Track Events", 101,-0.005,1.005); RecoEMFracTrack->SetXTitle("Reco EMFrac"); TH1F *RecoEMFracShower = new TH1F("RecoEMFracShower","Reco EMFrac for Shower Events", 101,-0.005,1.005); RecoEMFracShower->SetXTitle("Reco EMFrac"); TH1F *RecoEMFracTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoEMFracTrack->GetName(),SplitName[i].c_str()); RecoEMFracTrackSplit[i] = new TH1F(name,"Reco EMFrac for Track Events", 101,-0.005,1.005); } TH1F *RecoEMFracShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoEMFracShower->GetName(),SplitName[i].c_str()); RecoEMFracShowerSplit[i] = new TH1F(name,"Reco EMFrac for Shower Events", 101,-0.005,1.005); } TH1F *RecoEMFracProbTrack = new TH1F("RecoEMFracProbTrack","Reco EMFrac for Track Events", 101,-0.005,1.005); RecoEMFracProbTrack->SetXTitle("Reco EMFrac"); TH1F *RecoEMFracProbShower = new TH1F("RecoEMFracProbShower","Reco EMFrac for Shower Events", 101,-0.005,1.005); RecoEMFracProbShower->SetXTitle("Reco EMFrac"); TH1F *RecoEMFracProbTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoEMFracProbTrack->GetName(),SplitName[i].c_str()); RecoEMFracProbTrackSplit[i] = new TH1F(name,"Reco EMFrac for Track Events", 101,-0.005,1.005); } TH1F *RecoEMFracProbShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoEMFracProbShower->GetName(),SplitName[i].c_str()); RecoEMFracProbShowerSplit[i] = new TH1F(name,"Reco EMFrac for Shower Events", 101,-0.005,1.005); } //primary cluster IDs for NC and CC TH1F *UClusterIDTrack = new TH1F("UClusterIDTrack","U View Primary Cluster ID for Track Events",6,-0.5,5.5); UClusterIDTrack->SetXTitle("Cluster ID"); TH1F *UClusterIDShower = new TH1F("UClusterIDShower","U View Primary Cluster ID for Shower Events",6,-0.5,5.5); UClusterIDShower->SetXTitle("Cluster ID"); TH1F *VClusterIDTrack = new TH1F("VClusterIDTrack","V View Primary Cluster ID for Track Events",6,-0.5,5.5); VClusterIDTrack->SetXTitle("Cluster ID"); TH1F *VClusterIDShower = new TH1F("VClusterIDShower","V View Primary Cluster ID for Shower Events",6,-0.5,5.5); VClusterIDShower->SetXTitle("Cluster ID"); TH1F *UClusterIDTrackSplit[6] = {0}; TH1F *VClusterIDTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",UClusterIDTrack->GetName(),SplitName[i].c_str()); UClusterIDTrackSplit[i] = new TH1F(name,"U View Primary Cluster ID for Track Events",6,-0.5,5.5); sprintf(name,"%s_%s",VClusterIDTrack->GetName(),SplitName[i].c_str()); VClusterIDTrackSplit[i] = new TH1F(name,"V View Primary Cluster ID for Track Events",6,-0.5,5.5); } TH1F *UClusterIDShowerSplit[6] = {0}; TH1F *VClusterIDShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",UClusterIDShower->GetName(),SplitName[i].c_str()); UClusterIDShowerSplit[i] = new TH1F(name,"U View Primary Cluster ID for Shower Events",6,-0.5,5.5); sprintf(name,"%s_%s",VClusterIDShower->GetName(),SplitName[i].c_str()); VClusterIDShowerSplit[i] = new TH1F(name,"V View Primary Cluster ID for Shower Events",6,-0.5,5.5); } //primary cluster ProbEMs TH1F *UClusterProbEMTrack = new TH1F("UClusterProbEMTrack","U View Primary Cluster ProbEM for Track Events",100,0,1); UClusterProbEMTrack->SetXTitle("Cluster ProbEM"); TH1F *UClusterProbEMShower = new TH1F("UClusterProbEMShower","U View Primary Cluster ProbEM for Shower Events",100,0,1); UClusterProbEMShower->SetXTitle("Cluster ProbEM"); TH1F *VClusterProbEMTrack = new TH1F("VClusterProbEMTrack","V View Primary Cluster ProbEM for Track Events",100,0,1); VClusterProbEMTrack->SetXTitle("Cluster ProbEM"); TH1F *VClusterProbEMShower = new TH1F("VClusterProbEMShower","V View Primary Cluster ProbEM for Shower Events",100,0,1); VClusterProbEMShower->SetXTitle("Cluster ProbEM"); TH1F *UClusterProbEMTrackSplit[6] = {0}; TH1F *VClusterProbEMTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",UClusterProbEMTrack->GetName(),SplitName[i].c_str()); UClusterProbEMTrackSplit[i] = new TH1F(name,"U View Primary Cluster ProbEM for Track Events",100,0,1); sprintf(name,"%s_%s",VClusterProbEMTrack->GetName(),SplitName[i].c_str()); VClusterProbEMTrackSplit[i] = new TH1F(name,"V View Primary Cluster ProbEM for Track Events",100,0,1); } TH1F *UClusterProbEMShowerSplit[6] = {0}; TH1F *VClusterProbEMShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",UClusterProbEMShower->GetName(),SplitName[i].c_str()); UClusterProbEMShowerSplit[i] = new TH1F(name,"U View Primary Cluster ProbEM for Shower Events",100,0,1); sprintf(name,"%s_%s",VClusterProbEMShower->GetName(),SplitName[i].c_str()); VClusterProbEMShowerSplit[i] = new TH1F(name,"V View Primary Cluster ProbEM for Shower Events",100,0,1); } //Average #cluster with energy TH2F *AvgCluEnTrack = new TH2F("AvgCluEnTrack", "Average Number of Clusters with Reco Shower Energy", 100,0,25,20,0,20); AvgCluEnTrack->SetXTitle("Reco Shower Energy (GeV)"); AvgCluEnTrack->SetYTitle("# of Clusters"); TH2F *AvgCluEnShower = new TH2F("AvgCluEnShower", "Average Number of Clusters with Reco Shower Energy", 100,0,25,20,0,20); AvgCluEnShower->SetXTitle("Reco Shower Energy (GeV)"); AvgCluEnShower->SetYTitle("# of Clusters"); TH2F *AvgCluEnTracknox = new TH2F("AvgCluEnTracknox", "Average Number of Clusters with Reco Shower Energy (excl xtalk clu) - CC", 100,0,25,20,0,20); AvgCluEnTracknox->SetXTitle("Reco Shower Energy (GeV)"); AvgCluEnTracknox->SetYTitle("# of Clusters"); TH2F *AvgCluEnShowernox = new TH2F("AvgCluEnShowernox", "Average Number of Clusters with Reco Shower Energy (excl xtalk clu) - NC", 100,0,25,20,0,20); AvgCluEnShowernox->SetXTitle("Reco Shower Energy (GeV)"); AvgCluEnShowernox->SetYTitle("# of Clusters"); TH2F *AvgCluEnTrackSplit[6] = {0}; TH2F *AvgCluEnnoxTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",AvgCluEnTrack->GetName(),SplitName[i].c_str()); AvgCluEnTrackSplit[i] = new TH2F(name,"Average Number of Clusters with Reco Shower Energy", 100,0,25,20,0,20); sprintf(name,"%s_%s",AvgCluEnTracknox->GetName(),SplitName[i].c_str()); AvgCluEnnoxTrackSplit[i] = new TH2F(name, "Average Number of Clusters with Reco Shower Energy (excl xtalk clu)", 100,0,25,20,0,20); } TH2F *AvgCluEnShowerSplit[6] = {0}; TH2F *AvgCluEnnoxShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",AvgCluEnShower->GetName(),SplitName[i].c_str()); AvgCluEnShowerSplit[i] = new TH2F(name,"Average Number of Clusters with Reco Shower Energy", 100,0,25,20,0,20); sprintf(name,"%s_%s",AvgCluEnShowernox->GetName(),SplitName[i].c_str()); AvgCluEnnoxShowerSplit[i] = new TH2F(name, "Average Number of Clusters with Reco Shower Energy (excl xtalk clu)", 100,0,25,20,0,20); } TH1F *EMFracUVDiffTrack = new TH1F("EMFracUVDiffTrack","U/V Diff in Reco EMFrac",200,-2,2); EMFracUVDiffTrack->SetXTitle("U-V EMFrac"); TH1F *EMFracUVDiffShower = new TH1F("EMFracUVDiffShower","U/V Diff in Reco EMFrac",200,-2,2); EMFracUVDiffShower->SetXTitle("U-V EMFrac"); TH1F *EMFracUVDiffTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",EMFracUVDiffTrack->GetName(),SplitName[i].c_str()); EMFracUVDiffTrackSplit[i] = new TH1F(name,"U/V Diff in Reco EMFrac",200,-2,2); } TH1F *EMFracUVDiffShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",EMFracUVDiffShower->GetName(),SplitName[i].c_str()); EMFracUVDiffShowerSplit[i] = new TH1F(name,"U/V Diff in Reco EMFrac",200,-2,2); } for(int i=0;i<Nentries;i++){ if(i%1000==0) std::cout << float(i)*100./float(Nentries) << "% done" << std::endl; if(!this->GetEntry(i)) continue; //cout << ntpHeader->GetRun() << endl; if(eventSummary->nevent==0) continue; //if not MC: if(ntpHeader->GetVldContext().GetSimFlag()!=4) { Double_t snarltime = (ntpHeader->GetVldContext().GetTimeStamp().GetSec() + (ntpHeader->GetVldContext().GetTimeStamp().GetNanoSec()/1e9)); Int_t month = eventSummary->date.month; Int_t year = eventSummary->date.year; SpillInfoBlock spillinfo; pppinfo.GetSpillInfo(month,year,snarltime,spillinfo); if(spillinfo.GetTgtpos()>1000) continue; if(spillinfo.GetHorn()<16000) continue; if(spillinfo.GetHpos()<-2 || spillinfo.GetHpos()>-1.6) continue; if(spillinfo.GetVpos()<0.8 || spillinfo.GetVpos()>1.2) continue; if(spillinfo.GetMagnet()>-4800.) continue; datapot += spillinfo.GetPot(); } Int_t is_fid = 0; //fill tree once for each reconstructed event for(int j=0;j<eventSummary->nevent;j++){ if(!LoadEvent(j)) continue; //no event found if(ntpEvent->nshower==0) continue; //no shower found //get detector type: if(ntpHeader->GetVldContext().GetDetector()==Detector::kFar){ //fiducial criteria on vtx for far detector is_fid = 1; if(!IsFidVtxEvt(ntpEvent->index)) is_fid = 0; } else if(ntpHeader->GetVldContext().GetDetector()==Detector::kNear){ //fiducial criteria on vtx for near detector is_fid = 1; if(!IsFidVtxEvt(ntpEvent->index)) is_fid = 0; } if(is_fid==0) continue; Int_t shower = -1; if(!LoadLargestShowerFromEvent(ntpEvent->index,shower)) continue; if(!PassDataCleaningCuts(shower)) continue; Int_t noxClu = 0; double emfracU = 0; double emfracV = 0; double emfrac = RecoEMFrac(shower,0,emfracU,emfracV); double emfracprob = RecoEMFrac(shower,1,emfracU,emfracV); Int_t primaryUID = -1; Float_t primaryUprobem = 0; Double_t maxEnU = 0; Int_t primaryVID = -1; Float_t primaryVprobem = 0; Double_t maxEnV = 0; for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; if(!LoadCluster(index)) continue; if(ntpCluster->ph.gev>0.5 && (ntpCluster->id==0 || ntpCluster->id==1 || ntpCluster->id==3)) noxClu++; if(ntpCluster->planeview==2 && ntpCluster->ph.gev>maxEnU) { maxEnU = ntpCluster->ph.gev; primaryUID = ntpCluster->id; primaryUprobem = ntpCluster->probem; } if(ntpCluster->planeview==3 && ntpCluster->ph.gev>maxEnV) { maxEnV = ntpCluster->ph.gev; primaryVID = ntpCluster->id; primaryVprobem = ntpCluster->probem; } } Int_t track = -1; LoadLargestTrackFromEvent(ntpEvent->index,track); if(track>=0 && PassTrackCuts(track)) { //Fill track events: double y = 1. - RecoMuEnergy(0,track)/ntpEvent->ph.gev; for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; if(!LoadCluster(index)) continue; IDHistTrack->Fill(ntpCluster->id); } RecoyTrack->Fill(y); RecoEMFracTrack->Fill(emfrac); RecoEMFracProbTrack->Fill(emfracprob); UClusterIDTrack->Fill(primaryUID); VClusterIDTrack->Fill(primaryVID); if(primaryUID==0) UClusterProbEMTrack->Fill(primaryUprobem); if(primaryVID==0) VClusterProbEMTrack->Fill(primaryVprobem); AvgCluEnTrack->Fill(ntpShower->ph.gev,ntpShower->ncluster); AvgCluEnTracknox->Fill(ntpShower->ph.gev,noxClu); EMFracUVDiffTrack->Fill(emfracU-emfracV); Int_t mcevent = -1; Int_t hist_index = 5; if(LoadTruthForRecoTH(ntpEvent->index,mcevent)){ if(ntpTruth->iaction==1){ if(ntpTruth->inu==12 && ntpTruth->inunoosc==12) hist_index = 0; if(ntpTruth->inu==12 && ntpTruth->inunoosc==14) hist_index = 1; if(ntpTruth->inu==14 && ntpTruth->inunoosc==14) hist_index = 2; if(ntpTruth->inu==16 && ntpTruth->inunoosc==14) hist_index = 3; } else hist_index = 4; } for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; if(!LoadCluster(index)) continue; IDHistTrackSplit[hist_index]->Fill(ntpCluster->id); } RecoyTrackSplit[hist_index]->Fill(y); RecoEMFracTrackSplit[hist_index]->Fill(emfrac); RecoEMFracProbTrackSplit[hist_index]->Fill(emfracprob); UClusterIDTrackSplit[hist_index]->Fill(primaryUID); VClusterIDTrackSplit[hist_index]->Fill(primaryVID); if(primaryUID==0) UClusterProbEMTrackSplit[hist_index]->Fill(primaryUprobem); if(primaryVID==0) VClusterProbEMTrackSplit[hist_index]->Fill(primaryVprobem); AvgCluEnTrackSplit[hist_index]->Fill(ntpShower->ph.gev,ntpShower->ncluster); AvgCluEnnoxTrackSplit[hist_index]->Fill(ntpShower->ph.gev,noxClu); EMFracUVDiffTrackSplit[hist_index]->Fill(emfracU-emfracV); } else { //Fill shower events: for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; if(!LoadCluster(index)) continue; IDHistShower->Fill(ntpCluster->id); } if(primaryVprobem>0.2 && primaryUprobem>0.2 && primaryUID==0 && primaryVID==0) RecoyShower->Fill(emfracprob); RecoEMFracShower->Fill(emfrac); RecoEMFracProbShower->Fill(emfracprob); UClusterIDShower->Fill(primaryUID); VClusterIDShower->Fill(primaryVID); if(primaryUID==0) UClusterProbEMShower->Fill(primaryUprobem); if(primaryVID==0) VClusterProbEMShower->Fill(primaryVprobem); AvgCluEnShower->Fill(ntpShower->ph.gev,ntpShower->ncluster); AvgCluEnShowernox->Fill(ntpShower->ph.gev,noxClu); EMFracUVDiffShower->Fill(emfracU-emfracV); Int_t mcevent = -1; Int_t hist_index = 5; if(LoadTruthForRecoTH(ntpEvent->index,mcevent)){ if(ntpTruth->iaction==1){ if(ntpTruth->inu==12 && ntpTruth->inunoosc==12) hist_index = 0; if(ntpTruth->inu==12 && ntpTruth->inunoosc==14) hist_index = 1; if(ntpTruth->inu==14 && ntpTruth->inunoosc==14) hist_index = 2; if(ntpTruth->inu==16 && ntpTruth->inunoosc==14) hist_index = 3; } else hist_index = 4; } for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; if(!LoadCluster(index)) continue; IDHistShowerSplit[hist_index]->Fill(ntpCluster->id); } if(primaryVprobem>0.2 && primaryUprobem>0.2 && primaryUID==0 && primaryVID==0) { RecoyShowerSplit[hist_index]->Fill(emfrac); } RecoEMFracShowerSplit[hist_index]->Fill(emfrac); RecoEMFracProbShowerSplit[hist_index]->Fill(emfracprob); UClusterIDShowerSplit[hist_index]->Fill(primaryUID); VClusterIDShowerSplit[hist_index]->Fill(primaryVID); if(primaryUID==0) UClusterProbEMShowerSplit[hist_index]->Fill(primaryUprobem); if(primaryVID==0) VClusterProbEMShowerSplit[hist_index]->Fill(primaryVprobem); AvgCluEnShowerSplit[hist_index]->Fill(ntpShower->ph.gev,ntpShower->ncluster); AvgCluEnnoxShowerSplit[hist_index]->Fill(ntpShower->ph.gev,noxClu); EMFracUVDiffShowerSplit[hist_index]->Fill(emfracU-emfracV); } } } cout << "DATAPOT = " << datapot << endl; save->Write(); save->Close(); return; }

void MadCluAnalysis::DataMCComp (  char *  , char *  )

Definition at line 1280 of file MadCluAnalysis.cxx. { TFile *datafile = new TFile(dataname,"READ"); TH1F *DataIDHistTrack = (TH1F*) datafile->Get("IDHistTrack"); TH1F *DataRecoyTrack = (TH1F*) datafile->Get("RecoyTrack"); TH1F *DataRecoEMFracTrack = (TH1F*) datafile->Get("RecoEMFracTrack"); TH1F *DataRecoEMFracProbTrack = (TH1F*) datafile->Get("RecoEMFracProbTrack"); TH1F *DataUClusterIDTrack = (TH1F*) datafile->Get("UClusterIDTrack"); TH1F *DataVClusterIDTrack = (TH1F*) datafile->Get("VClusterIDTrack"); TH1F *DataUClusterProbEMTrack = (TH1F*) datafile->Get("UClusterProbEMTrack"); TH1F *DataVClusterProbEMTrack = (TH1F*) datafile->Get("VClusterProbEMTrack"); TH2F *DataAvgCluEnTrack = (TH2F*) datafile->Get("AvgCluEnTrack"); TH2F *DataAvgCluEnTracknox = (TH2F*) datafile->Get("AvgCluEnTracknox"); TH1F *DataEMFracUVDiffTrack = (TH1F*) datafile->Get("EMFracUVDiffTrack"); TH1F *DataIDHistShower = (TH1F*) datafile->Get("IDHistShower"); TH1F *DataRecoyShower = (TH1F*) datafile->Get("RecoyShower"); TH1F *DataRecoEMFracShower = (TH1F*) datafile->Get("RecoEMFracShower"); TH1F *DataRecoEMFracProbShower = (TH1F*) datafile->Get("RecoEMFracProbShower"); TH1F *DataUClusterIDShower = (TH1F*) datafile->Get("UClusterIDShower"); TH1F *DataVClusterIDShower = (TH1F*) datafile->Get("VClusterIDShower"); TH1F *DataUClusterProbEMShower = (TH1F*) datafile->Get("UClusterProbEMShower"); TH1F *DataVClusterProbEMShower = (TH1F*) datafile->Get("VClusterProbEMShower"); TH2F *DataAvgCluEnShower = (TH2F*) datafile->Get("AvgCluEnShower"); TH2F *DataAvgCluEnShowernox = (TH2F*) datafile->Get("AvgCluEnShowernox"); TH1F *DataEMFracUVDiffShower = (TH1F*) datafile->Get("EMFracUVDiffShower"); TFile *mcfile = new TFile(mcname,"READ"); TH1F *IDHistTrack = (TH1F*) mcfile->Get("IDHistTrack"); TH1F *RecoyTrack = (TH1F*) mcfile->Get("RecoyTrack"); TH1F *RecoEMFracTrack = (TH1F*) mcfile->Get("RecoEMFracTrack"); TH1F *RecoEMFracProbTrack = (TH1F*) mcfile->Get("RecoEMFracProbTrack"); TH1F *UClusterIDTrack = (TH1F*) mcfile->Get("UClusterIDTrack"); TH1F *VClusterIDTrack = (TH1F*) mcfile->Get("VClusterIDTrack"); TH1F *UClusterProbEMTrack = (TH1F*) mcfile->Get("UClusterProbEMTrack"); TH1F *VClusterProbEMTrack = (TH1F*) mcfile->Get("VClusterProbEMTrack"); TH2F *AvgCluEnTrack = (TH2F*) mcfile->Get("AvgCluEnTrack"); TH2F *AvgCluEnTracknox = (TH2F*) mcfile->Get("AvgCluEnTracknox"); TH1F *EMFracUVDiffTrack = (TH1F*) mcfile->Get("EMFracUVDiffTrack"); TH1F *IDHistShower = (TH1F*) mcfile->Get("IDHistShower"); TH1F *RecoyShower = (TH1F*) mcfile->Get("RecoyShower"); TH1F *RecoEMFracShower = (TH1F*) mcfile->Get("RecoEMFracShower"); TH1F *RecoEMFracProbShower = (TH1F*) mcfile->Get("RecoEMFracProbShower"); TH1F *UClusterIDShower = (TH1F*) mcfile->Get("UClusterIDShower"); TH1F *VClusterIDShower = (TH1F*) mcfile->Get("VClusterIDShower"); TH1F *UClusterProbEMShower = (TH1F*) mcfile->Get("UClusterProbEMShower"); TH1F *VClusterProbEMShower = (TH1F*) mcfile->Get("VClusterProbEMShower"); TH2F *AvgCluEnShower = (TH2F*) mcfile->Get("AvgCluEnShower"); TH2F *AvgCluEnShowernox = (TH2F*) mcfile->Get("AvgCluEnShowernox"); TH1F *EMFracUVDiffShower = (TH1F*) mcfile->Get("EMFracUVDiffShower"); string SplitName[6] = {"beamnue","nue","numu","nutau","nc","unknown"}; TH1F *IDHistTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",IDHistTrack->GetName(),SplitName[i].c_str()); IDHistTrackSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *IDHistShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",IDHistShower->GetName(),SplitName[i].c_str()); IDHistShowerSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *RecoyTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoyTrack->GetName(),SplitName[i].c_str()); RecoyTrackSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *RecoyShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoyShower->GetName(),SplitName[i].c_str()); RecoyShowerSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *RecoEMFracTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoEMFracTrack->GetName(),SplitName[i].c_str()); RecoEMFracTrackSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *RecoEMFracShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoEMFracShower->GetName(),SplitName[i].c_str()); RecoEMFracShowerSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *RecoEMFracProbTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoEMFracProbTrack->GetName(),SplitName[i].c_str()); RecoEMFracProbTrackSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *RecoEMFracProbShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",RecoEMFracProbShower->GetName(),SplitName[i].c_str()); RecoEMFracProbShowerSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *UClusterIDTrackSplit[6] = {0}; TH1F *VClusterIDTrackSplit[6] = {0}; TH1F *UClusterProbEMTrackSplit[6] = {0}; TH1F *VClusterProbEMTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",UClusterIDTrack->GetName(),SplitName[i].c_str()); UClusterIDTrackSplit[i] = (TH1F*) mcfile->Get(name); sprintf(name,"%s_%s",VClusterIDTrack->GetName(),SplitName[i].c_str()); VClusterIDTrackSplit[i] = (TH1F*) mcfile->Get(name); sprintf(name,"%s_%s",UClusterProbEMTrack->GetName(),SplitName[i].c_str()); UClusterProbEMTrackSplit[i] = (TH1F*) mcfile->Get(name); sprintf(name,"%s_%s",VClusterProbEMTrack->GetName(),SplitName[i].c_str()); VClusterProbEMTrackSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *UClusterIDShowerSplit[6] = {0}; TH1F *VClusterIDShowerSplit[6] = {0}; TH1F *UClusterProbEMShowerSplit[6] = {0}; TH1F *VClusterProbEMShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",UClusterIDShower->GetName(),SplitName[i].c_str()); UClusterIDShowerSplit[i] = (TH1F*) mcfile->Get(name); sprintf(name,"%s_%s",VClusterIDShower->GetName(),SplitName[i].c_str()); VClusterIDShowerSplit[i] = (TH1F*) mcfile->Get(name); sprintf(name,"%s_%s",UClusterProbEMShower->GetName(),SplitName[i].c_str()); UClusterProbEMShowerSplit[i] = (TH1F*) mcfile->Get(name); sprintf(name,"%s_%s",VClusterProbEMShower->GetName(),SplitName[i].c_str()); VClusterProbEMShowerSplit[i] = (TH1F*) mcfile->Get(name); } TH2F *AvgCluEnTrackSplit[6] = {0}; TH2F *AvgCluEnnoxTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",AvgCluEnTrack->GetName(),SplitName[i].c_str()); AvgCluEnTrackSplit[i] = (TH2F*) mcfile->Get(name); sprintf(name,"%s_%s",AvgCluEnTracknox->GetName(),SplitName[i].c_str()); AvgCluEnnoxTrackSplit[i] = (TH2F*) mcfile->Get(name); } TH2F *AvgCluEnShowerSplit[6] = {0}; TH2F *AvgCluEnnoxShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",AvgCluEnShower->GetName(),SplitName[i].c_str()); AvgCluEnShowerSplit[i] = (TH2F*) mcfile->Get(name); sprintf(name,"%s_%s",AvgCluEnShowernox->GetName(),SplitName[i].c_str()); AvgCluEnnoxShowerSplit[i] = (TH2F*) mcfile->Get(name); } TH1F *EMFracUVDiffTrackSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",EMFracUVDiffTrack->GetName(),SplitName[i].c_str()); EMFracUVDiffTrackSplit[i] = (TH1F*) mcfile->Get(name); } TH1F *EMFracUVDiffShowerSplit[6] = {0}; for(int i=0;i<6;i++) { char name[256]; sprintf(name,"%s_%s",EMFracUVDiffShower->GetName(),SplitName[i].c_str()); EMFracUVDiffShowerSplit[i] = (TH1F*) mcfile->Get(name); } //Double_t mc_scale = DataUClusterIDTrack->GetEntries() + DataUClusterIDShower->GetEntries(); //mc_scale /= (UClusterIDTrack->GetEntries() + UClusterIDShower->GetEntries()); //cout << "mc_scale = " << mc_scale << endl; cout << "MD: " << 550*2.4e13*40 << " POT" << endl; cout << "DATA: " << 354657e12 << " POT" << endl; Double_t mc_scale = 3.54657e17/(550.*2.4e13*40.); IDHistTrack->Scale(mc_scale); RecoyTrack->Scale(mc_scale); RecoEMFracTrack->Scale(mc_scale); RecoEMFracProbTrack->Scale(mc_scale); UClusterIDTrack->Scale(mc_scale); VClusterIDTrack->Scale(mc_scale); UClusterProbEMTrack->Scale(mc_scale); VClusterProbEMTrack->Scale(mc_scale); EMFracUVDiffTrack->Scale(mc_scale); IDHistTrack->SetLineColor(2); RecoyTrack->SetLineColor(2); RecoEMFracTrack->SetLineColor(2); RecoEMFracProbTrack->SetLineColor(2); UClusterIDTrack->SetLineColor(2); VClusterIDTrack->SetLineColor(2); UClusterProbEMTrack->SetLineColor(2); VClusterProbEMTrack->SetLineColor(2); EMFracUVDiffTrack->SetLineColor(2); for(int i=0;i<6;i++){ IDHistTrackSplit[i]->Scale(mc_scale); RecoyTrackSplit[i]->Scale(mc_scale); RecoEMFracTrackSplit[i]->Scale(mc_scale); RecoEMFracProbTrackSplit[i]->Scale(mc_scale); UClusterIDTrackSplit[i]->Scale(mc_scale); VClusterIDTrackSplit[i]->Scale(mc_scale); UClusterProbEMTrackSplit[i]->Scale(mc_scale); VClusterProbEMTrackSplit[i]->Scale(mc_scale); EMFracUVDiffTrackSplit[i]->Scale(mc_scale); Int_t col = i+3; if(col>=5) col+=1; IDHistTrackSplit[i]->SetLineColor(col); RecoyTrackSplit[i]->SetLineColor(col); RecoEMFracTrackSplit[i]->SetLineColor(col); RecoEMFracProbTrackSplit[i]->SetLineColor(col); UClusterIDTrackSplit[i]->SetLineColor(col); VClusterIDTrackSplit[i]->SetLineColor(col); UClusterProbEMTrackSplit[i]->SetLineColor(col); VClusterProbEMTrackSplit[i]->SetLineColor(col); EMFracUVDiffTrackSplit[i]->SetLineColor(col); } IDHistShower->Scale(mc_scale); RecoyShower->Scale(mc_scale); RecoEMFracShower->Scale(mc_scale); RecoEMFracProbShower->Scale(mc_scale); UClusterIDShower->Scale(mc_scale); VClusterIDShower->Scale(mc_scale); UClusterProbEMShower->Scale(mc_scale); VClusterProbEMShower->Scale(mc_scale); EMFracUVDiffShower->Scale(mc_scale); IDHistShower->SetLineColor(2); RecoyShower->SetLineColor(2); RecoEMFracShower->SetLineColor(2); RecoEMFracProbShower->SetLineColor(2); UClusterIDShower->SetLineColor(2); VClusterIDShower->SetLineColor(2); UClusterProbEMShower->SetLineColor(2); VClusterProbEMShower->SetLineColor(2); EMFracUVDiffShower->SetLineColor(2); for(int i=0;i<6;i++){ IDHistShowerSplit[i]->Scale(mc_scale); RecoyShowerSplit[i]->Scale(mc_scale); RecoEMFracShowerSplit[i]->Scale(mc_scale); RecoEMFracProbShowerSplit[i]->Scale(mc_scale); UClusterIDShowerSplit[i]->Scale(mc_scale); VClusterIDShowerSplit[i]->Scale(mc_scale); UClusterProbEMShowerSplit[i]->Scale(mc_scale); VClusterProbEMShowerSplit[i]->Scale(mc_scale); EMFracUVDiffShowerSplit[i]->Scale(mc_scale); Int_t col = i+3; if(col>=5) col+=1; IDHistShowerSplit[i]->SetLineColor(col); RecoyShowerSplit[i]->SetLineColor(col); RecoEMFracShowerSplit[i]->SetLineColor(col); RecoEMFracProbShowerSplit[i]->SetLineColor(col); UClusterIDShowerSplit[i]->SetLineColor(col); VClusterIDShowerSplit[i]->SetLineColor(col); UClusterProbEMShowerSplit[i]->SetLineColor(col); VClusterProbEMShowerSplit[i]->SetLineColor(col); EMFracUVDiffShowerSplit[i]->SetLineColor(col); } TCanvas *can1 = new TCanvas(); can1->Divide(1,2); can1->cd(1); DataIDHistTrack->Draw(); IDHistTrack->Draw("sames"); if(IDHistTrack->GetMaximum()>DataIDHistTrack->GetMaximum()){ DataIDHistTrack->SetMaximum(IDHistTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ IDHistTrackSplit[i]->Draw("same"); } can1->cd(2); DataIDHistShower->Draw(); IDHistShower->Draw("sames"); if(IDHistShower->GetMaximum()>DataIDHistShower->GetMaximum()){ DataIDHistShower->SetMaximum(IDHistShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ IDHistShowerSplit[i]->Draw("same"); } TCanvas *can2 = new TCanvas(); can2->Divide(1,2); can2->cd(1); DataRecoyTrack->Draw(); RecoyTrack->Draw("sames"); if(RecoyTrack->GetMaximum()>DataRecoyTrack->GetMaximum()){ DataRecoyTrack->SetMaximum(RecoyTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ RecoyTrackSplit[i]->Draw("same"); } can2->cd(2); DataRecoyShower->Draw(); RecoyShower->Draw("sames"); if(RecoyShower->GetMaximum()>DataRecoyShower->GetMaximum()){ DataRecoyShower->SetMaximum(RecoyShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ RecoyShowerSplit[i]->Draw("same"); } TCanvas *can3 = new TCanvas(); can3->Divide(1,2); can3->cd(1); DataRecoEMFracTrack->Draw(); RecoEMFracTrack->Draw("sames"); if(RecoEMFracTrack->GetMaximum()>DataRecoEMFracTrack->GetMaximum()){ DataRecoEMFracTrack->SetMaximum(RecoEMFracTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ RecoEMFracTrackSplit[i]->Draw("same"); } can3->cd(2); DataRecoEMFracShower->Draw(); RecoEMFracShower->Draw("sames"); if(RecoEMFracShower->GetMaximum()>DataRecoEMFracShower->GetMaximum()){ DataRecoEMFracShower->SetMaximum(RecoEMFracShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ RecoEMFracShowerSplit[i]->Draw("same"); } TCanvas *can4 = new TCanvas(); can4->Divide(1,2); can4->cd(1); DataRecoEMFracProbTrack->Draw(); RecoEMFracProbTrack->Draw("sames"); if(RecoEMFracProbTrack->GetMaximum()>DataRecoEMFracProbTrack->GetMaximum()){ DataRecoEMFracProbTrack->SetMaximum(RecoEMFracProbTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ RecoEMFracProbTrackSplit[i]->Draw("same"); } can4->cd(2); DataRecoEMFracProbShower->Draw(); RecoEMFracProbShower->Draw("sames"); if(RecoEMFracProbShower->GetMaximum()>DataRecoEMFracProbShower->GetMaximum()){ DataRecoEMFracProbShower->SetMaximum(RecoEMFracProbShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ RecoEMFracProbShowerSplit[i]->Draw("same"); } TCanvas *can5 = new TCanvas(); can5->Divide(1,2); can5->cd(1); DataUClusterIDTrack->Draw(); UClusterIDTrack->Draw("sames"); if(UClusterIDTrack->GetMaximum()>DataUClusterIDTrack->GetMaximum()){ DataUClusterIDTrack->SetMaximum(UClusterIDTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ UClusterIDTrackSplit[i]->Draw("same"); } can5->cd(2); DataUClusterIDShower->Draw(); UClusterIDShower->Draw("sames"); if(UClusterIDShower->GetMaximum()>DataUClusterIDShower->GetMaximum()){ DataUClusterIDShower->SetMaximum(UClusterIDShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ UClusterIDShowerSplit[i]->Draw("same"); } TCanvas *can5b = new TCanvas(); can5b->Divide(1,2); can5b->cd(1); DataUClusterProbEMTrack->Draw(); UClusterProbEMTrack->Draw("sames"); if(UClusterProbEMTrack->GetMaximum()>DataUClusterProbEMTrack->GetMaximum()){ DataUClusterProbEMTrack->SetMaximum(UClusterProbEMTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ UClusterProbEMTrackSplit[i]->Draw("same"); } can5b->cd(2); DataUClusterProbEMShower->Draw(); UClusterProbEMShower->Draw("sames"); if(UClusterProbEMShower->GetMaximum()>DataUClusterProbEMShower->GetMaximum()){ DataUClusterProbEMShower->SetMaximum(UClusterProbEMShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ UClusterProbEMShowerSplit[i]->Draw("same"); } TCanvas *can6 = new TCanvas(); can6->Divide(1,2); can6->cd(1); DataVClusterIDTrack->Draw(); VClusterIDTrack->Draw("sames"); if(VClusterIDTrack->GetMaximum()>DataVClusterIDTrack->GetMaximum()){ DataVClusterIDTrack->SetMaximum(VClusterIDTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ VClusterIDTrackSplit[i]->Draw("same"); } can6->cd(2); DataVClusterIDShower->Draw(); VClusterIDShower->Draw("sames"); if(VClusterIDShower->GetMaximum()>DataVClusterIDShower->GetMaximum()){ DataVClusterIDShower->SetMaximum(VClusterIDShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ VClusterIDShowerSplit[i]->Draw("same"); } TCanvas *can6b = new TCanvas(); can6b->Divide(1,2); can6b->cd(1); DataVClusterProbEMTrack->Draw(); VClusterProbEMTrack->Draw("sames"); if(VClusterProbEMTrack->GetMaximum()>DataVClusterProbEMTrack->GetMaximum()){ DataVClusterProbEMTrack->SetMaximum(VClusterProbEMTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ VClusterProbEMTrackSplit[i]->Draw("same"); } can6b->cd(2); DataVClusterProbEMShower->Draw(); VClusterProbEMShower->Draw("sames"); if(VClusterProbEMShower->GetMaximum()>DataVClusterProbEMShower->GetMaximum()){ DataVClusterProbEMShower->SetMaximum(VClusterProbEMShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ VClusterProbEMShowerSplit[i]->Draw("same"); } TCanvas *can7 = new TCanvas(); can7->Divide(3,2); can7->cd(1); DataAvgCluEnTrack->Draw("COLZ"); can7->cd(2); AvgCluEnTrack->Draw("COLZ"); can7->cd(3); TProfile *pDT = DataAvgCluEnTrack->ProfileX("pDT"); TProfile *pMT = AvgCluEnTrack->ProfileX("pMT"); pMT->SetLineColor(2); pMT->SetMarkerColor(2); pDT->Draw(); pMT->Draw("sames"); for(int i=0;i<6;i++){ TProfile *temp = AvgCluEnTrackSplit[i]->ProfileX(); Int_t col = i+3; if(col>=5) col+=1; temp->SetLineColor(col); temp->SetMarkerColor(col); //temp->Draw("same"); } can7->cd(4); DataAvgCluEnShower->Draw("COLZ"); can7->cd(5); AvgCluEnShower->Draw("COLZ"); can7->cd(6); TProfile *pDS = DataAvgCluEnShower->ProfileX("pDS"); TProfile *pMS = AvgCluEnShower->ProfileX("pMS"); pMS->SetLineColor(2); pMS->SetMarkerColor(2); pDS->Draw(); pMS->Draw("sames"); for(int i=0;i<6;i++){ TProfile *temp = AvgCluEnShowerSplit[i]->ProfileX(); Int_t col = i+3; if(col>=5) col+=1; temp->SetLineColor(col); temp->SetMarkerColor(col); //temp->Draw("same"); } TCanvas *can8 = new TCanvas(); can8->Divide(3,2); can8->cd(1); DataAvgCluEnTracknox->Draw("COLZ"); can8->cd(2); AvgCluEnTracknox->Draw("COLZ"); can8->cd(3); TProfile *pDTx = DataAvgCluEnTracknox->ProfileX("pDTx"); TProfile *pMTx = AvgCluEnTracknox->ProfileX("pMTx"); pMTx->SetLineColor(2); pMTx->SetMarkerColor(2); pDTx->Draw(); pMTx->Draw("sames"); for(int i=0;i<6;i++){ TProfile *temp = AvgCluEnnoxTrackSplit[i]->ProfileX(); Int_t col = i+3; if(col>=5) col+=1; temp->SetLineColor(col); temp->SetMarkerColor(col); //temp->Draw("same"); } can8->cd(4); DataAvgCluEnShowernox->Draw("COLZ"); can8->cd(5); AvgCluEnShowernox->Draw("COLZ"); can8->cd(6); TProfile *pDSx = DataAvgCluEnShowernox->ProfileX("pDSx"); TProfile *pMSx = AvgCluEnShowernox->ProfileX("pMSx"); pMSx->SetLineColor(2); pMSx->SetMarkerColor(2); pDSx->Draw(); pMSx->Draw("sames"); for(int i=0;i<6;i++){ TProfile *temp = AvgCluEnnoxShowerSplit[i]->ProfileX(); Int_t col = i+3; if(col>=5) col+=1; temp->SetLineColor(col); temp->SetMarkerColor(col); //temp->Draw("same"); } TCanvas *can9 = new TCanvas(); can9->Divide(1,2); can9->cd(1); DataEMFracUVDiffTrack->Draw(); EMFracUVDiffTrack->Draw("sames"); if(EMFracUVDiffTrack->GetMaximum()>DataEMFracUVDiffTrack->GetMaximum()){ DataEMFracUVDiffTrack->SetMaximum(EMFracUVDiffTrack->GetMaximum()*1.1); } for(int i=0;i<6;i++){ EMFracUVDiffTrackSplit[i]->Draw("same"); } can9->cd(2); DataEMFracUVDiffShower->Draw(); EMFracUVDiffShower->Draw("sames"); if(EMFracUVDiffShower->GetMaximum()>DataEMFracUVDiffShower->GetMaximum()){ DataEMFracUVDiffShower->SetMaximum(EMFracUVDiffShower->GetMaximum()*1.1); } for(int i=0;i<6;i++){ EMFracUVDiffShowerSplit[i]->Draw("same"); } }

Bool_t MadCluAnalysis::FillAveNumSSPlane (  Int_t  event  )

Definition at line 2703 of file MadCluAnalysis.cxx. References AveNumSSPlane(), PAN_AveNumSSPlaneU, and PAN_AveNumSSPlaneV. Referenced by CallUserFunctions(). { Double_t *aveNumSSPlane = AveNumSSPlane(event); PAN_AveNumSSPlaneU = aveNumSSPlane[0]; PAN_AveNumSSPlaneV = aveNumSSPlane[1]; delete [] aveNumSSPlane; return true; }

Bool_t MadCluAnalysis::FillBestProbEM (  Int_t  event  )

Definition at line 2950 of file MadCluAnalysis.cxx. References NtpSRShower::clu, NtpSRStripPulseHeight::gev, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), NtpSRShower::ncluster, PAN_BestProbEM, NtpSRCluster::ph, NtpSRCluster::planeview, and NtpSRCluster::probem. Referenced by CallUserFunctions(). { PAN_BestProbEM = 0; if(!LoadEvent(event)) return false; Int_t shower = -1; if(!LoadLargestShowerFromEvent(event,shower)) return false; Double_t bestEnU = 0.; Double_t bestEnV = 0.; Double_t bestProbU = 0; Double_t bestProbV = 0; Int_t *clusters = ntpShower->clu; for(int i=0;i<ntpShower->ncluster;i++){ if(!LoadCluster(clusters[i])) continue; if(ntpCluster->planeview==2&&ntpCluster->ph.gev>bestEnU) { bestProbU = ntpCluster->probem; bestEnU = ntpCluster->ph.gev; } else if(ntpCluster->planeview==3&&ntpCluster->ph.gev>bestEnV) { bestProbV += ntpCluster->probem; bestEnV = ntpCluster->ph.gev; } } PAN_BestProbEM = (bestProbU + bestProbV)/2.; return true; }

Bool_t MadCluAnalysis::FillChargeFracRMS (  Int_t  event, Int_t  method = 0 )

Definition at line 2655 of file MadCluAnalysis.cxx. References ChargeFracRMS(), PAN_ChargeFracRMSU, and PAN_ChargeFracRMSV. Referenced by CallUserFunctions(). { Double_t *chargeFracRMS = ChargeFracRMS(event,method); PAN_ChargeFracRMSU = chargeFracRMS[0]; PAN_ChargeFracRMSV = chargeFracRMS[1]; delete [] chargeFracRMS; return true; }

Bool_t MadCluAnalysis::FillNCluster (  Int_t  event  )

Definition at line 2925 of file MadCluAnalysis.cxx. References NtpSRShower::clu, NtpSRCluster::id, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), NtpSRShower::ncluster, NtpSRShower::nUcluster, NtpSRShower::nVcluster, PAN_NClusterU, PAN_NClusterV, PAN_NPhysClusterU, PAN_NPhysClusterV, and NtpSRCluster::planeview. Referenced by CallUserFunctions(). { PAN_NClusterU = 0; PAN_NClusterV = 0; PAN_NPhysClusterU = 0; PAN_NPhysClusterV = 0; if(!LoadEvent(event)) return false; Int_t shower = -1; if(!LoadLargestShowerFromEvent(event,shower)) return false; PAN_NClusterU = ntpShower->nUcluster; PAN_NClusterV = ntpShower->nVcluster; Int_t *clusters = ntpShower->clu; for(int i=0;i<ntpShower->ncluster;i++){ if(!LoadCluster(clusters[i])) continue; if(ntpCluster->id==0||ntpCluster->id==1||ntpCluster->id==3){ if(ntpCluster->planeview==2) PAN_NPhysClusterU+=1; else if(ntpCluster->planeview==3) PAN_NPhysClusterV+=1; } } return true; }

Bool_t MadCluAnalysis::FillPHWAvgDev (  Int_t  event  )

Definition at line 2601 of file MadCluAnalysis.cxx. References PAN_PHWAvgDevU, PAN_PHWAvgDevV, and PHWAvgDev(). Referenced by CallUserFunctions(). { Double_t *phwAvgDev = PHWAvgDev(event); PAN_PHWAvgDevU = phwAvgDev[0]; PAN_PHWAvgDevV = phwAvgDev[1]; delete [] phwAvgDev; return true; }

Bool_t MadCluAnalysis::FillPHWCluID (  Int_t  event  )

Definition at line 2557 of file MadCluAnalysis.cxx. References PAN_PHWCluIDU, PAN_PHWCluIDV, and PHWCluID(). Referenced by CallUserFunctions(). { Double_t *phwCluID = PHWCluID(event); PAN_PHWCluIDU = phwCluID[0]; PAN_PHWCluIDV = phwCluID[1]; delete [] phwCluID; return true; }

Bool_t MadCluAnalysis::FillPHWProbEM (  Int_t  event  )

Definition at line 2513 of file MadCluAnalysis.cxx. References PAN_PHWProbEMU, PAN_PHWProbEMV, and PHWProbEM(). Referenced by CallUserFunctions(). { Double_t *phwProbEM = PHWProbEM(event); PAN_PHWProbEMU = phwProbEM[0]; PAN_PHWProbEMV = phwProbEM[1]; delete [] phwProbEM; return true; }

Double_t MadCluAnalysis::FOM (   )

Definition at line 1860 of file MadCluAnalysis.cxx. References abs(), NtpSRVertex::dcosz, NtpMCTruth::emfrac, MadBase::GetEntry(), NtpSRStripPulseHeight::gev, NtpMCTruth::iaction, NtpSREvent::index, NtpMCTruth::inu, NtpMCTruth::inunoosc, NtpMCTruth::iresonance, MadQuantities::IsFidVtxEvt(), MadBase::LoadEvent(), MadBase::LoadTruthForRecoTH(), NtpSREventSummary::nevent, NtpSREvent::nshower, NtpSRShower::nstrip, PassBasicCuts(), NtpSRShower::ph, PHWAvgDev(), PHWCluID(), PHWProbEM(), NtpSRShower::vtx, and NtpMCTruth::y. { double NueCCn = 0.; double NueNCn = 0.; double NueBeamn = 0.; double NumuCCn = 0.; double NumuNCn = 0.; double NutauCCn = 0.; double NutauNCn = 0.; double NueQEn = 0.; double Unknown = 0; TH1F *totHist = new TH1F("totHist","totHist",10,0,10); totHist->SetLineColor(1); TH1F *beamnueHist = new TH1F("beamnueHist","beamnueHist",10,0,10); beamnueHist->SetLineColor(8); TH1F *nueCCHist = new TH1F("nueCCHist","nueCCHist",10,0,10); nueCCHist->SetLineColor(3); TH1F *numuCCHist = new TH1F("numuCCHist","numuCCHist",10,0,10); numuCCHist->SetLineColor(4); TH1F *nutauCCHist = new TH1F("nutauCCHist","nutauCCHist",10,0,10); nutauCCHist->SetLineColor(6); TH1F *NCHist = new TH1F("NCHist","NCHist",10,0,10); NCHist->SetLineColor(2); TH2F *totHist_nstp = new TH2F("totHist_nstp","totHist_nstp",10,0,10,300,0,300); totHist_nstp->SetLineColor(1); TH2F *beamnueHist_nstp = new TH2F("beamnueHist_nstp","beamnueHist_nstp",10,0,10,300,0,300); beamnueHist_nstp->SetLineColor(8); TH2F *nueCCHist_nstp = new TH2F("nueCCHist_nstp","nueCCHist_nstp",10,0,10,300,0,300); nueCCHist_nstp->SetLineColor(3); TH2F *numuCCHist_nstp = new TH2F("numuCCHist_nstp","numuCCHist_nstp",10,0,10,300,0,300); numuCCHist_nstp->SetLineColor(4); TH2F *nutauCCHist_nstp = new TH2F("nutauCCHist_nstp","nutauCCHist_nstp",10,0,10,300,0,300); nutauCCHist_nstp->SetLineColor(6); TH2F *NCHist_nstp = new TH2F("NCHist_nstp","NCHist_nstp",10,0,10,300,0,300); NCHist_nstp->SetLineColor(2); TH1F *totHist_dcz = new TH1F("totHist_dcz","totHist_dcz",100,0,1); totHist_dcz->SetLineColor(1); TH1F *beamnueHist_dcz = new TH1F("beamnueHist_dcz","beamnueHist_dcz",100,0,1); beamnueHist_dcz->SetLineColor(8); TH1F *nueCCHist_dcz = new TH1F("nueCCHist_dcz","nueCCHist_dcz",100,0,1); nueCCHist_dcz->SetLineColor(3); TH1F *numuCCHist_dcz = new TH1F("numuCCHist_dcz","numuCCHist_dcz",100,0,1); numuCCHist_dcz->SetLineColor(4); TH1F *nutauCCHist_dcz = new TH1F("nutauCCHist_dcz","nutauCCHist_dcz",100,0,1); nutauCCHist_dcz->SetLineColor(6); TH1F *NCHist_dcz = new TH1F("NCHist_dcz","NCHist_dcz",100,0,1); NCHist_dcz->SetLineColor(2); TH2F *totHist_phw = new TH2F("totHist_phw","totHist_phw",30,-0.5,5.5,20,0,1); totHist_phw->SetLineColor(1); TH2F *beamnueHist_phw = new TH2F("beamnueHist_phw","beamnueHist_phw",30,-0.5,5.5,20,0,1); beamnueHist_phw->SetLineColor(8); TH2F *nueCCHist_phw = new TH2F("nueCCHist_phw","nueCCHist_phw",30,-0.5,5.5,20,0,1); nueCCHist_phw->SetLineColor(3); TH2F *numuCCHist_phw = new TH2F("numuCCHist_phw","numuCCHist_phw",30,-0.5,5.5,20,0,1); numuCCHist_phw->SetLineColor(4); TH2F *nutauCCHist_phw = new TH2F("nutauCCHist_phw","nutauCCHist_phw",30,-0.5,5.5,20,0,1); nutauCCHist_phw->SetLineColor(6); TH2F *NCHist_phw = new TH2F("NCHist_phw","NCHist_phw",30,-0.5,5.5,20,0,1); NCHist_phw->SetLineColor(2); TH2F *highyCCHist_phw = new TH2F("highyHist_phw","highyHist_phw", 30,-0.5,5.5,20,0,1); TH2F *midyCCHist_phw = new TH2F("midyHist_phw","midyHist_phw", 30,-0.5,5.5,20,0,1); TH2F *lowyCCHist_phw = new TH2F("lowyHist_phw","lowyHist_phw", 30,-0.5,5.5,20,0,1); TH2F *qenueCCHist_phw = new TH2F("qenueCCHist_phw","qenueCCHist_phw", 30,-0.5,5.5,20,0,1); TH1F *totHist_avgdev = new TH1F("totHist_avgdev","totHist_avgdev",100,0,0.2); totHist_avgdev->SetLineColor(1); TH1F *beamnueHist_avgdev = new TH1F("beamnueHist_avgdev","beamnueHist_avgdev",100,0,0.2); beamnueHist_avgdev->SetLineColor(8); TH1F *nueCCHist_avgdev = new TH1F("nueCCHist_avgdev","nueCCHist_avgdev",100,0,0.2); nueCCHist_avgdev->SetLineColor(3); TH1F *numuCCHist_avgdev = new TH1F("numuCCHist_avgdev","numuCCHist_avgdev",100,0,0.2); numuCCHist_avgdev->SetLineColor(4); TH1F *nutauCCHist_avgdev = new TH1F("nutauCCHist_avgdev","nutauCCHist_avgdev",100,0,0.2); nutauCCHist_avgdev->SetLineColor(6); TH1F *NCHist_avgdev = new TH1F("NCHist_avgdev","NCHist_avgdev",100,0,0.2); NCHist_avgdev->SetLineColor(2); TH1F *highyCCHist_avgdev = new TH1F("highyHist_avgdev","highyHist_avgdev", 100,0,0.2); highyCCHist_avgdev->SetLineColor(5); TH1F *midyCCHist_avgdev = new TH1F("midyHist_avgdev","midyHist_avgdev", 100,0,0.2); midyCCHist_avgdev->SetLineColor(9); TH1F *lowyCCHist_avgdev = new TH1F("lowyHist_avgdev","lowyHist_avgdev", 100,0,0.2); lowyCCHist_avgdev->SetLineColor(3); TH1F *qenueCCHist_avgdev = new TH1F("qenueCCHist_avgdev", "qenueCCHist_avgdev", 100,0,0.2); qenueCCHist_avgdev->SetLineColor(4); TH2F *totHist_ratio = new TH2F("totHist_ratio","totHist_ratio",20,0,5,20,0,1); totHist_ratio->SetLineColor(1); TH2F *beamnueHist_ratio = new TH2F("beamnueHist_ratio","beamnueHist_ratio",20,0,5,20,0,1); beamnueHist_ratio->SetLineColor(8); TH2F *nueCCHist_ratio = new TH2F("nueCCHist_ratio","nueCCHist_ratio",20,0,5,20,0,1); nueCCHist_ratio->SetLineColor(3); TH2F *numuCCHist_ratio = new TH2F("numuCCHist_ratio","numuCCHist_ratio",20,0,5,20,0,1); numuCCHist_ratio->SetLineColor(4); TH2F *nutauCCHist_ratio = new TH2F("nutauCCHist_ratio","nutauCCHist_ratio",20,0,5,20,0,1); nutauCCHist_ratio->SetLineColor(6); TH2F *NCHist_ratio = new TH2F("NCHist_ratio","NCHist_ratio",20,0,5,20,0,1); NCHist_ratio->SetLineColor(2); TH2F *highyCCHist_ratio = new TH2F("highyHist_ratio","highyHist_ratio", 20,0,5,20,0,1); TH2F *midyCCHist_ratio = new TH2F("midyHist_ratio","midyHist_ratio", 20,0,5,20,0,1); TH2F *lowyCCHist_ratio = new TH2F("lowyHist_ratio","lowyHist_ratio", 20,0,5,20,0,1); TH2F *qenueCCHist_ratio = new TH2F("qenueCCHist_ratio","qenueCCHist_ratio", 20,0,5,20,0,1); Double_t nNuMuFiles = 10; Double_t nNuEFiles = 20; Double_t nNuTauFiles = 10; Double_t POT = 9.25; for(int i=0;i<Nentries;i++){ if(i%1000==0) std::cout << float(i)*100./float(Nentries) << "% done" << std::endl; if(!this->GetEntry(i)) continue; if(eventSummary->nevent==0) continue; //fill tree once for each reconstructed event for(int j=0;j<eventSummary->nevent;j++){ if(!LoadEvent(j)) continue; //no event found if(ntpEvent->nshower==0) continue; //no shower found //if(ntpEvent->ntrack!=0) continue; //no shower found if(!IsFidVtxEvt(ntpEvent->index)) continue; if(PassBasicCuts(ntpEvent->index)){ double prob = 1.; Int_t mcevent = -1; if(LoadTruthForRecoTH(ntpEvent->index,mcevent)) { //prob = Oscillate(ntpTruth->inu,ntpTruth->inunoosc, // ntpTruth->p4neu[3],735.,0.0025, // TMath::ASin(1.0)/2.,0.01); //prob = Oscillate(ntpTruth->inu,ntpTruth->inunoosc, // ntpTruth->p4neu[3],735.,0.002175, // TMath::ASin(TMath::Sqrt(0.925))/2., // 0.01); } //double emfracU = 0; //double emfracV = 0; //double emfracprob = RecoEMFrac(ntpShower->index,0,emfracU,emfracV); //if(emfracprob>0.75){ if(ntpTruth->emfrac>0.6){ //if(PassAnalysisCuts(ntpEvent->index)){ if( ntpShower->ph.gev>8. || ntpShower->ph.gev<2.) continue; if(mcevent>=0) { Double_t *phwcluid = this->PHWCluID(ntpEvent->index); Double_t *phwprobem = this->PHWProbEM(ntpEvent->index); Double_t rat_cluid = phwcluid[0]+phwcluid[1]; Double_t rat_probem = 0; if(phwprobem[0]>0 && phwprobem[1]>0) { if(phwprobem[0]>phwprobem[1]) rat_probem = phwprobem[1]/phwprobem[0]; else rat_probem = phwprobem[0]/phwprobem[1]; } Double_t *avgdev = this->PHWAvgDev(ntpEvent->index); double w1 = 0; if(ntpTruth->iaction==1){ if(abs(ntpTruth->inunoosc)==12 && abs(ntpTruth->inu)==12){ w1 = prob*(POT/((nNuEFiles+nNuMuFiles)*6.5)); NueBeamn+=w1; beamnueHist->Fill(ntpShower->ph.gev,w1); beamnueHist_nstp->Fill(ntpShower->ph.gev,ntpShower->nstrip,w1); beamnueHist_dcz->Fill(ntpShower->vtx.dcosz,w1); beamnueHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) beamnueHist_avgdev->Fill(avgdev[0],w1); beamnueHist_ratio->Fill(rat_cluid,rat_probem,w1); } else if(abs(ntpTruth->inu)==12) { w1 = prob*(POT/(nNuEFiles*6.5)); NueCCn+=w1; nueCCHist->Fill(ntpShower->ph.gev,w1); nueCCHist_nstp->Fill(ntpShower->ph.gev,ntpShower->nstrip,w1); nueCCHist_dcz->Fill(ntpShower->vtx.dcosz,w1); nueCCHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) nueCCHist_avgdev->Fill(avgdev[0],w1); nueCCHist_ratio->Fill(rat_cluid,rat_probem,w1); if(ntpTruth->iresonance==1001) { NueQEn+=w1; qenueCCHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) qenueCCHist_avgdev->Fill(avgdev[0],w1); qenueCCHist_ratio->Fill(rat_cluid,rat_probem,w1); } } else if(abs(ntpTruth->inu)==14) { w1 = prob*(POT/(nNuMuFiles*6.5)); NumuCCn+=w1; numuCCHist->Fill(ntpShower->ph.gev,w1); numuCCHist_nstp->Fill(ntpShower->ph.gev,ntpShower->nstrip,w1); numuCCHist_dcz->Fill(ntpShower->vtx.dcosz,w1); numuCCHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) numuCCHist_avgdev->Fill(avgdev[0],w1); numuCCHist_ratio->Fill(rat_cluid,rat_probem,w1); } else if(abs(ntpTruth->inu)==16) { w1 = prob*(POT/(nNuTauFiles*6.5)); NutauCCn+=w1; nutauCCHist->Fill(ntpShower->ph.gev,w1); nutauCCHist_nstp->Fill(ntpShower->ph.gev,ntpShower->nstrip,w1); nutauCCHist_dcz->Fill(ntpShower->vtx.dcosz,w1); nutauCCHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) nutauCCHist_avgdev->Fill(avgdev[0],w1); nutauCCHist_ratio->Fill(rat_cluid,rat_probem,w1); } if(ntpTruth->y<0.3) { lowyCCHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) lowyCCHist_avgdev->Fill(avgdev[0],w1); lowyCCHist_ratio->Fill(rat_cluid,rat_probem,w1); } else if(ntpTruth->y<0.6) { midyCCHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) midyCCHist_avgdev->Fill(avgdev[0],w1); midyCCHist_ratio->Fill(rat_cluid,rat_probem,w1); } else { highyCCHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) highyCCHist_avgdev->Fill(avgdev[0],w1); highyCCHist_ratio->Fill(rat_cluid,rat_probem,w1); } } else if(ntpTruth->iaction==0){ if(abs(ntpTruth->inu)==12) { w1 = prob*(POT/((nNuEFiles)*6.5)); NueNCn+=w1; } else if(abs(ntpTruth->inu)==14) { w1 = prob*(POT/((nNuMuFiles)*6.5)); NumuNCn+=w1; } else if(abs(ntpTruth->inu)==16) { w1 = prob*(POT/((nNuTauFiles)*6.5)); NutauNCn+=w1; } NCHist->Fill(ntpShower->ph.gev,w1); NCHist_nstp->Fill(ntpShower->ph.gev,ntpShower->nstrip,w1); NCHist_dcz->Fill(ntpShower->vtx.dcosz,w1); NCHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) NCHist_avgdev->Fill(avgdev[0],w1); NCHist_ratio->Fill(rat_cluid,rat_probem,w1); } totHist->Fill(ntpShower->ph.gev,w1); totHist_nstp->Fill(ntpShower->ph.gev,ntpShower->nstrip,w1); totHist_dcz->Fill(ntpShower->vtx.dcosz,w1); totHist_phw->Fill(phwcluid[0],phwprobem[0],w1); if(phwprobem[0]<0.05) totHist_avgdev->Fill(avgdev[0],w1); totHist_ratio->Fill(rat_cluid,rat_probem,w1); delete [] phwprobem; delete [] phwcluid; delete [] avgdev; } else Unknown+=prob*(POT/((nNuEFiles+nNuMuFiles+nNuTauFiles)*6.5)); } } } } Double_t fom = NueCCn/sqrt(NueBeamn+NumuCCn+NutauCCn+NueNCn+NumuNCn+NutauNCn+Unknown); cout<<"NueBeamn "<<NueBeamn<<" NueCCn "<<NueCCn<<" NumuCCn "<<NumuCCn<<" NutauCCn "<<NutauCCn <<" NueNCn "<<NueNCn<<" NumuNCn "<<NumuNCn<<" NutauNCn "<<NutauNCn<<" Unknown "<<Unknown << " fom " << fom<<" NueQEn "<<NueQEn<<endl; TCanvas *can = new TCanvas(); can->Divide(2,2); can->cd(1); totHist->Draw(); NCHist->Draw("sames"); numuCCHist->Draw("sames"); nutauCCHist->Draw("sames"); nueCCHist->Draw("sames"); beamnueHist->Draw("sames"); can->cd(2); totHist_dcz->Draw(); NCHist_dcz->Draw("sames"); numuCCHist_dcz->Draw("sames"); nutauCCHist_dcz->Draw("sames"); nueCCHist_dcz->Draw("sames"); beamnueHist_dcz->Draw("sames"); can->cd(3); NCHist_nstp->Draw("COLZ"); can->cd(4); nueCCHist_nstp->Draw("COLZ"); TCanvas *can2 = new TCanvas(); can2->Divide(2,3); can2->cd(1); totHist_phw->Draw("COLZ"); can2->cd(2); numuCCHist_phw->Draw("COLZ"); can2->cd(3); nutauCCHist_phw->Draw("COLZ"); can2->cd(4); nueCCHist_phw->Draw("COLZ"); can2->cd(5); beamnueHist_phw->Draw("COLZ"); can2->cd(6); NCHist_phw->Draw("COLZ"); TCanvas *can3 = new TCanvas(); can3->Divide(2,2); can3->cd(1); lowyCCHist_phw->Draw("LEGO2"); can3->cd(2); midyCCHist_phw->Draw("LEGO2"); can3->cd(3); highyCCHist_phw->Draw("LEGO2"); can3->cd(4); qenueCCHist_phw->Draw("LEGO2"); TCanvas *can4 = new TCanvas(); can4->Divide(1,2); can4->cd(1); totHist_avgdev->Draw(); numuCCHist_avgdev->Draw("sames"); nutauCCHist_avgdev->Draw("sames"); nueCCHist_avgdev->Draw("sames"); beamnueHist_avgdev->Draw("sames"); NCHist_avgdev->Draw("sames"); can4->cd(2); lowyCCHist_avgdev->Draw(); qenueCCHist_avgdev->Draw("sames"); midyCCHist_avgdev->Draw("sames"); highyCCHist_avgdev->Draw("sames"); NCHist_avgdev->Draw("sames"); TCanvas *can5 = new TCanvas(); can5->Divide(2,3); can5->cd(1); totHist_ratio->Draw("COLZ"); can5->cd(2); numuCCHist_ratio->Draw("COLZ"); can5->cd(3); nutauCCHist_ratio->Draw("COLZ"); can5->cd(4); nueCCHist_ratio->Draw("COLZ"); can5->cd(5); beamnueHist_ratio->Draw("COLZ"); can5->cd(6); NCHist_ratio->Draw("COLZ"); TCanvas *can6 = new TCanvas(); can6->Divide(2,2); can6->cd(1); lowyCCHist_ratio->Draw("LEGO2"); can6->cd(2); midyCCHist_ratio->Draw("LEGO2"); can6->cd(3); highyCCHist_ratio->Draw("LEGO2"); can6->cd(4); qenueCCHist_ratio->Draw("LEGO2"); return fom; }

Float_t MadCluAnalysis::GetWeight (  Int_t  event = 0  )  [protected, virtual]

Reimplemented from MadAnalysis. Definition at line 2471 of file MadCluAnalysis.cxx. { return 1; }

Bool_t MadCluAnalysis::IsNueFidVtxEvt (  Int_t  event  )

Definition at line 2908 of file MadCluAnalysis.cxx. References VldContext::GetDetector(), RecHeader::GetVldContext(), MadBase::LoadEvent(), NtpSREvent::vtx, NtpSRVertex::x, NtpSRVertex::y, and NtpSRVertex::z. Referenced by CallUserFunctions(). { if(!LoadEvent(ievt)) return false; if(ntpHeader->GetVldContext().GetDetector()==Detector::kFar){ if(ntpEvent->vtx.z<1 || ntpEvent->vtx.z>29 || //ends (ntpEvent->vtx.z<17&&ntpEvent->vtx.z>14) || //between SMs sqrt((ntpEvent->vtx.x*ntpEvent->vtx.x) //radial cut +(ntpEvent->vtx.y*ntpEvent->vtx.y))>3.87) return false; } else if(ntpHeader->GetVldContext().GetDetector()==Detector::kNear){ if(ntpEvent->vtx.z<1 || ntpEvent->vtx.z>5 || sqrt(((ntpEvent->vtx.x-1.4885)*(ntpEvent->vtx.x-1.4885)) + ((ntpEvent->vtx.y-0.1397)*(ntpEvent->vtx.y-0.1397)))>1) return false; } return true; }

Bool_t MadCluAnalysis::IsTrackInSlice (  Int_t  slice  )

Definition at line 1845 of file MadCluAnalysis.cxx. References NtpSRTrack::index, MadBase::LoadSlice(), MadBase::LoadTrack(), NtpSREventSummary::ntrack, and NtpSRTrack::slc. Referenced by PassDataCleaningCuts(). { Int_t cur_track = -1; if(ntpTrack) cur_track = ntpTrack->index; if(!LoadSlice(slice)) return false; for(int i=0;i<eventSummary->ntrack;i++){ if(!LoadTrack(i)) continue; if(ntpTrack->slc==slice) { LoadTrack(cur_track); return true; } } return false; }

Bool_t MadCluAnalysis::PassAnalysisCuts (  Int_t  event = 0  )  [protected, virtual]

Reimplemented from MadAnalysis. Definition at line 2283 of file MadCluAnalysis.cxx. References PassBasicCuts(), and PID(). { if(!PassBasicCuts(event)) return false; if(PID(event,1)>0. || PID(event,2)>0.) return true; return false; }

Bool_t MadCluAnalysis::PassBasicCuts (  Int_t  event = 0  )  [protected, virtual]

Reimplemented from MadAnalysis. Definition at line 2234 of file MadCluAnalysis.cxx. References NtpSRPlane::beg, NtpSRPlane::end, NtpSRMomentum::eqp, NtpSRTrack::fit, MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), MadBase::LoadLargestTrackFromEvent(), NtpSRTrack::momentum, NtpSRPlane::n, NtpSREvent::nshower, NtpSREvent::ntrack, NtpSRFitTrack::pass, NtpSRShower::plane, NtpSRTrack::plane, and NtpSRMomentum::qp. Referenced by FOM(), and PassAnalysisCuts(). { if(!LoadEvent(event)) return false; if(ntpEvent->nshower!=1) return false; Int_t shower = -1; if(!LoadLargestShowerFromEvent(event,shower)) return false; if(ntpEvent->ntrack==0) return true; Int_t ntrkplane = 0; Int_t trkpass = 0; Double_t trkerror = 0; Int_t track = -1; if(LoadLargestTrackFromEvent(event,track)) { ntrkplane = ntpTrack->plane.end-ntpTrack->plane.beg+1; if (ntrkplane<ntpTrack->plane.n) ntrkplane = ntpTrack->plane.n; if(ntpTrack->momentum.qp!=0) trkerror = ntpTrack->momentum.eqp/ntpTrack->momentum.qp; else trkerror = 1; trkpass = ntpTrack->fit.pass; } double trkshw_overlap = 0.; if (ntrkplane>0){ if(ntpShower->plane.beg<=ntpTrack->plane.beg && ntpShower->plane.end>=ntpTrack->plane.end) trkshw_overlap = 1.; else if(ntpShower->plane.beg>ntpTrack->plane.beg && ntpShower->plane.end<ntpTrack->plane.end) trkshw_overlap = double(ntpShower->plane.end - ntpShower->plane.beg+1)/double(ntrkplane); else if(ntpShower->plane.beg<=ntpTrack->plane.beg && ntpShower->plane.end>=ntpTrack->plane.beg) trkshw_overlap = double(ntpShower->plane.end - ntpTrack->plane.beg+1)/double(ntrkplane); else if(ntpShower->plane.beg<=ntpTrack->plane.end && ntpShower->plane.end>=ntpTrack->plane.end) trkshw_overlap = double(ntpTrack->plane.end - ntpShower->plane.beg+1)/double(ntrkplane); } if(trkshw_overlap>0.6 && (trkerror>0.4 || trkpass==0 || ntrkplane<20)) return true; return false; }

Bool_t MadCluAnalysis::PassDataCleaningCuts (  Int_t  shower  )

Definition at line 1820 of file MadCluAnalysis.cxx. References NtpSRShower::clu, NtpSRStripPulseHeight::gev, NtpSRCluster::id, IsTrackInSlice(), MadBase::LoadCluster(), MadBase::LoadShower(), NtpSRShower::ncluster, NtpSRShower::ph, NtpSRCluster::planeview, and NtpSRShower::slc. Referenced by DataDistributions(). { if(!LoadShower(shower)) return false; if(ntpShower->ph.gev<0.2) return false; Int_t nphyclusU = 0; Int_t nphyclusV = 0; for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; LoadCluster(index); if(ntpCluster->id==0 || ntpCluster->id==1 || ntpCluster->id==3) { if(ntpCluster->planeview==2) nphyclusU += 1; if(ntpCluster->planeview==3) nphyclusV += 1; } } if(nphyclusU==0&&nphyclusV==0) return 0; if(nphyclusU==0||nphyclusV==0) { if(IsTrackInSlice(ntpShower->slc) && ntpShower->ph.gev<0.5) return false; } return true; }

Bool_t MadCluAnalysis::PassTruthSignal (  Int_t  mcevent = 0  )  [protected, virtual]

Reimplemented from MadAnalysis. Definition at line 2289 of file MadCluAnalysis.cxx. References NtpMCTruth::iaction, NtpMCTruth::inu, and MadBase::LoadTruth(). { if(!LoadTruth(mcevent)) return false; if(ntpTruth->inu==12 && ntpTruth->iaction==1) return true; return false; }

Double_t * MadCluAnalysis::PHWAvgDev (  Int_t  event  )

Definition at line 2566 of file MadCluAnalysis.cxx. References NtpSRCluster::avgdev, NtpSRShower::clu, NtpSRStripPulseHeight::gev, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), NtpSRShower::ncluster, NtpSRCluster::ph, and NtpSRCluster::planeview. Referenced by FillPHWAvgDev(), and FOM(). { Double_t *phwAvgDev = new Double_t[2]; phwAvgDev[0] = 0; phwAvgDev[1] = 0; if(!LoadEvent(event)) { phwAvgDev[0] = -10.; phwAvgDev[1] = -10.; return phwAvgDev; } Int_t shower = -1; if(!LoadLargestShowerFromEvent(event,shower)) { phwAvgDev[0] = -10.; phwAvgDev[1] = -10.; return phwAvgDev; } Double_t totPH[2] = {0}; Int_t *clusters = ntpShower->clu; for(int k=0; k<ntpShower->ncluster; k++){ Int_t index = clusters[k]; if(!LoadCluster(index)) continue; if(ntpCluster->planeview==2){ phwAvgDev[0] += ntpCluster->ph.gev*ntpCluster->avgdev; totPH[0] += ntpCluster->ph.gev; } else if(ntpCluster->planeview==3){ phwAvgDev[1] += ntpCluster->ph.gev*ntpCluster->avgdev; totPH[1] += ntpCluster->ph.gev; } } phwAvgDev[0]/=totPH[0]; phwAvgDev[1]/=totPH[1]; return phwAvgDev; }

Double_t * MadCluAnalysis::PHWCluID (  Int_t  event  )

Definition at line 2522 of file MadCluAnalysis.cxx. References NtpSRShower::clu, NtpSRStripPulseHeight::gev, NtpSRCluster::id, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), NtpSRShower::ncluster, NtpSRCluster::ph, and NtpSRCluster::planeview. Referenced by FillPHWCluID(), and FOM(). { Double_t *phwCluID = new Double_t[2]; phwCluID[0] = 0; phwCluID[1] = 0; if(!LoadEvent(event)) { phwCluID[0] = -10.; phwCluID[1] = -10.; return phwCluID; } Int_t shower = -1; if(!LoadLargestShowerFromEvent(event,shower)) { phwCluID[0] = -10.; phwCluID[1] = -10.; return phwCluID; } Double_t totPH[2] = {0}; Int_t *clusters = ntpShower->clu; for(int k=0; k<ntpShower->ncluster; k++){ Int_t index = clusters[k]; if(!LoadCluster(index)) continue; if(ntpCluster->planeview==2){ phwCluID[0] += ntpCluster->ph.gev*ntpCluster->id; totPH[0] += ntpCluster->ph.gev; } else if(ntpCluster->planeview==3){ phwCluID[1] += ntpCluster->ph.gev*ntpCluster->id; totPH[1] += ntpCluster->ph.gev; } } phwCluID[0]/=totPH[0]; phwCluID[1]/=totPH[1]; return phwCluID; }

Double_t * MadCluAnalysis::PHWProbEM (  Int_t  event  )

Definition at line 2477 of file MadCluAnalysis.cxx. References NtpSRShower::clu, NtpSRStripPulseHeight::gev, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), NtpSRShower::ncluster, NtpSRCluster::ph, NtpSRCluster::planeview, and NtpSRCluster::probem. Referenced by FillPHWProbEM(), and FOM(). { Double_t *phwProbEM = new Double_t[2]; phwProbEM[0] = 0; phwProbEM[1] = 0; if(!LoadEvent(event)) { phwProbEM[0] = -10.; phwProbEM[1] = -10.; return phwProbEM; } Int_t shower = -1; if(!LoadLargestShowerFromEvent(event,shower)) { phwProbEM[0] = -10.; phwProbEM[1] = -10.; return phwProbEM; } Double_t totPH[2] = {0}; Int_t *clusters = ntpShower->clu; for(int k=0; k<ntpShower->ncluster; k++){ Int_t index = clusters[k]; if(!LoadCluster(index)) continue; if(ntpCluster->planeview==2){ if(ntpCluster->ph.gev>0.7) phwProbEM[0] += ntpCluster->ph.gev*ntpCluster->probem; totPH[0] += ntpCluster->ph.gev; } else if(ntpCluster->planeview==3){ if(ntpCluster->ph.gev>0.7) phwProbEM[1] += ntpCluster->ph.gev*ntpCluster->probem; totPH[1] += ntpCluster->ph.gev; } } phwProbEM[0]/=totPH[0]; phwProbEM[1]/=totPH[1]; return phwProbEM; }

Float_t MadCluAnalysis::PID (  Int_t  event = 0, Int_t  method = 0 )  [protected, virtual]

Reimplemented from MadAnalysis. Definition at line 2296 of file MadCluAnalysis.cxx. References NtpSRCluster::avgdev, NtpSRCluster::begplane, NtpSRPlane::begu, NtpSRPlane::begv, NtpSRShower::clu, flikelihoodDists, fneuralNet, NtpSRStripPulseHeight::gev, NtpSRCluster::id, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), MadBase::LoadLargestTrackFromEvent(), NtpSRShower::ncluster, PAN_ChargeFracRMSU, PAN_ChargeFracRMSV, PAN_PHWAvgDevU, PAN_PHWAvgDevV, PAN_PHWCluIDU, PAN_PHWProbEMU, PAN_PHWProbEMV, NtpSRCluster::ph, NtpSRShower::plane, NtpSRCluster::planeview, and NtpSRCluster::probem. Referenced by CallUserFunctions(), and PassAnalysisCuts(). { if(!LoadEvent(event)) return -10.; Int_t shower = -1; Int_t track = -1; LoadLargestTrackFromEvent(event,track); if(!LoadLargestShowerFromEvent(event,shower)) return -10.; if(method==0){ int goodU = 0; int goodV = 0; double maxU = 0; double maxV = 0; double avgdevU = 1.; double avgdevV = 1.; Int_t nemclu = 0; Int_t nhadclu = 0; Int_t nphysclu = 0; Int_t *clusters = ntpShower->clu; for(int k=0; k<ntpShower->ncluster; k++){ Int_t index = clusters[k]; if(!LoadCluster(index)) continue; if(ntpCluster->id==0) nemclu++; else if(ntpCluster->id==1) nhadclu++; else if(ntpCluster->id!=2 && ntpCluster->ph.gev>0.5) nphysclu++; if(ntpCluster->planeview==2){ if(ntpCluster->ph.gev>maxU){ if(ntpCluster->id==0 && TMath::Abs(ntpCluster->begplane-ntpShower->plane.begu)<=4){ goodU = 1; if(ntpCluster->probem>0.2) goodU=3; avgdevU = ntpCluster->avgdev; maxU = ntpCluster->ph.gev; } else if(ntpCluster->id==1 && ntpCluster->probem>0.5 && TMath::Abs(ntpCluster->begplane-ntpShower->plane.begu)<=4){ goodU = 2; avgdevU = ntpCluster->avgdev; maxU = ntpCluster->ph.gev; } else { goodU = 0; avgdevU = ntpCluster->avgdev; maxU = ntpCluster->ph.gev; } } } else if(ntpCluster->planeview==3){ if(ntpCluster->ph.gev>maxV){ if(ntpCluster->id==0 && TMath::Abs(ntpCluster->begplane-ntpShower->plane.begv)<=4){ goodV = 1; if(ntpCluster->probem>0.2) goodV=3; avgdevV = ntpCluster->avgdev; maxV = ntpCluster->ph.gev; } else if(ntpCluster->id==1 && ntpCluster->probem>0.5 && TMath::Abs(ntpCluster->begplane-ntpShower->plane.begv)<=4){ goodV = 2; avgdevV = ntpCluster->avgdev; maxV = ntpCluster->ph.gev; } else { goodV = 0; avgdevV = ntpCluster->avgdev; maxV = ntpCluster->ph.gev; } } } } if( ( (goodU==1 && avgdevU<0.05 && goodV==1 && avgdevV<0.05) || goodU==3 || goodV==3 ) && nemclu < 3){ return 1; } return 0.; } else if(method==1 || method==2){ Int_t index = 0; Float_t id = PAN_PHWCluIDU; Float_t em = PAN_PHWProbEMU; Float_t cf = PAN_ChargeFracRMSU; if(method==2) { index = 6; id = PAN_PHWCluIDV; em = PAN_PHWProbEMV; cf = PAN_ChargeFracRMSV; } Float_t ProbNC = 1.; Float_t ProbQE = 1.; Float_t PidParQENC = 0; //NC: int bin1 = flikelihoodDists[index+1]->FindBin(id); ProbNC *= flikelihoodDists[index+1]->GetBinContent(bin1); bin1 = flikelihoodDists[index+3]->FindBin(em); ProbNC *= flikelihoodDists[index+3]->GetBinContent(bin1); bin1 = flikelihoodDists[index+5]->FindBin(cf); ProbNC *= flikelihoodDists[index+5]->GetBinContent(bin1); //QE: bin1 = flikelihoodDists[index]->FindBin(id); ProbQE *= flikelihoodDists[index]->GetBinContent(bin1); bin1 = flikelihoodDists[index+2]->FindBin(em); ProbQE *= flikelihoodDists[index+2]->GetBinContent(bin1); bin1 = flikelihoodDists[index+4]->FindBin(cf); ProbQE *= flikelihoodDists[index+4]->GetBinContent(bin1); if(ProbQE!=0&&ProbNC!=0) PidParQENC = TMath::Sqrt(-TMath::Log(ProbNC)) - TMath::Sqrt(-TMath::Log(ProbQE)); else if(ProbQE==0 && ProbNC==0) PidParQENC=-10; else if(ProbQE==0) PidParQENC = TMath::Sqrt(-TMath::Log(ProbNC)) - 10; else if(ProbNC==0) PidParQENC = 10 - TMath::Sqrt( -TMath::Log(ProbQE)); return PidParQENC; } else if(method==3){ Double_t params[4]; params[0] = PAN_ChargeFracRMSU*2.; params[1] = PAN_PHWProbEMU*PAN_PHWProbEMU; params[2] = PAN_PHWCluIDU; params[3] = PAN_PHWAvgDevU*2.; return fneuralNet->Evaluate(0,params); } else if(method==4){ Double_t params[4]; params[0] = PAN_ChargeFracRMSV*2.; params[1] = PAN_PHWProbEMV*PAN_PHWProbEMV; params[2] = PAN_PHWCluIDV; params[3] = PAN_PHWAvgDevV*2.; return fneuralNet->Evaluate(0,params); } else if(method==5){ Double_t params[4]; params[0] = PAN_ChargeFracRMSU + PAN_ChargeFracRMSV; params[1] = (PAN_PHWProbEMU*PAN_PHWProbEMU + PAN_PHWProbEMV*PAN_PHWProbEMV)/2.; params[2] = (PAN_PHWCluIDU + PAN_PHWCluIDV)/2.; params[3] = PAN_PHWAvgDevU + PAN_PHWAvgDevV; return fneuralNet->Evaluate(0,params); } return 0; }

void MadCluAnalysis::Plot (  char *  fileName, Int_t  startEnt = 0 )

Definition at line 169 of file MadCluAnalysis.cxx. References abs(), NtpSRShower::clu, NtpMCTruth::emfrac, VldContext::GetDetector(), MadBase::GetEntry(), RecHeader::GetVldContext(), NtpSRStripPulseHeight::gev, NtpMCTruth::iaction, NtpSRCluster::id, NtpMCStdHep::IdHEP, NtpSREvent::index, MadQuantities::IsFidVtxEvt(), NtpMCStdHep::IstHEP, MadBase::LoadCluster(), MadBase::LoadEvent(), MadBase::LoadShower(), MadBase::LoadStdHep(), MadBase::LoadTruthForRecoTH(), NtpMCStdHep::mc, NtpSRShower::ncluster, NtpSREventSummary::nevent, NtpSREvent::nshower, NtpSREvent::ntrack, NtpMCStdHep::p4, NtpMCTruth::p4neu, NtpMCTruth::p4shw, NtpSRCluster::ph, NtpSRCluster::planeview, NtpSRCluster::probem, NtpSREvent::shw, NtpMCRecord::stdhep, NtpStRecord::stdhep, and NtpMCTruth::y. { char savename[256]; sprintf(savename,"Clu_%s.root",fileName); TFile *savefile = new TFile(savename,"RECREATE"); //all clusters ID TH1F *IDHist = new TH1F("IDHist","Cluster ID - CC",6,0,5); IDHist->SetXTitle("Cluster ID"); TH1F *IDHistNC = new TH1F("IDHistNC","Cluster ID - NC",6,0,5); IDHistNC->SetXTitle("Cluster ID"); //CC plots with y TH1F *yDist = new TH1F("yDist","y Distribution for CC Events",101,-0.005,1.005); TH1F *EMyDist = new TH1F("EMyDist","Fraction of CC Events with y",101,-0.005,1.005); EMyDist->SetXTitle("y"); TH1F *HADyDist = new TH1F("HADyDist","Fraction of CC Events with y",101,-0.005,1.005); HADyDist->SetXTitle("y"); TH1F *EMProbyDist = new TH1F("EMProbyDist","Fraction of CC Events with y",101,-0.005,1.005); EMProbyDist->SetXTitle("y"); TH1F *UVEMyDist = new TH1F("UVEMyDist","Fraction of CC Events with y",101,-0.005,1.005); UVEMyDist->SetXTitle("y"); TH1F *UVHADyDist = new TH1F("UVHADyDist","Fraction of CC Events with y",101,-0.005,1.005); UVHADyDist->SetXTitle("y"); TH1F *UVEMProbyDist = new TH1F("UVEMProbyDist","Fraction of CC Events with y",101,-0.005,1.005); UVEMProbyDist->SetXTitle("y"); //NC plots with y TH1F *yDistNC = new TH1F("yDistNC","y Distribution for NC Events",101,-0.005,1.005); TH1F *EMyDistNC = new TH1F("EMyDistNC","Fraction of NC Events with y",101,-0.005,1.005); EMyDistNC->SetXTitle("y"); TH1F *HADyDistNC = new TH1F("HADyDistNC","Fraction of NC Events with y", 101,-0.005,1.005); HADyDistNC->SetXTitle("y"); TH1F *EMProbyDistNC = new TH1F("EMProbyDistNC","Fraction of NC Events with y",101,-0.005,1.005); EMProbyDistNC->SetXTitle("y"); TH1F *UVEMyDistNC = new TH1F("UVEMyDistNC","Fraction of NC Events with y", 101,-0.005,1.005); UVEMyDistNC->SetXTitle("y"); TH1F *UVHADyDistNC = new TH1F("UVHADyDistNC","Fraction of NC Events with y", 101,-0.005,1.005); UVHADyDistNC->SetXTitle("y"); TH1F *UVEMProbyDistNC = new TH1F("UVEMProbyDistNC","Fraction of NC Events with y", 101,-0.005,1.005); UVEMProbyDistNC->SetXTitle("y"); //primary cluster IDs for NC and CC TH1F *UClusterIDCC = new TH1F("UClusterIDCC","U View Primary Cluster ID",6,-0.5,5.5); UClusterIDCC->SetXTitle("Cluster ID"); TH1F *UClusterIDNC = new TH1F("UClusterIDNC","U View Primary Cluster ID",6,-0.5,5.5); UClusterIDNC->SetXTitle("Cluster ID"); TH1F *VClusterIDCC = new TH1F("VClusterIDCC","V View Primary Cluster ID",6,-0.5,5.5); VClusterIDCC->SetXTitle("Cluster ID"); TH1F *VClusterIDNC = new TH1F("VClusterIDNC","V View Primary Cluster ID",6,-0.5,5.5); VClusterIDNC->SetXTitle("Cluster ID"); //Average #cluster with energy TH2F *AvgCluEnCC = new TH2F("AvgCluEnCC","Average Number of Clusters with Visible Energy - CC", 100,0,25,20,0,20); AvgCluEnCC->SetXTitle("Energy (GeV)"); AvgCluEnCC->SetYTitle("# of Clusters"); TH2F *AvgCluEnNC = new TH2F("AvgCluEnNC","Average Number of Clusters with Visible Energy - NC", 100,0,25,20,0,20); AvgCluEnNC->SetXTitle("Energy (GeV)"); AvgCluEnNC->SetYTitle("# of Clusters"); TH2F *AvgCluEnCCnox = new TH2F("AvgCluEnCCnox","Average Number of Clusters with Visible Energy (excl xtalk clu) - CC", 100,0,25,20,0,20); AvgCluEnCCnox->SetXTitle("Energy (GeV)"); AvgCluEnCCnox->SetYTitle("# of Clusters"); TH2F *AvgCluEnNCnox = new TH2F("AvgCluEnNCnox","Average Number of Clusters with Visible Energy (excl xtalk clu) - NC", 100,0,25,20,0,20); AvgCluEnNCnox->SetXTitle("Energy (GeV)"); AvgCluEnNCnox->SetYTitle("# of Clusters"); //Average #cluster with #stdhep final states TH2F *AvgCluSHCC = new TH2F("AvgCluSHCC","Average Number of Clusters Vs #StdHep Final States >0.1GeV - CC", 25,0,25,20,0,20); AvgCluSHCC->SetXTitle("#StdHep FS"); AvgCluSHCC->SetYTitle("# of Clusters"); TH2F *AvgCluSHNC = new TH2F("AvgCluSHNC","Average Number of Clusters Vs #StdHep Final States >0.1GeV - NC", 25,0,25,20,0,20); AvgCluSHNC->SetXTitle("#StdHep FS"); AvgCluSHNC->SetYTitle("# of Clusters"); TH2F *AvgCluSHCCnox = new TH2F("AvgCluSHCCnox","Average Number of Clusters Vs #StdHep Final States >0.5GeV (excl xtalk clu) - CC", 25,0,25,20,0,20); AvgCluSHCCnox->SetXTitle("#StdHep FS"); AvgCluSHCCnox->SetYTitle("# of Clusters"); TH2F *AvgCluSHNCnox = new TH2F("AvgCluSHNCnox","Average Number of Clusters Vs #StdHep Final States >0.5GeV (excl xtalk clu) - NC", 25,0,25,20,0,20); AvgCluSHNCnox->SetXTitle("#StdHep FS"); AvgCluSHNCnox->SetYTitle("# of Clusters"); //Average #cluster with #stdhep final states & energy TH3F *AvgCluSHEnergyCC = new TH3F("AvgCluSHEnergyCC","Average Number of Clusters Vs #StdHep Final States >0.1GeV vs Energy - CC", 100,0,25,25,0,25,50,0,50); AvgCluSHEnergyCC->SetXTitle("Energy (GeV)"); AvgCluSHEnergyCC->SetYTitle("#StdHep FS"); AvgCluSHEnergyCC->SetZTitle("# of Clusters"); TH3F *AvgCluSHEnergyNC = new TH3F("AvgCluSHEnergyNC","Average Number of Clusters Vs #StdHep Final States >0.1GeV vs Energy - NC", 100,0,25,25,0,25,50,0,50); AvgCluSHEnergyNC->SetXTitle("Energy (GeV)"); AvgCluSHEnergyNC->SetYTitle("#StdHep FS"); AvgCluSHEnergyNC->SetZTitle("# of Clusters"); TH3F *AvgCluSHEnergyCCnox = new TH3F("AvgCluSHEnergyCCnox","Average Number of Clusters Vs #StdHep Final States >0.5GeV (excl xtalk clu) vs Energy - CC", 100,0,25,25,0,25,50,0,50); AvgCluSHEnergyCCnox->SetXTitle("Energy (GeV)"); AvgCluSHEnergyCCnox->SetYTitle("#StdHep FS"); AvgCluSHEnergyCCnox->SetZTitle("# of Clusters"); TH3F *AvgCluSHEnergyNCnox = new TH3F("AvgCluSHEnergyNCnox","Average Number of Clusters Vs #StdHep Final States >0.5GeV (excl xtalk clu) vs Energy - NC", 100,0,25,25,0,25,50,0,50); AvgCluSHEnergyNCnox->SetXTitle("Energy (GeV)"); AvgCluSHEnergyNCnox->SetYTitle("#StdHep FS"); AvgCluSHEnergyNCnox->SetZTitle("# of Clusters"); //Reco'd EMFrac TH1F *EMFracCC = new TH1F("EMFracCC","Estimated EMFrac - CC Events",100,0,1); TH1F *EMFracNC = new TH1F("EMFracNC","Estimated EMFrac - NC Events",100,0,1); TH1F *EMFracCCDiff = new TH1F("EMFracCCDiff","(Truth - Estimated)/Truth EMFrac - CC Events",200,-2,2); TH1F *EMFracNCDiff = new TH1F("EMFracNCDiff","(Truth - Estimated)/Truth EMFrac - NC Events",200,-2,2); TH2F *EMFracCCY = new TH2F("EMFracCCY","Estimated EMFrac Vs Y - CC Events",100,0,1,100,0,1); TH2F *EMFracNCY = new TH2F("EMFracNCY","Estimated EMFrac Vs Y - NC Events",100,0,1,100,0,1); //2D reco vs true: TH2F *EMFrac_Truth_Reco = new TH2F("EMFracTrueReco","Reco vs Truth EMFrac", 100,0,1,100,0,1); EMFrac_Truth_Reco->SetXTitle("True EMFrac"); EMFrac_Truth_Reco->SetYTitle("Reco EMFrac"); TH2F *EMFrac_Truth_Reco_end0 = new TH2F("EMFracTrueReco_end0","Reco vs Truth EMFrac", 100,0,1,100,0,1); EMFrac_Truth_Reco_end0->SetXTitle("True EMFrac"); EMFrac_Truth_Reco_end0->SetYTitle("Reco EMFrac"); TH2F *EMFrac_Truth_Reco_end1 = new TH2F("EMFracTrueReco_end1","Reco vs Truth EMFrac", 100,0,1,100,0,1); EMFrac_Truth_Reco_end1->SetXTitle("True EMFrac"); EMFrac_Truth_Reco_end1->SetYTitle("Reco EMFrac"); TH2F *asymHistNC = new TH2F("asymHistNC","asymHistNC",1000,0,10,100,0,1); TH2F *asymHistCC = new TH2F("asymHistCC","asymHistCC",1000,0,10,100,0,1); TH2F *asymHistNC_NotPrime = new TH2F("asymHistNC_NotPrime","asymHistNC_NotPrime",1000,0,10,100,0,1); TH2F *asymHistCC_NotPrime = new TH2F("asymHistCC_NotPrime","asymHistCC_NotPrime",1000,0,10,100,0,1); TH2F *asymHist = new TH2F("asymHist","asymHist",1000,0,10,100,0,1); for(int i=startEnt;i<Nentries;i++){ if(i%1000==0) std::cout << float(i)*100./float(Nentries) << "% done" << std::endl; if(!this->GetEntry(i)) continue; if(eventSummary->nevent==0) continue; Int_t is_fid = 0; //fill tree once for each reconstructed event for(int j=0;j<eventSummary->nevent;j++){ if(!LoadEvent(j)) continue; //no event found if(ntpEvent->nshower==0) continue; //no shower found //get detector type: if(ntpHeader->GetVldContext().GetDetector()==Detector::kFar){ //fiducial criteria on vtx for far detector is_fid = 1; if(!IsFidVtxEvt(ntpEvent->index)) is_fid = 0; } else if(ntpHeader->GetVldContext().GetDetector()==Detector::kNear){ //fiducial criteria on vtx for near detector is_fid = 1; if(!IsFidVtxEvt(ntpEvent->index)) is_fid = 0; } //if(is_fid==0) continue; Int_t mcevent = -1; //check for a corresponding mc event if(LoadTruthForRecoTH(ntpEvent->index,mcevent)) { //loadtruth float nBigStdHep = 0; TClonesArray* pointStdhepArray = NULL; if(isST) pointStdhepArray = (strecord->stdhep); else pointStdhepArray = (mcrecord->stdhep); TClonesArray& stdhepArray = *pointStdhepArray; Int_t nStdHep = stdhepArray.GetEntries(); for(int k=0;k<nStdHep;k++){ LoadStdHep(k); if(ntpStdHep->mc==mcevent) { if(ntpStdHep->IstHEP==1){ //final state particles if(ntpStdHep->p4[3]>0.5&&!(abs(ntpStdHep->IdHEP)==12 || abs(ntpStdHep->IdHEP)==14 || abs(ntpStdHep->IdHEP)==16)) { nBigStdHep+=1; } } } } //Int_t shower = -1; //if(!LoadLargestShowerFromEvent(ntpEvent->index,shower)) continue; if(!LoadShower(ntpEvent->shw[0])) continue; int goodU = 0; int goodV = 0; int bestU = 0; int bestV = 0; double enU = 0; double enV = 0; double totU = 0; double totV = 0; for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; LoadCluster(index); if(ntpCluster->planeview==2){ totU+=ntpCluster->ph.gev; if(ntpCluster->ph.gev>enU) { enU = ntpCluster->ph.gev; bestU = index; if(ntpCluster->id==0) { if(ntpCluster->probem>0.2) goodU = 2; else goodU = 1; } else if(ntpCluster->id==1) goodU = -1; else goodU = 0; } } else if(ntpCluster->planeview==3){ totV+=ntpCluster->ph.gev; if(ntpCluster->ph.gev>enV) { enV = ntpCluster->ph.gev; bestV = index; if(ntpCluster->id==0){ if(ntpCluster->probem>0.2) goodV = 2; else goodV = 1; } else if(ntpCluster->id==1) goodV = -1; else goodV = 0; } } } double asymmetry = TMath::Abs(totU-totV); if((totU+totV)>0.) asymmetry /= (totU+totV); else asymmetry = 1; //cout << i << " " << j << " " << totU << " " << totV << " " << asymmetry << endl; if(ntpTruth->iaction==1){ yDist->Fill(ntpTruth->y); //if(totU<totV) asymHistCC->Fill(totU,asymmetry); //else asymHistCC->Fill(totV,asymmetry); if(asymmetry>0.8){ if(totU<totV) asymHistCC->Fill(totU+totV,totU); else asymHistCC->Fill(totU+totV,totV); } AvgCluEnCC->Fill(ntpTruth->p4neu[3],ntpShower->ncluster); //if(ntpTruth->p4neu[3]>2.5 && // ntpTruth->p4neu[3]<3.5) AvgCluSHCC->Fill(nBigStdHep,ntpShower->ncluster); AvgCluSHEnergyCC->Fill(ntpTruth->p4neu[3],nBigStdHep, ntpShower->ncluster); int cntU = 0; int cntV = 0; float EMFrac = 0; float EMFrac_end0 = 0; float EMFrac_end1 = 0; float EnTot = 0; float EnTot_end0 = 0; float EnTot_end1 = 0; for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; LoadCluster(index); IDHist->Fill(ntpCluster->id); if(ntpCluster->id==0 && ntpCluster->probem>0.2) { EMFrac+=ntpCluster->ph.gev; EnTot+=ntpCluster->ph.gev; if(ntpCluster->planeview==2){ EMFrac_end0+=ntpCluster->ph.gev; EnTot_end0+=ntpCluster->ph.gev; } if(ntpCluster->planeview==3){ EMFrac_end1+=ntpCluster->ph.gev; EnTot_end1+=ntpCluster->ph.gev; } } else { EnTot+=ntpCluster->ph.gev; if(ntpCluster->planeview==2){ EnTot_end0+=ntpCluster->ph.gev; } if(ntpCluster->planeview==3){ EnTot_end1+=ntpCluster->ph.gev; } } if(ntpCluster->ph.gev>0.25 && (ntpCluster->id==0 || ntpCluster->id==1 || ntpCluster->id==3)) { if(ntpCluster->planeview==2) cntU++; if(ntpCluster->planeview==3) cntV++; } if (index==bestU) UClusterIDCC->Fill(ntpCluster->id); if (index==bestV) VClusterIDCC->Fill(ntpCluster->id); } if(goodU>=1 || goodV>=1) EMyDist->Fill(ntpTruth->y); if(goodU==-1 || goodV==-1) HADyDist->Fill(ntpTruth->y); if(goodU==-1 && goodV==-1) UVHADyDist->Fill(ntpTruth->y); if(goodU>=1 && goodV>=1) UVEMyDist->Fill(ntpTruth->y); if(goodU==2 || goodV==2) EMProbyDist->Fill(ntpTruth->y); if(goodU==2 && goodV==2) UVEMProbyDist->Fill(ntpTruth->y); AvgCluEnCCnox->Fill(ntpTruth->p4neu[3],(cntU+cntV)/2.); //if(ntpTruth->p4neu[3]>2.5 && // ntpTruth->p4neu[3]<3.5) AvgCluSHCCnox->Fill(nBigStdHep,(cntU+cntV)/2.); AvgCluSHEnergyCCnox->Fill(ntpTruth->p4neu[3],nBigStdHep, (cntU+cntV)/2.); if(EnTot_end0!=0 || EnTot_end1!=0){ if(EnTot_end0==0) EMFrac = EMFrac_end1/EnTot_end1; else if(EnTot_end1==0) EMFrac = EMFrac_end0/EnTot_end0; else if(EMFrac_end0/EnTot_end0>EMFrac_end1/EnTot_end1) EMFrac = EMFrac_end0/EnTot_end0; else EMFrac = EMFrac_end1/EnTot_end1; //else EMFrac = (EMFrac_end0/EnTot_end0 + EMFrac_end1/EnTot_end1)/2.; EnTot = 1; } if(EnTot!=0) { EMFracCC->Fill(EMFrac/EnTot); EMFracCCY->Fill(ntpTruth->y,EMFrac/EnTot); if(ntpTruth->emfrac!=0) EMFracCCDiff->Fill((ntpTruth->emfrac - (EMFrac/EnTot))/ntpTruth->emfrac); else if(EMFrac==0) EMFracCCDiff->Fill(0); else EMFracCCDiff->Fill(2-EMFrac/EnTot); EMFrac_Truth_Reco->Fill(ntpTruth->emfrac,EMFrac/EnTot); if(EnTot_end0>0) EMFrac_Truth_Reco_end0->Fill(ntpTruth->emfrac,EMFrac_end0/EnTot_end0); if(EnTot_end1>0) EMFrac_Truth_Reco_end1->Fill(ntpTruth->emfrac,EMFrac_end1/EnTot_end1); } } else if(ntpTruth->iaction==0){ yDistNC->Fill(ntpTruth->y); if(totU<totV) asymHistNC->Fill(totU,asymmetry); else asymHistNC->Fill(totV,asymmetry); AvgCluEnNC->Fill(ntpTruth->p4shw[3],ntpShower->ncluster); //if(ntpTruth->p4neu[3]>3.0 && // ntpTruth->p4neu[3]<3.1) AvgCluSHNC->Fill(nBigStdHep,ntpShower->ncluster); AvgCluSHEnergyNC->Fill(ntpTruth->p4neu[3],nBigStdHep, ntpShower->ncluster); int cntU = 0; int cntV = 0; float EMFrac = 0; float EMFrac_end0 = 0; float EMFrac_end1 = 0; float EnTot = 0; float EnTot_end0 = 0; float EnTot_end1 = 0; for(int k=0; k<ntpShower->ncluster; k++){ int index = ntpShower->clu[k]; LoadCluster(index); IDHistNC->Fill(ntpCluster->id); if(ntpCluster->id==0 && ntpCluster->probem>0.2) { EMFrac+=ntpCluster->ph.gev; EnTot+=ntpCluster->ph.gev; if(ntpCluster->planeview==2){ EMFrac_end0+=ntpCluster->ph.gev; EnTot_end0+=ntpCluster->ph.gev; } if(ntpCluster->planeview==3){ EMFrac_end1+=ntpCluster->ph.gev; EnTot_end1+=ntpCluster->ph.gev; } } else { EnTot+=ntpCluster->ph.gev; if(ntpCluster->planeview==2){ EnTot_end0+=ntpCluster->ph.gev; } if(ntpCluster->planeview==3){ EnTot_end1+=ntpCluster->ph.gev; } } if(ntpCluster->ph.gev>0.5 && (ntpCluster->id==0 || ntpCluster->id==1 || ntpCluster->id==3)) { if(ntpCluster->planeview==2) cntU++; if(ntpCluster->planeview==3) cntV++; } if (index==bestU) UClusterIDNC->Fill(ntpCluster->id); if (index==bestV) VClusterIDNC->Fill(ntpCluster->id); } if(goodU>=1 || goodV>=1) EMyDistNC->Fill(ntpTruth->y); if(goodU==-1 || goodV==-1) HADyDistNC->Fill(ntpTruth->y); if(goodU==-1 && goodV==-1) UVHADyDistNC->Fill(ntpTruth->y); if(goodU>=1 && goodV>=1) UVEMyDistNC->Fill(ntpTruth->y); if(goodU==2 || goodV==2) EMProbyDistNC->Fill(ntpTruth->y); if(goodU==2 && goodV==2) UVEMProbyDistNC->Fill(ntpTruth->y); AvgCluEnNCnox->Fill(ntpTruth->p4shw[3],(cntU+cntV)/2.); //if(ntpTruth->p4neu[3]>3.0 && // ntpTruth->p4neu[3]<3.1) AvgCluSHNCnox->Fill(nBigStdHep,(cntU+cntV)/2.); AvgCluSHEnergyNCnox->Fill(ntpTruth->p4neu[3],nBigStdHep, (cntU+cntV)/2.); if(EnTot_end0!=0 || EnTot_end1!=0){ if(EnTot_end0==0) EMFrac = EMFrac_end1/EnTot_end1; else if(EnTot_end1==0) EMFrac = EMFrac_end0/EnTot_end0; else if(EMFrac_end0/EnTot_end0>EMFrac_end1/EnTot_end1) EMFrac = EMFrac_end0/EnTot_end0; else EMFrac = EMFrac_end1/EnTot_end1; //else EMFrac = (EMFrac_end0/EnTot_end0 + EMFrac_end1/EnTot_end1)/2.; EnTot = 1; } if(EnTot!=0) { EMFracNC->Fill(EMFrac/EnTot); EMFracNCY->Fill(ntpTruth->y,EMFrac/EnTot); if(ntpTruth->emfrac!=0) EMFracNCDiff->Fill((ntpTruth->emfrac - (EMFrac/EnTot))/ntpTruth->emfrac); else if(EMFrac==0) EMFracNCDiff->Fill(0); else EMFracNCDiff->Fill(2-EMFrac/EnTot); EMFrac_Truth_Reco->Fill(ntpTruth->emfrac,EMFrac/EnTot); if(EnTot_end0>0 && EnTot_end1==0) EMFrac_Truth_Reco_end0->Fill(ntpTruth->emfrac,EMFrac_end0/EnTot_end0); if(EnTot_end1>0 && EnTot_end0==0) EMFrac_Truth_Reco_end1->Fill(ntpTruth->emfrac,EMFrac_end1/EnTot_end1); } } //look at non primary showers for(int k=1;k<ntpEvent->nshower;k++){ totU = 0; totV = 0; Bool_t UHalo = true; Bool_t VHalo = true; if(LoadShower(ntpEvent->shw[k])){ for(int l=0; l<ntpShower->ncluster; l++){ if(LoadCluster(ntpShower->clu[l])){ if(ntpCluster->planeview==2){ totU+=ntpCluster->ph.gev; if(ntpCluster->id!=4) UHalo = false; } if(ntpCluster->planeview==3){ totV+=ntpCluster->ph.gev; if(ntpCluster->id!=4) VHalo = false; } } } } asymmetry = TMath::Abs(totU-totV); if((totU+totV)>0.) asymmetry /= (totU+totV); else asymmetry = 1; if(ntpTruth->iaction==1&&ntpEvent->ntrack==0){ if(UHalo || VHalo) { if(totU<totV) asymHistCC_NotPrime->Fill(totU,asymmetry); else asymHistCC_NotPrime->Fill(totV,asymmetry); } //if(asymmetry>0.8){ //if(totU<totV) asymHistCC_NotPrime->Fill(totU+totV,totU); //else asymHistCC_NotPrime->Fill(totU+totV,totV); //} } else if(ntpTruth->iaction==0){ if(totU<totV) asymHistNC_NotPrime->Fill(totU,asymmetry); else asymHistNC_NotPrime->Fill(totV,asymmetry); } } } else { //look at asymmetry for showers not associated with MC event Double_t totU = 0; Double_t totV = 0; if(LoadShower(ntpEvent->shw[0])){ for(int l=0; l<ntpShower->ncluster; l++){ if(LoadCluster(ntpShower->clu[l])){ if(ntpCluster->planeview==2){ totU+=ntpCluster->ph.gev; } if(ntpCluster->planeview==3){ totV+=ntpCluster->ph.gev; } } } } Double_t asymmetry = TMath::Abs(totU-totV); if((totU+totV)>0.) asymmetry /= (totU+totV); else asymmetry = 1; if(totU<totV) asymHist->Fill(totU,asymmetry); else asymHist->Fill(totV,asymmetry); } } } savefile->Write(); savefile->Close(); return; }

void MadCluAnalysis::Plotter (  char *  fileName  )

Definition at line 680 of file MadCluAnalysis.cxx. { TFile *file = new TFile(filename,"READ"); file->cd(); TH1F *IDHist = (TH1F*) file->Get("IDHist"); TH1F *IDHistNC = (TH1F*) file->Get("IDHistNC"); TH1F *yDist = (TH1F*) file->Get("yDist"); TH1F *EMyDist = (TH1F*) file->Get("EMyDist"); TH1F *HADyDist = (TH1F*) file->Get("HADyDist"); TH1F *EMProbyDist = (TH1F*) file->Get("EMProbyDist"); TH1F *UVHADyDist = (TH1F*) file->Get("UVHADyDist"); TH1F *UVEMyDist = (TH1F*) file->Get("UVEMyDist"); TH1F *UVEMProbyDist = (TH1F*) file->Get("UVEMProbyDist"); TH1F *yDistNC = (TH1F*) file->Get("yDistNC"); TH1F *EMyDistNC = (TH1F*) file->Get("EMyDistNC"); TH1F *HADyDistNC = (TH1F*) file->Get("HADyDistNC"); TH1F *EMProbyDistNC = (TH1F*) file->Get("EMProbyDistNC"); TH1F *UVHADyDistNC = (TH1F*) file->Get("UVHADyDistNC"); TH1F *UVEMyDistNC = (TH1F*) file->Get("UVEMyDistNC"); TH1F *UVEMProbyDistNC = (TH1F*) file->Get("UVEMProbyDistNC"); TH1F *UClusterIDCC = (TH1F*) file->Get("UClusterIDCC"); TH1F *UClusterIDNC = (TH1F*) file->Get("UClusterIDNC"); TH1F *VClusterIDCC = (TH1F*) file->Get("VClusterIDCC"); TH1F *VClusterIDNC = (TH1F*) file->Get("VClusterIDNC"); TH2F *AvgCluEnCC = (TH2F*) file->Get("AvgCluEnCC"); TH2F *AvgCluEnNC = (TH2F*) file->Get("AvgCluEnNC"); TH2F *AvgCluEnCCnox = (TH2F*) file->Get("AvgCluEnCCnox"); TH2F *AvgCluEnNCnox = (TH2F*) file->Get("AvgCluEnNCnox"); TH2F *AvgCluSHCC = (TH2F*) file->Get("AvgCluSHCC"); TH2F *AvgCluSHNC = (TH2F*) file->Get("AvgCluSHNC"); TH2F *AvgCluSHCCnox = (TH2F*) file->Get("AvgCluSHCCnox"); TH2F *AvgCluSHNCnox = (TH2F*) file->Get("AvgCluSHNCnox"); TH1F *EMFracCC = (TH1F*) file->Get("EMFracCC"); TH1F *EMFracNC = (TH1F*) file->Get("EMFracNC"); TH1F *EMFracCCDiff = (TH1F*) file->Get("EMFracCCDiff"); TH1F *EMFracNCDiff = (TH1F*) file->Get("EMFracNCDiff"); TH2F *EMFracCCY = (TH2F*) file->Get("EMFracCCY"); TH2F *EMFracNCY = (TH2F*) file->Get("EMFracNCY"); TH2F *EMFracTrueReco = (TH2F*) file->Get("EMFracTrueReco"); TH2F *EMFracTrueReco_end0 = (TH2F*) file->Get("EMFracTrueReco_end0"); TH2F *EMFracTrueReco_end1 = (TH2F*) file->Get("EMFracTrueReco_end1"); TCanvas *can = new TCanvas(); can->Divide(1,2); can->cd(1); EMyDist->Divide(yDist); HADyDist->Divide(yDist); EMProbyDist->Divide(yDist); UVEMyDist->Divide(yDist); UVHADyDist->Divide(yDist); UVEMProbyDist->Divide(yDist); EMyDist->SetLineColor(4); EMProbyDist->SetLineColor(6); HADyDist->SetLineColor(2); UVEMyDist->SetLineColor(3); UVHADyDist->SetLineColor(1); UVEMProbyDist->SetLineColor(7); EMyDist->Draw(); HADyDist->Draw("sames"); EMProbyDist->Draw("sames"); UVEMyDist->Draw("sames"); UVHADyDist->Draw("sames"); UVEMProbyDist->Draw("sames"); can->cd(2); EMyDistNC->Divide(yDistNC); HADyDistNC->Divide(yDistNC); EMProbyDistNC->Divide(yDistNC); UVEMyDistNC->Divide(yDistNC); UVHADyDistNC->Divide(yDistNC); UVEMProbyDistNC->Divide(yDistNC); EMyDistNC->SetLineColor(4); HADyDistNC->SetLineColor(2); EMProbyDistNC->SetLineColor(6); UVEMyDistNC->SetLineColor(3); UVHADyDistNC->SetLineColor(1); UVEMProbyDistNC->SetLineColor(7); EMyDistNC->Draw(); HADyDistNC->Draw("sames"); EMProbyDistNC->Draw("sames"); UVEMyDistNC->Draw("sames"); UVHADyDistNC->Draw("sames"); UVEMProbyDistNC->Draw("sames"); TCanvas *can1 = new TCanvas(); can1->Divide(1,2); can1->cd(1); UClusterIDCC->SetLineColor(4); float integ = UClusterIDCC->Integral(); UClusterIDCC->Scale(1./integ); UClusterIDCC->Draw(); UClusterIDNC->SetLineColor(2); integ = UClusterIDNC->Integral(); UClusterIDNC->Scale(1./integ); UClusterIDNC->Draw("sames"); can1->cd(2); VClusterIDCC->SetLineColor(4); integ = VClusterIDCC->Integral(); VClusterIDCC->Scale(1./integ); VClusterIDCC->Draw(); VClusterIDNC->SetLineColor(2); integ = VClusterIDNC->Integral(); VClusterIDNC->Scale(1./integ); VClusterIDNC->Draw("sames"); TCanvas *can2 = new TCanvas(); can2->Divide(2,2); can2->cd(1); AvgCluEnCC->Draw("COLZ"); TProfile *AvgCluEnCC_pfx = AvgCluEnCC->ProfileX(); AvgCluEnCC_pfx->Draw("same"); can2->cd(2); AvgCluEnNC->Draw("COLZ"); TProfile *AvgCluEnNC_pfx = AvgCluEnNC->ProfileX(); AvgCluEnNC_pfx->Draw("same"); can2->cd(3); AvgCluEnCCnox->Draw("COLZ"); TProfile *AvgCluEnCCnox_pfx = AvgCluEnCCnox->ProfileX(); AvgCluEnCCnox_pfx->Draw("same"); can2->cd(4); AvgCluEnNCnox->Draw("COLZ"); TProfile *AvgCluEnNCnox_pfx = AvgCluEnNCnox->ProfileX(); AvgCluEnNCnox_pfx->Draw("same"); TCanvas *can3 = new TCanvas(); can3->Divide(2,2); can3->cd(1); AvgCluSHCC->Draw("COLZ"); TProfile *AvgCluSHCC_pfx = AvgCluSHCC->ProfileX(); AvgCluSHCC_pfx->Draw("same"); can3->cd(2); AvgCluSHNC->Draw("COLZ"); TProfile *AvgCluSHNC_pfx = AvgCluSHNC->ProfileX(); AvgCluSHNC_pfx->Draw("same"); can3->cd(3); AvgCluSHCCnox->Draw("COLZ"); TProfile *AvgCluSHCCnox_pfx = AvgCluSHCCnox->ProfileX(); AvgCluSHCCnox_pfx->Draw("same"); can3->cd(4); AvgCluSHNCnox->Draw("COLZ"); TProfile *AvgCluSHNCnox_pfx = AvgCluSHNCnox->ProfileX(); AvgCluSHNCnox_pfx->Draw("same"); TCanvas *can4 = new TCanvas(); can4->cd(); IDHist->SetLineColor(4); IDHistNC->SetLineColor(2); IDHist->Draw(); IDHistNC->Draw("sames"); TCanvas *can5 = new TCanvas(); can5->Divide(2,2); can5->cd(1); EMFracCC->Draw(); can5->cd(2); EMFracNC->Draw(); can5->cd(3); EMFracCCDiff->Draw(); can5->cd(4); EMFracNCDiff->Draw(); TCanvas *can6 = new TCanvas(); can6->Divide(1,2); can6->cd(1); EMFracCCY->Draw("colz"); can6->cd(2); EMFracNCY->Draw("colz"); TCanvas *can7 = new TCanvas(); can7->Divide(1,3); can7->cd(1); EMFracTrueReco->Draw("COLZ"); TProfile *prof = EMFracTrueReco->ProfileX(); prof->Draw("sames"); can7->cd(2); EMFracTrueReco_end0->Draw("COLZ"); prof = EMFracTrueReco_end0->ProfileX(); prof->Draw("sames"); can7->cd(3); EMFracTrueReco_end1->Draw("COLZ"); prof = EMFracTrueReco_end1->ProfileX(); prof->Draw("sames"); return; }

Bool_t MadCluAnalysis::ReadNNFile (  char *  filename  )

Definition at line 2897 of file MadCluAnalysis.cxx. References fneuralNet. { TFile *f = new TFile(filename,"READ"); if(f->IsOpen()&&!f->IsZombie()){ fneuralNet = (TMultiLayerPerceptron*) f->Get("TMultiLayerPerceptron"); if(fneuralNet) cout << "Read NN file" << endl; } else return false; return true; }

Bool_t MadCluAnalysis::ReadPIDFile (  char *  filename  )

Definition at line 2875 of file MadCluAnalysis.cxx. References flikelihoodDists. { TFile *f = new TFile(filename,"READ"); if(f->IsOpen()&&!f->IsZombie()){ flikelihoodDists[0] = (TH1F*) f->Get("PhwIDU8"); flikelihoodDists[1] = (TH1F*) f->Get("PhwIDU10"); flikelihoodDists[2] = (TH1F*) f->Get("PhwProbEMU8"); flikelihoodDists[3] = (TH1F*) f->Get("PhwProbEMU10"); flikelihoodDists[4] = (TH1F*) f->Get("ChargeFracRMSU8"); flikelihoodDists[5] = (TH1F*) f->Get("ChargeFracRMSU10"); flikelihoodDists[6] = (TH1F*) f->Get("PhwIDV8"); flikelihoodDists[7] = (TH1F*) f->Get("PhwIDV10"); flikelihoodDists[8] = (TH1F*) f->Get("PhwProbEMV8"); flikelihoodDists[9] = (TH1F*) f->Get("PhwProbEMV10"); flikelihoodDists[10] = (TH1F*) f->Get("ChargeFracRMSV8"); flikelihoodDists[11] = (TH1F*) f->Get("ChargeFracRMSV10"); cout << "Read PID file" << endl; } else return false; return true; }

Float_t MadCluAnalysis::RecoAnalysisEnergy (  Int_t  event = 0  )  [protected, virtual]

Reimplemented from MadAnalysis. Definition at line 2456 of file MadCluAnalysis.cxx. References NtpSRStripPulseHeight::gev, MadBase::LoadEvent(), MadBase::LoadLargestShowerFromEvent(), MadBase::LoadLargestTrackFromEvent(), NtpSRShower::ph, and MadQuantities::RecoMuEnergy(). { if(!LoadEvent(event)) return 0; Int_t shower = -1; LoadLargestShowerFromEvent(event,shower); Int_t track = -1; LoadLargestTrackFromEvent(event,track); Double_t shw_en = 0; Double_t trk_en = 0; if(shower>=0) shw_en = ntpShower->ph.gev; if(track>=0) trk_en = RecoMuEnergy(0,track); return shw_en + trk_en; }

---

Member Data Documentation

Int_t MadCluAnalysis::fis_nue_pid

Definition at line 91 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

TH1F* MadCluAnalysis::flikelihoodDists[12]

Definition at line 98 of file MadCluAnalysis.h. Referenced by PID(), and ReadPIDFile().

TMultiLayerPerceptron* MadCluAnalysis::fneuralNet

Definition at line 99 of file MadCluAnalysis.h. Referenced by MadCluAnalysis(), PID(), and ReadNNFile().

Float_t MadCluAnalysis::fPID3

Definition at line 87 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Float_t MadCluAnalysis::fPID4

Definition at line 88 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Float_t MadCluAnalysis::fPID5

Definition at line 89 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_AveNumSSPlaneBest

Definition at line 74 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_AveNumSSPlaneU

Definition at line 72 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillAveNumSSPlane(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_AveNumSSPlaneV

Definition at line 73 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillAveNumSSPlane(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_BestProbEM

Definition at line 83 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillBestProbEM(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_ChargeFracRMSBest

Definition at line 70 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_ChargeFracRMSU

Definition at line 68 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillChargeFracRMS(), MadCluAnalysis(), and PID().

Double_t MadCluAnalysis::PAN_ChargeFracRMSV

Definition at line 69 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillChargeFracRMS(), MadCluAnalysis(), and PID().

Int_t MadCluAnalysis::PAN_NClusterU

Definition at line 93 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillNCluster(), and MadCluAnalysis().

Int_t MadCluAnalysis::PAN_NClusterV

Definition at line 94 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillNCluster(), and MadCluAnalysis().

Int_t MadCluAnalysis::PAN_NPhysClusterU

Definition at line 95 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillNCluster(), and MadCluAnalysis().

Int_t MadCluAnalysis::PAN_NPhysClusterV

Definition at line 96 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillNCluster(), and MadCluAnalysis().

Int_t MadCluAnalysis::PAN_NPhysShowerStripU

Definition at line 80 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Int_t MadCluAnalysis::PAN_NPhysShowerStripV

Definition at line 81 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Int_t MadCluAnalysis::PAN_NShowerStripBest

Definition at line 78 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Int_t MadCluAnalysis::PAN_NShowerStripU

Definition at line 76 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Int_t MadCluAnalysis::PAN_NShowerStripV

Definition at line 77 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_PHWAvgDevBest

Definition at line 66 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_PHWAvgDevU

Definition at line 64 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillPHWAvgDev(), MadCluAnalysis(), and PID().

Double_t MadCluAnalysis::PAN_PHWAvgDevV

Definition at line 65 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillPHWAvgDev(), MadCluAnalysis(), and PID().

Double_t MadCluAnalysis::PAN_PHWCluIDBest

Definition at line 58 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_PHWCluIDU

Definition at line 56 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillPHWCluID(), MadCluAnalysis(), and PID().

Double_t MadCluAnalysis::PAN_PHWCluIDV

Definition at line 57 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillPHWCluID(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_PHWProbEMBest

Definition at line 62 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_PHWProbEMU

Definition at line 60 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillPHWProbEM(), MadCluAnalysis(), and PID().

Double_t MadCluAnalysis::PAN_PHWProbEMV

Definition at line 61 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), FillPHWProbEM(), MadCluAnalysis(), and PID().

Double_t MadCluAnalysis::PAN_RatioProbEM

Definition at line 84 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

Double_t MadCluAnalysis::PAN_Weight

Definition at line 85 of file MadCluAnalysis.h. Referenced by AddUserBranches(), CallUserFunctions(), and MadCluAnalysis().

---

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

• MadCluAnalysis.h
• MadCluAnalysis.cxx

---